1
2#######################################################################
3#  TEST DESCRIPTION
4#######################################################################
5cat > cc3_noddy_hyp.info <<%EOF%
6   cc3_noddy_hyp
7   -------------
8   Molecule:         LIH
9   Theory Level:     CC3 / 6-31G
10   Test Purpose:     static and freq.-dep. first hyperpolarizability
11                     as computed with the noddy test code
12%EOF%
13
14#######################################################################
15#  INTEGRAL INPUT
16#######################################################################
17cat > cc3_noddy_hyp.mol <<%EOF%
18BASIS
196-31G
20blubber
21blubber
22    2    0  X  Y      1.00D-18
23        3.0   1    1    1
24Li    0.0000000000000000   0.000000000   1.30000000
25        1.0   1    1    1
26H     0.0000000000000000   0.000000000   0.00000000
27%EOF%
28
29#######################################################################
30#  DALTON INPUT
31#######################################################################
32cat > cc3_noddy_hyp.dal <<%EOF%
33**DALTON
34.RUN WAVE FUNCTIONS
35**INTEGRAL
36.DIPLEN
37**WAVE FUNCTION
38.CC
39*SCF INPUT
40.THRESH
41 1.0D-12
42.NODIIS
43*ORBITAL INPUT
44*CC INP
45.CC3
46.PRINT
47 1
48.THRENR
491.0D-14
50.THRLEQ
511.0D-14
52.MXLRV
53  1000
54.MAX IT
55  1000
56*CCLR
57.OPERATOR
58ZDIPLEN ZDIPLEN
59.FREQUE
60  3
610.00  0.10  0.20
62*CCQR
63.NOBMAT
64.OPERATOR
65ZDIPLEN ZDIPLEN ZDIPLEN
66.ORFREQ
67  2
68 0.00 0.10
69.SHGFRE
70  1
71 0.10
72*NODDY
73.XI
74.ETA
75.FMAT
76**END OF DALTON
77%EOF%
78#######################################################################
79
80
81#######################################################################
82#  CHECK SCRIPT
83#######################################################################
84echo '#!/bin/ksh
85log=$1
86
87if [ `uname` = Linux ]; then
88   GREP="egrep -a"
89else
90   GREP="egrep"
91fi
92
93# SCF, MP2 and CC3 energies:
94CRIT1=`$GREP "Total SCF   energy: * -7\.5744842494" $log | wc -l`
95CRIT2=`$GREP "Total MP2   energy: * -7\.5922238176" $log | wc -l`
96CRIT3=`$GREP "Total CC3   energy: * -7\.6000272811" $log | wc -l`
97TEST[1]=`expr $CRIT1 \+ $CRIT2 \+ $CRIT3`
98CTRL[1]=6
99ERROR[1]="ENERGIES NOT CORRECT"
100
101
102# CC3 polarizability :
103CRIT1=`$GREP  "ZDIPLEN *.unrel.. *..0000 *ZDIPLEN *.unrel.. *..0000 * 12\.684665" $log | wc -l`
104CRIT2=`$GREP  "...1000 * ..1000 * 230.56897" $log | wc -l`
105CRIT3=`$GREP  "...2000 * ..2000 * 8.8105639" $log | wc -l`
106TEST[2]=`expr $CRIT1  \+ $CRIT2 \+ $CRIT3`
107CTRL[2]=3
108ERROR[2]="CC3 POLARIZABILITY NOT CORRECT"
109
110# CC3 hyperpolarizability :
111CRIT1=`$GREP  "ZDIPLEN  .unrel.. *..0000   ZDIPLEN  .unrel.. *..0000   ZDIPLEN  .unrel.. *..0000 * -1050.3968" $log | wc -l`
112CRIT2=`$GREP  "..0000 * ..1000 * ...1000 * -401811.17" $log | wc -l`
113CRIT3=`$GREP  "..2000 * ..1000 * ...1000 * -115.90926" $log | wc -l`
114TEST[3]=`expr $CRIT1  \+ $CRIT2 \+ $CRIT3`
115CTRL[3]=3
116ERROR[3]="CC3 HYPERPOLARIZABILITY NOT CORRECT"
117
118
119PASSED=1
120for i in 1 2 3
121do
122   if [ ${TEST[i]} -ne ${CTRL[i]} ]; then
123     echo "${ERROR[i]} ( test = ${TEST[i]}; control = ${CTRL[i]} ); "
124     PASSED=0
125   fi
126done
127
128if [ $PASSED -eq 1 ]
129then
130  echo TEST ENDED PROPERLY
131  exit 0
132else
133  echo THERE IS A PROBLEM
134  exit 1
135fi
136
137' > cc3_noddy_hyp.check
138#######################################################################
139