From 10d25281d7a8801f0453819c7a46588c4866355f Mon Sep 17 00:00:00 2001 From: schnetter Date: Wed, 13 Apr 2005 17:35:09 +0000 Subject: Increase the fudge factor which avoids nans. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@158 6a3ddf76-46e1-4315-99d9-bc56cac1ef84 --- src/BrillLindquist.c | 13 ++++--------- 1 file changed, 4 insertions(+), 9 deletions(-) diff --git a/src/BrillLindquist.c b/src/BrillLindquist.c index c477690..279a40e 100644 --- a/src/BrillLindquist.c +++ b/src/BrillLindquist.c @@ -44,7 +44,6 @@ void BrillLindquist(CCTK_ARGUMENTS) CCTK_REAL hole_y0[MAX_HOLES], hole_z0[MAX_HOLES]; CCTK_REAL tmp1, tmp2, tmp3; CCTK_REAL xval, yval, zval; - CCTK_REAL x_2, y_2, z_2; int i, npoints; int make_conformal_derivs = 0; @@ -130,10 +129,6 @@ void BrillLindquist(CCTK_ARGUMENTS) psi[i] = 1.0; - x_2 = SQR(x[i]); - y_2 = SQR(y[i]); - z_2 = SQR(z[i]); - xval = x[i]; yval = y[i]; zval = z[i]; @@ -145,10 +140,10 @@ void BrillLindquist(CCTK_ARGUMENTS) * ------------ */ - tmp1 = sqrt(x_2+2.0*xval*hole_x0[n]+SQR(hole_x0[n]) - +y_2+2.0*yval*hole_y0[n]+SQR(hole_y0[n]) - +z_2+2.0*zval*hole_z0[n]+SQR(hole_z0[n]) - +1.0e-20); + tmp1 = sqrt(SQR(xval+hole_x0[n]) + +SQR(yval+hole_y0[n]) + +SQR(zval+hole_z0[n]) + +1.0e-10); psi[i] += hole_mass[n]/tmp1*0.5; -- cgit v1.2.3