blob: da94c88ed0cc5d5bb6abc61d3425d5db971f5b0c (
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
|
# Parameter definitions for thorn IDBrillData
# $Header$
shares: ADMBase
EXTENDS KEYWORD initial_data
{
"brilldata" :: "Brill wave initial data"
}
USES KEYWORD metric_type
shares: StaticConformal
USES KEYWORD conformal_storage
private:
# Parameters for elliptic solve
KEYWORD brill_solver "Which elliptic solver to use"
{
"sor" :: "Use SOR solver"
"petsc" :: "Use PETSc solver"
"bam" :: "Use bam solver"
} "sor"
KEYWORD brill_bound "Which boundary condition to use"
{
"const" :: "constant boundary: set brill_const_v0"
"robin" :: "Robin boundary: set brill_robin_falloff, brill_robin_inf"
} "const"
INT brill_robin_falloff "Fall-off of Robin BC"
{
0: :: "any positive integer value"
} 1
REAL brill_const_v0 "Value of constant BC"
{
: :: "anything goes"
} 1.0
REAL brill_robin_inf "Value at infinity of Robin BC"
{
: :: "anything goes"
} 1.0
REAL brill_thresh "How far (absolute norm) to go"
{
0.0: :: "Positive number please"
} 0.00001
# Brill wave parameters
BOOLEAN axisym "Axisymmetric initial data?"
{
} "yes"
INT brill_q "Form of function q [0,1,2]"
{
0:2 :: "Only cases 0,1,2 defined at the moment"
} 2
REAL brill_a "Brill wave: Amplitude"
{
: :: "Anything"
} 0.0
INT brill_b "Brill wave: rho^b"
{
: :: "Anything"
} 2
INT brill_c "Brill wave: (r^2 - r0^2)^(c/2)"
{
: :: "Anything"
} 2
REAL brill_rho0 "Brill wave: radius of torus in rho"
{
: :: "Anything"
} 0.0
REAL brill_r0 "Brill wave: radius of torus in r"
{
: :: "Anything"
} 0.0
REAL brill_srho "Brill wave: sigma in rho"
{
: :: "Anything"
} 1.0
REAL brill_sr "Brill wave: sigma in r"
{
: :: "Anything"
} 1.0
REAL brill_sz "Brill wave: sigma in z"
{
: :: "Anything"
} 1.0
# 3D Brill wave parameters
REAL brill_d "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: :: "Anything"
} 0.0
REAL brill_e "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: :: "Anything"
} 1.0
REAL brill_m "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: :: "Anything"
} 2.0
REAL brill_n "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: :: "Anything"
} 2.0
REAL brill_phi0 "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: :: "Anything"
} 0.0
# Additional parameters
REAL brill_rhofudge "delta rho for axis fudge"
{
0: :: "Positive please"
} 0.00001
INT sor_maxit "Maximum number of iterations"
{
0:* :: "Positive"
} 100
|
