← Back to branch summary

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

« back to all changes in this revision

Viewing changes to tests/cc4/cc4.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/cc4/cc4.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
 
 
19
$TOL = 10**-8;
 
20
if($SRC_PATH ne "") {
 
21
  $REF_FILE = "$SRC_PATH/output.ref";
 
22
}
 
23
else {
 
24
  $REF_FILE = "output.ref";
 
25
}
 
26
$TEST_FILE = "output.dat";
 
27
$RESULT = "cc4.test";
 
28
 
 
29
system ("$PSICMD");
 
30
 
 
31
$FAIL = 0;
 
32
 
 
33
open(RE, ">$RESULT") || die "cannot open $RESULT: $!";
 
34
select (RE);
 
35
printf "CC4:\n";
 
36
 
 
37
if (abs(seek_nuc($REF_FILE) - seek_nuc($TEST_FILE)) > $TOL) {
 
38
  fail_test("Nuclear Repulsion Energy"); $FAIL = 1;
 
39
}
 
40
else {
 
41
  pass_test("Nuclear Repulsion Energy");
 
42
}
 
43
 
 
44
if (abs(seek_scf($REF_FILE) - seek_scf($TEST_FILE)) > $TOL) {
 
45
  fail_test("RHF Energy"); $FAIL = 1;
 
46
}
 
47
else { 
 
48
  pass_test("RHF Energy");
 
49
}
 
50
 
 
51
if (abs(seek_ccsd($REF_FILE) - seek_ccsd($TEST_FILE)) > $TOL) {
 
52
  fail_test("CCSD Energy"); $FAIL = 1;
 
53
}
 
54
else { 
 
55
  pass_test("CCSD Energy");
 
56
}
 
57
 
 
58
if (abs(seek_ccsd_t($REF_FILE) - seek_ccsd_t($TEST_FILE)) > $TOL) {
 
59
  fail_test("CCSD(T) Energy"); $FAIL = 1;
 
60
}
 
61
else {
 
62
  pass_test("CCSD(T) Energy");
 
63
}
 
64
 
 
65
close (RE);
 
66
 
 
67
system("cat $RESULT");
 
68
 
 
69
exit($FAIL);
 
70