Skip to content

Commit

Permalink
TestRKAB: add unit test
Browse files Browse the repository at this point in the history
  • Loading branch information
@lwJi
lwJi committed Sep 20, 2024
1 parent 18c1dc9 commit a6682af
Show file tree
Hide file tree
Showing 19 changed files with 1,627 additions and 0 deletions.
79 changes: 79 additions & 0 deletions TestRKAB/test/gaussian.par
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,79 @@
ActiveThorns = "
BoxInBox
CarpetX
IOUtil
ODESolvers
TestRKAB
"

$out_every = 64
$nlevels = 3
$ncells = 32

TestRKAB::initial_condition = "Gaussian"
TestRKAB::gaussian_width = 0.17677669529 # sqrt(2)*W = 0.25

CarpetX::poison_undefined_values = yes
CarpetX::verbose = no

Cactus::cctk_show_schedule = yes
Cactus::presync_mode = "mixed-error"

CarpetX::xmin = -4.0
CarpetX::ymin = -4.0
CarpetX::zmin = -4.0

CarpetX::xmax = +4.0
CarpetX::ymax = +4.0
CarpetX::zmax = +4.0

CarpetX::ncells_x = $ncells
CarpetX::ncells_y = $ncells
CarpetX::ncells_z = $ncells

CarpetX::periodic_x = yes
CarpetX::periodic_y = yes
CarpetX::periodic_z = yes

CarpetX::max_num_levels = $nlevels
CarpetX::regrid_every = 128
CarpetX::regrid_error_threshold = 0.9

CarpetX::prolongation_type = "ddf"
#CarpetX::prolongation_order = 5
CarpetX::ghost_size = 3
CarpetX::dtfac = 0.25

CarpetX::blocking_factor_x = 2
CarpetX::blocking_factor_y = 2
CarpetX::blocking_factor_z = 2
#CarpetX::max_grid_size_x = 10000000
#CarpetX::max_grid_size_y = 10000000
#CarpetX::max_grid_size_z = 10000000
#CarpetX::max_tile_size_x = 10000000
#CarpetX::max_tile_size_y = 10000000
#CarpetX::max_tile_size_z = 10000000

BoxInBox::num_regions = 1
# Region 1
BoxInBox::shape_1 = "cube"
BoxInBox::num_levels_1 = $nlevels
BoxInBox::radius_1 = [ -1.0, 1.0, 0.5 ]


Cactus::terminate = "time"
Cactus::cctk_final_time = 2.0
#Cactus::cctk_itlast = $itlast

ODESolvers::method = "RKAB"
ODESolvers::verbose = no

IO::out_dir = $parfile
IO::out_every = $out_every

#CarpetX::out_norm_vars = "all"

CarpetX::out_tsv_vars = "
TestRKAB::state
TestRKAB::error
"
86 changes: 86 additions & 0 deletions TestRKAB/test/gaussian/testrkab-error.it000000.x.tsv
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,86 @@
# 1:iteration 2:time 3:patch 4:level 5:i 6:j 7:k 8:x 9:y 10:z 11:u_err 12:rho_err
0 0.0000000000000000e+00 0 0 -3 16 16 -4.7500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 -2 16 16 -4.5000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 -1 16 16 -4.2500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 0 16 16 -4.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 1 16 16 -3.7500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 2 16 16 -3.5000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 3 16 16 -3.2500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 4 16 16 -3.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 5 16 16 -2.7500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 6 16 16 -2.5000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 7 16 16 -2.2500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 8 16 16 -2.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 9 16 16 -1.7500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 10 16 16 -1.5000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 11 16 16 -1.2500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 12 16 16 -1.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 13 16 16 -7.5000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 14 16 16 -5.0000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 15 16 16 -2.5000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 16 16 16 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 17 16 16 2.5000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 18 16 16 5.0000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 19 16 16 7.5000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 20 16 16 1.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 21 16 16 1.2500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 22 16 16 1.5000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 23 16 16 1.7500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 24 16 16 2.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 25 16 16 2.2500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 26 16 16 2.5000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 27 16 16 2.7500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 28 16 16 3.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 29 16 16 3.2500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 30 16 16 3.5000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 31 16 16 3.7500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 32 16 16 4.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 33 16 16 4.2500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 34 16 16 4.5000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 0 35 16 16 4.7500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 21 32 32 -1.3750000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 22 32 32 -1.2500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 23 32 32 -1.1250000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 24 32 32 -1.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 25 32 32 -8.7500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 26 32 32 -7.5000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 27 32 32 -6.2500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 28 32 32 -5.0000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 29 32 32 -3.7500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 30 32 32 -2.5000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 31 32 32 -1.2500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 32 32 32 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 33 32 32 1.2500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 34 32 32 2.5000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 35 32 32 3.7500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 36 32 32 5.0000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 37 32 32 6.2500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 38 32 32 7.5000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 39 32 32 8.7500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 40 32 32 1.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 41 32 32 1.1250000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 42 32 32 1.2500000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 1 43 32 32 1.3750000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 53 64 64 -6.8750000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 54 64 64 -6.2500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 55 64 64 -5.6250000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 56 64 64 -5.0000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 57 64 64 -4.3750000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 58 64 64 -3.7500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 59 64 64 -3.1250000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 60 64 64 -2.5000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 61 64 64 -1.8750000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 62 64 64 -1.2500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 63 64 64 -6.2500000000000000e-02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 64 64 64 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 65 64 64 6.2500000000000000e-02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 66 64 64 1.2500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 67 64 64 1.8750000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 68 64 64 2.5000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 69 64 64 3.1250000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 70 64 64 3.7500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 71 64 64 4.3750000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 72 64 64 5.0000000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 73 64 64 5.6250000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 74 64 64 6.2500000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
0 0.0000000000000000e+00 0 2 75 64 64 6.8750000000000000e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
Loading

0 comments on commit a6682af

Please sign in to comment.