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
105
106
107
108
109
110
111
|
o1 = cos(q)
o2 = 2.00000000000000d0*q
o3 = cos(o2)
o4 = dnnnnngtil
o5 = sin(q)
o6 = o5**2
o7 = o5*o6
o8 = dnnngtil
o9 = dngtil
o10 = dnnnnnnngtil
o11 = o6**2
o12 = 4.0000000000000d0*q
o13 = cos(o12)
o14 = dnnnngtil
o15 = dnngtil
o16 = o1**2
o17 = dnnnnnngtil
o18 = o11*o6
o19 = gtil
o20 = byJ**2
o21 = o14**2
o22 = o3**2
o23 = o13**2
o24 = o15**2
o25 = o4**2
o26 = o8**2
o27 = o9**2
o28 = o16**2
o29 = o17**2
o30 = o10**2
o31 = o11**2
o32 = o31*o6
o33 = o19**2
c exc33(i,j,k) = 0
c exc32(i,j,k) = 1.33333333333333d0*o1*o10*o11*o5 - 2.083333333333
c & 33d-2*o1*(3.5200000000000d2 - 9.2800000000000d2*o3)*o4*o7 + 1.66
c & 666666666667d-1*o1*(1.15000000000000d2 - 4.1500000000000d2*o3)*o
c & 7*o8 - 1.56250000000000d0*o1*(8.0000000000000d0 - 4.000000000000
c & 0d1*o3)*o7*o9
c exc31(i,j,k) = -8.0000000000000d0*o11*o16*o17 + 1.33333333333333
c & d0*o17*o18 + 3.00000000000000d0*byJ*o6 - 2.08333333333333d-1*o15*(
c & -1.11000000000000d2 - 5.8100000000000d2*o13 - 2.68000000000000d2
c & *o3)*o6 - 5.2083333333333d-3*o14*(6.7200000000000d2 + 6.49600000
c & 00000d3*o13 + 2.04800000000000d3*o3)*o6 + o19*(5.0000000000000d1
c & *o11 + 7.5000000000000d2*o11*o16 - 1.25000000000000d2*o18 - 2.00
c & 000000000000d2*o16*o6)
c exc22(i,j,k) = 0
c exc21(i,j,k) = 0
c exc11(i,j,k) = 0
t1 = -9.6000000000000d1*byJ*o11*o16*o17 + 1.12000000000000d2*o11*o
& 14*o16*o17 + 1.08266666666667d3*o11*o13*o14*o16*o17 - 7.40000000
& 00000d2*o11*o15*o16*o17 - 3.8733333333333d3*o11*o13*o15*o16*o17
& + 1.60000000000000d1*byJ*o17*o18 - 1.86666666666667d1*o14*o17*o18
& - 1.80444444444444d2*o13*o14*o17*o18 + 1.23333333333333d2*o15*o1
& 7*o18 + 6.4555555555556d2*o13*o15*o17*o18 + 6.0000000000000d2*byJ*
& o11*o19 - 7.0000000000000d2*o11*o14*o19 - 6.7666666666667d3*o11*
& o13*o14*o19 + 4.6250000000000d3*o11*o15*o19 + 2.42083333333333d4
& *o11*o13*o15*o19 + 9.0000000000000d3*byJ*o11*o16*o19 - 1.050000000
& 00000d4*o11*o14*o16*o19 - 1.01500000000000d5*o11*o13*o14*o16*o19
& + 6.9375000000000d4*o11*o15*o16*o19 + 3.6312500000000d5*o11*o13
& *o15*o16*o19 - 1.50000000000000d3*byJ*o18*o19 + 1.75000000000000d3
& *o14*o18*o19 + 1.69166666666667d4*o13*o14*o18*o19 - 1.1562500000
& 0000d4*o15*o18*o19 - 6.0520833333333d4*o13*o15*o18*o19 - 2.66666
& 666666667d3*o16*o17*o18*o19 + 1.07555555555556d2*o11*o16*o25
t2 = 7.4755555555556d2*o11*o16*o22*o25 + 7.3472222222222d2*o11*o
& 16*o26 + 9.5680555555556d3*o11*o16*o22*o26 + 3.12500000000000d2*
& o11*o16*o27 + 7.8125000000000d3*o11*o16*o22*o27 + 6.400000000000
& 0d3*o11*o17*o19*o28 - 2.40000000000000d4*o17*o18*o19*o28 + 1.280
& 00000000000d2*o18*o28*o29 + 3.4133333333333d2*o11*o14*o16*o17*o3
& - 1.78666666666667d3*o11*o15*o16*o17*o3 - 5.6888888888889d1*o14
& *o17*o18*o3 + 2.97777777777778d2*o15*o17*o18*o3 - 2.133333333333
& 33d3*o11*o14*o19*o3 + 1.11666666666667d4*o11*o15*o19*o3 - 3.2000
& 000000000d4*o11*o14*o16*o19*o3 + 1.67500000000000d5*o11*o15*o16*
& o19*o3 + 5.3333333333333d3*o14*o18*o19*o3 - 2.79166666666667d4*o
& 15*o18*o19*o3 - 5.6711111111111d2*o11*o16*o25*o3 - 5.30277777777
& 78d3*o11*o16*o26*o3 - 3.12500000000000d3*o11*o16*o27*o3 + 2.6666
& 6666666667d2*o17*o19*o31 + 8.0000000000000d3*o16*o17*o19*o31 - 4
& .2666666666667d1*o16*o29*o31 + 3.5555555555556d0*o16*o30*o31 - 6
& .6666666666667d2*o17*o19*o32 + 3.5555555555556d0*o29*o32 - 4.000
& 0000000000d4*o11*o16*o33
t3 = 5.0000000000000d3*o18*o33 + 2.50000000000000d5*o16*o18*o33
& - 6.0000000000000d5*o11*o28*o33 + 1.12500000000000d6*o18*o28*o33
& - 2.50000000000000d4*o31*o33 - 3.7500000000000d5*o16*o31*o33 +
& 3.12500000000000d4*o32*o33 - 3.9111111111111d1*o10*o16*o18*o4 +
& 1.03111111111111d2*o10*o16*o18*o3*o4 - 4.2000000000000d1*byJ*o14*o
& 6 - 4.0600000000000d2*byJ*o13*o14*o6 + 2.77500000000000d2*byJ*o15*o6
& + 1.45250000000000d3*byJ*o13*o15*o6 - 3.2375000000000d2*o14*o15*o
& 6 - 4.8241666666667d3*o13*o14*o15*o6 - 2.40000000000000d3*byJ*o16*
& o19*o6 + 2.80000000000000d3*o14*o16*o19*o6 + 2.70666666666667d4*
& o13*o14*o16*o19*o6 - 1.85000000000000d4*o15*o16*o19*o6 - 9.68333
& 33333333d4*o13*o15*o16*o19*o6 + 1.80000000000000d1*o20*o6 + 2.45
& 000000000000d1*o21*o6 + 4.7366666666667d2*o13*o21*o6 - 2.3822222
& 2222222d3*o14*o15*o22*o6 + 2.27555555555556d2*o21*o22*o6 - 1.638
& 09722222222d4*o14*o15*o23*o6 + 2.28938888888889d3*o21*o23*o6
t4 = 1.06953125000000d3*o24*o6 + 1.11963541666667d4*o13*o24*o6 +
& 6.2347222222222d3*o22*o24*o6 + 2.93021701388889d4*o23*o24*o6 -
& 1.28000000000000d2*byJ*o14*o3*o6 + 6.7000000000000d2*byJ*o15*o3*o6 -
& 1.76833333333333d3*o14*o15*o3*o6 - 1.27205555555556d4*o13*o14*o
& 15*o3*o6 + 8.5333333333333d3*o14*o16*o19*o3*o6 - 4.4666666666667
& d4*o15*o16*o19*o3*o6 + 1.49333333333333d2*o21*o3*o6 + 1.44355555
& 555556d3*o13*o21*o3*o6 + 5.1645833333333d3*o24*o3*o6 + 2.7032638
& 8888889d4*o13*o24*o3*o6 + 8.0000000000000d4*o28*o33*o6 + 1.02222
& 222222222d2*o10*o16*o18*o8 - 3.6888888888889d2*o10*o16*o18*o3*o8
& - 5.6222222222222d2*o11*o16*o4*o8 - 5.3488888888889d3*o11*o16*o
& 22*o4*o8 + 3.5111111111111d3*o11*o16*o3*o4*o8 - 6.6666666666667d
& 1*o10*o16*o18*o9 + 3.3333333333333d2*o10*o16*o18*o3*o9 + 3.66666
& 66666667d2*o11*o16*o4*o9 + 4.8333333333333d3*o11*o16*o22*o4*o9 -
& 2.80000000000000d3*o11*o16*o3*o4*o9 - 9.5833333333333d2*o11*o16
& *o8*o9 - 1.72916666666667d4*o11*o16*o22*o8*o9 + 8.2500000000000d
& 3*o11*o16*o3*o8*o9
ksq(i,j,k) = t1 + t2 + t3 + t4
|
