blob: 8a6180b9c27018e3067f7e86f9f7bf9dbf835d1e (
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
|
o1 = a**2
o2 = -2.00000000000000d0*mass*rBL
o3 = rBL**2
o4 = o1 + o2 + o3
o5 = sqrt(o4)
o6 = cos(qgrd(j))
o7 = o6**2
o8 = o1*o7
o9 = o3 + o8
o10 = o9**2
o11 = 1/o10
o12 = 1.0!psi2dv(eta,q)
o13 = o12**2
o14 = 1/o13
o15 = sin(qgrd(j))
o16 = o15**2
o17 = -o1
o18 = mass**2
o19 = o17 + o18
o20 = sqrt(o19)
o21 = 1/o3
o22 = 2.50000000000000d-1*o19*o21
o23 = 1.00000000000000d0 + o22
o24 = o1 + o3
o25 = 2.00000000000000d0*o24*o3
o26 = o17 + o3
o27 = o26*o9
o28 = o25 + o27
o29 = o23**2
o30 = o10**2
o31 = 1/o30
o32 = o28**2
o33 = o13**2
o34 = o33**2
o35 = o33*o34
o36 = 1/o35
o37 = o1**2
o38 = o16**2
c exc33(i,j) = 0
c exc32(i,j) = -2.00000000000000d0*a*mass*o1*o11*o14*o15*o16*o5*o6
c & *rBL
c exc31(i,j) = 5.0000000000000d-1*a*mass*o11*o14*o16*o20*o23*o28*e
c & xp(eta)
c exc22(i,j) = 0
c exc21(i,j) = 0
c exc11(i,j) = 0
ksq(i,j) = 8.0000000000000d0*o1*o18*o3*o31*o36*o37*o38*o4*o7 + 5
& .0000000000000d-1*o1*o16*o18*o19*o29*o31*o32*o36*exp(2.000000000
& 00000d0*etagrd(i))
|