blob: 0a27f17fcf991e4cc0bacaff39991d13db2cdd22 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
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 "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_state .ne. 0) 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_state != 0)
{
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
|
