added the Einstien macros.

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

1 files changed, 48 insertions, 0 deletions

diff --git a/src/macro/DXDCG_guts.h b/src/macro/DXDCG_guts.h
new file mode 100644
index 0000000..8647d4f
--- /dev/null
+++ b/src/macro/DXDCG_guts.h
@@ -0,0 +1,48 @@
+/*@@
+  @header   DXDCG_guts.h
+  @date     Jun 98
+  @author   Gabrielle Allen
+  @desc
+  Macro to calculate the first derivatives of the 
+  conformal metric with respect to x
+
+  The macro is defined in terms of standard variables in
+
+  Requires: lower conformal metric at i+1,i-1 ; dx0
+
+  Provides: derivative of lower conformal metric wrt x
+
+  @seefile DXDCG_declare.h
+  @enddesc
+@@*/
+
+#ifndef DXDCG_GUTS
+#define DXDCG_GUTS
+
+#ifdef FCODE 
+
+      DXDCG_OO2DX = 1D0/(2D0*DXDCG_DX)
+    
+      DXDCG_DXDCGXX = DXDCG_OO2DX*(DXDCG_GXX_IP - DXDCG_GXX_IM)
+      DXDCG_DXDCGXY = DXDCG_OO2DX*(DXDCG_GXY_IP - DXDCG_GXY_IM)
+      DXDCG_DXDCGXZ = DXDCG_OO2DX*(DXDCG_GXZ_IP - DXDCG_GXZ_IM)
+      DXDCG_DXDCGYY = DXDCG_OO2DX*(DXDCG_GYY_IP - DXDCG_GYY_IM)
+      DXDCG_DXDCGYZ = DXDCG_OO2DX*(DXDCG_GYZ_IP - DXDCG_GYZ_IM)
+      DXDCG_DXDCGZZ = DXDCG_OO2DX*(DXDCG_GZZ_IP - DXDCG_GZZ_IM)
+
+#endif
+
+#ifdef CCODE
+
+      DXDCG_OO2DX = 1/(2*GH->dx0);
+    
+      DXDCG_DXDCGXX = DXDCG_OO2DX*(DXDCG_GXX_IP - DXDCG_GXX_IM);
+      DXDCG_DXDCGXY = DXDCG_OO2DX*(DXDCG_GXY_IP - DXDCG_GXY_IM);
+      DXDCG_DXDCGXZ = DXDCG_OO2DX*(DXDCG_GXZ_IP - DXDCG_GXZ_IM);
+      DXDCG_DXDCGYY = DXDCG_OO2DX*(DXDCG_GYY_IP - DXDCG_GYY_IM);
+      DXDCG_DXDCGYZ = DXDCG_OO2DX*(DXDCG_GYZ_IP - DXDCG_GYZ_IM);
+      DXDCG_DXDCGZZ = DXDCG_OO2DX*(DXDCG_GZZ_IP - DXDCG_GZZ_IM);
+
+#endif
+
+#endif