diff options
author | schnetter <schnetter@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2006-05-07 22:09:05 +0000 |
---|---|---|
committer | schnetter <schnetter@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2006-05-07 22:09:05 +0000 |
commit | f254b589c08bd1e846594a60e9043f39263ea9d4 (patch) | |
tree | f2f815ca18bbb64f70db51c45520bd13e0320c84 /src | |
parent | 7fd9ab1996d56de722d2945e298c85d3751ca5db (diff) |
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
Diffstat (limited to 'src')
-rw-r--r-- | src/Equations.c | 4 |
1 files changed, 2 insertions, 2 deletions
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) |