- changed status to resolved
- removed comment
IDAnalyticBH: Smooth out singularities
Issue #610
closed
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)
-
reporter -
- edited description
- changed status to closed
- Log in to comment
Applied.