diff options
| author | schnetter <schnetter@6a3ddf76-46e1-4315-99d9-bc56cac1ef84> | 2005-04-13 17:35:09 +0000 |
|---|---|---|
| committer | schnetter <schnetter@6a3ddf76-46e1-4315-99d9-bc56cac1ef84> | 2005-04-13 17:35:09 +0000 |
| commit | 10d25281d7a8801f0453819c7a46588c4866355f (patch) | |
| tree | 9b8d278adda64a5de746cc27992ba18d67bd14c1 | |
| parent | 94792b144fa5bf6359a4cf3755d96ff150fa315d (diff) | |
Increase the fudge factor which avoids nans.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@158 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
| -rw-r--r-- | src/BrillLindquist.c | 13 |
1 files 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; |