aboutsummaryrefslogtreecommitdiff
path: root/ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_constraints2.c
diff options
context:
space:
mode:
Diffstat (limited to 'ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_constraints2.c')
-rw-r--r--ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_constraints2.c285
1 files changed, 129 insertions, 156 deletions
diff --git a/ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_constraints2.c b/ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_constraints2.c
index fb5d2e5..5733cd9 100644
--- a/ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_constraints2.c
+++ b/ML_BSSN_MP_O8/src/ML_BSSN_MP_O8_constraints2.c
@@ -26,6 +26,12 @@ void ML_BSSN_MP_O8_constraints2_SelectBCs(CCTK_ARGUMENTS)
DECLARE_CCTK_PARAMETERS;
CCTK_INT ierr = 0;
+ ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN_MP_O8::ML_cons_detg","flat");
+ if (ierr < 0)
+ CCTK_WARN(1, "Failed to register flat BC for ML_BSSN_MP_O8::ML_cons_detg.");
+ ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN_MP_O8::ML_cons_Gamma","flat");
+ if (ierr < 0)
+ CCTK_WARN(1, "Failed to register flat BC for ML_BSSN_MP_O8::ML_cons_Gamma.");
ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN_MP_O8::ML_cons_traceA","flat");
if (ierr < 0)
CCTK_WARN(1, "Failed to register flat BC for ML_BSSN_MP_O8::ML_cons_traceA.");
@@ -53,9 +59,6 @@ void ML_BSSN_MP_O8_constraints2_Body(cGH const * restrict const cctkGH, int cons
return;
}
- const char *groups[] = {"Coordinates::jacobian","ML_BSSN_MP_O8::ML_cons_traceA","ML_BSSN_MP_O8::ML_curv","ML_BSSN_MP_O8::ML_lapse","ML_BSSN_MP_O8::ML_log_confac","ML_BSSN_MP_O8::ML_metric","ML_BSSN_MP_O8::ML_mom","ML_BSSN_MP_O8::ML_shift","ML_BSSN_MP_O8::ML_trace_curv"};
- GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_MP_O8_constraints2", 9, groups);
-
/* Include user-supplied include files */
/* Initialise finite differencing variables */
@@ -192,6 +195,9 @@ void ML_BSSN_MP_O8_constraints2_Body(cGH const * restrict const cctkGH, int cons
CCTK_REAL J33L = J33[index];
CCTK_REAL phiL = phi[index];
CCTK_REAL trKL = trK[index];
+ CCTK_REAL Xt1L = Xt1[index];
+ CCTK_REAL Xt2L = Xt2[index];
+ CCTK_REAL Xt3L = Xt3[index];
/* Include user supplied include files */
@@ -254,164 +260,116 @@ void ML_BSSN_MP_O8_constraints2_Body(cGH const * restrict const cctkGH, int cons
CCTK_REAL gtu33 = INV(detgt)*(gt11L*gt22L - SQR(gt12L));
- CCTK_REAL Gt111 = khalf*((gtu11*J11L - gtu12*J12L -
- gtu13*J13L)*PDstandardNth1gt11 + (gtu11*J21L - gtu12*J22L -
- gtu13*J23L)*PDstandardNth2gt11 + (gtu11*J31L - gtu12*J32L -
- gtu13*J33L)*PDstandardNth3gt11 + 2*(J11L*(gtu12*PDstandardNth1gt12 +
- gtu13*PDstandardNth1gt13) + J21L*(gtu12*PDstandardNth2gt12 +
- gtu13*PDstandardNth2gt13) + J31L*(gtu12*PDstandardNth3gt12 +
- gtu13*PDstandardNth3gt13)));
-
- CCTK_REAL Gt211 = khalf*((gtu12*J11L - gtu22*J12L -
- gtu23*J13L)*PDstandardNth1gt11 + (gtu12*J21L - gtu22*J22L -
- gtu23*J23L)*PDstandardNth2gt11 + (gtu12*J31L - gtu22*J32L -
- gtu23*J33L)*PDstandardNth3gt11 + 2*(J11L*(gtu22*PDstandardNth1gt12 +
- gtu23*PDstandardNth1gt13) + J21L*(gtu22*PDstandardNth2gt12 +
- gtu23*PDstandardNth2gt13) + J31L*(gtu22*PDstandardNth3gt12 +
- gtu23*PDstandardNth3gt13)));
-
- CCTK_REAL Gt311 = khalf*((gtu13*J11L - gtu23*J12L -
- gtu33*J13L)*PDstandardNth1gt11 + (gtu13*J21L - gtu23*J22L -
- gtu33*J23L)*PDstandardNth2gt11 + (gtu13*J31L - gtu23*J32L -
- gtu33*J33L)*PDstandardNth3gt11 + 2*(J11L*(gtu23*PDstandardNth1gt12 +
- gtu33*PDstandardNth1gt13) + J21L*(gtu23*PDstandardNth2gt12 +
- gtu33*PDstandardNth2gt13) + J31L*(gtu23*PDstandardNth3gt12 +
- gtu33*PDstandardNth3gt13)));
-
- CCTK_REAL Gt112 = khalf*(gtu11*(J12L*PDstandardNth1gt11 +
- J22L*PDstandardNth2gt11 + J32L*PDstandardNth3gt11) +
- gtu12*(J11L*PDstandardNth1gt22 + J21L*PDstandardNth2gt22 +
- J31L*PDstandardNth3gt22) + gtu13*(-(J13L*PDstandardNth1gt12) +
- J12L*PDstandardNth1gt13 + J11L*PDstandardNth1gt23 -
- J23L*PDstandardNth2gt12 + J22L*PDstandardNth2gt13 +
- J21L*PDstandardNth2gt23 - J33L*PDstandardNth3gt12 +
- J32L*PDstandardNth3gt13 + J31L*PDstandardNth3gt23));
+ CCTK_REAL Gtl111 = khalf*(J11L*PDstandardNth1gt11 +
+ J21L*PDstandardNth2gt11 + J31L*PDstandardNth3gt11);
- CCTK_REAL Gt212 = khalf*(gtu12*(J12L*PDstandardNth1gt11 +
- J22L*PDstandardNth2gt11 + J32L*PDstandardNth3gt11) +
- gtu22*(J11L*PDstandardNth1gt22 + J21L*PDstandardNth2gt22 +
- J31L*PDstandardNth3gt22) + gtu23*(-(J13L*PDstandardNth1gt12) +
- J12L*PDstandardNth1gt13 + J11L*PDstandardNth1gt23 -
- J23L*PDstandardNth2gt12 + J22L*PDstandardNth2gt13 +
- J21L*PDstandardNth2gt23 - J33L*PDstandardNth3gt12 +
- J32L*PDstandardNth3gt13 + J31L*PDstandardNth3gt23));
+ CCTK_REAL Gtl112 = khalf*(J12L*PDstandardNth1gt11 +
+ J22L*PDstandardNth2gt11 + J32L*PDstandardNth3gt11);
+
+ CCTK_REAL Gtl113 = khalf*(J13L*PDstandardNth1gt11 +
+ J23L*PDstandardNth2gt11 + J33L*PDstandardNth3gt11);
+
+ CCTK_REAL Gtl122 = J12L*PDstandardNth1gt12 -
+ J11L*khalf*PDstandardNth1gt22 + J22L*PDstandardNth2gt12 -
+ J21L*khalf*PDstandardNth2gt22 + J32L*PDstandardNth3gt12 -
+ J31L*khalf*PDstandardNth3gt22;
+
+ CCTK_REAL Gtl123 = khalf*(J13L*PDstandardNth1gt12 +
+ J12L*PDstandardNth1gt13 - J11L*PDstandardNth1gt23 +
+ J23L*PDstandardNth2gt12 + J22L*PDstandardNth2gt13 -
+ J21L*PDstandardNth2gt23 + J33L*PDstandardNth3gt12 +
+ J32L*PDstandardNth3gt13 - J31L*PDstandardNth3gt23);
- CCTK_REAL Gt312 = khalf*(gtu13*(J12L*PDstandardNth1gt11 +
- J22L*PDstandardNth2gt11 + J32L*PDstandardNth3gt11) +
- gtu23*(J11L*PDstandardNth1gt22 + J21L*PDstandardNth2gt22 +
- J31L*PDstandardNth3gt22) + gtu33*(-(J13L*PDstandardNth1gt12) +
+ CCTK_REAL Gtl133 = J13L*PDstandardNth1gt13 -
+ J11L*khalf*PDstandardNth1gt33 + J23L*PDstandardNth2gt13 -
+ J21L*khalf*PDstandardNth2gt33 + J33L*PDstandardNth3gt13 -
+ J31L*khalf*PDstandardNth3gt33;
+
+ 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);
+
+ CCTK_REAL Gtl213 = khalf*(J13L*PDstandardNth1gt12 -
+ J12L*PDstandardNth1gt13 + J11L*PDstandardNth1gt23 +
+ J23L*PDstandardNth2gt12 - J22L*PDstandardNth2gt13 +
+ J21L*PDstandardNth2gt23 + J33L*PDstandardNth3gt12 -
+ J32L*PDstandardNth3gt13 + J31L*PDstandardNth3gt23);
+
+ CCTK_REAL Gtl222 = khalf*(J12L*PDstandardNth1gt22 +
+ J22L*PDstandardNth2gt22 + J32L*PDstandardNth3gt22);
+
+ CCTK_REAL Gtl223 = khalf*(J13L*PDstandardNth1gt22 +
+ J23L*PDstandardNth2gt22 + J33L*PDstandardNth3gt22);
+
+ CCTK_REAL Gtl233 = J13L*PDstandardNth1gt23 -
+ J12L*khalf*PDstandardNth1gt33 + J23L*PDstandardNth2gt23 -
+ J22L*khalf*PDstandardNth2gt33 + J33L*PDstandardNth3gt23 -
+ J32L*khalf*PDstandardNth3gt33;
+
+ 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 -
J23L*PDstandardNth2gt12 + J22L*PDstandardNth2gt13 +
J21L*PDstandardNth2gt23 - J33L*PDstandardNth3gt12 +
- J32L*PDstandardNth3gt13 + J31L*PDstandardNth3gt23));
-
- CCTK_REAL Gt113 = khalf*(gtu11*(J13L*PDstandardNth1gt11 +
- J23L*PDstandardNth2gt11 + J33L*PDstandardNth3gt11) +
- gtu12*(J13L*PDstandardNth1gt12 - J12L*PDstandardNth1gt13 +
- J11L*PDstandardNth1gt23 + J23L*PDstandardNth2gt12 -
- J22L*PDstandardNth2gt13 + J21L*PDstandardNth2gt23 +
- J33L*PDstandardNth3gt12 - J32L*PDstandardNth3gt13 +
- J31L*PDstandardNth3gt23) + gtu13*(J11L*PDstandardNth1gt33 +
- J21L*PDstandardNth2gt33 + J31L*PDstandardNth3gt33));
-
- CCTK_REAL Gt213 = khalf*(gtu12*(J13L*PDstandardNth1gt11 +
- J23L*PDstandardNth2gt11 + J33L*PDstandardNth3gt11) +
- gtu22*(J13L*PDstandardNth1gt12 - J12L*PDstandardNth1gt13 +
- J11L*PDstandardNth1gt23 + J23L*PDstandardNth2gt12 -
- J22L*PDstandardNth2gt13 + J21L*PDstandardNth2gt23 +
- J33L*PDstandardNth3gt12 - J32L*PDstandardNth3gt13 +
- J31L*PDstandardNth3gt23) + gtu23*(J11L*PDstandardNth1gt33 +
- J21L*PDstandardNth2gt33 + J31L*PDstandardNth3gt33));
-
- CCTK_REAL Gt313 = khalf*(gtu13*(J13L*PDstandardNth1gt11 +
- J23L*PDstandardNth2gt11 + J33L*PDstandardNth3gt11) +
- gtu23*(J13L*PDstandardNth1gt12 - J12L*PDstandardNth1gt13 +
- J11L*PDstandardNth1gt23 + J23L*PDstandardNth2gt12 -
- J22L*PDstandardNth2gt13 + J21L*PDstandardNth2gt23 +
- J33L*PDstandardNth3gt12 - J32L*PDstandardNth3gt13 +
- J31L*PDstandardNth3gt23) + gtu33*(J11L*PDstandardNth1gt33 +
- J21L*PDstandardNth2gt33 + J31L*PDstandardNth3gt33));
-
- CCTK_REAL Gt122 = khalf*(gtu11*(-(J11L*PDstandardNth1gt22) +
- 2*(J12L*PDstandardNth1gt12 + J22L*PDstandardNth2gt12) -
- J21L*PDstandardNth2gt22 + 2*J32L*PDstandardNth3gt12 -
- J31L*PDstandardNth3gt22) + gtu12*(J12L*PDstandardNth1gt22 +
- J22L*PDstandardNth2gt22 + J32L*PDstandardNth3gt22) -
- gtu13*(J13L*PDstandardNth1gt22 + J23L*PDstandardNth2gt22 +
- J33L*PDstandardNth3gt22 - 2*(J12L*PDstandardNth1gt23 +
- J22L*PDstandardNth2gt23 + J32L*PDstandardNth3gt23)));
-
- CCTK_REAL Gt222 = khalf*(gtu12*(-(J11L*PDstandardNth1gt22) +
- 2*(J12L*PDstandardNth1gt12 + J22L*PDstandardNth2gt12) -
- J21L*PDstandardNth2gt22 + 2*J32L*PDstandardNth3gt12 -
- J31L*PDstandardNth3gt22) + gtu22*(J12L*PDstandardNth1gt22 +
- J22L*PDstandardNth2gt22 + J32L*PDstandardNth3gt22) -
- gtu23*(J13L*PDstandardNth1gt22 + J23L*PDstandardNth2gt22 +
- J33L*PDstandardNth3gt22 - 2*(J12L*PDstandardNth1gt23 +
- J22L*PDstandardNth2gt23 + J32L*PDstandardNth3gt23)));
-
- CCTK_REAL Gt322 = khalf*(gtu13*(-(J11L*PDstandardNth1gt22) +
- 2*(J12L*PDstandardNth1gt12 + J22L*PDstandardNth2gt12) -
- J21L*PDstandardNth2gt22 + 2*J32L*PDstandardNth3gt12 -
- J31L*PDstandardNth3gt22) + gtu23*(J12L*PDstandardNth1gt22 +
- J22L*PDstandardNth2gt22 + J32L*PDstandardNth3gt22) -
- gtu33*(J13L*PDstandardNth1gt22 + J23L*PDstandardNth2gt22 +
- J33L*PDstandardNth3gt22 - 2*(J12L*PDstandardNth1gt23 +
- J22L*PDstandardNth2gt23 + J32L*PDstandardNth3gt23)));
-
- CCTK_REAL Gt123 = khalf*(gtu12*(J13L*PDstandardNth1gt22 +
- J23L*PDstandardNth2gt22 + J33L*PDstandardNth3gt22) +
- gtu11*(J13L*PDstandardNth1gt12 + J12L*PDstandardNth1gt13 -
- J11L*PDstandardNth1gt23 + J23L*PDstandardNth2gt12 +
- J22L*PDstandardNth2gt13 - J21L*PDstandardNth2gt23 +
- J33L*PDstandardNth3gt12 + J32L*PDstandardNth3gt13 -
- J31L*PDstandardNth3gt23) + gtu13*(J12L*PDstandardNth1gt33 +
- J22L*PDstandardNth2gt33 + J32L*PDstandardNth3gt33));
-
- CCTK_REAL Gt223 = khalf*(gtu22*(J13L*PDstandardNth1gt22 +
- J23L*PDstandardNth2gt22 + J33L*PDstandardNth3gt22) +
- gtu12*(J13L*PDstandardNth1gt12 + J12L*PDstandardNth1gt13 -
- J11L*PDstandardNth1gt23 + J23L*PDstandardNth2gt12 +
- J22L*PDstandardNth2gt13 - J21L*PDstandardNth2gt23 +
- J33L*PDstandardNth3gt12 + J32L*PDstandardNth3gt13 -
- J31L*PDstandardNth3gt23) + gtu23*(J12L*PDstandardNth1gt33 +
- J22L*PDstandardNth2gt33 + J32L*PDstandardNth3gt33));
-
- CCTK_REAL Gt323 = khalf*(gtu23*(J13L*PDstandardNth1gt22 +
- J23L*PDstandardNth2gt22 + J33L*PDstandardNth3gt22) +
- gtu13*(J13L*PDstandardNth1gt12 + J12L*PDstandardNth1gt13 -
- J11L*PDstandardNth1gt23 + J23L*PDstandardNth2gt12 +
- J22L*PDstandardNth2gt13 - J21L*PDstandardNth2gt23 +
- J33L*PDstandardNth3gt12 + J32L*PDstandardNth3gt13 -
- J31L*PDstandardNth3gt23) + gtu33*(J12L*PDstandardNth1gt33 +
- J22L*PDstandardNth2gt33 + J32L*PDstandardNth3gt33));
-
- CCTK_REAL Gt133 = khalf*(gtu11*(-(J11L*PDstandardNth1gt33) +
- 2*(J13L*PDstandardNth1gt13 + J23L*PDstandardNth2gt13) -
- J21L*PDstandardNth2gt33 + 2*J33L*PDstandardNth3gt13 -
- J31L*PDstandardNth3gt33) + gtu12*(-(J12L*PDstandardNth1gt33) +
- 2*(J13L*PDstandardNth1gt23 + J23L*PDstandardNth2gt23) -
- J22L*PDstandardNth2gt33 + 2*J33L*PDstandardNth3gt23 -
- J32L*PDstandardNth3gt33) + gtu13*(J13L*PDstandardNth1gt33 +
- J23L*PDstandardNth2gt33 + J33L*PDstandardNth3gt33));
-
- CCTK_REAL Gt233 = khalf*(gtu12*(-(J11L*PDstandardNth1gt33) +
- 2*(J13L*PDstandardNth1gt13 + J23L*PDstandardNth2gt13) -
- J21L*PDstandardNth2gt33 + 2*J33L*PDstandardNth3gt13 -
- J31L*PDstandardNth3gt33) + gtu22*(-(J12L*PDstandardNth1gt33) +
- 2*(J13L*PDstandardNth1gt23 + J23L*PDstandardNth2gt23) -
- J22L*PDstandardNth2gt33 + 2*J33L*PDstandardNth3gt23 -
- J32L*PDstandardNth3gt33) + gtu23*(J13L*PDstandardNth1gt33 +
- J23L*PDstandardNth2gt33 + J33L*PDstandardNth3gt33));
-
- CCTK_REAL Gt333 = khalf*(gtu13*(-(J11L*PDstandardNth1gt33) +
- 2*(J13L*PDstandardNth1gt13 + J23L*PDstandardNth2gt13) -
- J21L*PDstandardNth2gt33 + 2*J33L*PDstandardNth3gt13 -
- J31L*PDstandardNth3gt33) + gtu23*(-(J12L*PDstandardNth1gt33) +
- 2*(J13L*PDstandardNth1gt23 + J23L*PDstandardNth2gt23) -
- J22L*PDstandardNth2gt33 + 2*J33L*PDstandardNth3gt23 -
- J32L*PDstandardNth3gt33) + gtu33*(J13L*PDstandardNth1gt33 +
- J23L*PDstandardNth2gt33 + J33L*PDstandardNth3gt33));
+ J32L*PDstandardNth3gt13 + J31L*PDstandardNth3gt23);
+
+ CCTK_REAL Gtl313 = khalf*(J11L*PDstandardNth1gt33 +
+ J21L*PDstandardNth2gt33 + J31L*PDstandardNth3gt33);
+
+ 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);
+
+ CCTK_REAL Gtl333 = khalf*(J13L*PDstandardNth1gt33 +
+ J23L*PDstandardNth2gt33 + J33L*PDstandardNth3gt33);
+
+ CCTK_REAL Gt111 = Gtl111*gtu11 + Gtl211*gtu12 + Gtl311*gtu13;
+
+ CCTK_REAL Gt211 = Gtl111*gtu12 + Gtl211*gtu22 + Gtl311*gtu23;
+
+ CCTK_REAL Gt311 = Gtl111*gtu13 + Gtl211*gtu23 + Gtl311*gtu33;
+
+ CCTK_REAL Gt112 = Gtl112*gtu11 + Gtl212*gtu12 + Gtl312*gtu13;
+
+ CCTK_REAL Gt212 = Gtl112*gtu12 + Gtl212*gtu22 + Gtl312*gtu23;
+
+ CCTK_REAL Gt312 = Gtl112*gtu13 + Gtl212*gtu23 + Gtl312*gtu33;
+
+ CCTK_REAL Gt113 = Gtl113*gtu11 + Gtl213*gtu12 + Gtl313*gtu13;
+
+ CCTK_REAL Gt213 = Gtl113*gtu12 + Gtl213*gtu22 + Gtl313*gtu23;
+
+ CCTK_REAL Gt313 = Gtl113*gtu13 + Gtl213*gtu23 + Gtl313*gtu33;
+
+ CCTK_REAL Gt122 = Gtl122*gtu11 + Gtl222*gtu12 + Gtl322*gtu13;
+
+ CCTK_REAL Gt222 = Gtl122*gtu12 + Gtl222*gtu22 + Gtl322*gtu23;
+
+ CCTK_REAL Gt322 = Gtl122*gtu13 + Gtl222*gtu23 + Gtl322*gtu33;
+
+ CCTK_REAL Gt123 = Gtl123*gtu11 + Gtl223*gtu12 + Gtl323*gtu13;
+
+ CCTK_REAL Gt223 = Gtl123*gtu12 + Gtl223*gtu22 + Gtl323*gtu23;
+
+ CCTK_REAL Gt323 = Gtl123*gtu13 + Gtl223*gtu23 + Gtl323*gtu33;
+
+ CCTK_REAL Gt133 = Gtl133*gtu11 + Gtl233*gtu12 + Gtl333*gtu13;
+
+ CCTK_REAL Gt233 = Gtl133*gtu12 + Gtl233*gtu22 + Gtl333*gtu23;
+
+ CCTK_REAL Gt333 = Gtl133*gtu13 + Gtl233*gtu23 + Gtl333*gtu33;
CCTK_REAL fac1 = IfThen(conformalMethod,-(khalf*INV(phiL)),1);
@@ -527,12 +485,27 @@ void ML_BSSN_MP_O8_constraints2_Body(cGH const * restrict const cctkGH, int cons
J33L*ktwothird*PDstandardNth3trK -
25.13274122871834590770114706623602307358*S3;
+ CCTK_REAL cSL = Log(detgt);
+
+ CCTK_REAL cXt1L = Gt111*gtu11 + Gt122*gtu22 + 2*(Gt112*gtu12 +
+ Gt113*gtu13 + Gt123*gtu23) + Gt133*gtu33 - Xt1L;
+
+ CCTK_REAL cXt2L = Gt211*gtu11 + Gt222*gtu22 + 2*(Gt212*gtu12 +
+ Gt213*gtu13 + Gt223*gtu23) + Gt233*gtu33 - Xt2L;
+
+ CCTK_REAL cXt3L = Gt311*gtu11 + Gt322*gtu22 + 2*(Gt312*gtu12 +
+ Gt313*gtu13 + Gt323*gtu23) + Gt333*gtu33 - Xt3L;
+
CCTK_REAL cAL = At11L*gtu11 + At22L*gtu22 + 2*(At12L*gtu12 +
At13L*gtu13 + At23L*gtu23) + At33L*gtu33;
/* Copy local copies back to grid functions */
cA[index] = cAL;
+ cS[index] = cSL;
+ cXt1[index] = cXt1L;
+ cXt2[index] = cXt2L;
+ cXt3[index] = cXt3L;
M1[index] = M1L;
M2[index] = M2L;
M3[index] = M3L;
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-03 06:29:56 +0000