aboutsummaryrefslogtreecommitdiff
path: root/src/macro/MOMYADM_guts.h
blob: 00e246decaa1931b1fba5e3e00b029cc3f7ab4b2 (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
/*@@
  @header   MOMYADM_guts.h
  @date     Aug 98
  @author   Gabrielle Allen
  @desc
  Macro to calculate the spacetime part of the 
  y-Momentum Constraint. That is:

       Del_j K_y^j - Del_y trK 

  @enddesc
@@*/

#ifndef MOMYADM_GUTS
#define MOMYADM_GUTS

#include "../../packages/CactusEinstein/Einstein/src/macro/UPPERMET_guts.h"
#include "../../packages/CactusEinstein/Einstein/src/macro/CDK_guts.h"

#ifdef FCODE 

      MOMYADM_MOMYADM = (CDXCDK_CDXCDKXY - CDYCDK_CDYCDKXX)*UPPERMET_UXX 
     &  + (CDXCDK_CDXCDKYY - CDYCDK_CDYCDKXY)*UPPERMET_UXY + 
     &  (CDXCDK_CDXCDKYZ - 2*CDYCDK_CDYCDKXZ + CDZCDK_CDZCDKXY)*UPPERMET_UXZ + 
     &  (-CDYCDK_CDYCDKYZ + CDZCDK_CDZCDKYY)*UPPERMET_UYZ + 
     &  (-CDYCDK_CDYCDKZZ + CDZCDK_CDZCDKYZ)*UPPERMET_UZZ

#endif

#ifdef CCODE

      MOMYADM_MOMYADM = (CDXCDK_CDXCDKXY - CDYCDK_CDYCDKXX)*UPPERMET_UXX 
       + (CDXCDK_CDXCDKYY - CDYCDK_CDYCDKXY)*UPPERMET_UXY + 
       (CDXCDK_CDXCDKYZ - 2*CDYCDK_CDYCDKXZ + CDZCDK_CDZCDKXY)*UPPERMET_UXZ + 
       (-CDYCDK_CDYCDKYZ + CDZCDK_CDZCDKYY)*UPPERMET_UYZ + 
       (-CDYCDK_CDYCDKZZ + CDZCDK_CDZCDKYZ)*UPPERMET_UZZ;

#endif

#endif