blob: e91a8e0080ef5d5f26433df7ca7a0a69a8924eb7 (
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
|
o1 = cos(phi)
o2 = cos(q)
o3 = o2**2
o4 = -3.00000000000000d0*eta
o5 = exp(o4)
o6 = sin(phi)
o7 = sin(q)
o8 = o7**2
o9 = dnnnnngtil
o10 = -3.6000000000000d1*o9
o11 = dnnnnnnngtil
o12 = -4.0000000000000d0*o11*o8
o13 = 1.16000000000000d2*o8*o9
o14 = dnnngtil
o15 = -4.1500000000000d2*o8
o16 = 1.50000000000000d2 + o15
o17 = o14*o16
o18 = dngtil
o19 = 5.0000000000000d0*o8
o20 = -2.00000000000000d0 + o19
o21 = 7.5000000000000d1*o18*o20
o22 = o10 + o12 + o13 + o17 + o21
o23 = dnnnnnngtil
o24 = 3.00000000000000d0*o23
o25 = 2.00000000000000d0*q
o26 = cos(o25)
o27 = 4.0000000000000d0*o23*o26
o28 = 4.0000000000000d0*q
o29 = cos(o28)
o30 = -7.0000000000000d0*o23*o29
o31 = dnnnngtil
o32 = 2.10000000000000d1*o31
o33 = 6.4000000000000d1*o26*o31
o34 = 2.03000000000000d2*o29*o31
o35 = o24 + o27 + o30 + o32 + o33 + o34
o36 = -4.0000000000000d0*o35
o37 = 2.68000000000000d2*o26
o38 = 5.8100000000000d2*o29
o39 = 1.11000000000000d2 + o37 + o38
o40 = dnngtil
o41 = 5.0000000000000d0*o39*o40
o42 = 2.00000000000000d1*o26
o43 = 3.5000000000000d1*o29
o44 = 9.0000000000000d0 + o42 + o43
o45 = gtil
o46 = -7.5000000000000d1*o44*o45
o47 = o36 + o41 + o46
o48 = o1**2
o49 = o6**2
o50 = o7*o8
hxx(i,j,k) = 6.6666666666667d-1*o1*o22*o3*o5*o6*o8 - 8.333333333
& 3333d-2*o1*o47*o5*o6*o8
hxy(i,j,k) = -3.3333333333333d-1*o22*o3*o48*o5*o8 + 4.1666666666
& 667d-2*o47*o48*o5*o8 + 3.3333333333333d-1*o22*o3*o49*o5*o8 - 4.1
& 666666666667d-2*o47*o49*o5*o8
hxz(i,j,k) = -3.3333333333333d-1*o2*o22*o5*o50*o6 - 4.1666666666
& 667d-2*o2*o47*o5*o6*o7
hyy(i,j,k) = -6.6666666666667d-1*o1*o22*o3*o5*o6*o8 + 8.33333333
& 33333d-2*o1*o47*o5*o6*o8
hyz(i,j,k) = 3.3333333333333d-1*o1*o2*o22*o5*o50 + 4.16666666666
& 67d-2*o1*o2*o47*o5*o7
hzz(i,j,k) = 0
|