diff options
Diffstat (limited to 'src/macro/STRESSENERGY_guts.h')
-rw-r--r-- | src/macro/STRESSENERGY_guts.h | 40 |
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 |