diff options
| author | Erik Schnetter <schnetter@cct.lsu.edu> | 2007-12-27 23:44:35 +0100 |
|---|---|---|
| committer | Erik Schnetter <schnetter@cct.lsu.edu> | 2007-12-27 23:44:35 +0100 |
| commit | 51c2468863e407195bd99c852071bdda74a889dd (patch) | |
| tree | 89b95b8d12a35614019133617be1d6bbc3489c8e /ML_WaveToy/src | |
| parent | 684fd4e0c6dd2c101967710e26c4281f2041396b (diff) | |
Regenerate code after removing leftover debugging code from Kranc
Diffstat (limited to 'ML_WaveToy/src')
| -rw-r--r-- | ML_WaveToy/src/Boundaries.c | 4 | ||||
| -rw-r--r-- | ML_WaveToy/src/Differencing.h | 78 | ||||
| -rw-r--r-- | ML_WaveToy/src/RegisterMoL.c | 4 | ||||
| -rw-r--r-- | ML_WaveToy/src/RegisterSymmetries.c | 4 | ||||
| -rw-r--r-- | ML_WaveToy/src/Startup.c | 4 | ||||
| -rw-r--r-- | ML_WaveToy/src/WT_Gaussian.c | 80 | ||||
| -rw-r--r-- | ML_WaveToy/src/WT_RHS.c | 80 | ||||
| -rw-r--r-- | ML_WaveToy/src/make.code.defn | 4 |
8 files changed, 129 insertions, 129 deletions
diff --git a/ML_WaveToy/src/Boundaries.c b/ML_WaveToy/src/Boundaries.c index e119d69..563f55d 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 (32-bit) (April 20, 2007) */ /* Mathematica script written by Ian Hinder and Sascha Husa */ diff --git a/ML_WaveToy/src/Differencing.h b/ML_WaveToy/src/Differencing.h index 6b60a8f..626c54f 100644 --- a/ML_WaveToy/src/Differencing.h +++ b/ML_WaveToy/src/Differencing.h @@ -1,39 +1,39 @@ -#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)])) +#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..b394a5d 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 (32-bit) (April 20, 2007) */ /* Mathematica script written by Ian Hinder and Sascha Husa */ diff --git a/ML_WaveToy/src/RegisterSymmetries.c b/ML_WaveToy/src/RegisterSymmetries.c index 82a4f12..cf9633e 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 (32-bit) (April 20, 2007) */ /* Mathematica script written by Ian Hinder and Sascha Husa */ diff --git a/ML_WaveToy/src/Startup.c b/ML_WaveToy/src/Startup.c index 49b9a0b..71be71f 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 (32-bit) (April 20, 2007) */ /* 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 aabb376..bbba1a1 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 (32-bit) (April 20, 2007) */ /* Mathematica script written by Ian Hinder and Sascha Husa */ @@ -34,25 +34,25 @@ void WT_Gaussian_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal /* Declare predefined quantities */ - 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; + CCTK_REAL p1o12dx = INITVALUE; + CCTK_REAL p1o12dy = INITVALUE; + CCTK_REAL p1o12dz = INITVALUE; + CCTK_REAL p1o144dx2dy = INITVALUE; + CCTK_REAL p1o144dx2dz = INITVALUE; + CCTK_REAL p1o144dxdy = INITVALUE; + CCTK_REAL p1o144dxdy2 = INITVALUE; + CCTK_REAL p1o144dxdz = INITVALUE; + CCTK_REAL p1o144dxdz2 = INITVALUE; + CCTK_REAL p1o144dy2dz = INITVALUE; + CCTK_REAL p1o144dydz = INITVALUE; + CCTK_REAL p1o144dydz2 = INITVALUE; + CCTK_REAL p1o1728dxdydz = INITVALUE; + CCTK_REAL p1o2dx3 = INITVALUE; + CCTK_REAL p1o2dy3 = INITVALUE; + CCTK_REAL p1o2dz3 = INITVALUE; + CCTK_REAL pm1o12dx2 = INITVALUE; + CCTK_REAL pm1o12dy2 = INITVALUE; + CCTK_REAL pm1o12dz2 = INITVALUE; if (verbose > 1) { @@ -83,25 +83,25 @@ void WT_Gaussian_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal hdzi = 0.5 * dzi; /* Initialize predefined quantities */ - 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.; + p1o12dx = INV(dx)/12.; + p1o12dy = INV(dy)/12.; + p1o12dz = INV(dz)/12.; + p1o144dx2dy = (INV(dy)*pow(dx,-2))/144.; + p1o144dx2dz = (INV(dz)*pow(dx,-2))/144.; + p1o144dxdy = (INV(dx)*INV(dy))/144.; + p1o144dxdy2 = (INV(dx)*pow(dy,-2))/144.; + p1o144dxdz = (INV(dx)*INV(dz))/144.; + p1o144dxdz2 = (INV(dx)*pow(dz,-2))/144.; + p1o144dy2dz = (INV(dz)*pow(dy,-2))/144.; + p1o144dydz = (INV(dy)*INV(dz))/144.; + p1o144dydz2 = (INV(dy)*pow(dz,-2))/144.; + p1o1728dxdydz = (INV(dx)*INV(dy)*INV(dz))/1728.; + p1o2dx3 = khalf*pow(dx,-3); + p1o2dy3 = khalf*pow(dy,-3); + p1o2dz3 = khalf*pow(dz,-3); + pm1o12dx2 = -pow(dx,-2)/12.; + pm1o12dy2 = -pow(dy,-2)/12.; + pm1o12dz2 = -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 2425966..f06f42a 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 (32-bit) (April 20, 2007) */ /* Mathematica script written by Ian Hinder and Sascha Husa */ @@ -34,25 +34,25 @@ void WT_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[3], /* Declare predefined quantities */ - 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; + CCTK_REAL p1o12dx = INITVALUE; + CCTK_REAL p1o12dy = INITVALUE; + CCTK_REAL p1o12dz = INITVALUE; + CCTK_REAL p1o144dx2dy = INITVALUE; + CCTK_REAL p1o144dx2dz = INITVALUE; + CCTK_REAL p1o144dxdy = INITVALUE; + CCTK_REAL p1o144dxdy2 = INITVALUE; + CCTK_REAL p1o144dxdz = INITVALUE; + CCTK_REAL p1o144dxdz2 = INITVALUE; + CCTK_REAL p1o144dy2dz = INITVALUE; + CCTK_REAL p1o144dydz = INITVALUE; + CCTK_REAL p1o144dydz2 = INITVALUE; + CCTK_REAL p1o1728dxdydz = INITVALUE; + CCTK_REAL p1o2dx3 = INITVALUE; + CCTK_REAL p1o2dy3 = INITVALUE; + CCTK_REAL p1o2dz3 = INITVALUE; + CCTK_REAL pm1o12dx2 = INITVALUE; + CCTK_REAL pm1o12dy2 = INITVALUE; + CCTK_REAL pm1o12dz2 = INITVALUE; if (verbose > 1) { @@ -83,25 +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 */ - 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.; + p1o12dx = INV(dx)/12.; + p1o12dy = INV(dy)/12.; + p1o12dz = INV(dz)/12.; + p1o144dx2dy = (INV(dy)*pow(dx,-2))/144.; + p1o144dx2dz = (INV(dz)*pow(dx,-2))/144.; + p1o144dxdy = (INV(dx)*INV(dy))/144.; + p1o144dxdy2 = (INV(dx)*pow(dy,-2))/144.; + p1o144dxdz = (INV(dx)*INV(dz))/144.; + p1o144dxdz2 = (INV(dx)*pow(dz,-2))/144.; + p1o144dy2dz = (INV(dz)*pow(dy,-2))/144.; + p1o144dydz = (INV(dy)*INV(dz))/144.; + p1o144dydz2 = (INV(dy)*pow(dz,-2))/144.; + p1o1728dxdydz = (INV(dx)*INV(dy)*INV(dz))/1728.; + p1o2dx3 = khalf*pow(dx,-3); + p1o2dy3 = khalf*pow(dy,-3); + p1o2dz3 = khalf*pow(dz,-3); + pm1o12dx2 = -pow(dx,-2)/12.; + pm1o12dy2 = -pow(dy,-2)/12.; + pm1o12dz2 = -pow(dz,-2)/12.; /* Loop over the grid points */ LC_LOOP3 (WT_RHS, diff --git a/ML_WaveToy/src/make.code.defn b/ML_WaveToy/src/make.code.defn index 0fa6a98..1cd6049 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 (32-bit) (April 20, 2007) # Mathematica script written by Ian Hinder and Sascha Husa |