diff options
Diffstat (limited to 'src/macro/DDA_guts.h')
-rw-r--r-- | src/macro/DDA_guts.h | 52 |
1 files changed, 52 insertions, 0 deletions
diff --git a/src/macro/DDA_guts.h b/src/macro/DDA_guts.h new file mode 100644 index 0000000..de4c849 --- /dev/null +++ b/src/macro/DDA_guts.h @@ -0,0 +1,52 @@ +/*@@ + @header DDA_guts.h + @date Jul 98 + @author Gabrielle Allen + @desc + Macro to calculate all second spatial derivative of lapse + @enddesc +@@*/ + +#ifndef DDA_GUTS +#define DDA_GUTS + +#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) + +#endif + +#ifdef CCODE + + DDA_OODX2 = 1/(DDA_DX*DDA_DX); + DDA_OODY2 = 1/(DDA_DY*DDA_DY); + DDA_OODZ2 = 1/(DDA_DZ*DDA_DZ); + DDA_OO4DXDY = 1/(4*DDA_DX*DDA_DY); + DDA_OO4DXDZ = 1/(4*DDA_DX*DDA_DZ); + DDA_OO4DYDZ = 1/(4*DDA_DY*DDA_DZ); + + DDA_DXXDA = DDA_OODX2*(DDA_A_IP - 2*DDA_A + DDA_A_IM); + DDA_DYYDA = DDA_OODY2*(DDA_A_JP - 2*DDA_A + DDA_A_JM); + DDA_DZZDA = DDA_OODZ2*(DDA_A_KP - 2*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); + +#endif + +#endif + |