added the Einstien macros.

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

1 files changed, 138 insertions, 0 deletions

diff --git a/src/macro/DXDK_declare.h b/src/macro/DXDK_declare.h
new file mode 100644
index 0000000..f31bb26
--- /dev/null
+++ b/src/macro/DXDK_declare.h
@@ -0,0 +1,138 @@
+/*@@
+  @header   DXDK_declare.h
+  @date     Jul 98
+  @author   Gabrielle Allen
+  @desc
+  Declarations for macro @seefile DXDK_guts.h to compute first 
+  derivatives of the extrinsic curvature with respect to x
+  @enddesc
+@@*/
+
+#ifndef DXDK_DECLARE
+#define DXDK_DECLARE
+
+#ifdef FCODE
+
+/* Input variables */
+#undef  DXDK_KXX_IP 
+#define DXDK_KXX_IP hxx(i+1,j,k)
+#undef  DXDK_KXX_IM
+#define DXDK_KXX_IM hxx(i-1,j,k)
+#undef  DXDK_KXY_IP 
+#define DXDK_KXY_IP hxy(i+1,j,k)
+#undef  DXDK_KXY_IM 
+#define DXDK_KXY_IM hxy(i-1,j,k)
+#undef  DXDK_KXZ_IP 
+#define DXDK_KXZ_IP hxz(i+1,j,k)
+#undef  DXDK_KXZ_IM
+#define DXDK_KXZ_IM hxz(i-1,j,k)
+#undef  DXDK_KYY_IP 
+#define DXDK_KYY_IP hyy(i+1,j,k)
+#undef  DXDK_KYY_IM 
+#define DXDK_KYY_IM hyy(i-1,j,k)
+#undef  DXDK_KYZ_IP 
+#define DXDK_KYZ_IP hyz(i+1,j,k)
+#undef  DXDK_KYZ_IM 
+#define DXDK_KYZ_IM hyz(i-1,j,k)
+#undef  DXDK_KZZ_IP 
+#define DXDK_KZZ_IP hzz(i+1,j,k)
+#undef  DXDK_KZZ_IM
+#define DXDK_KZZ_IM hzz(i-1,j,k)
+
+/* Output variables */ 
+#undef  DXDK_DXDKXX
+#define DXDK_DXDKXX  dxdk_dxdkxx
+#undef  DXDK_DXDKXY
+#define DXDK_DXDKXY  dxdk_dxdkxy
+#undef  DXDK_DXDKXZ
+#define DXDK_DXDKXZ  dxdk_dxdkxz
+#undef  DXDK_DXDKYY
+#define DXDK_DXDKYY  dxdk_dxdkyy
+#undef  DXDK_DXDKYZ
+#define DXDK_DXDKYZ  dxdk_dxdkyz
+#undef  DXDK_DXDKZZ
+#define DXDK_DXDKZZ  dxdk_dxdkzz
+
+/* Internal variables */
+#undef  DXDK_DX
+#define DXDK_DX    dx
+#undef  DXDK_OO2DX   
+#define DXDK_OO2DX dxdk_oo2dx
+
+/* Declare internal variables */
+      REAL DXDK_OO2DX
+
+/* Declare output variables */
+      REAL DXDK_DXDKXX
+      REAL DXDK_DXDKXY
+      REAL DXDK_DXDKXZ
+      REAL DXDK_DXDKYY
+      REAL DXDK_DXDKYZ
+      REAL DXDK_DXDKZZ
+
+#endif
+
+
+#ifdef CCODE
+
+/* Input variables */
+#undef  DXDK_KXX_IP 
+#define DXDK_KXX_IP hxx[ di+ijk]
+#undef  DXDK_KXX_IM
+#define DXDK_KXX_IM hxx[-di+ijk]
+#undef  DXDK_KXY_IP 
+#define DXDK_KXY_IP hxy[ di+ijk]
+#undef  DXDK_KXY_IM 
+#define DXDK_KXY_IM hxy[-di+ijk]
+#undef  DXDK_KXZ_IP 
+#define DXDK_KXZ_IP hxz[ di+ijk]
+#undef  DXDK_KXZ_IM
+#define DXDK_KXZ_IM hxz[-di+ijk]
+#undef  DXDK_KYY_IP 
+#define DXDK_KYY_IP hyy[ di+ijk]
+#undef  DXDK_KYY_IM 
+#define DXDK_KYY_IM hyy[-di+ijk]
+#undef  DXDK_KYZ_IP 
+#define DXDK_KYZ_IP hyz[ di+ijk]
+#undef  DXDK_KYZ_IM 
+#define DXDK_KYZ_IM hyz[-di+ijk]
+#undef  DXDK_KZZ_IP 
+#define DXDK_KZZ_IP hzz[ di+ijk]
+#undef  DXDK_KZZ_IM
+#define DXDK_KZZ_IM hzz[-di+ijk]
+
+/* Output variables */ 
+#undef  DXDK_DXDKXX
+#define DXDK_DXDKXX  dxdk_dxdkxx
+#undef  DXDK_DXDKXY
+#define DXDK_DXDKXY  dxdk_dxdkxy
+#undef  DXDK_DXDKXZ
+#define DXDK_DXDKXZ  dxdk_dxdkxz
+#undef  DXDK_DXDKYY
+#define DXDK_DXDKYY  dxdk_dxdkyy
+#undef  DXDK_DXDKYZ
+#define DXDK_DXDKYZ  dxdk_dxdkyz
+#undef  DXDK_DXDKZZ
+#define DXDK_DXDKZZ  dxdk_dxdkzz
+
+/* Internal variables */
+#undef  DXDK_DX
+#define DXDK_DX    dx0
+#undef  DXDK_OO2DX   
+#define DXDK_OO2DX dxdk_oo2dx
+
+/* Declare internal variables */
+double DXDK_DX;
+double DXDK_OO2DX;
+
+/* Declare output variables */
+double DXDK_DXDKXX;
+double DXDK_DXDKXY;
+double DXDK_DXDKXZ;
+double DXDK_DXDKYY;
+double DXDK_DXDKYZ;
+double DXDK_DXDKZZ;
+
+#endif
+
+#endif