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/DZDCG_guts.h b/src/macro/DZDCG_guts.h
new file mode 100644
index 0000000..6f4858a
--- /dev/null
+++ b/src/macro/DZDCG_guts.h
@@ -0,0 +1,48 @@
+/*@@
+  @header   DZDCG_guts.h
+  @date     Jun 98
+  @author   Gabrielle Allen
+  @desc
+  Macro to calculate the first derivatives of the 
+  conformal metric with respect to z
+
+  The macro is defined in terms of standard variables in
+
+  Requires: lower conformal metric at k+1,k-1 ; dz0
+
+  Provides: derivative of lower conformal metric wrt z
+
+  @seefile DZDCG_declare.h
+  @enddesc
+@@*/
+
+#ifndef DZDCG_GUTS
+#define DZDCG_GUTS
+
+#ifdef FCODE 
+
+      DZDCG_OO2DZ = 1D0/(2D0*DZDCG_DZ)
+    
+      DZDCG_DZDCGXX = DZDCG_OO2DZ*(DZDCG_GXX_KP - DZDCG_GXX_KM)
+      DZDCG_DZDCGXY = DZDCG_OO2DZ*(DZDCG_GXY_KP - DZDCG_GXY_KM)
+      DZDCG_DZDCGXZ = DZDCG_OO2DZ*(DZDCG_GXZ_KP - DZDCG_GXZ_KM)
+      DZDCG_DZDCGYY = DZDCG_OO2DZ*(DZDCG_GYY_KP - DZDCG_GYY_KM)
+      DZDCG_DZDCGYZ = DZDCG_OO2DZ*(DZDCG_GYZ_KP - DZDCG_GYZ_KM)
+      DZDCG_DZDCGZZ = DZDCG_OO2DZ*(DZDCG_GZZ_KP - DZDCG_GZZ_KM)
+
+#endif
+
+#ifdef CCODE
+
+      DZDCG_OO2DZ = 1/(2*GH->dz0);
+    
+      DZDCG_DZDCGXX = DZDCG_OO2DZ*(DZDCG_GXX_KP - DZDCG_GXX_KM);
+      DZDCG_DZDCGXY = DZDCG_OO2DZ*(DZDCG_GXY_KP - DZDCG_GXY_KM);
+      DZDCG_DZDCGXZ = DZDCG_OO2DZ*(DZDCG_GXZ_KP - DZDCG_GXZ_KM);
+      DZDCG_DZDCGYY = DZDCG_OO2DZ*(DZDCG_GYY_KP - DZDCG_GYY_KM);
+      DZDCG_DZDCGYZ = DZDCG_OO2DZ*(DZDCG_GYZ_KP - DZDCG_GYZ_KM);
+      DZDCG_DZDCGZZ = DZDCG_OO2DZ*(DZDCG_GZZ_KP - DZDCG_GZZ_KM);
+
+#endif
+
+#endif