aboutsummaryrefslogtreecommitdiff
path: root/ML_ADM
diff options
context:
space:
mode:
authorPeter Diener <diener@10-4-34-85.(none)>2007-11-09 15:19:26 -0600
committerPeter Diener <diener@10-4-34-85.(none)>2007-11-09 15:19:26 -0600
commitb08330b05518b8083dc368d012010868786e813c (patch)
treeca77546416cfdd73c9f0c79d483ab1f11141386b /ML_ADM
parent901aee61dc2a0b884b5982162a93d3023bfcb80e (diff)
parent3d4d942646f4628487e8786df216dafd766eccca (diff)
Merge branch 'master' of ssh://relativity.phys.lsu.edu/home/perturbed/gitroot/McLachlan
Conflicts: ML_BSSN/src/ML_BSSN_constraints.c m/McLachlan.m m/McLachlan.sh
Diffstat (limited to 'ML_ADM')
-rw-r--r--ML_ADM/src/Differencing.h24
-rw-r--r--ML_ADM/src/ML_ADM_RHS.c552
-rw-r--r--ML_ADM/src/ML_ADM_constraints.c438
3 files changed, 507 insertions, 507 deletions
diff --git a/ML_ADM/src/Differencing.h b/ML_ADM/src/Differencing.h
index 7ada6a7..876855b 100644
--- a/ML_ADM/src/Differencing.h
+++ b/ML_ADM/src/Differencing.h
@@ -1,12 +1,12 @@
-#define PDstandard4th1(u,i,j,k) (p1o12dx*(-8*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] + 8*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
-#define PDstandard4th2(u,i,j,k) (p1o12dy*(-8*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] + 8*u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))
-#define PDstandard4th3(u,i,j,k) (p1o12dz*(-8*u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] + 8*u[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))
-#define PDstandard4th11(u,i,j,k) (pm1o12dx2*(-16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,j,k)]))
-#define PDstandard4th22(u,i,j,k) (pm1o12dy2*(-16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)]) + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,j,k)]))
-#define PDstandard4th33(u,i,j,k) (pm1o12dz2*(-16*(u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,j,k)]))
-#define PDstandard4th12(u,i,j,k) (p1o144dxdy*(-64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)]))
-#define PDstandard4th13(u,i,j,k) (p1o144dxdz*(-64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)]))
-#define PDstandard4th21(u,i,j,k) (p1o144dxdy*(-64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)]))
-#define PDstandard4th23(u,i,j,k) (p1o144dydz*(-64*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)]))
-#define PDstandard4th31(u,i,j,k) (p1o144dxdz*(-64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)]))
-#define PDstandard4th32(u,i,j,k) (p1o144dydz*(-64*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)]))
+#define PDstandardNth1(u,i,j,k) (p1o12dx*(-8*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] + 8*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
+#define PDstandardNth2(u,i,j,k) (p1o12dy*(-8*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] + 8*u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))
+#define PDstandardNth3(u,i,j,k) (p1o12dz*(-8*u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] + 8*u[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))
+#define PDstandardNth11(u,i,j,k) (pm1o12dx2*(-16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,j,k)]))
+#define PDstandardNth22(u,i,j,k) (pm1o12dy2*(-16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)]) + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,j,k)]))
+#define PDstandardNth33(u,i,j,k) (pm1o12dz2*(-16*(u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,j,k)]))
+#define PDstandardNth12(u,i,j,k) (p1o144dxdy*(-64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)]))
+#define PDstandardNth13(u,i,j,k) (p1o144dxdz*(-64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)]))
+#define PDstandardNth21(u,i,j,k) (p1o144dxdy*(-64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)]))
+#define PDstandardNth23(u,i,j,k) (p1o144dydz*(-64*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)]))
+#define PDstandardNth31(u,i,j,k) (p1o144dxdz*(-64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)]))
+#define PDstandardNth32(u,i,j,k) (p1o144dydz*(-64*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)]))
diff --git a/ML_ADM/src/ML_ADM_RHS.c b/ML_ADM/src/ML_ADM_RHS.c
index cc16858..20b9f39 100644
--- a/ML_ADM/src/ML_ADM_RHS.c
+++ b/ML_ADM/src/ML_ADM_RHS.c
@@ -114,105 +114,105 @@ void ML_ADM_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[
/* Declare precomputed derivatives*/
/* Declare derivatives */
- CCTK_REAL PDstandard4th1alpha = INITVALUE;
- CCTK_REAL PDstandard4th2alpha = INITVALUE;
- CCTK_REAL PDstandard4th3alpha = INITVALUE;
- CCTK_REAL PDstandard4th11alpha = INITVALUE;
- CCTK_REAL PDstandard4th22alpha = INITVALUE;
- CCTK_REAL PDstandard4th33alpha = INITVALUE;
- CCTK_REAL PDstandard4th12alpha = INITVALUE;
- CCTK_REAL PDstandard4th13alpha = INITVALUE;
- CCTK_REAL PDstandard4th23alpha = INITVALUE;
- CCTK_REAL PDstandard4th1beta1 = INITVALUE;
- CCTK_REAL PDstandard4th2beta1 = INITVALUE;
- CCTK_REAL PDstandard4th3beta1 = INITVALUE;
- CCTK_REAL PDstandard4th1beta2 = INITVALUE;
- CCTK_REAL PDstandard4th2beta2 = INITVALUE;
- CCTK_REAL PDstandard4th3beta2 = INITVALUE;
- CCTK_REAL PDstandard4th1beta3 = INITVALUE;
- CCTK_REAL PDstandard4th2beta3 = INITVALUE;
- CCTK_REAL PDstandard4th3beta3 = INITVALUE;
- CCTK_REAL PDstandard4th1g11 = INITVALUE;
- CCTK_REAL PDstandard4th2g11 = INITVALUE;
- CCTK_REAL PDstandard4th3g11 = INITVALUE;
- CCTK_REAL PDstandard4th22g11 = INITVALUE;
- CCTK_REAL PDstandard4th33g11 = INITVALUE;
- CCTK_REAL PDstandard4th12g11 = INITVALUE;
- CCTK_REAL PDstandard4th13g11 = INITVALUE;
- CCTK_REAL PDstandard4th21g11 = INITVALUE;
- CCTK_REAL PDstandard4th23g11 = INITVALUE;
- CCTK_REAL PDstandard4th31g11 = INITVALUE;
- CCTK_REAL PDstandard4th32g11 = INITVALUE;
- CCTK_REAL PDstandard4th1g12 = INITVALUE;
- CCTK_REAL PDstandard4th2g12 = INITVALUE;
- CCTK_REAL PDstandard4th3g12 = INITVALUE;
- CCTK_REAL PDstandard4th33g12 = INITVALUE;
- CCTK_REAL PDstandard4th12g12 = INITVALUE;
- CCTK_REAL PDstandard4th13g12 = INITVALUE;
- CCTK_REAL PDstandard4th21g12 = INITVALUE;
- CCTK_REAL PDstandard4th23g12 = INITVALUE;
- CCTK_REAL PDstandard4th31g12 = INITVALUE;
- CCTK_REAL PDstandard4th32g12 = INITVALUE;
- CCTK_REAL PDstandard4th1g13 = INITVALUE;
- CCTK_REAL PDstandard4th2g13 = INITVALUE;
- CCTK_REAL PDstandard4th3g13 = INITVALUE;
- CCTK_REAL PDstandard4th22g13 = INITVALUE;
- CCTK_REAL PDstandard4th12g13 = INITVALUE;
- CCTK_REAL PDstandard4th13g13 = INITVALUE;
- CCTK_REAL PDstandard4th21g13 = INITVALUE;
- CCTK_REAL PDstandard4th23g13 = INITVALUE;
- CCTK_REAL PDstandard4th31g13 = INITVALUE;
- CCTK_REAL PDstandard4th32g13 = INITVALUE;
- CCTK_REAL PDstandard4th1g22 = INITVALUE;
- CCTK_REAL PDstandard4th2g22 = INITVALUE;
- CCTK_REAL PDstandard4th3g22 = INITVALUE;
- CCTK_REAL PDstandard4th11g22 = INITVALUE;
- CCTK_REAL PDstandard4th33g22 = INITVALUE;
- CCTK_REAL PDstandard4th12g22 = INITVALUE;
- CCTK_REAL PDstandard4th13g22 = INITVALUE;
- CCTK_REAL PDstandard4th21g22 = INITVALUE;
- CCTK_REAL PDstandard4th23g22 = INITVALUE;
- CCTK_REAL PDstandard4th31g22 = INITVALUE;
- CCTK_REAL PDstandard4th32g22 = INITVALUE;
- CCTK_REAL PDstandard4th1g23 = INITVALUE;
- CCTK_REAL PDstandard4th2g23 = INITVALUE;
- CCTK_REAL PDstandard4th3g23 = INITVALUE;
- CCTK_REAL PDstandard4th11g23 = INITVALUE;
- CCTK_REAL PDstandard4th12g23 = INITVALUE;
- CCTK_REAL PDstandard4th13g23 = INITVALUE;
- CCTK_REAL PDstandard4th21g23 = INITVALUE;
- CCTK_REAL PDstandard4th23g23 = INITVALUE;
- CCTK_REAL PDstandard4th31g23 = INITVALUE;
- CCTK_REAL PDstandard4th32g23 = INITVALUE;
- CCTK_REAL PDstandard4th1g33 = INITVALUE;
- CCTK_REAL PDstandard4th2g33 = INITVALUE;
- CCTK_REAL PDstandard4th3g33 = INITVALUE;
- CCTK_REAL PDstandard4th11g33 = INITVALUE;
- CCTK_REAL PDstandard4th22g33 = INITVALUE;
- CCTK_REAL PDstandard4th12g33 = INITVALUE;
- CCTK_REAL PDstandard4th13g33 = INITVALUE;
- CCTK_REAL PDstandard4th21g33 = INITVALUE;
- CCTK_REAL PDstandard4th23g33 = INITVALUE;
- CCTK_REAL PDstandard4th31g33 = INITVALUE;
- CCTK_REAL PDstandard4th32g33 = INITVALUE;
- CCTK_REAL PDstandard4th1K11 = INITVALUE;
- CCTK_REAL PDstandard4th2K11 = INITVALUE;
- CCTK_REAL PDstandard4th3K11 = INITVALUE;
- CCTK_REAL PDstandard4th1K12 = INITVALUE;
- CCTK_REAL PDstandard4th2K12 = INITVALUE;
- CCTK_REAL PDstandard4th3K12 = INITVALUE;
- CCTK_REAL PDstandard4th1K13 = INITVALUE;
- CCTK_REAL PDstandard4th2K13 = INITVALUE;
- CCTK_REAL PDstandard4th3K13 = INITVALUE;
- CCTK_REAL PDstandard4th1K22 = INITVALUE;
- CCTK_REAL PDstandard4th2K22 = INITVALUE;
- CCTK_REAL PDstandard4th3K22 = INITVALUE;
- CCTK_REAL PDstandard4th1K23 = INITVALUE;
- CCTK_REAL PDstandard4th2K23 = INITVALUE;
- CCTK_REAL PDstandard4th3K23 = INITVALUE;
- CCTK_REAL PDstandard4th1K33 = INITVALUE;
- CCTK_REAL PDstandard4th2K33 = INITVALUE;
- CCTK_REAL PDstandard4th3K33 = INITVALUE;
+ CCTK_REAL PDstandardNth1alpha = INITVALUE;
+ CCTK_REAL PDstandardNth2alpha = INITVALUE;
+ CCTK_REAL PDstandardNth3alpha = INITVALUE;
+ CCTK_REAL PDstandardNth11alpha = INITVALUE;
+ CCTK_REAL PDstandardNth22alpha = INITVALUE;
+ CCTK_REAL PDstandardNth33alpha = INITVALUE;
+ CCTK_REAL PDstandardNth12alpha = INITVALUE;
+ CCTK_REAL PDstandardNth13alpha = INITVALUE;
+ CCTK_REAL PDstandardNth23alpha = INITVALUE;
+ CCTK_REAL PDstandardNth1beta1 = INITVALUE;
+ CCTK_REAL PDstandardNth2beta1 = INITVALUE;
+ CCTK_REAL PDstandardNth3beta1 = INITVALUE;
+ CCTK_REAL PDstandardNth1beta2 = INITVALUE;
+ CCTK_REAL PDstandardNth2beta2 = INITVALUE;
+ CCTK_REAL PDstandardNth3beta2 = INITVALUE;
+ CCTK_REAL PDstandardNth1beta3 = INITVALUE;
+ CCTK_REAL PDstandardNth2beta3 = INITVALUE;
+ CCTK_REAL PDstandardNth3beta3 = INITVALUE;
+ CCTK_REAL PDstandardNth1g11 = INITVALUE;
+ CCTK_REAL PDstandardNth2g11 = INITVALUE;
+ CCTK_REAL PDstandardNth3g11 = INITVALUE;
+ CCTK_REAL PDstandardNth22g11 = INITVALUE;
+ CCTK_REAL PDstandardNth33g11 = INITVALUE;
+ CCTK_REAL PDstandardNth12g11 = INITVALUE;
+ CCTK_REAL PDstandardNth13g11 = INITVALUE;
+ CCTK_REAL PDstandardNth21g11 = INITVALUE;
+ CCTK_REAL PDstandardNth23g11 = INITVALUE;
+ CCTK_REAL PDstandardNth31g11 = INITVALUE;
+ CCTK_REAL PDstandardNth32g11 = INITVALUE;
+ CCTK_REAL PDstandardNth1g12 = INITVALUE;
+ CCTK_REAL PDstandardNth2g12 = INITVALUE;
+ CCTK_REAL PDstandardNth3g12 = INITVALUE;
+ CCTK_REAL PDstandardNth33g12 = INITVALUE;
+ CCTK_REAL PDstandardNth12g12 = INITVALUE;
+ CCTK_REAL PDstandardNth13g12 = INITVALUE;
+ CCTK_REAL PDstandardNth21g12 = INITVALUE;
+ CCTK_REAL PDstandardNth23g12 = INITVALUE;
+ CCTK_REAL PDstandardNth31g12 = INITVALUE;
+ CCTK_REAL PDstandardNth32g12 = INITVALUE;
+ CCTK_REAL PDstandardNth1g13 = INITVALUE;
+ CCTK_REAL PDstandardNth2g13 = INITVALUE;
+ CCTK_REAL PDstandardNth3g13 = INITVALUE;
+ CCTK_REAL PDstandardNth22g13 = INITVALUE;
+ CCTK_REAL PDstandardNth12g13 = INITVALUE;
+ CCTK_REAL PDstandardNth13g13 = INITVALUE;
+ CCTK_REAL PDstandardNth21g13 = INITVALUE;
+ CCTK_REAL PDstandardNth23g13 = INITVALUE;
+ CCTK_REAL PDstandardNth31g13 = INITVALUE;
+ CCTK_REAL PDstandardNth32g13 = INITVALUE;
+ CCTK_REAL PDstandardNth1g22 = INITVALUE;
+ CCTK_REAL PDstandardNth2g22 = INITVALUE;
+ CCTK_REAL PDstandardNth3g22 = INITVALUE;
+ CCTK_REAL PDstandardNth11g22 = INITVALUE;
+ CCTK_REAL PDstandardNth33g22 = INITVALUE;
+ CCTK_REAL PDstandardNth12g22 = INITVALUE;
+ CCTK_REAL PDstandardNth13g22 = INITVALUE;
+ CCTK_REAL PDstandardNth21g22 = INITVALUE;
+ CCTK_REAL PDstandardNth23g22 = INITVALUE;
+ CCTK_REAL PDstandardNth31g22 = INITVALUE;
+ CCTK_REAL PDstandardNth32g22 = INITVALUE;
+ CCTK_REAL PDstandardNth1g23 = INITVALUE;
+ CCTK_REAL PDstandardNth2g23 = INITVALUE;
+ CCTK_REAL PDstandardNth3g23 = INITVALUE;
+ CCTK_REAL PDstandardNth11g23 = INITVALUE;
+ CCTK_REAL PDstandardNth12g23 = INITVALUE;
+ CCTK_REAL PDstandardNth13g23 = INITVALUE;
+ CCTK_REAL PDstandardNth21g23 = INITVALUE;
+ CCTK_REAL PDstandardNth23g23 = INITVALUE;
+ CCTK_REAL PDstandardNth31g23 = INITVALUE;
+ CCTK_REAL PDstandardNth32g23 = INITVALUE;
+ CCTK_REAL PDstandardNth1g33 = INITVALUE;
+ CCTK_REAL PDstandardNth2g33 = INITVALUE;
+ CCTK_REAL PDstandardNth3g33 = INITVALUE;
+ CCTK_REAL PDstandardNth11g33 = INITVALUE;
+ CCTK_REAL PDstandardNth22g33 = INITVALUE;
+ CCTK_REAL PDstandardNth12g33 = INITVALUE;
+ CCTK_REAL PDstandardNth13g33 = INITVALUE;
+ CCTK_REAL PDstandardNth21g33 = INITVALUE;
+ CCTK_REAL PDstandardNth23g33 = INITVALUE;
+ CCTK_REAL PDstandardNth31g33 = INITVALUE;
+ CCTK_REAL PDstandardNth32g33 = INITVALUE;
+ CCTK_REAL PDstandardNth1K11 = INITVALUE;
+ CCTK_REAL PDstandardNth2K11 = INITVALUE;
+ CCTK_REAL PDstandardNth3K11 = INITVALUE;
+ CCTK_REAL PDstandardNth1K12 = INITVALUE;
+ CCTK_REAL PDstandardNth2K12 = INITVALUE;
+ CCTK_REAL PDstandardNth3K12 = INITVALUE;
+ CCTK_REAL PDstandardNth1K13 = INITVALUE;
+ CCTK_REAL PDstandardNth2K13 = INITVALUE;
+ CCTK_REAL PDstandardNth3K13 = INITVALUE;
+ CCTK_REAL PDstandardNth1K22 = INITVALUE;
+ CCTK_REAL PDstandardNth2K22 = INITVALUE;
+ CCTK_REAL PDstandardNth3K22 = INITVALUE;
+ CCTK_REAL PDstandardNth1K23 = INITVALUE;
+ CCTK_REAL PDstandardNth2K23 = INITVALUE;
+ CCTK_REAL PDstandardNth3K23 = INITVALUE;
+ CCTK_REAL PDstandardNth1K33 = INITVALUE;
+ CCTK_REAL PDstandardNth2K33 = INITVALUE;
+ CCTK_REAL PDstandardNth3K33 = INITVALUE;
/* Assign local copies of grid functions */
alphaL = alpha[index];
@@ -237,81 +237,81 @@ void ML_ADM_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[
/* Include user supplied include files */
/* Precompute derivatives (new style) */
- PDstandard4th1alpha = PDstandard4th1(alpha, i, j, k);
- PDstandard4th2alpha = PDstandard4th2(alpha, i, j, k);
- PDstandard4th3alpha = PDstandard4th3(alpha, i, j, k);
- PDstandard4th11alpha = PDstandard4th11(alpha, i, j, k);
- PDstandard4th22alpha = PDstandard4th22(alpha, i, j, k);
- PDstandard4th33alpha = PDstandard4th33(alpha, i, j, k);
- PDstandard4th12alpha = PDstandard4th12(alpha, i, j, k);
- PDstandard4th13alpha = PDstandard4th13(alpha, i, j, k);
- PDstandard4th23alpha = PDstandard4th23(alpha, i, j, k);
- PDstandard4th1beta1 = PDstandard4th1(beta1, i, j, k);
- PDstandard4th2beta1 = PDstandard4th2(beta1, i, j, k);
- PDstandard4th3beta1 = PDstandard4th3(beta1, i, j, k);
- PDstandard4th1beta2 = PDstandard4th1(beta2, i, j, k);
- PDstandard4th2beta2 = PDstandard4th2(beta2, i, j, k);
- PDstandard4th3beta2 = PDstandard4th3(beta2, i, j, k);
- PDstandard4th1beta3 = PDstandard4th1(beta3, i, j, k);
- PDstandard4th2beta3 = PDstandard4th2(beta3, i, j, k);
- PDstandard4th3beta3 = PDstandard4th3(beta3, i, j, k);
- PDstandard4th1g11 = PDstandard4th1(g11, i, j, k);
- PDstandard4th2g11 = PDstandard4th2(g11, i, j, k);
- PDstandard4th3g11 = PDstandard4th3(g11, i, j, k);
- PDstandard4th22g11 = PDstandard4th22(g11, i, j, k);
- PDstandard4th33g11 = PDstandard4th33(g11, i, j, k);
- PDstandard4th23g11 = PDstandard4th23(g11, i, j, k);
- PDstandard4th1g12 = PDstandard4th1(g12, i, j, k);
- PDstandard4th2g12 = PDstandard4th2(g12, i, j, k);
- PDstandard4th3g12 = PDstandard4th3(g12, i, j, k);
- PDstandard4th33g12 = PDstandard4th33(g12, i, j, k);
- PDstandard4th12g12 = PDstandard4th12(g12, i, j, k);
- PDstandard4th13g12 = PDstandard4th13(g12, i, j, k);
- PDstandard4th23g12 = PDstandard4th23(g12, i, j, k);
- PDstandard4th1g13 = PDstandard4th1(g13, i, j, k);
- PDstandard4th2g13 = PDstandard4th2(g13, i, j, k);
- PDstandard4th3g13 = PDstandard4th3(g13, i, j, k);
- PDstandard4th22g13 = PDstandard4th22(g13, i, j, k);
- PDstandard4th12g13 = PDstandard4th12(g13, i, j, k);
- PDstandard4th13g13 = PDstandard4th13(g13, i, j, k);
- PDstandard4th23g13 = PDstandard4th23(g13, i, j, k);
- PDstandard4th1g22 = PDstandard4th1(g22, i, j, k);
- PDstandard4th2g22 = PDstandard4th2(g22, i, j, k);
- PDstandard4th3g22 = PDstandard4th3(g22, i, j, k);
- PDstandard4th11g22 = PDstandard4th11(g22, i, j, k);
- PDstandard4th33g22 = PDstandard4th33(g22, i, j, k);
- PDstandard4th13g22 = PDstandard4th13(g22, i, j, k);
- PDstandard4th1g23 = PDstandard4th1(g23, i, j, k);
- PDstandard4th2g23 = PDstandard4th2(g23, i, j, k);
- PDstandard4th3g23 = PDstandard4th3(g23, i, j, k);
- PDstandard4th11g23 = PDstandard4th11(g23, i, j, k);
- PDstandard4th12g23 = PDstandard4th12(g23, i, j, k);
- PDstandard4th13g23 = PDstandard4th13(g23, i, j, k);
- PDstandard4th23g23 = PDstandard4th23(g23, i, j, k);
- PDstandard4th1g33 = PDstandard4th1(g33, i, j, k);
- PDstandard4th2g33 = PDstandard4th2(g33, i, j, k);
- PDstandard4th3g33 = PDstandard4th3(g33, i, j, k);
- PDstandard4th11g33 = PDstandard4th11(g33, i, j, k);
- PDstandard4th22g33 = PDstandard4th22(g33, i, j, k);
- PDstandard4th12g33 = PDstandard4th12(g33, i, j, k);
- PDstandard4th1K11 = PDstandard4th1(K11, i, j, k);
- PDstandard4th2K11 = PDstandard4th2(K11, i, j, k);
- PDstandard4th3K11 = PDstandard4th3(K11, i, j, k);
- PDstandard4th1K12 = PDstandard4th1(K12, i, j, k);
- PDstandard4th2K12 = PDstandard4th2(K12, i, j, k);
- PDstandard4th3K12 = PDstandard4th3(K12, i, j, k);
- PDstandard4th1K13 = PDstandard4th1(K13, i, j, k);
- PDstandard4th2K13 = PDstandard4th2(K13, i, j, k);
- PDstandard4th3K13 = PDstandard4th3(K13, i, j, k);
- PDstandard4th1K22 = PDstandard4th1(K22, i, j, k);
- PDstandard4th2K22 = PDstandard4th2(K22, i, j, k);
- PDstandard4th3K22 = PDstandard4th3(K22, i, j, k);
- PDstandard4th1K23 = PDstandard4th1(K23, i, j, k);
- PDstandard4th2K23 = PDstandard4th2(K23, i, j, k);
- PDstandard4th3K23 = PDstandard4th3(K23, i, j, k);
- PDstandard4th1K33 = PDstandard4th1(K33, i, j, k);
- PDstandard4th2K33 = PDstandard4th2(K33, i, j, k);
- PDstandard4th3K33 = PDstandard4th3(K33, i, j, k);
+ PDstandardNth1alpha = PDstandardNth1(alpha, i, j, k);
+ PDstandardNth2alpha = PDstandardNth2(alpha, i, j, k);
+ PDstandardNth3alpha = PDstandardNth3(alpha, i, j, k);
+ PDstandardNth11alpha = PDstandardNth11(alpha, i, j, k);
+ PDstandardNth22alpha = PDstandardNth22(alpha, i, j, k);
+ PDstandardNth33alpha = PDstandardNth33(alpha, i, j, k);
+ PDstandardNth12alpha = PDstandardNth12(alpha, i, j, k);
+ PDstandardNth13alpha = PDstandardNth13(alpha, i, j, k);
+ PDstandardNth23alpha = PDstandardNth23(alpha, i, j, k);
+ PDstandardNth1beta1 = PDstandardNth1(beta1, i, j, k);
+ PDstandardNth2beta1 = PDstandardNth2(beta1, i, j, k);
+ PDstandardNth3beta1 = PDstandardNth3(beta1, i, j, k);
+ PDstandardNth1beta2 = PDstandardNth1(beta2, i, j, k);
+ PDstandardNth2beta2 = PDstandardNth2(beta2, i, j, k);
+ PDstandardNth3beta2 = PDstandardNth3(beta2, i, j, k);
+ PDstandardNth1beta3 = PDstandardNth1(beta3, i, j, k);
+ PDstandardNth2beta3 = PDstandardNth2(beta3, i, j, k);
+ PDstandardNth3beta3 = PDstandardNth3(beta3, i, j, k);
+ PDstandardNth1g11 = PDstandardNth1(g11, i, j, k);
+ PDstandardNth2g11 = PDstandardNth2(g11, i, j, k);
+ PDstandardNth3g11 = PDstandardNth3(g11, i, j, k);
+ PDstandardNth22g11 = PDstandardNth22(g11, i, j, k);
+ PDstandardNth33g11 = PDstandardNth33(g11, i, j, k);
+ PDstandardNth23g11 = PDstandardNth23(g11, i, j, k);
+ PDstandardNth1g12 = PDstandardNth1(g12, i, j, k);
+ PDstandardNth2g12 = PDstandardNth2(g12, i, j, k);
+ PDstandardNth3g12 = PDstandardNth3(g12, i, j, k);
+ PDstandardNth33g12 = PDstandardNth33(g12, i, j, k);
+ PDstandardNth12g12 = PDstandardNth12(g12, i, j, k);
+ PDstandardNth13g12 = PDstandardNth13(g12, i, j, k);
+ PDstandardNth23g12 = PDstandardNth23(g12, i, j, k);
+ PDstandardNth1g13 = PDstandardNth1(g13, i, j, k);
+ PDstandardNth2g13 = PDstandardNth2(g13, i, j, k);
+ PDstandardNth3g13 = PDstandardNth3(g13, i, j, k);
+ PDstandardNth22g13 = PDstandardNth22(g13, i, j, k);
+ PDstandardNth12g13 = PDstandardNth12(g13, i, j, k);
+ PDstandardNth13g13 = PDstandardNth13(g13, i, j, k);
+ PDstandardNth23g13 = PDstandardNth23(g13, i, j, k);
+ PDstandardNth1g22 = PDstandardNth1(g22, i, j, k);
+ PDstandardNth2g22 = PDstandardNth2(g22, i, j, k);
+ PDstandardNth3g22 = PDstandardNth3(g22, i, j, k);
+ PDstandardNth11g22 = PDstandardNth11(g22, i, j, k);
+ PDstandardNth33g22 = PDstandardNth33(g22, i, j, k);
+ PDstandardNth13g22 = PDstandardNth13(g22, i, j, k);
+ PDstandardNth1g23 = PDstandardNth1(g23, i, j, k);
+ PDstandardNth2g23 = PDstandardNth2(g23, i, j, k);
+ PDstandardNth3g23 = PDstandardNth3(g23, i, j, k);
+ PDstandardNth11g23 = PDstandardNth11(g23, i, j, k);
+ PDstandardNth12g23 = PDstandardNth12(g23, i, j, k);
+ PDstandardNth13g23 = PDstandardNth13(g23, i, j, k);
+ PDstandardNth23g23 = PDstandardNth23(g23, i, j, k);
+ PDstandardNth1g33 = PDstandardNth1(g33, i, j, k);
+ PDstandardNth2g33 = PDstandardNth2(g33, i, j, k);
+ PDstandardNth3g33 = PDstandardNth3(g33, i, j, k);
+ PDstandardNth11g33 = PDstandardNth11(g33, i, j, k);
+ PDstandardNth22g33 = PDstandardNth22(g33, i, j, k);
+ PDstandardNth12g33 = PDstandardNth12(g33, i, j, k);
+ PDstandardNth1K11 = PDstandardNth1(K11, i, j, k);
+ PDstandardNth2K11 = PDstandardNth2(K11, i, j, k);
+ PDstandardNth3K11 = PDstandardNth3(K11, i, j, k);
+ PDstandardNth1K12 = PDstandardNth1(K12, i, j, k);
+ PDstandardNth2K12 = PDstandardNth2(K12, i, j, k);
+ PDstandardNth3K12 = PDstandardNth3(K12, i, j, k);
+ PDstandardNth1K13 = PDstandardNth1(K13, i, j, k);
+ PDstandardNth2K13 = PDstandardNth2(K13, i, j, k);
+ PDstandardNth3K13 = PDstandardNth3(K13, i, j, k);
+ PDstandardNth1K22 = PDstandardNth1(K22, i, j, k);
+ PDstandardNth2K22 = PDstandardNth2(K22, i, j, k);
+ PDstandardNth3K22 = PDstandardNth3(K22, i, j, k);
+ PDstandardNth1K23 = PDstandardNth1(K23, i, j, k);
+ PDstandardNth2K23 = PDstandardNth2(K23, i, j, k);
+ PDstandardNth3K23 = PDstandardNth3(K23, i, j, k);
+ PDstandardNth1K33 = PDstandardNth1(K33, i, j, k);
+ PDstandardNth2K33 = PDstandardNth2(K33, i, j, k);
+ PDstandardNth3K33 = PDstandardNth3(K33, i, j, k);
/* Precompute derivatives (old style) */
@@ -330,97 +330,97 @@ void ML_ADM_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[
gu33 = INV(detg)*(g11L*g22L - SQR(g12L));
- G111 = khalf*(gu11*PDstandard4th1g11 + 2*(gu21*PDstandard4th1g12 + gu31*PDstandard4th1g13) - gu21*PDstandard4th2g11 -
- gu31*PDstandard4th3g11);
+ G111 = khalf*(gu11*PDstandardNth1g11 + 2*(gu21*PDstandardNth1g12 + gu31*PDstandardNth1g13) - gu21*PDstandardNth2g11 -
+ gu31*PDstandardNth3g11);
- G211 = khalf*(gu21*PDstandard4th1g11 + 2*(gu22*PDstandard4th1g12 + gu32*PDstandard4th1g13) - gu22*PDstandard4th2g11 -
- gu32*PDstandard4th3g11);
+ G211 = khalf*(gu21*PDstandardNth1g11 + 2*(gu22*PDstandardNth1g12 + gu32*PDstandardNth1g13) - gu22*PDstandardNth2g11 -
+ gu32*PDstandardNth3g11);
- G311 = khalf*(gu31*PDstandard4th1g11 + 2*(gu32*PDstandard4th1g12 + gu33*PDstandard4th1g13) - gu32*PDstandard4th2g11 -
- gu33*PDstandard4th3g11);
+ G311 = khalf*(gu31*PDstandardNth1g11 + 2*(gu32*PDstandardNth1g12 + gu33*PDstandardNth1g13) - gu32*PDstandardNth2g11 -
+ gu33*PDstandardNth3g11);
- G112 = khalf*(gu21*PDstandard4th1g22 + gu11*PDstandard4th2g11 +
- gu31*(PDstandard4th1g23 + PDstandard4th2g13 - PDstandard4th3g12));
+ G112 = khalf*(gu21*PDstandardNth1g22 + gu11*PDstandardNth2g11 +
+ gu31*(PDstandardNth1g23 + PDstandardNth2g13 - PDstandardNth3g12));
- G212 = khalf*(gu22*PDstandard4th1g22 + gu21*PDstandard4th2g11 +
- gu32*(PDstandard4th1g23 + PDstandard4th2g13 - PDstandard4th3g12));
+ G212 = khalf*(gu22*PDstandardNth1g22 + gu21*PDstandardNth2g11 +
+ gu32*(PDstandardNth1g23 + PDstandardNth2g13 - PDstandardNth3g12));
- G312 = khalf*(gu32*PDstandard4th1g22 + gu31*PDstandard4th2g11 +
- gu33*(PDstandard4th1g23 + PDstandard4th2g13 - PDstandard4th3g12));
+ G312 = khalf*(gu32*PDstandardNth1g22 + gu31*PDstandardNth2g11 +
+ gu33*(PDstandardNth1g23 + PDstandardNth2g13 - PDstandardNth3g12));
- G113 = khalf*(gu31*PDstandard4th1g33 + gu11*PDstandard4th3g11 +
- gu21*(PDstandard4th1g23 - PDstandard4th2g13 + PDstandard4th3g12));
+ G113 = khalf*(gu31*PDstandardNth1g33 + gu11*PDstandardNth3g11 +
+ gu21*(PDstandardNth1g23 - PDstandardNth2g13 + PDstandardNth3g12));
- G213 = khalf*(gu32*PDstandard4th1g33 + gu21*PDstandard4th3g11 +
- gu22*(PDstandard4th1g23 - PDstandard4th2g13 + PDstandard4th3g12));
+ G213 = khalf*(gu32*PDstandardNth1g33 + gu21*PDstandardNth3g11 +
+ gu22*(PDstandardNth1g23 - PDstandardNth2g13 + PDstandardNth3g12));
- G313 = khalf*(gu33*PDstandard4th1g33 + gu31*PDstandard4th3g11 +
- gu32*(PDstandard4th1g23 - PDstandard4th2g13 + PDstandard4th3g12));
+ G313 = khalf*(gu33*PDstandardNth1g33 + gu31*PDstandardNth3g11 +
+ gu32*(PDstandardNth1g23 - PDstandardNth2g13 + PDstandardNth3g12));
- G122 = khalf*(gu11*(-PDstandard4th1g22 + 2*PDstandard4th2g12) + gu21*PDstandard4th2g22 +
- gu31*(2*PDstandard4th2g23 - PDstandard4th3g22));
+ G122 = khalf*(gu11*(-PDstandardNth1g22 + 2*PDstandardNth2g12) + gu21*PDstandardNth2g22 +
+ gu31*(2*PDstandardNth2g23 - PDstandardNth3g22));
- G222 = khalf*(gu21*(-PDstandard4th1g22 + 2*PDstandard4th2g12) + gu22*PDstandard4th2g22 +
- gu32*(2*PDstandard4th2g23 - PDstandard4th3g22));
+ G222 = khalf*(gu21*(-PDstandardNth1g22 + 2*PDstandardNth2g12) + gu22*PDstandardNth2g22 +
+ gu32*(2*PDstandardNth2g23 - PDstandardNth3g22));
- G322 = khalf*(gu31*(-PDstandard4th1g22 + 2*PDstandard4th2g12) + gu32*PDstandard4th2g22 +
- gu33*(2*PDstandard4th2g23 - PDstandard4th3g22));
+ G322 = khalf*(gu31*(-PDstandardNth1g22 + 2*PDstandardNth2g12) + gu32*PDstandardNth2g22 +
+ gu33*(2*PDstandardNth2g23 - PDstandardNth3g22));
- G123 = khalf*(gu31*PDstandard4th2g33 + gu11*(-PDstandard4th1g23 + PDstandard4th2g13 + PDstandard4th3g12) +
- gu21*PDstandard4th3g22);
+ G123 = khalf*(gu31*PDstandardNth2g33 + gu11*(-PDstandardNth1g23 + PDstandardNth2g13 + PDstandardNth3g12) +
+ gu21*PDstandardNth3g22);
- G223 = khalf*(gu32*PDstandard4th2g33 + gu21*(-PDstandard4th1g23 + PDstandard4th2g13 + PDstandard4th3g12) +
- gu22*PDstandard4th3g22);
+ G223 = khalf*(gu32*PDstandardNth2g33 + gu21*(-PDstandardNth1g23 + PDstandardNth2g13 + PDstandardNth3g12) +
+ gu22*PDstandardNth3g22);
- G323 = khalf*(gu33*PDstandard4th2g33 + gu31*(-PDstandard4th1g23 + PDstandard4th2g13 + PDstandard4th3g12) +
- gu32*PDstandard4th3g22);
+ G323 = khalf*(gu33*PDstandardNth2g33 + gu31*(-PDstandardNth1g23 + PDstandardNth2g13 + PDstandardNth3g12) +
+ gu32*PDstandardNth3g22);
- G133 = khalf*(-(gu11*PDstandard4th1g33) - gu21*PDstandard4th2g33 + 2*gu11*PDstandard4th3g13 +
- 2*gu21*PDstandard4th3g23 + gu31*PDstandard4th3g33);
+ G133 = khalf*(-(gu11*PDstandardNth1g33) - gu21*PDstandardNth2g33 + 2*gu11*PDstandardNth3g13 +
+ 2*gu21*PDstandardNth3g23 + gu31*PDstandardNth3g33);
- G233 = khalf*(-(gu21*PDstandard4th1g33) - gu22*PDstandard4th2g33 + 2*gu21*PDstandard4th3g13 +
- 2*gu22*PDstandard4th3g23 + gu32*PDstandard4th3g33);
+ G233 = khalf*(-(gu21*PDstandardNth1g33) - gu22*PDstandardNth2g33 + 2*gu21*PDstandardNth3g13 +
+ 2*gu22*PDstandardNth3g23 + gu32*PDstandardNth3g33);
- G333 = khalf*(-(gu31*PDstandard4th1g33) - gu32*PDstandard4th2g33 + 2*gu31*PDstandard4th3g13 +
- 2*gu32*PDstandard4th3g23 + gu33*PDstandard4th3g33);
+ G333 = khalf*(-(gu31*PDstandardNth1g33) - gu32*PDstandardNth2g33 + 2*gu31*PDstandardNth3g13 +
+ 2*gu32*PDstandardNth3g23 + gu33*PDstandardNth3g33);
- R11 = khalf*(4*G213*G312 + G211*(2*G112 - 2*G222 - 2*G323) + G311*(2*G113 - 2*G333) - gu22*PDstandard4th11g22 -
- 2*(G111*G212 + G223*G311 + G111*G313 + gu32*PDstandard4th11g23) - gu33*PDstandard4th11g33 +
- 2*gu22*PDstandard4th12g12 + 2*gu32*PDstandard4th12g13 + 2*gu32*PDstandard4th13g12 + 2*gu33*PDstandard4th13g13 -
- gu22*PDstandard4th22g11 - 2*gu32*PDstandard4th23g11 - gu33*PDstandard4th33g11 + 2*SQR(G212) + 2*SQR(G313));
+ R11 = khalf*(4*G213*G312 + G211*(2*G112 - 2*G222 - 2*G323) + G311*(2*G113 - 2*G333) - gu22*PDstandardNth11g22 -
+ 2*(G111*G212 + G223*G311 + G111*G313 + gu32*PDstandardNth11g23) - gu33*PDstandardNth11g33 +
+ 2*gu22*PDstandardNth12g12 + 2*gu32*PDstandardNth12g13 + 2*gu32*PDstandardNth13g12 + 2*gu33*PDstandardNth13g13 -
+ gu22*PDstandardNth22g11 - 2*gu32*PDstandardNth23g11 - gu33*PDstandardNth33g11 + 2*SQR(G212) + 2*SQR(G313));
R12 = khalf*(2*(G122*G211 + G123*G311 + G213*G322 + G313*G323) -
- 2*(G112*G212 + G112*G313 + G212*G323 + G312*G333 + gu21*PDstandard4th12g12) - gu32*PDstandard4th12g23 -
- gu33*PDstandard4th12g33 + gu31*(PDstandard4th11g23 - PDstandard4th12g13 - PDstandard4th13g12) +
- gu32*PDstandard4th13g22 + gu33*PDstandard4th13g23 + gu21*(PDstandard4th11g22 + PDstandard4th22g11) +
- gu32*PDstandard4th22g13 + gu31*PDstandard4th23g11 - gu32*PDstandard4th23g12 + gu33*PDstandard4th23g13 -
- gu33*PDstandard4th33g12);
+ 2*(G112*G212 + G112*G313 + G212*G323 + G312*G333 + gu21*PDstandardNth12g12) - gu32*PDstandardNth12g23 -
+ gu33*PDstandardNth12g33 + gu31*(PDstandardNth11g23 - PDstandardNth12g13 - PDstandardNth13g12) +
+ gu32*PDstandardNth13g22 + gu33*PDstandardNth13g23 + gu21*(PDstandardNth11g22 + PDstandardNth22g11) +
+ gu32*PDstandardNth22g13 + gu31*PDstandardNth23g11 - gu32*PDstandardNth23g12 + gu33*PDstandardNth23g13 -
+ gu33*PDstandardNth33g12);
R13 = khalf*(2*(G123*G211 + G212*G223 + G133*G311 + G233*G312) -
- 2*(G213*G222 + G223*G313 + G113*(G212 + G313) + gu31*PDstandard4th13g13) +
- gu21*(PDstandard4th11g23 - PDstandard4th12g13 - PDstandard4th13g12 + PDstandard4th23g11) +
- gu22*(PDstandard4th12g23 - PDstandard4th13g22 - PDstandard4th22g13 + PDstandard4th23g12) +
- gu31*(PDstandard4th11g33 + PDstandard4th33g11) +
- gu32*(PDstandard4th12g33 - PDstandard4th13g23 - PDstandard4th23g13 + PDstandard4th33g12));
+ 2*(G213*G222 + G223*G313 + G113*(G212 + G313) + gu31*PDstandardNth13g13) +
+ gu21*(PDstandardNth11g23 - PDstandardNth12g13 - PDstandardNth13g12 + PDstandardNth23g11) +
+ gu22*(PDstandardNth12g23 - PDstandardNth13g22 - PDstandardNth22g13 + PDstandardNth23g12) +
+ gu31*(PDstandardNth11g33 + PDstandardNth33g11) +
+ gu32*(PDstandardNth12g33 - PDstandardNth13g23 - PDstandardNth23g13 + PDstandardNth33g12));
R22 = khalf*(4*G123*G312 + G122*(2*G212 - 2*(G111 + G313)) + G322*(2*G223 - 2*G333) -
- 2*(G113*G322 + G222*(G112 + G323) + gu31*PDstandard4th13g22) +
- gu11*(-PDstandard4th11g22 + 2*PDstandard4th12g12 - PDstandard4th22g11) +
- gu31*(-2*PDstandard4th22g13 + 2*(PDstandard4th12g23 + PDstandard4th23g12)) +
- gu33*(-PDstandard4th22g33 + 2*PDstandard4th23g23 - PDstandard4th33g22) + 2*(SQR(G112) + SQR(G323)));
+ 2*(G113*G322 + G222*(G112 + G323) + gu31*PDstandardNth13g22) +
+ gu11*(-PDstandardNth11g22 + 2*PDstandardNth12g12 - PDstandardNth22g11) +
+ gu31*(-2*PDstandardNth22g13 + 2*(PDstandardNth12g23 + PDstandardNth23g12)) +
+ gu33*(-PDstandardNth22g33 + 2*PDstandardNth23g23 - PDstandardNth33g22) + 2*(SQR(G112) + SQR(G323)));
R23 = khalf*(2*(G112*G113 + G122*G213 + G133*G312 + G233*G322) +
- gu11*(-PDstandard4th11g23 + PDstandard4th12g13 + PDstandard4th13g12 - PDstandard4th23g11) +
- gu21*(-PDstandard4th12g23 + PDstandard4th13g22 + PDstandard4th22g13 - PDstandard4th23g12) -
- 2*(G111*G123 + G113*G323 + G223*(G112 + G323) + gu32*PDstandard4th23g23) +
- gu31*(PDstandard4th12g33 - PDstandard4th13g23 - PDstandard4th23g13 + PDstandard4th33g12) +
- gu32*(PDstandard4th22g33 + PDstandard4th33g22));
-
- R33 = khalf*(4*G123*G213 - gu11*PDstandard4th11g33 -
- 2*(G111*G133 + G133*G212 + G112*G233 + G222*G233 + G113*G333 + G223*G333 + gu21*PDstandard4th12g33) +
- 2*(G133*G313 + G233*G323 + gu11*PDstandard4th13g13) + 2*gu21*PDstandard4th13g23 - gu22*PDstandard4th22g33 +
- 2*gu21*PDstandard4th23g13 + 2*gu22*PDstandard4th23g23 - gu11*PDstandard4th33g11 - 2*gu21*PDstandard4th33g12 -
- gu22*PDstandard4th33g22 + 2*SQR(G113) + 2*SQR(G223));
+ gu11*(-PDstandardNth11g23 + PDstandardNth12g13 + PDstandardNth13g12 - PDstandardNth23g11) +
+ gu21*(-PDstandardNth12g23 + PDstandardNth13g22 + PDstandardNth22g13 - PDstandardNth23g12) -
+ 2*(G111*G123 + G113*G323 + G223*(G112 + G323) + gu32*PDstandardNth23g23) +
+ gu31*(PDstandardNth12g33 - PDstandardNth13g23 - PDstandardNth23g13 + PDstandardNth33g12) +
+ gu32*(PDstandardNth22g33 + PDstandardNth33g22));
+
+ R33 = khalf*(4*G123*G213 - gu11*PDstandardNth11g33 -
+ 2*(G111*G133 + G133*G212 + G112*G233 + G222*G233 + G113*G333 + G223*G333 + gu21*PDstandardNth12g33) +
+ 2*(G133*G313 + G233*G323 + gu11*PDstandardNth13g13) + 2*gu21*PDstandardNth13g23 - gu22*PDstandardNth22g33 +
+ 2*gu21*PDstandardNth23g13 + 2*gu22*PDstandardNth23g23 - gu11*PDstandardNth33g11 - 2*gu21*PDstandardNth33g12 -
+ gu22*PDstandardNth33g22 + 2*SQR(G113) + 2*SQR(G223));
Km11 = gu11*K11L + gu21*K12L + gu31*K13L;
@@ -442,57 +442,57 @@ void ML_ADM_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[
trK = Km11 + Km22 + Km33;
- g11rhsL = -2*alphaL*K11L + 2*(g11L*PDstandard4th1beta1 + g12L*PDstandard4th1beta2 + g13L*PDstandard4th1beta3) +
- beta1L*PDstandard4th1g11 + beta2L*PDstandard4th2g11 + beta3L*PDstandard4th3g11;
+ g11rhsL = -2*alphaL*K11L + 2*(g11L*PDstandardNth1beta1 + g12L*PDstandardNth1beta2 + g13L*PDstandardNth1beta3) +
+ beta1L*PDstandardNth1g11 + beta2L*PDstandardNth2g11 + beta3L*PDstandardNth3g11;
- g12rhsL = -2*alphaL*K12L + g22L*PDstandard4th1beta2 + g23L*PDstandard4th1beta3 + beta1L*PDstandard4th1g12 +
- g11L*PDstandard4th2beta1 + g12L*(PDstandard4th1beta1 + PDstandard4th2beta2) + g13L*PDstandard4th2beta3 +
- beta2L*PDstandard4th2g12 + beta3L*PDstandard4th3g12;
+ g12rhsL = -2*alphaL*K12L + g22L*PDstandardNth1beta2 + g23L*PDstandardNth1beta3 + beta1L*PDstandardNth1g12 +
+ g11L*PDstandardNth2beta1 + g12L*(PDstandardNth1beta1 + PDstandardNth2beta2) + g13L*PDstandardNth2beta3 +
+ beta2L*PDstandardNth2g12 + beta3L*PDstandardNth3g12;
- g13rhsL = -2*alphaL*K13L + g23L*PDstandard4th1beta2 + g33L*PDstandard4th1beta3 + beta1L*PDstandard4th1g13 +
- beta2L*PDstandard4th2g13 + g11L*PDstandard4th3beta1 + g12L*PDstandard4th3beta2 +
- g13L*(PDstandard4th1beta1 + PDstandard4th3beta3) + beta3L*PDstandard4th3g13;
+ g13rhsL = -2*alphaL*K13L + g23L*PDstandardNth1beta2 + g33L*PDstandardNth1beta3 + beta1L*PDstandardNth1g13 +
+ beta2L*PDstandardNth2g13 + g11L*PDstandardNth3beta1 + g12L*PDstandardNth3beta2 +
+ g13L*(PDstandardNth1beta1 + PDstandardNth3beta3) + beta3L*PDstandardNth3g13;
- g22rhsL = -2*alphaL*K22L + beta1L*PDstandard4th1g22 +
- 2*(g12L*PDstandard4th2beta1 + g22L*PDstandard4th2beta2 + g23L*PDstandard4th2beta3) + beta2L*PDstandard4th2g22 +
- beta3L*PDstandard4th3g22;
+ g22rhsL = -2*alphaL*K22L + beta1L*PDstandardNth1g22 +
+ 2*(g12L*PDstandardNth2beta1 + g22L*PDstandardNth2beta2 + g23L*PDstandardNth2beta3) + beta2L*PDstandardNth2g22 +
+ beta3L*PDstandardNth3g22;
- g23rhsL = -2*alphaL*K23L + beta1L*PDstandard4th1g23 + g13L*PDstandard4th2beta1 + g33L*PDstandard4th2beta3 +
- beta2L*PDstandard4th2g23 + g12L*PDstandard4th3beta1 + g22L*PDstandard4th3beta2 +
- g23L*(PDstandard4th2beta2 + PDstandard4th3beta3) + beta3L*PDstandard4th3g23;
+ g23rhsL = -2*alphaL*K23L + beta1L*PDstandardNth1g23 + g13L*PDstandardNth2beta1 + g33L*PDstandardNth2beta3 +
+ beta2L*PDstandardNth2g23 + g12L*PDstandardNth3beta1 + g22L*PDstandardNth3beta2 +
+ g23L*(PDstandardNth2beta2 + PDstandardNth3beta3) + beta3L*PDstandardNth3g23;
- g33rhsL = -2*alphaL*K33L + beta1L*PDstandard4th1g33 + beta2L*PDstandard4th2g33 +
- 2*(g13L*PDstandard4th3beta1 + g23L*PDstandard4th3beta2 + g33L*PDstandard4th3beta3) + beta3L*PDstandard4th3g33;
+ g33rhsL = -2*alphaL*K33L + beta1L*PDstandardNth1g33 + beta2L*PDstandardNth2g33 +
+ 2*(g13L*PDstandardNth3beta1 + g23L*PDstandardNth3beta2 + g33L*PDstandardNth3beta3) + beta3L*PDstandardNth3g33;
- K11rhsL = -PDstandard4th11alpha + G111*PDstandard4th1alpha +
- 2*(K11L*PDstandard4th1beta1 + K12L*PDstandard4th1beta2 + K13L*PDstandard4th1beta3) + beta1L*PDstandard4th1K11 +
- G211*PDstandard4th2alpha + beta2L*PDstandard4th2K11 + G311*PDstandard4th3alpha + beta3L*PDstandard4th3K11 +
+ K11rhsL = -PDstandardNth11alpha + G111*PDstandardNth1alpha +
+ 2*(K11L*PDstandardNth1beta1 + K12L*PDstandardNth1beta2 + K13L*PDstandardNth1beta3) + beta1L*PDstandardNth1K11 +
+ G211*PDstandardNth2alpha + beta2L*PDstandardNth2K11 + G311*PDstandardNth3alpha + beta3L*PDstandardNth3K11 +
alphaL*(-2*(K11L*Km11 + K12L*Km21 + K13L*Km31) + R11 + K11L*trK);
- K12rhsL = -PDstandard4th12alpha + G112*PDstandard4th1alpha + K22L*PDstandard4th1beta2 + K23L*PDstandard4th1beta3 +
- beta1L*PDstandard4th1K12 + G212*PDstandard4th2alpha + K11L*PDstandard4th2beta1 +
- K12L*(PDstandard4th1beta1 + PDstandard4th2beta2) + K13L*PDstandard4th2beta3 + beta2L*PDstandard4th2K12 +
- G312*PDstandard4th3alpha + beta3L*PDstandard4th3K12 +
+ K12rhsL = -PDstandardNth12alpha + G112*PDstandardNth1alpha + K22L*PDstandardNth1beta2 + K23L*PDstandardNth1beta3 +
+ beta1L*PDstandardNth1K12 + G212*PDstandardNth2alpha + K11L*PDstandardNth2beta1 +
+ K12L*(PDstandardNth1beta1 + PDstandardNth2beta2) + K13L*PDstandardNth2beta3 + beta2L*PDstandardNth2K12 +
+ G312*PDstandardNth3alpha + beta3L*PDstandardNth3K12 +
alphaL*(-2*(K11L*Km12 + K12L*Km22 + K13L*Km32) + R12 + K12L*trK);
- K13rhsL = -PDstandard4th13alpha + G113*PDstandard4th1alpha + K23L*PDstandard4th1beta2 + K33L*PDstandard4th1beta3 +
- beta1L*PDstandard4th1K13 + G213*PDstandard4th2alpha + beta2L*PDstandard4th2K13 + G313*PDstandard4th3alpha +
- K11L*PDstandard4th3beta1 + K12L*PDstandard4th3beta2 + K13L*(PDstandard4th1beta1 + PDstandard4th3beta3) +
- beta3L*PDstandard4th3K13 + alphaL*(-2*(K11L*Km13 + K12L*Km23 + K13L*Km33) + R13 + K13L*trK);
+ K13rhsL = -PDstandardNth13alpha + G113*PDstandardNth1alpha + K23L*PDstandardNth1beta2 + K33L*PDstandardNth1beta3 +
+ beta1L*PDstandardNth1K13 + G213*PDstandardNth2alpha + beta2L*PDstandardNth2K13 + G313*PDstandardNth3alpha +
+ K11L*PDstandardNth3beta1 + K12L*PDstandardNth3beta2 + K13L*(PDstandardNth1beta1 + PDstandardNth3beta3) +
+ beta3L*PDstandardNth3K13 + alphaL*(-2*(K11L*Km13 + K12L*Km23 + K13L*Km33) + R13 + K13L*trK);
- K22rhsL = G122*PDstandard4th1alpha + beta1L*PDstandard4th1K22 - PDstandard4th22alpha + G222*PDstandard4th2alpha +
- 2*(K12L*PDstandard4th2beta1 + K22L*PDstandard4th2beta2 + K23L*PDstandard4th2beta3) + beta2L*PDstandard4th2K22 +
- G322*PDstandard4th3alpha + beta3L*PDstandard4th3K22 +
+ K22rhsL = G122*PDstandardNth1alpha + beta1L*PDstandardNth1K22 - PDstandardNth22alpha + G222*PDstandardNth2alpha +
+ 2*(K12L*PDstandardNth2beta1 + K22L*PDstandardNth2beta2 + K23L*PDstandardNth2beta3) + beta2L*PDstandardNth2K22 +
+ G322*PDstandardNth3alpha + beta3L*PDstandardNth3K22 +
alphaL*(-2*(K12L*Km12 + K22L*Km22 + K23L*Km32) + R22 + K22L*trK);
- K23rhsL = G123*PDstandard4th1alpha + beta1L*PDstandard4th1K23 - PDstandard4th23alpha + G223*PDstandard4th2alpha +
- K13L*PDstandard4th2beta1 + K33L*PDstandard4th2beta3 + beta2L*PDstandard4th2K23 + G323*PDstandard4th3alpha +
- K12L*PDstandard4th3beta1 + K22L*PDstandard4th3beta2 + K23L*(PDstandard4th2beta2 + PDstandard4th3beta3) +
- beta3L*PDstandard4th3K23 + alphaL*(-2*(K12L*Km13 + K22L*Km23 + K23L*Km33) + R23 + K23L*trK);
+ K23rhsL = G123*PDstandardNth1alpha + beta1L*PDstandardNth1K23 - PDstandardNth23alpha + G223*PDstandardNth2alpha +
+ K13L*PDstandardNth2beta1 + K33L*PDstandardNth2beta3 + beta2L*PDstandardNth2K23 + G323*PDstandardNth3alpha +
+ K12L*PDstandardNth3beta1 + K22L*PDstandardNth3beta2 + K23L*(PDstandardNth2beta2 + PDstandardNth3beta3) +
+ beta3L*PDstandardNth3K23 + alphaL*(-2*(K12L*Km13 + K22L*Km23 + K23L*Km33) + R23 + K23L*trK);
- K33rhsL = G133*PDstandard4th1alpha + beta1L*PDstandard4th1K33 + G233*PDstandard4th2alpha + beta2L*PDstandard4th2K33 -
- PDstandard4th33alpha + G333*PDstandard4th3alpha +
- 2*(K13L*PDstandard4th3beta1 + K23L*PDstandard4th3beta2 + K33L*PDstandard4th3beta3) + beta3L*PDstandard4th3K33 +
+ K33rhsL = G133*PDstandardNth1alpha + beta1L*PDstandardNth1K33 + G233*PDstandardNth2alpha + beta2L*PDstandardNth2K33 -
+ PDstandardNth33alpha + G333*PDstandardNth3alpha +
+ 2*(K13L*PDstandardNth3beta1 + K23L*PDstandardNth3beta2 + K33L*PDstandardNth3beta3) + beta3L*PDstandardNth3K33 +
alphaL*(-2*(K13L*Km13 + K23L*Km23 + K33L*Km33) + R33 + K33L*trK);
alpharhsL = 0;
diff --git a/ML_ADM/src/ML_ADM_constraints.c b/ML_ADM/src/ML_ADM_constraints.c
index f373870..389568d 100644
--- a/ML_ADM/src/ML_ADM_constraints.c
+++ b/ML_ADM/src/ML_ADM_constraints.c
@@ -113,84 +113,84 @@ void ML_ADM_constraints_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL
/* Declare precomputed derivatives*/
/* Declare derivatives */
- CCTK_REAL PDstandard4th1g11 = INITVALUE;
- CCTK_REAL PDstandard4th2g11 = INITVALUE;
- CCTK_REAL PDstandard4th3g11 = INITVALUE;
- CCTK_REAL PDstandard4th22g11 = INITVALUE;
- CCTK_REAL PDstandard4th33g11 = INITVALUE;
- CCTK_REAL PDstandard4th12g11 = INITVALUE;
- CCTK_REAL PDstandard4th13g11 = INITVALUE;
- CCTK_REAL PDstandard4th21g11 = INITVALUE;
- CCTK_REAL PDstandard4th23g11 = INITVALUE;
- CCTK_REAL PDstandard4th31g11 = INITVALUE;
- CCTK_REAL PDstandard4th32g11 = INITVALUE;
- CCTK_REAL PDstandard4th1g12 = INITVALUE;
- CCTK_REAL PDstandard4th2g12 = INITVALUE;
- CCTK_REAL PDstandard4th3g12 = INITVALUE;
- CCTK_REAL PDstandard4th33g12 = INITVALUE;
- CCTK_REAL PDstandard4th12g12 = INITVALUE;
- CCTK_REAL PDstandard4th13g12 = INITVALUE;
- CCTK_REAL PDstandard4th21g12 = INITVALUE;
- CCTK_REAL PDstandard4th23g12 = INITVALUE;
- CCTK_REAL PDstandard4th31g12 = INITVALUE;
- CCTK_REAL PDstandard4th32g12 = INITVALUE;
- CCTK_REAL PDstandard4th1g13 = INITVALUE;
- CCTK_REAL PDstandard4th2g13 = INITVALUE;
- CCTK_REAL PDstandard4th3g13 = INITVALUE;
- CCTK_REAL PDstandard4th22g13 = INITVALUE;
- CCTK_REAL PDstandard4th12g13 = INITVALUE;
- CCTK_REAL PDstandard4th13g13 = INITVALUE;
- CCTK_REAL PDstandard4th21g13 = INITVALUE;
- CCTK_REAL PDstandard4th23g13 = INITVALUE;
- CCTK_REAL PDstandard4th31g13 = INITVALUE;
- CCTK_REAL PDstandard4th32g13 = INITVALUE;
- CCTK_REAL PDstandard4th1g22 = INITVALUE;
- CCTK_REAL PDstandard4th2g22 = INITVALUE;
- CCTK_REAL PDstandard4th3g22 = INITVALUE;
- CCTK_REAL PDstandard4th11g22 = INITVALUE;
- CCTK_REAL PDstandard4th33g22 = INITVALUE;
- CCTK_REAL PDstandard4th12g22 = INITVALUE;
- CCTK_REAL PDstandard4th13g22 = INITVALUE;
- CCTK_REAL PDstandard4th21g22 = INITVALUE;
- CCTK_REAL PDstandard4th23g22 = INITVALUE;
- CCTK_REAL PDstandard4th31g22 = INITVALUE;
- CCTK_REAL PDstandard4th32g22 = INITVALUE;
- CCTK_REAL PDstandard4th1g23 = INITVALUE;
- CCTK_REAL PDstandard4th2g23 = INITVALUE;
- CCTK_REAL PDstandard4th3g23 = INITVALUE;
- CCTK_REAL PDstandard4th11g23 = INITVALUE;
- CCTK_REAL PDstandard4th12g23 = INITVALUE;
- CCTK_REAL PDstandard4th13g23 = INITVALUE;
- CCTK_REAL PDstandard4th21g23 = INITVALUE;
- CCTK_REAL PDstandard4th23g23 = INITVALUE;
- CCTK_REAL PDstandard4th31g23 = INITVALUE;
- CCTK_REAL PDstandard4th32g23 = INITVALUE;
- CCTK_REAL PDstandard4th1g33 = INITVALUE;
- CCTK_REAL PDstandard4th2g33 = INITVALUE;
- CCTK_REAL PDstandard4th3g33 = INITVALUE;
- CCTK_REAL PDstandard4th11g33 = INITVALUE;
- CCTK_REAL PDstandard4th22g33 = INITVALUE;
- CCTK_REAL PDstandard4th12g33 = INITVALUE;
- CCTK_REAL PDstandard4th13g33 = INITVALUE;
- CCTK_REAL PDstandard4th21g33 = INITVALUE;
- CCTK_REAL PDstandard4th23g33 = INITVALUE;
- CCTK_REAL PDstandard4th31g33 = INITVALUE;
- CCTK_REAL PDstandard4th32g33 = INITVALUE;
- CCTK_REAL PDstandard4th2K11 = INITVALUE;
- CCTK_REAL PDstandard4th3K11 = INITVALUE;
- CCTK_REAL PDstandard4th1K12 = INITVALUE;
- CCTK_REAL PDstandard4th2K12 = INITVALUE;
- CCTK_REAL PDstandard4th3K12 = INITVALUE;
- CCTK_REAL PDstandard4th1K13 = INITVALUE;
- CCTK_REAL PDstandard4th2K13 = INITVALUE;
- CCTK_REAL PDstandard4th3K13 = INITVALUE;
- CCTK_REAL PDstandard4th1K22 = INITVALUE;
- CCTK_REAL PDstandard4th3K22 = INITVALUE;
- CCTK_REAL PDstandard4th1K23 = INITVALUE;
- CCTK_REAL PDstandard4th2K23 = INITVALUE;
- CCTK_REAL PDstandard4th3K23 = INITVALUE;
- CCTK_REAL PDstandard4th1K33 = INITVALUE;
- CCTK_REAL PDstandard4th2K33 = INITVALUE;
+ CCTK_REAL PDstandardNth1g11 = INITVALUE;
+ CCTK_REAL PDstandardNth2g11 = INITVALUE;
+ CCTK_REAL PDstandardNth3g11 = INITVALUE;
+ CCTK_REAL PDstandardNth22g11 = INITVALUE;
+ CCTK_REAL PDstandardNth33g11 = INITVALUE;
+ CCTK_REAL PDstandardNth12g11 = INITVALUE;
+ CCTK_REAL PDstandardNth13g11 = INITVALUE;
+ CCTK_REAL PDstandardNth21g11 = INITVALUE;
+ CCTK_REAL PDstandardNth23g11 = INITVALUE;
+ CCTK_REAL PDstandardNth31g11 = INITVALUE;
+ CCTK_REAL PDstandardNth32g11 = INITVALUE;
+ CCTK_REAL PDstandardNth1g12 = INITVALUE;
+ CCTK_REAL PDstandardNth2g12 = INITVALUE;
+ CCTK_REAL PDstandardNth3g12 = INITVALUE;
+ CCTK_REAL PDstandardNth33g12 = INITVALUE;
+ CCTK_REAL PDstandardNth12g12 = INITVALUE;
+ CCTK_REAL PDstandardNth13g12 = INITVALUE;
+ CCTK_REAL PDstandardNth21g12 = INITVALUE;
+ CCTK_REAL PDstandardNth23g12 = INITVALUE;
+ CCTK_REAL PDstandardNth31g12 = INITVALUE;
+ CCTK_REAL PDstandardNth32g12 = INITVALUE;
+ CCTK_REAL PDstandardNth1g13 = INITVALUE;
+ CCTK_REAL PDstandardNth2g13 = INITVALUE;
+ CCTK_REAL PDstandardNth3g13 = INITVALUE;
+ CCTK_REAL PDstandardNth22g13 = INITVALUE;
+ CCTK_REAL PDstandardNth12g13 = INITVALUE;
+ CCTK_REAL PDstandardNth13g13 = INITVALUE;
+ CCTK_REAL PDstandardNth21g13 = INITVALUE;
+ CCTK_REAL PDstandardNth23g13 = INITVALUE;
+ CCTK_REAL PDstandardNth31g13 = INITVALUE;
+ CCTK_REAL PDstandardNth32g13 = INITVALUE;
+ CCTK_REAL PDstandardNth1g22 = INITVALUE;
+ CCTK_REAL PDstandardNth2g22 = INITVALUE;
+ CCTK_REAL PDstandardNth3g22 = INITVALUE;
+ CCTK_REAL PDstandardNth11g22 = INITVALUE;
+ CCTK_REAL PDstandardNth33g22 = INITVALUE;
+ CCTK_REAL PDstandardNth12g22 = INITVALUE;
+ CCTK_REAL PDstandardNth13g22 = INITVALUE;
+ CCTK_REAL PDstandardNth21g22 = INITVALUE;
+ CCTK_REAL PDstandardNth23g22 = INITVALUE;
+ CCTK_REAL PDstandardNth31g22 = INITVALUE;
+ CCTK_REAL PDstandardNth32g22 = INITVALUE;
+ CCTK_REAL PDstandardNth1g23 = INITVALUE;
+ CCTK_REAL PDstandardNth2g23 = INITVALUE;
+ CCTK_REAL PDstandardNth3g23 = INITVALUE;
+ CCTK_REAL PDstandardNth11g23 = INITVALUE;
+ CCTK_REAL PDstandardNth12g23 = INITVALUE;
+ CCTK_REAL PDstandardNth13g23 = INITVALUE;
+ CCTK_REAL PDstandardNth21g23 = INITVALUE;
+ CCTK_REAL PDstandardNth23g23 = INITVALUE;
+ CCTK_REAL PDstandardNth31g23 = INITVALUE;
+ CCTK_REAL PDstandardNth32g23 = INITVALUE;
+ CCTK_REAL PDstandardNth1g33 = INITVALUE;
+ CCTK_REAL PDstandardNth2g33 = INITVALUE;
+ CCTK_REAL PDstandardNth3g33 = INITVALUE;
+ CCTK_REAL PDstandardNth11g33 = INITVALUE;
+ CCTK_REAL PDstandardNth22g33 = INITVALUE;
+ CCTK_REAL PDstandardNth12g33 = INITVALUE;
+ CCTK_REAL PDstandardNth13g33 = INITVALUE;
+ CCTK_REAL PDstandardNth21g33 = INITVALUE;
+ CCTK_REAL PDstandardNth23g33 = INITVALUE;
+ CCTK_REAL PDstandardNth31g33 = INITVALUE;
+ CCTK_REAL PDstandardNth32g33 = INITVALUE;
+ CCTK_REAL PDstandardNth2K11 = INITVALUE;
+ CCTK_REAL PDstandardNth3K11 = INITVALUE;
+ CCTK_REAL PDstandardNth1K12 = INITVALUE;
+ CCTK_REAL PDstandardNth2K12 = INITVALUE;
+ CCTK_REAL PDstandardNth3K12 = INITVALUE;
+ CCTK_REAL PDstandardNth1K13 = INITVALUE;
+ CCTK_REAL PDstandardNth2K13 = INITVALUE;
+ CCTK_REAL PDstandardNth3K13 = INITVALUE;
+ CCTK_REAL PDstandardNth1K22 = INITVALUE;
+ CCTK_REAL PDstandardNth3K22 = INITVALUE;
+ CCTK_REAL PDstandardNth1K23 = INITVALUE;
+ CCTK_REAL PDstandardNth2K23 = INITVALUE;
+ CCTK_REAL PDstandardNth3K23 = INITVALUE;
+ CCTK_REAL PDstandardNth1K33 = INITVALUE;
+ CCTK_REAL PDstandardNth2K33 = INITVALUE;
/* Assign local copies of grid functions */
g11L = g11[index];
@@ -211,60 +211,60 @@ void ML_ADM_constraints_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL
/* Include user supplied include files */
/* Precompute derivatives (new style) */
- PDstandard4th1g11 = PDstandard4th1(g11, i, j, k);
- PDstandard4th2g11 = PDstandard4th2(g11, i, j, k);
- PDstandard4th3g11 = PDstandard4th3(g11, i, j, k);
- PDstandard4th22g11 = PDstandard4th22(g11, i, j, k);
- PDstandard4th33g11 = PDstandard4th33(g11, i, j, k);
- PDstandard4th23g11 = PDstandard4th23(g11, i, j, k);
- PDstandard4th1g12 = PDstandard4th1(g12, i, j, k);
- PDstandard4th2g12 = PDstandard4th2(g12, i, j, k);
- PDstandard4th3g12 = PDstandard4th3(g12, i, j, k);
- PDstandard4th33g12 = PDstandard4th33(g12, i, j, k);
- PDstandard4th12g12 = PDstandard4th12(g12, i, j, k);
- PDstandard4th13g12 = PDstandard4th13(g12, i, j, k);
- PDstandard4th23g12 = PDstandard4th23(g12, i, j, k);
- PDstandard4th1g13 = PDstandard4th1(g13, i, j, k);
- PDstandard4th2g13 = PDstandard4th2(g13, i, j, k);
- PDstandard4th3g13 = PDstandard4th3(g13, i, j, k);
- PDstandard4th22g13 = PDstandard4th22(g13, i, j, k);
- PDstandard4th12g13 = PDstandard4th12(g13, i, j, k);
- PDstandard4th13g13 = PDstandard4th13(g13, i, j, k);
- PDstandard4th23g13 = PDstandard4th23(g13, i, j, k);
- PDstandard4th1g22 = PDstandard4th1(g22, i, j, k);
- PDstandard4th2g22 = PDstandard4th2(g22, i, j, k);
- PDstandard4th3g22 = PDstandard4th3(g22, i, j, k);
- PDstandard4th11g22 = PDstandard4th11(g22, i, j, k);
- PDstandard4th33g22 = PDstandard4th33(g22, i, j, k);
- PDstandard4th13g22 = PDstandard4th13(g22, i, j, k);
- PDstandard4th1g23 = PDstandard4th1(g23, i, j, k);
- PDstandard4th2g23 = PDstandard4th2(g23, i, j, k);
- PDstandard4th3g23 = PDstandard4th3(g23, i, j, k);
- PDstandard4th11g23 = PDstandard4th11(g23, i, j, k);
- PDstandard4th12g23 = PDstandard4th12(g23, i, j, k);
- PDstandard4th13g23 = PDstandard4th13(g23, i, j, k);
- PDstandard4th23g23 = PDstandard4th23(g23, i, j, k);
- PDstandard4th1g33 = PDstandard4th1(g33, i, j, k);
- PDstandard4th2g33 = PDstandard4th2(g33, i, j, k);
- PDstandard4th3g33 = PDstandard4th3(g33, i, j, k);
- PDstandard4th11g33 = PDstandard4th11(g33, i, j, k);
- PDstandard4th22g33 = PDstandard4th22(g33, i, j, k);
- PDstandard4th12g33 = PDstandard4th12(g33, i, j, k);
- PDstandard4th2K11 = PDstandard4th2(K11, i, j, k);
- PDstandard4th3K11 = PDstandard4th3(K11, i, j, k);
- PDstandard4th1K12 = PDstandard4th1(K12, i, j, k);
- PDstandard4th2K12 = PDstandard4th2(K12, i, j, k);
- PDstandard4th3K12 = PDstandard4th3(K12, i, j, k);
- PDstandard4th1K13 = PDstandard4th1(K13, i, j, k);
- PDstandard4th2K13 = PDstandard4th2(K13, i, j, k);
- PDstandard4th3K13 = PDstandard4th3(K13, i, j, k);
- PDstandard4th1K22 = PDstandard4th1(K22, i, j, k);
- PDstandard4th3K22 = PDstandard4th3(K22, i, j, k);
- PDstandard4th1K23 = PDstandard4th1(K23, i, j, k);
- PDstandard4th2K23 = PDstandard4th2(K23, i, j, k);
- PDstandard4th3K23 = PDstandard4th3(K23, i, j, k);
- PDstandard4th1K33 = PDstandard4th1(K33, i, j, k);
- PDstandard4th2K33 = PDstandard4th2(K33, i, j, k);
+ PDstandardNth1g11 = PDstandardNth1(g11, i, j, k);
+ PDstandardNth2g11 = PDstandardNth2(g11, i, j, k);
+ PDstandardNth3g11 = PDstandardNth3(g11, i, j, k);
+ PDstandardNth22g11 = PDstandardNth22(g11, i, j, k);
+ PDstandardNth33g11 = PDstandardNth33(g11, i, j, k);
+ PDstandardNth23g11 = PDstandardNth23(g11, i, j, k);
+ PDstandardNth1g12 = PDstandardNth1(g12, i, j, k);
+ PDstandardNth2g12 = PDstandardNth2(g12, i, j, k);
+ PDstandardNth3g12 = PDstandardNth3(g12, i, j, k);
+ PDstandardNth33g12 = PDstandardNth33(g12, i, j, k);
+ PDstandardNth12g12 = PDstandardNth12(g12, i, j, k);
+ PDstandardNth13g12 = PDstandardNth13(g12, i, j, k);
+ PDstandardNth23g12 = PDstandardNth23(g12, i, j, k);
+ PDstandardNth1g13 = PDstandardNth1(g13, i, j, k);
+ PDstandardNth2g13 = PDstandardNth2(g13, i, j, k);
+ PDstandardNth3g13 = PDstandardNth3(g13, i, j, k);
+ PDstandardNth22g13 = PDstandardNth22(g13, i, j, k);
+ PDstandardNth12g13 = PDstandardNth12(g13, i, j, k);
+ PDstandardNth13g13 = PDstandardNth13(g13, i, j, k);
+ PDstandardNth23g13 = PDstandardNth23(g13, i, j, k);
+ PDstandardNth1g22 = PDstandardNth1(g22, i, j, k);
+ PDstandardNth2g22 = PDstandardNth2(g22, i, j, k);
+ PDstandardNth3g22 = PDstandardNth3(g22, i, j, k);
+ PDstandardNth11g22 = PDstandardNth11(g22, i, j, k);
+ PDstandardNth33g22 = PDstandardNth33(g22, i, j, k);
+ PDstandardNth13g22 = PDstandardNth13(g22, i, j, k);
+ PDstandardNth1g23 = PDstandardNth1(g23, i, j, k);
+ PDstandardNth2g23 = PDstandardNth2(g23, i, j, k);
+ PDstandardNth3g23 = PDstandardNth3(g23, i, j, k);
+ PDstandardNth11g23 = PDstandardNth11(g23, i, j, k);
+ PDstandardNth12g23 = PDstandardNth12(g23, i, j, k);
+ PDstandardNth13g23 = PDstandardNth13(g23, i, j, k);
+ PDstandardNth23g23 = PDstandardNth23(g23, i, j, k);
+ PDstandardNth1g33 = PDstandardNth1(g33, i, j, k);
+ PDstandardNth2g33 = PDstandardNth2(g33, i, j, k);
+ PDstandardNth3g33 = PDstandardNth3(g33, i, j, k);
+ PDstandardNth11g33 = PDstandardNth11(g33, i, j, k);
+ PDstandardNth22g33 = PDstandardNth22(g33, i, j, k);
+ PDstandardNth12g33 = PDstandardNth12(g33, i, j, k);
+ PDstandardNth2K11 = PDstandardNth2(K11, i, j, k);
+ PDstandardNth3K11 = PDstandardNth3(K11, i, j, k);
+ PDstandardNth1K12 = PDstandardNth1(K12, i, j, k);
+ PDstandardNth2K12 = PDstandardNth2(K12, i, j, k);
+ PDstandardNth3K12 = PDstandardNth3(K12, i, j, k);
+ PDstandardNth1K13 = PDstandardNth1(K13, i, j, k);
+ PDstandardNth2K13 = PDstandardNth2(K13, i, j, k);
+ PDstandardNth3K13 = PDstandardNth3(K13, i, j, k);
+ PDstandardNth1K22 = PDstandardNth1(K22, i, j, k);
+ PDstandardNth3K22 = PDstandardNth3(K22, i, j, k);
+ PDstandardNth1K23 = PDstandardNth1(K23, i, j, k);
+ PDstandardNth2K23 = PDstandardNth2(K23, i, j, k);
+ PDstandardNth3K23 = PDstandardNth3(K23, i, j, k);
+ PDstandardNth1K33 = PDstandardNth1(K33, i, j, k);
+ PDstandardNth2K33 = PDstandardNth2(K33, i, j, k);
/* Precompute derivatives (old style) */
@@ -283,95 +283,95 @@ void ML_ADM_constraints_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL
gu33 = INV(detg)*(g11L*g22L - SQR(g12L));
- G111 = khalf*(gu11*PDstandard4th1g11 + 2*(gu21*PDstandard4th1g12 + gu31*PDstandard4th1g13) - gu21*PDstandard4th2g11 -
- gu31*PDstandard4th3g11);
+ G111 = khalf*(gu11*PDstandardNth1g11 + 2*(gu21*PDstandardNth1g12 + gu31*PDstandardNth1g13) - gu21*PDstandardNth2g11 -
+ gu31*PDstandardNth3g11);
- G211 = khalf*(gu21*PDstandard4th1g11 + 2*(gu22*PDstandard4th1g12 + gu32*PDstandard4th1g13) - gu22*PDstandard4th2g11 -
- gu32*PDstandard4th3g11);
+ G211 = khalf*(gu21*PDstandardNth1g11 + 2*(gu22*PDstandardNth1g12 + gu32*PDstandardNth1g13) - gu22*PDstandardNth2g11 -
+ gu32*PDstandardNth3g11);
- G311 = khalf*(gu31*PDstandard4th1g11 + 2*(gu32*PDstandard4th1g12 + gu33*PDstandard4th1g13) - gu32*PDstandard4th2g11 -
- gu33*PDstandard4th3g11);
+ G311 = khalf*(gu31*PDstandardNth1g11 + 2*(gu32*PDstandardNth1g12 + gu33*PDstandardNth1g13) - gu32*PDstandardNth2g11 -
+ gu33*PDstandardNth3g11);
- G112 = khalf*(gu21*PDstandard4th1g22 + gu11*PDstandard4th2g11 +
- gu31*(PDstandard4th1g23 + PDstandard4th2g13 - PDstandard4th3g12));
+ G112 = khalf*(gu21*PDstandardNth1g22 + gu11*PDstandardNth2g11 +
+ gu31*(PDstandardNth1g23 + PDstandardNth2g13 - PDstandardNth3g12));
- G212 = khalf*(gu22*PDstandard4th1g22 + gu21*PDstandard4th2g11 +
- gu32*(PDstandard4th1g23 + PDstandard4th2g13 - PDstandard4th3g12));
+ G212 = khalf*(gu22*PDstandardNth1g22 + gu21*PDstandardNth2g11 +
+ gu32*(PDstandardNth1g23 + PDstandardNth2g13 - PDstandardNth3g12));
- G312 = khalf*(gu32*PDstandard4th1g22 + gu31*PDstandard4th2g11 +
- gu33*(PDstandard4th1g23 + PDstandard4th2g13 - PDstandard4th3g12));
+ G312 = khalf*(gu32*PDstandardNth1g22 + gu31*PDstandardNth2g11 +
+ gu33*(PDstandardNth1g23 + PDstandardNth2g13 - PDstandardNth3g12));
- G113 = khalf*(gu31*PDstandard4th1g33 + gu11*PDstandard4th3g11 +
- gu21*(PDstandard4th1g23 - PDstandard4th2g13 + PDstandard4th3g12));
+ G113 = khalf*(gu31*PDstandardNth1g33 + gu11*PDstandardNth3g11 +
+ gu21*(PDstandardNth1g23 - PDstandardNth2g13 + PDstandardNth3g12));
- G213 = khalf*(gu32*PDstandard4th1g33 + gu21*PDstandard4th3g11 +
- gu22*(PDstandard4th1g23 - PDstandard4th2g13 + PDstandard4th3g12));
+ G213 = khalf*(gu32*PDstandardNth1g33 + gu21*PDstandardNth3g11 +
+ gu22*(PDstandardNth1g23 - PDstandardNth2g13 + PDstandardNth3g12));
- G313 = khalf*(gu33*PDstandard4th1g33 + gu31*PDstandard4th3g11 +
- gu32*(PDstandard4th1g23 - PDstandard4th2g13 + PDstandard4th3g12));
+ G313 = khalf*(gu33*PDstandardNth1g33 + gu31*PDstandardNth3g11 +
+ gu32*(PDstandardNth1g23 - PDstandardNth2g13 + PDstandardNth3g12));
- G122 = khalf*(gu11*(-PDstandard4th1g22 + 2*PDstandard4th2g12) + gu21*PDstandard4th2g22 +
- gu31*(2*PDstandard4th2g23 - PDstandard4th3g22));
+ G122 = khalf*(gu11*(-PDstandardNth1g22 + 2*PDstandardNth2g12) + gu21*PDstandardNth2g22 +
+ gu31*(2*PDstandardNth2g23 - PDstandardNth3g22));
- G222 = khalf*(gu21*(-PDstandard4th1g22 + 2*PDstandard4th2g12) + gu22*PDstandard4th2g22 +
- gu32*(2*PDstandard4th2g23 - PDstandard4th3g22));
+ G222 = khalf*(gu21*(-PDstandardNth1g22 + 2*PDstandardNth2g12) + gu22*PDstandardNth2g22 +
+ gu32*(2*PDstandardNth2g23 - PDstandardNth3g22));
- G322 = khalf*(gu31*(-PDstandard4th1g22 + 2*PDstandard4th2g12) + gu32*PDstandard4th2g22 +
- gu33*(2*PDstandard4th2g23 - PDstandard4th3g22));
+ G322 = khalf*(gu31*(-PDstandardNth1g22 + 2*PDstandardNth2g12) + gu32*PDstandardNth2g22 +
+ gu33*(2*PDstandardNth2g23 - PDstandardNth3g22));
- G123 = khalf*(gu31*PDstandard4th2g33 + gu11*(-PDstandard4th1g23 + PDstandard4th2g13 + PDstandard4th3g12) +
- gu21*PDstandard4th3g22);
+ G123 = khalf*(gu31*PDstandardNth2g33 + gu11*(-PDstandardNth1g23 + PDstandardNth2g13 + PDstandardNth3g12) +
+ gu21*PDstandardNth3g22);
- G223 = khalf*(gu32*PDstandard4th2g33 + gu21*(-PDstandard4th1g23 + PDstandard4th2g13 + PDstandard4th3g12) +
- gu22*PDstandard4th3g22);
+ G223 = khalf*(gu32*PDstandardNth2g33 + gu21*(-PDstandardNth1g23 + PDstandardNth2g13 + PDstandardNth3g12) +
+ gu22*PDstandardNth3g22);
- G323 = khalf*(gu33*PDstandard4th2g33 + gu31*(-PDstandard4th1g23 + PDstandard4th2g13 + PDstandard4th3g12) +
- gu32*PDstandard4th3g22);
+ G323 = khalf*(gu33*PDstandardNth2g33 + gu31*(-PDstandardNth1g23 + PDstandardNth2g13 + PDstandardNth3g12) +
+ gu32*PDstandardNth3g22);
- G133 = khalf*(-(gu11*PDstandard4th1g33) - gu21*PDstandard4th2g33 + 2*gu11*PDstandard4th3g13 +
- 2*gu21*PDstandard4th3g23 + gu31*PDstandard4th3g33);
+ G133 = khalf*(-(gu11*PDstandardNth1g33) - gu21*PDstandardNth2g33 + 2*gu11*PDstandardNth3g13 +
+ 2*gu21*PDstandardNth3g23 + gu31*PDstandardNth3g33);
- G233 = khalf*(-(gu21*PDstandard4th1g33) - gu22*PDstandard4th2g33 + 2*gu21*PDstandard4th3g13 +
- 2*gu22*PDstandard4th3g23 + gu32*PDstandard4th3g33);
+ G233 = khalf*(-(gu21*PDstandardNth1g33) - gu22*PDstandardNth2g33 + 2*gu21*PDstandardNth3g13 +
+ 2*gu22*PDstandardNth3g23 + gu32*PDstandardNth3g33);
- G333 = khalf*(-(gu31*PDstandard4th1g33) - gu32*PDstandard4th2g33 + 2*gu31*PDstandard4th3g13 +
- 2*gu32*PDstandard4th3g23 + gu33*PDstandard4th3g33);
+ G333 = khalf*(-(gu31*PDstandardNth1g33) - gu32*PDstandardNth2g33 + 2*gu31*PDstandardNth3g13 +
+ 2*gu32*PDstandardNth3g23 + gu33*PDstandardNth3g33);
- R11 = khalf*(-(gu22*PDstandard4th11g22) - 2*(G111*G122 + G111*G133 + G211*G222 + G211*G233 + G311*G322 + G311*G333 +
- gu32*PDstandard4th11g23) - gu33*PDstandard4th11g33 + 2*gu22*PDstandard4th12g12 + 2*gu32*PDstandard4th12g13 +
- 2*gu32*PDstandard4th13g12 + 2*gu33*PDstandard4th13g13 - gu22*PDstandard4th22g11 - 2*gu32*PDstandard4th23g11 -
- gu33*PDstandard4th33g11 + 2*SQR(G112) + 2*SQR(G113) + 2*SQR(G212) + 2*SQR(G213) + 2*SQR(G312) + 2*SQR(G313));
+ R11 = khalf*(-(gu22*PDstandardNth11g22) - 2*(G111*G122 + G111*G133 + G211*G222 + G211*G233 + G311*G322 + G311*G333 +
+ gu32*PDstandardNth11g23) - gu33*PDstandardNth11g33 + 2*gu22*PDstandardNth12g12 + 2*gu32*PDstandardNth12g13 +
+ 2*gu32*PDstandardNth13g12 + 2*gu33*PDstandardNth13g13 - gu22*PDstandardNth22g11 - 2*gu32*PDstandardNth23g11 -
+ gu33*PDstandardNth33g11 + 2*SQR(G112) + 2*SQR(G113) + 2*SQR(G212) + 2*SQR(G213) + 2*SQR(G312) + 2*SQR(G313));
- R12 = khalf*(2*(G113*G123 + G213*G223 + G313*G323) - 2*(G112*G133 + G212*G233 + G312*G333 + gu21*PDstandard4th12g12) -
- gu32*PDstandard4th12g23 - gu33*PDstandard4th12g33 +
- gu31*(PDstandard4th11g23 - PDstandard4th12g13 - PDstandard4th13g12) + gu32*PDstandard4th13g22 +
- gu33*PDstandard4th13g23 + gu21*(PDstandard4th11g22 + PDstandard4th22g11) + gu32*PDstandard4th22g13 +
- gu31*PDstandard4th23g11 - gu32*PDstandard4th23g12 + gu33*PDstandard4th23g13 - gu33*PDstandard4th33g12);
+ R12 = khalf*(2*(G113*G123 + G213*G223 + G313*G323) - 2*(G112*G133 + G212*G233 + G312*G333 + gu21*PDstandardNth12g12) -
+ gu32*PDstandardNth12g23 - gu33*PDstandardNth12g33 +
+ gu31*(PDstandardNth11g23 - PDstandardNth12g13 - PDstandardNth13g12) + gu32*PDstandardNth13g22 +
+ gu33*PDstandardNth13g23 + gu21*(PDstandardNth11g22 + PDstandardNth22g11) + gu32*PDstandardNth22g13 +
+ gu31*PDstandardNth23g11 - gu32*PDstandardNth23g12 + gu33*PDstandardNth23g13 - gu33*PDstandardNth33g12);
- R13 = khalf*(2*(G112*G123 + G212*G223 + G312*G323) - 2*(G113*G122 + G213*G222 + G313*G322 + gu31*PDstandard4th13g13) +
- gu21*(PDstandard4th11g23 - PDstandard4th12g13 - PDstandard4th13g12 + PDstandard4th23g11) +
- gu22*(PDstandard4th12g23 - PDstandard4th13g22 - PDstandard4th22g13 + PDstandard4th23g12) +
- gu31*(PDstandard4th11g33 + PDstandard4th33g11) +
- gu32*(PDstandard4th12g33 - PDstandard4th13g23 - PDstandard4th23g13 + PDstandard4th33g12));
+ R13 = khalf*(2*(G112*G123 + G212*G223 + G312*G323) - 2*(G113*G122 + G213*G222 + G313*G322 + gu31*PDstandardNth13g13) +
+ gu21*(PDstandardNth11g23 - PDstandardNth12g13 - PDstandardNth13g12 + PDstandardNth23g11) +
+ gu22*(PDstandardNth12g23 - PDstandardNth13g22 - PDstandardNth22g13 + PDstandardNth23g12) +
+ gu31*(PDstandardNth11g33 + PDstandardNth33g11) +
+ gu32*(PDstandardNth12g33 - PDstandardNth13g23 - PDstandardNth23g13 + PDstandardNth33g12));
- R22 = khalf*(-2*(G122*(G111 + G133) + G222*(G211 + G233) + G322*(G311 + G333) + gu31*PDstandard4th13g22) +
- gu11*(-PDstandard4th11g22 + 2*PDstandard4th12g12 - PDstandard4th22g11) +
- gu31*(-2*PDstandard4th22g13 + 2*(PDstandard4th12g23 + PDstandard4th23g12)) +
- gu33*(-PDstandard4th22g33 + 2*PDstandard4th23g23 - PDstandard4th33g22) +
+ R22 = khalf*(-2*(G122*(G111 + G133) + G222*(G211 + G233) + G322*(G311 + G333) + gu31*PDstandardNth13g22) +
+ gu11*(-PDstandardNth11g22 + 2*PDstandardNth12g12 - PDstandardNth22g11) +
+ gu31*(-2*PDstandardNth22g13 + 2*(PDstandardNth12g23 + PDstandardNth23g12)) +
+ gu33*(-PDstandardNth22g33 + 2*PDstandardNth23g23 - PDstandardNth33g22) +
2*(SQR(G112) + SQR(G123) + SQR(G212) + SQR(G223) + SQR(G312) + SQR(G323)));
R23 = khalf*(2*(G112*G113 + G212*G213 + G312*G313) +
- gu11*(-PDstandard4th11g23 + PDstandard4th12g13 + PDstandard4th13g12 - PDstandard4th23g11) +
- gu21*(-PDstandard4th12g23 + PDstandard4th13g22 + PDstandard4th22g13 - PDstandard4th23g12) -
- 2*(G111*G123 + G211*G223 + G311*G323 + gu32*PDstandard4th23g23) +
- gu31*(PDstandard4th12g33 - PDstandard4th13g23 - PDstandard4th23g13 + PDstandard4th33g12) +
- gu32*(PDstandard4th22g33 + PDstandard4th33g22));
+ gu11*(-PDstandardNth11g23 + PDstandardNth12g13 + PDstandardNth13g12 - PDstandardNth23g11) +
+ gu21*(-PDstandardNth12g23 + PDstandardNth13g22 + PDstandardNth22g13 - PDstandardNth23g12) -
+ 2*(G111*G123 + G211*G223 + G311*G323 + gu32*PDstandardNth23g23) +
+ gu31*(PDstandardNth12g33 - PDstandardNth13g23 - PDstandardNth23g13 + PDstandardNth33g12) +
+ gu32*(PDstandardNth22g33 + PDstandardNth33g22));
- R33 = khalf*(gu11*(-PDstandard4th11g33 + 2*PDstandard4th13g13 - PDstandard4th33g11) -
+ R33 = khalf*(gu11*(-PDstandardNth11g33 + 2*PDstandardNth13g13 - PDstandardNth33g11) -
2*((G111 + G122)*G133 + (G211 + G222)*G233 + (G311 + G322)*G333 +
- gu21*(PDstandard4th12g33 + PDstandard4th33g12)) +
- gu22*(-PDstandard4th22g33 + 2*PDstandard4th23g23 - PDstandard4th33g22) +
- 2*(gu21*(PDstandard4th13g23 + PDstandard4th23g13) + SQR(G113) + SQR(G123) + SQR(G213) + SQR(G223) + SQR(G313) +
+ gu21*(PDstandardNth12g33 + PDstandardNth33g12)) +
+ gu22*(-PDstandardNth22g33 + 2*PDstandardNth23g23 - PDstandardNth33g22) +
+ 2*(gu21*(PDstandardNth13g23 + PDstandardNth23g13) + SQR(G113) + SQR(G123) + SQR(G213) + SQR(G223) + SQR(G313) +
SQR(G323)));
trR = gu11*R11 + gu22*R22 + 2*(gu21*R12 + gu31*R13 + gu32*R23) + gu33*R33;
@@ -398,31 +398,31 @@ void ML_ADM_constraints_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL
HL = -2*(Km12*Km21 + Km13*Km31 + Km23*Km32) + trR - SQR(Km11) - SQR(Km22) - SQR(Km33) + SQR(trK);
- M1L = gu21*(-(G112*K11L) + G111*K12L - G212*K12L - G312*K13L + G211*K22L + G311*K23L - PDstandard4th1K12 +
- PDstandard4th2K11) + gu22*(-(G122*K11L) + G112*K12L - G222*K12L - G322*K13L + G212*K22L + G312*K23L -
- PDstandard4th1K22 + PDstandard4th2K12) + gu31*
- (-(G113*K11L) - G213*K12L + G111*K13L - G313*K13L + G211*K23L + G311*K33L - PDstandard4th1K13 + PDstandard4th3K11)\
+ M1L = gu21*(-(G112*K11L) + G111*K12L - G212*K12L - G312*K13L + G211*K22L + G311*K23L - PDstandardNth1K12 +
+ PDstandardNth2K11) + gu22*(-(G122*K11L) + G112*K12L - G222*K12L - G322*K13L + G212*K22L + G312*K23L -
+ PDstandardNth1K22 + PDstandardNth2K12) + gu31*
+ (-(G113*K11L) - G213*K12L + G111*K13L - G313*K13L + G211*K23L + G311*K33L - PDstandardNth1K13 + PDstandardNth3K11)\
+ gu32*(G113*K12L + G112*K13L + G213*K22L + (G212 + G313)*K23L + G312*K33L -
- 2*(G123*K11L + G223*K12L + G323*K13L + PDstandard4th1K23) + PDstandard4th2K13 + PDstandard4th3K12) +
- gu33*(-(G133*K11L) - G233*K12L + G113*K13L - G333*K13L + G213*K23L + G313*K33L - PDstandard4th1K33 +
- PDstandard4th3K13);
-
- M2L = gu11*(G112*K11L + (-G111 + G212)*K12L + G312*K13L - G211*K22L - G311*K23L + PDstandard4th1K12 -
- PDstandard4th2K11) + gu21*(G122*K11L + (-G112 + G222)*K12L + G322*K13L - G212*K22L - G312*K23L +
- PDstandard4th1K22 - PDstandard4th2K12) + gu31*
- (G123*K11L + (-2*G113 + G223)*K12L + (G112 + G323)*K13L + G212*K23L + G312*K33L + PDstandard4th1K23 -
- 2*(G213*K22L + G313*K23L + PDstandard4th2K13) + PDstandard4th3K12) +
- gu32*(-(G123*K12L) + G122*K13L - G223*K22L + G222*K23L - G323*K23L + G322*K33L - PDstandard4th2K23 +
- PDstandard4th3K22) + gu33*(-(G133*K12L) + G123*K13L - G233*K22L + G223*K23L - G333*K23L + G323*K33L -
- PDstandard4th2K33 + PDstandard4th3K23);
-
- M3L = gu11*(G113*K11L + G213*K12L + (-G111 + G313)*K13L - G211*K23L - G311*K33L + PDstandard4th1K13 -
- PDstandard4th3K11) + gu21*(G123*K11L + (G113 + G223)*K12L + (-2*G112 + G323)*K13L + G213*K22L +
- (-2*G212 + G313)*K23L + PDstandard4th1K23 + PDstandard4th2K13 - 2*(G312*K33L + PDstandard4th3K12)) +
- gu31*(G133*K11L + G233*K12L + (-G113 + G333)*K13L - G213*K23L - G313*K33L + PDstandard4th1K33 -
- PDstandard4th3K13) + gu22*(G123*K12L - G122*K13L + G223*K22L - G222*K23L + G323*K23L - G322*K33L +
- PDstandard4th2K23 - PDstandard4th3K22) + gu32*
- (G133*K12L - G123*K13L + G233*K22L - G223*K23L + G333*K23L - G323*K33L + PDstandard4th2K33 - PDstandard4th3K23);
+ 2*(G123*K11L + G223*K12L + G323*K13L + PDstandardNth1K23) + PDstandardNth2K13 + PDstandardNth3K12) +
+ gu33*(-(G133*K11L) - G233*K12L + G113*K13L - G333*K13L + G213*K23L + G313*K33L - PDstandardNth1K33 +
+ PDstandardNth3K13);
+
+ M2L = gu11*(G112*K11L + (-G111 + G212)*K12L + G312*K13L - G211*K22L - G311*K23L + PDstandardNth1K12 -
+ PDstandardNth2K11) + gu21*(G122*K11L + (-G112 + G222)*K12L + G322*K13L - G212*K22L - G312*K23L +
+ PDstandardNth1K22 - PDstandardNth2K12) + gu31*
+ (G123*K11L + (-2*G113 + G223)*K12L + (G112 + G323)*K13L + G212*K23L + G312*K33L + PDstandardNth1K23 -
+ 2*(G213*K22L + G313*K23L + PDstandardNth2K13) + PDstandardNth3K12) +
+ gu32*(-(G123*K12L) + G122*K13L - G223*K22L + G222*K23L - G323*K23L + G322*K33L - PDstandardNth2K23 +
+ PDstandardNth3K22) + gu33*(-(G133*K12L) + G123*K13L - G233*K22L + G223*K23L - G333*K23L + G323*K33L -
+ PDstandardNth2K33 + PDstandardNth3K23);
+
+ M3L = gu11*(G113*K11L + G213*K12L + (-G111 + G313)*K13L - G211*K23L - G311*K33L + PDstandardNth1K13 -
+ PDstandardNth3K11) + gu21*(G123*K11L + (G113 + G223)*K12L + (-2*G112 + G323)*K13L + G213*K22L +
+ (-2*G212 + G313)*K23L + PDstandardNth1K23 + PDstandardNth2K13 - 2*(G312*K33L + PDstandardNth3K12)) +
+ gu31*(G133*K11L + G233*K12L + (-G113 + G333)*K13L - G213*K23L - G313*K33L + PDstandardNth1K33 -
+ PDstandardNth3K13) + gu22*(G123*K12L - G122*K13L + G223*K22L - G222*K23L + G323*K23L - G322*K33L +
+ PDstandardNth2K23 - PDstandardNth3K22) + gu32*
+ (G133*K12L - G123*K13L + G233*K22L - G223*K23L + G333*K23L - G323*K33L + PDstandardNth2K33 - PDstandardNth3K23);
/* Copy local copies back to grid functions */
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-29 10:00:48 +0000