Use "epsilon" also for Brill-LindquistHEADmaster

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.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@183 6a3ddf76-46e1-4315-99d9-bc56cac1ef84

1 files changed, 1 insertions, 1 deletions

diff --git a/src/BrillLindquist.c b/src/BrillLindquist.c
index eef51b5..0323ab4 100644
--- a/src/BrillLindquist.c
+++ b/src/BrillLindquist.c
@@ -143,7 +143,7 @@ void BrillLindquist(CCTK_ARGUMENTS)
       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;