Update generated WaveToy

12 files changed, 220 insertions, 189 deletions

diff --git a/ML_WaveToy/configuration.ccl b/ML_WaveToy/configuration.ccl
index b343426..c4e15c5 100644
--- a/ML_WaveToy/configuration.ccl
+++ b/ML_WaveToy/configuration.ccl
@@ -1,5 +1,6 @@
-# File produced by user diener
-# Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007)
+# File produced by user eschnett
+# Produced with Mathematica Version 6.0 for Mac OS X x86 (64-bit) (May 21, 2008)
 
 # Mathematica script written by Ian Hinder and Sascha Husa
 
+REQUIRES LoopControl
diff --git a/ML_WaveToy/interface.ccl b/ML_WaveToy/interface.ccl
index fae6a52..6037ad5 100644
--- a/ML_WaveToy/interface.ccl
+++ b/ML_WaveToy/interface.ccl
@@ -1,5 +1,5 @@
-# File produced by user diener
-# Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007)
+# File produced by user eschnett
+# Produced with Mathematica Version 6.0 for Mac OS X x86 (64-bit) (May 21, 2008)
 
 # Mathematica script written by Ian Hinder and Sascha Husa
 
@@ -13,13 +13,11 @@ USES INCLUDE: GenericFD.h
 USES INCLUDE: Symmetry.h
 USES INCLUDE: sbp_calc_coeffs.h
 USES INCLUDE: Boundary.h
+USES INCLUDE: loopcontrol.h
 
 CCTK_INT FUNCTION MoLRegisterEvolved(CCTK_INT IN EvolvedIndex, CCTK_INT IN RHSIndex)
 USES FUNCTION MoLRegisterEvolved
 
-CCTK_INT FUNCTION MoLRegisterConstrained(CCTK_INT IN ConstrainedIndex)
-USES FUNCTION MoLRegisterConstrained
-
 SUBROUTINE Diff_coeff(CCTK_POINTER_TO_CONST IN cctkGH, CCTK_INT IN dir, CCTK_INT IN nsize, CCTK_INT OUT ARRAY imin, CCTK_INT OUT ARRAY imax, CCTK_REAL OUT ARRAY q, CCTK_INT IN table_handle)
 USES FUNCTION Diff_coeff
 
diff --git a/ML_WaveToy/param.ccl b/ML_WaveToy/param.ccl
index d341963..745aaa1 100644
--- a/ML_WaveToy/param.ccl
+++ b/ML_WaveToy/param.ccl
@@ -1,5 +1,5 @@
-# File produced by user diener
-# Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007)
+# File produced by user eschnett
+# Produced with Mathematica Version 6.0 for Mac OS X x86 (64-bit) (May 21, 2008)
 
 # Mathematica script written by Ian Hinder and Sascha Husa
 
@@ -16,7 +16,6 @@ USES CCTK_INT boundary_width
 shares: MethodOfLines
 
 USES CCTK_INT MoL_Num_Evolved_Vars
-USES CCTK_INT MoL_Num_Constrained_Vars
 
 restricted:
 CCTK_INT verbose "verbose"
@@ -31,12 +30,6 @@ CCTK_INT ML_WaveToy_MaxNumEvolvedVars "Number of evolved variables used by this
 } 2
 
 restricted:
-CCTK_INT ML_WaveToy_MaxNumConstrainedVars "Number of constrained variables used by this thorn" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Constrained_Vars
-{
-  58:58 :: "Number of constrained variables used by this thorn"
-} 58
-
-restricted:
 CCTK_INT timelevels "Number of active timelevels"
 {
   0:2 :: ""
diff --git a/ML_WaveToy/schedule.ccl b/ML_WaveToy/schedule.ccl
index a02f75f..9b0da13 100644
--- a/ML_WaveToy/schedule.ccl
+++ b/ML_WaveToy/schedule.ccl
@@ -1,5 +1,5 @@
-# File produced by user diener
-# Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007)
+# File produced by user eschnett
+# Produced with Mathematica Version 6.0 for Mac OS X x86 (64-bit) (May 21, 2008)
 
 # Mathematica script written by Ian Hinder and Sascha Husa
 
@@ -38,7 +38,7 @@ schedule ML_WaveToy_RegisterVars in MoL_Register
   OPTIONS: meta
 } "Register Variables for MoL"
 
