Regenerate WaveToy thorns

Correct thorn name: Use ML_FOWaveToy instead of ML_FOWavetoy (upper case T).

Regenerate thorns after correcting error in Kranc.

5 files changed, 119 insertions, 52 deletions

diff --git a/ML_WaveToy/interface.ccl b/ML_WaveToy/interface.ccl
index be30c79..515e050 100644
--- a/ML_WaveToy/interface.ccl
+++ b/ML_WaveToy/interface.ccl
@@ -43,13 +43,13 @@ CCTK_REAL WT_urhs type=GF timelevels=1 tags='tensortypealias="Scalar" tensorweig
 } "WT_urhs"
 
 public:
-CCTK_REAL WT_rho type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
+CCTK_REAL WT_rho type=GF timelevels=2 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
 {
   rho
 } "WT_rho"
 
 public:
-CCTK_REAL WT_u type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
+CCTK_REAL WT_u type=GF timelevels=2 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
 {
   u
 } "WT_u"
diff --git a/ML_WaveToy/schedule.ccl b/ML_WaveToy/schedule.ccl
index f02078a..5308b4c 100644
--- a/ML_WaveToy/schedule.ccl
+++ b/ML_WaveToy/schedule.ccl
@@ -8,9 +8,9 @@ STORAGE: WT_rhorhs[1]
 
 STORAGE: WT_urhs[1]
 
-STORAGE: WT_rho[3]
+STORAGE: WT_rho[2]
 
