diff options
| author | herrmann <herrmann@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2004-06-10 13:55:25 +0000 |
|---|---|---|
| committer | herrmann <herrmann@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2004-06-10 13:55:25 +0000 |
| commit | 29fbb0c8978abc82e4443dbf1a0c7b3cb52a82ba (patch) | |
| tree | 0a3ee981760bb039cee5ac2553bed6d1c44cd9ca /src | |
| parent | 4e9143deab72ea0213cdf8de81c10341fdead108 (diff) | |
changing computation of R for the "evaluation" case. without this change R
is NaN for x close to zero. was already done by Erik for the Taylor case.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@15 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
Diffstat (limited to 'src')
| -rw-r--r-- | src/FuncAndJacobian.c | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/src/FuncAndJacobian.c b/src/FuncAndJacobian.c index 96adeb4..7c30da0 100644 --- a/src/FuncAndJacobian.c +++ b/src/FuncAndJacobian.c @@ -725,7 +725,7 @@ PunctIntPolAtArbitPosition (int ivar, int nvar, int n1, aux1 = 0.5 * (xs * xs + rs2 - 1); aux2 = sqrt (aux1 * aux1 + rs2); X = asinh (sqrt (aux1 + aux2)); - R = asin (sqrt (-aux1 + aux2)); + R = asin (min(1.0, sqrt (-aux1 + aux2))); if (x < 0) R = Pi - R; |