blob: 39a175f7af4b50ac789aa932c1b0e36f9e2fcf0c (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
o1 = cos(q)
o2 = 1.
o3 = o2**2
o4 = 1/o3
o5 = sin(q)
o6 = 2.00000000000000d0 + 3.0d0
o7 = o5**2
o8 = -(o6*o7)
o9 = 1.00000000000000d0 + 3.0d0 + o8
o10 = o3**2
o11 = o10**2
o12 = dngtil**2
o13 = o1**2
o14 = gtil**2
o15 = o9**2
c exc33(i,j,k) = 0
c exc32(i,j,k) = -(dngtil*o1*o4*o5**n)
c exc31(i,j,k) = gtil*o4*o5**(-1.00000000000000d0 + n)*o9
c exc22(i,j,k) = 0
c exc21(i,j,k) = 0
c exc11(i,j,k) = 0
ksq(i,j,k) = (2.00000000000000d0*o5**(2.00000000000000d0*(-2.000
& 00000000000d0 + 3.d0))*(o14*o15 + o12*o13*o7))/(o10*o11)
|