diff options
author | allen <allen@0a4070d5-58f5-498f-b6c0-2693e757fa0f> | 2000-01-24 08:31:55 +0000 |
---|---|---|
committer | allen <allen@0a4070d5-58f5-498f-b6c0-2693e757fa0f> | 2000-01-24 08:31:55 +0000 |
commit | 1ce7be9da747a06a27ed8c880c4be740d27c5331 (patch) | |
tree | f40ff221e8f51315ff8e0a0a7b1d7d42cafb52f2 /src/IDAxiBrillBH.F | |
parent | cd8420110973d389f85d83bc0fff1b4dab495d95 (diff) |
Fix from Thomas Radke for floating point exception on the Alphas caused
by exponention by a real number .... here's Thomas's description:
The include file generated by Mathematica contains an exponentiation
operation
with a REAL exponent.
My Fortran book says that
y**x = exp(x * ln(y))
On O2K this intrinsic function checks the exponent against 0,
on Alphas it doesn't and you get a FP exception if y <= 0.
I just changed the type of the exponent to be INTEGER instead of REAL
which was probably also the intention of the programmer.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@8 0a4070d5-58f5-498f-b6c0-2693e757fa0f
Diffstat (limited to 'src/IDAxiBrillBH.F')
-rw-r--r-- | src/IDAxiBrillBH.F | 1 |
1 files changed, 1 insertions, 0 deletions
diff --git a/src/IDAxiBrillBH.F b/src/IDAxiBrillBH.F index d064efc..01000ef 100644 --- a/src/IDAxiBrillBH.F +++ b/src/IDAxiBrillBH.F @@ -62,6 +62,7 @@ c Perhaps this and others should go into cctk.h real*8 o70,o71,o72,o73,o74,o75,o76,o77,o78,o79 real*8 o80,o81,o82,o83,o84,o85,o86,o87,o88,o89 real*8 o90,o91,o92,o93,o94,o95,o96,o97,o98,o99 + integer i22 real*8 pi real*8 adm integer :: nx,ny,nz |