diff options
author | Erik Schnetter <schnetter@cct.lsu.edu> | 2008-01-02 15:20:16 +0100 |
---|---|---|
committer | Erik Schnetter <schnetter@cct.lsu.edu> | 2008-01-02 15:20:16 +0100 |
commit | 09fff22bc4f2e1c7115822f86f859230ff22b85d (patch) | |
tree | 312a17e792189f222c49d5ad56aff0fafc483880 /ML_WaveToy/src | |
parent | aa0b728b9fa5c76f81f2c981be18a9a233246ed4 (diff) |
Update generated code
Update the generated code after the recent changes to Kranc.
Diffstat (limited to 'ML_WaveToy/src')
-rw-r--r-- | ML_WaveToy/src/Differencing.h | 78 | ||||
-rw-r--r-- | ML_WaveToy/src/WT_Gaussian.c | 1 | ||||
-rw-r--r-- | ML_WaveToy/src/WT_RHS.c | 1 |
3 files changed, 41 insertions, 39 deletions
diff --git a/ML_WaveToy/src/Differencing.h b/ML_WaveToy/src/Differencing.h index 626c54f..fd9befd 100644 --- a/ML_WaveToy/src/Differencing.h +++ b/ML_WaveToy/src/Differencing.h @@ -1,39 +1,39 @@ -#define PDstandardNth1(u,i,j,k) (p1o12dx*(-8*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] + 8*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) -#define PDstandardNth2(u,i,j,k) (p1o12dy*(-8*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] + 8*u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)])) -#define PDstandardNth3(u,i,j,k) (p1o12dz*(-8*u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] + 8*u[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)])) -#define PDstandardNth11(u,i,j,k) (pm1o12dx2*(-16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,j,k)])) -#define PDstandardNth22(u,i,j,k) (pm1o12dy2*(-16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)]) + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,j,k)])) -#define PDstandardNth33(u,i,j,k) (pm1o12dz2*(-16*(u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,j,k)])) -#define PDstandardNth12(u,i,j,k) (p1o144dxdy*(-64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)])) -#define PDstandardNth13(u,i,j,k) (p1o144dxdz*(-64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)])) -#define PDstandardNth21(u,i,j,k) (p1o144dxdy*(-64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)])) -#define PDstandardNth23(u,i,j,k) (p1o144dydz*(-64*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)])) -#define PDstandardNth31(u,i,j,k) (p1o144dxdz*(-64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)])) -#define PDstandardNth32(u,i,j,k) (p1o144dydz*(-64*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)])) -#define PDstandardNth111(u,i,j,k) (p1o2dx3*(2*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 2*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) -#define PDstandardNth222(u,i,j,k) (p1o2dy3*(2*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 2*u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)])) -#define PDstandardNth333(u,i,j,k) (p1o2dz3*(2*u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 2*u[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)])) -#define PDstandardNth112(u,i,j,k) (p1o144dx2dy*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)] + 30*(8*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 8*u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))) -#define PDstandardNth113(u,i,j,k) (p1o144dx2dz*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)] + 30*(8*u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 8*u[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))) -#define PDstandardNth221(u,i,j,k) (p1o144dxdy2*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)] - 30*u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) -#define PDstandardNth223(u,i,j,k) (p1o144dy2dz*(-128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)] + 30*(8*u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 8*u[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))) -#define PDstandardNth331(u,i,j,k) (p1o144dxdz2*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)] - 30*u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) -#define PDstandardNth332(u,i,j,k) (p1o144dydz2*(-128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)] - 30*u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)])) -#define PDstandardNth121(u,i,j,k) (p1o144dx2dy*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)] + 30*(8*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 8*u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))) -#define PDstandardNth131(u,i,j,k) (p1o144dx2dz*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)] + 30*(8*u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 8*u[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))) -#define PDstandardNth212(u,i,j,k) (p1o144dxdy2*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)] - 30*u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) -#define PDstandardNth232(u,i,j,k) (p1o144dy2dz*(-128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)] + 30*(8*u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 8*u[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))) -#define PDstandardNth313(u,i,j,k) (p1o144dxdz2*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)] - 30*u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) -#define PDstandardNth323(u,i,j,k) (p1o144dydz2*(-128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)] - 30*u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)])) -#define PDstandardNth211(u,i,j,k) (p1o144dx2dy*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)] + 30*(8*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 8*u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))) -#define PDstandardNth311(u,i,j,k) (p1o144dx2dz*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)] + 30*(8*u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 8*u[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))) -#define PDstandardNth122(u,i,j,k) (p1o144dxdy2*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)] - 30*u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) -#define PDstandardNth322(u,i,j,k) (p1o144dy2dz*(-128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)] + 30*(8*u[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 8*u[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))) -#define PDstandardNth133(u,i,j,k) (p1o144dxdz2*(-128*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)] - 30*u[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) -#define PDstandardNth233(u,i,j,k) (p1o144dydz2*(-128*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)]) + 128*(u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 240*u[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 240*u[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + 16*(u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)]) - 16*(u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)] - 30*u[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + 30*u[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)])) -#define PDstandardNth123(u,i,j,k) (p1o1728dxdydz*(-512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)])) -#define PDstandardNth132(u,i,j,k) (p1o1728dxdydz*(-512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)])) -#define PDstandardNth213(u,i,j,k) (p1o1728dxdydz*(-512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)])) -#define PDstandardNth231(u,i,j,k) (p1o1728dxdydz*(-512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)])) -#define PDstandardNth312(u,i,j,k) (p1o1728dxdydz*(-512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)])) -#define PDstandardNth321(u,i,j,k) (p1o1728dxdydz*(-512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*(u[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - u[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)])) +#define PDstandardNth1(u,i,j,k) (p1o12dx*(-8*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] + 8*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) +#define PDstandardNth2(u,i,j,k) (p1o12dy*(-8*(u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] + 8*(u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)])) +#define PDstandardNth3(u,i,j,k) (p1o12dz*(-8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] + 8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)])) +#define PDstandardNth11(u,i,j,k) (pm1o12dx2*(-16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,i,j,k)])) +#define PDstandardNth22(u,i,j,k) (pm1o12dy2*(-16*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)]) + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,i,j,k)])) +#define PDstandardNth33(u,i,j,k) (pm1o12dz2*(-16*((u)[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,i,j,k)])) +#define PDstandardNth12(u,i,j,k) (p1o144dxdy*(-64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)])) +#define PDstandardNth13(u,i,j,k) (p1o144dxdz*(-64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)])) +#define PDstandardNth21(u,i,j,k) (p1o144dxdy*(-64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)])) +#define PDstandardNth23(u,i,j,k) (p1o144dydz*(-64*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)])) +#define PDstandardNth31(u,i,j,k) (p1o144dxdz*(-64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)])) +#define PDstandardNth32(u,i,j,k) (p1o144dydz*(-64*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)])) +#define PDstandardNth111(u,i,j,k) (p1o2dx3*(2*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 2*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) +#define PDstandardNth222(u,i,j,k) (p1o2dy3*(2*(u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 2*(u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)])) +#define PDstandardNth333(u,i,j,k) (p1o2dz3*(2*(u)[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 2*(u)[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)])) +#define PDstandardNth112(u,i,j,k) (p1o144dx2dy*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)] + 30*(8*(u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 8*(u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))) +#define PDstandardNth113(u,i,j,k) (p1o144dx2dz*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)] + 30*(8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))) +#define PDstandardNth221(u,i,j,k) (p1o144dxdy2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) +#define PDstandardNth223(u,i,j,k) (p1o144dy2dz*(-128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)] + 30*(8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))) +#define PDstandardNth331(u,i,j,k) (p1o144dxdz2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) +#define PDstandardNth332(u,i,j,k) (p1o144dydz2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)])) +#define PDstandardNth121(u,i,j,k) (p1o144dx2dy*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)] + 30*(8*(u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 8*(u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))) +#define PDstandardNth131(u,i,j,k) (p1o144dx2dz*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)] + 30*(8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))) +#define PDstandardNth212(u,i,j,k) (p1o144dxdy2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) +#define PDstandardNth232(u,i,j,k) (p1o144dy2dz*(-128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)] + 30*(8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))) +#define PDstandardNth313(u,i,j,k) (p1o144dxdz2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) +#define PDstandardNth323(u,i,j,k) (p1o144dydz2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)])) +#define PDstandardNth211(u,i,j,k) (p1o144dx2dy*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)] + 30*(8*(u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 8*(u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)]))) +#define PDstandardNth311(u,i,j,k) (p1o144dx2dz*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)] + 30*(8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))) +#define PDstandardNth122(u,i,j,k) (p1o144dxdy2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,k)]) - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,k)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) +#define PDstandardNth322(u,i,j,k) (p1o144dy2dz*(-128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)] + 30*(8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,-1 + k)] - 8*(u)[CCTK_GFINDEX3D(cctkGH,i,j,1 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,j,2 + k)]))) +#define PDstandardNth133(u,i,j,k) (p1o144dxdz2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,2 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,2 + k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,-1 + i,j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,1 + i,j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,1 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,1 + k)]) - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,2 + k)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,-2 + i,j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,2 + i,j,k)])) +#define PDstandardNth233(u,i,j,k) (p1o144dydz2*(-128*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,1 + k)]) + 128*((u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,2 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,2 + k)]) + 240*(u)[CCTK_GFINDEX3D(cctkGH,i,-1 + j,k)] - 240*(u)[CCTK_GFINDEX3D(cctkGH,i,1 + j,k)] + 16*((u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,1 + k)]) - 16*((u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,1 + k)]) - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,2 + k)] - 30*(u)[CCTK_GFINDEX3D(cctkGH,i,-2 + j,k)] + 30*(u)[CCTK_GFINDEX3D(cctkGH,i,2 + j,k)])) +#define PDstandardNth123(u,i,j,k) (p1o1728dxdydz*(-512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)])) +#define PDstandardNth132(u,i,j,k) (p1o1728dxdydz*(-512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)])) +#define PDstandardNth213(u,i,j,k) (p1o1728dxdydz*(-512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)])) +#define PDstandardNth231(u,i,j,k) (p1o1728dxdydz*(-512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)])) +#define PDstandardNth312(u,i,j,k) (p1o1728dxdydz*(-512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)])) +#define PDstandardNth321(u,i,j,k) (p1o1728dxdydz*(-512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-1 + k)]) + 512*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,1 + k)]) + 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-1 + k)]) - 64*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,1 + k)]) - 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-1 + k)]) + 8*((u)[CCTK_GFINDEX3D(cctkGH,-1 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-1 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,1 + i,2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-1 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,1 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-1 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,1 + k)]) + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,-2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,-2 + i,2 + j,2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,-2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,-2 + j,2 + k)] + (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,-2 + k)] - (u)[CCTK_GFINDEX3D(cctkGH,2 + i,2 + j,2 + k)])) diff --git a/ML_WaveToy/src/WT_Gaussian.c b/ML_WaveToy/src/WT_Gaussian.c index bbba1a1..5caf9ca 100644 --- a/ML_WaveToy/src/WT_Gaussian.c +++ b/ML_WaveToy/src/WT_Gaussian.c @@ -104,6 +104,7 @@ void WT_Gaussian_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal pm1o12dz2 = -pow(dz,-2)/12.; /* Loop over the grid points */ + _Pragma ("omp parallel") LC_LOOP3 (WT_Gaussian, i,j,k, min[0],min[1],min[2], max[0],max[1],max[2], cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) diff --git a/ML_WaveToy/src/WT_RHS.c b/ML_WaveToy/src/WT_RHS.c index f06f42a..41d3939 100644 --- a/ML_WaveToy/src/WT_RHS.c +++ b/ML_WaveToy/src/WT_RHS.c @@ -104,6 +104,7 @@ void WT_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[3], pm1o12dz2 = -pow(dz,-2)/12.; /* Loop over the grid points */ + _Pragma ("omp parallel") LC_LOOP3 (WT_RHS, i,j,k, min[0],min[1],min[2], max[0],max[1],max[2], cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) |