blob: e624a2ccd5da4da6b72e2cc253bda12a69893d89 (
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
|
# Parameter definitions for thorn IDBrillData
shares:einstein
EXTENDS KEYWORD initial_data ""
{
"brilldata" :: "Brill wave initial data"
} ""
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"
REAL 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"
} 1
REAL brill_a "Brill wave: Amplitude"
{
: ::
} 0.0
REAL brill_b "Brill wave: rho^b"
{
: ::
} 2.0
REAL brill_c "Brill wave: (r^2 - r0^2)^(c/2)"
{
: ::
} 2.0
REAL brill_rho0 "Brill wave: radius of torus in rho"
{
: ::
} 0.0
REAL brill_r0 "Brill wave: radius of torus in r"
{
: ::
} 0.0
REAL brill_srho "Brill wave: sigma in rho"
{
: ::
} 1.0
REAL brill_sr "Brill wave: sigma in r"
{
: ::
} 1.0
# 3D Brill wave parameters
REAL brill_d "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: ::
} 0.0
REAL brill_e "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: ::
} 1.0
REAL brill_m "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: ::
} 2.0
REAL brill_n "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: ::
} 2.0
REAL brill_phi0 "3D Brill wave: d rho^m cos^2(n (phi + phi0))"
{
: ::
} 0.0
KEYWORD outputpsi "Output conformal factor?"
{
"yes" :: "Output brillpsi"
"no" :: "Don't output brillpsi"
}"yes"
# Additional parameters
REAL brill_eps "epsilon for finite differencing"
{
0: :: "Positive please"
} 1.e-6
REAL brill_rhofudge "delta rho for axis fudge"
{
0: :: "Positive please"
} 1.e-6
|