-schedule ML_WaveToy_RegisterSymmetries at BASEGRID
+schedule ML_WaveToy_RegisterSymmetries in SymmetryRegister
 {
   LANG: C
   OPTIONS: meta
diff --git a/ML_WaveToy/src/Boundaries.c b/ML_WaveToy/src/Boundaries.c
index e119d69..294f430 100644
--- a/ML_WaveToy/src/Boundaries.c
+++ b/ML_WaveToy/src/Boundaries.c
@@ -1,5 +1,5 @@
-/*  File produced by user diener */
-/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
+/*  File produced by user eschnett */
+/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (64-bit) (May 21, 2008) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
@@ -22,16 +22,16 @@
 
 void ML_WaveToy_CheckBoundaries(CCTK_ARGUMENTS)
 {
-  DECLARE_CCTK_ARGUMENTS
-  DECLARE_CCTK_PARAMETERS
+  DECLARE_CCTK_ARGUMENTS;
+  DECLARE_CCTK_PARAMETERS;
   
   return;
 }
 
 void ML_WaveToy_ApplyBoundConds(CCTK_ARGUMENTS)
 {
-  DECLARE_CCTK_ARGUMENTS
-  DECLARE_CCTK_PARAMETERS
+  DECLARE_CCTK_ARGUMENTS;
+  DECLARE_CCTK_PARAMETERS;
   
   CCTK_INT ierr = 0;
   
@@ -43,7 +43,7 @@ void ML_WaveToy_ApplyBoundConds(CCTK_ARGUMENTS)
     ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, 
                       "ML_WaveToy::WT_rho", WT_rho_bound);
     if (ierr < 0)
-       CCTK_WARN(-1, "Failed to register WT_rho_bound BC for ML_WaveToy::WT_rho!");
+       CCTK_WARN(0, "Failed to register WT_rho_bound BC for ML_WaveToy::WT_rho!");
   }
   
   if (CCTK_EQUALS(WT_u_bound, "none"  ) ||
@@ -54,7 +54,7 @@ void ML_WaveToy_ApplyBoundConds(CCTK_ARGUMENTS)
     ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, 
                       "ML_WaveToy::WT_u", WT_u_bound);
     if (ierr < 0)
-       CCTK_WARN(-1, "Failed to register WT_u_bound BC for ML_WaveToy::WT_u!");
+       CCTK_WARN(0, "Failed to register WT_u_bound BC for ML_WaveToy::WT_u!");
   }
   
   if (CCTK_EQUALS(rho_bound, "none"  ) ||
@@ -65,7 +65,7 @@ void ML_WaveToy_ApplyBoundConds(CCTK_ARGUMENTS)
     ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, 
                       "ML_WaveToy::rho", rho_bound);
     if (ierr < 0)
-       CCTK_WARN(-1, "Failed to register rho_bound BC for ML_WaveToy::rho!");
+       CCTK_WARN(0, "Failed to register rho_bound BC for ML_WaveToy::rho!");
   }
   
   if (CCTK_EQUALS(u_bound, "none"  ) ||
@@ -76,142 +76,150 @@ void ML_WaveToy_ApplyBoundConds(CCTK_ARGUMENTS)
     ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, 
                       "ML_WaveToy::u", u_bound);
     if (ierr < 0)
