blob: 6aa344fe93fd5abdf62b4b599c402fa74e961716 (
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 "CactusEinstein/Einstein/src/macro/UPPERMET_guts.h"
#include "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
|