src/macro/TRKK_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

/*@@
  @header   TRKK_guts.h
  @date     Aug 98
  @author   Gabrielle Allen
  @desc

  Macro to calculate the trace of the K_lj K^j_i

  @enddesc
@@*/

#ifndef TRKK_GUTS
#define TRKK_GUTS

#include "UPPERMET_guts.h"
#include "KK_guts.h"

#ifdef FCODE

      TRKK_TRKK = UPPERMET_UXX*KK_KKXX + UPPERMET_UYY*KK_KKYY \
           + UPPERMET_UZZ*KK_KKZZ + 2D0*(UPPERMET_UXY*KK_KKXY \
           + UPPERMET_UXZ*KK_KKXZ + UPPERMET_UYZ*KK_KKYZ)
 
#endif

#ifdef CCODE
 
      TRKK_TRKK = UPPERMET_UXX*KK_KKXX + UPPERMET_UYY*KK_KKYY
           + UPPERMET_UZZ*KK_KKZZ + 2.0*(UPPERMET_UXY*KK_KKXY
           + UPPERMET_UXZ*KK_KKXZ + UPPERMET_UYZ*KK_KKYZ);
 
#endif

#endif