From f254b589c08bd1e846594a60e9043f39263ea9d4 Mon Sep 17 00:00:00 2001 From: schnetter Date: Sun, 7 May 2006 22:09:05 +0000 Subject: Use the correct power of TP_epsilon for the extrinsic curvature git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@55 b2a53a04-0f4f-0410-87ed-f9f25ced00cf --- src/Equations.c | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'src') diff --git a/src/Equations.c b/src/Equations.c index 17032c4..42a2a08 100644 --- a/src/Equations.c +++ b/src/Equations.c @@ -88,8 +88,8 @@ BY_Aijofxyz (double x, double y, double z, double Aij[3][3]) r2_plus = (x - par_b) * (x - par_b) + y * y + z * z; r2_minus = (x + par_b) * (x + par_b) + y * y + z * z; - r2_plus = sqrt (pow (r2_plus, 2) + pow (TP_epsilon, 2)); - r2_minus = sqrt (pow (r2_minus, 2) + pow (TP_epsilon, 2)); + r2_plus = sqrt (pow (r2_plus, 2) + pow (TP_epsilon, 4)); + r2_minus = sqrt (pow (r2_minus, 2) + pow (TP_epsilon, 4)); if (r2_plus < TP_Tiny) r2_plus = TP_Tiny; if (r2_minus < TP_Tiny) -- cgit v1.2.3