Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Add EinsteinExact version of Exact's KS-tilted test
The output is sufficiently similar to that of Exact (which uses finite differencing to compute K_ij, rather than analytic expressions as in EinsteinExact) that I believe the test is working, but I have not performed a convergence test in Exact's "eps" parameter.
- Loading branch information
Showing
236 changed files
with
21,323 additions
and
0 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,113 @@ | ||
Cactus::cctk_run_title = "Tilted Kerr/Kerr-Schild data" | ||
Cactus::cctk_full_warnings = yes | ||
|
||
Cactus::cctk_itlast = 2 | ||
|
||
ActiveThorns = "IOUtil" | ||
IO::out_dir = $parfile | ||
IO::out_fileinfo = "axis labels" | ||
IO::parfile_write = "no" | ||
|
||
ActiveThorns = "AEILocalInterp LocalInterp LocalReduce" | ||
|
||
ActiveThorns = "PUGH PUGHInterp PUGHReduce PUGHSlab" | ||
driver::global_nsize = 49 | ||
driver::ghost_size = 2 | ||
|
||
ActiveThorns = "NaNChecker" | ||
|
||
ActiveThorns = "Boundary CartGrid3D CoordBase SymBase SymBase SphericalSurface CoordBase" | ||
CartGrid3D::type = "byspacing" | ||
CartGrid3D::domain = "full" | ||
CartGrid3D::dxyz = 0.125 | ||
CartGrid3D::avoid_origin = no | ||
CoordBase::boundary_size_x_upper = 2 | ||
CoordBase::boundary_size_y_upper = 2 | ||
CoordBase::boundary_size_z_upper = 2 | ||
CoordBase::boundary_size_x_lower = 2 | ||
CoordBase::boundary_size_y_lower = 2 | ||
CoordBase::boundary_size_z_lower = 2 | ||
|
||
ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase" | ||
SpaceMask::use_mask = yes | ||
|
||
ActiveThorns = "KerrSchild" | ||
|
||
ADMBase::initial_data = "KerrSchild" | ||
ADMBase::initial_lapse = "KerrSchild" | ||
ADMBase::initial_shift = "KerrSchild" | ||
ADMBase::initial_dtlapse = "KerrSchild" | ||
ADMBase::initial_dtshift = "KerrSchild" | ||
|
||
KerrSchild::theta = 0.4 # in radians, approx. 23 degrees | ||
KerrSchild::M = 1 | ||
KerrSchild::a = 0.6 | ||
|
||
ActiveThorns = "NoExcision" | ||
NoExcision::num_regions = 1 | ||
NoExcision::radius[0] = 0.8 | ||
NoExcision::Minkowski_scale[0] = 2.6637047122099 # value found at r=0.875 | ||
NoExcision::reduce_rhs[0] = yes | ||
|
||
ActiveThorns = "Time" | ||
Time::dtfac = 0.25 | ||
|
||
ActiveThorns = "MoL" | ||
MoL::ODE_Method = "RK3" | ||
|
||
ActiveThorns = "ML_BSSN ML_BSSN_Helper GenericFD LoopControl TmunuBase NewRad" | ||
ADMBase::evolution_method = ML_BSSN | ||
ML_BSSN::fdOrder = 2 | ||
ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" | ||
ML_BSSN::my_rhs_boundary_condition = "NewRad" | ||
# These are necessary | ||
ML_BSSN::ML_log_confac_bound = "none" | ||
ML_BSSN::ML_metric_bound = "none" | ||
ML_BSSN::ML_Gamma_bound = "none" | ||
ML_BSSN::ML_trace_curv_bound = "none" | ||
ML_BSSN::ML_curv_bound = "none" | ||
ML_BSSN::ML_lapse_bound = "none" | ||
ML_BSSN::ML_dtlapse_bound = "none" | ||
ML_BSSN::ML_shift_bound = "none" | ||
ML_BSSN::ML_dtshift_bound = "none" | ||
|
||
ActiveThorns = "ADMConstraints" | ||
|
||
ActiveThorns = "AHFinderDirect" | ||
AHFinderDirect::output_OpenDX_control_files = no | ||
AHFinderDirect::geometry_interpolator_name = "Lagrange polynomial interpolation" | ||
AHFinderDirect::geometry_interpolator_pars = "order=2" | ||
AHFinderDirect::surface_interpolator_name = "Lagrange polynomial interpolation" | ||
AHFinderDirect::surface_interpolator_pars = "order=2" | ||
|
||
ActiveThorns = "IOBasic" | ||
IOBasic::outInfo_every = 1 | ||
IOBasic::outInfo_reductions = "norm2" | ||
IOBasic::outInfo_vars = " | ||
ADMConstraints::ham | ||
ADMBase::gxx | ||
" | ||
|
||
IOBasic::outScalar_every = 1 | ||
IOBasic::outScalar_reductions = "average count maximum minimum norm1 norm2 norm_inf" | ||
IOBasic::outScalar_vars = " | ||
ADMBase::metric | ||
ADMBase::curv | ||
ADMBase::lapse | ||
ADMBase::shift | ||
ADMConstraints::hamiltonian | ||
ADMConstraints::momentum | ||
SpaceMask::mask | ||
" | ||
|
||
ActiveThorns = "IOASCII" | ||
IOASCII::out1D_every = 1 | ||
IOASCII::out1D_vars = " | ||
ADMBase::metric | ||
ADMBase::curv | ||
ADMBase::lapse | ||
ADMBase::shift | ||
ADMConstraints::hamiltonian | ||
ADMConstraints::momentum | ||
SpaceMask::mask | ||
" |
45 changes: 45 additions & 0 deletions
45
EinsteinExact_Test/test/KS-tilted-EE/BH_diagnostics.ah1.gp
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,45 @@ | ||
# apparent horizon 1/1 | ||
# | ||
# column 1 = cctk_iteration | ||
# column 2 = cctk_time | ||
# column 3 = centroid_x | ||
# column 4 = centroid_y | ||
# column 5 = centroid_z | ||
# column 6 = min radius | ||
# column 7 = max radius | ||
# column 8 = mean radius | ||
# column 9 = quadrupole_xx | ||
# column 10 = quadrupole_xy | ||
# column 11 = quadrupole_xz | ||
# column 12 = quadrupole_yy | ||
# column 13 = quadrupole_yz | ||
# column 14 = quadrupole_zz | ||
# column 15 = min x | ||
# column 16 = max x | ||
# column 17 = min y | ||
# column 18 = max y | ||
# column 19 = min z | ||
# column 20 = max z | ||
# column 21 = xy-plane circumference | ||
# column 22 = xz-plane circumference | ||
# column 23 = yz-plane circumference | ||
# column 24 = ratio of xz/xy-plane circumferences | ||
# column 25 = ratio of yz/xy-plane circumferences | ||
# column 26 = area | ||
# column 27 = m_irreducible | ||
# column 28 = areal radius | ||
# column 29 = expansion Theta_(l) | ||
# column 30 = inner expansion Theta_(n) | ||
# column 31 = product of the expansions | ||
# column 32 = mean curvature | ||
# column 33 = gradient of the areal radius | ||
# column 34 = gradient of the expansion Theta_(l) | ||
# column 35 = gradient of the inner expansion Theta_(n) | ||
# column 36 = gradient of the product of the expansions | ||
# column 37 = gradient of the mean curvature | ||
# column 38 = minimum of the mean curvature | ||
# column 39 = maximum of the mean curvature | ||
# column 40 = integral of the mean curvature | ||
0 0.000 -0.000000 -0.000000 0.000000 1.797436675 1.895003149 1.862690223 1.196781968 -4.924942002e-06 2.312431388e-07 1.178602681 0.04303625084 1.095066992 -1.894581621 1.894581621 -1.878779470 1.878779470 -1.810946770 1.810946770 12.39444852 11.73317682 11.60097103 0.9466477508 0.9359812184 45.12053438 0.9474410326 1.894882065 1.577582837e-15 -1.469351101 -2.315736692e-15 0.7346755503 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.7086609142 0.7502906623 33.17672887 | ||
1 0.031 -0.000000 -0.000000 -0.000000 1.797372448 1.894934281 1.862626252 1.196702580 -5.225525238e-06 -1.878491706e-08 1.178523290 0.04303781555 1.094987529 -1.894534939 1.894534939 -1.878730681 1.878730681 -1.810890209 1.810890209 12.39460725 11.73324945 11.60103583 0.9466414875 0.9359744600 45.12029054 0.9474384725 1.894876945 -2.909184128e-15 -1.469190774 4.167231013e-15 0.7345953872 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.7086091918 0.7501493126 33.17302269 | ||
2 0.062 -0.000000 0.000000 0.000000 1.797232278 1.894796545 1.862490582 1.196534815 -5.332127631e-06 -5.600195251e-07 1.178355203 0.04303981998 1.094818168 -1.894418286 1.894418286 -1.878610525 1.878610525 -1.810757789 1.810757789 12.39476525 11.73332824 11.60110855 0.9466357775 0.9359683964 45.12015820 0.9474370831 1.894874166 -7.304396136e-14 -1.468873330 1.025105290e-13 0.7344366648 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.7084954921 0.7499218727 33.16583812 |
158 changes: 158 additions & 0 deletions
158
EinsteinExact_Test/test/KS-tilted-EE/alp_3D_diagonal.xg
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,158 @@ | ||
"x-label diagonal | ||
"y-label ADMBASE::alp | ||
|
||
|
||
"Time = 0.0000000000000 | ||
-5.1961524227066 0.8491838137798 | ||
-4.9796460717605 0.8440302314853 | ||
-4.7631397208144 0.8385063860414 | ||
-4.5466333698683 0.8325701064365 | ||
-4.3301270189222 0.8261724128828 | ||
-4.1136206679761 0.8192560464056 | ||
-3.8971143170300 0.8117535880699 | ||
-3.6806079660839 0.8035850231432 | ||
-3.4641016151378 0.7946545420772 | ||
-3.2475952641916 0.7848462727130 | ||
-3.0310889132455 0.7740184844577 | ||
-2.8145825622994 0.7619955561369 | ||
-2.5980762113533 0.7485565831352 | ||
-2.3815698604072 0.7334187813939 | ||
-2.1650635094611 0.7162125660409 | ||
-1.9485571585150 0.6964428409450 | ||
-1.7320508075689 0.6734267500570 | ||
-1.5155444566228 0.6461911061836 | ||
-1.2990381056767 0.6133094408108 | ||
-1.0825317547305 0.5727502701744 | ||
-0.8660254037844 0.5231578523833 | ||
-0.6495190528383 1.0000000000000 | ||
-0.4330127018922 1.0000000000000 | ||
-0.2165063509461 1.0000000000000 | ||
0.0000000000000 1.0000000000000 | ||
0.2165063509461 1.0000000000000 | ||
0.4330127018922 1.0000000000000 | ||
0.6495190528383 1.0000000000000 | ||
0.8660254037844 0.5231578523833 | ||
1.0825317547305 0.5727502701744 | ||
1.2990381056767 0.6133094408108 | ||
1.5155444566228 0.6461911061836 | ||
1.7320508075689 0.6734267500570 | ||
1.9485571585150 0.6964428409450 | ||
2.1650635094611 0.7162125660409 | ||
2.3815698604072 0.7334187813939 | ||
2.5980762113533 0.7485565831352 | ||
2.8145825622994 0.7619955561369 | ||
3.0310889132455 0.7740184844577 | ||
3.2475952641916 0.7848462727130 | ||
3.4641016151378 0.7946545420772 | ||
3.6806079660839 0.8035850231432 | ||
3.8971143170300 0.8117535880699 | ||
4.1136206679761 0.8192560464056 | ||
4.3301270189222 0.8261724128828 | ||
4.5466333698683 0.8325701064365 | ||
4.7631397208144 0.8385063860414 | ||
4.9796460717605 0.8440302314853 | ||
5.1961524227066 0.8491838137798 | ||
|
||
|
||
"Time = 0.0312500000000 | ||
-5.1961524227066 0.8493721158841 | ||
-4.9796460717605 0.8442397867995 | ||
-4.7631397208144 0.8370507547443 | ||
-4.5466333698683 0.8307560075103 | ||
-4.3301270189222 0.8242677587729 | ||
-4.1136206679761 0.8172002921121 | ||
-3.8971143170300 0.8095213566283 | ||
-3.6806079660839 0.8011528398115 | ||
-3.4641016151378 0.7919935664094 | ||
-3.2475952641916 0.7819216486263 | ||
-3.0310889132455 0.7707876181682 | ||
-2.8145825622994 0.7584057354483 | ||
-2.5980762113533 0.7445416057973 | ||
-2.3815698604072 0.7288940401085 | ||
-2.1650635094611 0.7110676256931 | ||
-1.9485571585150 0.6905297180928 | ||
-1.7320508075689 0.6665402985916 | ||
-1.5155444566228 0.6380334670640 | ||
-1.2990381056767 0.6034175782281 | ||
-1.0825317547305 0.5603066871570 | ||
-0.8660254037844 0.5231578523833 | ||
-0.6495190528383 1.0000000000000 | ||
-0.4330127018922 1.0000000000000 | ||
-0.2165063509461 1.0000000000000 | ||
0.0000000000000 1.0000000000000 | ||
0.2165063509461 1.0000000000000 | ||
0.4330127018922 1.0000000000000 | ||
0.6495190528383 1.0000000000000 | ||
0.8660254037844 0.5231578523833 | ||
1.0825317547305 0.5603066871570 | ||
1.2990381056767 0.6034175782281 | ||
1.5155444566228 0.6380334670640 | ||
1.7320508075689 0.6665402985916 | ||
1.9485571585150 0.6905297180928 | ||
2.1650635094611 0.7110676256931 | ||
2.3815698604072 0.7288940401085 | ||
2.5980762113533 0.7445416057973 | ||
2.8145825622994 0.7584057354483 | ||
3.0310889132455 0.7707876181682 | ||
3.2475952641916 0.7819216486263 | ||
3.4641016151378 0.7919935664094 | ||
3.6806079660839 0.8011528398115 | ||
3.8971143170300 0.8095213566283 | ||
4.1136206679761 0.8172002921121 | ||
4.3301270189222 0.8242677587729 | ||
4.5466333698683 0.8307560075103 | ||
4.7631397208144 0.8370507547443 | ||
4.9796460717605 0.8442397867995 | ||
5.1961524227066 0.8493721158841 | ||
|
||
|
||
"Time = 0.0625000000000 | ||
-5.1961524227066 0.8495625621014 | ||
-4.9796460717605 0.8444066077885 | ||
-4.7631397208144 0.8359422474452 | ||
-4.5466333698683 0.8289201450046 | ||
-4.3301270189222 0.8223583165762 | ||
-4.1136206679761 0.8151637556547 | ||
-3.8971143170300 0.8073123821167 | ||
-3.6806079660839 0.7987477552586 | ||
-3.4641016151378 0.7893652929364 | ||
-3.2475952641916 0.7790369050674 | ||
-3.0310889132455 0.7676059290648 | ||
-2.8145825622994 0.7548773465080 | ||
-2.5980762113533 0.7406045202802 | ||
-2.3815698604072 0.7244697899727 | ||
-2.1650635094611 0.7060550269194 | ||
-1.9485571585150 0.6847952775622 | ||
-1.7320508075689 0.6599026981887 | ||
-1.5155444566228 0.6302365982193 | ||
-1.2990381056767 0.5940782165569 | ||
-1.0825317547305 0.5487494070793 | ||
-0.8660254037844 0.5231578523833 | ||
-0.6495190528383 1.0000000000000 | ||
-0.4330127018922 1.0000000000000 | ||
-0.2165063509461 1.0000000000000 | ||
0.0000000000000 1.0000000000000 | ||
0.2165063509461 1.0000000000000 | ||
0.4330127018922 1.0000000000000 | ||
0.6495190528383 1.0000000000000 | ||
0.8660254037844 0.5231578523833 | ||
1.0825317547305 0.5487494070793 | ||
1.2990381056767 0.5940782165569 | ||
1.5155444566228 0.6302365982193 | ||
1.7320508075689 0.6599026981887 | ||
1.9485571585150 0.6847952775622 | ||
2.1650635094611 0.7060550269194 | ||
2.3815698604072 0.7244697899727 | ||
2.5980762113533 0.7406045202802 | ||
2.8145825622994 0.7548773465080 | ||
3.0310889132455 0.7676059290648 | ||
3.2475952641916 0.7790369050674 | ||
3.4641016151378 0.7893652929364 | ||
3.6806079660839 0.7987477552586 | ||
3.8971143170300 0.8073123821167 | ||
4.1136206679761 0.8151637556547 | ||
4.3301270189222 0.8223583165762 | ||
4.5466333698683 0.8289201450046 | ||
4.7631397208144 0.8359422474452 | ||
4.9796460717605 0.8444066077885 | ||
5.1961524227066 0.8495625621014 |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,6 @@ | ||
"x-label time | ||
"y-label ADMBASE::alp | ||
"alp v time | ||
0.0000000000000 0.7647839801295 | ||
0.0312500000000 0.7617064217053 | ||
0.0625000000000 0.7586685846462 |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,6 @@ | ||
"x-label time | ||
"y-label ADMBASE::alp | ||
"alp v time | ||
0.0000000000000 117649.0000000000000 | ||
0.0312500000000 117649.0000000000000 | ||
0.0625000000000 117649.0000000000000 |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,6 @@ | ||
"x-label time | ||
"y-label ADMBASE::alp | ||
"alp v time | ||
0.0000000000000 1.0000000000000 | ||
0.0312500000000 1.0000000000000 | ||
0.0625000000000 1.0000000000000 |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,6 @@ | ||
"x-label time | ||
"y-label ADMBASE::alp | ||
"alp v time | ||
0.0000000000000 0.4579988863699 | ||
0.0312500000000 0.4579988863699 | ||
0.0625000000000 0.4579988863699 |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,6 @@ | ||
"x-label time | ||
"y-label ADMBASE::alp | ||
"alp v time | ||
0.0000000000000 0.7647839801295 | ||
0.0312500000000 0.7617064217053 | ||
0.0625000000000 0.7586685846462 |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,6 @@ | ||
"x-label time | ||
"y-label ADMBASE::alp | ||
"alp v time | ||
0.0000000000000 0.7669778982599 | ||
0.0312500000000 0.7640654264959 | ||
0.0625000000000 0.7611917276668 |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,6 @@ | ||
"x-label time | ||
"y-label ADMBASE::alp | ||
"alp v time | ||
0.0000000000000 1.0000000000000 | ||
0.0312500000000 1.0000000000000 | ||
0.0625000000000 1.0000000000000 |
Oops, something went wrong.