aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorallen <allen@6a3ddf76-46e1-4315-99d9-bc56cac1ef84>2002-11-01 08:31:35 +0000
committerallen <allen@6a3ddf76-46e1-4315-99d9-bc56cac1ef84>2002-11-01 08:31:35 +0000
commit03633b15df26458293cd074b2804bc6e199b6cfb (patch)
tree298e2ec0c5aebc78abb2a3d2efa5f851d215dfd5
parent924216071673e4de7f7de2468b974bcc43a3837c (diff)
Fix to a bug in the Schwarzschild metric for the physical case (Note: this is just
in the physical metric, the conformal (usual) metric is fine!). Spotted and fixed by Frank Herrmann. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@123 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
-rw-r--r--src/Schwarzschild.c4
1 files changed, 2 insertions, 2 deletions
diff --git a/src/Schwarzschild.c b/src/Schwarzschild.c
index fa442e8..4fc33dd 100644
--- a/src/Schwarzschild.c
+++ b/src/Schwarzschild.c
@@ -91,8 +91,8 @@ void Schwarzschild(CCTK_ARGUMENTS)
{
for (i = 0; i < npoints; i++)
{
- r_squared = r[i] * r[i];
- gxx[i] = one + mass/two/(r_squared * r_squared);
+ tmp = one + mass/two/r[i];
+ gxx[i] = tmp*tmp*tmp*tmp;
gyy[i] = gxx[i];
gzz[i] = gxx[i];
gxy[i] = zero;