diff options
| author | schnetter <schnetter@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2006-03-22 11:18:14 +0000 |
|---|---|---|
| committer | schnetter <schnetter@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2006-03-22 11:18:14 +0000 |
| commit | a1a9c98fc4f1e907341844611d37fbddf7aca23a (patch) | |
| tree | 7279f1cc8465f1c2c32c1d4aff61a5da6e91bbf0 /src | |
| parent | 4c9ade9003c1c8b4ef7784350956d0fb3766853e (diff) | |
Smooth out Aij as well
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@53 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
Diffstat (limited to 'src')
| -rw-r--r-- | src/Equations.c | 2 |
1 files changed, 2 insertions, 0 deletions
diff --git a/src/Equations.c b/src/Equations.c index 4bbcbef..17032c4 100644 --- a/src/Equations.c +++ b/src/Equations.c @@ -88,6 +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)); if (r2_plus < TP_Tiny) r2_plus = TP_Tiny; if (r2_minus < TP_Tiny) |