added the Einstien macros.

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

1 files changed, 80 insertions, 0 deletions

diff --git a/src/macro/LIEG_guts.h b/src/macro/LIEG_guts.h
new file mode 100644
index 0000000..a769079
--- /dev/null
+++ b/src/macro/LIEG_guts.h
@@ -0,0 +1,80 @@
+/*@@
+  @header   LIEG_guts.h
+  @date     Jun 98
+  @author   Gabrielle Allen
+  @desc
+  Macro to calculate the Lie derivative of the lower 
+  conformal metric
+
+  @enddesc
+@@*/
+
+#ifndef LIEG_GUTS
+#define LIEG_GUTS
+
+#include "../../packages/CactusEinstein/Einstein/src/macro/DB_guts.h"
+#include "../../packages/CactusEinstein/Einstein/src/macro/DG_guts.h"
+
+#ifdef FCODE
+
+      LIEG_LGXX = DXDCG_DXDCGXX*LIEG_BX + DYDCG_DYDCGXX*LIEG_BY +
+     & DZDCG_DZDCGXX*LIEG_BZ + 2*(DXDB_DXDBX*LIEG_GXX + DXDB_DXDBY*LIEG_GXY +
+     & DXDB_DXDBZ*LIEG_GXZ)
+
+      LIEG_LGXY = DXDCG_DXDCGXY*LIEG_BX + DYDCG_DYDCGXY*LIEG_BY +
+     & DZDCG_DZDCGXY*LIEG_BZ + DYDB_DYDBX*LIEG_GXX + (DXDB_DXDBX +
+     & DYDB_DYDBY)*LIEG_GXY + DYDB_DYDBZ*LIEG_GXZ + DXDB_DXDBY*LIEG_GYY +
+     & DXDB_DXDBZ*LIEG_GYZ
+
+      LIEG_LGXZ = DXDCG_DXDCGXZ*LIEG_BX + DYDCG_DYDCGXZ*LIEG_BY +
+     & DZDCG_DZDCGXZ*LIEG_BZ + DZDB_DZDBX*LIEG_GXX + DZDB_DZDBY*LIEG_GXY +
+     & (DXDB_DXDBX + DZDB_DZDBZ)*LIEG_GXZ + DXDB_DXDBY*LIEG_GYZ +
+     & DXDB_DXDBZ*LIEG_GZZ
+
+      LIEG_LGYY = DXDCG_DXDCGYY*LIEG_BX +DYDCG_DYDCGYY*LIEG_BY +
+     & DZDCG_DZDCGYY*LIEG_BZ + 2*(DYDB_DYDBX*LIEG_GXY + DYDB_DYDBY*LIEG_GYY +
+     & DYDB_DYDBZ*LIEG_GYZ)
+
+      LIEG_LGYZ = DXDCG_DXDCGYZ*LIEG_BX + DYDCG_DYDCGYZ*LIEG_BY +
+     & DZDCG_DZDCGYZ*LIEG_BZ + DZDB_DZDBX*LIEG_GXY + DYDB_DYDBX*LIEG_GXZ +
+     & DZDB_DZDBY*LIEG_GYY + (DYDB_DYDBY + DZDB_DZDBZ)*LIEG_GYZ +
+     & DYDB_DYDBZ*LIEG_GZZ
+
+      LIEG_LGZZ = DXDCG_DXDCGZZ*LIEG_BX + DYDCG_DYDCGZZ*LIEG_BY +
+     &DZDCG_DZDCGZZ*LIEG_BZ + 2*(DZDB_DZDBX*LIEG_GXZ + DZDB_DZDBY*LIEG_GYZ +
+     &DZDB_DZDBZ*LIEG_GZZ)
+
+#endif
+
+#ifdef CCODE
+
+      LIEG_LGXX = DXDCG_DXDCGXX*LIEG_BX + DYDCG_DYDCGXX*LIEG_BY +
+      DZDCG_DZDCGXX*LIEG_BZ + 2*(DXDB_DXDBX*LIEG_GXX + DXDB_DXDBY*LIEG_GXY +
+      DXDB_DXDBZ*LIEG_GZ) ;
+
+      LIEG_LGXY = DXDCG_DXDCGXY*LIEG_BX + DYDCG_DYDCGXY*LIEG_BY +
+      DZDCG_DZDCGXY*LIEG_BZ + CDYDB_CDYDBX*LIEG_GXX + (DXDB_DXDBX +
+      DYDB_DYDBY)*LIEG_GXY + DYDB_DYDBZ*LIEG_GXZ+ DXDB_DXDBY*LIEG_GYY +
+      DXDB_DXDBZ*LIEG_GYZ ;
+
+      LIEG_LGXZ = DXDCG_DXDCGXZ*LIEG_BX + DYDCG_DYDCGXZ*LIEG_BY +
+      DZDCG_DZDCGXZ*LIEG_BZ + DZDB_DZDBX*LIEG_GXX + DZDB_DZDBY*LIEG_GXY +
+      (DXDB_DXDBX + DZDB_DZDBZ)*LIEG_GXZ+ DXDB_DXDBY*LIEG_GYZ +
+      DXDB_DXDBZ*LIEG_GZZ ;
+
+      LIEG_LGYY = DXDCG_DXDCGYY*LIEG_BX + DYDCG_DYDCGYY*LIEG_BY +
+      DZDCG_DZDCGYY*LIEG_BZ + 2*(DYDB_DYDBX*LIEG_GXY + DYDB_DYDBY*LIEG_GYY +
+      DYDB_DYDBZ*LIEG_GYZ) ;
+
+      LIEG_LGYZ = DXDCG_DXDCGYZ*LIEG_BX + DYDCG_DYDCGYZ*LIEG_BY +
+      DZDCG_DZDCGYZ*LIEG_BZ + DZDB_DZDBX*LIEG_GXY + DYDB_DYDBX*LIEG_GXZ+
+      DZDB_DZDBY*LIEG_GYY + (DYDB_DYDBY + DZDB_DZDBZ)*LIEG_GYZ +
+      DYDB_DYDBZ*LIEG_GZZ ;
+
+      LIEG_LGZZ = DXDCG_DXDCGZZ*LIEG_BX + DYDCG_DYDCGZZ*LIEG_BY +
+      DZDCG_DZDCGZZ*LIEG_BZ + 2*(DZDB_DZDBX*LIEG_GXZ + DZDB_DZDBY*LIEG_GYZ +
+      DZDB_DZDBZ*LIEG_GZZ) ;
+
+#endif
+
+#endif