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
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
|
\documentclass{article}
% Use the Cactus ThornGuide style file
% (Automatically used from Cactus distribution, if you have a
% thorn without the Cactus Flesh download this from the Cactus
% homepage at www.cactuscode.org)
\usepackage{../../../../doc/ThornGuide/cactus}
\begin{document}
\title{StaticConformal}
\author{Tom Goodale et al}
\date{$ $Date$ $}
\maketitle
% Do not delete next line
% START CACTUS THORNGUIDE
\begin{abstract}
Base thorn to provide the variables for the static conformal factor
\end{abstract}
\section{Purpose}
This thorn provides the variables defining a static conformal factor
which is used to transform the physical metric. If this thorn is
active and the {\tt ADMBase::metric\_type} parameter is set to
``{\tt static conformal}'' then the {\tt ADMBase::g...} variables are the
conformal values as opposed to the physical values.
The transformation is
$$ g_{ij}^{\mbox{physical}} = \psi^4 g_{ij}^{\mbox{conformal}} $$
The extrinsic curvature is not transformed.
Memory is provided for the conformal factor {\tt psi}, its first
derivatives {\tt psix}, {\tt psiy}, {\tt psiz}, and its second
derivatives {\tt psixx}, {\tt psixy}, {\tt psixz}, {\tt psiyy}, {\tt
psiyz}, and {\tt psizz} depending on the setting of the {\tt
conformal\_storage} parameter.
Note that the first and second ``derivative'' grid functions have an
additional factor of $1 / \psi$ normalisation since this is the most
common use of the derivative. I.e., the grid functions are
\begin{eqnarray*}
{\tt psi} &=& \psi, \\
{\tt psix} &=& \psi_x/\psi, \qquad \mbox{etc}\\
{\tt psixx} &=& \psi_{ij}/\psi \qquad \mbox{etc}
\end{eqnarray*}
Thorns need to check the value of the grid scalar
{\tt conformal\_state} to determine how many levels of these variables have
actually been calculated before using the conformal factor:
\begin{description}
\item[{\tt conformal\_state=0}]
No conformal factor has been calculated --- thorns may
assume the conformal factor is 1 at all points. (I.e., the metric is physical.)
\item[{\tt conformal\_state=1}]
The conformal factor has been calulated, but no derivatives.
\item[{\tt conformal\_state=2}]
The conformal factor and its first derivatives have been calculated.
\item[{\tt conformal\_state=3}]
The conformal factor and its first and second derivatives have been calculated.
\end{description}
Note that this means that if you only want to know whether {\tt psi} contains
the values for the conformal factor you can check for {\tt conformal\_state > 0}.
\section{Utilities}
{\tt StaticConformal} provides functions to convert between physical and conformal 3-metric values. It is very important to understand that these functions
apply the conversion {\it in place}. That is, if {\tt gxx} contains the conformal metric value, when the routine is exited it will now contain the physical metric value. These functions {\it do not} change the value of {\tt conformal\_state} and should be used with due care. (These functions are for example used
by some analysis thorns who work only with the physical metric, they apply the
transformation on entry to the analysis routine and switch it back on exit).
\begin{description}
\item[Convert from conformal to physical:]
{\tt
\begin{verbatim}
StaticConf_ConfToPhysInPlace (cctk_lsh(1),
cctk_lsh(2),
cctk_lsh(3),
psi,
gxx,
gxy,
gxz,
gyy,
gyz,
gzz)
\end{verbatim}
}
\item[Convert from physical to conformal:]
{\tt
\begin{verbatim}
StaticConf_PhysToConfInPlace (cctk_lsh(1),
cctk_lsh(2),
cctk_lsh(3),
psi,
gxx,
gxy,
gxz,
gyy,
gyz,
gzz)
\end{verbatim}
}
\end{description}
\section{Comments}
The {\tt StaticConformal} thorn itself does not calculate any conformal
factor, but does initialise the {\tt conformal\_state} variable to 0.
Please note, no thorn should use the {\tt conformal\_state} variable
unless the parameter {\tt metric\_type} is {\tt "static conformal"}.
The {\tt conformal\_state} variable is not assigned storage or
initialised by the \texttt{StaticConformal} thorn in any other case.
However, thorns are free to themselves assign storage for {\tt
conformal\_state} and initialise it to zero if {\tt metric\_type} is
``physical''. In this one case it is safe for them to use the
\texttt{conformal\_state} variable if \texttt{metric\_type} is not
``static conformal''. This method allows them to use just one set of
ifs rather than first checking the metric\_type and then the
conformal\_state variable if the metric\_type is ``static conformal''.
% Do not delete next line
% END CACTUS THORNGUIDE
\end{document}
|
