MPI initialized with 4 MPI processes MPI initialized with thread support level 3 OMP initialized with 4 OMP threads AMReX (23.10-13-gecaa14456e3f) initialized Starting run at 06:46:27 UTC on 2023-10-13. Successfully read inputs file ... Castro git describe: 23.10-35-gf8446cd94 AMReX git describe: 23.10-13-gecaa14456 Microphysics git describe: 23.10-6-g657406dc reading extern runtime parameters ... 13 Species: He4 C12 O16 Ne20 Mg24 Si28 S32 Ar36 Ca40 Ti44 Cr48 Fe52 Ni56 Successfully read inputs file ... Generated initial model for primary WD of mass 0.9 solar masses, central density 1.574e+07 g cm**-3, and radius 6.747e+08 cm. Generated initial model for secondary WD of mass 0.9 solar masses, central density 1.574e+07 g cm**-3, and radius 6.747e+08 cm. Generated binary orbit of distance 1.781e+09 cm = 1.190e-04 AU. The primary orbits the center of mass at distance -8.903e+08 cm = -5.951e-05 AU. The secondary orbits the center of mass at distance 8.903e+08 cm = 5.951e-05 AU. The initial orbital period is 30.543687 s. The initial orbital speed of the primary is -1.831e+08 cm/s. The initial orbital speed of the secondary is 1.831e+08 cm/s. Castro::numpts_1d at level 0 is 174 Setting the gravity type to PoissonGrav Initializing the data at level 0 Done initializing the level 0 data ...estimated hydro-limited timestep at level 0: 0.05426072637 ...hydro-limited CFL timestep constrained at (i,j,k) = (64,0,0) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.05426072637 grid_places() time: 0.03482254 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.034776853 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.369164632 Gravity::multilevel_solve_for_new_phi() time = 1.806306752 Gravity::fill_multipole_BCs() time = 0.370127889 ... norm of composite residual at level 0 3.509358608e-10 INITIAL GRIDS Level 0 27 grids 884736 cells 100 % of domain smallest grid: 32 x 32 x 32 biggest grid: 32 x 32 x 32 PLOTFILE: file = wdmerger_3D_plt00000 Write plotfile time = 1.086336246 seconds ...estimated hydro-limited timestep at level 0: 0.05426072637 ...hydro-limited CFL timestep constrained at (i,j,k) = (64,0,0) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.05426072637 [Level 0 step 1] ADVANCE at time 0 with dt = 0.005426072637  Beginning subcycle 1 starting at time 0 with dt = 0.005426072637 Estimated number of subcycles remaining: 1 Castro::construct_old_gravity() time = 0.869907 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -3.287624606e+47 ener added: 6.129070894e+49 Contributions to the state from old-time sources: mass added: 0 xmom added: 1.449362815e+25 ymom added: -1.046793351e+26 zmom added: -4.133074948e+24 eint added: 0 ener added: -1.270533196e+50 Castro::do_old_sources() time = 2.341768878 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 9.00547245e+47 ener added: 1.651139801e+46 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 20.71990059 on level 0 ... new-time composite Poisson gravity solve from level 0 to level 0 Gravity::fill_multipole_BCs() time = 0.650990975 Gravity::multilevel_solve_for_new_phi() time = 1.719004261 Castro::construct_new_gravity() time = 1.865160207 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.440206272e+49 ener added: 2.911791643e+45 Contributions to the state from new-time sources: mass added: 0 xmom added: 1.451207927e+26 ymom added: 1.075144428e+26 zmom added: -5.613375295e+26 eint added: 0 ener added: 1.260457978e+50 Castro::do_new_sources() time = 2.396836246 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.184461852e+27 ener added: 2.243834889e+27 ...estimated hydro-limited timestep at level 0: 0.05461684804 ...hydro-limited CFL timestep constrained at (i,j,k) = (0,29,94) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.05461684804  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.02466204421 [Level 0 step 1] Advanced 884736 cells Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 5.190352593e+26 ener added: 6.213783704e+26 Updating the rotational period from 30.54368673 s to 30.45939751 s. [STEP 1] Coarse TimeStep time: 38.04053681 [STEP 1] FAB kilobyte spread across MPI nodes: [679617 ... 761594] STEP = 1 TIME = 0.005426072637 DT = 0.005426072637 ...estimated hydro-limited timestep at level 0: 0.05461684804 ...hydro-limited CFL timestep constrained at (i,j,k) = (0,29,94) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.05461684804 Castro::compute_new_dt : limiting dt at level 0 ... new dt computed: 0.05461684804 ... but limiting to: 0.0059686799 = 1.1 * 0.005426072637 Now regridding at level lbase = 0 grid_places() time: 0.033338256 new finest: 1 Castro::numpts_1d at level 1 is 340 Gravity::fill_multipole_BCs() time = 0.417467468 Gravity::multilevel_solve_for_new_phi() time = 1.985879734 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.306135866e+46 ener added: 9.211440366e+51 STEP = 1 TIME = 0.005426072637 : REGRID with lbase = 0 Level 1 4 grids 65536 cells 0.9259259259 % of domain smallest grid: 32 x 32 x 16 biggest grid: 32 x 32 x 16 ...estimated hydro-limited timestep at level 0: 0.05461684804 ...hydro-limited CFL timestep constrained at (i,j,k) = (0,29,94) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.05461684804 ...estimated hydro-limited timestep at level 1: 2.526135286e-09 ...hydro-limited CFL timestep constrained at (i,j,k) = (87,103,92) Castro::estTimeStep (hydro-limited) at level 1: estdt = 2.526135286e-09 [Level 0 step 2] ADVANCE at time 0.005426072637 with dt = 2.526135286e-09 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.849569753e+31 ener added: 2.785937322e+31  Beginning subcycle 1 starting at time 0.005426072637 with dt = 2.526135286e-09 Estimated number of subcycles remaining: 1 Castro::construct_old_gravity() time = 0.789683972 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.200514105e+40 ener added: 6.433158191e+33 Contributions to the state from old-time sources: mass added: 0 xmom added: 5.533896722e+24 ymom added: 4.756222073e+26 zmom added: -2.866993882e+23 eint added: 0 ener added: -1.721443002e+45 Castro::do_old_sources() time = 2.05677388 on level 0 ... Entering construct_ctu_hydro_source() on level 0 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 21.72245456 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.147330623e+31 ener added: 1.147330701e+31 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 4.551111111e+23 ener added: 6.281481481e+23 Castro::construct_old_gravity() time = 0.042881074 on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.197253775e+40 ener added: 1.79168512e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: 7.571938195e+22 ymom added: 1.564283897e+22 zmom added: 6.441948907e+21 eint added: 0 ener added: -2.586789965e+47 Castro::do_old_sources() time = 0.195928339 on level 1 ... Entering construct_ctu_hydro_source() on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 9.204912576e+51 ener added: 6.747117479e+48 ... Leaving construct_ctu_hydro_source() on level 1 Castro::construct_ctu_hydro_source() time = 2.853483402 on level 1 Castro::reflux() at level 0 : time = 2.36064763 ... new-time composite Poisson gravity solve from level 0 to level 1 Gravity::fill_multipole_BCs() time = 0.660239211 Gravity::multilevel_solve_for_new_phi() time = 1.642955777 Castro::construct_new_gravity() time = 1.757616765 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -3.46375568e+41 ener added: 7.485558923e+35 Contributions to the state from new-time sources: mass added: 0 xmom added: -4.864072325e+21 ymom added: -6.244635709e+21 zmom added: -1.246148803e+20 eint added: 0 ener added: -3.223710603e+41 Castro::do_new_sources() time = 1.98990893 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.14873037e+26 ener added: 4.251851852e+22 ...estimated hydro-limited timestep at level 0: 0.05461716025 ...hydro-limited CFL timestep constrained at (i,j,k) = (0,29,94) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.05461716025 Castro::construct_new_gravity() time = 0.030323508 on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.564837764e+47 ener added: 4.559472528e+41 Contributions to the state from new-time sources: mass added: 0 xmom added: 5.11224561e+18 ymom added: -1.461418079e+19 zmom added: -8.220975268e+19 eint added: 0 ener added: 2.564611685e+47 Castro::do_new_sources() time = 0.146946879 on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 3.281876859e+29 ener added: 3.281875295e+29 ...estimated hydro-limited timestep at level 1: 6.287037368e-07 ...hydro-limited CFL timestep constrained at (i,j,k) = (121,90,90) Castro::estTimeStep (hydro-limited) at level 1: estdt = 6.287037368e-07  Subcycle completed Subcycling complete Zones advanced per microsecond from level 0 to level 1: 0.02249906495 [Level 0 step 2] Advanced 884736 cells [Level 1 step 2] ADVANCE at time 0.005426072637 with dt = 2.526135286e-09 Advance at this level has already been completed. [Level 1 step 2] Advanced 65536 cells Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.036240741e+24 ener added: 5.917722222e+24 Updating the rotational period from 30.45939751 s to 30.52749735 s. [STEP 2] Coarse TimeStep time: 50.501152 [STEP 2] FAB kilobyte spread across MPI nodes: [833440 ... 934735] STEP = 2 TIME = 0.005426075163 DT = 2.526135286e-09 ...estimated hydro-limited timestep at level 0: 0.0008131491259 ...hydro-limited CFL timestep constrained at (i,j,k) = (37,49,42) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.0008131491259 ...estimated hydro-limited timestep at level 1: 6.287037368e-07 ...hydro-limited CFL timestep constrained at (i,j,k) = (121,90,90) Castro::estTimeStep (hydro-limited) at level 1: estdt = 6.287037368e-07 Castro::compute_new_dt : limiting dt at level 0 ... new dt computed: 0.0008131491259 ... but limiting to: 2.778748814e-09 = 1.1 * 2.526135286e-09 Castro::compute_new_dt : limiting dt at level 1 ... new dt computed: 6.287037368e-07 ... but limiting to: 2.778748814e-09 = 1.1 * 2.526135286e-09 Now regridding at level lbase = 0 grid_places() time: 0.032005373 new finest: 1 Castro::numpts_1d at level 1 is 340 Gravity::fill_multipole_BCs() time = 0.405919721 Gravity::multilevel_solve_for_new_phi() time = 0.535042895 STEP = 2 TIME = 0.005426075163 : REGRID with lbase = 0 Level 1 4 grids 65536 cells 0.9259259259 % of domain smallest grid: 32 x 32 x 16 biggest grid: 32 x 32 x 16 ...estimated hydro-limited timestep at level 0: 0.0008131491259 ...hydro-limited CFL timestep constrained at (i,j,k) = (37,49,42) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.0008131491259 ...estimated hydro-limited timestep at level 1: 6.287037368e-07 ...hydro-limited CFL timestep constrained at (i,j,k) = (121,90,90) Castro::estTimeStep (hydro-limited) at level 1: estdt = 6.287037368e-07 [Level 0 step 3] ADVANCE at time 0.005426075163 with dt = 2.778748814e-09  Beginning subcycle 1 starting at time 0.005426075163 with dt = 2.778748814e-09 Estimated number of subcycles remaining: 1 Castro::construct_old_gravity() time = 0.830848507 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: 6.083120164e+24 ymom added: 5.231862796e+26 zmom added: -3.157842738e+23 eint added: 0 ener added: -1.937110674e+45 Castro::do_old_sources() time = 2.068196309 on level 0 ... Entering construct_ctu_hydro_source() on level 0 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 21.46089152 on level 0 Castro::construct_old_gravity() time = 0.044124163 on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.818691488e+40 ener added: 1.212217726e+36 Contributions to the state from old-time sources: mass added: 0 xmom added: 8.443638501e+22 ymom added: 4.427218578e+21 zmom added: 7.10314939e+21 eint added: 0 ener added: -2.344787748e+45 Castro::do_old_sources() time = 0.196115381 on level 1 ... Entering construct_ctu_hydro_source() on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.83673455e+48 ener added: 2.077900047e+42 ... Leaving construct_ctu_hydro_source() on level 1 Castro::construct_ctu_hydro_source() time = 2.851676053 on level 1 Castro::reflux() at level 0 : time = 2.265534176 ... new-time composite Poisson gravity solve from level 0 to level 1 Gravity::fill_multipole_BCs() time = 0.574712249 Gravity::multilevel_solve_for_new_phi() time = 1.1283843 Castro::construct_new_gravity() time = 1.243122136 on level 0 Contributions to the state from new-time sources: mass added: 0 xmom added: -2.108598872e+21 ymom added: -7.90284783e+21 zmom added: 1.10682106e+20 eint added: 0 ener added: -3.427512073e+41 Castro::do_new_sources() time = 1.981166194 on level 0 ...estimated hydro-limited timestep at level 0: 0.0008131520194 ...hydro-limited CFL timestep constrained at (i,j,k) = (37,49,42) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.0008131520194 Castro::construct_new_gravity() time = 0.031023647 on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 5.244373444e+43 ener added: 1.74625952e+38 Contributions to the state from new-time sources: mass added: 0 xmom added: -5.605243277e+17 ymom added: 2.1743379e+19 zmom added: 9.538641885e+19 eint added: 0 ener added: 5.303047905e+43 Castro::do_new_sources() time = 0.14023339 on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.495038611e+26 ener added: 1.500090463e+26 ...estimated hydro-limited timestep at level 1: 6.564309766e-07 ...hydro-limited CFL timestep constrained at (i,j,k) = (121,90,90) Castro::estTimeStep (hydro-limited) at level 1: estdt = 6.564309766e-07  Subcycle completed Subcycling complete Zones advanced per microsecond from level 0 to level 1: 0.02307708085 [Level 0 step 3] Advanced 884736 cells [Level 1 step 3] ADVANCE at time 0.005426075163 with dt = 2.778748814e-09 Advance at this level has already been completed. [Level 1 step 3] Advanced 65536 cells Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.283851852e+24 ener added: 6.674962963e+24 Updating the rotational period from 30.52749735 s to 30.52751553 s. [STEP 3] Coarse TimeStep time: 47.38416933 [STEP 3] FAB kilobyte spread across MPI nodes: [833440 ... 934735] STEP = 3 TIME = 0.005426077942 DT = 2.778748814e-09 ...estimated hydro-limited timestep at level 0: 0.0008130329735 ...hydro-limited CFL timestep constrained at (i,j,k) = (37,49,42) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.0008130329735 ...estimated hydro-limited timestep at level 1: 6.564309766e-07 ...hydro-limited CFL timestep constrained at (i,j,k) = (121,90,90) Castro::estTimeStep (hydro-limited) at level 1: estdt = 6.564309766e-07 Castro::compute_new_dt : limiting dt at level 0 ... new dt computed: 0.0008130329735 ... but limiting to: 3.056623696e-09 = 1.1 * 2.778748814e-09 Castro::compute_new_dt : limiting dt at level 1 ... new dt computed: 6.564309766e-07 ... but limiting to: 3.056623696e-09 = 1.1 * 2.778748814e-09 Now regridding at level lbase = 0 grid_places() time: 0.031934916 new finest: 1 Castro::numpts_1d at level 1 is 340 Gravity::fill_multipole_BCs() time = 0.409385671 Gravity::multilevel_solve_for_new_phi() time = 0.540210844 STEP = 3 TIME = 0.005426077942 : REGRID with lbase = 0 Level 1 4 grids 65536 cells 0.9259259259 % of domain smallest grid: 32 x 32 x 16 biggest grid: 32 x 32 x 16 ...estimated hydro-limited timestep at level 0: 0.0008130329735 ...hydro-limited CFL timestep constrained at (i,j,k) = (37,49,42) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.0008130329735 ...estimated hydro-limited timestep at level 1: 6.564309766e-07 ...hydro-limited CFL timestep constrained at (i,j,k) = (121,90,90) Castro::estTimeStep (hydro-limited) at level 1: estdt = 6.564309766e-07 [Level 0 step 4] ADVANCE at time 0.005426077942 with dt = 3.056623696e-09  Beginning subcycle 1 starting at time 0.005426077942 with dt = 3.056623695e-09 Estimated number of subcycles remaining: 1 Castro::construct_old_gravity() time = 0.92888448 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: 6.694655893e+24 ymom added: 5.755038681e+26 zmom added: -3.471615689e+23 eint added: 0 ener added: -2.130751807e+45 Castro::do_old_sources() time = 2.084383934 on level 0 ... Entering construct_ctu_hydro_source() on level 0 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 21.64894574 on level 0 Castro::construct_old_gravity() time = 0.080990796 on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -2.150425501e+40 ener added: 1.373884257e+36 Contributions to the state from old-time sources: mass added: 0 xmom added: 9.203810089e+22 ymom added: 1.239621202e+22 zmom added: 8.023901326e+21 eint added: 0 ener added: -2.520808161e+45 Castro::do_old_sources() time = 0.193814721 on level 1 ... Entering construct_ctu_hydro_source() on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.318881554e+48 ener added: 1.675623034e+42 ... Leaving construct_ctu_hydro_source() on level 1 Castro::construct_ctu_hydro_source() time = 2.855559479 on level 1 Castro::reflux() at level 0 : time = 2.290938893 ... new-time composite Poisson gravity solve from level 0 to level 1 Gravity::fill_multipole_BCs() time = 0.581292015 Gravity::multilevel_solve_for_new_phi() time = 1.133942925 Castro::construct_new_gravity() time = 1.248612599 on level 0 Contributions to the state from new-time sources: mass added: 0 xmom added: -2.405469507e+21 ymom added: -9.735672131e+21 zmom added: -7.437548959e+18 eint added: 0 ener added: -3.723007452e+41 Castro::do_new_sources() time = 1.968822289 on level 0 ...estimated hydro-limited timestep at level 0: 0.000813036171 ...hydro-limited CFL timestep constrained at (i,j,k) = (37,49,42) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.000813036171 Castro::construct_new_gravity() time = 0.030137785 on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 3.983601425e+43 ener added: 1.617265925e+38 Contributions to the state from new-time sources: mass added: 0 xmom added: 3.521439958e+18 ymom added: -2.667031699e+19 zmom added: -1.338680215e+19 eint added: 0 ener added: 4.042918443e+43 Castro::do_new_sources() time = 0.144179011 on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.851958796e+25 ener added: 2.794684722e+25 ...estimated hydro-limited timestep at level 1: 6.863488975e-07 ...hydro-limited CFL timestep constrained at (i,j,k) = (121,90,90) Castro::estTimeStep (hydro-limited) at level 1: estdt = 6.863488975e-07  Subcycle completed Subcycling complete Zones advanced per microsecond from level 0 to level 1: 0.02282630291 [Level 0 step 4] Advanced 884736 cells [Level 1 step 4] ADVANCE at time 0.005426077942 with dt = 3.056623696e-09 Advance at this level has already been completed. [Level 1 step 4] Advanced 65536 cells Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -3.365462963e+23 ener added: 4.62962963e+19 Updating the rotational period from 30.52751553 s to 30.52751552 s. [STEP 4] Coarse TimeStep time: 47.80455623 [STEP 4] FAB kilobyte spread across MPI nodes: [833440 ... 934735] STEP = 4 TIME = 0.005426080998 DT = 3.056623696e-09 ...estimated hydro-limited timestep at level 0: 0.0002506596503 ...hydro-limited CFL timestep constrained at (i,j,k) = (57,41,46) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.0002506596503 ...estimated hydro-limited timestep at level 1: 6.863488975e-07 ...hydro-limited CFL timestep constrained at (i,j,k) = (121,90,90) Castro::estTimeStep (hydro-limited) at level 1: estdt = 6.863488975e-07 Castro::compute_new_dt : limiting dt at level 0 ... new dt computed: 0.0002506596503 ... but limiting to: 3.362286065e-09 = 1.1 * 3.056623696e-09 Castro::compute_new_dt : limiting dt at level 1 ... new dt computed: 6.863488975e-07 ... but limiting to: 3.362286065e-09 = 1.1 * 3.056623696e-09 Now regridding at level lbase = 0 grid_places() time: 0.032820052 new finest: 1 Castro::numpts_1d at level 1 is 340 Gravity::fill_multipole_BCs() time = 0.412781294 Gravity::multilevel_solve_for_new_phi() time = 0.54809848 STEP = 4 TIME = 0.005426080998 : REGRID with lbase = 0 Level 1 4 grids 65536 cells 0.9259259259 % of domain smallest grid: 32 x 32 x 16 biggest grid: 32 x 32 x 16 ...estimated hydro-limited timestep at level 0: 0.0002506596503 ...hydro-limited CFL timestep constrained at (i,j,k) = (57,41,46) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.0002506596503 ...estimated hydro-limited timestep at level 1: 6.863488975e-07 ...hydro-limited CFL timestep constrained at (i,j,k) = (121,90,90) Castro::estTimeStep (hydro-limited) at level 1: estdt = 6.863488975e-07 [Level 0 step 5] ADVANCE at time 0.005426080998 with dt = 3.362286065e-09  Beginning subcycle 1 starting at time 0.005426080998 with dt = 3.362286065e-09 Estimated number of subcycles remaining: 1 Castro::construct_old_gravity() time = 0.76660303 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: 7.366517551e+24 ymom added: 6.330531998e+26 zmom added: -3.819051789e+23 eint added: 0 ener added: -2.343742535e+45 Castro::do_old_sources() time = 2.100421477 on level 0 ... Entering construct_ctu_hydro_source() on level 0 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 21.46414375 on level 0 Castro::construct_old_gravity() time = 0.068067397 on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -2.560109508e+40 ener added: 1.662186972e+36 Contributions to the state from old-time sources: mass added: 0 xmom added: 1.013699038e+23 ymom added: 1.829917012e+22 zmom added: 8.797223426e+21 eint added: 0 ener added: -2.728291723e+45 Castro::do_old_sources() time = 0.211818603 on level 1 ... Entering construct_ctu_hydro_source() on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.05099824e+48 ener added: 1.448876286e+42 ... Leaving construct_ctu_hydro_source() on level 1 Castro::construct_ctu_hydro_source() time = 2.885581485 on level 1 Castro::reflux() at level 0 : time = 2.422013428 ... new-time composite Poisson gravity solve from level 0 to level 1 Gravity::fill_multipole_BCs() time = 0.791457641 Gravity::multilevel_solve_for_new_phi() time = 1.381133492 Castro::construct_new_gravity() time = 1.497034785 on level 0 Contributions to the state from new-time sources: mass added: 0 xmom added: -2.798425632e+21 ymom added: -1.178643014e+22 zmom added: -5.143756067e+18 eint added: 0 ener added: -4.148475158e+41 Castro::do_new_sources() time = 2.113301696 on level 0 ...estimated hydro-limited timestep at level 0: 0.0006828000901 ...hydro-limited CFL timestep constrained at (i,j,k) = (57,41,46) Castro::estTimeStep (hydro-limited) at level 0: estdt = 0.0006828000901 Castro::construct_new_gravity() time = 0.032567028 on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 3.326552265e+43 ener added: 1.577214716e+38 Contributions to the state from new-time sources: mass added: 0 xmom added: 2.01306964e+18 ymom added: -2.459866116e+19 zmom added: -1.244730808e+19 eint added: 0 ener added: 3.390623254e+43 Castro::do_new_sources() time = 0.149106734 on level 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 4.005277778e+24 ener added: 3.882013889e+24 ...estimated hydro-limited timestep at level 1: 7.169281148e-07 ...hydro-limited CFL timestep constrained at (i,j,k) = (121,90,90) Castro::estTimeStep (hydro-limited) at level 1: estdt = 7.169281148e-07  Subcycle completed Subcycling complete Zones advanced per microsecond from level 0 to level 1: 0.02267348736 [Level 0 step 5] Advanced 884736 cells [Level 1 step 5] ADVANCE at time 0.005426080998 with dt = 3.362286065e-09 Advance at this level has already been completed. [Level 1 step 5] Advanced 65536 cells Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -2.814813006e+22 ener added: 2.962962963e+20 Updating the rotational period from 30.52751552 s to 30.52751551 s. [STEP 5] Coarse TimeStep time: 48.21134807 [STEP 5] FAB kilobyte spread across MPI nodes: [833440 ... 934735] STEP = 5 TIME = 0.005426084361 DT = 3.362286065e-09 PLOTFILE: file = wdmerger_3D_plt00005 Write plotfile time = 0.566007806 seconds Ending run at 06:50:33 UTC on 2023-10-13. Run time = 245.3335691 Run time without initialization = 232.509689 Average number of zones advanced per microsecond: 0.021 Average number of zones advanced per microsecond per rank: 0.005 TinyProfiler total time across processes [min...avg...max]: 245.3 ... 245.3 ... 245.3 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 9 95.9178 104.1255 107.5382 43.83% Castro::computeTemp() 88 15 16.9 17.68 7.21% amrex::Copy() 736 12.89 14.41 15.11 6.16% Castro::normalize_species() 87 10.1 11.83 12.89 5.25% check_for_negative_density() 9 0.6786 4.134 12.43 5.07% MultiFab::AddProduct() 4205 10.44 11.68 12.18 4.97% FabArray::ParallelCopy_finish() 830 1.578 4.283 10.88 4.43% Castro::reset_internal_energy(Fab) 1476 8.446 9.386 9.714 3.96% Castro::estTimeStep() 26 4.044 5.777 9.563 3.90% FabArray::setVal() 5078 7.552 8.422 8.99 3.66% MultiFab::Subtract() 175 6.636 7.433 7.778 3.17% Gravity::fill_multipole_BCs() 11 4.371 4.787 5.571 2.27% Castro::reset_internal_energy(MultiFab) 88 3.546 4.388 5.435 2.22% FabArray::sum() 4205 3.5 3.908 4.076 1.66% FluxRegister::Reflux() 24 2.993 3.344 3.508 1.43% Castro::enforce_min_density() 87 1.638 2.062 2.696 1.10% FabArray::ParallelCopy_nowait() 830 2.051 2.274 2.673 1.09% MLPoisson::Fsmooth() 1376 2.157 2.423 2.561 1.04% Castro::wd_update() 6 0.5184 0.9752 1.939 0.79% Castro::reflux() 4 1.513 1.696 1.788 0.73% Castro::FluxRegCrseInit 9 1.488 1.692 1.773 0.72% FillPatchInterp(MF) 30 1.208 1.314 1.421 0.58% Castro::construct_old_gravity_source() 18 0.8988 1.099 1.263 0.51% Amr::InitAmr() 1 1.25 1.25 1.25 0.51% amrex::Add() 143 0.9959 1.116 1.194 0.49% Castro::initMFs() 5 1.108 1.142 1.164 0.47% FillBoundary_finish() 1561 0.4248 0.6463 1.091 0.44% MLPoisson::Fapply() 326 0.9089 1.008 1.056 0.43% Castro::problem_post_timestep() 10 0.7775 0.8808 1.043 0.42% Castro::construct_new_gravity_source() 9 0.741 0.8374 0.8799 0.36% FillPatchIterator::Initialize 62 0.7595 0.8435 0.8797 0.36% FabArray::Saxpy() 125 0.6942 0.7781 0.8141 0.33% Amr::writePlotFile() 2 0.333 0.5989 0.7892 0.32% VisMF::Write(FabArray) 3 0.6205 0.7127 0.7744 0.32% Castro::volInBoundary() 42 0.6212 0.6902 0.7707 0.31% StateData::FillBoundary(geom) 50 0.6213 0.6304 0.6406 0.26% Castro::construct_new_source() 63 0.4919 0.5379 0.563 0.23% MLCellLinOp::applyBC() 1750 0.3724 0.4345 0.4726 0.19% FillBoundary_nowait() 1561 0.3321 0.3907 0.4533 0.18% MultiFab::contains_nan() 18 0.3962 0.4333 0.4522 0.18% MultiFab::max() 18 0.02724 0.1509 0.435 0.18% FabArray::Xpay() 291 0.3413 0.3823 0.4028 0.16% Castro::initialize_do_advance() 9 0.2352 0.2661 0.2832 0.12% FabArray::setDomainBndry() 92 0.1987 0.227 0.2467 0.10% Castro::construct_old_source() 63 0.1885 0.2022 0.2111 0.09% amrex::average_down 233 0.1928 0.2024 0.2081 0.08% Castro::fill_hybrid_hydro_source() 27 0.1669 0.1861 0.1943 0.08% MLCellLinOp::defineAuxData() 11 0.1537 0.1677 0.1744 0.07% MLPoisson::FFlux() 150 0.004573 0.04201 0.154 0.06% MLCellLinOp::compGrad() 18 0.1293 0.1439 0.1495 0.06% Castro::initData() 1 0.1307 0.1399 0.1472 0.06% FabArray::mult() 56 0.1194 0.1351 0.1456 0.06% Castro::swap_state_time_levels() 9 0.09331 0.1183 0.1338 0.05% MLMG::addInterpCorrection() 155 0.1083 0.1208 0.1281 0.05% StateData::define() 25 0.1028 0.116 0.1247 0.05% Castro::evaluate_source_change() 193 0.1066 0.1166 0.1236 0.05% Castro::initialize_advance() 9 0.07909 0.08888 0.1036 0.04% Gravity::get_new_grav_vector() 10 0.08691 0.0968 0.1017 0.04% Gravity::get_old_grav_vector() 9 0.08523 0.0927 0.09602 0.04% MLMG::prepareForSolve() 10 0.07027 0.078 0.08825 0.04% MLCellLinOp::updateSolBC() 25 0.06684 0.06752 0.06796 0.03% MLMG::solve() 10 0.001472 0.04273 0.06734 0.03% AmrLevel::FillPatch() 62 0.06018 0.06367 0.06593 0.03% FabArray::norminf(mask) 21 0.05513 0.06208 0.06526 0.03% Castro::apply_tagging_restrictions() 6 0.04795 0.05413 0.05822 0.02% AmrMesh::MakeNewGrids() 6 0.04741 0.04903 0.05338 0.02% Gravity::multilevel_solve_for_new_phi() 10 0.03446 0.04131 0.04918 0.02% MLCellLinOp::updateCorBC() 17 0.04556 0.04585 0.04596 0.02% BndryData::define() 27 0.04352 0.0448 0.04579 0.02% FabArray::norminf() 117 0.03911 0.04092 0.04245 0.02% MLMG::ResNormInf() 54 0.004035 0.02164 0.04143 0.02% Castro::do_new_sources() 9 0.0005473 0.009952 0.03576 0.01% MLMG::MLResNormInf() 15 9.9e-05 0.02432 0.0353 0.01% MLCellLinOp::reflux() 30 0.004094 0.01154 0.03354 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 622 0.01844 0.02275 0.02893 0.01% MLCellLinOp::defineBC() 11 0.02596 0.0262 0.02652 0.01% MultiFab::Multiply() 8 0.02182 0.02438 0.0261 0.01% StateDataPhysBCFunct::() 95 0.002031 0.01475 0.02606 0.01% Castro::FluxRegFineAdd() 9 0.02558 0.02564 0.02574 0.01% MLLinOp::makeSubCommunicator() 11 0.0004699 0.008624 0.02012 0.01% Castro::build_fine_mask() 44 0.01804 0.01928 0.02007 0.01% Castro::do_old_sources() 9 0.0006151 0.005555 0.01991 0.01% MLMG::interpCorrection_1 17 0.01808 0.01859 0.0193 0.01% FabArrayBase::CPC::define() 416 0.01395 0.01561 0.01924 0.01% DistributionMapping::LeastUsedCPUs() 5 0.0001491 0.007358 0.01864 0.01% Castro::construct_old_gravity() 9 0.0003687 0.005133 0.01836 0.01% Castro::Castro() 5 0.01484 0.01653 0.01726 0.01% Castro::buildMetrics() 5 0.01418 0.01602 0.01684 0.01% Castro::linear_to_hybrid_momentum() 1 0.01516 0.01583 0.01622 0.01% FillPatchTwoLevels 31 0.01167 0.01367 0.0151 0.01% Gravity::actual_multilevel_solve() 10 0.009924 0.0117 0.01452 0.01% MLCellLinOp::prepareForSolve() 11 0.01165 0.01191 0.0122 0.00% Castro::construct_new_gravity() 9 0.0004369 0.003728 0.01204 0.00% Castro::enforce_consistent_e() 1 0.009507 0.01042 0.01086 0.00% FabArrayBase::getCPC() 1440 0.009892 0.01036 0.01075 0.00% Gravity::update_max_rhs() 10 0.007531 0.008806 0.01048 0.00% AMRErrorTag::operator() 6 0.008596 0.009721 0.0104 0.00% Gravity::actual_solve_with_mlmg() 11 0.009247 0.009751 0.01019 0.00% Castro::finalize_do_advance() 9 0.006034 0.007071 0.008398 0.00% Castro::apply_problem_tags() 6 0.006575 0.007342 0.007623 0.00% AmrLevel::derive() 6 0.004661 0.00529 0.006117 0.00% average_down_faces 48 0.005725 0.005825 0.006026 0.00% MLCellLinOp::setLevelBC() 19 0.003784 0.004622 0.005388 0.00% MLCGSolver::bicgstab 7 0 0.001334 0.005335 0.00% FabArray::FillBoundary() 1561 0.004424 0.004743 0.005181 0.00% MLCellLinOp::smooth() 688 0.004056 0.004248 0.00439 0.00% FabArrayBase::getFB() 1561 0.003238 0.003432 0.003866 0.00% amrex::Dot() 115 0 0.0009072 0.003629 0.00% FabArrayBase::FB::FB() 19 0.002625 0.00297 0.003293 0.00% FabArray::ParallelCopy() 830 0.003078 0.003189 0.003279 0.00% Gravity::test_composite_phi() 1 0.00215 0.002413 0.00258 0.00% FabArray::LinComb() 57 0 0.0004904 0.001962 0.00% main() 1 0.001454 0.001619 0.001805 0.00% MLCellLinOp::apply() 326 0.0009208 0.001143 0.00168 0.00% MLLinOp::defineGrids() 11 0.00133 0.001483 0.001617 0.00% Amr::regrid() 5 0.0009731 0.001198 0.001573 0.00% Castro::clean_state() 87 0.001235 0.001323 0.001471 0.00% Amr::coarseTimeStep() 5 0.000545 0.001123 0.001317 0.00% TagBoxArray::collate() 6 0.0003083 0.0009875 0.001309 0.00% AmrLevel::FillCoarsePatch() 5 0.0009522 0.001038 0.001238 0.00% MLCGSolver::ParallelAllReduce 153 0 0.0002879 0.001151 0.00% Gravity::average_fine_ec_onto_crse_ec() 8 0.0009173 0.0009314 0.0009409 0.00% OwnerMask() 6 0.0007333 0.0008402 0.0008969 0.00% MLCellLinOp::correctionResidual() 196 0.0007084 0.0007635 0.000847 0.00% MLMG::mgVcycle() 65 0.0006012 0.0006291 0.0006437 0.00% Castro::do_advance_ctu() 5 0.0005805 0.0006039 0.0006282 0.00% FillPatchSingleLevel 92 0.0005658 0.0005774 0.0005927 0.00% Castro::finalize_advance() 9 0.000396 0.000515 0.0005898 0.00% Amr::grid_places() 6 0.0004461 0.0004891 0.0005366 0.00% FillPatchIterator::FillFromLevel0() 31 0.0004154 0.0004431 0.000489 0.00% TagBoxArray::mapPRD 6 0.0004427 0.0004599 0.0004776 0.00% MLMG:computeResOfCorrection() 155 0.000407 0.0004294 0.0004528 0.00% Castro::subcycle_advance_ctu() 5 7.2e-05 0.0001786 0.0004334 0.00% MLMG::actualBottomSolve() 31 4.025e-05 0.0001313 0.0004008 0.00% MLMG::oneIter() 31 0.0003585 0.000367 0.000378 0.00% MLCellLinOp::solutionResidual() 72 0.0003582 0.0003673 0.0003745 0.00% Castro::computeNewDt() 8 0.0001105 0.0001765 0.0003461 0.00% FillPatchIterator::FillFromTwoLevels() 31 0.0003038 0.0003166 0.0003246 0.00% MLMG::mgVcycle_down::0 31 0.0002791 0.0002964 0.0003071 0.00% Castro::enforce_speed_limit() 87 0.0002585 0.0002633 0.0002741 0.00% Gravity::solve_phi_with_mlmg() 11 0.0002389 0.0002571 0.0002666 0.00% MLMG::mgVcycle_down::1 31 0.0002448 0.0002546 0.0002607 0.00% Amr::defBaseLevel() 1 0.0002468 0.0002498 0.0002555 0.00% Amr::timeStep() 9 0.0001858 0.0002115 0.0002417 0.00% MLMG::mgVcycle_bottom 65 0.000217 0.0002207 0.0002239 0.00% MLMG::mgVcycle_down::2 31 0.0002028 0.0002117 0.0002211 0.00% makeSFC 33 0.0001726 0.0001945 0.0002197 0.00% MLMG::mgVcycle_down::3 31 0.0002063 0.0002128 0.0002189 0.00% Castro::post_regrid() 10 0.000182 0.0001989 0.0002128 0.00% AmrLevel::AmrLevel(dm) 5 0.0001909 0.0002024 0.0002082 0.00% MLMG::mgVcycle_up::0 31 0.0001887 0.0001945 0.0002021 0.00% Castro::check_for_nan() 18 0.0001347 0.0001673 0.0001974 0.00% FabArrayBase::TheFPinfo() 25 0.0001673 0.0001784 0.0001943 0.00% MLMG::mgVcycle_down::4 31 0.0001831 0.000188 0.000194 0.00% Amr::FinalizeInit() 1 3.651e-05 7.29e-05 0.0001749 0.00% Castro::react_state() 18 0.0001384 0.0001526 0.0001662 0.00% MLMG::mgVcycle_up::2 31 0.0001392 0.0001452 0.0001509 0.00% MLMG::mgVcycle_up::1 31 0.0001441 0.0001472 0.0001508 0.00% MLMG::computeResWithCrseSolFineCor() 17 0.0001377 0.0001426 0.0001467 0.00% MLMG::mgVcycle_up::4 31 0.0001184 0.0001254 0.0001393 0.00% MLMG::mgVcycle_up::3 31 0.0001301 0.0001336 0.0001375 0.00% FPinfo::FPinfo() 3 0.0001153 0.0001287 0.0001364 0.00% MLLinOp::makeAgglomeratedDMap 11 0.0001142 0.0001211 0.0001338 0.00% DistributionMapping::Distribute() 35 0.0001108 0.0001229 0.000132 0.00% Castro::advance() 5 9.334e-05 0.0001122 0.0001265 0.00% FabArray::RecvLayoutMask() 8 9.569e-05 0.0001058 0.0001262 0.00% Castro::expand_state() 13 0.0001079 0.000116 0.0001227 0.00% Castro::apply_source_to_state() 18 9.465e-05 0.0001071 0.0001198 0.00% Amr::writeSmallPlotFile() 1 9.477e-05 0.0001059 0.0001166 0.00% MLPoisson::define() 11 0.0001042 0.0001102 0.0001163 0.00% Castro::errorEst() 6 0.0001055 0.0001106 0.0001161 0.00% Gravity::swapTimeLevels() 9 8.482e-05 9.269e-05 0.0001009 0.00% ClusterList::intersect() 1 0 2.435e-05 9.741e-05 0.00% MLMG::computeResidual() 31 8.41e-05 8.894e-05 9.528e-05 0.00% Castro::create_source_corrector() 9 5.52e-05 7.543e-05 8.883e-05 0.00% MLLinOp::define() 11 6.517e-05 7.292e-05 8.827e-05 0.00% Castro::post_timestep() 9 7.441e-05 8.004e-05 8.783e-05 0.00% Castro::avgDown(state_indx) 20 6.975e-05 7.698e-05 8.683e-05 0.00% MLMG::computeMLResidual() 15 7.835e-05 8.093e-05 8.31e-05 0.00% knapsack() 3 6.554e-05 6.924e-05 7.508e-05 0.00% MLMG::miniCycle() 34 7.07e-05 7.208e-05 7.331e-05 0.00% Castro::retry_advance_ctu() 5 2.33e-05 3.829e-05 7.222e-05 0.00% MLMG::computeResWithCrseCorFineCor() 17 6.575e-05 6.721e-05 6.91e-05 0.00% MLMG::getGradSolution() 10 5.744e-05 6.26e-05 6.612e-05 0.00% amrex::average_down_w_geom 20 5.084e-05 5.424e-05 5.76e-05 0.00% Castro::construct_new_hybrid_source() 9 4.918e-05 5.143e-05 5.363e-05 0.00% Castro::init(old) 3 4.427e-05 4.793e-05 5.259e-05 0.00% Amr::initSubcycle() 1 5.109e-05 5.165e-05 5.212e-05 0.00% MLPoisson::prepareForSolve() 11 4.766e-05 4.978e-05 5.188e-05 0.00% DistributionMapping::SFCProcessorMapDoIt() 2 4.696e-05 4.83e-05 4.957e-05 0.00% Gravity::set_mass_offset() 6 3.503e-05 4.118e-05 4.911e-05 0.00% AmrMesh-cluster 1 0 1.125e-05 4.499e-05 0.00% Amr::bldFineLevels() 1 1.901e-05 2.792e-05 4.382e-05 0.00% Castro::construct_old_hybrid_source() 9 3.377e-05 3.82e-05 4.104e-05 0.00% Castro::post_init() 1 2.696e-05 3.256e-05 3.95e-05 0.00% MLMG::MLRhsNormInf() 10 3.445e-05 3.525e-05 3.605e-05 0.00% Castro::avgDown() 4 3.03e-05 3.238e-05 3.58e-05 0.00% DistributionMapping::KnapSackDoIt() 3 3.046e-05 3.347e-05 3.501e-05 0.00% Castro::construct_old_rotation_source() 9 2.227e-05 2.598e-05 3.046e-05 0.00% BoxList::complementIn 12 2.205e-05 2.488e-05 3.023e-05 0.00% MLMG::compResidual() 1 2.856e-05 2.922e-05 3.012e-05 0.00% BoxList::parallelComplementIn() 12 2.609e-05 2.724e-05 2.878e-05 0.00% Castro::construct_new_rotation_source() 9 2.122e-05 2.419e-05 2.819e-05 0.00% Castro::derive() 6 1.788e-05 2.085e-05 2.79e-05 0.00% Castro::computeInitialDt() 2 1.779e-05 2.054e-05 2.748e-05 0.00% Castro::init() 1 1.68e-05 1.744e-05 1.818e-05 0.00% Amr::InitializeInit() 1 7.231e-06 1.182e-05 1.772e-05 0.00% ClusterList::chop() 1 0 2.8e-06 1.12e-05 0.00% Amr::initialInit() 1 4.722e-06 7.691e-06 1.09e-05 0.00% Amr::init() 1 6.288e-06 6.994e-06 7.732e-06 0.00% knapsack()swap 1 3.386e-06 3.635e-06 3.816e-06 0.00% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 245.3 245.3 245.3 100.00% Amr::coarseTimeStep() 5 231.9 231.9 231.9 94.54% Amr::timeStep() 9 230.6 230.6 230.6 94.00% Castro::advance() 5 202.8 202.8 202.8 82.68% Castro::subcycle_advance_ctu() 5 194 194.4 195.4 79.64% Castro::do_advance_ctu() 5 194 194.4 195.4 79.64% Castro::construct_ctu_hydro_source() 9 118.2 118.3 118.3 48.24% Castro::clean_state() 87 71.9 81.65 85.82 34.98% Castro::computeTemp() 88 45.83 51.7 54.69 22.29% Castro::reset_internal_energy(MultiFab) 88 30.83 34.81 37.01 15.09% Castro::evaluate_source_change() 193 17.85 19.95 20.82 8.49% Castro::enforce_min_density() 87 16.78 19.06 20.49 8.35% Castro::post_regrid() 10 16.14 17.67 18.28 7.45% Amr::regrid() 5 15.93 17.36 17.95 7.32% Castro::initialize_do_advance() 9 14.9 15.39 15.63 6.37% amrex::Copy() 736 12.89 14.41 15.11 6.16% AmrLevel::FillPatch() 62 8.027 9.965 14.73 6.01% FillPatchIterator::Initialize 62 6.29 8.286 13.2 5.38% FabArray::ParallelCopy() 830 3.86 6.622 13.1 5.34% Castro::normalize_species() 87 10.1 11.83 12.89 5.25% Gravity::solve_phi_with_mlmg() 11 11.52 11.94 12.72 5.19% check_for_negative_density() 9 0.6786 4.134 12.43 5.07% Gravity::multilevel_solve_for_new_phi() 10 11.24 11.65 12.42 5.06% Gravity::actual_multilevel_solve() 10 11.19 11.6 12.38 5.05% MultiFab::AddProduct() 4205 10.44 11.68 12.18 4.97% Amr::init() 1 11.57 11.57 11.57 4.72% Castro::do_old_sources() 9 11.32 11.37 11.45 4.67% FillPatchSingleLevel 92 4.008 6.09 11.17 4.55% Castro::do_new_sources() 9 9.34 10.53 11.01 4.49% FabArray::ParallelCopy_finish() 830 1.578 4.283 10.88 4.43% Amr::initialInit() 1 10.49 10.55 10.74 4.38% FillPatchIterator::FillFromLevel0() 31 3.614 5.592 10.53 4.29% Castro::post_timestep() 9 10.36 10.36 10.36 4.22% Castro::reset_internal_energy(Fab) 1476 8.446 9.386 9.714 3.96% Castro::estTimeStep() 26 4.044 5.777 9.563 3.90% Castro::reflux() 4 8.46 9.039 9.335 3.81% FabArray::setVal() 5078 7.552 8.422 8.99 3.66% Castro::initialize_advance() 9 7.452 8.406 8.805 3.59% MultiFab::Subtract() 175 6.636 7.433 7.778 3.17% Castro::construct_new_gravity() 9 6.603 6.992 7.734 3.15% Gravity::actual_solve_with_mlmg() 11 7.049 7.056 7.064 2.88% Amr::FinalizeInit() 1 6.149 6.379 6.759 2.75% Castro::expand_state() 13 5.256 5.719 6.598 2.69% MLMG::solve() 10 6.527 6.54 6.558 2.67% Gravity::fill_multipole_BCs() 11 4.452 4.865 5.642 2.30% MLMG::oneIter() 31 5.305 5.321 5.338 2.18% FluxRegister::Reflux() 24 3.843 4.382 4.686 1.91% MLMG::mgVcycle() 65 4.523 4.577 4.606 1.88% Amr::InitializeInit() 1 3.728 4.173 4.595 1.87% Amr::defBaseLevel() 1 3.728 4.173 4.595 1.87% Castro::construct_old_gravity() 9 0.7337 1.725 4.422 1.80% Gravity::get_old_grav_vector() 9 0.7201 1.72 4.422 1.80% Castro::finalize_do_advance() 9 1.369 2.228 4.303 1.75% Castro::problem_post_timestep() 10 2.812 3.253 4.083 1.66% FabArray::sum() 4205 3.5 3.908 4.076 1.66% Castro::computeNewDt() 8 2.149 2.628 3.894 1.59% MLCellLinOp::smooth() 688 3.454 3.513 3.572 1.46% Castro::initData() 1 2.326 2.732 3.17 1.29% Castro::post_init() 1 2.618 2.765 2.976 1.21% Castro::FluxRegCrseInit 9 1.819 2.148 2.94 1.20% FabArray::ParallelCopy_nowait() 830 2.11 2.336 2.74 1.12% MLPoisson::Fsmooth() 1376 2.157 2.423 2.561 1.04% Castro::wd_update() 6 1.201 1.665 2.56 1.04% MLMG::mgVcycle_down::0 31 2.059 2.104 2.161 0.88% Castro::construct_new_source() 63 1.778 1.985 2.081 0.85% MLCellLinOp::applyBC() 1750 1.27 1.483 1.928 0.79% FillPatchIterator::FillFromTwoLevels() 31 1.519 1.624 1.72 0.70% FillPatchTwoLevels 31 1.518 1.623 1.719 0.70% Amr::writePlotFile() 2 1.396 1.588 1.653 0.67% Amr::bldFineLevels() 1 1.36 1.522 1.636 0.67% FabArray::FillBoundary() 1561 0.7972 1.048 1.555 0.63% MLMG::mgVcycle_up::0 31 1.44 1.462 1.488 0.61% Castro::Castro() 5 1.376 1.404 1.457 0.59% FillPatchInterp(MF) 30 1.208 1.314 1.421 0.58% Castro::computeInitialDt() 2 0.5065 0.9286 1.373 0.56% MLCellLinOp::apply() 326 1.315 1.333 1.359 0.55% Castro::construct_old_gravity_source() 18 0.8988 1.099 1.263 0.51% Castro::construct_old_source() 63 0.9126 1.101 1.262 0.51% Amr::InitAmr() 1 1.25 1.25 1.25 0.51% amrex::Add() 143 0.9959 1.116 1.194 0.49% Castro::initMFs() 5 1.113 1.148 1.17 0.48% FillBoundary_finish() 1561 0.4248 0.6463 1.091 0.44% MLPoisson::Fapply() 326 0.9089 1.008 1.056 0.43% MLCellLinOp::solutionResidual() 72 1.02 1.036 1.048 0.43% Castro::construct_new_gravity_source() 9 0.741 0.8374 0.8799 0.36% Gravity::get_new_grav_vector() 10 0.8222 0.8485 0.8715 0.36% MLCellLinOp::correctionResidual() 196 0.821 0.827 0.8358 0.34% FabArray::Saxpy() 125 0.6942 0.7781 0.8141 0.33% VisMF::Write(FabArray) 3 0.6205 0.7127 0.7744 0.32% Castro::volInBoundary() 42 0.6213 0.6903 0.7707 0.31% MLMG:computeResOfCorrection() 155 0.6654 0.6973 0.7184 0.29% StateDataPhysBCFunct::() 95 0.6234 0.6452 0.6586 0.27% StateData::FillBoundary(geom) 50 0.6213 0.6304 0.6406 0.26% MLMG::computeResWithCrseSolFineCor() 17 0.4648 0.4967 0.578 0.24% MLMG::computeMLResidual() 15 0.4141 0.4514 0.5006 0.20% Gravity::update_max_rhs() 10 0.08739 0.2098 0.4888 0.20% Castro::init() 1 0.479 0.4804 0.4829 0.20% AmrLevel::FillCoarsePatch() 5 0.479 0.4804 0.4828 0.20% FillBoundary_nowait() 1561 0.3381 0.3971 0.4598 0.19% Castro::check_for_nan() 18 0.3964 0.4335 0.4524 0.18% MultiFab::contains_nan() 18 0.3962 0.4333 0.4522 0.18% MultiFab::max() 18 0.02724 0.1509 0.435 0.18% MLMG::mgVcycle_down::1 31 0.3302 0.3624 0.4184 0.17% Gravity::test_composite_phi() 1 0.4166 0.4171 0.4177 0.17% FabArray::Xpay() 291 0.3413 0.3823 0.4028 0.16% Castro::init(old) 3 0.3668 0.3703 0.3761 0.15% MLCellLinOp::reflux() 30 0.1833 0.2321 0.363 0.15% MLMG::computeResidual() 31 0.3587 0.3599 0.3613 0.15% Castro::apply_source_to_state() 18 0.2773 0.3111 0.3239 0.13% MLPoisson::define() 11 0.2897 0.2977 0.3057 0.12% Castro::react_state() 18 0.2575 0.2753 0.3053 0.12% amrex::average_down 233 0.2585 0.2666 0.2756 0.11% FabArray::setDomainBndry() 92 0.1987 0.227 0.2467 0.10% MLMG::prepareForSolve() 10 0.2005 0.2204 0.2348 0.10% MLMG::mgVcycle_up::1 31 0.2207 0.226 0.2298 0.09% MLCellLinOp::defineAuxData() 11 0.1824 0.1982 0.2093 0.09% Amr::grid_places() 6 0.1995 0.1997 0.1998 0.08% AmrMesh::MakeNewGrids() 6 0.1991 0.1992 0.1993 0.08% Castro::fill_hybrid_hydro_source() 27 0.1669 0.1861 0.1943 0.08% MLMG::addInterpCorrection() 155 0.1306 0.172 0.1936 0.08% MLMG::mgVcycle_bottom 65 0.1576 0.1663 0.1918 0.08% MLMG::getGradSolution() 10 0.182 0.1834 0.1846 0.08% MLCellLinOp::compGrad() 18 0.1819 0.1834 0.1845 0.08% MLMG::miniCycle() 34 0.1577 0.1579 0.1581 0.06% MLPoisson::FFlux() 150 0.004573 0.04201 0.154 0.06% FabArray::mult() 56 0.1194 0.1351 0.1456 0.06% Castro::avgDown() 4 0.1165 0.1269 0.1339 0.05% Castro::avgDown(state_indx) 20 0.1165 0.1268 0.1338 0.05% amrex::average_down_w_geom 20 0.1164 0.1268 0.1338 0.05% Castro::swap_state_time_levels() 9 0.09331 0.1183 0.1338 0.05% MLCellLinOp::updateSolBC() 25 0.07755 0.09776 0.133 0.05% Castro::construct_new_hybrid_source() 9 0.1108 0.1242 0.1302 0.05% Castro::errorEst() 6 0.1065 0.1193 0.1271 0.05% AmrLevel::AmrLevel(dm) 5 0.103 0.1162 0.1249 0.05% StateData::define() 25 0.1028 0.116 0.1247 0.05% MLMG::ResNormInf() 54 0.0767 0.09223 0.1045 0.04% MLCellLinOp::defineBC() 11 0.08634 0.08878 0.09172 0.04% MLMG::mgVcycle_down::2 31 0.07705 0.08231 0.09115 0.04% MLMG::computeResWithCrseCorFineCor() 17 0.08688 0.08802 0.0887 0.04% MLMG::mgVcycle_up::3 31 0.02592 0.07044 0.08621 0.04% MLMG::MLResNormInf() 15 0.04452 0.06767 0.07682 0.03% FabArray::norminf(mask) 21 0.05513 0.06208 0.06526 0.03% BndryData::define() 27 0.05897 0.0613 0.06449 0.03% MLMG::interpCorrection_1 17 0.02629 0.04256 0.06432 0.03% Castro::construct_old_hybrid_source() 9 0.05618 0.06192 0.06419 0.03% Castro::apply_tagging_restrictions() 6 0.04795 0.05413 0.05822 0.02% MLMG::mgVcycle_up::2 31 0.05244 0.05334 0.05414 0.02% MLCellLinOp::updateCorBC() 17 0.05112 0.05146 0.05182 0.02% Castro::derive() 6 0.04324 0.048 0.05073 0.02% AmrLevel::derive() 6 0.04322 0.04797 0.0507 0.02% Gravity::swapTimeLevels() 9 0.03511 0.03844 0.04511 0.02% FabArray::norminf() 117 0.03911 0.04092 0.04245 0.02% Castro::buildMetrics() 5 0.03627 0.03939 0.04194 0.02% MLMG::mgVcycle_down::3 31 0.03086 0.03214 0.03544 0.01% MLMG::actualBottomSolve() 31 4.025e-05 0.00855 0.03408 0.01% MLCGSolver::bicgstab 7 0 0.008329 0.03332 0.01% TagBoxArray::mapPRD 6 0.004005 0.01748 0.03329 0.01% FabArrayBase::getCPC() 1440 0.02427 0.02597 0.02967 0.01% MLMG::MLRhsNormInf() 10 0.02537 0.02819 0.02925 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 622 0.01844 0.02275 0.02893 0.01% MultiFab::Multiply() 8 0.02182 0.02438 0.0261 0.01% Castro::FluxRegFineAdd() 9 0.02558 0.02564 0.02574 0.01% MLMG::compResidual() 1 0.0189 0.02088 0.02254 0.01% Gravity::average_fine_ec_onto_crse_ec() 8 0.0148 0.01716 0.02228 0.01% MLLinOp::define() 11 0.002402 0.01061 0.02218 0.01% MLLinOp::defineGrids() 11 0.002337 0.01053 0.02209 0.01% MLLinOp::makeSubCommunicator() 11 0.0004699 0.008624 0.02012 0.01% Castro::build_fine_mask() 44 0.01804 0.01928 0.02007 0.01% FabArrayBase::CPC::define() 416 0.01395 0.01561 0.01924 0.01% MLMG::mgVcycle_down::4 31 0.005334 0.009136 0.0191 0.01% FabArrayBase::TheFPinfo() 25 0.0005364 0.007725 0.01901 0.01% FPinfo::FPinfo() 3 0.0003421 0.007547 0.01883 0.01% DistributionMapping::KnapSackDoIt() 3 0.0002118 0.007418 0.0187 0.01% DistributionMapping::LeastUsedCPUs() 5 0.0001491 0.007358 0.01864 0.01% Castro::linear_to_hybrid_momentum() 1 0.01516 0.01583 0.01622 0.01% MLMG::mgVcycle_up::4 31 0.004397 0.007169 0.01449 0.01% MLPoisson::prepareForSolve() 11 0.0117 0.01196 0.01224 0.00% MLCellLinOp::prepareForSolve() 11 0.01165 0.01191 0.0122 0.00% Castro::enforce_consistent_e() 1 0.009507 0.01042 0.01086 0.00% average_down_faces 48 0.009316 0.009916 0.01063 0.00% AMRErrorTag::operator() 6 0.008596 0.009721 0.0104 0.00% Castro::apply_problem_tags() 6 0.006575 0.007342 0.007623 0.00% FabArrayBase::getFB() 1561 0.005971 0.006402 0.006982 0.00% MLCellLinOp::setLevelBC() 19 0.003784 0.004622 0.005388 0.00% amrex::Dot() 115 0 0.0009072 0.003629 0.00% FabArrayBase::FB::FB() 19 0.002625 0.00297 0.003293 0.00% MLCGSolver::ParallelAllReduce 153 0 0.0007727 0.003091 0.00% FabArray::LinComb() 57 0 0.0004904 0.001962 0.00% TagBoxArray::collate() 6 0.0003083 0.0009875 0.001309 0.00% OwnerMask() 6 0.0007333 0.0008402 0.0008969 0.00% Castro::finalize_advance() 9 0.000396 0.000515 0.0005898 0.00% MLLinOp::makeAgglomeratedDMap 11 0.000385 0.0004264 0.0004719 0.00% makeSFC 33 0.0002708 0.0003052 0.0003382 0.00% Castro::enforce_speed_limit() 87 0.0002585 0.0002633 0.0002741 0.00% AmrMesh-cluster 1 0 3.84e-05 0.0001536 0.00% DistributionMapping::Distribute() 35 0.0001108 0.0001229 0.000132 0.00% FabArray::RecvLayoutMask() 8 9.569e-05 0.0001058 0.0001262 0.00% DistributionMapping::SFCProcessorMapDoIt() 2 9.129e-05 0.0001069 0.000118 0.00% Amr::writeSmallPlotFile() 1 9.477e-05 0.0001059 0.0001166 0.00% ClusterList::intersect() 1 0 2.435e-05 9.741e-05 0.00% Castro::create_source_corrector() 9 5.52e-05 7.543e-05 8.883e-05 0.00% knapsack() 3 6.893e-05 7.287e-05 7.89e-05 0.00% Castro::retry_advance_ctu() 5 2.33e-05 3.829e-05 7.222e-05 0.00% BoxList::parallelComplementIn() 12 4.822e-05 5.211e-05 5.773e-05 0.00% Amr::initSubcycle() 1 5.109e-05 5.165e-05 5.212e-05 0.00% Gravity::set_mass_offset() 6 3.503e-05 4.118e-05 4.911e-05 0.00% Castro::construct_old_rotation_source() 9 2.227e-05 2.598e-05 3.046e-05 0.00% BoxList::complementIn 12 2.205e-05 2.488e-05 3.023e-05 0.00% Castro::construct_new_rotation_source() 9 2.122e-05 2.419e-05 2.819e-05 0.00% ClusterList::chop() 1 0 2.8e-06 1.12e-05 0.00% knapsack()swap 1 3.386e-06 3.635e-06 3.816e-06 0.00% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [2] Pinned Memory Usage: --------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max --------------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 4 4 0 B 0 B 0 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (23.10-13-gecaa14456e3f) finalized