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/DZDK_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/DZDK_guts.h')
-rw-r--r-- | src/macro/DZDK_guts.h | 26 |
1 files changed, 17 insertions, 9 deletions
diff --git a/src/macro/DZDK_guts.h b/src/macro/DZDK_guts.h index 1e205ab..cc8e5e6 100644 --- a/src/macro/DZDK_guts.h +++ b/src/macro/DZDK_guts.h @@ -13,15 +13,23 @@ #ifdef FCODE - DZDK_OO2DZ = 1D0/(2D0*DZDK_DZ) - - DZDK_DZDKXX = DZDK_OO2DZ*(DZDK_KXX_KP - DZDK_KXX_KM) - DZDK_DZDKXY = DZDK_OO2DZ*(DZDK_KXY_KP - DZDK_KXY_KM) - DZDK_DZDKXZ = DZDK_OO2DZ*(DZDK_KXZ_KP - DZDK_KXZ_KM) - DZDK_DZDKYY = DZDK_OO2DZ*(DZDK_KYY_KP - DZDK_KYY_KM) - DZDK_DZDKYZ = DZDK_OO2DZ*(DZDK_KYZ_KP - DZDK_KYZ_KM) - DZDK_DZDKZZ = DZDK_OO2DZ*(DZDK_KZZ_KP - DZDK_KZZ_KM) - +#include "ADM_Derivative.h" + + if (spatial_order.eq.2) then + DZDK_DZDKXX = ADM_DZ_2(kxx,i,j,k) + DZDK_DZDKXY = ADM_DZ_2(kxy,i,j,k) + DZDK_DZDKXZ = ADM_DZ_2(kxz,i,j,k) + DZDK_DZDKYY = ADM_DZ_2(kyy,i,j,k) + DZDK_DZDKYZ = ADM_DZ_2(kyz,i,j,k) + DZDK_DZDKZZ = ADM_DZ_2(kzz,i,j,k) + else + DZDK_DZDKXX = ADM_DZ_4(kxx,i,j,k) + DZDK_DZDKXY = ADM_DZ_4(kxy,i,j,k) + DZDK_DZDKXZ = ADM_DZ_4(kxz,i,j,k) + DZDK_DZDKYY = ADM_DZ_4(kyy,i,j,k) + DZDK_DZDKYZ = ADM_DZ_4(kyz,i,j,k) + DZDK_DZDKZZ = ADM_DZ_4(kzz,i,j,k) + end if #endif #ifdef CCODE |