Update generated code

Update the generated code after the recent changes to Kranc.

4 files changed, 46 insertions, 44 deletions

diff --git a/ML_WaveToy/param.ccl b/ML_WaveToy/param.ccl
index 17883ec..5964f1e 100644
--- a/ML_WaveToy/param.ccl
+++ b/ML_WaveToy/param.ccl
@@ -21,7 +21,7 @@ USES CCTK_INT MoL_Num_Constrained_Vars
 restricted:
 CCTK_REAL verbose "verbose"
 {
-  "*:*" :: "no restrictions"
+  "*:*" :: ""
 } 0
 
 restricted:
@@ -39,25 +39,25 @@ CCTK_INT ML_WaveToy_MaxNumConstrainedVars "Number of constrained variables used
 restricted:
 CCTK_INT WT_Gaussian_calc_every "WT_Gaussian_calc_every"
 {
-  *:* :: "no restrictions"
+  *:* :: ""
 } 1
 
 restricted:
 CCTK_INT WT_RHS_calc_every "WT_RHS_calc_every"
 {
-  *:* :: "no restrictions"
+  *:* :: ""
 } 1
 
 restricted:
 CCTK_INT WT_Gaussian_calc_offset "WT_Gaussian_calc_offset"
 {
-  *:* :: "no restrictions"
+  *:* :: ""
 } 0
 
 restricted:
 CCTK_INT WT_RHS_calc_offset "WT_RHS_calc_offset"
 {
-  *:* :: "no restrictions"
+  *:* :: ""
 } 0
 
 private:
