blob: d351f25d2170d0993e5d5b6db7fcdd57d5c94089 (
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
|
o1 = -xi
o2 = cos(qgrd(j))
o3 = xi*o2
o4 = 1.d0 + o1 + o3
o5 = -eta0_me
o6 = etagrd(i) + o5
o7 = o6**2
o8 = sigma_me**2
o9 = o8**2
o10 = 1/o9
o11 = 1/o8
o12 = -(o11*o7)
o13 = exp(o12)
o14 = eta0_me + etagrd(i)
o15 = o14**2
o16 = -(o11*o15)
o17 = exp(o16)
o18 = sin(qgrd(j))
o19 = o18**n_me
o20 = o2**2
o21 = o13 + o17
qfetaeta(i,j) = amp_me*o19*o4*(-2.d0*o11*o13 - 2.d0*o11*o17 + 4.d0*o10*o15*o
& 17 + 4.d0*o10*o13*o7)
qfqq(i,j) = -(amp_me*xi*n_me*o19*o2*o21) - amp_me*xi*(1.d0 + n_me)*o19*o2*o21 -
& amp_me*n_me*o19*o21*o4 + amp_me*(-1.d0 + n_me)*n_me*o18**(-2 + n)*o20*o21*o
& 4
|