added the Einstien macros.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinBase/ADMMacros/trunk@2 b1d164ef-f17a-46e7-89d4-021c7118ef4e

1 files changed, 197 insertions, 0 deletions

diff --git a/src/macro/DDA_declare.h b/src/macro/DDA_declare.h
new file mode 100644
index 0000000..b9d7a29
--- /dev/null
+++ b/src/macro/DDA_declare.h
@@ -0,0 +1,197 @@
+/*@@
+  @header   DDA_declare.h
+  @date     Jul 98
+  @author   Gabrielle Allen
+  @desc
+  Declarations for macro @seefile DDA_guts.h to compute second
+  spatial derivatives of the lapse
+  @enddesc
+@@*/
+
+#ifndef DDA_DECLARE
+#define DDA_DECLARE
+
+#ifdef FCODE
+
+/* Input variables */
+#undef  DDA_A
+#define DDA_A alp(i,j,k)
+#undef  DDA_A_IP 
+#define DDA_A_IP alp(i+1,j,k)
+#undef  DDA_A_IM
+#define DDA_A_IM alp(i-1,j,k)
+#undef  DDA_A_JP 
+#define DDA_A_JP alp(i,j+1,k)
+#undef  DDA_A_JM 
+#define DDA_A_JM alp(i,j-1,k)
+#undef  DDA_A_KP 
+#define DDA_A_KP alp(i,j,k+1)
+#undef  DDA_A_KM
+#define DDA_A_KM alp(i,j,k-1)
+#undef  DDA_A_IPJP 
+#define DDA_A_IPJP alp(i+1,j+1,k)
+#undef  DDA_A_IPJM 
+#define DDA_A_IPJM alp(i+1,j-1,k)
+#undef  DDA_A_IMJP 
+#define DDA_A_IMJP alp(i-1,j+1,k)
+#undef  DDA_A_IMJM 
+#define DDA_A_IMJM alp(i-1,j-1,k)
+#undef  DDA_A_IPKP 
+#define DDA_A_IPKP alp(i+1,j,k+1)
+#undef  DDA_A_IPKM 
+#define DDA_A_IPKM alp(i+1,j,k-1)
+#undef  DDA_A_IMKP 
+#define DDA_A_IMKP alp(i-1,j,k+1)
+#undef  DDA_A_IMKM 
+#define DDA_A_IMKM alp(i-1,j,k-1)
+#undef  DDA_A_JPKP 
+#define DDA_A_JPKP alp(i,j+1,k+1)
+#undef  DDA_A_JPKM 
+#define DDA_A_JPKM alp(i,j+1,k-1)
+#undef  DDA_A_JMKP 
+#define DDA_A_JMKP alp(i,j-1,k+1)
+#undef  DDA_A_JMKM 
+#define DDA_A_JMKM alp(i,j-1,k-1)
+
+/* Output variables */ 
+#undef  DDA_DXXDA
+#define DDA_DXXDA dda_dxxda
+#undef  DDA_DXYDA
+#define DDA_DXYDA dda_dxyda
+#undef  DDA_DXZDA
+#define DDA_DXZDA dda_dxzda
+#undef  DDA_DYYDA
+#define DDA_DYYDA dda_dyyda
+#undef  DDA_DYZDA
+#define DDA_DYZDA dda_dyzda
+#undef  DDA_DZZDA
+#define DDA_DZZDA dda_dzzda
+
+/* Internal variables */
+#undef  DDA_OODX2   
+#define DDA_OODX2 dda_oodx2
+#undef  DDA_OODY2   
+#define DDA_OODY2 dda_oody2
+#undef  DDA_OODZ2   
+#define DDA_OODZ2 dda_oodz2
+#undef  DDA_OO4DXDY
+#define DDA_OO4DXDY dda_oo4dxdy
+#undef  DDA_OO4DXDZ
+#define DDA_OO4DXDZ dda_oo4dxdz
+#undef  DDA_OO4DYDZ
+#define DDA_OO4DYDZ dda_oo4dydz
+#undef  DDA_DX
+#define DDA_DX dx
+#undef  DDA_DY
+#define DDA_DY dy
+#undef  DDA_DZ
+#define DDA_DZ dz
+
+/* Declare internal variables */
+      REAL DDA_OODX2
+      REAL DDA_OODY2
+      REAL DDA_OODZ2
+      REAL DDA_OO4DXDY
+      REAL DDA_OO4DXDZ
+      REAL DDA_OO4DYDZ
+
+/* Declare output variables */
+      REAL DDA_DXXDA
+      REAL DDA_DXYDA
+      REAL DDA_DXZDA
+      REAL DDA_DYYDA
+      REAL DDA_DYZDA
+      REAL DDA_DZZDA
+
+#endif
+
+
+#ifdef CCODE
+
+/* Input variables */
+#undef  DDA_A
+#define DDA_A alp[ijk]
+#undef  DDA_A_IP 
+#define DDA_A_IP alp[ di + ijk]
+#undef  DDA_A_IM
+#define DDA_A_IM alp[-di + ijk]
+#undef  DDA_A_JP 
+#define DDA_A_JP alp[ dj + ijk]
+#undef  DDA_A_JM 
+#define DDA_A_JM alp[-dj + ijk]
+#undef  DDA_A_KP 
+#define DDA_A_KP alp[ dk + ijk]
+#undef  DDA_A_KM
+#define DDA_A_KM alp[-dk + ijk]
+#undef  DDA_A_IPJP 
+#define DDA_A_IPJP alp[ di + dj + ijk]
+#undef  DDA_A_IPJM 
+#define DDA_A_IPJM alp[ di - dj + ijk]
+#undef  DDA_A_IMJP 
+#define DDA_A_IMJP alp[-di + dj + ijk]
+#undef  DDA_A_IMJM 
+#define DDA_A_IMJM alp[-di - dj + ijk]
+#undef  DDA_A_IPKP 
+#define DDA_A_IPKP alp[ di + dk + ijk]
+#undef  DDA_A_IPKM 
+#define DDA_A_IPKM alp[ di - dk + ijk]
+#undef  DDA_A_IMKP 
+#define DDA_A_IMKP alp[-di + dk + ijk]
+#undef  DDA_A_IMKM 
+#define DDA_A_IMKM alp[-di - dk + ijk]
+#undef  DDA_A_JPKP 
+#define DDA_A_JPKP alp[ dj + dk + ijk]
+#undef  DDA_A_JPKM 
+#define DDA_A_JPKM alp[ dj - dk + ijk]
+#undef  DDA_A_JMKP 
+#define DDA_A_JMKP alp[-dj + dk + ijk]
+#undef  DDA_A_JMKM 
+#define DDA_A_JMKM alp[-dj - dk + ijk]
+
+/* Output variables */ 
+#undef  DDA_DXXDA
+#define DDA_DXXDA dda_dxxda
+#undef  DDA_DXYDA
+#define DDA_DXYDA dda_dxyda
+#undef  DDA_DXZDA
+#define DDA_DXZDA dda_dxzda
+#undef  DDA_DYYDA
+#define DDA_DYYDA dda_dyyda
+#undef  DDA_DYZDA
+#define DDA_DYZDA dda_dyzda
+#undef  DDA_DZZDA
+#define DDA_DZZDA dda_dzzda
+
+/* Internal variables */
+#undef  DDA_OODX2   
+#define DDA_OODX2 dda_oodx2
+#undef  DDA_OODY2   
+#define DDA_OODY2 dda_oody2
+#undef  DDA_OODZ2   
+#define DDA_OODZ2 dda_oodz2
+#undef  DDA_OO4DXDY
+#define DDA_OO4DXDY dda_oo4dxdy
+#undef  DDA_OO4DXDZ
+#define DDA_OO4DXDZ dda_oo4dxdz
+#undef  DDA_OO4DYDZ
+#define DDA_OO4DYDZ dda_oo4dydz
+
+/* Declare internal variables */
+      double DDA_OODX2;
+      double DDA_OODY2;
+      double DDA_OODZ2;
+      double DDA_OO4DXDY;
+      double DDA_OO4DXDZ;
+      double DDA_OO4DYDZ;
+
+/* Declare output variables */
+      double DDA_DXXDA;
+      double DDA_DXYDA;
+      double DDA_DXZDA;
+      double DDA_DYYDA;
+      double DDA_DYZDA;
+      double DDA_DZZDA;
+
+#endif
+
+#endif