From 8dac0f5f80630eb093da9f5183aa5143382184a4 Mon Sep 17 00:00:00 2001 From: Peter Diener Date: Tue, 29 Jan 2013 16:03:19 -0600 Subject: 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 --- m/McLachlan_BSSN.m | 1 + 1 file changed, 1 insertion(+) 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], -- cgit v1.2.3