aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorPeter Diener <diener@cct.lsu.edu>2013-01-29 16:03:19 -0600
committerPeter Diener <diener@cct.lsu.edu>2013-01-29 16:03:19 -0600
commit8dac0f5f80630eb093da9f5183aa5143382184a4 (patch)
treebb35f634c30dc6a4a1ee133a280dd90b436f1b29
parent3c5b9316279ea056c467d607558e4c16f7ab25e0 (diff)
Fix the sign for 2nd and 6th order dissipation operators.
Fix the sign for 2nd and 6th order dissipation operators. Signed-off-by: Peter Diener <diener@cct.lsu.edu>
-rw-r--r--m/McLachlan_BSSN.m1
1 files changed, 1 insertions, 0 deletions
diff --git a/m/McLachlan_BSSN.m b/m/McLachlan_BSSN.m
index 20ff871..a835560 100644
--- a/m/McLachlan_BSSN.m
+++ b/m/McLachlan_BSSN.m
@@ -37,6 +37,7 @@ derivatives =
PDstandardNth[i_,j_] -> StandardCenteredDifferenceOperator[1,fdOrder/2,i] *
StandardCenteredDifferenceOperator[1,fdOrder/2,j],
PDdissipationNth[i_] ->
+ (-1)^(fdOrder/2) *
spacing[i]^(fdOrder+1) / 2^(fdOrder+2) *
StandardCenteredDifferenceOperator[fdOrder+2,fdOrder/2+1,i],