diff options
| author | Erik Schnetter <schnetter@gmail.com> | 2012-02-05 21:54:48 -0500 |
|---|---|---|
| committer | Erik Schnetter <schnetter@gmail.com> | 2012-02-05 21:54:48 -0500 |
| commit | 45ec267196a8f999bd0bc0b72a18ce2bbcdbf4af (patch) | |
| tree | 4491a43bf978eae93c5f633e5dbadceff3c0517a /m | |
| parent | 21a06b511674e1b75e05b83981fb7ad71d083ae2 (diff) | |
Correct etaExpr to avoid (inconsequential) division by zero
Diffstat (limited to 'm')
| -rw-r--r-- | m/McLachlan_BSSN.m | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/m/McLachlan_BSSN.m b/m/McLachlan_BSSN.m index a8f6984..2d4eeb3 100644 --- a/m/McLachlan_BSSN.m +++ b/m/McLachlan_BSSN.m @@ -169,7 +169,7 @@ ddetgtExpr[la_] = Sum [D[Det[MatrixOfComponents[gt[la, lb]]], X] PD[X, la], {X, Union[Flatten[MatrixOfComponents[gt[la, lb]]]]}]; -etaExpr = Min [SpatialBetaDriverRadius / r, 1]; +etaExpr = SpatialBetaDriverRadius / Max [r, SpatialBetaDriverRadius]; thetaExpr = Min [Exp [1 - r / SpatialShiftGammaCoeffRadius], 1]; |