diff options
author | schnetter <schnetter@e296648e-0e4f-0410-bd07-d597d9acff87> | 2009-08-27 19:53:23 +0000 |
---|---|---|
committer | schnetter <schnetter@e296648e-0e4f-0410-bd07-d597d9acff87> | 2009-08-27 19:53:23 +0000 |
commit | 86ce0be086a38b55d07ec8fdfeb48402eaf66639 (patch) | |
tree | e4e32e351a12e9451f80a8115dc18d1022b6f6b8 | |
parent | 768b3046a61d7c0c00719c24db81f0a65e91f69a (diff) |
Protect sqrt argument with max, ensuring it cannot be negative.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/Exact/trunk@251 e296648e-0e4f-0410-bd07-d597d9acff87
-rw-r--r-- | src/metrics/Kerr_KerrSchild.F77 | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/src/metrics/Kerr_KerrSchild.F77 b/src/metrics/Kerr_KerrSchild.F77 index 139f271..ee1e077 100644 --- a/src/metrics/Kerr_KerrSchild.F77 +++ b/src/metrics/Kerr_KerrSchild.F77 @@ -81,7 +81,7 @@ C Spherical auxiliary coordinate r and angle theta in BH rest frame. r02 = 0.5d0 * (rho02 - a**2) $ + sqrt(0.25d0 * (rho02 - a**2)**2 + a**2 * z0**2) - r0 = sqrt(r02) + r0 = sqrt(max(0.0d0,r02)) if (Kerr_KerrSchild__parabolic .eq. 0) then C Use a power law to avoid the singularity r0 = (r0**power + eps**power)**(1.0d0/power) |