blob: aa565d6ccb56cf4edcf24733486be2b4d63604c4 (
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
52
53
54
55
56
57
58
59
60
61
62
63
64
65
|
/*@@
@header CHR1_guts.h
@date Jun 98
@author Gabrielle Allen
@desc
Macro to calculate Christoffel symbols of the first kind.
That is CHR1_cab = 1/2 (g_ac,b + gbc,a - gab,c)
@enddesc
@@*/
#ifndef CHR1_GUTS
#define CHR1_GUTS
#include "DG_guts.h"
#ifdef FCODE
CHR1_XXX = DXDG_DXDGXX/2D0
CHR1_XXY = DYDG_DYDGXX/2D0
CHR1_XXZ = DZDG_DZDGXX/2D0
CHR1_XYY = -DXDG_DXDGYY/2D0 + DYDG_DYDGXY
CHR1_XYZ = (-DXDG_DXDGYZ + DYDG_DYDGXZ + DZDG_DZDGXY)/2D0
CHR1_XZZ = -DXDG_DXDGZZ/2D0 + DZDG_DZDGXZ
CHR1_YXX = DXDG_DXDGXY - DYDG_DYDGXX/2D0
CHR1_YXY = DXDG_DXDGYY/2D0
CHR1_YXZ = ( DXDG_DXDGYZ - DYDG_DYDGXZ + DZDG_DZDGXY)/2D0
CHR1_YYY = DYDG_DYDGYY/2D0
CHR1_YYZ = DZDG_DZDGYY/2D0
CHR1_YZZ = -DYDG_DYDGZZ/2D0 + DZDG_DZDGYZ
CHR1_ZXX = DXDG_DXDGXZ - DZDG_DZDGXX/2D0
CHR1_ZXY = ( DXDG_DXDGYZ + DYDG_DYDGXZ - DZDG_DZDGXY)/2D0
CHR1_ZXZ = DXDG_DXDGZZ/2D0
CHR1_ZYY = DYDG_DYDGYZ - DZDG_DZDGYY/2D0
CHR1_ZYZ = DYDG_DYDGZZ/2D0
CHR1_ZZZ = DZDG_DZDGZZ/2D0
#endif
#ifdef CCODE
CHR1_XXX = DXDG_DXDGXX/2;
CHR1_XXY = DYDG_DYDGXX/2;
CHR1_XXZ = DZDG_DZDGXX/2;
CHR1_XYY = -DXDG_DXDGYY/2 + DYDG_DYDGXY;
CHR1_XYZ = (-DXDG_DXDGYZ + DYDG_DYDGXZ + DZDG_DZDGXY)/2;
CHR1_XZZ = -DXDG_DXDGZZ/2 + DZDG_DZDGXZ;
CHR1_YXX = DXDG_DXDGXY - DYDG_DYDGXX/2;
CHR1_YXY = DXDG_DXDGYY/2;
CHR1_YXZ = ( DXDG_DXDGYZ - DYDG_DYDGXZ + DZDG_DZDGXY)/2;
CHR1_YYY = DYDG_DYDGYY/2;
CHR1_YYZ = DZDG_DZDGYY/2;
CHR1_YZZ = -DYDG_DYDGZZ/2 + DZDG_DZDGYZ;
CHR1_ZXX = DXDG_DXDGXZ - DZDG_DZDGXX/2;
CHR1_ZXY = ( DXDG_DXDGYZ + DYDG_DYDGXZ - DZDG_DZDGXY)/2;
CHR1_ZXZ = DXDG_DXDGZZ/2;
CHR1_ZYY = DYDG_DYDGYZ - DZDG_DZDGYY/2;
CHR1_ZYZ = DYDG_DYDGZZ/2;
CHR1_ZZZ = DZDG_DZDGZZ/2;
#endif
#endif
|
