diff options
Diffstat (limited to 'src/macro/ADM_Derivative.h')
-rw-r--r-- | src/macro/ADM_Derivative.h | 39 |
1 files changed, 21 insertions, 18 deletions
diff --git a/src/macro/ADM_Derivative.h b/src/macro/ADM_Derivative.h index c9b4fb4..8e62310 100644 --- a/src/macro/ADM_Derivative.h +++ b/src/macro/ADM_Derivative.h @@ -28,40 +28,43 @@ + var(i,j,k-1)) #define ADM_DXY_2(var,i,j,k) IDXY*(var(i+1,j+1,k) - var(i+1,j-1,k) \ - - var(i-1,j+1,k) + var(i-1,j-1,k)) + - var(i-1,j+1,k) + var(i-1,j-1,k)) #define ADM_DXZ_2(var,i,j,k) IDXZ*(var(i+1,j,k+1) - var(i+1,j,k-1) \ - - var(i-1,j,k+1) + var(i-1,j,k-1)) + - var(i-1,j,k+1) + var(i-1,j,k-1)) #define ADM_DYZ_2(var,i,j,k) IDYZ*(var(i,j+1,k+1) - var(i,j+1,k-1) \ - - var(i,j-1,k+1) + var(i,j-1,k-1)) + - var(i,j-1,k+1) + var(i,j-1,k-1)) /* * 4th order operators. */ -#define ADM_DX_4(var,i,j,k) (I2DX/6.d0)*(-var(i+2,j,k) + 8.d0*var(i+1,j,k) \ - - 8.d0*var(i-1,j,k) + var(i-2,j,k)) -#define ADM_DY_4(var,i,j,k) (I2DY/6.d0)*(-var(i,j+2,k) + 8.d0*var(i,j+1,k) \ - - 8.d0*var(i,j-1,k) + var(i,j-2,k)) -#define ADM_DZ_4(var,i,j,k) (I2DZ/6.d0)*(-var(i,j,k+2) + 8.d0*var(i,j,k+1) \ - - 8.d0*var(i,j,k-1) + var(i,j,k-2)) +#define ADM_DX_4(var,i,j,k) I12DX*(-var(i+2,j,k) + var(i-2,j,k) \ + + 8.d0*(var(i+1,j,k) - var(i-1,j,k))) +#define ADM_DY_4(var,i,j,k) I12DY*(-var(i,j+2,k) + var(i,j-2,k) \ + + 8.d0*(var(i,j+1,k) - var(i,j-1,k))) +#define ADM_DZ_4(var,i,j,k) I12DZ*(-var(i,j,k+2) + var(i,j,k-2) \ + + 8.d0*(var(i,j,k+1) - var(i,j,k-1))) -#define ADM_DXX_4(var,i,j,k) (IDXX/12.d0)*(-var(i+2,j,k) + 16.d0*var(i+1,j,k)\ - - 30.d0*var(i,j,k) + 16.d0*var(i-1,j,k) - var(i-2,j,k)) -#define ADM_DYY_4(var,i,j,k) (IDYY/12.d0)*(-var(i,j+2,k) + 16.d0*var(i,j+1,k)\ - - 30.d0*var(i,j,k) + 16.d0*var(i,j-1,k) - var(i,j-2,k)) -#define ADM_DZZ_4(var,i,j,k) (IDZZ/12.d0)*(-var(i,j,k+2) + 16.d0*var(i,j,k+1)\ - - 30.d0*var(i,j,k) + 16.d0*var(i,j,k-1) - var(i,j,k-2)) +#define ADM_DXX_4(var,i,j,k) I12DXX*(-var(i+2,j,k) - var(i-2,j,k) \ + + 16.d0*(var(i+1,j,k) + var(i-1,j,k)) \ + - 30.d0*var(i,j,k)) +#define ADM_DYY_4(var,i,j,k) I12DYY*(-var(i,j+2,k) - var(i,j-2,k) \ + + 16.d0*(var(i,j+1,k) + var(i,j-1,k)) \ + - 30.d0*var(i,j,k)) +#define ADM_DZZ_4(var,i,j,k) I12DZZ*(-var(i,j,k+2) - var(i,j,k-2) \ + + 16.d0*(var(i,j,k+1) + var(i,j,k-1)) \ + - 30.d0*var(i,j,k)) -#define ADM_DXY_4(var,i,j,k) (IDXY/36.d0)* \ +#define ADM_DXY_4(var,i,j,k) I36DXY* \ (var(i+2,j+2,k) - var(i+2,j-2,k) - var(i-2,j+2,k) + var(i-2,j-2,k) \ + 8.d0*(-var(i+2,j+1,k) + var(i+2,j-1,k) - var(i+1,j+2,k) + var(i+1,j-2,k) \ + var(i-2,j+1,k) - var(i-2,j-1,k) - var(i-1,j-2,k) + var(i-1,j+2,k)) \ + 64.d0*(var(i+1,j+1,k) - var(i+1,j-1,k) - var(i-1,j+1,k) + var(i-1,j-1,k))) -#define ADM_DXZ_4(var,i,j,k) (IDXZ/36.d0)* \ +#define ADM_DXZ_4(var,i,j,k) I36DXZ* \ (var(i+2,j,k+2) - var(i+2,j,k-2) - var(i-2,j,k+2) + var(i-2,j,k-2) \ + 8.d0*(-var(i+2,j,k+1) + var(i+2,j,k-1) - var(i+1,j,k+2) + var(i+1,j,k-2) \ + var(i-2,j,k+1) - var(i-2,j,k-1) - var(i-1,j,k-2) + var(i-1,j,k+2)) \ + 64.d0*(var(i+1,j,k+1) - var(i+1,j,k-1) - var(i-1,j,k+1) + var(i-1,j,k-1))) -#define ADM_DYZ_4(var,i,j,k) (IDYZ/36.d0)* \ +#define ADM_DYZ_4(var,i,j,k) I36DYZ* \ (var(i,j+2,k+2) - var(i,j+2,k-2) - var(i,j-2,k+2) + var(i,j-2,k-2) \ + 8.d0*(-var(i,j+2,k+1) + var(i,j+2,k-1) - var(i,j+1,k+2) + var(i,j+1,k-2) \ + var(i,j-2,k+1) - var(i,j-2,k-1) - var(i,j-1,k-2) + var(i,j-1,k+2)) \ |