diff --git a/ML_WaveToy/src/Differencing.h b/ML_WaveToy/src/Differencing.h
index 626c54f..fd9befd 100644
--- a/ML_WaveToy/src/Differencing.h
+++ b/ML_WaveToy/src/Differencing.h
@@ -1,39 +1,39 @@
-#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)]))
-#define PDstandardNth111(u,i,j,k) (p1o2dx3*(2*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 2*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 PDstandardNth222(u,i,j,k) (p1o2dy3*(2*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 2*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 PDstandardNth333(u,i,j,k) (p1o2dz3*(2*u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 2*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 PDstandardNth112(u,i,j,k) (p1o144dx2dy*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*(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)] + 30*(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 PDstandardNth113(u,i,j,k) (p1o144dx2dz*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*(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)] + 30*(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 PDstandardNth221(u,i,j,k) (p1o144dxdy2*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)]) + 128*(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)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)]) - 16*(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)] - 30*u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
-#define PDstandardNth223(u,i,j,k) (p1o144dy2dz*(-128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*(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)] + 30*(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 PDstandardNth331(u,i,j,k) (p1o144dxdz2*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)]) + 128*(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)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)]) - 16*(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)] - 30*u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
-#define PDstandardNth332(u,i,j,k) (p1o144dydz2*(-128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)]) + 128*(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)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)]) - 16*(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)] - 30*u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))
-#define PDstandardNth121(u,i,j,k) (p1o144dx2dy*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*(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)] + 30*(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 PDstandardNth131(u,i,j,k) (p1o144dx2dz*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*(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)] + 30*(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 PDstandardNth212(u,i,j,k) (p1o144dxdy2*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)]) + 128*(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)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)]) - 16*(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)] - 30*u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
-#define PDstandardNth232(u,i,j,k) (p1o144dy2dz*(-128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*(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)] + 30*(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 PDstandardNth313(u,i,j,k) (p1o144dxdz2*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)]) + 128*(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)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)]) - 16*(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)] - 30*u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
-#define PDstandardNth323(u,i,j,k) (p1o144dydz2*(-128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)]) + 128*(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)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)]) - 16*(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)] - 30*u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))
-#define PDstandardNth211(u,i,j,k) (p1o144dx2dy*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*(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)] + 30*(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 PDstandardNth311(u,i,j,k) (p1o144dx2dz*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*(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)] + 30*(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 PDstandardNth122(u,i,j,k) (p1o144dxdy2*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)]) + 128*(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)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)]) - 16*(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)] - 30*u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
-#define PDstandardNth322(u,i,j,k) (p1o144dy2dz*(-128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*(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)] + 30*(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 PDstandardNth133(u,i,j,k) (p1o144dxdz2*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)]) + 128*(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)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)]) - 16*(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)] - 30*u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
-#define PDstandardNth233(u,i,j,k) (p1o144dydz2*(-128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)]) + 128*(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)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)]) - 16*(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)] - 30*u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))
-#define PDstandardNth123(u,i,j,k) (p1o1728dxdydz*(-512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)]))
-#define PDstandardNth132(u,i,j,k) (p1o1728dxdydz*(-512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)]))
-#define PDstandardNth213(u,i,j,k) (p1o1728dxdydz*(-512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)]))
-#define PDstandardNth231(u,i,j,k) (p1o1728dxdydz*(-512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)]))
-#define PDstandardNth312(u,i,j,k) (p1o1728dxdydz*(-512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)]))
-#define PDstandardNth321(u,i,j,k) (p1o1728dxdydz*(-512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + 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)]))
+#define PDstandardNth111(u,i,j,k) (p1o2dx3*(2*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 2*(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 PDstandardNth222(u,i,j,k) (p1o2dy3*(2*(u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 2*(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 PDstandardNth333(u,i,j,k) (p1o2dz3*(2*(u)[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 2*(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 PDstandardNth112(u,i,j,k) (p1o144dx2dy*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*((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)] + 30*(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 PDstandardNth113(u,i,j,k) (p1o144dx2dz*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*((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)] + 30*(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 PDstandardNth221(u,i,j,k) (p1o144dxdy2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)]) + 128*((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)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)]) - 16*((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)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
+#define PDstandardNth223(u,i,j,k) (p1o144dy2dz*(-128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*((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)] + 30*(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 PDstandardNth331(u,i,j,k) (p1o144dxdz2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)]) + 128*((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)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)]) - 16*((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)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
+#define PDstandardNth332(u,i,j,k) (p1o144dydz2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)]) + 128*((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)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)]) - 16*((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)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))
+#define PDstandardNth121(u,i,j,k) (p1o144dx2dy*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*((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)] + 30*(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 PDstandardNth131(u,i,j,k) (p1o144dx2dz*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*((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)] + 30*(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 PDstandardNth212(u,i,j,k) (p1o144dxdy2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)]) + 128*((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)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)]) - 16*((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)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
+#define PDstandardNth232(u,i,j,k) (p1o144dy2dz*(-128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*((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)] + 30*(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 PDstandardNth313(u,i,j,k) (p1o144dxdz2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)]) + 128*((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)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)]) - 16*((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)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
+#define PDstandardNth323(u,i,j,k) (p1o144dydz2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)]) + 128*((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)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)]) - 16*((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)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))
+#define PDstandardNth211(u,i,j,k) (p1o144dx2dy*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*((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)] + 30*(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 PDstandardNth311(u,i,j,k) (p1o144dx2dz*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*((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)] + 30*(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 PDstandardNth122(u,i,j,k) (p1o144dxdy2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)]) + 128*((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)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)]) - 16*((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)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
+#define PDstandardNth322(u,i,j,k) (p1o144dy2dz*(-128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*((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)] + 30*(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 PDstandardNth133(u,i,j,k) (p1o144dxdz2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)]) + 128*((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)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)]) - 16*((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)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)]))
+#define PDstandardNth233(u,i,j,k) (p1o144dydz2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)]) + 128*((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)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)]) - 16*((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)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))
+#define PDstandardNth123(u,i,j,k) (p1o1728dxdydz*(-512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)]))
+#define PDstandardNth132(u,i,j,k) (p1o1728dxdydz*(-512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)]))
+#define PDstandardNth213(u,i,j,k) (p1o1728dxdydz*(-512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)]))
+#define PDstandardNth231(u,i,j,k) (p1o1728dxdydz*(-512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)]))
+#define PDstandardNth312(u,i,j,k) (p1o1728dxdydz*(-512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)]))
+#define PDstandardNth321(u,i,j,k) (p1o1728dxdydz*(-512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)]))
diff --git a/ML_WaveToy/src/WT_Gaussian.c b/ML_WaveToy/src/WT_Gaussian.c
index bbba1a1..5caf9ca 100644
--- a/ML_WaveToy/src/WT_Gaussian.c
+++ b/ML_WaveToy/src/WT_Gaussian.c
@@ -104,6 +104,7 @@ void WT_Gaussian_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal
   pm1o12dz2 = -pow(dz,-2)/12.;
   
   /* Loop over the grid points */
+  _Pragma ("omp parallel")
   LC_LOOP3 (WT_Gaussian,
             i,j,k, min[0],min[1],min[2], max[0],max[1],max[2],
             cctk_lsh[0],cctk_lsh[1],cctk_lsh[2])
diff --git a/ML_WaveToy/src/WT_RHS.c b/ML_WaveToy/src/WT_RHS.c
index f06f42a..41d3939 100644
--- a/ML_WaveToy/src/WT_RHS.c
+++ b/ML_WaveToy/src/WT_RHS.c
@@ -104,6 +104,7 @@ void WT_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[3],
   pm1o12dz2 = -pow(dz,-2)/12.;
   
   /* Loop over the grid points */
+  _Pragma ("omp parallel")
   LC_LOOP3 (WT_RHS,
             i,j,k, min[0],min[1],min[2], max[0],max[1],max[2],
             cctk_lsh[0],cctk_lsh[1],cctk_lsh[2])