diff options
author | schnetter <schnetter@6a3ddf76-46e1-4315-99d9-bc56cac1ef84> | 2004-06-23 07:13:24 +0000 |
---|---|---|
committer | schnetter <schnetter@6a3ddf76-46e1-4315-99d9-bc56cac1ef84> | 2004-06-23 07:13:24 +0000 |
commit | e340f99ace935245fb42ee42a8e2225bce351d65 (patch) | |
tree | 8214c75c0ac20a0b9902f989401ab710f9f674ba /src | |
parent | 374fbfbc1f5ece3be910bfd1359fc9f595f24ccf (diff) |
Add lapse "schwarzschild" as an alias for "schwarz".
Do not produce nans on the origin.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@150 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
Diffstat (limited to 'src')
-rw-r--r-- | src/Schwarzschild.c | 12 |
1 files changed, 7 insertions, 5 deletions
diff --git a/src/Schwarzschild.c b/src/Schwarzschild.c index 5eb4c6f..850bd9a 100644 --- a/src/Schwarzschild.c +++ b/src/Schwarzschild.c @@ -1,3 +1,4 @@ + /*@@ @file Schwarzschild.c @date Sun Oct 17 10:35:41 1999 @@ -67,14 +68,14 @@ void Schwarzschild(CCTK_ARGUMENTS) for (i = 0; i < npoints; i++) { /* Compute conformal factor */ - psi[i] = ( one + mass/two/r[i]); + psi[i] = ( one + mass/(two*r[i]+1e-20)); if(make_conformal_derivs) { /* derivatives of psi / psi */ r_squared = r[i]*r[i]; r_cubed = r[i]*r_squared; - tmp = mass/two/r_cubed/psi[i]; + tmp = mass/(two*r_cubed*psi[i]+1e-20); psix[i] = -x[i]*tmp; @@ -83,7 +84,7 @@ void Schwarzschild(CCTK_ARGUMENTS) if(*conformal_state > 2) { - tmp = mass/two/(r_squared*r_cubed)/psi[i]; + tmp = mass/(two*r_squared*r_cubed*psi[i]+1e-20); psixy[i] = three*x[i]*y[i]*tmp; psixz[i] = three*x[i]*z[i]*tmp; psiyz[i] = three*y[i]*z[i]*tmp; @@ -105,7 +106,7 @@ void Schwarzschild(CCTK_ARGUMENTS) { for (i = 0; i < npoints; i++) { - tmp = one + mass/two/r[i]; + tmp = one + mass/(two*r[i]+1e-20); gxx[i] = tmp*tmp*tmp*tmp; gyy[i] = gxx[i]; gzz[i] = gxx[i]; @@ -116,7 +117,8 @@ void Schwarzschild(CCTK_ARGUMENTS) } /* If the initial lapse is not one ... */ - if (CCTK_Equals(initial_lapse,"schwarz")) + if (CCTK_Equals(initial_lapse,"schwarzschild") + || CCTK_Equals(initial_lapse,"schwarz")) { CCTK_INFO("Initialise with Schwarzschild lapse"); |