blob: a590721c25e0fbd58957f84b44ef8cf33d157c66 (
plain)
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
|
o1 = cos(qgrd(j))
o2 = 2.00000000000000d0*qgrd(j)
o3 = cos(o2)
o4 = dnnnnngtil
o5 = sin(qgrd(j))
o6 = o5**2
o7 = o5*o6
o8 = dnnngtil
o9 = dngtil
o10 = dnnnnnnngtil
o11 = o6**2
o12 = 4.0000000000000d0*qgrd(j)
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
t1 = 2.88000000000000d1*byJ*o11*o16*o17 + 1.00800000000000d1*o11*o
& 14*o16*o17 + 9.7440000000000d1*o11*o13*o14*o16*o17 - 6.660000000
& 0000d1*o11*o15*o16*o17 - 3.4860000000000d2*o11*o13*o15*o16*o17 -
& 4.8000000000000d0*byJ*o17*o18 - 1.68000000000000d0*o14*o17*o18 -
& 1.62400000000000d1*o13*o14*o17*o18 + 1.11000000000000d1*o15*o17*
& o18 + 5.8100000000000d1*o13*o15*o17*o18 - 1.80000000000000d2*byJ*o
& 11*o19 - 6.3000000000000d1*o11*o14*o19 - 6.0900000000000d2*o11*o
& 13*o14*o19 + 4.1625000000000d2*o11*o15*o19 + 2.17875000000000d3*
& o11*o13*o15*o19 - 2.70000000000000d3*byJ*o11*o16*o19 - 9.450000000
& 0000d2*o11*o14*o16*o19 - 9.1350000000000d3*o11*o13*o14*o16*o19 +
& 6.2437500000000d3*o11*o15*o16*o19 + 3.2681250000000d4*o11*o13*o
& 15*o16*o19 + 4.5000000000000d2*byJ*o18*o19 + 1.57500000000000d2*o1
& 4*o18*o19 + 1.52250000000000d3*o13*o14*o18*o19 - 1.0406250000000
& 0d3*o15*o18*o19 - 5.4468750000000d3*o13*o15*o18*o19 - 2.40000000
& 000000d2*o16*o17*o18*o19 + 9.6800000000000d0*o11*o16*o25
t2 = 6.7280000000000d1*o11*o16*o22*o25 + 6.6125000000000d1*o11*o
& 16*o26 + 8.6112500000000d2*o11*o16*o22*o26 + 2.81250000000000d1*
& o11*o16*o27 + 7.0312500000000d2*o11*o16*o22*o27 + 5.760000000000
& 0d2*o11*o17*o19*o28 - 2.16000000000000d3*o17*o18*o19*o28 + 1.152
& 00000000000d1*o18*o28*o29 + 3.07200000000000d1*o11*o14*o16*o17*o
& 3 - 1.60800000000000d2*o11*o15*o16*o17*o3 - 5.1200000000000d0*o1
& 4*o17*o18*o3 + 2.68000000000000d1*o15*o17*o18*o3 - 1.92000000000
& 000d2*o11*o14*o19*o3 + 1.00500000000000d3*o11*o15*o19*o3 - 2.880
& 00000000000d3*o11*o14*o16*o19*o3 + 1.50750000000000d4*o11*o15*o1
& 6*o19*o3 + 4.8000000000000d2*o14*o18*o19*o3 - 2.51250000000000d3
& *o15*o18*o19*o3 - 5.1040000000000d1*o11*o16*o25*o3 - 4.772500000
& 0000d2*o11*o16*o26*o3 - 2.81250000000000d2*o11*o16*o27*o3 + 2.40
& 000000000000d1*o17*o19*o31 + 7.2000000000000d2*o16*o17*o19*o31 -
& 3.8400000000000d0*o16*o29*o31 + 3.2000000000000d-1*o16*o30*o31
& - 6.0000000000000d1*o17*o19*o32 + 3.2000000000000d-1*o29*o32 - 3
& .6000000000000d3*o11*o16*o33
t3 = 4.5000000000000d2*o18*o33 + 2.25000000000000d4*o16*o18*o33
& - 5.4000000000000d4*o11*o28*o33 + 1.01250000000000d5*o18*o28*o33
& - 2.25000000000000d3*o31*o33 - 3.3750000000000d4*o16*o31*o33 +
& 2.81250000000000d3*o32*o33 - 3.5200000000000d0*o10*o16*o18*o4 +
& 9.2800000000000d0*o10*o16*o18*o3*o4 + 1.26000000000000d1*byJ*o14*o
& 6 + 1.21800000000000d2*byJ*o13*o14*o6 - 8.3250000000000d1*byJ*o15*o6
& - 4.3575000000000d2*byJ*o13*o15*o6 - 2.91375000000000d1*o14*o15*o
& 6 - 4.3417500000000d2*o13*o14*o15*o6 + 7.2000000000000d2*byJ*o16*o
& 19*o6 + 2.52000000000000d2*o14*o16*o19*o6 + 2.43600000000000d3*o
& 13*o14*o16*o19*o6 - 1.66500000000000d3*o15*o16*o19*o6 - 8.715000
& 0000000d3*o13*o15*o16*o19*o6 + 1.80000000000000d1*o20*o6 + 2.205
& 00000000000d0*o21*o6 + 4.2630000000000d1*o13*o21*o6 - 2.14400000
& 000000d2*o14*o15*o22*o6 + 2.04800000000000d1*o21*o22*o6 - 1.4742
& 8750000000d3*o14*o15*o23*o6 + 2.06045000000000d2*o21*o23*o6
t4 = 9.6257812500000d1*o24*o6 + 1.00767187500000d3*o13*o24*o6 +
& 5.6112500000000d2*o22*o24*o6 + 2.63719531250000d3*o23*o24*o6 + 3
& .8400000000000d1*byJ*o14*o3*o6 - 2.01000000000000d2*byJ*o15*o3*o6 -
& 1.59150000000000d2*o14*o15*o3*o6 - 1.14485000000000d3*o13*o14*o1
& 5*o3*o6 + 7.6800000000000d2*o14*o16*o19*o3*o6 - 4.0200000000000d
& 3*o15*o16*o19*o3*o6 + 1.34400000000000d1*o21*o3*o6 + 1.299200000
& 00000d2*o13*o21*o3*o6 + 4.6481250000000d2*o24*o3*o6 + 2.43293750
& 000000d3*o13*o24*o3*o6 + 7.2000000000000d3*o28*o33*o6 + 9.200000
& 0000000d0*o10*o16*o18*o8 - 3.3200000000000d1*o10*o16*o18*o3*o8 -
& 5.0600000000000d1*o11*o16*o4*o8 - 4.8140000000000d2*o11*o16*o22
& *o4*o8 + 3.16000000000000d2*o11*o16*o3*o4*o8 - 6.0000000000000d0
& *o10*o16*o18*o9 + 3.00000000000000d1*o10*o16*o18*o3*o9 + 3.30000
& 00000000d1*o11*o16*o4*o9 + 4.3500000000000d2*o11*o16*o22*o4*o9 -
& 2.52000000000000d2*o11*o16*o3*o4*o9 - 8.6250000000000d1*o11*o16
& *o8*o9 - 1.55625000000000d3*o11*o16*o22*o8*o9 + 7.4250000000000d
& 2*o11*o16*o3*o8*o9
ksq(i,j) = t1 + t2 + t3 + t4
|