#610: IDAnalyticBH: Smooth out singularities -----------------------------------+---------------------------------------- Reporter: eschnett | Owner: Type: defect | Status: new Priority: major | Milestone: Component: EinsteinToolkit thorn | Version: Keywords: | -----------------------------------+---------------------------------------- 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;
}}}