1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
|
#!/bin/sh
currtestdirname=tmpdir_test_fmm
if test "$top_builddir" = ""; then
top_builddir=..
fi
if test "$top_srcdir" = ""; then
top_srcdir=..
fi
. "$top_srcdir"/test/functions
# Run each test in a separate directory, to allow "make check -j" to work properly.
currdir=`pwd` ; cd $top_builddir ; top_builddir_pwd=`pwd` ; cd $currdir
rm -rf $currtestdirname ; mkdir $currtestdirname ; cd $currtestdirname
ln -s "$top_builddir_pwd"/source/ergo ./ergo
if test `./ergo -e precision` = 'single'; then
echo SKIPPED
cd .. ; rm -r $currtestdirname
exit 0
fi
errorfilename=ergoscf.out.error.fmm
echo
echo Testing four_h2o HF/6-31G using FMM
rm -f ergoscf.out
./ergo <<EOINPUT > /dev/null
molecule_inline Angstrom
H 8.675000 0.000000 -16.640000
H 9.579936 0.000000 -17.920455
H 7.203213 0.037757 -17.673336
H 8.108149 0.037757 -18.953791
H 9.241330 0.571501 -2.493601
H 10.146266 0.571501 -3.774055
H 7.134485 0.659897 -2.982496
H 8.039421 0.659897 -4.262951
O 8.675000 0.000000 -17.600000
O 7.203213 0.037757 -18.633336
O 9.241330 0.571501 -3.453601
O 7.134485 0.659897 -3.942496
EOF
basis = "6-31G"
run "HF"
EOINPUT
if
check_final_energy -303.7606521 1e-7 ;
then
echo OK
else
echo ERROR
mv ergoscf.out $errorfilename
echo output file saved as $errorfilename
exit 1
fi
rm ergoscf.out
echo Testing four_h2o HF/6-31G* using FMM
rm -f ergoscf.out
./ergo <<EOINPUT > /dev/null
molecule_inline Angstrom
H 8.675000 0.000000 -16.640000
H 9.579936 0.000000 -17.920455
H 7.203213 0.037757 -17.673336
H 8.108149 0.037757 -18.953791
H 9.241330 0.571501 -2.493601
H 10.146266 0.571501 -3.774055
H 7.134485 0.659897 -2.982496
H 8.039421 0.659897 -4.262951
O 8.675000 0.000000 -17.600000
O 7.203213 0.037757 -18.633336
O 9.241330 0.571501 -3.453601
O 7.134485 0.659897 -3.942496
EOF
basis = "6-31Gs"
initial_density = "density.bin"
run "HF"
EOINPUT
if
check_final_energy -303.8352704 1e-7 ;
then
echo OK
else
echo ERROR
mv ergoscf.out $errorfilename
echo output file saved as $errorfilename
exit 1
fi
rm ergoscf.out
rm density.bin
cd ..
rm -r $currtestdirname
echo
echo FMM tests completed successfully!
echo
|
