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/DDA_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/DDA_guts.h')
| -rw-r--r-- | src/macro/DDA_guts.h | 31 |
1 files changed, 17 insertions, 14 deletions
diff --git a/src/macro/DDA_guts.h b/src/macro/DDA_guts.h index de4c849..7f115f0 100644 --- a/src/macro/DDA_guts.h +++ b/src/macro/DDA_guts.h @@ -12,20 +12,23 @@ #ifdef FCODE - DDA_OODX2 = 1D0/(DDA_DX*DDA_DX) - DDA_OODY2 = 1D0/(DDA_DY*DDA_DY) - DDA_OODZ2 = 1D0/(DDA_DZ*DDA_DZ) - DDA_OO4DXDY = 1D0/(4D0*DDA_DX*DDA_DY) - DDA_OO4DXDZ = 1D0/(4D0*DDA_DX*DDA_DZ) - DDA_OO4DYDZ = 1D0/(4D0*DDA_DY*DDA_DZ) - - DDA_DXXDA = DDA_OODX2*(DDA_A_IP - 2D0*DDA_A + DDA_A_IM) - DDA_DYYDA = DDA_OODY2*(DDA_A_JP - 2D0*DDA_A + DDA_A_JM) - DDA_DZZDA = DDA_OODZ2*(DDA_A_KP - 2D0*DDA_A + DDA_A_KM) - - DDA_DXYDA = DDA_OO4DXDY*(DDA_A_IPJP-DDA_A_IPJM-DDA_A_IMJP+DDA_A_IMJM) - DDA_DXZDA = DDA_OO4DXDZ*(DDA_A_IPKP-DDA_A_IPKM-DDA_A_IMKP+DDA_A_IMKM) - DDA_DYZDA = DDA_OO4DYDZ*(DDA_A_JPKP-DDA_A_JPKM-DDA_A_JMKP+DDA_A_JMKM) +#include "ADM_Derivative.h" + + if (spatial_order.eq.2) then + DDA_DXXDA = ADM_DXX_2(alp,i,j,k) + DDA_DXYDA = ADM_DXY_2(alp,i,j,k) + DDA_DXZDA = ADM_DXZ_2(alp,i,j,k) + DDA_DYYDA = ADM_DYY_2(alp,i,j,k) + DDA_DYZDA = ADM_DYZ_2(alp,i,j,k) + DDA_DZZDA = ADM_DZZ_2(alp,i,j,k) + else + DDA_DXXDA = ADM_DXX_4(alp,i,j,k) + DDA_DXYDA = ADM_DXY_4(alp,i,j,k) + DDA_DXZDA = ADM_DXZ_4(alp,i,j,k) + DDA_DYYDA = ADM_DYY_4(alp,i,j,k) + DDA_DYZDA = ADM_DYZ_4(alp,i,j,k) + DDA_DZZDA = ADM_DZZ_4(alp,i,j,k) + end if #endif |