aboutsummaryrefslogtreecommitdiff
path: root/src/ckerr.x
diff options
context:
space:
mode:
authorryoji <ryoji@b6f3ac56-194f-0410-8878-cdf6079d7f1b>2001-09-11 05:47:49 +0000
committerryoji <ryoji@b6f3ac56-194f-0410-8878-cdf6079d7f1b>2001-09-11 05:47:49 +0000
commit5dae8919dc27e53baf338422ee1906a7389fd2b6 (patch)
tree3870f164c28aba522bdc637d94bb1caea0d244af /src/ckerr.x
parent266e41ac2faa66e868afb82d510f58fe551c4d81 (diff)
sergio ivp
bugs are fixed git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiOddBrillBH/trunk@35 b6f3ac56-194f-0410-8878-cdf6079d7f1b
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))