blob: 9ec04f9c9f2a1a9e73a52a10f59d9543876fee63 (
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
|
/*@@
@header TRK_guts.h
@date Jun 98
@author Gabrielle Allen
@desc
Macro to calculate the trace of the extrinsic curvature
The macro is defined in terms of standard variables in
TRK_declare.h
Requires: Upper physical metric
Lower extrainsic curvature
Provides: Trace of the extrinsic curvature
@enddesc
@@*/
#ifndef TRK_GUTS
#define TRK_GUTS
#include "CactusEinstein/Einstein/src/macro/UPPERMET_guts.h"
#ifdef FCODE
TRK_TRK = UPPERMET_UXX*TRK_HXX + UPPERMET_UYY*TRK_HYY
& + UPPERMET_UZZ*TRK_HZZ + 2D0*(UPPERMET_UXY*TRK_HXY
& + UPPERMET_UXZ*TRK_HXZ + UPPERMET_UYZ*TRK_HYZ)
#endif
#ifdef CCODE
TRK_TRK = UPPERMET_UXX*TRK_HXX + UPPERMET_UYY*TRK_HYY
+ UPPERMET_UZZ*TRK_HZZ + 2*(UPPERMET_UXY*TRK_HXY
+ UPPERMET_UXZ*TRK_HXZ + UPPERMET_UYZ*TRK_HYZ);
#endif
#endif
|