-       CCTK_WARN(-1, "Failed to register u_bound BC for ML_WaveToy::u!");
+       CCTK_WARN(0, "Failed to register u_bound BC for ML_WaveToy::u!");
   }
   
   if (CCTK_EQUALS(WT_rho_bound, "radiative"))
   {
    /* apply radiation boundary condition */
-    CCTK_INT handle_WT_rho_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
-    if (handle_WT_rho_bound < 0) CCTK_WARN(-1, "could not create table!");
+    static CCTK_INT handle_WT_rho_bound = -1;
+    if (handle_WT_rho_bound < 0) handle_WT_rho_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+    if (handle_WT_rho_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_WT_rho_bound , WT_rho_bound_limit, "LIMIT") < 0)
-       CCTK_WARN(-1, "could not set LIMIT value in table!");
+       CCTK_WARN(0, "could not set LIMIT value in table!");
     if (Util_TableSetReal(handle_WT_rho_bound ,WT_rho_bound_speed, "SPEED") < 0)
-       CCTK_WARN(-1, "could not set SPEED value in table!");
+       CCTK_WARN(0, "could not set SPEED value in table!");
   
     ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_WT_rho_bound, 
                       "ML_WaveToy::WT_rho", "Radiation");
   
     if (ierr < 0)
-       CCTK_WARN(-1, "Failed to register Radiation BC for ML_WaveToy::WT_rho!");
+       CCTK_WARN(0, "Failed to register Radiation BC for ML_WaveToy::WT_rho!");
   
   }
   
   if (CCTK_EQUALS(WT_u_bound, "radiative"))
   {
    /* apply radiation boundary condition */
-    CCTK_INT handle_WT_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
-    if (handle_WT_u_bound < 0) CCTK_WARN(-1, "could not create table!");
+    static CCTK_INT handle_WT_u_bound = -1;
+    if (handle_WT_u_bound < 0) handle_WT_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+    if (handle_WT_u_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_WT_u_bound , WT_u_bound_limit, "LIMIT") < 0)
-       CCTK_WARN(-1, "could not set LIMIT value in table!");
+       CCTK_WARN(0, "could not set LIMIT value in table!");
     if (Util_TableSetReal(handle_WT_u_bound ,WT_u_bound_speed, "SPEED") < 0)
-       CCTK_WARN(-1, "could not set SPEED value in table!");
+       CCTK_WARN(0, "could not set SPEED value in table!");
   
     ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_WT_u_bound, 
                       "ML_WaveToy::WT_u", "Radiation");
   
     if (ierr < 0)
-       CCTK_WARN(-1, "Failed to register Radiation BC for ML_WaveToy::WT_u!");
+       CCTK_WARN(0, "Failed to register Radiation BC for ML_WaveToy::WT_u!");
   
   }
   
   if (CCTK_EQUALS(rho_bound, "radiative"))
   {
    /* apply radiation boundary condition */
-    CCTK_INT handle_rho_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
-    if (handle_rho_bound < 0) CCTK_WARN(-1, "could not create table!");
+    static CCTK_INT handle_rho_bound = -1;
+    if (handle_rho_bound < 0) handle_rho_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+    if (handle_rho_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_rho_bound , rho_bound_limit, "LIMIT") < 0)
-       CCTK_WARN(-1, "could not set LIMIT value in table!");
+       CCTK_WARN(0, "could not set LIMIT value in table!");
     if (Util_TableSetReal(handle_rho_bound ,rho_bound_speed, "SPEED") < 0)
-        CCTK_WARN(-1, "could not set SPEED value in table!");
+        CCTK_WARN(0, "could not set SPEED value in table!");
   
     ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_rho_bound, 
                       "ML_WaveToy::rho", "Radiation");
   
     if (ierr < 0)
-       CCTK_WARN(-1, "Failed to register Radiation BC for ML_WaveToy::rho!");
+       CCTK_WARN(0, "Failed to register Radiation BC for ML_WaveToy::rho!");
   
   }
   
   if (CCTK_EQUALS(u_bound, "radiative"))
   {
    /* apply radiation boundary condition */
-    CCTK_INT handle_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
-    if (handle_u_bound < 0) CCTK_WARN(-1, "could not create table!");
+    static CCTK_INT handle_u_bound = -1;
+    if (handle_u_bound < 0) handle_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+    if (handle_u_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_u_bound , u_bound_limit, "LIMIT") < 0)
-       CCTK_WARN(-1, "could not set LIMIT value in table!");
+       CCTK_WARN(0, "could not set LIMIT value in table!");
     if (Util_TableSetReal(handle_u_bound ,u_bound_speed, "SPEED") < 0)
