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
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
|
# Parameter definitions for thorn AHFinderDirect
# $Header$
################################################################################
#
# parameters to define the patch system
#
private:
real origin_x "global x coordinate of patch system origin"
{
*:* :: "any real number"
} 0.0
real origin_y "global y coordinate of patch system origin"
{
*:* :: "any real number"
} 0.0
real origin_z "global z coordinate of patch system origin"
{
*:* :: "any real number"
} 0.0
keyword patch_system_type "what type of patch system should we use?"
{
"full sphere" :: "full sphere, no symmetries"
"+z hemisphere" :: "mirror symmetry across z=0 plane"
"+xy quadrant" :: "90 degree periodic rotation symmetry about z axis"
"+xz quadrant" :: "mirror symmetry across z=0 plane *and* \
180 degree periodic rotation symmetry about z axis"
"+xyz octant" :: "mirror symmetry across z=0 plane *and* \
90 degree periodic rotation symmetry about z axis"
} "full sphere"
int N_ghost_points "number of ghost zones on each side of a patch"
{
0:* :: "any integer >= 0"
} 2
int N_overlap_points \
"number of grid points that nominally-just-touching patches should overlap"
{
1:*:2 :: "any integer >= 0; current implementation requires that it be odd"
} 1
real delta_drho_dsigma "angular grid spacing of patches, in degrees"
{
(0.0:* :: "any real number > 0.0"
} 5.0
################################################################################
#
# parameters for the interpolator used to interpolate the "geometry"
# $g_{ij}$ and $K_{ij}$ to the apparent horizon surface position
#
string geometry_interpolator_name \
"name under which the geometry interpolation operator is registered in Cactus"
{
.* :: "any string"
} "generalized polynomial interpolation"
string geometry_interpolator_pars \
"parameters for the geometry interpolator"
{
.* :: "any string acceptable to Util_TableSetFromString()"
} ""
################################################################################
#
# parameters for the interpatch interpolator
#
string interpatch_interpolator_name \
"name under which the interpatch interpolation operator is registered in Cactus"
{
.* :: "any string"
} "generalized polynomial interpolation"
string interpatch_interpolator_pars \
"parameters for the interpatch interpolator"
{
.* :: "any string acceptable to Util_TableSetFromString()"
} ""
################################################################################
#
# parameters specifying the initial guess for the apparent horizon shape
#
keyword initial_guess_method \
"method used to set up initial guess for apparent horizon shape"
{
"read from file":: "read from input file"
"ellipsoid" :: "ellipsoid specified by initial_guess__ellipsoid__* parameters"
} "read from file"
# parameters for initial_guess_method = "read from file"
string initial_guess__read_from_file__file_name \
"file name from which to read initial guess for apparent horizon shape"
{
.* :: "any string"
} "h.dat"
# parameters for initial_guess_method = "ellipsoid"
real initial_guess__ellipsoid__center_global_x \
"global x coordinate of ellipsoid center"
{
*:* :: "any real number"
} 0.0
real initial_guess__ellipsoid__center_global_y \
"global y coordinate of ellipsoid center"
{
*:* :: "any real number"
} 0.0
real initial_guess__ellipsoid__center_global_z \
"global z coordinate of ellipsoid center"
{
*:* :: "any real number"
} 0.0
real initial_guess__ellipsoid__radius_x "x coordinate radius of ellipsoid"
{
*:* :: "any real number"
} 2.0
real initial_guess__ellipsoid__radius_y "y coordinate radius of ellipsoid"
{
*:* :: "any real number"
} 2.0
real initial_guess__ellipsoid__radius_z "z coordinate radius of ellipsoid"
{
*:* :: "any real number"
} 2.0
################################################################################
#
# overall parameters for the apparent horizon finder
#
keyword method "top-level method used to find the apparent horizon"
{
"horizon" :: \
"don't find the apparent horizon, just evaluate LHS function H(h)"
} "horizon"
################################################################################
#
# parameters for the test driver src/util/test_patch_system.cc
#
keyword which_test "which test should we do?"
{
"gridfn" :: "set up test fn(x,y,z), print it"
"read gridfn" :: "read in ghosted test fn(x,y,z), print it"
"synchronize" :: "set up test fn(x,y,z), synchronize it, print errors"
"synchronize Jacobian":: \
"set up test fn(x,y,z), compute Jacobian of synchronize(), compare with NP"
"derivatives" :: "set up test fn(rho,sigma), take derivs, print errors"
} "gridfn"
int which_derivs "bit flags to specify which derivatives to test"
{
0:63 :: "any set of bit flags"
} 63
#
# Don't set this too small or roundoff errors will become large.
# Don't set this too large or finite differencing errors will become large.
# In practice the default value should be fine.
#
real NP_Jacobian__perturbation_amplitude \
"NP Jacobian: numerical-perturbation amplitude"
{
(0.0:* :: "any real number > 0"
} 1.0e-6
string Jacobian_file_name "name of the Jacobian output file"
{
.+ :: "any valid file name"
} "Jacobian.dat"
|