-STORAGE: WT_u[3]
+STORAGE: WT_u[2]
 
 schedule ML_WaveToy_Startup at STARTUP
 {
diff --git a/ML_WaveToy/src/Differencing.h b/ML_WaveToy/src/Differencing.h
index 876855b..6b60a8f 100644
--- a/ML_WaveToy/src/Differencing.h
+++ b/ML_WaveToy/src/Differencing.h
@@ -1,12 +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 PDstandardNth1(u,i,j,k) (qp1o12dx*(-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) (qp1o12dy*(-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) (qp1o12dz*(-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) (qpm1o12dx2*(-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) (qpm1o12dy2*(-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) (qpm1o12dz2*(-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) (qp1o144dxdy*(-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) (qp1o144dxdz*(-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) (qp1o144dxdy*(-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) (qp1o144dydz*(-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) (qp1o144dxdz*(-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) (qp1o144dydz*(-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) (qp1o2dx3*(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) (qp1o2dy3*(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) (qp1o2dz3*(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) (qp1o144dx2dy*(-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) (qp1o144dx2dz*(-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) (qp1o144dxdy2*(-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) (qp1o144dy2dz*(-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) (qp1o144dxdz2*(-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) (qp1o144dydz2*(-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) (qp1o144dx2dy*(-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) (qp1o144dx2dz*(-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) (qp1o144dxdy2*(-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) (qp1o144dy2dz*(-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) (qp1o144dxdz2*(-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) (qp1o144dydz2*(-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) (qp1o144dx2dy*(-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) (qp1o144dx2dz*(-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) (qp1o144dxdy2*(-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) (qp1o144dy2dz*(-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) (qp1o144dxdz2*(-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) (qp1o144dydz2*(-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) (qp1o1728dxdydz*(-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) (qp1o1728dxdydz*(-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) (qp1o1728dxdydz*(-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) (qp1o1728dxdydz*(-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) (qp1o1728dxdydz*(-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) (qp1o1728dxdydz*(-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 2b1e3b7..ffac818 100644
--- a/ML_WaveToy/src/WT_Gaussian.c
+++ b/ML_WaveToy/src/WT_Gaussian.c
@@ -34,15 +34,25 @@ void WT_Gaussian_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal
   
   
   /* Declare predefined quantities */
-  CCTK_REAL p1o12dx = INITVALUE;
-  CCTK_REAL p1o12dy = INITVALUE;
-  CCTK_REAL p1o12dz = INITVALUE;
-  CCTK_REAL p1o144dxdy = INITVALUE;
-  CCTK_REAL p1o144dxdz = INITVALUE;
-  CCTK_REAL p1o144dydz = INITVALUE;
-  CCTK_REAL pm1o12dx2 = INITVALUE;
-  CCTK_REAL pm1o12dy2 = INITVALUE;
-  CCTK_REAL pm1o12dz2 = INITVALUE;
+  CCTK_REAL qp1o12dx = INITVALUE;
+  CCTK_REAL qp1o12dy = INITVALUE;
+  CCTK_REAL qp1o12dz = INITVALUE;
+  CCTK_REAL qp1o144dx2dy = INITVALUE;
+  CCTK_REAL qp1o144dx2dz = INITVALUE;
+  CCTK_REAL qp1o144dxdy = INITVALUE;
+  CCTK_REAL qp1o144dxdy2 = INITVALUE;
+  CCTK_REAL qp1o144dxdz = INITVALUE;
+  CCTK_REAL qp1o144dxdz2 = INITVALUE;
+  CCTK_REAL qp1o144dy2dz = INITVALUE;
+  CCTK_REAL qp1o144dydz = INITVALUE;
+  CCTK_REAL qp1o144dydz2 = INITVALUE;
+  CCTK_REAL qp1o1728dxdydz = INITVALUE;
+  CCTK_REAL qp1o2dx3 = INITVALUE;
+  CCTK_REAL qp1o2dy3 = INITVALUE;
+  CCTK_REAL qp1o2dz3 = INITVALUE;
+  CCTK_REAL qpm1o12dx2 = INITVALUE;
+  CCTK_REAL qpm1o12dy2 = INITVALUE;
+  CCTK_REAL qpm1o12dz2 = INITVALUE;
   
   if (verbose > 1)
   {
@@ -73,15 +83,25 @@ void WT_Gaussian_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal
   hdzi = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  p1o12dx = INV(dx)/12.;
-  p1o12dy = INV(dy)/12.;
-  p1o12dz = INV(dz)/12.;
-  p1o144dxdy = (INV(dx)*INV(dy))/144.;
-  p1o144dxdz = (INV(dx)*INV(dz))/144.;
-  p1o144dydz = (INV(dy)*INV(dz))/144.;
-  pm1o12dx2 = -pow(dx,-2)/12.;
-  pm1o12dy2 = -pow(dy,-2)/12.;
-  pm1o12dz2 = -pow(dz,-2)/12.;
+  qp1o12dx = INV(dx)/12.;
+  qp1o12dy = INV(dy)/12.;
+  qp1o12dz = INV(dz)/12.;
+  qp1o144dx2dy = (INV(dy)*pow(dx,-2))/144.;
+  qp1o144dx2dz = (INV(dz)*pow(dx,-2))/144.;
+  qp1o144dxdy = (INV(dx)*INV(dy))/144.;
+  qp1o144dxdy2 = (INV(dx)*pow(dy,-2))/144.;
+  qp1o144dxdz = (INV(dx)*INV(dz))/144.;
+  qp1o144dxdz2 = (INV(dx)*pow(dz,-2))/144.;
+  qp1o144dy2dz = (INV(dz)*pow(dy,-2))/144.;
+  qp1o144dydz = (INV(dy)*INV(dz))/144.;
+  qp1o144dydz2 = (INV(dy)*pow(dz,-2))/144.;
+  qp1o1728dxdydz = (INV(dx)*INV(dy)*INV(dz))/1728.;
+  qp1o2dx3 = khalf*pow(dx,-3);
+  qp1o2dy3 = khalf*pow(dy,-3);
+  qp1o2dz3 = khalf*pow(dz,-3);
+  qpm1o12dx2 = -pow(dx,-2)/12.;
+  qpm1o12dy2 = -pow(dy,-2)/12.;
+  qpm1o12dz2 = -pow(dz,-2)/12.;
   
   /* Loop over the grid points */
   LC_LOOP3 (WT_Gaussian,
diff --git a/ML_WaveToy/src/WT_RHS.c b/ML_WaveToy/src/WT_RHS.c
index f9ff504..81d2144 100644
--- a/ML_WaveToy/src/WT_RHS.c
+++ b/ML_WaveToy/src/WT_RHS.c
@@ -34,15 +34,25 @@ void WT_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[3],
   
   
   /* Declare predefined quantities */
-  CCTK_REAL p1o12dx = INITVALUE;
-  CCTK_REAL p1o12dy = INITVALUE;
-  CCTK_REAL p1o12dz = INITVALUE;
-  CCTK_REAL p1o144dxdy = INITVALUE;
-  CCTK_REAL p1o144dxdz = INITVALUE;
-  CCTK_REAL p1o144dydz = INITVALUE;
-  CCTK_REAL pm1o12dx2 = INITVALUE;
-  CCTK_REAL pm1o12dy2 = INITVALUE;
-  CCTK_REAL pm1o12dz2 = INITVALUE;
+  CCTK_REAL qp1o12dx = INITVALUE;
+  CCTK_REAL qp1o12dy = INITVALUE;
+  CCTK_REAL qp1o12dz = INITVALUE;
+  CCTK_REAL qp1o144dx2dy = INITVALUE;
+  CCTK_REAL qp1o144dx2dz = INITVALUE;
+  CCTK_REAL qp1o144dxdy = INITVALUE;
+  CCTK_REAL qp1o144dxdy2 = INITVALUE;
+  CCTK_REAL qp1o144dxdz = INITVALUE;
+  CCTK_REAL qp1o144dxdz2 = INITVALUE;
+  CCTK_REAL qp1o144dy2dz = INITVALUE;
+  CCTK_REAL qp1o144dydz = INITVALUE;
+  CCTK_REAL qp1o144dydz2 = INITVALUE;
+  CCTK_REAL qp1o1728dxdydz = INITVALUE;
+  CCTK_REAL qp1o2dx3 = INITVALUE;
+  CCTK_REAL qp1o2dy3 = INITVALUE;
+  CCTK_REAL qp1o2dz3 = INITVALUE;
+  CCTK_REAL qpm1o12dx2 = INITVALUE;
+  CCTK_REAL qpm1o12dy2 = INITVALUE;
+  CCTK_REAL qpm1o12dz2 = INITVALUE;
   
   if (verbose > 1)
   {
@@ -73,15 +83,25 @@ void WT_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[3],
   hdzi = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  p1o12dx = INV(dx)/12.;
-  p1o12dy = INV(dy)/12.;
-  p1o12dz = INV(dz)/12.;
-  p1o144dxdy = (INV(dx)*INV(dy))/144.;
-  p1o144dxdz = (INV(dx)*INV(dz))/144.;
-  p1o144dydz = (INV(dy)*INV(dz))/144.;
-  pm1o12dx2 = -pow(dx,-2)/12.;
-  pm1o12dy2 = -pow(dy,-2)/12.;
-  pm1o12dz2 = -pow(dz,-2)/12.;
+  qp1o12dx = INV(dx)/12.;
+  qp1o12dy = INV(dy)/12.;
+  qp1o12dz = INV(dz)/12.;
+  qp1o144dx2dy = (INV(dy)*pow(dx,-2))/144.;
+  qp1o144dx2dz = (INV(dz)*pow(dx,-2))/144.;
+  qp1o144dxdy = (INV(dx)*INV(dy))/144.;
+  qp1o144dxdy2 = (INV(dx)*pow(dy,-2))/144.;
+  qp1o144dxdz = (INV(dx)*INV(dz))/144.;
+  qp1o144dxdz2 = (INV(dx)*pow(dz,-2))/144.;
+  qp1o144dy2dz = (INV(dz)*pow(dy,-2))/144.;
+  qp1o144dydz = (INV(dy)*INV(dz))/144.;
+  qp1o144dydz2 = (INV(dy)*pow(dz,-2))/144.;
+  qp1o1728dxdydz = (INV(dx)*INV(dy)*INV(dz))/1728.;
+  qp1o2dx3 = khalf*pow(dx,-3);
+  qp1o2dy3 = khalf*pow(dy,-3);
+  qp1o2dz3 = khalf*pow(dz,-3);
+  qpm1o12dx2 = -pow(dx,-2)/12.;
+  qpm1o12dy2 = -pow(dy,-2)/12.;
+  qpm1o12dz2 = -pow(dz,-2)/12.;
   
   /* Loop over the grid points */
   LC_LOOP3 (WT_RHS,