From 03633b15df26458293cd074b2804bc6e199b6cfb Mon Sep 17 00:00:00 2001 From: allen Date: Fri, 1 Nov 2002 08:31:35 +0000 Subject: 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 --- src/Schwarzschild.c | 4 ++-- 1 file 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; -- cgit v1.2.3