diff options
Diffstat (limited to 'src/ksq_sergio.x')
-rw-r--r-- | src/ksq_sergio.x | 46 |
1 files changed, 19 insertions, 27 deletions
diff --git a/src/ksq_sergio.x b/src/ksq_sergio.x index 8a6180b..76a947f 100644 --- a/src/ksq_sergio.x +++ b/src/ksq_sergio.x @@ -14,36 +14,28 @@ 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 + o17 = o1 + o3 + o18 = 2.00000000000000d0*o17*o3 + o19 = -o1 + o20 = o19 + o3 + o21 = o20*o9 + o22 = o18 + o21 + o23 = mass**2 + o24 = o10**2 + o25 = 1/o24 + o26 = o22**2 + o27 = o13**2 + o28 = o27**2 + o29 = o27*o28 + o30 = 1/o29 + o31 = o1**2 + o32 = 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 exc31(i,j) = a*mass*o11*o14*o16*o22 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)) + ksq(i,j) = 2.00000000000000d0*o1*o16*o23*o25*o26*o30 + 8.0000000 + & 000000d0*o1*o23*o25*o3*o30*o31*o32*o4*o7 |