aboutsummaryrefslogtreecommitdiff
path: root/ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_RHS2.c
diff options
context:
space:
mode:
authorErik Schnetter <schnetter@cct.lsu.edu>2011-01-03 21:01:50 -0500
committerErik Schnetter <schnetter@cct.lsu.edu>2011-01-03 21:01:50 -0500
commitc14ee6f3d20ffb6e9099de25f03a5c2f49544dab (patch)
treefe1e4d9818d677b88e1d4258027a25e580d2fa69 /ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_RHS2.c
parentfcd47f9d3637004987ac59562574f9e3d09e1a5c (diff)
Re-generate code
Diffstat (limited to 'ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_RHS2.c')
-rw-r--r--ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_RHS2.c27
1 files changed, 12 insertions, 15 deletions
diff --git a/ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_RHS2.c b/ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_RHS2.c
index 526974b..4b8a2d2 100644
--- a/ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_RHS2.c
+++ b/ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_RHS2.c
@@ -50,9 +50,6 @@ void ML_BSSN_MP_O8_RHS2_Body(cGH const * restrict const cctkGH, int const dir, i
return;
}
- const char *groups[] = {"Coordinates::jacobian","Coordinates::jacobian2","ML_BSSN_MP_O8::ML_curv","ML_BSSN_MP_O8::ML_curvrhs","ML_BSSN_MP_O8::ML_Gamma","ML_BSSN_MP_O8::ML_lapse","ML_BSSN_MP_O8::ML_log_confac","ML_BSSN_MP_O8::ML_metric","ML_BSSN_MP_O8::ML_shift","ML_BSSN_MP_O8::ML_trace_curv"};
- GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_MP_O8_RHS2", 10, groups);
-
/* Include user-supplied include files */
/* Initialise finite differencing variables */
@@ -510,10 +507,10 @@ void ML_BSSN_MP_O8_RHS2_Body(cGH const * restrict const cctkGH, int const dir, i
J21L*khalf*PDstandardNth2gt33 + J33L*PDstandardNth3gt13 -
J31L*khalf*PDstandardNth3gt33;
- CCTK_REAL Gtl211 = khalf*(-(J12L*PDstandardNth1gt11) +
- 2*J11L*PDstandardNth1gt12 - J22L*PDstandardNth2gt11 +
- 2*J21L*PDstandardNth2gt12 - J32L*PDstandardNth3gt11 +
- 2*J31L*PDstandardNth3gt12);
+ CCTK_REAL Gtl211 = -(J12L*khalf*PDstandardNth1gt11) +
+ J11L*PDstandardNth1gt12 - J22L*khalf*PDstandardNth2gt11 +
+ J21L*PDstandardNth2gt12 - J32L*khalf*PDstandardNth3gt11 +
+ J31L*PDstandardNth3gt12;
CCTK_REAL Gtl212 = khalf*(J11L*PDstandardNth1gt22 +
J21L*PDstandardNth2gt22 + J31L*PDstandardNth3gt22);
@@ -535,10 +532,10 @@ void ML_BSSN_MP_O8_RHS2_Body(cGH const * restrict const cctkGH, int const dir, i
J22L*khalf*PDstandardNth2gt33 + J33L*PDstandardNth3gt23 -
J32L*khalf*PDstandardNth3gt33;
- CCTK_REAL Gtl311 = khalf*(-(J13L*PDstandardNth1gt11) +
- 2*J11L*PDstandardNth1gt13 - J23L*PDstandardNth2gt11 +
- 2*J21L*PDstandardNth2gt13 - J33L*PDstandardNth3gt11 +
- 2*J31L*PDstandardNth3gt13);
+ CCTK_REAL Gtl311 = -(J13L*khalf*PDstandardNth1gt11) +
+ J11L*PDstandardNth1gt13 - J23L*khalf*PDstandardNth2gt11 +
+ J21L*PDstandardNth2gt13 - J33L*khalf*PDstandardNth3gt11 +
+ J31L*PDstandardNth3gt13;
CCTK_REAL Gtl312 = khalf*(-(J13L*PDstandardNth1gt12) +
J12L*PDstandardNth1gt13 + J11L*PDstandardNth1gt23 -
@@ -549,10 +546,10 @@ void ML_BSSN_MP_O8_RHS2_Body(cGH const * restrict const cctkGH, int const dir, i
CCTK_REAL Gtl313 = khalf*(J11L*PDstandardNth1gt33 +
J21L*PDstandardNth2gt33 + J31L*PDstandardNth3gt33);
- CCTK_REAL Gtl322 = khalf*(-(J13L*PDstandardNth1gt22) +
- 2*J12L*PDstandardNth1gt23 - J23L*PDstandardNth2gt22 +
- 2*J22L*PDstandardNth2gt23 - J33L*PDstandardNth3gt22 +
- 2*J32L*PDstandardNth3gt23);
+ CCTK_REAL Gtl322 = -(J13L*khalf*PDstandardNth1gt22) +
+ J12L*PDstandardNth1gt23 - J23L*khalf*PDstandardNth2gt22 +
+ J22L*PDstandardNth2gt23 - J33L*khalf*PDstandardNth3gt22 +
+ J32L*PDstandardNth3gt23;
CCTK_REAL Gtl323 = khalf*(J12L*PDstandardNth1gt33 +
J22L*PDstandardNth2gt33 + J32L*PDstandardNth3gt33);