blob: c7d5b27e6b495aa5baeede5a1f7f7d1dc99d34ab (
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
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
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
|
# Parameter definitions for thorn IDBrillData
# $Header$
shares: ADMBase
EXTENDS KEYWORD initial_data
{
"brilldata" :: "Brill wave initial data"
"brilldata2D" :: "Brill wave initial data assuming axisymmetry"
}
USES KEYWORD metric_type
shares: StaticConformal
USES KEYWORD conformal_storage
private:
# Parameters for elliptic solve
KEYWORD solver "Which elliptic solver to use"
{
"sor" :: "Use SOR solver"
"petsc" :: "Use PETSc solver"
"bam" :: "Use bam solver"
} "sor"
KEYWORD bound "Which boundary condition to use"
{
"const" :: "constant boundary: set const_v0"
"robin" :: "Robin boundary: set robin_falloff, robin_inf"
} "robin"
INT robin_falloff "Fall-off of Robin BC"
{
0: :: "any positive integer value"
} 1
REAL const_v0 "Value of constant BC"
{
: :: "anything goes"
} 1.0
REAL robin_inf "Value at infinity of Robin BC"
{
: :: "anything goes"
} 1.0
REAL thresh "How far (absolute norm) to go"
{
0.0: :: "Positive number please"
} 0.00001
# Brill wave parameters
KEYWORD q_function "Form of function q [0,1,2]"
{
"exp" :: "contains e^{-z^2} factor"
"eppley" :: "contains 1/(1+r^2) factor"
"gundlach" :: "contains e^{-r^2} factor"
} "gundlach"
# -------------------------------------------------------------
REAL exp_a "Exp Brill wave: Amplitude"
{
: :: "Anything"
} 0.0
INT exp_b "Exp Brill wave: used in exponent in rho: rho^(2+b)"
{
: :: "Anything"
} 2
REAL exp_rho0 "Exp Brill wave: radius of torus in rho"
{
0:* :: "Positive"
} 0.0
REAL exp_sigmaz "Exp Brill wave: sigma in z"
{
(0:* :: "Positive"
} 1.0
# -------------------------------------------------------------
REAL eppley_a "Eppley Brill wave: Amplitude"
{
: :: "Anything"
} 0.0
INT eppley_b "Eppley Brill wave: used in exponent in rho: rho^b"
{
: :: "Anything"
} 2
INT eppley_c "Eppley Brill wave: (r^2 - r0^2)^(c/2)"
{
: :: "Anything"
} 2
REAL eppley_r0 "Eppley Brill wave: radius of torus in r"
{
0:* :: "Positive"
} 0.0
REAL eppley_sigmarho "Eppley Brill wave: sigma in rho"
{
: :: "Anything"
} 1.0
REAL eppley_sigmar "Eppley Brill wave: sigma in r"
{
: :: "Anything"
} 1.0
# -------------------------------------------------------------
REAL gundlach_a "Gundlach Brill wave: Amplitude"
{
: :: "Anything"
} 0.0
INT gundlach_b "Gundlach Brill wave: used in exponent in rho: rho^b"
{
: :: "Anything"
} 2
REAL gundlach_sigmarho "Gundlach Brill wave: sigma in rho"
{
: :: "Anything"
} 1.0
REAL gundlach_r0 "Gundlach Brill wave: radius of torus in r"
{
0:* :: "Positive"
} 0.0
INT gundlach_c "Gundlach Brill wave: (r^2 - r0^2)^(c/2)"
{
: :: "Anything"
} 2
REAL gundlach_sigmar "Gundlach Brill wave: sigma in r"
{
: :: "Anything"
} 1.0
# -------------------------------------------------------------
# 3D Brill wave parameters
REAL brill3d_d "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: :: "Anything"
} 0.0
REAL brill3d_e "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: :: "Anything"
} 1.0
REAL brill3d_m "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: :: "Anything"
} 2.0
REAL brill3d_n "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: :: "Anything"
} 2.0
REAL brill3d_phi0 "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: :: "Anything"
} 0.0
# Additional parameters
REAL rhofudge "delta rho for axis fudge"
{
0: :: "Positive please"
} 0.00001
INT sor_maxit "Maximum number of iterations"
{
0:* :: "Positive"
} 10000
BOOLEAN output_coeffs "output coefficients for elliptic solve"
{
} "no"
|
