aboutsummaryrefslogtreecommitdiff
path: root/src/macro/STRESSENERGY_guts.h
diff options
context:
space:
mode:
Diffstat (limited to 'src/macro/STRESSENERGY_guts.h')
-rw-r--r--src/macro/STRESSENERGY_guts.h40
1 files changed, 40 insertions, 0 deletions
diff --git a/src/macro/STRESSENERGY_guts.h b/src/macro/STRESSENERGY_guts.h
new file mode 100644
index 0000000..cffa42f
--- /dev/null
+++ b/src/macro/STRESSENERGY_guts.h
@@ -0,0 +1,40 @@
+/*@@
+ @header STRESSENERGY_guts.h
+ @date Nov 98
+ @author Gabrielle Allen
+ @desc
+ Macro to calculate the hydro quantities
+ scheme
+ @enddesc
+@@*/
+
+#ifndef STRESSENERGY_GUTS
+#define STRESSENERGY_GUTS
+
+c These geometrical quantities are supplied for CalcTmunu
+
+#include "../../packages/CactusEinstein/Einstein/src/macro/DETG_guts.h"
+#include "../../packages/CactusEinstein/Einstein/src/macro/UPPERMET_guts.h"
+
+ uxx = UPPERMET_UXX; uyy = UPPERMET_UYY; uzz = UPPERMET_UZZ
+ uxy = UPPERMET_UXY; uxz = UPPERMET_UXZ; uyz = UPPERMET_UYZ
+
+ detg = DETG_DETCG
+
+c Zero the stress-energy tensor
+
+ STRESSENERGY_TTT = 0.0D0; STRESSENERGY_TTX = 0.0D0
+ STRESSENERGY_TTY = 0.0D0; STRESSENERGY_TTZ = 0.0D0
+ STRESSENERGY_TXX = 0.0D0; STRESSENERGY_TXY = 0.0D0
+ STRESSENERGY_TXZ = 0.0D0; STRESSENERGY_TYY = 0.0D0
+ STRESSENERGY_TYZ = 0.0D0; STRESSENERGY_TZZ = 0.0D0
+
+
+c Calculate the stress-energy tensor
+#include "./CalcTmunu.inc"
+
+
+c Calculate the trace of the stress-energy tensor
+#include "../../packages/CactusEinstein/Einstein/src/macro/TRT_guts.h"
+
+#endif