aboutsummaryrefslogtreecommitdiff
path: root/src/macro/LIEG_declare.h
diff options
context:
space:
mode:
Diffstat (limited to 'src/macro/LIEG_declare.h')
-rw-r--r--src/macro/LIEG_declare.h110
1 files changed, 110 insertions, 0 deletions
diff --git a/src/macro/LIEG_declare.h b/src/macro/LIEG_declare.h
new file mode 100644
index 0000000..3f52f68
--- /dev/null
+++ b/src/macro/LIEG_declare.h
@@ -0,0 +1,110 @@
+/*@@
+ @header LIEG_declare.h
+ @date Jun 98
+ @author Gabrielle Allen
+ @desc
+ Macro to calculate the Lie derivative of the lower
+ physical metric
+
+ @enddesc
+@@*/
+
+#ifndef LIEG_DECLARE
+#define LIEG_DECLARE
+
+#include "../../packages/CactusEinstein/Einstein/src/macro/DB_declare.h"
+#include "../../packages/CactusEinstein/Einstein/src/macro/DG_declare.h"
+
+#ifdef FCODE
+
+/* Input variables */
+#undef LIEG_BX
+#define LIEG_BX betax(i,j,k)
+#undef LIEG_BY
+#define LIEG_BY betay(i,j,k)
+#undef LIEG_BZ
+#define LIEG_BZ betaz(i,j,k)
+#undef LIEG_GXX
+#define LIEG_GXX gxx(i,j,k)
+#undef LIEG_GXY
+#define LIEG_GXY gxy(i,j,k)
+#undef LIEG_GXZ
+#define LIEG_GXZ gxz(i,j,k)
+#undef LIEG_GYY
+#define LIEG_GYY gyy(i,j,k)
+#undef LIEG_GYZ
+#define LIEG_GYZ gyz(i,j,k)
+#undef LIEG_GZZ
+#define LIEG_GZZ gzz(i,j,k)
+
+/* Output variables */
+#undef LIEG_LGXX
+#define LIEG_LGXX lieg_lgxx
+#undef LIEG_LGXY
+#define LIEG_LGXY lieg_lgxy
+#undef LIEG_LGXZ
+#define LIEG_LGXZ lieg_lgxz
+#undef LIEG_LGYY
+#define LIEG_LGYY lieg_lgyy
+#undef LIEG_LGYZ
+#define LIEG_LGYZ lieg_lgyz
+#undef LIEG_LGZZ
+#define LIEG_LGZZ lieg_lgzz
+
+/* Declare output variables */
+ REAL LIEG_LGXX
+ REAL LIEG_LGXY
+ REAL LIEG_LGXZ
+ REAL LIEG_LGYY
+ REAL LIEG_LGYZ
+ REAL LIEG_LGZZ
+
+#endif
+
+#ifdef CCODE
+
+/* Input variables */
+#undef LIEG_BX
+#define LIEG_BX betax[ijk]
+#undef LIEG_BY
+#define LIEG_BY betay[ijk]
+#undef LIEG_BZ
+#define LIEG_BZ betaz[ijk]
+#undef LIEG_GXX
+#define LIEG_GXX gxx[ijk]
+#undef LIEG_GXY
+#define LIEG_GXY gxy[ijk]
+#undef LIEG_GXZ
+#define LIEG_GXZ gxz[ijk]
+#undef LIEG_GYY
+#define LIEG_GYY gyy[ijk]
+#undef LIEG_GYZ
+#define LIEG_GYZ gyz[ijk]
+#undef LIEG_GZZ
+#define LIEG_GZZ gzz[ijk]
+
+/* Output variables */
+#undef LIEG_LGXX
+#define LIEG_LGXX lieg_lgxx
+#undef LIEG_LGXY
+#define LIEG_LGXY lieg_lgxy
+#undef LIEG_LGXZ
+#define LIEG_LGXZ lieg_lgxz
+#undef LIEG_LGYY
+#define LIEG_LGYY lieg_lgyy
+#undef LIEG_LGYZ
+#define LIEG_LGYZ lieg_lgyz
+#undef LIEG_LGZZ
+#define LIEG_LGZZ lieg_lgzz
+
+/* Declare output variables */
+ double LIEG_LGXX;
+ double LIEG_LGXY;
+ double LIEG_LGXZ;
+ double LIEG_LGYY;
+ double LIEG_LGYZ;
+ double LIEG_LGZZ;
+
+#endif
+
+#endif