diff options
author | evans <evans@b1d164ef-f17a-46e7-89d4-021c7118ef4e> | 1999-03-09 21:04:02 +0000 |
---|---|---|
committer | evans <evans@b1d164ef-f17a-46e7-89d4-021c7118ef4e> | 1999-03-09 21:04:02 +0000 |
commit | 0a13cbb3d0f5ebf208d962d374a197c66bd13b1b (patch) | |
tree | 19582fa856a1193d17170ad971353f58e7be36c9 /src/macro/TRT_guts.h | |
parent | c58897b8676212b9aee16d6d2f80ca2adeb50c6c (diff) |
added the Einstien macros.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinBase/ADMMacros/trunk@2 b1d164ef-f17a-46e7-89d4-021c7118ef4e
Diffstat (limited to 'src/macro/TRT_guts.h')
-rw-r--r-- | src/macro/TRT_guts.h | 62 |
1 files changed, 62 insertions, 0 deletions
diff --git a/src/macro/TRT_guts.h b/src/macro/TRT_guts.h new file mode 100644 index 0000000..efda382 --- /dev/null +++ b/src/macro/TRT_guts.h @@ -0,0 +1,62 @@ +/*@@ + @header TRT_guts.h + @date Jun 98 + @author Gabrielle Allen + @desc + + Macro to calculate the trace of the (4-)stress energy tensor + + @enddesc +@@*/ + +#ifndef TRT_GUTS +#define TRT_GUTS + +#include "../../packages/CactusEinstein/Einstein/src/macro/UPPERMET_guts.h" + +#ifdef FCODE + + TRT_IALP2 = 1D0/TRT_ALP**2 + + TRT_TRT = -TRT_TTT*TRT_IALP2+TRT_TXX*UPPERMET_UXX+ + & TRT_TYY*UPPERMET_UYY+TRT_TZZ*UPPERMET_UZZ+2D0*( + & TRT_TXY*UPPERMET_UXY+TRT_TXZ*UPPERMET_UXZ+TRT_TYZ* + & UPPERMET_UYZ) + + + if (shift_active /= SHIFT_INACTIVE) then + + TRT_TRT = TRT_TRT - TRT_IALP2*(TRT_TXX*TRT_BX*TRT_BX + & + TRT_TYY*TRT_BY*TRT_BY + TRT_TZZ*TRT_BZ*TRT_BZ + & + 2D0*(TRT_TXY*TRT_BX*TRT_BY + TRT_TXZ*TRT_BX*TRT_BZ + & + TRT_TYZ*TRT_BY*TRT_BZ -TRT_TTX*TRT_BX -TRT_TTY*TRT_BY + & -TRT_TTZ*TRT_BZ)) + + endif + +#endif + +#ifdef CCODE + + TRT_IALP2 = 1D0/TRT_ALP**2; + + TRT_TRT = -TRT_TTT*TRT_IALP2+TRT_TXX*UPPERMET_UXX+ + TRT_TYY*UPPERMET_UYY+TRT_TZZ*UPPERMET_UZZ+2D0*( + TRT_TXY*UPPERMET_UXY+TRT_TXZ*UPPERMET_UXZ+TRT_TYZ* + UPPERMET_UYZ) + + + if (shift) + { + TRT_TRT = TRT_TRT - TRT_IALP2*(TRT_TXX*TRT_BX*TRT_BX + + TRT_TYY*TRT_BY*TRT_BY + TRT_TZZ*TRT_BZ*TRT_BZ + + 2D0*(TRT_TXY*TRT_BX*TRT_BY + TRT_TXZ*TRT_BX*TRT_BZ + + TRT_TYZ*TRT_BY*TRT_BZ -TRT_TTX*TRT_BX -TRT_TTY*TRT_BY + -TRT_TTZ*TRT_BZ)) + + } + +#endif + +#endif + |