diff options
Diffstat (limited to 'src/qfunc_even.x')
| -rw-r--r-- | src/qfunc_even.x | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/src/qfunc_even.x b/src/qfunc_even.x index d351f25..00f9872 100644 --- a/src/qfunc_even.x +++ b/src/qfunc_even.x @@ -1,7 +1,7 @@ o1 = -xi o2 = cos(qgrd(j)) o3 = xi*o2 - o4 = 1.d0 + o1 + o3 + o4 = 1.0 + o1 + o3 o5 = -eta0_me o6 = etagrd(i) + o5 o7 = o6**2 @@ -19,8 +19,8 @@ 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 + qfetaeta(i,j) = amp_me*o19*o4*(-2.0*o11*o13 - 2.0*o11*o17 + 4.0*o10*o15*o + & 17 + 4.0*o10*o13*o7) + qfqq(i,j) = -(amp_me*xi*n_me*o19*o2*o21) - amp_me*xi*(1.0 + n_me)*o19*o2*o21 - + & amp_me*n_me*o19*o21*o4 + amp_me*(-1.0 + n_me)*n_me*o18**(-2 + n)*o20*o21*o & 4 |