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
|
/*@@
@header DZZDG_guts.h
@date Jun 98
@author Gabrielle Allen
@desc
Macro to calculate the (first and) second derivatives of the
physical metric with respect to z
The macro is defined in terms of standard variables in
@seefile DZZDG_declare.h
The macro uses @seefile DZDG_guts.h and @seefile DZDG_declare.h
@enddesc
@@*/
#ifndef DZZDG_GUTS
#define DZZDG_GUTS
#include "DZDG_guts.h"
#ifdef FCODE
#include "ADM_Derivative.h"
/* Factor involving 2nd derivative of conformal factor */
IF (conformal_state .eq. 0) THEN
DZZDG_FAC = 0
ELSE
DZZDG_FAC = DZDG_PSI4*(4*DZZDG_DZZDPSI_O_PSI \
+ 12*DZDG_DZDPSI_O_PSI*DZDG_DZDPSI_O_PSI)
ENDIF
if (local_spatial_order.eq.2) then
DZZDG_DZZDGXX = 2*DZDCG_DZDCGXX*DZDG_FAC + DZZDG_FAC*DZDG_GXX \
+ DZDG_PSI4*ADM_DZZ_2(gxx,i,j,k)
DZZDG_DZZDGXY = 2*DZDCG_DZDCGXY*DZDG_FAC + DZZDG_FAC*DZDG_GXY \
+ DZDG_PSI4*ADM_DZZ_2(gxy,i,j,k)
DZZDG_DZZDGXZ = 2*DZDCG_DZDCGXZ*DZDG_FAC + DZZDG_FAC*DZDG_GXZ \
+ DZDG_PSI4*ADM_DZZ_2(gxz,i,j,k)
DZZDG_DZZDGYY = 2*DZDCG_DZDCGYY*DZDG_FAC + DZZDG_FAC*DZDG_GYY \
+ DZDG_PSI4*ADM_DZZ_2(gyy,i,j,k)
DZZDG_DZZDGYZ = 2*DZDCG_DZDCGYZ*DZDG_FAC + DZZDG_FAC*DZDG_GYZ \
+ DZDG_PSI4*ADM_DZZ_2(gyz,i,j,k)
DZZDG_DZZDGZZ = 2*DZDCG_DZDCGZZ*DZDG_FAC + DZZDG_FAC*DZDG_GZZ \
+ DZDG_PSI4*ADM_DZZ_2(gzz,i,j,k)
else
DZZDG_DZZDGXX = 2*DZDCG_DZDCGXX*DZDG_FAC + DZZDG_FAC*DZDG_GXX \
+ DZDG_PSI4*ADM_DZZ_4(gxx,i,j,k)
DZZDG_DZZDGXY = 2*DZDCG_DZDCGXY*DZDG_FAC + DZZDG_FAC*DZDG_GXY \
+ DZDG_PSI4*ADM_DZZ_4(gxy,i,j,k)
DZZDG_DZZDGXZ = 2*DZDCG_DZDCGXZ*DZDG_FAC + DZZDG_FAC*DZDG_GXZ \
+ DZDG_PSI4*ADM_DZZ_4(gxz,i,j,k)
DZZDG_DZZDGYY = 2*DZDCG_DZDCGYY*DZDG_FAC + DZZDG_FAC*DZDG_GYY \
+ DZDG_PSI4*ADM_DZZ_4(gyy,i,j,k)
DZZDG_DZZDGYZ = 2*DZDCG_DZDCGYZ*DZDG_FAC + DZZDG_FAC*DZDG_GYZ \
+ DZDG_PSI4*ADM_DZZ_4(gyz,i,j,k)
DZZDG_DZZDGZZ = 2*DZDCG_DZDCGZZ*DZDG_FAC + DZZDG_FAC*DZDG_GZZ \
+ DZDG_PSI4*ADM_DZZ_4(gzz,i,j,k)
end if
#endif
#ifdef CCODE
/* Factor involving 2nd derivative of conformal factor */
DZZDG_FAC = ((*conformal_state == 0)?0:
DZDG_PSI4*(4*DZZDG_DZZDPSI_O_PSI + 12*DZDG_DZDPSI_O_PSI*DZDG_DZDPSI_O_PSI));
/* Now calculate the second deriatives */
DZZDG_DZZDGXX = 2*DZDCG_DZDCGXX*DZDG_FAC+DZZDG_FAC*DZDG_GXX+DZDG_PSI4
*DZZDG_OODZ2*(DZDCG_GXX_KP-2*DZDG_GXX+DZDCG_GXX_KM);
DZZDG_DZZDGXY = 2*DZDCG_DZDCGXY*DZDG_FAC+DZZDG_FAC*DZDG_GXY+DZDG_PSI4
*DZZDG_OODZ2*(DZDCG_GXY_KP-2*DZDG_GXY+DZDCG_GXY_KM);
DZZDG_DZZDGXZ = 2*DZDCG_DZDCGXZ*DZDG_FAC+DZZDG_FAC*DZDG_GXZ+DZDG_PSI4
*DZZDG_OODZ2*(DZDCG_GXZ_KP-2*DZDG_GXZ+DZDCG_GXZ_KM);
DZZDG_DZZDGYY = 2*DZDCG_DZDCGYY*DZDG_FAC+DZZDG_FAC*DZDG_GYY+DZDG_PSI4
*DZZDG_OODZ2*(DZDCG_GYY_KP-2*DZDG_GYY+DZDCG_GYY_KM);
DZZDG_DZZDGYZ = 2*DZDCG_DZDCGYZ*DZDG_FAC+DZZDG_FAC*DZDG_GYZ+DZDG_PSI4
*DZZDG_OODZ2*(DZDCG_GYZ_KP-2*DZDG_GYZ+DZDCG_GYZ_KM);
DZZDG_DZZDGZZ = 2*DZDCG_DZDCGZZ*DZDG_FAC+DZZDG_FAC*DZDG_GZZ+DZDG_PSI4
*DZZDG_OODZ2*(DZDCG_GZZ_KP-2*DZDG_GZZ+DZDCG_GZZ_KM);
#endif
#endif
|
