1
# Make sure color.rb is included
4
PSITEST_ETOL = 10e-8; # Default test criterion for energies
5
PSITEST_ENUCTOL = 10e-10; # Check nuclear repulsion energy tighter than other energies
7
def test_scf_energy(expected)
10
# Make sure SCF is in the checkpoint
11
if Psi::Chkpt::exist?("SCF energy")
12
if (Psi::Chkpt::escf - expected).abs < PSITEST_ETOL
13
puts "SCF final energy: " + green("PASSED")
15
puts "SCF final energy: " + red("FAILED")
16
puts "Obtained: #{Psi::Chkpt::escf}"
17
puts "Expected: #{expected}"
21
puts red("Error: ") + "SCF energy not found in checkpoint."
28
def test_ccsd_energy(expected)
31
# Make sure CCSD is in the checkpoint
32
if Psi::Chkpt::exist?("CCSD Energy")
33
if (Psi::Chkpt::eccsd - expected).abs < PSITEST_ETOL
34
puts "CCSD final energy: " + green("PASSED")
36
puts "CCSD final energy: " + red("FAILED")
37
puts "Obtained: #{Psi::Chkpt::eccsd}"
38
puts "Expected: #{expected}"
42
puts red("Error: ") + "CCSD energy not found in checkpoint."
49
def test_t_energy(expected)
52
# Make sure (T) is in the checkpoint
53
if Psi::Chkpt::exist?("(T) Energy")
54
if (Psi::Chkpt::e_t - expected).abs < PSITEST_ETOL
55
puts "(T) final energy: " + green("PASSED")
57
puts "(T) final energy: " + red("FAILED")
58
puts "Obtained: #{Psi::Chkpt::e_t}"
59
puts "Expected: #{expected}"
63
puts red("Error: ") + "(T) energy not found in checkpoint."
70
def test_nuclear_repulsion(expected)
73
# Make sure it is found in the checkpoint
74
if Psi::Chkpt::exist?("Nuclear rep. energy")
75
if (Psi::Chkpt::enuc - expected).abs < PSITEST_ENUCTOL
76
puts "Nuclear repulsion energy: " + green("PASSED")
78
puts "Nuclear repulsion energy: " + red("FAILED")
79
puts "Obtained: #{Psi::Chkpt::enuc}"
80
puts "Expected: #{expected}"
84
puts red("Error:") + "Nuclear repulsion energy not found in checkpoint."
91
def test_total_energy(expected)
94
# Make sure it is found in the checkpoint
95
if Psi::Chkpt::exist?("Total energy")
96
if (Psi::Chkpt::etot - expected).abs < PSITEST_ETOL
97
puts "Total energy: " + green("PASSED")
99
puts "Total energy: " + red("FAILED")
100
puts "Obtained: #{Psi::Chkpt::etot}"
101
puts "Expected: #{expected}"
105
puts red("Error:") + "Total energy not found in checkpoint."