blob: eaa690893ab6fb52c734e6ed8c824c58c6e873c7 (
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
43
44
45
46
47
48
49
50
51
|
/*@@
@header GAMMA_guts.h
@date Oct 98
@author Gabrielle Allen
@desc
Macro to calculate Gauge Source Functions (Gammas)
That is GAMMA_Ga = g^bc CHR2_abc
@enddesc
@@*/
#ifndef GAMMA_GUTS
#define GAMMA_GUTS
#include "CactusEinstein/Einstein/src/macro/CHR2_guts.h"
#include "CactusEinstein/Einstein/src/macro/UPPERMET_guts.h"
#ifdef FCODE
GAMMA_GAMMAX = UPPERMET_UXX*CHR2_XXX + 2D0*UPPERMET_UXY*CHR2_XXY
& + 2D0*UPPERMET_UXZ*CHR2_XXZ + UPPERMET_UYY*CHR2_XYY
& + 2D0*UPPERMET_UYZ*CHR2_XYZ + UPPERMET_UZZ*CHR2_XZZ
GAMMA_GAMMAY = UPPERMET_UXX*CHR2_YXX + 2D0*UPPERMET_UXY*CHR2_YXY
& + 2D0*UPPERMET_UXZ*CHR2_YXZ + UPPERMET_UYY*CHR2_YYY
& + 2D0*UPPERMET_UYZ*CHR2_YYZ + UPPERMET_UZZ*CHR2_YZZ
GAMMA_GAMMAZ = UPPERMET_UXX*CHR2_ZXX + 2D0*UPPERMET_UXY*CHR2_ZXY
& + 2D0*UPPERMET_UXZ*CHR2_ZXZ + UPPERMET_UYY*CHR2_ZYY
& + 2D0*UPPERMET_UYZ*CHR2_ZYZ + UPPERMET_UZZ*CHR2_ZZZ
#endif
#ifdef CCODE
GAMMA_GAMMAX = UPPERMET_UXX*CHR2_XXX + 2*UPPERMET_UXY*CHR2_XXY
+ 2*UPPERMET_UXZ*CHR2_XXZ + UPPERMET_UYY*CHR2_XYY
+ 2*UPPERMET_UYZ*CHR2_XYZ + UPPERMET_UZZ*CHR2_XZZ;
GAMMA_GAMMAY = UPPERMET_UXX*CHR2_YXX + 2*UPPERMET_UXY*CHR2_YXY
+ 2*UPPERMET_UXZ*CHR2_YXZ + UPPERMET_UYY*CHR2_YYY
+ 2*UPPERMET_UYZ*CHR2_YYZ + UPPERMET_UZZ*CHR2_YZZ;
GAMMA_GAMMAZ = UPPERMET_UXX*CHR2_ZXX + 2*UPPERMET_UXY*CHR2_ZXY
+ 2*UPPERMET_UXZ*CHR2_ZXZ + UPPERMET_UYY*CHR2_ZYY
+ 2*UPPERMET_UYZ*CHR2_ZYZ + UPPERMET_UZZ*CHR2_ZZZ;
#endif
#endif
|