-        CCTK_WARN(-1, "could not set SPEED value in table!");
+        CCTK_WARN(0, "could not set SPEED value in table!");
   
     ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_u_bound, 
                       "ML_WaveToy::u", "Radiation");
   
     if (ierr < 0)
-       CCTK_WARN(-1, "Failed to register Radiation BC for ML_WaveToy::u!");
+       CCTK_WARN(0, "Failed to register Radiation BC for ML_WaveToy::u!");
   
   }
   
   if (CCTK_EQUALS(WT_rho_bound, "scalar"))
   {
    /* apply scalar boundary condition */
-    CCTK_INT handle_WT_rho_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
-    if (handle_WT_rho_bound < 0) CCTK_WARN(-1, "could not create table!");
+    static CCTK_INT handle_WT_rho_bound = -1;
+    if (handle_WT_rho_bound < 0) handle_WT_rho_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+    if (handle_WT_rho_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_WT_rho_bound ,WT_rho_bound_scalar, "SCALAR") < 0)
-        CCTK_WARN(-1, "could not set SCALAR value in table!");
+        CCTK_WARN(0, "could not set SCALAR value in table!");
   
     ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_WT_rho_bound, 
                       "ML_WaveToy::WT_rho", "scalar");
   
     if (ierr < 0)
-       CCTK_WARN(-1, "Failed to register Scalar BC for ML_WaveToy::WT_rho!");
+       CCTK_WARN(0, "Failed to register Scalar BC for ML_WaveToy::WT_rho!");
   
   }
   
   if (CCTK_EQUALS(WT_u_bound, "scalar"))
   {
    /* apply scalar boundary condition */
-    CCTK_INT handle_WT_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
-    if (handle_WT_u_bound < 0) CCTK_WARN(-1, "could not create table!");
+    static CCTK_INT handle_WT_u_bound = -1;
+    if (handle_WT_u_bound < 0) handle_WT_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+    if (handle_WT_u_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_WT_u_bound ,WT_u_bound_scalar, "SCALAR") < 0)
-        CCTK_WARN(-1, "could not set SCALAR value in table!");
+        CCTK_WARN(0, "could not set SCALAR value in table!");
   
     ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_WT_u_bound, 
                       "ML_WaveToy::WT_u", "scalar");
   
     if (ierr < 0)
-       CCTK_WARN(-1, "Failed to register Scalar BC for ML_WaveToy::WT_u!");
+       CCTK_WARN(0, "Failed to register Scalar BC for ML_WaveToy::WT_u!");
   
   }
   
   if (CCTK_EQUALS(rho_bound, "scalar"))
   {
    /* apply scalar boundary condition */
-    CCTK_INT handle_rho_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
-    if (handle_rho_bound < 0) CCTK_WARN(-1, "could not create table!");
+    static CCTK_INT handle_rho_bound = -1;
+    if (handle_rho_bound < 0) handle_rho_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+    if (handle_rho_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_rho_bound ,rho_bound_scalar, "SCALAR") < 0)
-      CCTK_WARN(-1, "could not set SCALAR value in table!");
+      CCTK_WARN(0, "could not set SCALAR value in table!");
   
     ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_rho_bound, 
                       "ML_WaveToy::rho", "scalar");
   
     if (ierr < 0)
-       CCTK_WARN(-1, "Error in registering Scalar BC for ML_WaveToy::rho!");
+       CCTK_WARN(0, "Error in registering Scalar BC for ML_WaveToy::rho!");
   
   }
   
   if (CCTK_EQUALS(u_bound, "scalar"))
   {
    /* apply scalar boundary condition */
-    CCTK_INT handle_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
-    if (handle_u_bound < 0) CCTK_WARN(-1, "could not create table!");
+    static CCTK_INT handle_u_bound = -1;
+    if (handle_u_bound < 0) handle_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+    if (handle_u_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_u_bound ,u_bound_scalar, "SCALAR") < 0)
-      CCTK_WARN(-1, "could not set SCALAR value in table!");
+      CCTK_WARN(0, "could not set SCALAR value in table!");
   
     ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_u_bound, 
                       "ML_WaveToy::u", "scalar");
   
     if (ierr < 0)
