IDAnalyticBH: Smooth out singularities

IDAnalyticBH has an "epsilon" parameter to smooth out singularites. This parameter should also be used in Brill-Lindquist initial data. Otherwise, moving puncture evolution of these data is not possible without an explicit smoothing step.

Index: src/BrillLindquist.c
===================================================================
--- src/BrillLindquist.c    (revision 180)
+++ src/BrillLindquist.c    (working copy)
@@ -143,7 +143,7 @@
       tmp1 = sqrt(SQR(xval+hole_x0[n])
                  +SQR(yval+hole_y0[n])
                  +SQR(zval+hole_z0[n])
-                 +1.0e-20);
+                 +SQR(epsilon));

       psi[i] += hole_mass[n]/tmp1*0.5;

Keyword:

Comments (2)