aboutsummaryrefslogtreecommitdiff
path: root/src/ckerr.x
diff options
context:
space:
mode:
Diffstat (limited to 'src/ckerr.x')
-rw-r--r--src/ckerr.x9
1 files changed, 6 insertions, 3 deletions
diff --git a/src/ckerr.x b/src/ckerr.x
index b6dbed5..7d7e754 100644
--- a/src/ckerr.x
+++ b/src/ckerr.x
@@ -1,4 +1,7 @@
- mass = sqrt(2+sqrt(4+byJ**2))
+ mass = sqrt(2.00000000000000d0 + sqrt(4.0000000000000d0 + byJ**2
+ & ))
a = byJ/mass
- rBL = sqrt(mass**2-a**2)/2*exp(eta(i,j,k))*(1+(mass+a)/(sqrt(mass**2-a**2
- & )*exp(eta(i,j,k))))*(1+(mass-a)/(sqrt(mass**2-a**2)*exp(eta(i,j,k))))
+ rBL = 5.0000000000000d-1*exp(eta(i,j,k))*(1.00000000000000d0 + (1.d0*(a
+ & + mass)*exp(-eta(i,j,k)))/sqrt(-a**2 + mass**2))*sqrt(-a**2 + mass**2)
+ & *(1.00000000000000d0 + (1.d0*(a + mass)*exp(-eta(i,j,k)))/sqrt(a**2 + m
+ & ass**2))