-       CCTK_WARN(-1, "Error in registering Scalar BC for ML_WaveToy::u!");
+       CCTK_WARN(0, "Error in registering Scalar BC for ML_WaveToy::u!");
   
   }
   return;
diff --git a/ML_WaveToy/src/Differencing.h b/ML_WaveToy/src/Differencing.h
new file mode 100644
index 0000000..fd9befd
--- /dev/null
+++ b/ML_WaveToy/src/Differencing.h
@@ -0,0 +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)]))
diff --git a/ML_WaveToy/src/RegisterMoL.c b/ML_WaveToy/src/RegisterMoL.c
index 86ccf50..7d7a356 100644
--- a/ML_WaveToy/src/RegisterMoL.c
+++ b/ML_WaveToy/src/RegisterMoL.c
@@ -1,5 +1,5 @@
-/*  File produced by user diener */
-/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
+/*  File produced by user eschnett */
+/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (64-bit) (May 21, 2008) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
@@ -9,8 +9,8 @@
 
 void ML_WaveToy_RegisterVars(CCTK_ARGUMENTS)
 {
-  DECLARE_CCTK_ARGUMENTS
-  DECLARE_CCTK_PARAMETERS
+  DECLARE_CCTK_ARGUMENTS;
+  DECLARE_CCTK_PARAMETERS;
   
   CCTK_INT ierr = 0;
   
diff --git a/ML_WaveToy/src/RegisterSymmetries.c b/ML_WaveToy/src/RegisterSymmetries.c
index 82a4f12..77696e4 100644
--- a/ML_WaveToy/src/RegisterSymmetries.c
+++ b/ML_WaveToy/src/RegisterSymmetries.c
@@ -1,5 +1,5 @@
-/*  File produced by user diener */
-/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
+/*  File produced by user eschnett */
+/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (64-bit) (May 21, 2008) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
@@ -10,8 +10,8 @@
 
 void ML_WaveToy_RegisterSymmetries(CCTK_ARGUMENTS)
 {
-  DECLARE_CCTK_ARGUMENTS
-  DECLARE_CCTK_PARAMETERS
+  DECLARE_CCTK_ARGUMENTS;
+  DECLARE_CCTK_PARAMETERS;
   
   
   /* array holding symmetry definitions */
diff --git a/ML_WaveToy/src/Startup.c b/ML_WaveToy/src/Startup.c
index 49b9a0b..7016d26 100644
--- a/ML_WaveToy/src/Startup.c
+++ b/ML_WaveToy/src/Startup.c
@@ -1,5 +1,5 @@
-/*  File produced by user diener */
-/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
+/*  File produced by user eschnett */
+/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (64-bit) (May 21, 2008) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
diff --git a/ML_WaveToy/src/WT_Gaussian.c b/ML_WaveToy/src/WT_Gaussian.c
index 075b622..bec7489 100644
--- a/ML_WaveToy/src/WT_Gaussian.c
+++ b/ML_WaveToy/src/WT_Gaussian.c
@@ -1,5 +1,5 @@
-/*  File produced by user diener */
-/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
+/*  File produced by user eschnett */
+/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (64-bit) (May 21, 2008) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
@@ -14,6 +14,7 @@
 #include "cctk_Parameters.h"
 #include "GenericFD.h"
 #include "Differencing.h"
+#include "loopcontrol.h"
 
 /* Define macros used in calculations */
 #define INITVALUE  (42)
@@ -22,17 +23,12 @@
 #define CUB(x) ((x) * (x) * (x))
 #define QAD(x) ((x) * (x) * (x) * (x))
 
-void WT_Gaussian_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[3], CCTK_REAL tangentA[3], CCTK_REAL tangentB[3], CCTK_INT min[3], CCTK_INT max[3], CCTK_INT n_subblock_gfs, CCTK_REAL *subblock_gfs[])
+void WT_Gaussian_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT const face, CCTK_REAL const normal[3], CCTK_REAL const tangentA[3], CCTK_REAL const tangentB[3], CCTK_INT const min[3], CCTK_INT const max[3], CCTK_INT const n_subblock_gfs, CCTK_REAL * const subblock_gfs[])
 {
-  DECLARE_CCTK_ARGUMENTS
-  DECLARE_CCTK_PARAMETERS
+  DECLARE_CCTK_ARGUMENTS;
+  DECLARE_CCTK_PARAMETERS;
   
   
-  /* Declare the variables used for looping over grid points */
-  CCTK_INT i = INITVALUE, j = INITVALUE, k = INITVALUE;
-  CCTK_INT index = INITVALUE;
-  CCTK_INT subblock_index = INITVALUE;
-  
   /* Declare finite differencing variables */
   CCTK_REAL dx = INITVALUE, dy = INITVALUE, dz = INITVALUE;
   CCTK_REAL dxi = INITVALUE, dyi = INITVALUE, dzi = INITVALUE;
@@ -111,54 +107,54 @@ 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 */
-  for (k = min[2]; k < max[2]; k++)
+  #pragma omp parallel
+  LC_LOOP3 (WT_Gaussian,
+            i,j,k, min[0],min[1],min[2], max[0],max[1],max[2],
+            cctk_lssh[CCTK_LSSH_IDX(0,0)],cctk_lssh[CCTK_LSSH_IDX(0,1)],cctk_lssh[CCTK_LSSH_IDX(0,2)])
   {
-    for (j = min[1]; j < max[1]; j++)
-    {
-      for (i = min[0]; i < max[0]; i++)
-      {
-         index  =  CCTK_GFINDEX3D(cctkGH,i,j,k) ;
-         subblock_index  =  i - min[0] + (max[0] - min[0]) * (j - min[1] + (max[1]-min[1]) * (k - min[2])) ;
-        
-        /* Declare shorthands */
-        
-        /* Declare local copies of grid functions */
-        CCTK_REAL rhoL = INITVALUE;
-        CCTK_REAL uL = INITVALUE;
-        /* Declare precomputed derivatives*/
-        
-        /* Declare derivatives */
-        
-        /* Assign local copies of grid functions */
-        
-        /* Assign local copies of subblock grid functions */
-        
-        /* Include user supplied include files */
-        
-        /* Precompute derivatives (new style) */
-        
-        /* Precompute derivatives (old style) */
-        
-        /* Calculate temporaries and grid functions */
-        uL  =  0;
-        
-        rhoL  =  0;
-        
-        
-        /* Copy local copies back to grid functions */
-        rho[index] = rhoL;
-        u[index] = uL;
-        
-        /* Copy local copies back to subblock grid functions */
-      }
-    }
+    int index = INITVALUE;
+    int subblock_index = INITVALUE;
+    index = CCTK_GFINDEX3D(cctkGH,i,j,k);
+    subblock_index = i - min[0] + (max[0] - min[0]) * (j - min[1] + (max[1]-min[1]) * (k - min[2]));
+    
+    /* Declare shorthands */
+    
+    /* Declare local copies of grid functions */
+    CCTK_REAL rhoL = INITVALUE;
+    CCTK_REAL uL = INITVALUE;
+    /* Declare precomputed derivatives*/
+    
+    /* Declare derivatives */
+    
+    /* Assign local copies of grid functions */
+    
+    /* Assign local copies of subblock grid functions */
+    
+    /* Include user supplied include files */
+    
+    /* Precompute derivatives (new style) */
+    
+    /* Precompute derivatives (old style) */
+    
+    /* Calculate temporaries and grid functions */
+    uL  =  0;
+    
+    rhoL  =  0;
+    
+    
+    /* Copy local copies back to grid functions */
+    rho[index] = rhoL;
+    u[index] = uL;
+    
+    /* Copy local copies back to subblock grid functions */
   }
+  LC_ENDLOOP3 (WT_Gaussian);
 }
 
 void WT_Gaussian(CCTK_ARGUMENTS)
 {
-  DECLARE_CCTK_ARGUMENTS
-  DECLARE_CCTK_PARAMETERS
+  DECLARE_CCTK_ARGUMENTS;
+  DECLARE_CCTK_PARAMETERS;
   
   GenericFD_LoopOverEverything(cctkGH, &WT_Gaussian_Body);
 }
