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
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
|
/*@@
@header DYYDG_guts.h
@date Jun 98
@author Gabrielle Allen
@desc
Macro to calculate the (first and) second derivatives of the
physical metric with respect to y
The macro is defined in terms of standard variables in
@seefile DYYDG_declare.h
The macro uses @seefile DXDG_guts.h and @seefile DXDG_declare.h
@enddesc
@@*/
#ifndef DYYDG_GUTS
#define DYYDG_GUTS
#include "DYDG_guts.h"
#ifdef FCODE
#include "ADM_Derivative.h"
/* Factor involving 2nd derivative of conformal factor */
IF (conformal_state .eq. 0) THEN
DYYDG_FAC = 0
ELSE
DYYDG_FAC = DYDG_PSI4*(4*DYYDG_DYYDPSI_O_PSI \
+ 12*DYDG_DYDPSI_O_PSI*DYDG_DYDPSI_O_PSI)
ENDIF
/* Now calculate the second deriatives */
if (spatial_order.eq.2) then
DYYDG_DYYDGXX = 2*DYDCG_DYDCGXX*DYDG_FAC + DYYDG_FAC*DYDG_GXX \
+ DYDG_PSI4*ADM_DYY_2(gxx,i,j,k)
DYYDG_DYYDGXY = 2*DYDCG_DYDCGXY*DYDG_FAC + DYYDG_FAC*DYDG_GXY \
+ DYDG_PSI4*ADM_DYY_2(gxy,i,j,k)
DYYDG_DYYDGXZ = 2*DYDCG_DYDCGXZ*DYDG_FAC + DYYDG_FAC*DYDG_GXZ \
+ DYDG_PSI4*ADM_DYY_2(gxz,i,j,k)
DYYDG_DYYDGYY = 2*DYDCG_DYDCGYY*DYDG_FAC + DYYDG_FAC*DYDG_GYY \
+ DYDG_PSI4*ADM_DYY_2(gyy,i,j,k)
DYYDG_DYYDGYZ = 2*DYDCG_DYDCGYZ*DYDG_FAC + DYYDG_FAC*DYDG_GYZ \
+ DYDG_PSI4*ADM_DYY_2(gyz,i,j,k)
DYYDG_DYYDGZZ = 2*DYDCG_DYDCGZZ*DYDG_FAC + DYYDG_FAC*DYDG_GZZ \
+ DYDG_PSI4*ADM_DYY_2(gzz,i,j,k)
else
DYYDG_DYYDGXX = 2*DYDCG_DYDCGXX*DYDG_FAC + DYYDG_FAC*DYDG_GXX \
+ DYDG_PSI4*ADM_DYY_4(gxx,i,j,k)
DYYDG_DYYDGXY = 2*DYDCG_DYDCGXY*DYDG_FAC + DYYDG_FAC*DYDG_GXY \
+ DYDG_PSI4*ADM_DYY_4(gxy,i,j,k)
DYYDG_DYYDGXZ = 2*DYDCG_DYDCGXZ*DYDG_FAC + DYYDG_FAC*DYDG_GXZ \
+ DYDG_PSI4*ADM_DYY_4(gxz,i,j,k)
DYYDG_DYYDGYY = 2*DYDCG_DYDCGYY*DYDG_FAC + DYYDG_FAC*DYDG_GYY \
+ DYDG_PSI4*ADM_DYY_4(gyy,i,j,k)
DYYDG_DYYDGYZ = 2*DYDCG_DYDCGYZ*DYDG_FAC + DYYDG_FAC*DYDG_GYZ \
+ DYDG_PSI4*ADM_DYY_4(gyz,i,j,k)
DYYDG_DYYDGZZ = 2*DYDCG_DYDCGZZ*DYDG_FAC + DYYDG_FAC*DYDG_GZZ \
+ DYDG_PSI4*ADM_DYY_4(gzz,i,j,k)
end if
#endif
#ifdef CCODE
/* Factor involving 2nd derivative of conformal factor */
DYYDG_FAC = ((*conformal_state == 0)?0:
DYDG_PSI4*(4*DYYDG_DYYDPSI_O_PSI + 12*DYDG_DYDPSI_O_PSI*DYDG_DYDPSI_O_PSI));
/* Now calculate the second deriatives */
DYYDG_DYYDGXX = 2*DYDCG_DYDCGXX*DYDG_FAC+DYYDG_FAC*DYDG_GXX+DYDG_PSI4
*DYYDG_OODY2*(DYDCG_GXX_JP-2*DYDG_GXX+DYDCG_GXX_JM);
DYYDG_DYYDGXY = 2*DYDCG_DYDCGXY*DYDG_FAC+DYYDG_FAC*DYDG_GXY+DYDG_PSI4
*DYYDG_OODY2*(DYDCG_GXY_JP-2*DYDG_GXY+DYDCG_GXY_JM);
DYYDG_DYYDGXZ = 2*DYDCG_DYDCGXZ*DYDG_FAC+DYYDG_FAC*DYDG_GXZ+DYDG_PSI4
*DYYDG_OODY2*(DYDCG_GXZ_JP-2*DYDG_GXZ+DYDCG_GXZ_JM);
DYYDG_DYYDGYY = 2*DYDCG_DYDCGYY*DYDG_FAC+DYYDG_FAC*DYDG_GYY+DYDG_PSI4
*DYYDG_OODY2*(DYDCG_GYY_JP-2*DYDG_GYY+DYDCG_GYY_JM);
DYYDG_DYYDGYZ = 2*DYDCG_DYDCGYZ*DYDG_FAC+DYYDG_FAC*DYDG_GYZ+DYDG_PSI4
*DYYDG_OODY2*(DYDCG_GYZ_JP-2*DYDG_GYZ+DYDCG_GYZ_JM);
DYYDG_DYYDGZZ = 2*DYDCG_DYDCGZZ*DYDG_FAC+DYYDG_FAC*DYDG_GZZ+DYDG_PSI4
*DYYDG_OODY2*(DYDCG_GZZ_JP-2*DYDG_GZZ+DYDCG_GZZ_JM);
#endif
#endif
|