User: eschnett Date: 2011/10/27 10:30 AM
Modified: /trunk/src/ BrillLindquist.c
Log: Use "epsilon" also for Brill-Lindquist
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.
File Changes:
Directory: /trunk/src/ ======================
File [modified]: BrillLindquist.c Delta lines: +1 -1 =================================================================== --- trunk/src/BrillLindquist.c 2011-10-25 18:45:16 UTC (rev 182) +++ trunk/src/BrillLindquist.c 2011-10-27 15:30:44 UTC (rev 183) @@ -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;