new parameter to be able to avoid nans at or near the puncture locations

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@44 b2a53a04-0f4f-0410-87ed-f9f25ced00cf

1 files changed, 4 insertions, 0 deletions

diff --git a/src/Equations.c b/src/Equations.c
index 1c9c5f4..4bbcbef 100644
--- a/src/Equations.c
+++ b/src/Equations.c
@@ -88,6 +88,10 @@ BY_Aijofxyz (double x, double y, double z, double Aij[3][3])
 
   r2_plus = (x - par_b) * (x - par_b) + y * y + z * z;
   r2_minus = (x + par_b) * (x + par_b) + y * y + z * z;
+  if (r2_plus < TP_Tiny)
+    r2_plus = TP_Tiny;
+  if (r2_minus < TP_Tiny)
+    r2_minus = TP_Tiny;
   r_plus = sqrt (r2_plus);
   r_minus = sqrt (r2_minus);
   r3_plus = r_plus * r2_plus;