diff options
| author | schnetter <schnetter@6a3ddf76-46e1-4315-99d9-bc56cac1ef84> | 2004-06-23 07:13:51 +0000 |
|---|---|---|
| committer | schnetter <schnetter@6a3ddf76-46e1-4315-99d9-bc56cac1ef84> | 2004-06-23 07:13:51 +0000 |
| commit | 8e9335ee54ecc33e8b07be4ce6fd815b316f768c (patch) | |
| tree | 9489129c55587d8d379000c6f4b798ef76b7358d | |
| parent | e340f99ace935245fb42ee42a8e2225bce351d65 (diff) | |
Do not produce nans on the origin.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@151 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
| -rw-r--r-- | src/BrillLindquist.c | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/src/BrillLindquist.c b/src/BrillLindquist.c index e0ce8ac..c477690 100644 --- a/src/BrillLindquist.c +++ b/src/BrillLindquist.c @@ -147,7 +147,8 @@ 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])); + +z_2+2.0*zval*hole_z0[n]+SQR(hole_z0[n]) + +1.0e-20); psi[i] += hole_mass[n]/tmp1*0.5; |