← Back to branch summary

~ubuntu-branches/ubuntu/quantal/psicode/quantal

« back to all changes in this revision

Viewing changes to tests/cc13/cc13.test.pl

  • Committer: Bazaar Package Importer
  • Author(s): Michael Banck
  • Date: 2006-09-10 14:01:33 UTC
  • Revision ID: james.westby@ubuntu.com-20060910140133-ib2j86trekykfsfv
Tags: upstream-3.2.3
Import upstream version 3.2.3

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/cc13/cc13.test.pl
 
1
#!/usr/bin/perl  
 
2
 
 
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; }
 
7
}
 
8
 
 
9
if($SRC_PATH ne "") {
 
10
  require($SRC_PATH . "/../psitest.pl");
 
11
}
 
12
else {
 
13
  require("../psitest.pl");
 
14
}
 
15
 
 
16
# build the command for the psi3 driver
 
17
$PSICMD = build_psi_cmd($QUIET, $SRC_PATH, $EXEC_PATH);
 
18
$TOL = 10**-8;
 
19
$GTOL = 10**-8;
 
20
if($SRC_PATH ne "") {
 
21
  $REF_FILE = "$SRC_PATH/file11.ref";
 
22
}
 
23
else {
 
24
  $REF_FILE = "file11.ref";
 
25
}
 
26
$TEST_FILE = "psi.file11.dat";
 
27
$RESULT = "cc13.test";
 
28
 
 
29
system ("$PSICMD");
 
30
 
 
31
$FAIL = 0;
 
32
 
 
33
$natom = seek_natom_file11($REF_FILE,"CCSD");
 
34
 
 
35
open(RE, ">$RESULT") || die "cannot open $RESULT $!"; 
 
36
select (RE);
 
37
printf "CC13:\n";
 
38
 
 
39
if(abs(seek_energy_file11($REF_FILE,"CCSD") - seek_energy_file11($TEST_FILE,"CCSD")) > $TOL) {
 
40
  fail_test("CCSD energy"); $FAIL = 1;
 
41
}
 
42
else {
 
43
  pass_test("CCSD energy");
 
44
}
 
45
 
 
46
@geom_ref = seek_geom_file11($REF_FILE, "CCSD");
 
47
@geom_test = seek_geom_file11($TEST_FILE, "CCSD");
 
48
if(!compare_arrays(\@geom_ref, \@geom_test, $natom, 3, $GTOL)) {
 
49
  fail_test("CCSD Geometry"); $FAIL = 1;
 
50
}
 
51
else {
 
52
  pass_test("CCSD Geometry");
 
53
}
 
54
 
 
55
@grad_ref = seek_grad_file11($REF_FILE, "CCSD");
 
56
@grad_test = seek_grad_file11($TEST_FILE, "CCSD");
 
57
 
 
58
if(!compare_arrays(\@grad_ref, \@grad_test, $natom, 3, $GTOL)) {
 
59
  fail_test("CCSD Gradient"); $FAIL = 1;
 
60
}
 
61
else {
 
62
  pass_test("CCSD Gradient");
 
63
}
 
64
 
 
65
close (RE);
 
66
 
 
67
system("cat $RESULT");
 
68
 
 
69
exit($FAIL);
 
70