src/macro/TRK_guts.h


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