aboutsummaryrefslogtreecommitdiff
path: root/src/ksq_odd.x
blob: 3f116b934083d66a1287dfdfd266d1d4f1f962da (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
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
        o1 = mm**2
        o2 = 5.0000000000000d1*o1
        o3 = -8.0000000000000d1 + o2
        o4 = 1/o3
        o5 = 2.00000000000000d0*qgrd(j)
        o6 = cos(o5)
        o7 = 5.0000000000000d0*o6
        o8 = 3.00000000000000d0 + o7
        o9 = dnnnnngtil
        o10 = mm*phigrd(k)
        o11 = sin(o10)
        o12 = sin(qgrd(j))
        o13 = o12**2
        o14 = o13**2
        o15 = 2.50000000000000d1*o1
        o16 = -4.0000000000000d1 + o15
        o17 = 1/o16
        o18 = dnnngtil
        o19 = 2.00000000000000d1*o1
        o20 = -3.2000000000000d1 + o19
        o21 = 1/o20
        o22 = dngtil
        o23 = 5.0000000000000d0*o1
        o24 = -8.0000000000000d0 + o23
        o25 = 1/o24
        o26 = cos(o10)
        o27 = cos(qgrd(j))
        o28 = 1.03000000000000d2*o1
        o29 = -9.2800000000000d2 + o28
        o30 = o12*o13
        o31 = 3.5000000000000d1*o1
        o32 = -8.3000000000000d1 + o31
        o33 = -1.60000000000000d1 + o1
        o34 = dnnnnnnngtil
        o35 = 4.0000000000000d0*qgrd(j)
        o36 = cos(o35)
        o37 = dnnnngtil
        o38 = 1.00000000000000d1*o1
        o39 = -1.60000000000000d1 + o38
        o40 = 1/o39
        o41 = dnngtil
        o42 = dnnnnnngtil
        o43 = o27**2
        o44 = gtil
        o45 = o13*o14
        o46 = 1.07000000000000d2*o1
        o47 = 5.3000000000000d1*o1
        o48 = sin(o5)
        o49 = byJ**2
        o50 = o24**2
        o51 = 1/o50
        o52 = o26**2
        o53 = o37**2
        o54 = o1**2
        o55 = o6**2
        o56 = o36**2
        o57 = o39**2
        o58 = 1/o57
        o59 = o41**2
        o60 = o9**2
        o61 = o16**2
        o62 = 1/o61
        o63 = o18**2
        o64 = o20**2
        o65 = 1/o64
        o66 = o22**2
        o67 = o11**2
        o68 = o1*o54
        o69 = o3**2
        o70 = 1/o69
        o71 = o42**2
        o72 = o34**2
        o73 = o14**2
        o74 = o48**2
        o75 = o43**2
        o76 = o44**2
        t1 = -1.00800000000000d2*byJ*o13*o25*o26*o37 - 1.39500000000000d1*
     &  byJ*o1*o13*o25*o26*o37 - 9.7440000000000d2*byJ*o13*o25*o26*o36*o37 +
     &   1.08150000000000d2*byJ*o1*o13*o25*o26*o36*o37 + 1.33200000000000d
     &  3*byJ*o13*o26*o40*o41 - 4.6800000000000d2*byJ*o1*o13*o26*o40*o41 + 6
     &  .9720000000000d3*byJ*o13*o26*o36*o40*o41 - 2.94000000000000d3*byJ*o1
     &  *o13*o26*o36*o40*o41 - 9.6000000000000d2*byJ*o14*o26*o4*o42 + 6.00
     &  00000000000d1*byJ*o1*o14*o26*o4*o42 - 1.80000000000000d2*byJ*o14*o26
     &  *o44 + 7.2000000000000d2*byJ*o13*o26*o43*o44 - 2.70000000000000d3*
     &  byJ*o14*o26*o43*o44 + 4.5000000000000d2*byJ*o26*o44*o45 + 1.80000000
     &  000000d1*o13*o49 - 3.7296000000000d3*o13*o25*o37*o40*o41*o52 + 7
     &  .9425000000000d2*o1*o13*o25*o37*o40*o41*o52 - 5.5574400000000d4*
     &  o13*o25*o36*o37*o40*o41*o52 + 2.21991000000000d4*o1*o13*o25*o36*
     &  o37*o40*o41*o52 + 2.68800000000000d3*o14*o25*o37*o4*o42*o52 + 2.
     &  04000000000000d2*o1*o14*o25*o37*o4*o42*o52 + 2.59840000000000d4*
     &  o14*o25*o36*o37*o4*o42*o52 - 4.5080000000000d3*o1*o14*o25*o36*o3
     &  7*o4*o42*o52 - 3.5520000000000d4*o14*o4*o40*o41*o42*o52 + 1.4700
     &  0000000000d4*o1*o14*o4*o40*o41*o42*o52 - 1.85920000000000d5*o14*
     &  o36*o4*o40*o41*o42*o52 + 9.0020000000000d4*o1*o14*o36*o4*o40*o41
     &  *o42*o52
        t2 = -1.10400000000000d5*o14*o17*o18*o21*o22*o43*o52 + 2.9160000
     &  0000000d4*o1*o14*o17*o18*o21*o22*o43*o52 + 5.0400000000000d2*o14
     &  *o25*o37*o44*o52 + 6.9750000000000d1*o1*o14*o25*o37*o44*o52 + 4.
     &  8720000000000d3*o14*o25*o36*o37*o44*o52 - 5.4075000000000d2*o1*o
     &  14*o25*o36*o37*o44*o52 - 6.6600000000000d3*o14*o40*o41*o44*o52 +
     &   2.34000000000000d3*o1*o14*o40*o41*o44*o52 - 3.4860000000000d4*o
     &  14*o36*o40*o41*o44*o52 + 1.47000000000000d4*o1*o14*o36*o40*o41*o
     &  44*o52 - 2.01600000000000d3*o13*o25*o37*o43*o44*o52 - 2.79000000
     &  000000d2*o1*o13*o25*o37*o43*o44*o52 + 7.5600000000000d3*o14*o25*
     &  o37*o43*o44*o52 + 1.04625000000000d3*o1*o14*o25*o37*o43*o44*o52 
     &  - 1.94880000000000d4*o13*o25*o36*o37*o43*o44*o52 + 2.16300000000
     &  000d3*o1*o13*o25*o36*o37*o43*o44*o52 + 7.3080000000000d4*o14*o25
     &  *o36*o37*o43*o44*o52 - 8.1112500000000d3*o1*o14*o25*o36*o37*o43*
     &  o44*o52 + 2.66400000000000d4*o13*o40*o41*o43*o44*o52 - 9.3600000
     &  000000d3*o1*o13*o40*o41*o43*o44*o52 - 9.9900000000000d4*o14*o40*
     &  o41*o43*o44*o52 + 3.5100000000000d4*o1*o14*o40*o41*o43*o44*o52 +
     &   1.39440000000000d5*o13*o36*o40*o41*o43*o44*o52 - 5.880000000000
     &  0d4*o1*o13*o36*o40*o41*o43*o44*o52 - 5.2290000000000d5*o14*o36*o
     &  40*o41*o43*o44*o52 + 2.20500000000000d5*o1*o14*o36*o40*o41*o43*o
     &  44*o52 - 1.92000000000000d4*o14*o4*o42*o43*o44*o52
        t3 = 1.20000000000000d3*o1*o14*o4*o42*o43*o44*o52 - 7.6800000000
     &  000d3*o17*o21*o22*o34*o43*o45*o52 + 1.44000000000000d3*o1*o17*o2
     &  1*o22*o34*o43*o45*o52 - 1.26000000000000d3*o25*o37*o44*o45*o52 -
     &   1.74375000000000d2*o1*o25*o37*o44*o45*o52 - 1.21800000000000d4*
     &  o25*o36*o37*o44*o45*o52 + 1.35187500000000d3*o1*o25*o36*o37*o44*
     &  o45*o52 + 1.66500000000000d4*o40*o41*o44*o45*o52 - 5.85000000000
     &  00d3*o1*o40*o41*o44*o45*o52 + 8.7150000000000d4*o36*o40*o41*o44*
     &  o45*o52 - 3.6750000000000d4*o1*o36*o40*o41*o44*o45*o52 + 4.80000
     &  00000000d3*o4*o42*o44*o45*o52 - 3.00000000000000d2*o1*o4*o42*o44
     &  *o45*o52 + 7.2000000000000d4*o4*o42*o43*o44*o45*o52 - 4.50000000
     &  00000d3*o1*o4*o42*o43*o44*o45*o52 + 1.41120000000000d2*o13*o51*o
     &  52*o53 + 3.9060000000000d1*o1*o13*o51*o52*o53 + 2.72832000000000
     &  d3*o13*o36*o51*o52*o53 + 7.4760000000000d1*o1*o13*o36*o51*o52*o5
     &  3 + 1.81350000000000d2*o13*o25*o37*o40*o41*o52*o54 - 2.667000000
     &  00000d2*o13*o25*o36*o37*o40*o41*o52*o54 - 2.32500000000000d1*o14
     &  *o25*o37*o4*o42*o52*o54 + 1.80250000000000d2*o14*o25*o36*o37*o4*
     &  o42*o52*o54 - 7.8000000000000d2*o14*o4*o40*o41*o42*o52*o54 - 4.9
     &  000000000000d3*o14*o36*o4*o40*o41*o42*o52*o54 - 1.92000000000000
     &  d3*o14*o17*o18*o21*o22*o43*o52*o54 - 6.0000000000000d1*o17*o21*o
     &  22*o34*o43*o45*o52*o54
        t4 = 2.70281250000000d0*o13*o51*o52*o53*o54 - 4.1908125000000d1*
     &  o13*o36*o51*o52*o53*o54 - 2.74432000000000d4*o13*o25*o37*o40*o41
     &  *o52*o55 + 9.5432000000000d3*o1*o13*o25*o37*o40*o41*o52*o55 + 1.
     &  14688000000000d4*o14*o25*o37*o4*o42*o52*o55 - 4.2560000000000d2*
     &  o1*o14*o25*o37*o4*o42*o52*o55 - 1.20064000000000d5*o14*o4*o40*o4
     &  1*o42*o52*o55 + 5.2304000000000d4*o1*o14*o4*o40*o41*o42*o52*o55 
     &  - 1.99200000000000d6*o14*o17*o18*o21*o22*o43*o52*o55 + 2.0850000
     &  0000000d6*o1*o14*o17*o18*o21*o22*o43*o52*o55 + 1.31072000000000d
     &  3*o13*o51*o52*o53*o55 + 6.6560000000000d1*o1*o13*o51*o52*o53*o55
     &   + 2.60000000000000d2*o13*o25*o37*o40*o41*o52*o54*o55 - 1.820000
     &  00000000d1*o14*o25*o37*o4*o42*o52*o54*o55 - 2.80000000000000d3*o
     &  14*o4*o40*o41*o42*o52*o54*o55 - 5.2500000000000d5*o14*o17*o18*o2
     &  1*o22*o43*o52*o54*o55 + 8.4500000000000d-1*o13*o51*o52*o53*o54*o
     &  55 - 1.88708800000000d5*o13*o25*o37*o40*o41*o52*o56 + 1.00521050
     &  000000d5*o1*o13*o25*o37*o40*o41*o52*o56 + 1.31868800000000d4*o13
     &  *o51*o52*o53*o56 - 2.92726000000000d3*o1*o13*o51*o52*o53*o56 - 8
     &  .8322500000000d3*o13*o25*o37*o40*o41*o52*o54*o56 + 1.62450312500
     &  000d2*o13*o51*o52*o53*o54*o56 + 2.46420000000000d4*o13*o52*o58*o
     &  59 - 1.73160000000000d4*o1*o13*o52*o58*o59 + 2.57964000000000d5*
     &  o13*o36*o52*o58*o59 - 1.99416000000000d5*o1*o13*o36*o52*o58*o59
        t5 = 3.04200000000000d3*o13*o52*o54*o58*o59 + 3.8220000000000d4*
     &  o13*o36*o52*o54*o58*o59 + 1.43648000000000d5*o13*o52*o55*o58*o59
     &   - 1.07200000000000d5*o1*o13*o52*o55*o58*o59 + 2.00000000000000d
     &  4*o13*o52*o54*o55*o58*o59 + 6.7512200000000d5*o13*o52*o56*o58*o5
     &  9 - 5.6938000000000d5*o1*o13*o52*o56*o58*o59 + 1.20050000000000d
     &  5*o13*o52*o54*o56*o58*o59 - 3.07200000000000d2*byJ*o13*o25*o26*o37
     &  *o6 - 7.8000000000000d0*byJ*o1*o13*o25*o26*o37*o6 + 3.216000000000
     &  0d3*byJ*o13*o26*o40*o41*o6 - 1.20000000000000d3*byJ*o1*o13*o26*o40*o
     &  41*o6 - 1.34400000000000d3*byJ*o14*o26*o4*o42*o6 + 8.4000000000000
     &  d1*byJ*o1*o14*o26*o4*o42*o6 - 2.03712000000000d4*o13*o25*o37*o40*o
     &  41*o52*o6 + 5.8188000000000d3*o1*o13*o25*o37*o40*o41*o52*o6 - 1.
     &  46540800000000d5*o13*o25*o36*o37*o40*o41*o52*o6 + 6.571880000000
     &  0d4*o1*o13*o25*o36*o37*o40*o41*o52*o6 + 1.19552000000000d4*o14*o
     &  25*o37*o4*o42*o52*o6 - 1.84000000000000d1*o1*o14*o25*o37*o4*o42*
     &  o52*o6 + 3.6377600000000d4*o14*o25*o36*o37*o4*o42*o52*o6 - 6.311
     &  2000000000d3*o1*o14*o25*o36*o37*o4*o42*o52*o6 - 1.35488000000000
     &  d5*o14*o4*o40*o41*o42*o52*o6 + 5.7940000000000d4*o1*o14*o4*o40*o
     &  41*o42*o52*o6 - 2.60288000000000d5*o14*o36*o4*o40*o41*o42*o52*o6
     &   + 1.26028000000000d5*o1*o14*o36*o4*o40*o41*o42*o52*o6 + 9.50400
     &  00000000d5*o14*o17*o18*o21*o22*o43*o52*o6
        t6 = -6.3960000000000d5*o1*o14*o17*o18*o21*o22*o43*o52*o6 + 1.53
     &  600000000000d3*o14*o25*o37*o44*o52*o6 + 3.9000000000000d1*o1*o14
     &  *o25*o37*o44*o52*o6 - 1.60800000000000d4*o14*o40*o41*o44*o52*o6 
     &  + 6.0000000000000d3*o1*o14*o40*o41*o44*o52*o6 - 6.1440000000000d
     &  3*o13*o25*o37*o43*o44*o52*o6 - 1.56000000000000d2*o1*o13*o25*o37
     &  *o43*o44*o52*o6 + 2.30400000000000d4*o14*o25*o37*o43*o44*o52*o6 
     &  + 5.8500000000000d2*o1*o14*o25*o37*o43*o44*o52*o6 + 6.4320000000
     &  000d4*o13*o40*o41*o43*o44*o52*o6 - 2.40000000000000d4*o1*o13*o40
     &  *o41*o43*o44*o52*o6 - 2.41200000000000d5*o14*o40*o41*o43*o44*o52
     &  *o6 + 9.0000000000000d4*o1*o14*o40*o41*o43*o44*o52*o6 - 2.688000
     &  00000000d4*o14*o4*o42*o43*o44*o52*o6 + 1.68000000000000d3*o1*o14
     &  *o4*o42*o43*o44*o52*o6 + 3.8400000000000d4*o17*o21*o22*o34*o43*o
     &  45*o52*o6 - 2.64000000000000d4*o1*o17*o21*o22*o34*o43*o45*o52*o6
     &   - 3.8400000000000d3*o25*o37*o44*o45*o52*o6 - 9.7500000000000d1*
     &  o1*o25*o37*o44*o45*o52*o6 + 4.0200000000000d4*o40*o41*o44*o45*o5
     &  2*o6 - 1.50000000000000d4*o1*o40*o41*o44*o45*o52*o6 + 6.72000000
     &  00000d3*o4*o42*o44*o45*o52*o6 - 4.2000000000000d2*o1*o4*o42*o44*
     &  o45*o52*o6 + 1.00800000000000d5*o4*o42*o43*o44*o45*o52*o6 - 6.30
     &  00000000000d3*o1*o4*o42*o43*o44*o45*o52*o6 + 8.6016000000000d2*o
     &  13*o51*o52*o53*o6 + 1.40880000000000d2*o1*o13*o51*o52*o53*o6
        t7 = 8.3148800000000d3*o13*o36*o51*o52*o53*o6 - 7.1176000000000d
     &  2*o1*o13*o36*o51*o52*o53*o6 + 5.6640000000000d2*o13*o25*o37*o40*
     &  o41*o52*o54*o6 - 2.96800000000000d3*o13*o25*o36*o37*o40*o41*o52*
     &  o54*o6 - 4.5550000000000d1*o14*o25*o37*o4*o42*o52*o54*o6 + 2.523
     &  50000000000d2*o14*o25*o36*o37*o4*o42*o52*o54*o6 - 3.092000000000
     &  00d3*o14*o4*o40*o41*o42*o52*o54*o6 - 6.8600000000000d3*o14*o36*o
     &  4*o40*o41*o42*o52*o54*o6 + 6.9000000000000d4*o14*o17*o18*o21*o22
     &  *o43*o52*o54*o6 + 1.50000000000000d3*o17*o21*o22*o34*o43*o45*o52
     &  *o54*o6 + 3.02250000000000d0*o13*o51*o52*o53*o54*o6 - 2.34325000
     &  000000d1*o13*o36*o51*o52*o53*o54*o6 + 1.18992000000000d5*o13*o52
     &  *o58*o59*o6 - 8.6208000000000d4*o1*o13*o52*o58*o59*o6 + 6.228320
     &  0000000d5*o13*o36*o52*o58*o59*o6 - 4.9504000000000d5*o1*o13*o36*
     &  o52*o58*o59*o6 + 1.56000000000000d4*o13*o52*o54*o58*o59*o6 + 9.8
     &  000000000000d4*o13*o36*o52*o54*o58*o59*o6 + 6.1952000000000d2*o1
     &  4*o43*o51*o52*o60 - 1.93600000000000d2*o1*o14*o43*o51*o52*o60 + 
     &  1.51250000000000d1*o14*o43*o51*o52*o54*o60 + 4.3059200000000d3*o
     &  14*o43*o51*o52*o55*o60 - 9.5584000000000d2*o1*o14*o43*o51*o52*o5
     &  5*o60 + 5.3045000000000d1*o14*o43*o51*o52*o54*o55*o60 - 3.266560
     &  0000000d3*o14*o43*o51*o52*o6*o60 + 8.7296000000000d2*o1*o14*o43*
     &  o51*o52*o6*o60 - 5.6650000000000d1*o14*o43*o51*o52*o54*o6*o60
        t8 = 1.47200000000000d4*o18*o34*o43*o45*o52*o62 - 2.968000000000
     &  00d3*o1*o18*o34*o43*o45*o52*o62 + 1.28000000000000d2*o18*o34*o43
     &  *o45*o52*o54*o62 - 5.3120000000000d4*o18*o34*o43*o45*o52*o6*o62 
     &  + 2.57200000000000d4*o1*o18*o34*o43*o45*o52*o6*o62 - 1.400000000
     &  00000d3*o18*o34*o43*o45*o52*o54*o6*o62 + 1.05800000000000d5*o14*
     &  o43*o52*o62*o63 - 2.94400000000000d4*o1*o14*o43*o52*o62*o63 + 2.
     &  04800000000000d3*o14*o43*o52*o54*o62*o63 + 1.37780000000000d6*o1
     &  4*o43*o52*o55*o62*o63 - 1.16200000000000d6*o1*o14*o43*o52*o55*o6
     &  2*o63 + 2.45000000000000d5*o14*o43*o52*o54*o55*o62*o63 - 7.63600
     &  00000000d5*o14*o43*o52*o6*o62*o63 + 4.2824000000000d5*o1*o14*o43
     &  *o52*o6*o62*o63 - 4.4800000000000d4*o14*o43*o52*o54*o6*o62*o63 +
     &   2.88000000000000d4*o14*o43*o52*o65*o66 - 7.2000000000000d3*o1*o
     &  14*o43*o52*o65*o66 + 4.5000000000000d2*o14*o43*o52*o54*o65*o66 +
     &   7.2000000000000d5*o14*o43*o52*o55*o65*o66 - 9.0000000000000d5*o
     &  1*o14*o43*o52*o55*o65*o66 + 2.81250000000000d5*o14*o43*o52*o54*o
     &  55*o65*o66 - 2.88000000000000d5*o14*o43*o52*o6*o65*o66 + 2.16000
     &  000000000d5*o1*o14*o43*o52*o6*o65*o66 - 2.25000000000000d4*o14*o
     &  43*o52*o54*o6*o65*o66 - 2.10600000000000d4*o1*o14*o17*o18*o21*o2
     &  2*o67 - 6.7562000000000d4*o1*o14*o17*o18*o22*o25*o67 - 7.0180000
     &  000000d4*o1*o14*o17*o18*o22*o40*o67 + 1.92000000000000d3*o1*o14*
     &  o17*o42*o43*o44*o67
        t9 = -3.3600000000000d3*o1*o17*o42*o43*o44*o45*o67 - 5.440000000
     &  0000d2*o1*o17*o25*o27*o30*o37*o42*o48*o67 + 1.98400000000000d3*o
     &  1*o17*o27*o30*o40*o41*o42*o48*o67 - 1.02000000000000d3*o1*o12*o2
     &  5*o27*o37*o44*o48*o67 + 1.78500000000000d3*o1*o25*o27*o30*o37*o4
     &  4*o48*o67 + 3.7200000000000d3*o1*o12*o27*o40*o41*o44*o48*o67 - 6
     &  .5100000000000d3*o1*o27*o30*o40*o41*o44*o48*o67 + 4.625200000000
     &  0d4*o14*o17*o18*o22*o25*o54*o67 + 6.6494000000000d4*o14*o17*o18*
     &  o22*o40*o54*o67 - 1.20000000000000d2*o14*o17*o42*o43*o44*o54*o67
     &   + 2.10000000000000d2*o17*o42*o43*o44*o45*o54*o67 + 1.2200000000
     &  0000d2*o17*o25*o27*o30*o37*o42*o48*o54*o67 - 5.7200000000000d2*o
     &  17*o27*o30*o40*o41*o42*o48*o54*o67 + 1.65000000000000d2*o12*o25*
     &  o27*o37*o44*o48*o54*o67 - 2.88750000000000d2*o25*o27*o30*o37*o44
     &  *o48*o54*o67 - 8.4000000000000d2*o12*o27*o40*o41*o44*o48*o54*o67
     &   + 1.47000000000000d3*o27*o30*o40*o41*o44*o48*o54*o67 - 5.850000
     &  0000000d4*o1*o14*o17*o18*o21*o22*o55*o67 - 1.42506000000000d5*o1
     &  *o14*o17*o18*o22*o25*o55*o67 - 1.31652000000000d5*o1*o14*o17*o18
     &  *o22*o40*o55*o67 + 9.2016000000000d4*o14*o17*o18*o22*o25*o54*o55
     &  *o67 + 1.25082000000000d5*o14*o17*o18*o22*o40*o54*o55*o67 - 7.02
     &  00000000000d4*o1*o14*o17*o18*o21*o22*o6*o67 - 1.96332000000000d5
     &  *o1*o14*o17*o18*o22*o25*o6*o67 - 1.92408000000000d5*o1*o14*o17*o
     &  18*o22*o40*o6*o67 - 1.29920000000000d3*o1*o17*o25*o27*o30*o37*o4
     &  2*o48*o6*o67 + 1.33760000000000d4*o1*o17*o27*o30*o40*o41*o42*o48
     &  *o6*o67
        t10 = -2.43600000000000d3*o1*o12*o25*o27*o37*o44*o48*o6*o67 + 4.
     &  2630000000000d3*o1*o25*o27*o30*o37*o44*o48*o6*o67 + 2.5080000000
     &  0000d4*o1*o12*o27*o40*o41*o44*o48*o6*o67 - 4.3890000000000d4*o1*
     &  o27*o30*o40*o41*o44*o48*o6*o67 + 1.30692000000000d5*o14*o17*o18*
     &  o22*o25*o54*o6*o67 + 1.82664000000000d5*o14*o17*o18*o22*o40*o54*
     &  o6*o67 + 1.08400000000000d2*o17*o25*o27*o30*o37*o42*o48*o54*o6*o
     &  67 - 6.5320000000000d3*o17*o27*o30*o40*o41*o42*o48*o54*o6*o67 + 
     &  5.1000000000000d1*o12*o25*o27*o37*o44*o48*o54*o6*o67 - 8.9250000
     &  000000d1*o25*o27*o30*o37*o44*o48*o54*o6*o67 - 1.06800000000000d4
     &  *o12*o27*o40*o41*o44*o48*o54*o6*o67 + 1.86900000000000d4*o27*o30
     &  *o40*o41*o44*o48*o54*o6*o67 + 4.0960000000000d3*o1*o14*o60*o62*o
     &  67 - 5.1200000000000d2*o14*o54*o60*o62*o67 + 9.2160000000000d3*o
     &  1*o14*o55*o60*o62*o67 - 1.15200000000000d3*o14*o54*o55*o60*o62*o
     &  67 + 1.22880000000000d4*o1*o14*o6*o60*o62*o67 - 1.53600000000000
     &  d3*o14*o54*o6*o60*o62*o67 + 2.63438000000000d5*o1*o14*o62*o63*o6
     &  7 - 1.49158000000000d5*o14*o54*o62*o63*o67 + 5.8030200000000d5*o
     &  1*o14*o55*o62*o63*o67 - 3.2531400000000d5*o14*o54*o55*o62*o63*o6
     &  7 + 7.8147600000000d5*o1*o14*o6*o62*o63*o67 - 4.4056800000000d5*
     &  o14*o54*o6*o62*o63*o67 + 6.8890000000000d3*o1*o14*o51*o66*o67 - 
     &  5.8100000000000d3*o14*o51*o54*o66*o67 + 1.36890000000000d4*o1*o1
     &  4*o51*o55*o66*o67
        t11 = -1.05300000000000d4*o14*o51*o54*o55*o66*o67 + 1.4641000000
     &  0000d4*o1*o14*o58*o66*o67 - 1.69400000000000d4*o14*o54*o58*o66*o
     &  67 + 2.52810000000000d4*o1*o14*o55*o58*o66*o67 - 2.8620000000000
     &  0d4*o14*o54*o55*o58*o66*o67 + 1.94220000000000d4*o1*o14*o51*o6*o
     &  66*o67 - 1.56600000000000d4*o14*o51*o54*o6*o66*o67 + 3.847800000
     &  0000d4*o1*o14*o58*o6*o66*o67 - 4.4040000000000d4*o14*o54*o58*o6*
     &  o66*o67 + 8.1000000000000d3*o1*o14*o65*o66*o67 + 2.2500000000000
     &  0d4*o1*o14*o55*o65*o66*o67 + 2.70000000000000d4*o1*o14*o6*o65*o6
     &  6*o67 - 7.4900000000000d3*o14*o17*o18*o22*o25*o67*o68 - 1.498000
     &  00000000d4*o14*o17*o18*o22*o40*o67*o68 - 5.5000000000000d0*o17*o
     &  25*o27*o30*o37*o42*o48*o67*o68 + 2.80000000000000d1*o17*o27*o30*
     &  o40*o41*o42*o48*o67*o68 - 1.43100000000000d4*o14*o17*o18*o22*o25
     &  *o55*o67*o68 - 2.86200000000000d4*o14*o17*o18*o22*o40*o55*o67*o6
     &  8 - 2.07600000000000d4*o14*o17*o18*o22*o25*o6*o67*o68 - 4.152000
     &  0000000d4*o14*o17*o18*o22*o40*o6*o67*o68 - 1.70000000000000d0*o1
     &  7*o25*o27*o30*o37*o42*o48*o6*o67*o68 + 3.5600000000000d2*o17*o27
     &  *o30*o40*o41*o42*o48*o6*o67*o68 + 1.60000000000000d1*o14*o60*o62
     &  *o67*o68 + 3.6000000000000d1*o14*o55*o60*o62*o67*o68 + 4.8000000
     &  000000d1*o14*o6*o60*o62*o67*o68 + 2.28980000000000d4*o14*o62*o63
     &  *o67*o68 + 5.0562000000000d4*o14*o55*o62*o63*o67*o68
        t12 = 6.8052000000000d4*o14*o6*o62*o63*o67*o68 + 1.2250000000000
     &  0d3*o14*o51*o66*o67*o68 + 2.02500000000000d3*o14*o51*o55*o66*o67
     &  *o68 + 4.9000000000000d3*o14*o58*o66*o67*o68 + 8.1000000000000d3
     &  *o14*o55*o58*o66*o67*o68 + 3.15000000000000d3*o14*o51*o6*o66*o67
     &  *o68 + 1.26000000000000d4*o14*o58*o6*o66*o67*o68 + 1.42420000000
     &  000d4*o1*o14*o60*o67*o70 - 1.90400000000000d3*o14*o54*o60*o67*o7
     &  0 + 3.15540000000000d4*o1*o14*o55*o60*o67*o70 - 4.2480000000000d
     &  3*o14*o54*o55*o60*o67*o70 + 4.2396000000000d4*o1*o14*o6*o60*o67*
     &  o70 - 5.6880000000000d3*o14*o54*o6*o60*o67*o70 + 6.4000000000000
     &  d1*o14*o60*o67*o68*o70 + 1.44000000000000d2*o14*o55*o60*o67*o68*
     &  o70 + 1.92000000000000d2*o14*o6*o60*o67*o68*o70 + 5.120000000000
     &  0d2*o1*o43*o45*o62*o67*o71 - 6.4000000000000d1*o43*o45*o54*o62*o
     &  67*o71 + 2.00000000000000d0*o43*o45*o62*o67*o68*o71 + 1.28000000
     &  000000d4*o45*o52*o70*o71 - 1.60000000000000d3*o1*o45*o52*o70*o71
     &   + 5.0000000000000d1*o45*o52*o54*o70*o71 + 2.50880000000000d4*o4
     &  5*o52*o55*o70*o71 - 3.13600000000000d3*o1*o45*o52*o55*o70*o71 + 
     &  9.8000000000000d1*o45*o52*o54*o55*o70*o71 + 3.5840000000000d4*o4
     &  5*o52*o6*o70*o71 - 4.4800000000000d3*o1*o45*o52*o6*o70*o71 + 1.4
     &  0000000000000d2*o45*o52*o54*o6*o70*o71
        t13 = -1.20000000000000d4*o4*o42*o44*o52*o73 + 7.5000000000000d2
     &  *o1*o4*o42*o44*o52*o73 - 1.68000000000000d4*o4*o42*o44*o52*o6*o7
     &  3 + 1.05000000000000d3*o1*o4*o42*o44*o52*o6*o73 + 5.120000000000
     &  0d2*o43*o52*o62*o72*o73 - 6.4000000000000d1*o1*o43*o52*o62*o72*o
     &  73 + 2.00000000000000d0*o43*o52*o54*o62*o72*o73 - 1.054000000000
     &  00d3*o1*o25*o37*o40*o41*o67*o74 + 1.44500000000000d2*o1*o51*o53*
     &  o67*o74 + 4.0850000000000d2*o25*o37*o40*o41*o54*o67*o74 - 4.6750
     &  000000000d1*o51*o53*o54*o67*o74 - 1.69708000000000d4*o1*o25*o37*
     &  o40*o41*o55*o67*o74 + 8.2418000000000d2*o1*o51*o53*o55*o67*o74 +
     &   7.5821000000000d3*o25*o37*o40*o41*o54*o55*o67*o74 - 3.451000000
     &  0000d1*o51*o53*o54*o55*o67*o74 + 1.92200000000000d3*o1*o58*o59*o
     &  67*o74 - 8.6800000000000d2*o54*o58*o59*o67*o74 + 8.7362000000000
     &  d4*o1*o55*o58*o59*o67*o74 - 7.4404000000000d4*o54*o55*o58*o59*o6
     &  7*o74 - 9.6232000000000d3*o1*o25*o37*o40*o41*o6*o67*o74 + 6.9020
     &  000000000d2*o1*o51*o53*o6*o67*o74 + 4.7966000000000d3*o25*o37*o4
     &  0*o41*o54*o6*o67*o74 - 1.26100000000000d2*o51*o53*o54*o6*o67*o74
     &   + 2.59160000000000d4*o1*o58*o59*o6*o67*o74 - 1.68880000000000d4
     &  *o54*o58*o59*o6*o67*o74 - 3.8500000000000d1*o25*o37*o40*o41*o67*
     &  o68*o74 + 3.7812500000000d0*o51*o53*o67*o68*o74
        t14 = -1.51300000000000d2*o25*o37*o40*o41*o55*o67*o68*o74 + 3.61
     &  25000000000d-1*o51*o53*o55*o67*o68*o74 + 9.8000000000000d1*o58*o
     &  59*o67*o68*o74 + 1.58420000000000d4*o55*o58*o59*o67*o68*o74 - 5.
     &  0140000000000d2*o25*o37*o40*o41*o6*o67*o68*o74 + 2.3375000000000
     &  0d0*o51*o53*o6*o67*o68*o74 + 2.49200000000000d3*o58*o59*o6*o67*o
     &  68*o74 - 3.6000000000000d3*o14*o43*o52*o76 + 4.5000000000000d2*o
     &  45*o52*o76 + 2.25000000000000d4*o43*o45*o52*o76 + 1.800000000000
     &  00d3*o1*o13*o43*o67*o76 - 6.3000000000000d3*o1*o14*o43*o67*o76 +
     &   5.5125000000000d3*o1*o43*o45*o67*o76 - 2.25000000000000d3*o52*o
     &  73*o76 + 2.81250000000000d3*o13*o52*o73*o76 - 3.3750000000000d4*
     &  o43*o52*o73*o76 + 7.2000000000000d3*o13*o52*o75*o76 - 5.40000000
     &  00000d4*o14*o52*o75*o76 + 1.01250000000000d5*o45*o52*o75*o76 - 1
     &  .61920000000000d4*o14*o17*o18*o25*o43*o52*o9 + 4.7828000000000d3
     &  *o1*o14*o17*o18*o25*o43*o52*o9 + 8.4480000000000d3*o14*o21*o22*o
     &  25*o43*o52*o9 - 2.37600000000000d3*o1*o14*o21*o22*o25*o43*o52*o9
     &   - 1.12640000000000d3*o17*o25*o34*o43*o45*o52*o9 + 2.46400000000
     &  000d2*o1*o17*o25*o34*o43*o45*o52*o9 - 3.5200000000000d2*o14*o17*
     &  o18*o25*o43*o52*o54*o9 + 1.65000000000000d2*o14*o21*o22*o25*o43*
     &  o52*o54*o9
        t15 = -1.10000000000000d1*o17*o25*o34*o43*o45*o52*o54*o9 - 1.540
     &  48000000000d5*o14*o17*o18*o25*o43*o52*o55*o9 + 8.2058000000000d4
     &  *o1*o14*o17*o18*o25*o43*o52*o55*o9 + 1.11360000000000d5*o14*o21*
     &  o22*o25*o43*o52*o55*o9 - 8.1960000000000d4*o1*o14*o21*o22*o25*o4
     &  3*o52*o55*o9 - 7.2100000000000d3*o14*o17*o18*o25*o43*o52*o54*o55
     &  *o9 + 7.7250000000000d3*o14*o21*o22*o25*o43*o52*o54*o55*o9 + 1.0
     &  1120000000000d5*o14*o17*o18*o25*o43*o52*o6*o9 - 4.4447200000000d
     &  4*o1*o14*o17*o18*o25*o43*o52*o6*o9 - 6.4512000000000d4*o14*o21*o
     &  22*o25*o43*o52*o6*o9 + 3.8256000000000d4*o1*o14*o21*o22*o25*o43*
     &  o52*o6*o9 + 2.96960000000000d3*o17*o25*o34*o43*o45*o52*o6*o9 - 5
     &  .1520000000000d2*o1*o17*o25*o34*o43*o45*o52*o6*o9 + 4.5092000000
     &  000d3*o14*o17*o18*o25*o43*o52*o54*o6*o9 - 4.4340000000000d3*o14*
     &  o21*o22*o25*o43*o52*o54*o6*o9 + 2.06000000000000d1*o17*o25*o34*o
     &  43*o45*o52*o54*o6*o9 + 1.06240000000000d4*o1*o14*o17*o22*o25*o67
     &  *o9 - 7.1126000000000d4*o1*o14*o17*o18*o4*o67*o9 + 1.62000000000
     &  000d3*o1*o14*o21*o22*o4*o67*o9 + 2.87980000000000d4*o1*o14*o22*o
     &  4*o40*o67*o9 - 5.1440000000000d3*o14*o17*o22*o25*o54*o67*o9 + 3.
     &  01060000000000d4*o14*o17*o18*o4*o54*o67*o9 - 1.85960000000000d4*
     &  o14*o22*o4*o40*o54*o67*o9 + 2.24640000000000d4*o1*o14*o17*o22*o2
     &  5*o55*o67*o9 - 1.52406000000000d5*o1*o14*o17*o18*o4*o55*o67*o9 +
     &   4.5000000000000d3*o1*o14*o21*o22*o4*o55*o67*o9 + 5.628600000000
     &  0d4*o1*o14*o22*o4*o40*o55*o67*o9
        t16 = -1.00440000000000d4*o14*o17*o22*o25*o54*o55*o67*o9 + 6.622
     &  2000000000d4*o14*o17*o18*o4*o54*o55*o67*o9 - 3.5676000000000d4*o
     &  14*o22*o4*o40*o54*o55*o67*o9 + 3.09120000000000d4*o1*o14*o17*o22
     &  *o25*o6*o67*o9 - 2.08212000000000d5*o1*o14*o17*o18*o4*o6*o67*o9 
     &  + 5.4000000000000d3*o1*o14*o21*o22*o4*o6*o67*o9 + 8.067600000000
     &  0d4*o1*o14*o22*o4*o40*o6*o67*o9 - 1.44120000000000d4*o14*o17*o22
     &  *o25*o54*o6*o67*o9 + 8.9304000000000d4*o14*o17*o18*o4*o54*o6*o67
     &  *o9 - 5.1648000000000d4*o14*o22*o4*o40*o54*o6*o67*o9 - 5.2096000
     &  000000d4*o1*o14*o18*o62*o67*o9 + 1.69520000000000d4*o14*o18*o54*
     &  o62*o67*o9 - 1.16928000000000d5*o1*o14*o18*o55*o62*o67*o9 + 3.78
     &  36000000000d4*o14*o18*o54*o55*o62*o67*o9 - 1.56096000000000d5*o1
     &  *o14*o18*o6*o62*o67*o9 + 5.0652000000000d4*o14*o18*o54*o6*o62*o6
     &  7*o9 + 2.80000000000000d2*o14*o17*o22*o25*o67*o68*o9 - 1.7120000
     &  0000000d3*o14*o17*o18*o4*o67*o68*o9 + 1.12000000000000d3*o14*o22
     &  *o4*o40*o67*o68*o9 + 5.4000000000000d2*o14*o17*o22*o25*o55*o67*o
     &  68*o9 - 3.8160000000000d3*o14*o17*o18*o4*o55*o67*o68*o9 + 2.1600
     &  0000000000d3*o14*o22*o4*o40*o55*o67*o68*o9 + 7.8000000000000d2*o
     &  14*o17*o22*o25*o6*o67*o68*o9 - 5.1120000000000d3*o14*o17*o18*o4*
     &  o6*o67*o68*o9 + 3.12000000000000d3*o14*o22*o4*o40*o6*o67*o68*o9 
     &  - 8.5600000000000d2*o14*o18*o62*o67*o68*o9 - 1.90800000000000d3*
     &  o14*o18*o55*o62*o67*o68*o9 - 2.55600000000000d3*o14*o18*o6*o62*o
     &  67*o68*o9
        ksq(i,j,k) = t1 + t2 + t3 + t4 + t5 + t6 + t7 + t8 + t9 + t10 + 
     &  t11 + t12 + t13 + t14 + t15 + t16