Initializing CUDA...
CUDA initialized with 1 GPU
AMReX (21.07-8-gfffd3d976c7e) initialized
Calling Setup()
Calling ReadParameters()
WARNING: burning_cutoff_density_lo not supplied in the inputs file
WARNING: setting burning_cutoff_density_lo = base_cutoff_density
reading extern runtime parameters ...
Calling VariableSetup()
Initializing Helmholtz EOS and using Coulomb corrections.
Calling BCSetup()
Calling BaseStateGeometry::Init()
Calling Init()
Calling InitData()
initdata model_File = model.hse.cool.coulomb
model file = model.hse.cool.coulomb
reading initial model
640 points found in the initial model
6 variables found in the initial model
model file mapping, level: 0
dr of MAESTRO base state = 562500.000000
dr of input file data = 562500.000000
maximum radius (cell-centered) of input model = 359718750.000000
setting r_cutoff to 365
radius at r_cutoff 205593750
Maximum HSE Error = 0.000061
(after putting initial model into base state arrays, and
for density < base_cutoff_density)
Writing plotfile reacting_bubble-3d_pltInitData after InitData
inner sponge: r_sp , r_tp : 186468750, 224718750
Time to write plotfile: 0.524370425
inner sponge: r_sp , r_tp : 186468750, 224718750
Doing initial projection
Calling nodal solver
MLMG: Initial rhs = 0
MLMG: Initial residual (resid0) = 0
MLMG: No iterations needed
MLMG: Timers: Solve = 0.077667208 Iter = 0 Bottom = 0
Done calling nodal solver
Writing plotfile reacting_bubble-3d_pltafter_InitProj after InitProj
inner sponge: r_sp , r_tp : 186468750, 224718750
Time to write plotfile: 0.427095466
Call to firstdt for level 0 gives dt_lev = 0.0004741795139
Multiplying dt_lev by init_shrink; dt_lev = 0.0004741795139
Minimum firstdt over all levels = 0.0004741795139
Doing initial divu iteration #1
Calling nodal solver
MLMG: Initial rhs = 2054.815795
MLMG: Initial residual (resid0) = 2054.815795
MLMG: Final Iter. 7 resid, resid/bnorm = 3.94901889e-10, 1.921835961e-13
MLMG: Timers: Solve = 3.400218369 Iter = 3.327084164 Bottom = 0.024300881
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 0.1719927642
Minimum estdt over all levels = 0.1719927642
Call to estdt at end of istep_divu_iter = 1 gives dt = 0.1719927642
Multiplying dt by init_shrink; dt = 0.1719927642
Ignoring this new dt since it's larger than the previous dt = 0.0004741795139
Writing plotfile reacting_bubble-3d_pltafter_DivuIter after final DivuIter
inner sponge: r_sp , r_tp : 186468750, 224718750
Time to write plotfile: 0.400389066
Doing initial pressure iteration #1
Timestep 0 starts with TIME = 0 DT = 0.0004741795139
Cell Count:
Level 0, 1327104 cells
inner sponge: r_sp , r_tp : 186468750, 224718750
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 2728631.974
MLMG: Initial residual (resid0) = 2728631.974
MLMG: Final Iter. 9 resid, resid/bnorm = 1.603597775e-05, 5.876929504e-12
MLMG: Timers: Solve = 0.201021686 Iter = 0.195757424 Bottom = 0.143742222
<<< STEP 4 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs = 2728627.802
MLMG: Initial residual (resid0) = 8728.675618
MLMG: Final Iter. 5 resid, resid/bnorm = 0.0002170299413, 7.953812575e-11
MLMG: Timers: Solve = 0.118763655 Iter = 0.116533357 Bottom = 0.087834516
<<< STEP 8 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs = 1.15755262e+10
MLMG: Initial residual (resid0) = 1.15755262e+10
MLMG: Final Iter. 7 resid, resid/bnorm = 0.01014206602, 8.761645772e-13
MLMG: Timers: Solve = 3.42500241 Iter = 3.351229354 Bottom = 0.024139712
Done calling nodal solver
Timestep 0 ends with TIME = 0.0004741795139 DT = 0.0004741795139
Timing summary:
Advection :0.793498306 seconds
MAC Proj :0.339609989 seconds
Nodal Proj :3.553173624 seconds
Reactions :0.112502932 seconds
Misc :0.04427762 seconds
Base State :0.012134834 seconds
Time to advance time step: 4.843296088
Writing plotfile 0 after all initialization
inner sponge: r_sp , r_tp : 186468750, 224718750
Time to write plotfile: 0.392625235
Calling Evolve()
Timestep 1 starts with TIME = 0 DT = 0.0004741795139
Cell Count:
Level 0, 1327104 cells
inner sponge: r_sp , r_tp : 186468750, 224718750
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 1308380.202
MLMG: Initial residual (resid0) = 1308380.202
MLMG: Final Iter. 8 resid, resid/bnorm = 0.0001284286409, 9.815850219e-11
MLMG: Timers: Solve = 0.172806654 Iter = 0.170436386 Bottom = 0.124745051
<<< STEP 4 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs = 1308380.239
MLMG: Initial residual (resid0) = 7.635486734
MLMG: Final Iter. 5 resid, resid/bnorm = 4.092841049e-05, 3.128174002e-11
MLMG: Timers: Solve = 0.104695687 Iter = 0.102396745 Bottom = 0.073726052
<<< STEP 8 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs = 5488883.997
MLMG: Initial residual (resid0) = 5488883.997
MLMG: Final Iter. 7 resid, resid/bnorm = 5.941594922e-06, 1.082477772e-12
MLMG: Timers: Solve = 3.441144315 Iter = 3.370856267 Bottom = 0.024120968
Done calling nodal solver
Timestep 1 ends with TIME = 0.0004741795139 DT = 0.0004741795139
Timing summary:
Advection :0.754304619 seconds
MAC Proj :0.289273118 seconds
Nodal Proj :3.56947723 seconds
Reactions :0.110956487 seconds
Misc :0.023009024 seconds
Base State :0.013949667 seconds
Time to advance time step: 4.747300254
Call to estdt for level 0 gives dt_lev = 0.2012294727
Minimum estdt over all levels = 0.2012294727
Call to estdt at beginning of step 2 gives dt =0.2012294727
dt_growth factor limits the new dt = 0.0005215974653
Timestep 2 starts with TIME = 0.0004741795139 DT = 0.0005215974653
Cell Count:
Level 0, 1327104 cells
inner sponge: r_sp , r_tp : 186468750, 224718750
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 1426181.41
MLMG: Initial residual (resid0) = 1426181.41
MLMG: Final Iter. 9 resid, resid/bnorm = 3.271783318e-05, 2.294086359e-11
MLMG: Timers: Solve = 0.212795059 Iter = 0.21055547 Bottom = 0.158938737
<<< STEP 4 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs = 1426291.688
MLMG: Initial residual (resid0) = 36463.44575
MLMG: Final Iter. 8 resid, resid/bnorm = 0.0001356700103, 9.512080274e-11
MLMG: Timers: Solve = 0.14329961 Iter = 0.141072147 Bottom = 0.095287996
<<< STEP 8 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs = 6063485.842
MLMG: Initial residual (resid0) = 6063485.842
MLMG: Final Iter. 8 resid, resid/bnorm = 4.934472649e-05, 8.138012981e-12
MLMG: Timers: Solve = 3.940807219 Iter = 3.867259032 Bottom = 0.027098307
Done calling nodal solver
Timestep 2 ends with TIME = 0.0009957769792 DT = 0.0005215974653
Timing summary:
Advection :0.753095152 seconds
MAC Proj :0.367825727 seconds
Nodal Proj :4.068817655 seconds
Reactions :0.111004408 seconds
Misc :0.023217658 seconds
Base State :0.012865646 seconds
Time to advance time step: 5.324126178
Call to estdt for level 0 gives dt_lev = 0.201266633
Minimum estdt over all levels = 0.201266633
Call to estdt at beginning of step 3 gives dt =0.201266633
dt_growth factor limits the new dt = 0.0005737572118
Timestep 3 starts with TIME = 0.0009957769792 DT = 0.0005737572118
Cell Count:
Level 0, 1327104 cells
inner sponge: r_sp , r_tp : 186468750, 224718750
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 1597824.953
MLMG: Initial residual (resid0) = 1597824.953
MLMG: Final Iter. 9 resid, resid/bnorm = 0.0001106822459, 6.927057038e-11
MLMG: Timers: Solve = 0.181484424 Iter = 0.179214219 Bottom = 0.12795068
<<< STEP 4 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs = 1597946.234
MLMG: Initial residual (resid0) = 56931.55491
MLMG: Final Iter. 9 resid, resid/bnorm = 2.066239267e-05, 1.293059318e-11
MLMG: Timers: Solve = 0.186451517 Iter = 0.184232483 Bottom = 0.132989139
<<< STEP 8 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs = 6669975.866
MLMG: Initial residual (resid0) = 6669975.866
MLMG: Final Iter. 8 resid, resid/bnorm = 2.929759992e-05, 4.392459659e-12
MLMG: Timers: Solve = 3.86898938 Iter = 3.7986089 Bottom = 0.027428549
Done calling nodal solver
Timestep 3 ends with TIME = 0.001569534191 DT = 0.0005737572118
Timing summary:
Advection :0.742030893 seconds
MAC Proj :0.37981964 seconds
Nodal Proj :3.994088222 seconds
Reactions :0.111227177 seconds
Misc :0.02451876 seconds
Base State :0.013462769 seconds
Time to advance time step: 5.251847567
Writing plotfile 3
inner sponge: r_sp , r_tp : 186468750, 224718750
Time to write plotfile: 0.390275262
Total Time: 29.43911162
Unused ParmParse Variables:
[TOP]::maestro.v(nvals = 1) :: [1]
[TOP]::amr.ref_ratio(nvals = 4) :: [2, 2, 2, 2]
[TOP]::amr.check_file(nvals = 1) :: [reacting_bubble-3d_chk]
[TOP]::amr.checkpoint_files_output(nvals = 1) :: [0]
Total GPU global memory (MB): 12066
Free GPU global memory (MB): 8718
[The Arena] space allocated (MB): 9049
[The Arena] space used (MB): 0
[The Pinned Arena] space allocated (MB): 132
[The Pinned Arena] space used (MB): 0
AMReX (21.07-8-gfffd3d976c7e) finalized