diff options
| author | pollney <pollney@b1d164ef-f17a-46e7-89d4-021c7118ef4e> | 2003-06-26 11:37:17 +0000 |
|---|---|---|
| committer | pollney <pollney@b1d164ef-f17a-46e7-89d4-021c7118ef4e> | 2003-06-26 11:37:17 +0000 |
| commit | f6be10b987d8db12434c2607dcd6d434a277157b (patch) | |
| tree | bdfcb7c3758a15d492447b34e2fd727e2ab56aa5 /src/macro/DZDCG_guts.h | |
| parent | 948ce995e5d8976df0b5a69d7f4c3cb63a3448da (diff) | |
Update of macros for 4th order differencing.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinBase/ADMMacros/trunk@53 b1d164ef-f17a-46e7-89d4-021c7118ef4e
Diffstat (limited to 'src/macro/DZDCG_guts.h')
| -rw-r--r-- | src/macro/DZDCG_guts.h | 25 |
1 files changed, 16 insertions, 9 deletions
diff --git a/src/macro/DZDCG_guts.h b/src/macro/DZDCG_guts.h index 82f3a09..08664d0 100644 --- a/src/macro/DZDCG_guts.h +++ b/src/macro/DZDCG_guts.h @@ -21,15 +21,22 @@ #ifdef FCODE - DZDCG_OO2DZ = 1D0/(2D0*DZDCG_DZ) - - DZDCG_DZDCGXX = DZDCG_OO2DZ*(DZDCG_GXX_KP - DZDCG_GXX_KM) - DZDCG_DZDCGXY = DZDCG_OO2DZ*(DZDCG_GXY_KP - DZDCG_GXY_KM) - DZDCG_DZDCGXZ = DZDCG_OO2DZ*(DZDCG_GXZ_KP - DZDCG_GXZ_KM) - DZDCG_DZDCGYY = DZDCG_OO2DZ*(DZDCG_GYY_KP - DZDCG_GYY_KM) - DZDCG_DZDCGYZ = DZDCG_OO2DZ*(DZDCG_GYZ_KP - DZDCG_GYZ_KM) - DZDCG_DZDCGZZ = DZDCG_OO2DZ*(DZDCG_GZZ_KP - DZDCG_GZZ_KM) - +#include "ADM_Derivative.h" + if (spatial_order.eq.2) then + DZDCG_DZDCGXX = ADM_DZ_2(gxx,i,j,k) + DZDCG_DZDCGXY = ADM_DZ_2(gxy,i,j,k) + DZDCG_DZDCGXZ = ADM_DZ_2(gxz,i,j,k) + DZDCG_DZDCGYY = ADM_DZ_2(gyy,i,j,k) + DZDCG_DZDCGYZ = ADM_DZ_2(gyz,i,j,k) + DZDCG_DZDCGZZ = ADM_DZ_2(gzz,i,j,k) + else + DZDCG_DZDCGXX = ADM_DZ_4(gxx,i,j,k) + DZDCG_DZDCGXY = ADM_DZ_4(gxy,i,j,k) + DZDCG_DZDCGXZ = ADM_DZ_4(gxz,i,j,k) + DZDCG_DZDCGYY = ADM_DZ_4(gyy,i,j,k) + DZDCG_DZDCGYZ = ADM_DZ_4(gyz,i,j,k) + DZDCG_DZDCGZZ = ADM_DZ_4(gzz,i,j,k) + end if #endif #ifdef CCODE |