diff --git a/ML_WaveToy/src/WT_RHS.c b/ML_WaveToy/src/WT_RHS.c
index 2409a58..2b96bb3 100644
--- a/ML_WaveToy/src/WT_RHS.c
+++ b/ML_WaveToy/src/WT_RHS.c
@@ -1,5 +1,5 @@
-/*  File produced by user diener */
-/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
+/*  File produced by user eschnett */
+/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (64-bit) (May 21, 2008) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
@@ -14,6 +14,7 @@
 #include "cctk_Parameters.h"
 #include "GenericFD.h"
 #include "Differencing.h"
+#include "loopcontrol.h"
 
 /* Define macros used in calculations */
 #define INITVALUE  (42)
@@ -22,17 +23,12 @@
 #define CUB(x) ((x) * (x) * (x))
 #define QAD(x) ((x) * (x) * (x) * (x))
 
-void WT_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[3], CCTK_REAL tangentA[3], CCTK_REAL tangentB[3], CCTK_INT min[3], CCTK_INT max[3], CCTK_INT n_subblock_gfs, CCTK_REAL *subblock_gfs[])
+void WT_RHS_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT const face, CCTK_REAL const normal[3], CCTK_REAL const tangentA[3], CCTK_REAL const tangentB[3], CCTK_INT const min[3], CCTK_INT const max[3], CCTK_INT const n_subblock_gfs, CCTK_REAL * const subblock_gfs[])
 {
-  DECLARE_CCTK_ARGUMENTS
-  DECLARE_CCTK_PARAMETERS
+  DECLARE_CCTK_ARGUMENTS;
+  DECLARE_CCTK_PARAMETERS;
   
   
-  /* Declare the variables used for looping over grid points */
-  CCTK_INT i = INITVALUE, j = INITVALUE, k = INITVALUE;
-  CCTK_INT index = INITVALUE;
-  CCTK_INT subblock_index = INITVALUE;
-  
   /* Declare finite differencing variables */
   CCTK_REAL dx = INITVALUE, dy = INITVALUE, dz = INITVALUE;
   CCTK_REAL dxi = INITVALUE, dyi = INITVALUE, dzi = INITVALUE;
@@ -111,62 +107,62 @@ 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 */
-  for (k = min[2]; k < max[2]; k++)
+  #pragma omp parallel
+  LC_LOOP3 (WT_RHS,
+            i,j,k, min[0],min[1],min[2], max[0],max[1],max[2],
+            cctk_lssh[CCTK_LSSH_IDX(0,0)],cctk_lssh[CCTK_LSSH_IDX(0,1)],cctk_lssh[CCTK_LSSH_IDX(0,2)])
   {
-    for (j = min[1]; j < max[1]; j++)
-    {
-      for (i = min[0]; i < max[0]; i++)
-      {
-         index  =  CCTK_GFINDEX3D(cctkGH,i,j,k) ;
-         subblock_index  =  i - min[0] + (max[0] - min[0]) * (j - min[1] + (max[1]-min[1]) * (k - min[2])) ;
-        
-        /* Declare shorthands */
-        
-        /* Declare local copies of grid functions */
-        CCTK_REAL rhoL = INITVALUE, rhorhsL = INITVALUE;
-        CCTK_REAL uL = INITVALUE, urhsL = INITVALUE;
-        /* Declare precomputed derivatives*/
-        
-        /* Declare derivatives */
-        CCTK_REAL PDstandardNth11u = INITVALUE;
-        CCTK_REAL PDstandardNth22u = INITVALUE;
-        CCTK_REAL PDstandardNth33u = INITVALUE;
-        
-        /* Assign local copies of grid functions */
-        rhoL = rho[index];
-        uL = u[index];
-        
-        /* Assign local copies of subblock grid functions */
-        
-        /* Include user supplied include files */
-        
-        /* Precompute derivatives (new style) */
-        PDstandardNth11u = PDstandardNth11(u, i, j, k);
-        PDstandardNth22u = PDstandardNth22(u, i, j, k);
-        PDstandardNth33u = PDstandardNth33(u, i, j, k);
-        
-        /* Precompute derivatives (old style) */
-        
-        /* Calculate temporaries and grid functions */
-        urhsL  =  rhoL;
-        
-        rhorhsL  =  PDstandardNth11u + PDstandardNth22u + PDstandardNth33u;
-        
-        
-        /* Copy local copies back to grid functions */
-        rhorhs[index] = rhorhsL;
-        urhs[index] = urhsL;
-        
-        /* Copy local copies back to subblock grid functions */
-      }
-    }
+    int index = INITVALUE;
+    int subblock_index = INITVALUE;
+    index = CCTK_GFINDEX3D(cctkGH,i,j,k);
+    subblock_index = i - min[0] + (max[0] - min[0]) * (j - min[1] + (max[1]-min[1]) * (k - min[2]));
+    
+    /* Declare shorthands */
+    
+    /* Declare local copies of grid functions */
+    CCTK_REAL rhoL = INITVALUE, rhorhsL = INITVALUE;
+    CCTK_REAL uL = INITVALUE, urhsL = INITVALUE;
+    /* Declare precomputed derivatives*/
+    
+    /* Declare derivatives */
+    CCTK_REAL PDstandardNth11u = INITVALUE;
+    CCTK_REAL PDstandardNth22u = INITVALUE;
+    CCTK_REAL PDstandardNth33u = INITVALUE;
+    
+    /* Assign local copies of grid functions */
+    rhoL = rho[index];
+    uL = u[index];
+    
+    /* Assign local copies of subblock grid functions */
+    
+    /* Include user supplied include files */
+    
+    /* Precompute derivatives (new style) */
+    PDstandardNth11u = PDstandardNth11(u, i, j, k);
+    PDstandardNth22u = PDstandardNth22(u, i, j, k);
+    PDstandardNth33u = PDstandardNth33(u, i, j, k);
+    
+    /* Precompute derivatives (old style) */
+    
+    /* Calculate temporaries and grid functions */
+    urhsL  =  rhoL;
+    
+    rhorhsL  =  PDstandardNth11u + PDstandardNth22u + PDstandardNth33u;
+    
+    
+    /* Copy local copies back to grid functions */
+    rhorhs[index] = rhorhsL;
+    urhs[index] = urhsL;
+    
+    /* Copy local copies back to subblock grid functions */
   }
+  LC_ENDLOOP3 (WT_RHS);
 }
 
 void WT_RHS(CCTK_ARGUMENTS)
 {
-  DECLARE_CCTK_ARGUMENTS
-  DECLARE_CCTK_PARAMETERS
+  DECLARE_CCTK_ARGUMENTS;
+  DECLARE_CCTK_PARAMETERS;
   
   GenericFD_LoopOverInterior(cctkGH, &WT_RHS_Body);
 }
diff --git a/ML_WaveToy/src/make.code.defn b/ML_WaveToy/src/make.code.defn
index 0fa6a98..e88f77e 100644
--- a/ML_WaveToy/src/make.code.defn
+++ b/ML_WaveToy/src/make.code.defn
@@ -1,5 +1,5 @@
-# File produced by user diener
-# Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007)
+# File produced by user eschnett
+# Produced with Mathematica Version 6.0 for Mac OS X x86 (64-bit) (May 21, 2008)
 
 # Mathematica script written by Ian Hinder and Sascha Husa