IDAnalyticBH: Smooth out singularities

[Erik Schnetter]

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:

[Erik Schnetter]

Applied.