From 0a625978969719a333695aa42e280eeecfeb55ef Mon Sep 17 00:00:00 2001
From: eschnett <eschnett@6a3ddf76-46e1-4315-99d9-bc56cac1ef84>
Date: Thu, 27 Oct 2011 15:30:44 +0000
Subject: 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.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@183 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
---
 src/BrillLindquist.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

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;
 
-- 
cgit v1.2.3