diff options
author | eschnett <eschnett@6a3ddf76-46e1-4315-99d9-bc56cac1ef84> | 2011-10-27 15:30:44 +0000 |
---|---|---|
committer | eschnett <eschnett@6a3ddf76-46e1-4315-99d9-bc56cac1ef84> | 2011-10-27 15:30:44 +0000 |
commit | 0a625978969719a333695aa42e280eeecfeb55ef (patch) | |
tree | a6639f89e2aed54b5c8a6ea2b0695d834c422cd3 | |
parent | f300f50f897ae891c414a829b47e77f9ddf8ef5c (diff) |
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
-rw-r--r-- | src/BrillLindquist.c | 2 |
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; |