diff options
Diffstat (limited to 'ML_BSSN_MP_O8/src/Differencing.h')
| -rw-r--r-- | ML_BSSN_MP_O8/src/Differencing.h | 54 |
1 files changed, 27 insertions, 27 deletions
diff --git a/ML_BSSN_MP_O8/src/Differencing.h b/ML_BSSN_MP_O8/src/Differencing.h index 450d07a..5cd43c9 100644 --- a/ML_BSSN_MP_O8/src/Differencing.h +++ b/ML_BSSN_MP_O8/src/Differencing.h @@ -1,27 +1,27 @@ -#define PDstandardNth1(u,i,j,k) (p1o840dx*(-672*(u)[index+di*(-1)+dj*(0)+dk*(0)] + 672*(u)[index+di*(1)+dj*(0)+dk*(0)] + 168*(u)[index+di*(-2)+dj*(0)+dk*(0)] - 168*(u)[index+di*(2)+dj*(0)+dk*(0)] - 32*(u)[index+di*(-3)+dj*(0)+dk*(0)] + 32*(u)[index+di*(3)+dj*(0)+dk*(0)] + 3*(u)[index+di*(-4)+dj*(0)+dk*(0)] - 3*(u)[index+di*(4)+dj*(0)+dk*(0)])) -#define PDstandardNth2(u,i,j,k) (p1o840dy*(-672*(u)[index+di*(0)+dj*(-1)+dk*(0)] + 672*(u)[index+di*(0)+dj*(1)+dk*(0)] + 168*(u)[index+di*(0)+dj*(-2)+dk*(0)] - 168*(u)[index+di*(0)+dj*(2)+dk*(0)] - 32*(u)[index+di*(0)+dj*(-3)+dk*(0)] + 32*(u)[index+di*(0)+dj*(3)+dk*(0)] + 3*(u)[index+di*(0)+dj*(-4)+dk*(0)] - 3*(u)[index+di*(0)+dj*(4)+dk*(0)])) -#define PDstandardNth3(u,i,j,k) (p1o840dz*(-672*(u)[index+di*(0)+dj*(0)+dk*(-1)] + 672*(u)[index+di*(0)+dj*(0)+dk*(1)] + 168*(u)[index+di*(0)+dj*(0)+dk*(-2)] - 168*(u)[index+di*(0)+dj*(0)+dk*(2)] - 32*(u)[index+di*(0)+dj*(0)+dk*(-3)] + 32*(u)[index+di*(0)+dj*(0)+dk*(3)] + 3*(u)[index+di*(0)+dj*(0)+dk*(-4)] - 3*(u)[index+di*(0)+dj*(0)+dk*(4)])) -#define PDstandardNth11(u,i,j,k) (p1o5040dx2*(-14350*(u)[index+di*(0)+dj*(0)+dk*(0)] + 8064*((u)[index+di*(-1)+dj*(0)+dk*(0)] + (u)[index+di*(1)+dj*(0)+dk*(0)]) - 1008*((u)[index+di*(-2)+dj*(0)+dk*(0)] + (u)[index+di*(2)+dj*(0)+dk*(0)]) + 128*((u)[index+di*(-3)+dj*(0)+dk*(0)] + (u)[index+di*(3)+dj*(0)+dk*(0)]) - 9*((u)[index+di*(-4)+dj*(0)+dk*(0)] + (u)[index+di*(4)+dj*(0)+dk*(0)]))) -#define PDstandardNth22(u,i,j,k) (p1o5040dy2*(-14350*(u)[index+di*(0)+dj*(0)+dk*(0)] + 8064*((u)[index+di*(0)+dj*(-1)+dk*(0)] + (u)[index+di*(0)+dj*(1)+dk*(0)]) - 1008*((u)[index+di*(0)+dj*(-2)+dk*(0)] + (u)[index+di*(0)+dj*(2)+dk*(0)]) + 128*((u)[index+di*(0)+dj*(-3)+dk*(0)] + (u)[index+di*(0)+dj*(3)+dk*(0)]) - 9*((u)[index+di*(0)+dj*(-4)+dk*(0)] + (u)[index+di*(0)+dj*(4)+dk*(0)]))) -#define PDstandardNth33(u,i,j,k) (p1o5040dz2*(-14350*(u)[index+di*(0)+dj*(0)+dk*(0)] + 8064*((u)[index+di*(0)+dj*(0)+dk*(-1)] + (u)[index+di*(0)+dj*(0)+dk*(1)]) - 1008*((u)[index+di*(0)+dj*(0)+dk*(-2)] + (u)[index+di*(0)+dj*(0)+dk*(2)]) + 128*((u)[index+di*(0)+dj*(0)+dk*(-3)] + (u)[index+di*(0)+dj*(0)+dk*(3)]) - 9*((u)[index+di*(0)+dj*(0)+dk*(-4)] + (u)[index+di*(0)+dj*(0)+dk*(4)]))) -#define PDstandardNth12(u,i,j,k) (p1o705600dxdy*(-451584*((u)[index+di*(-1)+dj*(1)+dk*(0)] + (u)[index+di*(1)+dj*(-1)+dk*(0)]) + 451584*((u)[index+di*(-1)+dj*(-1)+dk*(0)] + (u)[index+di*(1)+dj*(1)+dk*(0)]) + 112896*((u)[index+di*(-1)+dj*(2)+dk*(0)] + (u)[index+di*(1)+dj*(-2)+dk*(0)] + (u)[index+di*(-2)+dj*(1)+dk*(0)] + (u)[index+di*(2)+dj*(-1)+dk*(0)]) - 112896*((u)[index+di*(-1)+dj*(-2)+dk*(0)] + (u)[index+di*(1)+dj*(2)+dk*(0)] + (u)[index+di*(-2)+dj*(-1)+dk*(0)] + (u)[index+di*(2)+dj*(1)+dk*(0)]) - 28224*((u)[index+di*(-2)+dj*(2)+dk*(0)] + (u)[index+di*(2)+dj*(-2)+dk*(0)]) + 28224*((u)[index+di*(-2)+dj*(-2)+dk*(0)] + (u)[index+di*(2)+dj*(2)+dk*(0)]) - 21504*((u)[index+di*(-1)+dj*(3)+dk*(0)] + (u)[index+di*(1)+dj*(-3)+dk*(0)] + (u)[index+di*(-3)+dj*(1)+dk*(0)] + (u)[index+di*(3)+dj*(-1)+dk*(0)]) + 21504*((u)[index+di*(-1)+dj*(-3)+dk*(0)] + (u)[index+di*(1)+dj*(3)+dk*(0)] + (u)[index+di*(-3)+dj*(-1)+dk*(0)] + (u)[index+di*(3)+dj*(1)+dk*(0)]) + 5376*((u)[index+di*(-2)+dj*(3)+dk*(0)] + (u)[index+di*(2)+dj*(-3)+dk*(0)] + (u)[index+di*(-3)+dj*(2)+dk*(0)] + (u)[index+di*(3)+dj*(-2)+dk*(0)]) - 5376*((u)[index+di*(-2)+dj*(-3)+dk*(0)] + (u)[index+di*(2)+dj*(3)+dk*(0)] + (u)[index+di*(-3)+dj*(-2)+dk*(0)] + (u)[index+di*(3)+dj*(2)+dk*(0)]) - 1024*((u)[index+di*(-3)+dj*(3)+dk*(0)] + (u)[index+di*(3)+dj*(-3)+dk*(0)]) + 1024*((u)[index+di*(-3)+dj*(-3)+dk*(0)] + (u)[index+di*(3)+dj*(3)+dk*(0)]) + 2016*((u)[index+di*(-1)+dj*(4)+dk*(0)] + (u)[index+di*(1)+dj*(-4)+dk*(0)] + (u)[index+di*(-4)+dj*(1)+dk*(0)] + (u)[index+di*(4)+dj*(-1)+dk*(0)]) - 2016*((u)[index+di*(-1)+dj*(-4)+dk*(0)] + (u)[index+di*(1)+dj*(4)+dk*(0)] + (u)[index+di*(-4)+dj*(-1)+dk*(0)] + (u)[index+di*(4)+dj*(1)+dk*(0)]) - 504*((u)[index+di*(-2)+dj*(4)+dk*(0)] + (u)[index+di*(2)+dj*(-4)+dk*(0)] + (u)[index+di*(-4)+dj*(2)+dk*(0)] + (u)[index+di*(4)+dj*(-2)+dk*(0)]) + 504*((u)[index+di*(-2)+dj*(-4)+dk*(0)] + (u)[index+di*(2)+dj*(4)+dk*(0)] + (u)[index+di*(-4)+dj*(-2)+dk*(0)] + (u)[index+di*(4)+dj*(2)+dk*(0)]) + 96*((u)[index+di*(-3)+dj*(4)+dk*(0)] + (u)[index+di*(3)+dj*(-4)+dk*(0)] + (u)[index+di*(-4)+dj*(3)+dk*(0)] + (u)[index+di*(4)+dj*(-3)+dk*(0)]) - 96*((u)[index+di*(-3)+dj*(-4)+dk*(0)] + (u)[index+di*(3)+dj*(4)+dk*(0)] + (u)[index+di*(-4)+dj*(-3)+dk*(0)] + (u)[index+di*(4)+dj*(3)+dk*(0)]) - 9*((u)[index+di*(-4)+dj*(4)+dk*(0)] + (u)[index+di*(4)+dj*(-4)+dk*(0)]) + 9*((u)[index+di*(-4)+dj*(-4)+dk*(0)] + (u)[index+di*(4)+dj*(4)+dk*(0)]))) -#define PDstandardNth13(u,i,j,k) (p1o705600dxdz*(-451584*((u)[index+di*(-1)+dj*(0)+dk*(1)] + (u)[index+di*(1)+dj*(0)+dk*(-1)]) + 451584*((u)[index+di*(-1)+dj*(0)+dk*(-1)] + (u)[index+di*(1)+dj*(0)+dk*(1)]) + 112896*((u)[index+di*(-1)+dj*(0)+dk*(2)] + (u)[index+di*(1)+dj*(0)+dk*(-2)] + (u)[index+di*(-2)+dj*(0)+dk*(1)] + (u)[index+di*(2)+dj*(0)+dk*(-1)]) - 112896*((u)[index+di*(-1)+dj*(0)+dk*(-2)] + (u)[index+di*(1)+dj*(0)+dk*(2)] + (u)[index+di*(-2)+dj*(0)+dk*(-1)] + (u)[index+di*(2)+dj*(0)+dk*(1)]) - 28224*((u)[index+di*(-2)+dj*(0)+dk*(2)] + (u)[index+di*(2)+dj*(0)+dk*(-2)]) + 28224*((u)[index+di*(-2)+dj*(0)+dk*(-2)] + (u)[index+di*(2)+dj*(0)+dk*(2)]) - 21504*((u)[index+di*(-1)+dj*(0)+dk*(3)] + (u)[index+di*(1)+dj*(0)+dk*(-3)] + (u)[index+di*(-3)+dj*(0)+dk*(1)] + (u)[index+di*(3)+dj*(0)+dk*(-1)]) + 21504*((u)[index+di*(-1)+dj*(0)+dk*(-3)] + (u)[index+di*(1)+dj*(0)+dk*(3)] + (u)[index+di*(-3)+dj*(0)+dk*(-1)] + (u)[index+di*(3)+dj*(0)+dk*(1)]) + 5376*((u)[index+di*(-2)+dj*(0)+dk*(3)] + (u)[index+di*(2)+dj*(0)+dk*(-3)] + (u)[index+di*(-3)+dj*(0)+dk*(2)] + (u)[index+di*(3)+dj*(0)+dk*(-2)]) - 5376*((u)[index+di*(-2)+dj*(0)+dk*(-3)] + (u)[index+di*(2)+dj*(0)+dk*(3)] + (u)[index+di*(-3)+dj*(0)+dk*(-2)] + (u)[index+di*(3)+dj*(0)+dk*(2)]) - 1024*((u)[index+di*(-3)+dj*(0)+dk*(3)] + (u)[index+di*(3)+dj*(0)+dk*(-3)]) + 1024*((u)[index+di*(-3)+dj*(0)+dk*(-3)] + (u)[index+di*(3)+dj*(0)+dk*(3)]) + 2016*((u)[index+di*(-1)+dj*(0)+dk*(4)] + (u)[index+di*(1)+dj*(0)+dk*(-4)] + (u)[index+di*(-4)+dj*(0)+dk*(1)] + (u)[index+di*(4)+dj*(0)+dk*(-1)]) - 2016*((u)[index+di*(-1)+dj*(0)+dk*(-4)] + (u)[index+di*(1)+dj*(0)+dk*(4)] + (u)[index+di*(-4)+dj*(0)+dk*(-1)] + (u)[index+di*(4)+dj*(0)+dk*(1)]) - 504*((u)[index+di*(-2)+dj*(0)+dk*(4)] + (u)[index+di*(2)+dj*(0)+dk*(-4)] + (u)[index+di*(-4)+dj*(0)+dk*(2)] + (u)[index+di*(4)+dj*(0)+dk*(-2)]) + 504*((u)[index+di*(-2)+dj*(0)+dk*(-4)] + (u)[index+di*(2)+dj*(0)+dk*(4)] + (u)[index+di*(-4)+dj*(0)+dk*(-2)] + (u)[index+di*(4)+dj*(0)+dk*(2)]) + 96*((u)[index+di*(-3)+dj*(0)+dk*(4)] + (u)[index+di*(3)+dj*(0)+dk*(-4)] + (u)[index+di*(-4)+dj*(0)+dk*(3)] + (u)[index+di*(4)+dj*(0)+dk*(-3)]) - 96*((u)[index+di*(-3)+dj*(0)+dk*(-4)] + (u)[index+di*(3)+dj*(0)+dk*(4)] + (u)[index+di*(-4)+dj*(0)+dk*(-3)] + (u)[index+di*(4)+dj*(0)+dk*(3)]) - 9*((u)[index+di*(-4)+dj*(0)+dk*(4)] + (u)[index+di*(4)+dj*(0)+dk*(-4)]) + 9*((u)[index+di*(-4)+dj*(0)+dk*(-4)] + (u)[index+di*(4)+dj*(0)+dk*(4)]))) -#define PDstandardNth21(u,i,j,k) (p1o705600dxdy*(-451584*((u)[index+di*(-1)+dj*(1)+dk*(0)] + (u)[index+di*(1)+dj*(-1)+dk*(0)]) + 451584*((u)[index+di*(-1)+dj*(-1)+dk*(0)] + (u)[index+di*(1)+dj*(1)+dk*(0)]) + 112896*((u)[index+di*(-1)+dj*(2)+dk*(0)] + (u)[index+di*(1)+dj*(-2)+dk*(0)] + (u)[index+di*(-2)+dj*(1)+dk*(0)] + (u)[index+di*(2)+dj*(-1)+dk*(0)]) - 112896*((u)[index+di*(-1)+dj*(-2)+dk*(0)] + (u)[index+di*(1)+dj*(2)+dk*(0)] + (u)[index+di*(-2)+dj*(-1)+dk*(0)] + (u)[index+di*(2)+dj*(1)+dk*(0)]) - 28224*((u)[index+di*(-2)+dj*(2)+dk*(0)] + (u)[index+di*(2)+dj*(-2)+dk*(0)]) + 28224*((u)[index+di*(-2)+dj*(-2)+dk*(0)] + (u)[index+di*(2)+dj*(2)+dk*(0)]) - 21504*((u)[index+di*(-1)+dj*(3)+dk*(0)] + (u)[index+di*(1)+dj*(-3)+dk*(0)] + (u)[index+di*(-3)+dj*(1)+dk*(0)] + (u)[index+di*(3)+dj*(-1)+dk*(0)]) + 21504*((u)[index+di*(-1)+dj*(-3)+dk*(0)] + (u)[index+di*(1)+dj*(3)+dk*(0)] + (u)[index+di*(-3)+dj*(-1)+dk*(0)] + (u)[index+di*(3)+dj*(1)+dk*(0)]) + 5376*((u)[index+di*(-2)+dj*(3)+dk*(0)] + (u)[index+di*(2)+dj*(-3)+dk*(0)] + (u)[index+di*(-3)+dj*(2)+dk*(0)] + (u)[index+di*(3)+dj*(-2)+dk*(0)]) - 5376*((u)[index+di*(-2)+dj*(-3)+dk*(0)] + (u)[index+di*(2)+dj*(3)+dk*(0)] + (u)[index+di*(-3)+dj*(-2)+dk*(0)] + (u)[index+di*(3)+dj*(2)+dk*(0)]) - 1024*((u)[index+di*(-3)+dj*(3)+dk*(0)] + (u)[index+di*(3)+dj*(-3)+dk*(0)]) + 1024*((u)[index+di*(-3)+dj*(-3)+dk*(0)] + (u)[index+di*(3)+dj*(3)+dk*(0)]) + 2016*((u)[index+di*(-1)+dj*(4)+dk*(0)] + (u)[index+di*(1)+dj*(-4)+dk*(0)] + (u)[index+di*(-4)+dj*(1)+dk*(0)] + (u)[index+di*(4)+dj*(-1)+dk*(0)]) - 2016*((u)[index+di*(-1)+dj*(-4)+dk*(0)] + (u)[index+di*(1)+dj*(4)+dk*(0)] + (u)[index+di*(-4)+dj*(-1)+dk*(0)] + (u)[index+di*(4)+dj*(1)+dk*(0)]) - 504*((u)[index+di*(-2)+dj*(4)+dk*(0)] + (u)[index+di*(2)+dj*(-4)+dk*(0)] + (u)[index+di*(-4)+dj*(2)+dk*(0)] + (u)[index+di*(4)+dj*(-2)+dk*(0)]) + 504*((u)[index+di*(-2)+dj*(-4)+dk*(0)] + (u)[index+di*(2)+dj*(4)+dk*(0)] + (u)[index+di*(-4)+dj*(-2)+dk*(0)] + (u)[index+di*(4)+dj*(2)+dk*(0)]) + 96*((u)[index+di*(-3)+dj*(4)+dk*(0)] + (u)[index+di*(3)+dj*(-4)+dk*(0)] + (u)[index+di*(-4)+dj*(3)+dk*(0)] + (u)[index+di*(4)+dj*(-3)+dk*(0)]) - 96*((u)[index+di*(-3)+dj*(-4)+dk*(0)] + (u)[index+di*(3)+dj*(4)+dk*(0)] + (u)[index+di*(-4)+dj*(-3)+dk*(0)] + (u)[index+di*(4)+dj*(3)+dk*(0)]) - 9*((u)[index+di*(-4)+dj*(4)+dk*(0)] + (u)[index+di*(4)+dj*(-4)+dk*(0)]) + 9*((u)[index+di*(-4)+dj*(-4)+dk*(0)] + (u)[index+di*(4)+dj*(4)+dk*(0)]))) -#define PDstandardNth23(u,i,j,k) (p1o705600dydz*(-451584*((u)[index+di*(0)+dj*(-1)+dk*(1)] + (u)[index+di*(0)+dj*(1)+dk*(-1)]) + 451584*((u)[index+di*(0)+dj*(-1)+dk*(-1)] + (u)[index+di*(0)+dj*(1)+dk*(1)]) + 112896*((u)[index+di*(0)+dj*(-1)+dk*(2)] + (u)[index+di*(0)+dj*(1)+dk*(-2)] + (u)[index+di*(0)+dj*(-2)+dk*(1)] + (u)[index+di*(0)+dj*(2)+dk*(-1)]) - 112896*((u)[index+di*(0)+dj*(-1)+dk*(-2)] + (u)[index+di*(0)+dj*(1)+dk*(2)] + (u)[index+di*(0)+dj*(-2)+dk*(-1)] + (u)[index+di*(0)+dj*(2)+dk*(1)]) - 28224*((u)[index+di*(0)+dj*(-2)+dk*(2)] + (u)[index+di*(0)+dj*(2)+dk*(-2)]) + 28224*((u)[index+di*(0)+dj*(-2)+dk*(-2)] + (u)[index+di*(0)+dj*(2)+dk*(2)]) - 21504*((u)[index+di*(0)+dj*(-1)+dk*(3)] + (u)[index+di*(0)+dj*(1)+dk*(-3)] + (u)[index+di*(0)+dj*(-3)+dk*(1)] + (u)[index+di*(0)+dj*(3)+dk*(-1)]) + 21504*((u)[index+di*(0)+dj*(-1)+dk*(-3)] + (u)[index+di*(0)+dj*(1)+dk*(3)] + (u)[index+di*(0)+dj*(-3)+dk*(-1)] + (u)[index+di*(0)+dj*(3)+dk*(1)]) + 5376*((u)[index+di*(0)+dj*(-2)+dk*(3)] + (u)[index+di*(0)+dj*(2)+dk*(-3)] + (u)[index+di*(0)+dj*(-3)+dk*(2)] + (u)[index+di*(0)+dj*(3)+dk*(-2)]) - 5376*((u)[index+di*(0)+dj*(-2)+dk*(-3)] + (u)[index+di*(0)+dj*(2)+dk*(3)] + (u)[index+di*(0)+dj*(-3)+dk*(-2)] + (u)[index+di*(0)+dj*(3)+dk*(2)]) - 1024*((u)[index+di*(0)+dj*(-3)+dk*(3)] + (u)[index+di*(0)+dj*(3)+dk*(-3)]) + 1024*((u)[index+di*(0)+dj*(-3)+dk*(-3)] + (u)[index+di*(0)+dj*(3)+dk*(3)]) + 2016*((u)[index+di*(0)+dj*(-1)+dk*(4)] + (u)[index+di*(0)+dj*(1)+dk*(-4)] + (u)[index+di*(0)+dj*(-4)+dk*(1)] + (u)[index+di*(0)+dj*(4)+dk*(-1)]) - 2016*((u)[index+di*(0)+dj*(-1)+dk*(-4)] + (u)[index+di*(0)+dj*(1)+dk*(4)] + (u)[index+di*(0)+dj*(-4)+dk*(-1)] + (u)[index+di*(0)+dj*(4)+dk*(1)]) - 504*((u)[index+di*(0)+dj*(-2)+dk*(4)] + (u)[index+di*(0)+dj*(2)+dk*(-4)] + (u)[index+di*(0)+dj*(-4)+dk*(2)] + (u)[index+di*(0)+dj*(4)+dk*(-2)]) + 504*((u)[index+di*(0)+dj*(-2)+dk*(-4)] + (u)[index+di*(0)+dj*(2)+dk*(4)] + (u)[index+di*(0)+dj*(-4)+dk*(-2)] + (u)[index+di*(0)+dj*(4)+dk*(2)]) + 96*((u)[index+di*(0)+dj*(-3)+dk*(4)] + (u)[index+di*(0)+dj*(3)+dk*(-4)] + (u)[index+di*(0)+dj*(-4)+dk*(3)] + (u)[index+di*(0)+dj*(4)+dk*(-3)]) - 96*((u)[index+di*(0)+dj*(-3)+dk*(-4)] + (u)[index+di*(0)+dj*(3)+dk*(4)] + (u)[index+di*(0)+dj*(-4)+dk*(-3)] + (u)[index+di*(0)+dj*(4)+dk*(3)]) - 9*((u)[index+di*(0)+dj*(-4)+dk*(4)] + (u)[index+di*(0)+dj*(4)+dk*(-4)]) + 9*((u)[index+di*(0)+dj*(-4)+dk*(-4)] + (u)[index+di*(0)+dj*(4)+dk*(4)]))) -#define PDstandardNth31(u,i,j,k) (p1o705600dxdz*(-451584*((u)[index+di*(-1)+dj*(0)+dk*(1)] + (u)[index+di*(1)+dj*(0)+dk*(-1)]) + 451584*((u)[index+di*(-1)+dj*(0)+dk*(-1)] + (u)[index+di*(1)+dj*(0)+dk*(1)]) + 112896*((u)[index+di*(-1)+dj*(0)+dk*(2)] + (u)[index+di*(1)+dj*(0)+dk*(-2)] + (u)[index+di*(-2)+dj*(0)+dk*(1)] + (u)[index+di*(2)+dj*(0)+dk*(-1)]) - 112896*((u)[index+di*(-1)+dj*(0)+dk*(-2)] + (u)[index+di*(1)+dj*(0)+dk*(2)] + (u)[index+di*(-2)+dj*(0)+dk*(-1)] + (u)[index+di*(2)+dj*(0)+dk*(1)]) - 28224*((u)[index+di*(-2)+dj*(0)+dk*(2)] + (u)[index+di*(2)+dj*(0)+dk*(-2)]) + 28224*((u)[index+di*(-2)+dj*(0)+dk*(-2)] + (u)[index+di*(2)+dj*(0)+dk*(2)]) - 21504*((u)[index+di*(-1)+dj*(0)+dk*(3)] + (u)[index+di*(1)+dj*(0)+dk*(-3)] + (u)[index+di*(-3)+dj*(0)+dk*(1)] + (u)[index+di*(3)+dj*(0)+dk*(-1)]) + 21504*((u)[index+di*(-1)+dj*(0)+dk*(-3)] + (u)[index+di*(1)+dj*(0)+dk*(3)] + (u)[index+di*(-3)+dj*(0)+dk*(-1)] + (u)[index+di*(3)+dj*(0)+dk*(1)]) + 5376*((u)[index+di*(-2)+dj*(0)+dk*(3)] + (u)[index+di*(2)+dj*(0)+dk*(-3)] + (u)[index+di*(-3)+dj*(0)+dk*(2)] + (u)[index+di*(3)+dj*(0)+dk*(-2)]) - 5376*((u)[index+di*(-2)+dj*(0)+dk*(-3)] + (u)[index+di*(2)+dj*(0)+dk*(3)] + (u)[index+di*(-3)+dj*(0)+dk*(-2)] + (u)[index+di*(3)+dj*(0)+dk*(2)]) - 1024*((u)[index+di*(-3)+dj*(0)+dk*(3)] + (u)[index+di*(3)+dj*(0)+dk*(-3)]) + 1024*((u)[index+di*(-3)+dj*(0)+dk*(-3)] + (u)[index+di*(3)+dj*(0)+dk*(3)]) + 2016*((u)[index+di*(-1)+dj*(0)+dk*(4)] + (u)[index+di*(1)+dj*(0)+dk*(-4)] + (u)[index+di*(-4)+dj*(0)+dk*(1)] + (u)[index+di*(4)+dj*(0)+dk*(-1)]) - 2016*((u)[index+di*(-1)+dj*(0)+dk*(-4)] + (u)[index+di*(1)+dj*(0)+dk*(4)] + (u)[index+di*(-4)+dj*(0)+dk*(-1)] + (u)[index+di*(4)+dj*(0)+dk*(1)]) - 504*((u)[index+di*(-2)+dj*(0)+dk*(4)] + (u)[index+di*(2)+dj*(0)+dk*(-4)] + (u)[index+di*(-4)+dj*(0)+dk*(2)] + (u)[index+di*(4)+dj*(0)+dk*(-2)]) + 504*((u)[index+di*(-2)+dj*(0)+dk*(-4)] + (u)[index+di*(2)+dj*(0)+dk*(4)] + (u)[index+di*(-4)+dj*(0)+dk*(-2)] + (u)[index+di*(4)+dj*(0)+dk*(2)]) + 96*((u)[index+di*(-3)+dj*(0)+dk*(4)] + (u)[index+di*(3)+dj*(0)+dk*(-4)] + (u)[index+di*(-4)+dj*(0)+dk*(3)] + (u)[index+di*(4)+dj*(0)+dk*(-3)]) - 96*((u)[index+di*(-3)+dj*(0)+dk*(-4)] + (u)[index+di*(3)+dj*(0)+dk*(4)] + (u)[index+di*(-4)+dj*(0)+dk*(-3)] + (u)[index+di*(4)+dj*(0)+dk*(3)]) - 9*((u)[index+di*(-4)+dj*(0)+dk*(4)] + (u)[index+di*(4)+dj*(0)+dk*(-4)]) + 9*((u)[index+di*(-4)+dj*(0)+dk*(-4)] + (u)[index+di*(4)+dj*(0)+dk*(4)]))) -#define PDstandardNth32(u,i,j,k) (p1o705600dydz*(-451584*((u)[index+di*(0)+dj*(-1)+dk*(1)] + (u)[index+di*(0)+dj*(1)+dk*(-1)]) + 451584*((u)[index+di*(0)+dj*(-1)+dk*(-1)] + (u)[index+di*(0)+dj*(1)+dk*(1)]) + 112896*((u)[index+di*(0)+dj*(-1)+dk*(2)] + (u)[index+di*(0)+dj*(1)+dk*(-2)] + (u)[index+di*(0)+dj*(-2)+dk*(1)] + (u)[index+di*(0)+dj*(2)+dk*(-1)]) - 112896*((u)[index+di*(0)+dj*(-1)+dk*(-2)] + (u)[index+di*(0)+dj*(1)+dk*(2)] + (u)[index+di*(0)+dj*(-2)+dk*(-1)] + (u)[index+di*(0)+dj*(2)+dk*(1)]) - 28224*((u)[index+di*(0)+dj*(-2)+dk*(2)] + (u)[index+di*(0)+dj*(2)+dk*(-2)]) + 28224*((u)[index+di*(0)+dj*(-2)+dk*(-2)] + (u)[index+di*(0)+dj*(2)+dk*(2)]) - 21504*((u)[index+di*(0)+dj*(-1)+dk*(3)] + (u)[index+di*(0)+dj*(1)+dk*(-3)] + (u)[index+di*(0)+dj*(-3)+dk*(1)] + (u)[index+di*(0)+dj*(3)+dk*(-1)]) + 21504*((u)[index+di*(0)+dj*(-1)+dk*(-3)] + (u)[index+di*(0)+dj*(1)+dk*(3)] + (u)[index+di*(0)+dj*(-3)+dk*(-1)] + (u)[index+di*(0)+dj*(3)+dk*(1)]) + 5376*((u)[index+di*(0)+dj*(-2)+dk*(3)] + (u)[index+di*(0)+dj*(2)+dk*(-3)] + (u)[index+di*(0)+dj*(-3)+dk*(2)] + (u)[index+di*(0)+dj*(3)+dk*(-2)]) - 5376*((u)[index+di*(0)+dj*(-2)+dk*(-3)] + (u)[index+di*(0)+dj*(2)+dk*(3)] + (u)[index+di*(0)+dj*(-3)+dk*(-2)] + (u)[index+di*(0)+dj*(3)+dk*(2)]) - 1024*((u)[index+di*(0)+dj*(-3)+dk*(3)] + (u)[index+di*(0)+dj*(3)+dk*(-3)]) + 1024*((u)[index+di*(0)+dj*(-3)+dk*(-3)] + (u)[index+di*(0)+dj*(3)+dk*(3)]) + 2016*((u)[index+di*(0)+dj*(-1)+dk*(4)] + (u)[index+di*(0)+dj*(1)+dk*(-4)] + (u)[index+di*(0)+dj*(-4)+dk*(1)] + (u)[index+di*(0)+dj*(4)+dk*(-1)]) - 2016*((u)[index+di*(0)+dj*(-1)+dk*(-4)] + (u)[index+di*(0)+dj*(1)+dk*(4)] + (u)[index+di*(0)+dj*(-4)+dk*(-1)] + (u)[index+di*(0)+dj*(4)+dk*(1)]) - 504*((u)[index+di*(0)+dj*(-2)+dk*(4)] + (u)[index+di*(0)+dj*(2)+dk*(-4)] + (u)[index+di*(0)+dj*(-4)+dk*(2)] + (u)[index+di*(0)+dj*(4)+dk*(-2)]) + 504*((u)[index+di*(0)+dj*(-2)+dk*(-4)] + (u)[index+di*(0)+dj*(2)+dk*(4)] + (u)[index+di*(0)+dj*(-4)+dk*(-2)] + (u)[index+di*(0)+dj*(4)+dk*(2)]) + 96*((u)[index+di*(0)+dj*(-3)+dk*(4)] + (u)[index+di*(0)+dj*(3)+dk*(-4)] + (u)[index+di*(0)+dj*(-4)+dk*(3)] + (u)[index+di*(0)+dj*(4)+dk*(-3)]) - 96*((u)[index+di*(0)+dj*(-3)+dk*(-4)] + (u)[index+di*(0)+dj*(3)+dk*(4)] + (u)[index+di*(0)+dj*(-4)+dk*(-3)] + (u)[index+di*(0)+dj*(4)+dk*(3)]) - 9*((u)[index+di*(0)+dj*(-4)+dk*(4)] + (u)[index+di*(0)+dj*(4)+dk*(-4)]) + 9*((u)[index+di*(0)+dj*(-4)+dk*(-4)] + (u)[index+di*(0)+dj*(4)+dk*(4)]))) -#define PDdissipationNth1(u,i,j,k) (p1o1024dx*(-252*(u)[index+di*(0)+dj*(0)+dk*(0)] + 210*((u)[index+di*(-1)+dj*(0)+dk*(0)] + (u)[index+di*(1)+dj*(0)+dk*(0)]) - 120*((u)[index+di*(-2)+dj*(0)+dk*(0)] + (u)[index+di*(2)+dj*(0)+dk*(0)]) + 45*((u)[index+di*(-3)+dj*(0)+dk*(0)] + (u)[index+di*(3)+dj*(0)+dk*(0)]) - 10*((u)[index+di*(-4)+dj*(0)+dk*(0)] + (u)[index+di*(4)+dj*(0)+dk*(0)]) + (u)[index+di*(-5)+dj*(0)+dk*(0)] + (u)[index+di*(5)+dj*(0)+dk*(0)])) -#define PDdissipationNth2(u,i,j,k) (p1o1024dy*(-252*(u)[index+di*(0)+dj*(0)+dk*(0)] + 210*((u)[index+di*(0)+dj*(-1)+dk*(0)] + (u)[index+di*(0)+dj*(1)+dk*(0)]) - 120*((u)[index+di*(0)+dj*(-2)+dk*(0)] + (u)[index+di*(0)+dj*(2)+dk*(0)]) + 45*((u)[index+di*(0)+dj*(-3)+dk*(0)] + (u)[index+di*(0)+dj*(3)+dk*(0)]) - 10*((u)[index+di*(0)+dj*(-4)+dk*(0)] + (u)[index+di*(0)+dj*(4)+dk*(0)]) + (u)[index+di*(0)+dj*(-5)+dk*(0)] + (u)[index+di*(0)+dj*(5)+dk*(0)])) -#define PDdissipationNth3(u,i,j,k) (p1o1024dz*(-252*(u)[index+di*(0)+dj*(0)+dk*(0)] + 210*((u)[index+di*(0)+dj*(0)+dk*(-1)] + (u)[index+di*(0)+dj*(0)+dk*(1)]) - 120*((u)[index+di*(0)+dj*(0)+dk*(-2)] + (u)[index+di*(0)+dj*(0)+dk*(2)]) + 45*((u)[index+di*(0)+dj*(0)+dk*(-3)] + (u)[index+di*(0)+dj*(0)+dk*(3)]) - 10*((u)[index+di*(0)+dj*(0)+dk*(-4)] + (u)[index+di*(0)+dj*(0)+dk*(4)]) + (u)[index+di*(0)+dj*(0)+dk*(-5)] + (u)[index+di*(0)+dj*(0)+dk*(5)])) -#define PDupwindNth1(u,i,j,k) (pm1o840dx*(378*(u)[index+di*(0)+dj*(0)+dk*(0)] - 60*(u)[index+di*(-2*dir1)+dj*(0)+dk*(0)] + 5*(u)[index+di*(-3*dir1)+dj*(0)+dk*(0)] - 140*(u)[index+di*(3*dir1)+dj*(0)+dk*(0)] + 30*(u)[index+di*(4*dir1)+dj*(0)+dk*(0)] - 3*(u)[index+di*(5*dir1)+dj*(0)+dk*(0)] + 420*((u)[index+di*(2*dir1)+dj*(0)+dk*(0)] + (u)[index+di*(-dir1)+dj*(0)+dk*(0)]) - 1050*(u)[index+di*(dir1)+dj*(0)+dk*(0)])*dir1) -#define PDupwindNth2(u,i,j,k) (pm1o840dy*(378*(u)[index+di*(0)+dj*(0)+dk*(0)] - 60*(u)[index+di*(0)+dj*(-2*dir2)+dk*(0)] + 5*(u)[index+di*(0)+dj*(-3*dir2)+dk*(0)] - 140*(u)[index+di*(0)+dj*(3*dir2)+dk*(0)] + 30*(u)[index+di*(0)+dj*(4*dir2)+dk*(0)] - 3*(u)[index+di*(0)+dj*(5*dir2)+dk*(0)] + 420*((u)[index+di*(0)+dj*(2*dir2)+dk*(0)] + (u)[index+di*(0)+dj*(-dir2)+dk*(0)]) - 1050*(u)[index+di*(0)+dj*(dir2)+dk*(0)])*dir2) -#define PDupwindNth3(u,i,j,k) (pm1o840dz*(378*(u)[index+di*(0)+dj*(0)+dk*(0)] - 60*(u)[index+di*(0)+dj*(0)+dk*(-2*dir3)] + 5*(u)[index+di*(0)+dj*(0)+dk*(-3*dir3)] - 140*(u)[index+di*(0)+dj*(0)+dk*(3*dir3)] + 30*(u)[index+di*(0)+dj*(0)+dk*(4*dir3)] - 3*(u)[index+di*(0)+dj*(0)+dk*(5*dir3)] + 420*((u)[index+di*(0)+dj*(0)+dk*(2*dir3)] + (u)[index+di*(0)+dj*(0)+dk*(-dir3)]) - 1050*(u)[index+di*(0)+dj*(0)+dk*(dir3)])*dir3) -#define PDonesided1(u,i,j,k) (p1odx*(-(u)[index+di*(0)+dj*(0)+dk*(0)] + (u)[index+di*(dir1)+dj*(0)+dk*(0)])*dir1) -#define PDonesided2(u,i,j,k) (p1ody*(-(u)[index+di*(0)+dj*(0)+dk*(0)] + (u)[index+di*(0)+dj*(dir2)+dk*(0)])*dir2) -#define PDonesided3(u,i,j,k) (p1odz*(-(u)[index+di*(0)+dj*(0)+dk*(0)] + (u)[index+di*(0)+dj*(0)+dk*(dir3)])*dir3) -#define PDupwindNthAnti1(u,i,j,k) (p1o1680dx*(-1470*(u)[index+di*(-1)+dj*(0)+dk*(0)] + 1470*(u)[index+di*(1)+dj*(0)+dk*(0)] + 480*(u)[index+di*(-2)+dj*(0)+dk*(0)] - 480*(u)[index+di*(2)+dj*(0)+dk*(0)] - 145*(u)[index+di*(-3)+dj*(0)+dk*(0)] + 145*(u)[index+di*(3)+dj*(0)+dk*(0)] + 30*(u)[index+di*(-4)+dj*(0)+dk*(0)] - 30*(u)[index+di*(4)+dj*(0)+dk*(0)] - 3*(u)[index+di*(-5)+dj*(0)+dk*(0)] + 3*(u)[index+di*(5)+dj*(0)+dk*(0)])) -#define PDupwindNthSymm1(u,i,j,k) (p1o560dx*(-252*(u)[index+di*(0)+dj*(0)+dk*(0)] + 210*((u)[index+di*(-1)+dj*(0)+dk*(0)] + (u)[index+di*(1)+dj*(0)+dk*(0)]) - 120*((u)[index+di*(-2)+dj*(0)+dk*(0)] + (u)[index+di*(2)+dj*(0)+dk*(0)]) + 45*((u)[index+di*(-3)+dj*(0)+dk*(0)] + (u)[index+di*(3)+dj*(0)+dk*(0)]) - 10*((u)[index+di*(-4)+dj*(0)+dk*(0)] + (u)[index+di*(4)+dj*(0)+dk*(0)]) + (u)[index+di*(-5)+dj*(0)+dk*(0)] + (u)[index+di*(5)+dj*(0)+dk*(0)])) -#define PDupwindNthAnti2(u,i,j,k) (p1o1680dy*(-1470*(u)[index+di*(0)+dj*(-1)+dk*(0)] + 1470*(u)[index+di*(0)+dj*(1)+dk*(0)] + 480*(u)[index+di*(0)+dj*(-2)+dk*(0)] - 480*(u)[index+di*(0)+dj*(2)+dk*(0)] - 145*(u)[index+di*(0)+dj*(-3)+dk*(0)] + 145*(u)[index+di*(0)+dj*(3)+dk*(0)] + 30*(u)[index+di*(0)+dj*(-4)+dk*(0)] - 30*(u)[index+di*(0)+dj*(4)+dk*(0)] - 3*(u)[index+di*(0)+dj*(-5)+dk*(0)] + 3*(u)[index+di*(0)+dj*(5)+dk*(0)])) -#define PDupwindNthSymm2(u,i,j,k) (p1o560dy*(-252*(u)[index+di*(0)+dj*(0)+dk*(0)] + 210*((u)[index+di*(0)+dj*(-1)+dk*(0)] + (u)[index+di*(0)+dj*(1)+dk*(0)]) - 120*((u)[index+di*(0)+dj*(-2)+dk*(0)] + (u)[index+di*(0)+dj*(2)+dk*(0)]) + 45*((u)[index+di*(0)+dj*(-3)+dk*(0)] + (u)[index+di*(0)+dj*(3)+dk*(0)]) - 10*((u)[index+di*(0)+dj*(-4)+dk*(0)] + (u)[index+di*(0)+dj*(4)+dk*(0)]) + (u)[index+di*(0)+dj*(-5)+dk*(0)] + (u)[index+di*(0)+dj*(5)+dk*(0)])) -#define PDupwindNthAnti3(u,i,j,k) (p1o1680dz*(-1470*(u)[index+di*(0)+dj*(0)+dk*(-1)] + 1470*(u)[index+di*(0)+dj*(0)+dk*(1)] + 480*(u)[index+di*(0)+dj*(0)+dk*(-2)] - 480*(u)[index+di*(0)+dj*(0)+dk*(2)] - 145*(u)[index+di*(0)+dj*(0)+dk*(-3)] + 145*(u)[index+di*(0)+dj*(0)+dk*(3)] + 30*(u)[index+di*(0)+dj*(0)+dk*(-4)] - 30*(u)[index+di*(0)+dj*(0)+dk*(4)] - 3*(u)[index+di*(0)+dj*(0)+dk*(-5)] + 3*(u)[index+di*(0)+dj*(0)+dk*(5)])) -#define PDupwindNthSymm3(u,i,j,k) (p1o560dz*(-252*(u)[index+di*(0)+dj*(0)+dk*(0)] + 210*((u)[index+di*(0)+dj*(0)+dk*(-1)] + (u)[index+di*(0)+dj*(0)+dk*(1)]) - 120*((u)[index+di*(0)+dj*(0)+dk*(-2)] + (u)[index+di*(0)+dj*(0)+dk*(2)]) + 45*((u)[index+di*(0)+dj*(0)+dk*(-3)] + (u)[index+di*(0)+dj*(0)+dk*(3)]) - 10*((u)[index+di*(0)+dj*(0)+dk*(-4)] + (u)[index+di*(0)+dj*(0)+dk*(4)]) + (u)[index+di*(0)+dj*(0)+dk*(-5)] + (u)[index+di*(0)+dj*(0)+dk*(5)])) +#define PDstandardNth1(u) (p1o840dx*(-672*(u)[di*(-1)+dj*(0)+dk*(0)] + 672*(u)[di*(1)+dj*(0)+dk*(0)] + 168*(u)[di*(-2)+dj*(0)+dk*(0)] - 168*(u)[di*(2)+dj*(0)+dk*(0)] - 32*(u)[di*(-3)+dj*(0)+dk*(0)] + 32*(u)[di*(3)+dj*(0)+dk*(0)] + 3*(u)[di*(-4)+dj*(0)+dk*(0)] - 3*(u)[di*(4)+dj*(0)+dk*(0)])) +#define PDstandardNth2(u) (p1o840dy*(-672*(u)[di*(0)+dj*(-1)+dk*(0)] + 672*(u)[di*(0)+dj*(1)+dk*(0)] + 168*(u)[di*(0)+dj*(-2)+dk*(0)] - 168*(u)[di*(0)+dj*(2)+dk*(0)] - 32*(u)[di*(0)+dj*(-3)+dk*(0)] + 32*(u)[di*(0)+dj*(3)+dk*(0)] + 3*(u)[di*(0)+dj*(-4)+dk*(0)] - 3*(u)[di*(0)+dj*(4)+dk*(0)])) +#define PDstandardNth3(u) (p1o840dz*(-672*(u)[di*(0)+dj*(0)+dk*(-1)] + 672*(u)[di*(0)+dj*(0)+dk*(1)] + 168*(u)[di*(0)+dj*(0)+dk*(-2)] - 168*(u)[di*(0)+dj*(0)+dk*(2)] - 32*(u)[di*(0)+dj*(0)+dk*(-3)] + 32*(u)[di*(0)+dj*(0)+dk*(3)] + 3*(u)[di*(0)+dj*(0)+dk*(-4)] - 3*(u)[di*(0)+dj*(0)+dk*(4)])) +#define PDstandardNth11(u) (p1o5040dx2*(-14350*(u)[di*(0)+dj*(0)+dk*(0)] + 8064*((u)[di*(-1)+dj*(0)+dk*(0)] + (u)[di*(1)+dj*(0)+dk*(0)]) - 1008*((u)[di*(-2)+dj*(0)+dk*(0)] + (u)[di*(2)+dj*(0)+dk*(0)]) + 128*((u)[di*(-3)+dj*(0)+dk*(0)] + (u)[di*(3)+dj*(0)+dk*(0)]) - 9*((u)[di*(-4)+dj*(0)+dk*(0)] + (u)[di*(4)+dj*(0)+dk*(0)]))) +#define PDstandardNth22(u) (p1o5040dy2*(-14350*(u)[di*(0)+dj*(0)+dk*(0)] + 8064*((u)[di*(0)+dj*(-1)+dk*(0)] + (u)[di*(0)+dj*(1)+dk*(0)]) - 1008*((u)[di*(0)+dj*(-2)+dk*(0)] + (u)[di*(0)+dj*(2)+dk*(0)]) + 128*((u)[di*(0)+dj*(-3)+dk*(0)] + (u)[di*(0)+dj*(3)+dk*(0)]) - 9*((u)[di*(0)+dj*(-4)+dk*(0)] + (u)[di*(0)+dj*(4)+dk*(0)]))) +#define PDstandardNth33(u) (p1o5040dz2*(-14350*(u)[di*(0)+dj*(0)+dk*(0)] + 8064*((u)[di*(0)+dj*(0)+dk*(-1)] + (u)[di*(0)+dj*(0)+dk*(1)]) - 1008*((u)[di*(0)+dj*(0)+dk*(-2)] + (u)[di*(0)+dj*(0)+dk*(2)]) + 128*((u)[di*(0)+dj*(0)+dk*(-3)] + (u)[di*(0)+dj*(0)+dk*(3)]) - 9*((u)[di*(0)+dj*(0)+dk*(-4)] + (u)[di*(0)+dj*(0)+dk*(4)]))) +#define PDstandardNth12(u) (p1o705600dxdy*(-451584*((u)[di*(-1)+dj*(1)+dk*(0)] + (u)[di*(1)+dj*(-1)+dk*(0)]) + 451584*((u)[di*(-1)+dj*(-1)+dk*(0)] + (u)[di*(1)+dj*(1)+dk*(0)]) + 112896*((u)[di*(-1)+dj*(2)+dk*(0)] + (u)[di*(1)+dj*(-2)+dk*(0)] + (u)[di*(-2)+dj*(1)+dk*(0)] + (u)[di*(2)+dj*(-1)+dk*(0)]) - 112896*((u)[di*(-1)+dj*(-2)+dk*(0)] + (u)[di*(1)+dj*(2)+dk*(0)] + (u)[di*(-2)+dj*(-1)+dk*(0)] + (u)[di*(2)+dj*(1)+dk*(0)]) - 28224*((u)[di*(-2)+dj*(2)+dk*(0)] + (u)[di*(2)+dj*(-2)+dk*(0)]) + 28224*((u)[di*(-2)+dj*(-2)+dk*(0)] + (u)[di*(2)+dj*(2)+dk*(0)]) - 21504*((u)[di*(-1)+dj*(3)+dk*(0)] + (u)[di*(1)+dj*(-3)+dk*(0)] + (u)[di*(-3)+dj*(1)+dk*(0)] + (u)[di*(3)+dj*(-1)+dk*(0)]) + 21504*((u)[di*(-1)+dj*(-3)+dk*(0)] + (u)[di*(1)+dj*(3)+dk*(0)] + (u)[di*(-3)+dj*(-1)+dk*(0)] + (u)[di*(3)+dj*(1)+dk*(0)]) + 5376*((u)[di*(-2)+dj*(3)+dk*(0)] + (u)[di*(2)+dj*(-3)+dk*(0)] + (u)[di*(-3)+dj*(2)+dk*(0)] + (u)[di*(3)+dj*(-2)+dk*(0)]) - 5376*((u)[di*(-2)+dj*(-3)+dk*(0)] + (u)[di*(2)+dj*(3)+dk*(0)] + (u)[di*(-3)+dj*(-2)+dk*(0)] + (u)[di*(3)+dj*(2)+dk*(0)]) - 1024*((u)[di*(-3)+dj*(3)+dk*(0)] + (u)[di*(3)+dj*(-3)+dk*(0)]) + 1024*((u)[di*(-3)+dj*(-3)+dk*(0)] + (u)[di*(3)+dj*(3)+dk*(0)]) + 2016*((u)[di*(-1)+dj*(4)+dk*(0)] + (u)[di*(1)+dj*(-4)+dk*(0)] + (u)[di*(-4)+dj*(1)+dk*(0)] + (u)[di*(4)+dj*(-1)+dk*(0)]) - 2016*((u)[di*(-1)+dj*(-4)+dk*(0)] + (u)[di*(1)+dj*(4)+dk*(0)] + (u)[di*(-4)+dj*(-1)+dk*(0)] + (u)[di*(4)+dj*(1)+dk*(0)]) - 504*((u)[di*(-2)+dj*(4)+dk*(0)] + (u)[di*(2)+dj*(-4)+dk*(0)] + (u)[di*(-4)+dj*(2)+dk*(0)] + (u)[di*(4)+dj*(-2)+dk*(0)]) + 504*((u)[di*(-2)+dj*(-4)+dk*(0)] + (u)[di*(2)+dj*(4)+dk*(0)] + (u)[di*(-4)+dj*(-2)+dk*(0)] + (u)[di*(4)+dj*(2)+dk*(0)]) + 96*((u)[di*(-3)+dj*(4)+dk*(0)] + (u)[di*(3)+dj*(-4)+dk*(0)] + (u)[di*(-4)+dj*(3)+dk*(0)] + (u)[di*(4)+dj*(-3)+dk*(0)]) - 96*((u)[di*(-3)+dj*(-4)+dk*(0)] + (u)[di*(3)+dj*(4)+dk*(0)] + (u)[di*(-4)+dj*(-3)+dk*(0)] + (u)[di*(4)+dj*(3)+dk*(0)]) - 9*((u)[di*(-4)+dj*(4)+dk*(0)] + (u)[di*(4)+dj*(-4)+dk*(0)]) + 9*((u)[di*(-4)+dj*(-4)+dk*(0)] + (u)[di*(4)+dj*(4)+dk*(0)]))) +#define PDstandardNth13(u) (p1o705600dxdz*(-451584*((u)[di*(-1)+dj*(0)+dk*(1)] + (u)[di*(1)+dj*(0)+dk*(-1)]) + 451584*((u)[di*(-1)+dj*(0)+dk*(-1)] + (u)[di*(1)+dj*(0)+dk*(1)]) + 112896*((u)[di*(-1)+dj*(0)+dk*(2)] + (u)[di*(1)+dj*(0)+dk*(-2)] + (u)[di*(-2)+dj*(0)+dk*(1)] + (u)[di*(2)+dj*(0)+dk*(-1)]) - 112896*((u)[di*(-1)+dj*(0)+dk*(-2)] + (u)[di*(1)+dj*(0)+dk*(2)] + (u)[di*(-2)+dj*(0)+dk*(-1)] + (u)[di*(2)+dj*(0)+dk*(1)]) - 28224*((u)[di*(-2)+dj*(0)+dk*(2)] + (u)[di*(2)+dj*(0)+dk*(-2)]) + 28224*((u)[di*(-2)+dj*(0)+dk*(-2)] + (u)[di*(2)+dj*(0)+dk*(2)]) - 21504*((u)[di*(-1)+dj*(0)+dk*(3)] + (u)[di*(1)+dj*(0)+dk*(-3)] + (u)[di*(-3)+dj*(0)+dk*(1)] + (u)[di*(3)+dj*(0)+dk*(-1)]) + 21504*((u)[di*(-1)+dj*(0)+dk*(-3)] + (u)[di*(1)+dj*(0)+dk*(3)] + (u)[di*(-3)+dj*(0)+dk*(-1)] + (u)[di*(3)+dj*(0)+dk*(1)]) + 5376*((u)[di*(-2)+dj*(0)+dk*(3)] + (u)[di*(2)+dj*(0)+dk*(-3)] + (u)[di*(-3)+dj*(0)+dk*(2)] + (u)[di*(3)+dj*(0)+dk*(-2)]) - 5376*((u)[di*(-2)+dj*(0)+dk*(-3)] + (u)[di*(2)+dj*(0)+dk*(3)] + (u)[di*(-3)+dj*(0)+dk*(-2)] + (u)[di*(3)+dj*(0)+dk*(2)]) - 1024*((u)[di*(-3)+dj*(0)+dk*(3)] + (u)[di*(3)+dj*(0)+dk*(-3)]) + 1024*((u)[di*(-3)+dj*(0)+dk*(-3)] + (u)[di*(3)+dj*(0)+dk*(3)]) + 2016*((u)[di*(-1)+dj*(0)+dk*(4)] + (u)[di*(1)+dj*(0)+dk*(-4)] + (u)[di*(-4)+dj*(0)+dk*(1)] + (u)[di*(4)+dj*(0)+dk*(-1)]) - 2016*((u)[di*(-1)+dj*(0)+dk*(-4)] + (u)[di*(1)+dj*(0)+dk*(4)] + (u)[di*(-4)+dj*(0)+dk*(-1)] + (u)[di*(4)+dj*(0)+dk*(1)]) - 504*((u)[di*(-2)+dj*(0)+dk*(4)] + (u)[di*(2)+dj*(0)+dk*(-4)] + (u)[di*(-4)+dj*(0)+dk*(2)] + (u)[di*(4)+dj*(0)+dk*(-2)]) + 504*((u)[di*(-2)+dj*(0)+dk*(-4)] + (u)[di*(2)+dj*(0)+dk*(4)] + (u)[di*(-4)+dj*(0)+dk*(-2)] + (u)[di*(4)+dj*(0)+dk*(2)]) + 96*((u)[di*(-3)+dj*(0)+dk*(4)] + (u)[di*(3)+dj*(0)+dk*(-4)] + (u)[di*(-4)+dj*(0)+dk*(3)] + (u)[di*(4)+dj*(0)+dk*(-3)]) - 96*((u)[di*(-3)+dj*(0)+dk*(-4)] + (u)[di*(3)+dj*(0)+dk*(4)] + (u)[di*(-4)+dj*(0)+dk*(-3)] + (u)[di*(4)+dj*(0)+dk*(3)]) - 9*((u)[di*(-4)+dj*(0)+dk*(4)] + (u)[di*(4)+dj*(0)+dk*(-4)]) + 9*((u)[di*(-4)+dj*(0)+dk*(-4)] + (u)[di*(4)+dj*(0)+dk*(4)]))) +#define PDstandardNth21(u) (p1o705600dxdy*(-451584*((u)[di*(-1)+dj*(1)+dk*(0)] + (u)[di*(1)+dj*(-1)+dk*(0)]) + 451584*((u)[di*(-1)+dj*(-1)+dk*(0)] + (u)[di*(1)+dj*(1)+dk*(0)]) + 112896*((u)[di*(-1)+dj*(2)+dk*(0)] + (u)[di*(1)+dj*(-2)+dk*(0)] + (u)[di*(-2)+dj*(1)+dk*(0)] + (u)[di*(2)+dj*(-1)+dk*(0)]) - 112896*((u)[di*(-1)+dj*(-2)+dk*(0)] + (u)[di*(1)+dj*(2)+dk*(0)] + (u)[di*(-2)+dj*(-1)+dk*(0)] + (u)[di*(2)+dj*(1)+dk*(0)]) - 28224*((u)[di*(-2)+dj*(2)+dk*(0)] + (u)[di*(2)+dj*(-2)+dk*(0)]) + 28224*((u)[di*(-2)+dj*(-2)+dk*(0)] + (u)[di*(2)+dj*(2)+dk*(0)]) - 21504*((u)[di*(-1)+dj*(3)+dk*(0)] + (u)[di*(1)+dj*(-3)+dk*(0)] + (u)[di*(-3)+dj*(1)+dk*(0)] + (u)[di*(3)+dj*(-1)+dk*(0)]) + 21504*((u)[di*(-1)+dj*(-3)+dk*(0)] + (u)[di*(1)+dj*(3)+dk*(0)] + (u)[di*(-3)+dj*(-1)+dk*(0)] + (u)[di*(3)+dj*(1)+dk*(0)]) + 5376*((u)[di*(-2)+dj*(3)+dk*(0)] + (u)[di*(2)+dj*(-3)+dk*(0)] + (u)[di*(-3)+dj*(2)+dk*(0)] + (u)[di*(3)+dj*(-2)+dk*(0)]) - 5376*((u)[di*(-2)+dj*(-3)+dk*(0)] + (u)[di*(2)+dj*(3)+dk*(0)] + (u)[di*(-3)+dj*(-2)+dk*(0)] + (u)[di*(3)+dj*(2)+dk*(0)]) - 1024*((u)[di*(-3)+dj*(3)+dk*(0)] + (u)[di*(3)+dj*(-3)+dk*(0)]) + 1024*((u)[di*(-3)+dj*(-3)+dk*(0)] + (u)[di*(3)+dj*(3)+dk*(0)]) + 2016*((u)[di*(-1)+dj*(4)+dk*(0)] + (u)[di*(1)+dj*(-4)+dk*(0)] + (u)[di*(-4)+dj*(1)+dk*(0)] + (u)[di*(4)+dj*(-1)+dk*(0)]) - 2016*((u)[di*(-1)+dj*(-4)+dk*(0)] + (u)[di*(1)+dj*(4)+dk*(0)] + (u)[di*(-4)+dj*(-1)+dk*(0)] + (u)[di*(4)+dj*(1)+dk*(0)]) - 504*((u)[di*(-2)+dj*(4)+dk*(0)] + (u)[di*(2)+dj*(-4)+dk*(0)] + (u)[di*(-4)+dj*(2)+dk*(0)] + (u)[di*(4)+dj*(-2)+dk*(0)]) + 504*((u)[di*(-2)+dj*(-4)+dk*(0)] + (u)[di*(2)+dj*(4)+dk*(0)] + (u)[di*(-4)+dj*(-2)+dk*(0)] + (u)[di*(4)+dj*(2)+dk*(0)]) + 96*((u)[di*(-3)+dj*(4)+dk*(0)] + (u)[di*(3)+dj*(-4)+dk*(0)] + (u)[di*(-4)+dj*(3)+dk*(0)] + (u)[di*(4)+dj*(-3)+dk*(0)]) - 96*((u)[di*(-3)+dj*(-4)+dk*(0)] + (u)[di*(3)+dj*(4)+dk*(0)] + (u)[di*(-4)+dj*(-3)+dk*(0)] + (u)[di*(4)+dj*(3)+dk*(0)]) - 9*((u)[di*(-4)+dj*(4)+dk*(0)] + (u)[di*(4)+dj*(-4)+dk*(0)]) + 9*((u)[di*(-4)+dj*(-4)+dk*(0)] + (u)[di*(4)+dj*(4)+dk*(0)]))) +#define PDstandardNth23(u) (p1o705600dydz*(-451584*((u)[di*(0)+dj*(-1)+dk*(1)] + (u)[di*(0)+dj*(1)+dk*(-1)]) + 451584*((u)[di*(0)+dj*(-1)+dk*(-1)] + (u)[di*(0)+dj*(1)+dk*(1)]) + 112896*((u)[di*(0)+dj*(-1)+dk*(2)] + (u)[di*(0)+dj*(1)+dk*(-2)] + (u)[di*(0)+dj*(-2)+dk*(1)] + (u)[di*(0)+dj*(2)+dk*(-1)]) - 112896*((u)[di*(0)+dj*(-1)+dk*(-2)] + (u)[di*(0)+dj*(1)+dk*(2)] + (u)[di*(0)+dj*(-2)+dk*(-1)] + (u)[di*(0)+dj*(2)+dk*(1)]) - 28224*((u)[di*(0)+dj*(-2)+dk*(2)] + (u)[di*(0)+dj*(2)+dk*(-2)]) + 28224*((u)[di*(0)+dj*(-2)+dk*(-2)] + (u)[di*(0)+dj*(2)+dk*(2)]) - 21504*((u)[di*(0)+dj*(-1)+dk*(3)] + (u)[di*(0)+dj*(1)+dk*(-3)] + (u)[di*(0)+dj*(-3)+dk*(1)] + (u)[di*(0)+dj*(3)+dk*(-1)]) + 21504*((u)[di*(0)+dj*(-1)+dk*(-3)] + (u)[di*(0)+dj*(1)+dk*(3)] + (u)[di*(0)+dj*(-3)+dk*(-1)] + (u)[di*(0)+dj*(3)+dk*(1)]) + 5376*((u)[di*(0)+dj*(-2)+dk*(3)] + (u)[di*(0)+dj*(2)+dk*(-3)] + (u)[di*(0)+dj*(-3)+dk*(2)] + (u)[di*(0)+dj*(3)+dk*(-2)]) - 5376*((u)[di*(0)+dj*(-2)+dk*(-3)] + (u)[di*(0)+dj*(2)+dk*(3)] + (u)[di*(0)+dj*(-3)+dk*(-2)] + (u)[di*(0)+dj*(3)+dk*(2)]) - 1024*((u)[di*(0)+dj*(-3)+dk*(3)] + (u)[di*(0)+dj*(3)+dk*(-3)]) + 1024*((u)[di*(0)+dj*(-3)+dk*(-3)] + (u)[di*(0)+dj*(3)+dk*(3)]) + 2016*((u)[di*(0)+dj*(-1)+dk*(4)] + (u)[di*(0)+dj*(1)+dk*(-4)] + (u)[di*(0)+dj*(-4)+dk*(1)] + (u)[di*(0)+dj*(4)+dk*(-1)]) - 2016*((u)[di*(0)+dj*(-1)+dk*(-4)] + (u)[di*(0)+dj*(1)+dk*(4)] + (u)[di*(0)+dj*(-4)+dk*(-1)] + (u)[di*(0)+dj*(4)+dk*(1)]) - 504*((u)[di*(0)+dj*(-2)+dk*(4)] + (u)[di*(0)+dj*(2)+dk*(-4)] + (u)[di*(0)+dj*(-4)+dk*(2)] + (u)[di*(0)+dj*(4)+dk*(-2)]) + 504*((u)[di*(0)+dj*(-2)+dk*(-4)] + (u)[di*(0)+dj*(2)+dk*(4)] + (u)[di*(0)+dj*(-4)+dk*(-2)] + (u)[di*(0)+dj*(4)+dk*(2)]) + 96*((u)[di*(0)+dj*(-3)+dk*(4)] + (u)[di*(0)+dj*(3)+dk*(-4)] + (u)[di*(0)+dj*(-4)+dk*(3)] + (u)[di*(0)+dj*(4)+dk*(-3)]) - 96*((u)[di*(0)+dj*(-3)+dk*(-4)] + (u)[di*(0)+dj*(3)+dk*(4)] + (u)[di*(0)+dj*(-4)+dk*(-3)] + (u)[di*(0)+dj*(4)+dk*(3)]) - 9*((u)[di*(0)+dj*(-4)+dk*(4)] + (u)[di*(0)+dj*(4)+dk*(-4)]) + 9*((u)[di*(0)+dj*(-4)+dk*(-4)] + (u)[di*(0)+dj*(4)+dk*(4)]))) +#define PDstandardNth31(u) (p1o705600dxdz*(-451584*((u)[di*(-1)+dj*(0)+dk*(1)] + (u)[di*(1)+dj*(0)+dk*(-1)]) + 451584*((u)[di*(-1)+dj*(0)+dk*(-1)] + (u)[di*(1)+dj*(0)+dk*(1)]) + 112896*((u)[di*(-1)+dj*(0)+dk*(2)] + (u)[di*(1)+dj*(0)+dk*(-2)] + (u)[di*(-2)+dj*(0)+dk*(1)] + (u)[di*(2)+dj*(0)+dk*(-1)]) - 112896*((u)[di*(-1)+dj*(0)+dk*(-2)] + (u)[di*(1)+dj*(0)+dk*(2)] + (u)[di*(-2)+dj*(0)+dk*(-1)] + (u)[di*(2)+dj*(0)+dk*(1)]) - 28224*((u)[di*(-2)+dj*(0)+dk*(2)] + (u)[di*(2)+dj*(0)+dk*(-2)]) + 28224*((u)[di*(-2)+dj*(0)+dk*(-2)] + (u)[di*(2)+dj*(0)+dk*(2)]) - 21504*((u)[di*(-1)+dj*(0)+dk*(3)] + (u)[di*(1)+dj*(0)+dk*(-3)] + (u)[di*(-3)+dj*(0)+dk*(1)] + (u)[di*(3)+dj*(0)+dk*(-1)]) + 21504*((u)[di*(-1)+dj*(0)+dk*(-3)] + (u)[di*(1)+dj*(0)+dk*(3)] + (u)[di*(-3)+dj*(0)+dk*(-1)] + (u)[di*(3)+dj*(0)+dk*(1)]) + 5376*((u)[di*(-2)+dj*(0)+dk*(3)] + (u)[di*(2)+dj*(0)+dk*(-3)] + (u)[di*(-3)+dj*(0)+dk*(2)] + (u)[di*(3)+dj*(0)+dk*(-2)]) - 5376*((u)[di*(-2)+dj*(0)+dk*(-3)] + (u)[di*(2)+dj*(0)+dk*(3)] + (u)[di*(-3)+dj*(0)+dk*(-2)] + (u)[di*(3)+dj*(0)+dk*(2)]) - 1024*((u)[di*(-3)+dj*(0)+dk*(3)] + (u)[di*(3)+dj*(0)+dk*(-3)]) + 1024*((u)[di*(-3)+dj*(0)+dk*(-3)] + (u)[di*(3)+dj*(0)+dk*(3)]) + 2016*((u)[di*(-1)+dj*(0)+dk*(4)] + (u)[di*(1)+dj*(0)+dk*(-4)] + (u)[di*(-4)+dj*(0)+dk*(1)] + (u)[di*(4)+dj*(0)+dk*(-1)]) - 2016*((u)[di*(-1)+dj*(0)+dk*(-4)] + (u)[di*(1)+dj*(0)+dk*(4)] + (u)[di*(-4)+dj*(0)+dk*(-1)] + (u)[di*(4)+dj*(0)+dk*(1)]) - 504*((u)[di*(-2)+dj*(0)+dk*(4)] + (u)[di*(2)+dj*(0)+dk*(-4)] + (u)[di*(-4)+dj*(0)+dk*(2)] + (u)[di*(4)+dj*(0)+dk*(-2)]) + 504*((u)[di*(-2)+dj*(0)+dk*(-4)] + (u)[di*(2)+dj*(0)+dk*(4)] + (u)[di*(-4)+dj*(0)+dk*(-2)] + (u)[di*(4)+dj*(0)+dk*(2)]) + 96*((u)[di*(-3)+dj*(0)+dk*(4)] + (u)[di*(3)+dj*(0)+dk*(-4)] + (u)[di*(-4)+dj*(0)+dk*(3)] + (u)[di*(4)+dj*(0)+dk*(-3)]) - 96*((u)[di*(-3)+dj*(0)+dk*(-4)] + (u)[di*(3)+dj*(0)+dk*(4)] + (u)[di*(-4)+dj*(0)+dk*(-3)] + (u)[di*(4)+dj*(0)+dk*(3)]) - 9*((u)[di*(-4)+dj*(0)+dk*(4)] + (u)[di*(4)+dj*(0)+dk*(-4)]) + 9*((u)[di*(-4)+dj*(0)+dk*(-4)] + (u)[di*(4)+dj*(0)+dk*(4)]))) +#define PDstandardNth32(u) (p1o705600dydz*(-451584*((u)[di*(0)+dj*(-1)+dk*(1)] + (u)[di*(0)+dj*(1)+dk*(-1)]) + 451584*((u)[di*(0)+dj*(-1)+dk*(-1)] + (u)[di*(0)+dj*(1)+dk*(1)]) + 112896*((u)[di*(0)+dj*(-1)+dk*(2)] + (u)[di*(0)+dj*(1)+dk*(-2)] + (u)[di*(0)+dj*(-2)+dk*(1)] + (u)[di*(0)+dj*(2)+dk*(-1)]) - 112896*((u)[di*(0)+dj*(-1)+dk*(-2)] + (u)[di*(0)+dj*(1)+dk*(2)] + (u)[di*(0)+dj*(-2)+dk*(-1)] + (u)[di*(0)+dj*(2)+dk*(1)]) - 28224*((u)[di*(0)+dj*(-2)+dk*(2)] + (u)[di*(0)+dj*(2)+dk*(-2)]) + 28224*((u)[di*(0)+dj*(-2)+dk*(-2)] + (u)[di*(0)+dj*(2)+dk*(2)]) - 21504*((u)[di*(0)+dj*(-1)+dk*(3)] + (u)[di*(0)+dj*(1)+dk*(-3)] + (u)[di*(0)+dj*(-3)+dk*(1)] + (u)[di*(0)+dj*(3)+dk*(-1)]) + 21504*((u)[di*(0)+dj*(-1)+dk*(-3)] + (u)[di*(0)+dj*(1)+dk*(3)] + (u)[di*(0)+dj*(-3)+dk*(-1)] + (u)[di*(0)+dj*(3)+dk*(1)]) + 5376*((u)[di*(0)+dj*(-2)+dk*(3)] + (u)[di*(0)+dj*(2)+dk*(-3)] + (u)[di*(0)+dj*(-3)+dk*(2)] + (u)[di*(0)+dj*(3)+dk*(-2)]) - 5376*((u)[di*(0)+dj*(-2)+dk*(-3)] + (u)[di*(0)+dj*(2)+dk*(3)] + (u)[di*(0)+dj*(-3)+dk*(-2)] + (u)[di*(0)+dj*(3)+dk*(2)]) - 1024*((u)[di*(0)+dj*(-3)+dk*(3)] + (u)[di*(0)+dj*(3)+dk*(-3)]) + 1024*((u)[di*(0)+dj*(-3)+dk*(-3)] + (u)[di*(0)+dj*(3)+dk*(3)]) + 2016*((u)[di*(0)+dj*(-1)+dk*(4)] + (u)[di*(0)+dj*(1)+dk*(-4)] + (u)[di*(0)+dj*(-4)+dk*(1)] + (u)[di*(0)+dj*(4)+dk*(-1)]) - 2016*((u)[di*(0)+dj*(-1)+dk*(-4)] + (u)[di*(0)+dj*(1)+dk*(4)] + (u)[di*(0)+dj*(-4)+dk*(-1)] + (u)[di*(0)+dj*(4)+dk*(1)]) - 504*((u)[di*(0)+dj*(-2)+dk*(4)] + (u)[di*(0)+dj*(2)+dk*(-4)] + (u)[di*(0)+dj*(-4)+dk*(2)] + (u)[di*(0)+dj*(4)+dk*(-2)]) + 504*((u)[di*(0)+dj*(-2)+dk*(-4)] + (u)[di*(0)+dj*(2)+dk*(4)] + (u)[di*(0)+dj*(-4)+dk*(-2)] + (u)[di*(0)+dj*(4)+dk*(2)]) + 96*((u)[di*(0)+dj*(-3)+dk*(4)] + (u)[di*(0)+dj*(3)+dk*(-4)] + (u)[di*(0)+dj*(-4)+dk*(3)] + (u)[di*(0)+dj*(4)+dk*(-3)]) - 96*((u)[di*(0)+dj*(-3)+dk*(-4)] + (u)[di*(0)+dj*(3)+dk*(4)] + (u)[di*(0)+dj*(-4)+dk*(-3)] + (u)[di*(0)+dj*(4)+dk*(3)]) - 9*((u)[di*(0)+dj*(-4)+dk*(4)] + (u)[di*(0)+dj*(4)+dk*(-4)]) + 9*((u)[di*(0)+dj*(-4)+dk*(-4)] + (u)[di*(0)+dj*(4)+dk*(4)]))) +#define PDdissipationNth1(u) (p1o1024dx*(-252*(u)[di*(0)+dj*(0)+dk*(0)] + 210*((u)[di*(-1)+dj*(0)+dk*(0)] + (u)[di*(1)+dj*(0)+dk*(0)]) - 120*((u)[di*(-2)+dj*(0)+dk*(0)] + (u)[di*(2)+dj*(0)+dk*(0)]) + 45*((u)[di*(-3)+dj*(0)+dk*(0)] + (u)[di*(3)+dj*(0)+dk*(0)]) - 10*((u)[di*(-4)+dj*(0)+dk*(0)] + (u)[di*(4)+dj*(0)+dk*(0)]) + (u)[di*(-5)+dj*(0)+dk*(0)] + (u)[di*(5)+dj*(0)+dk*(0)])) +#define PDdissipationNth2(u) (p1o1024dy*(-252*(u)[di*(0)+dj*(0)+dk*(0)] + 210*((u)[di*(0)+dj*(-1)+dk*(0)] + (u)[di*(0)+dj*(1)+dk*(0)]) - 120*((u)[di*(0)+dj*(-2)+dk*(0)] + (u)[di*(0)+dj*(2)+dk*(0)]) + 45*((u)[di*(0)+dj*(-3)+dk*(0)] + (u)[di*(0)+dj*(3)+dk*(0)]) - 10*((u)[di*(0)+dj*(-4)+dk*(0)] + (u)[di*(0)+dj*(4)+dk*(0)]) + (u)[di*(0)+dj*(-5)+dk*(0)] + (u)[di*(0)+dj*(5)+dk*(0)])) +#define PDdissipationNth3(u) (p1o1024dz*(-252*(u)[di*(0)+dj*(0)+dk*(0)] + 210*((u)[di*(0)+dj*(0)+dk*(-1)] + (u)[di*(0)+dj*(0)+dk*(1)]) - 120*((u)[di*(0)+dj*(0)+dk*(-2)] + (u)[di*(0)+dj*(0)+dk*(2)]) + 45*((u)[di*(0)+dj*(0)+dk*(-3)] + (u)[di*(0)+dj*(0)+dk*(3)]) - 10*((u)[di*(0)+dj*(0)+dk*(-4)] + (u)[di*(0)+dj*(0)+dk*(4)]) + (u)[di*(0)+dj*(0)+dk*(-5)] + (u)[di*(0)+dj*(0)+dk*(5)])) +#define PDupwindNth1(u) (pm1o840dx*(378*(u)[di*(0)+dj*(0)+dk*(0)] - 60*(u)[di*(-2*dir1)+dj*(0)+dk*(0)] + 5*(u)[di*(-3*dir1)+dj*(0)+dk*(0)] - 140*(u)[di*(3*dir1)+dj*(0)+dk*(0)] + 30*(u)[di*(4*dir1)+dj*(0)+dk*(0)] - 3*(u)[di*(5*dir1)+dj*(0)+dk*(0)] + 420*((u)[di*(2*dir1)+dj*(0)+dk*(0)] + (u)[di*(-dir1)+dj*(0)+dk*(0)]) - 1050*(u)[di*(dir1)+dj*(0)+dk*(0)])*dir1) +#define PDupwindNth2(u) (pm1o840dy*(378*(u)[di*(0)+dj*(0)+dk*(0)] - 60*(u)[di*(0)+dj*(-2*dir2)+dk*(0)] + 5*(u)[di*(0)+dj*(-3*dir2)+dk*(0)] - 140*(u)[di*(0)+dj*(3*dir2)+dk*(0)] + 30*(u)[di*(0)+dj*(4*dir2)+dk*(0)] - 3*(u)[di*(0)+dj*(5*dir2)+dk*(0)] + 420*((u)[di*(0)+dj*(2*dir2)+dk*(0)] + (u)[di*(0)+dj*(-dir2)+dk*(0)]) - 1050*(u)[di*(0)+dj*(dir2)+dk*(0)])*dir2) +#define PDupwindNth3(u) (pm1o840dz*(378*(u)[di*(0)+dj*(0)+dk*(0)] - 60*(u)[di*(0)+dj*(0)+dk*(-2*dir3)] + 5*(u)[di*(0)+dj*(0)+dk*(-3*dir3)] - 140*(u)[di*(0)+dj*(0)+dk*(3*dir3)] + 30*(u)[di*(0)+dj*(0)+dk*(4*dir3)] - 3*(u)[di*(0)+dj*(0)+dk*(5*dir3)] + 420*((u)[di*(0)+dj*(0)+dk*(2*dir3)] + (u)[di*(0)+dj*(0)+dk*(-dir3)]) - 1050*(u)[di*(0)+dj*(0)+dk*(dir3)])*dir3) +#define PDonesided1(u) (p1odx*(-(u)[di*(0)+dj*(0)+dk*(0)] + (u)[di*(dir1)+dj*(0)+dk*(0)])*dir1) +#define PDonesided2(u) (p1ody*(-(u)[di*(0)+dj*(0)+dk*(0)] + (u)[di*(0)+dj*(dir2)+dk*(0)])*dir2) +#define PDonesided3(u) (p1odz*(-(u)[di*(0)+dj*(0)+dk*(0)] + (u)[di*(0)+dj*(0)+dk*(dir3)])*dir3) +#define PDupwindNthAnti1(u) (p1o1680dx*(-1470*(u)[di*(-1)+dj*(0)+dk*(0)] + 1470*(u)[di*(1)+dj*(0)+dk*(0)] + 480*(u)[di*(-2)+dj*(0)+dk*(0)] - 480*(u)[di*(2)+dj*(0)+dk*(0)] - 145*(u)[di*(-3)+dj*(0)+dk*(0)] + 145*(u)[di*(3)+dj*(0)+dk*(0)] + 30*(u)[di*(-4)+dj*(0)+dk*(0)] - 30*(u)[di*(4)+dj*(0)+dk*(0)] - 3*(u)[di*(-5)+dj*(0)+dk*(0)] + 3*(u)[di*(5)+dj*(0)+dk*(0)])) +#define PDupwindNthSymm1(u) (p1o560dx*(-252*(u)[di*(0)+dj*(0)+dk*(0)] + 210*((u)[di*(-1)+dj*(0)+dk*(0)] + (u)[di*(1)+dj*(0)+dk*(0)]) - 120*((u)[di*(-2)+dj*(0)+dk*(0)] + (u)[di*(2)+dj*(0)+dk*(0)]) + 45*((u)[di*(-3)+dj*(0)+dk*(0)] + (u)[di*(3)+dj*(0)+dk*(0)]) - 10*((u)[di*(-4)+dj*(0)+dk*(0)] + (u)[di*(4)+dj*(0)+dk*(0)]) + (u)[di*(-5)+dj*(0)+dk*(0)] + (u)[di*(5)+dj*(0)+dk*(0)])) +#define PDupwindNthAnti2(u) (p1o1680dy*(-1470*(u)[di*(0)+dj*(-1)+dk*(0)] + 1470*(u)[di*(0)+dj*(1)+dk*(0)] + 480*(u)[di*(0)+dj*(-2)+dk*(0)] - 480*(u)[di*(0)+dj*(2)+dk*(0)] - 145*(u)[di*(0)+dj*(-3)+dk*(0)] + 145*(u)[di*(0)+dj*(3)+dk*(0)] + 30*(u)[di*(0)+dj*(-4)+dk*(0)] - 30*(u)[di*(0)+dj*(4)+dk*(0)] - 3*(u)[di*(0)+dj*(-5)+dk*(0)] + 3*(u)[di*(0)+dj*(5)+dk*(0)])) +#define PDupwindNthSymm2(u) (p1o560dy*(-252*(u)[di*(0)+dj*(0)+dk*(0)] + 210*((u)[di*(0)+dj*(-1)+dk*(0)] + (u)[di*(0)+dj*(1)+dk*(0)]) - 120*((u)[di*(0)+dj*(-2)+dk*(0)] + (u)[di*(0)+dj*(2)+dk*(0)]) + 45*((u)[di*(0)+dj*(-3)+dk*(0)] + (u)[di*(0)+dj*(3)+dk*(0)]) - 10*((u)[di*(0)+dj*(-4)+dk*(0)] + (u)[di*(0)+dj*(4)+dk*(0)]) + (u)[di*(0)+dj*(-5)+dk*(0)] + (u)[di*(0)+dj*(5)+dk*(0)])) +#define PDupwindNthAnti3(u) (p1o1680dz*(-1470*(u)[di*(0)+dj*(0)+dk*(-1)] + 1470*(u)[di*(0)+dj*(0)+dk*(1)] + 480*(u)[di*(0)+dj*(0)+dk*(-2)] - 480*(u)[di*(0)+dj*(0)+dk*(2)] - 145*(u)[di*(0)+dj*(0)+dk*(-3)] + 145*(u)[di*(0)+dj*(0)+dk*(3)] + 30*(u)[di*(0)+dj*(0)+dk*(-4)] - 30*(u)[di*(0)+dj*(0)+dk*(4)] - 3*(u)[di*(0)+dj*(0)+dk*(-5)] + 3*(u)[di*(0)+dj*(0)+dk*(5)])) +#define PDupwindNthSymm3(u) (p1o560dz*(-252*(u)[di*(0)+dj*(0)+dk*(0)] + 210*((u)[di*(0)+dj*(0)+dk*(-1)] + (u)[di*(0)+dj*(0)+dk*(1)]) - 120*((u)[di*(0)+dj*(0)+dk*(-2)] + (u)[di*(0)+dj*(0)+dk*(2)]) + 45*((u)[di*(0)+dj*(0)+dk*(-3)] + (u)[di*(0)+dj*(0)+dk*(3)]) - 10*((u)[di*(0)+dj*(0)+dk*(-4)] + (u)[di*(0)+dj*(0)+dk*(4)]) + (u)[di*(0)+dj*(0)+dk*(-5)] + (u)[di*(0)+dj*(0)+dk*(5)])) |