aboutsummaryrefslogtreecommitdiff
path: root/src/macro/CHR1_guts.h
blob: 834957b521f2cbec79c38c911fb18f0b1bbca68d (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*0.5
      CHR1_XXY = DYDG_DYDGXX*0.5
      CHR1_XXZ = DZDG_DZDGXX*0.5
      CHR1_XYY = -DXDG_DXDGYY*0.5 + DYDG_DYDGXY
      CHR1_XYZ = (-DXDG_DXDGYZ + DYDG_DYDGXZ + DZDG_DZDGXY)*0.5
      CHR1_XZZ = -DXDG_DXDGZZ*0.5 + DZDG_DZDGXZ
      CHR1_YXX = DXDG_DXDGXY - DYDG_DYDGXX*0.5
      CHR1_YXY = DXDG_DXDGYY*0.5
      CHR1_YXZ = ( DXDG_DXDGYZ - DYDG_DYDGXZ + DZDG_DZDGXY)*0.5
      CHR1_YYY = DYDG_DYDGYY*0.5
      CHR1_YYZ = DZDG_DZDGYY*0.5
      CHR1_YZZ = -DYDG_DYDGZZ*0.5 + DZDG_DZDGYZ
      CHR1_ZXX =  DXDG_DXDGXZ - DZDG_DZDGXX*0.5
      CHR1_ZXY = ( DXDG_DXDGYZ + DYDG_DYDGXZ - DZDG_DZDGXY)*0.5
      CHR1_ZXZ = DXDG_DXDGZZ*0.5
      CHR1_ZYY = DYDG_DYDGYZ - DZDG_DZDGYY*0.5
      CHR1_ZYZ = DYDG_DYDGZZ*0.5
      CHR1_ZZZ = DZDG_DZDGZZ*0.5

#endif

#ifdef CCODE

      CHR1_XXX = DXDG_DXDGXX*0.5;
      CHR1_XXY = DYDG_DYDGXX*0.5;
      CHR1_XXZ = DZDG_DZDGXX*0.5;
      CHR1_XYY = -DXDG_DXDGYY*0.5 + DYDG_DYDGXY;
      CHR1_XYZ = (-DXDG_DXDGYZ + DYDG_DYDGXZ + DZDG_DZDGXY)*0.5;
      CHR1_XZZ = -DXDG_DXDGZZ*0.5 + DZDG_DZDGXZ;
      CHR1_YXX = DXDG_DXDGXY - DYDG_DYDGXX*0.5;
      CHR1_YXY = DXDG_DXDGYY*0.5;
      CHR1_YXZ = ( DXDG_DXDGYZ - DYDG_DYDGXZ + DZDG_DZDGXY)*0.5;
      CHR1_YYY = DYDG_DYDGYY*0.5;
      CHR1_YYZ = DZDG_DZDGYY*0.5;
      CHR1_YZZ = -DYDG_DYDGZZ*0.5 + DZDG_DZDGYZ;
      CHR1_ZXX =  DXDG_DXDGXZ - DZDG_DZDGXX*0.5;
      CHR1_ZXY = ( DXDG_DXDGYZ + DYDG_DYDGXZ - DZDG_DZDGXY)*0.5;
      CHR1_ZXZ = DXDG_DXDGZZ*0.5;
      CHR1_ZYY = DYDG_DYDGYZ - DZDG_DZDGYY*0.5;
      CHR1_ZYZ = DYDG_DYDGZZ*0.5;
      CHR1_ZZZ = DZDG_DZDGZZ*0.5;

#endif

#endif