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);
24
$REF_FILE = "$SRC_PATH/output.ref";
27
$REF_FILE = "output.ref";
29
$TEST_FILE = "output.dat";
30
$RESULT = "scf-freq.test";
36
open(RE, ">$RESULT") || die "cannot open $RESULT: $!";
40
if (abs(seek_nuc($REF_FILE) - seek_nuc($TEST_FILE)) > $TOL) {
41
fail_test("Nuclear Repulsion Energy"); $FAIL = 1;
44
pass_test("Nuclear Repulsion Energy");
47
if (abs(seek_scf($REF_FILE) - seek_scf($TEST_FILE)) > $TOL) {
48
fail_test("RHF Energy"); $FAIL = 1;
51
pass_test("RHF Energy");
54
@freq_ref = seek_anal_freq($REF_FILE,"Harmonic Frequency",$NDOF);
55
@freq_test = seek_anal_freq($TEST_FILE,"Harmonic Frequency",$NDOF);
57
if(!compare_arrays(\@freq_ref, \@freq_test, $NDOF, 1, $HTOL)) {
58
fail_test("SCF Frequencies"); $FAIL = 1;
61
pass_test("SCF Frequencies");
64
@int_ref = seek_int($REF_FILE,"Harmonic Frequency",$NDOF);
65
@int_test = seek_int($TEST_FILE,"Harmonic Frequency",$NDOF);
67
if(!compare_arrays(\@int_ref, \@int_test, $NDOF, 1, $ITOL)) {
68
fail_test("SCF Intensities"); $FAIL = 1;
71
pass_test("SCF Intensities");
76
system("cat $RESULT");