3
while ($ARGV = shift) {
4
if ("$ARGV" eq "-q") { $QUIET = 1; }
5
elsif("$ARGV" eq "-i") { $SRC_PATH = shift; }
6
elsif("$ARGV" eq "-x") { $EXEC_PATH = shift; }
10
require($SRC_PATH . "/../psitest.pl");
13
require("../psitest.pl");
16
# build the command for the psi3 driver
17
$PSICMD = build_psi_cmd($QUIET, $SRC_PATH, $EXEC_PATH);
23
$REF_FILE = "$SRC_PATH/file11.ref";
26
$REF_FILE = "file11.ref";
28
$TEST_FILE = "psi.file11.dat";
31
$REF_OUT = "$SRC_PATH/output.ref";
34
$REF_OUT = "output.ref";
36
$TEST_OUT = "output.dat";
37
$RESULT = "scf-freq-numer.test";
44
$natom = seek_natom_file11($REF_FILE,"SCF");
46
open(RE, ">$RESULT") || die "cannot open $RESULT $!";
48
printf "SCF-FREQ-NUMER:\n";
50
if(abs(seek_energy_file11($REF_FILE,"SCF") - seek_energy_file11($TEST_FILE,"SCF")) > $TOL) {
51
fail_test("SCF energy"); $FAIL = 1;
54
pass_test("SCF energy");
57
@geom_ref = seek_geom_file11($REF_FILE, "SCF");
58
@geom_test = seek_geom_file11($TEST_FILE, "SCF");
59
if(!compare_arrays(\@geom_ref, \@geom_test, $natom, 3, $GTOL)) {
60
fail_test("SCF Geometry"); $FAIL = 1;
63
pass_test("SCF Geometry");
66
@grad_ref = seek_grad_file11($REF_FILE, "SCF");
67
@grad_test = seek_grad_file11($TEST_FILE, "SCF");
69
if(!compare_arrays(\@grad_ref, \@grad_test, $natom, 3, $GTOL)) {
70
fail_test("SCF Gradient"); $FAIL = 1;
73
pass_test("SCF Gradient");
76
@freq_ref = seek_findif_freq($REF_OUT,"Harmonic Vibrational Frequencies",$NDOF);
77
@freq_test = seek_findif_freq($TEST_OUT,"Harmonic Vibrational Frequencies",$NDOF);
79
if(!compare_arrays(\@freq_ref, \@freq_test, $NDOF, 1, $HTOL)) {
80
fail_test("SCF Frequencies"); $FAIL = 1;
83
pass_test("SCF Frequencies");
88
system("cat $RESULT");