MPI initialized with 4 MPI processes MPI initialized with thread support level 3 OMP initialized with 4 OMP threads AMReX (23.07-13-gb41f85ab36ac) initialized Starting run at 19:16:14 UTC on 2023-07-08. Successfully read inputs file ... Castro git describe: 23.07-24-g8ffe34660 AMReX git describe: 23.07-13-gb41f85ab3 Microphysics git describe: 23.07-19-g4bea489a 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.054341025 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.085479651 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.394967769 Gravity::multilevel_solve_for_new_phi() time = 1.748493126 Gravity::fill_multipole_BCs() time = 0.362083211 ... 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 = 0.667840335 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 [1m[32m Beginning subcycle 1 starting at time 0 with dt = 0.005426072637[0m Estimated number of subcycles remaining: 1 Castro::construct_old_gravity() time = 0.822174867 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.482623477 ... Entering construct_ctu_hydro_source() 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() Castro::construct_ctu_hydro_source() time = 20.70372343 ... new-time composite Poisson gravity solve from level 0 to level 0 Gravity::fill_multipole_BCs() time = 0.655753968 Gravity::multilevel_solve_for_new_phi() time = 1.688452085 Castro::construct_new_gravity() time = 1.80497966 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.451207933e+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.375502272 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 [1m[32m Subcycle completed[0m Subcycling complete Zones advanced per microsecond at this level: 0.0247787717 [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: 37.92963669 [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.032731138 new finest: 1 Castro::numpts_1d at level 1 is 340 Gravity::fill_multipole_BCs() time = 0.401655752 Gravity::multilevel_solve_for_new_phi() time = 1.91372847 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 [1m[32m Beginning subcycle 1 starting at time 0.005426072637 with dt = 2.526135286e-09[0m Estimated number of subcycles remaining: 1 Castro::construct_old_gravity() time = 0.822573681 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.15757758 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 21.29807673 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.043713898 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.575026313e+22 ymom added: 1.534769107e+22 zmom added: 6.442093022e+21 eint added: 0 ener added: -2.586789965e+47 Castro::do_old_sources() time = 0.210043023 ... Entering construct_ctu_hydro_source() 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() Castro::construct_ctu_hydro_source() time = 2.840443681 Castro::reflux() at level 0 : time = 2.326072588 ... new-time composite Poisson gravity solve from level 0 to level 1 Gravity::fill_multipole_BCs() time = 0.57323856 Gravity::multilevel_solve_for_new_phi() time = 1.54819991 Castro::construct_new_gravity() time = 1.665844311 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 = 2.133008906 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.029213964 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.112156155e+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.148002301 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 [1m[32m Subcycle completed[0m Subcycling complete Zones advanced per microsecond from level 0 to level 1: 0.02263485551 [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.1446632 [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.031775967 new finest: 1 Castro::numpts_1d at level 1 is 340 Gravity::fill_multipole_BCs() time = 0.403441501 Gravity::multilevel_solve_for_new_phi() time = 0.533406829 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 [1m[32m Beginning subcycle 1 starting at time 0.005426075163 with dt = 2.778748814e-09[0m Estimated number of subcycles remaining: 1 Castro::construct_old_gravity() time = 0.796110039 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.10782819 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 21.42028495 Castro::construct_old_gravity() time = 0.08035403 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.444240294e+22 ymom added: 4.132070673e+21 zmom added: 7.10314939e+21 eint added: 0 ener added: -2.344787748e+45 Castro::do_old_sources() time = 0.230942135 ... Entering construct_ctu_hydro_source() 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() Castro::construct_ctu_hydro_source() time = 2.827351232 Castro::reflux() at level 0 : time = 2.375965254 ... new-time composite Poisson gravity solve from level 0 to level 1 Gravity::fill_multipole_BCs() time = 0.633270325 Gravity::multilevel_solve_for_new_phi() time = 1.186860443 Castro::construct_new_gravity() time = 1.303306625 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 = 2.099601079 ...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.029913307 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.605965175e+17 ymom added: 2.17388754e+19 zmom added: 9.538642311e+19 eint added: 0 ener added: 5.303047905e+43 Castro::do_new_sources() time = 0.151043907 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 [1m[32m Subcycle completed[0m Subcycling complete Zones advanced per microsecond from level 0 to level 1: 0.02284421498 [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.82220562 [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.03168022 new finest: 1 Castro::numpts_1d at level 1 is 340 Gravity::fill_multipole_BCs() time = 0.402437283 Gravity::multilevel_solve_for_new_phi() time = 0.531447917 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 [1m[32m Beginning subcycle 1 starting at time 0.005426077942 with dt = 3.056623695e-09[0m Estimated number of subcycles remaining: 1 Castro::construct_old_gravity() time = 0.842792722 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.141650124 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 21.48797085 Castro::construct_old_gravity() time = 0.038641807 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.202438518e+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.207236337 ... Entering construct_ctu_hydro_source() 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() Castro::construct_ctu_hydro_source() time = 2.860943536 Castro::reflux() at level 0 : time = 2.368928527 ... new-time composite Poisson gravity solve from level 0 to level 1 Gravity::fill_multipole_BCs() time = 0.611697053 Gravity::multilevel_solve_for_new_phi() time = 1.196980771 Castro::construct_new_gravity() time = 1.312770244 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 = 2.141029472 ...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.032927969 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.521555441e+18 ymom added: -2.667031699e+19 zmom added: -1.338680256e+19 eint added: 0 ener added: 4.042918443e+43 Castro::do_new_sources() time = 0.193991412 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 [1m[32m Subcycle completed[0m Subcycling complete Zones advanced per microsecond from level 0 to level 1: 0.02278038609 [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.90087505 [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.032639137 new finest: 1 Castro::numpts_1d at level 1 is 340 Gravity::fill_multipole_BCs() time = 0.407219914 Gravity::multilevel_solve_for_new_phi() time = 0.539990764 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 [1m[32m Beginning subcycle 1 starting at time 0.005426080998 with dt = 3.362286065e-09[0m Estimated number of subcycles remaining: 1 Castro::construct_old_gravity() time = 0.841016309 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.186750776 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 21.53039053 Castro::construct_old_gravity() time = 0.040488603 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.013350437e+23 ymom added: 1.770887431e+22 zmom added: 8.796935195e+21 eint added: 0 ener added: -2.728291723e+45 Castro::do_old_sources() time = 0.201312521 ... Entering construct_ctu_hydro_source() 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() Castro::construct_ctu_hydro_source() time = 2.786587475 Castro::reflux() at level 0 : time = 2.356188846 ... new-time composite Poisson gravity solve from level 0 to level 1 Gravity::fill_multipole_BCs() time = 0.6532852 Gravity::multilevel_solve_for_new_phi() time = 1.226054569 Castro::construct_new_gravity() time = 1.344106986 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.156995459 ...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.029926106 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.013069657e+18 ymom added: -2.460766836e+19 zmom added: -1.244731042e+19 eint added: 0 ener added: 3.390623254e+43 Castro::do_new_sources() time = 0.153817636 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 [1m[32m Subcycle completed[0m Subcycling complete Zones advanced per microsecond from level 0 to level 1: 0.02271825251 [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.115437 [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.586641268 seconds Ending run at 19:20:19 UTC on 2023-07-08. Run time = 244.4993795 Run time without initialization = 232.5014351 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]: 244.5 ... 244.5 ... 244.5 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 9 92.1761 102.7615 106.9464 43.74% Castro::computeTemp() 88 14.7 16.62 17.37 7.10% check_for_negative_density() 9 0.3131 4.529 15.17 6.21% amrex::Copy() 736 12.86 14.31 15.06 6.16% Castro::normalize_species() 87 10.08 11.41 12.72 5.20% MultiFab::AddProduct() 4205 10.46 11.67 12.23 5.00% FabArray::ParallelCopy_finish() 830 1.481 4.412 10.78 4.41% Castro::reset_internal_energy(Fab) 1476 8.284 9.247 9.601 3.93% Castro::estTimeStep() 26 4.022 5.691 9.534 3.90% FabArray::setVal() 5078 7.347 8.215 8.742 3.58% MultiFab::Subtract() 175 6.646 7.433 7.77 3.18% Castro::reset_internal_energy(MultiFab) 88 3.808 4.757 6 2.45% Gravity::fill_multipole_BCs() 11 4.268 4.64 5.415 2.21% FabArray::sum() 4205 3.446 3.857 4.06 1.66% FluxRegister::Reflux() 24 3.043 3.421 3.568 1.46% Castro::enforce_min_density() 87 1.791 2.346 3.307 1.35% FabArray::ParallelCopy_nowait() 830 2.148 2.321 2.705 1.11% MLPoisson::Fsmooth() 1376 2.067 2.328 2.429 0.99% Castro::wd_update() 6 0.5081 1.016 1.989 0.81% Castro::FluxRegCrseInit 9 1.565 1.794 1.947 0.80% Castro::reflux() 4 1.577 1.775 1.89 0.77% FabArray::Saxpy() 125 1.471 1.649 1.728 0.71% FillPatchInterp(MF) 30 1.217 1.328 1.428 0.58% amrex::Add() 143 1.063 1.194 1.3 0.53% Amr::InitAmr() 1 1.23 1.23 1.23 0.50% Castro::construct_old_gravity_source() 18 0.9225 1.045 1.209 0.49% Castro::initMFs() 5 1.101 1.147 1.172 0.48% FillBoundary_finish() 1561 0.4048 0.6259 1.092 0.45% MLPoisson::Fapply() 326 0.885 0.9854 1.024 0.42% FillPatchIterator::Initialize 62 0.8058 0.8924 0.9387 0.38% Castro::construct_new_gravity_source() 9 0.7489 0.8431 0.913 0.37% Castro::problem_post_timestep() 10 0.7974 0.8475 0.9015 0.37% VisMF::Write(FabArray) 3 0.6387 0.7443 0.8239 0.34% Castro::volInBoundary() 42 0.5746 0.6511 0.7695 0.31% StateData::FillBoundary(geom) 50 0.6105 0.6224 0.639 0.26% Castro::construct_new_source() 63 0.492 0.5447 0.5738 0.23% MultiFab::max() 18 0.01806 0.1569 0.4997 0.20% MLCellLinOp::applyBC() 1750 0.3561 0.4198 0.4509 0.18% MultiFab::contains_nan() 18 0.375 0.4199 0.44 0.18% FabArray::Xpay() 291 0.3428 0.3844 0.4001 0.16% FillBoundary_nowait() 1561 0.3061 0.3557 0.3893 0.16% Amr::writePlotFile() 2 0.08498 0.2167 0.3641 0.15% Castro::initialize_do_advance() 9 0.256 0.2836 0.3133 0.13% FabArray::setDomainBndry() 92 0.2052 0.2199 0.2305 0.09% Castro::construct_old_source() 63 0.1906 0.205 0.221 0.09% amrex::average_down 233 0.1986 0.208 0.2124 0.09% Castro::fill_hybrid_hydro_source() 27 0.1699 0.1908 0.2069 0.08% Castro::initData() 1 0.1244 0.151 0.1758 0.07% MLCellLinOp::defineAuxData() 11 0.1562 0.1691 0.1745 0.07% MLPoisson::FFlux() 150 0.004671 0.04248 0.1558 0.06% MLCellLinOp::compGrad() 18 0.1201 0.1346 0.1412 0.06% FabArray::mult() 56 0.1167 0.1314 0.1383 0.06% StateData::define() 25 0.1232 0.1325 0.1366 0.06% Castro::swap_state_time_levels() 9 0.1015 0.1147 0.1304 0.05% Castro::evaluate_source_change() 193 0.1139 0.1197 0.1248 0.05% MLMG::addInterpCorrection() 155 0.1059 0.1173 0.1218 0.05% Gravity::get_new_grav_vector() 10 0.08517 0.09614 0.1051 0.04% Gravity::get_old_grav_vector() 9 0.0774 0.09171 0.1046 0.04% Castro::initialize_advance() 9 0.08094 0.09109 0.09891 0.04% AmrLevel::FillPatch() 62 0.08321 0.08847 0.09491 0.04% MLMG::prepareForSolve() 10 0.07199 0.08186 0.09386 0.04% Castro::apply_tagging_restrictions() 6 0.05505 0.06461 0.09138 0.04% MLCellLinOp::updateSolBC() 25 0.06765 0.06832 0.06941 0.03% MLMG::solve() 10 0.001051 0.04303 0.06753 0.03% DistributionMapping::LeastUsedCPUs() 5 0.004813 0.02641 0.06615 0.03% FabArray::norminf(mask) 21 0.05607 0.06312 0.06584 0.03% AmrMesh::MakeNewGrids() 6 0.04715 0.0541 0.06296 0.03% Gravity::multilevel_solve_for_new_phi() 10 0.04088 0.0446 0.04903 0.02% MLCellLinOp::updateCorBC() 17 0.04602 0.04645 0.04731 0.02% BndryData::define() 27 0.04349 0.04457 0.0457 0.02% MLMG::ResNormInf() 54 0.0009835 0.01849 0.04476 0.02% FabArray::norminf() 117 0.03652 0.03856 0.03973 0.02% MLMG::MLResNormInf() 15 9.926e-05 0.02478 0.03619 0.01% MLCellLinOp::reflux() 30 0.004287 0.01177 0.03368 0.01% MLLinOp::makeSubCommunicator() 11 0.0003215 0.01445 0.03367 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 622 0.01912 0.02253 0.02868 0.01% MLCellLinOp::defineBC() 11 0.02636 0.02655 0.0268 0.01% MultiFab::Multiply() 8 0.02233 0.02509 0.02613 0.01% Castro::FluxRegFineAdd() 9 0.0255 0.02573 0.02594 0.01% Castro::build_fine_mask() 44 0.01747 0.01974 0.02248 0.01% Castro::do_old_sources() 9 0.0006041 0.005786 0.0212 0.01% FabArrayBase::CPC::define() 416 0.01442 0.01632 0.01976 0.01% StateDataPhysBCFunct::() 95 0.00226 0.007886 0.01953 0.01% MLMG::interpCorrection_1 17 0.01842 0.01886 0.01948 0.01% Castro::buildMetrics() 5 0.0155 0.01779 0.01932 0.01% Castro::Castro() 5 0.01672 0.01804 0.01853 0.01% AMRErrorTag::operator() 6 0.00999 0.01261 0.01825 0.01% Castro::linear_to_hybrid_momentum() 1 0.01371 0.01573 0.01688 0.01% FillPatchTwoLevels 31 0.0135 0.01411 0.01557 0.01% Castro::do_new_sources() 9 0.0005264 0.0041 0.01469 0.01% Castro::construct_old_gravity() 9 0.0003928 0.005376 0.01382 0.01% Gravity::actual_multilevel_solve() 10 0.01203 0.01269 0.01312 0.01% Castro::finalize_do_advance() 9 0.008564 0.01094 0.01303 0.01% MLCellLinOp::prepareForSolve() 11 0.01165 0.01176 0.01198 0.00% Gravity::update_max_rhs() 10 0.00774 0.009631 0.01158 0.00% FabArrayBase::getCPC() 1440 0.01025 0.01062 0.01078 0.00% Castro::enforce_consistent_e() 1 0.009071 0.01024 0.01075 0.00% Castro::apply_problem_tags() 6 0.007394 0.008625 0.01038 0.00% Gravity::actual_solve_with_mlmg() 11 0.009441 0.009665 0.009949 0.00% AmrLevel::derive() 6 0.00432 0.005327 0.006359 0.00% average_down_faces 48 0.005725 0.005811 0.006039 0.00% MLCellLinOp::setLevelBC() 19 0.003679 0.004576 0.005499 0.00% MLCGSolver::bicgstab 7 0 0.001292 0.005169 0.00% FabArray::FillBoundary() 1561 0.00441 0.004587 0.005022 0.00% Castro::construct_new_gravity() 9 0.0004401 0.001622 0.003885 0.00% MLCellLinOp::smooth() 688 0.003784 0.00381 0.00383 0.00% FabArrayBase::getFB() 1561 0.003096 0.003241 0.003634 0.00% FabArrayBase::FB::FB() 19 0.002806 0.003157 0.003514 0.00% amrex::Dot() 115 0 0.0008461 0.003384 0.00% FabArray::ParallelCopy() 830 0.003079 0.003174 0.003288 0.00% Gravity::test_composite_phi() 1 0.002081 0.00237 0.002588 0.00% main() 1 0.001653 0.001777 0.001923 0.00% FabArray::LinComb() 57 0 0.0004704 0.001882 0.00% MLLinOp::defineGrids() 11 0.001412 0.001574 0.001842 0.00% MLCellLinOp::apply() 326 0.0008911 0.001064 0.001567 0.00% Amr::regrid() 5 0.001119 0.001251 0.001492 0.00% TagBoxArray::collate() 6 0.00059 0.001005 0.00147 0.00% MLCGSolver::ParallelAllReduce 153 0 0.000365 0.00146 0.00% Castro::clean_state() 87 0.001244 0.001282 0.001323 0.00% AmrLevel::FillCoarsePatch() 5 0.0009092 0.001088 0.001285 0.00% Amr::coarseTimeStep() 5 0.0005557 0.001033 0.001225 0.00% OwnerMask() 6 0.0008155 0.000943 0.001014 0.00% Gravity::average_fine_ec_onto_crse_ec() 8 0.0009087 0.0009335 0.0009462 0.00% MLCellLinOp::correctionResidual() 196 0.0007316 0.0007756 0.0008497 0.00% Castro::do_advance_ctu() 5 0.0006487 0.0006692 0.0007057 0.00% Castro::finalize_advance() 9 0.0005355 0.0006233 0.0006926 0.00% MLMG::mgVcycle() 65 0.0005858 0.0006053 0.0006223 0.00% FillPatchSingleLevel 92 0.000538 0.0005614 0.0005884 0.00% Amr::grid_places() 6 0.0004962 0.0005445 0.0005851 0.00% TagBoxArray::mapPRD 6 0.0004795 0.0005099 0.0005434 0.00% FillPatchIterator::FillFromLevel0() 31 0.0004075 0.0004523 0.0005014 0.00% Castro::subcycle_advance_ctu() 5 7.479e-05 0.0002243 0.0004556 0.00% MLMG:computeResOfCorrection() 155 0.0003906 0.0003967 0.000414 0.00% MLCellLinOp::solutionResidual() 72 0.0003843 0.0003904 0.0003993 0.00% MLMG::actualBottomSolve() 31 4.15e-05 0.0001278 0.0003837 0.00% MLMG::oneIter() 31 0.0003505 0.0003593 0.0003682 0.00% FillPatchIterator::FillFromTwoLevels() 31 0.0003091 0.0003249 0.0003457 0.00% Castro::computeNewDt() 8 0.0001223 0.0001773 0.0003348 0.00% MLMG::mgVcycle_down::0 31 0.0002815 0.0002929 0.0003062 0.00% Gravity::solve_phi_with_mlmg() 11 0.0002405 0.0002623 0.0002769 0.00% Amr::timeStep() 9 0.0002151 0.0002379 0.000271 0.00% Castro::enforce_speed_limit() 87 0.0002386 0.0002492 0.0002627 0.00% Amr::defBaseLevel() 1 0.0002063 0.000235 0.0002568 0.00% MLMG::mgVcycle_down::1 31 0.0002351 0.0002426 0.0002508 0.00% Castro::post_regrid() 10 0.0001948 0.0002142 0.0002504 0.00% MLMG::mgVcycle_bottom 65 0.0002146 0.0002181 0.0002224 0.00% AmrLevel::AmrLevel(dm) 5 0.0002014 0.0002124 0.0002223 0.00% makeSFC 33 0.0001797 0.0001948 0.0002143 0.00% MLMG::mgVcycle_up::0 31 0.0001826 0.0001918 0.0002044 0.00% MLMG::mgVcycle_down::2 31 0.0001937 0.0001971 0.0002022 0.00% MLMG::mgVcycle_down::3 31 0.0001958 0.0001974 0.0001991 0.00% FabArrayBase::TheFPinfo() 25 0.0001753 0.0001844 0.0001897 0.00% MLMG::mgVcycle_down::4 31 0.0001733 0.0001786 0.0001886 0.00% Castro::apply_source_to_state() 18 0.0001034 0.0001217 0.0001693 0.00% Amr::FinalizeInit() 1 3.326e-05 7.124e-05 0.0001683 0.00% Castro::check_for_nan() 18 0.0001265 0.0001413 0.0001545 0.00% Castro::react_state() 18 0.0001294 0.0001387 0.0001523 0.00% MLMG::mgVcycle_up::1 31 0.0001371 0.0001422 0.0001498 0.00% FPinfo::FPinfo() 3 0.0001259 0.0001331 0.0001467 0.00% MLMG::computeResWithCrseSolFineCor() 17 0.0001352 0.0001378 0.0001404 0.00% MLMG::mgVcycle_up::2 31 0.0001315 0.0001337 0.0001361 0.00% MLMG::mgVcycle_up::4 31 0.0001139 0.0001201 0.000132 0.00% DistributionMapping::Distribute() 35 0.0001148 0.0001214 0.0001319 0.00% MLMG::mgVcycle_up::3 31 0.000125 0.0001287 0.0001301 0.00% Gravity::swapTimeLevels() 9 9.496e-05 0.0001072 0.0001279 0.00% MLPoisson::define() 11 0.0001063 0.0001157 0.0001274 0.00% MLLinOp::makeAgglomeratedDMap 11 0.0001143 0.0001186 0.0001255 0.00% Castro::expand_state() 13 0.0001128 0.0001172 0.0001229 0.00% FabArray::RecvLayoutMask() 8 9.756e-05 0.0001055 0.0001218 0.00% Amr::writeSmallPlotFile() 1 9.556e-05 0.0001072 0.0001207 0.00% Castro::errorEst() 6 0.0001049 0.0001129 0.0001187 0.00% Castro::advance() 5 0.0001029 0.0001078 0.0001127 0.00% MLLinOp::define() 11 6.886e-05 8.061e-05 9.739e-05 0.00% MLMG::computeMLResidual() 15 7.786e-05 8.349e-05 9.288e-05 0.00% ClusterList::intersect() 1 0 2.288e-05 9.154e-05 0.00% MLMG::computeResidual() 31 8.137e-05 8.471e-05 8.87e-05 0.00% Castro::post_timestep() 9 8.455e-05 8.595e-05 8.711e-05 0.00% Castro::avgDown(state_indx) 20 7.06e-05 7.743e-05 8.377e-05 0.00% knapsack() 3 6.68e-05 7.02e-05 7.522e-05 0.00% MLMG::computeResWithCrseCorFineCor() 17 6.558e-05 6.745e-05 7.194e-05 0.00% MLMG::miniCycle() 34 6.808e-05 6.925e-05 7.108e-05 0.00% Castro::create_source_corrector() 9 5.596e-05 6.23e-05 6.684e-05 0.00% MLMG::getGradSolution() 10 5.943e-05 6.25e-05 6.587e-05 0.00% DistributionMapping::SFCProcessorMapDoIt() 2 4.911e-05 5.505e-05 6.169e-05 0.00% amrex::average_down_w_geom 20 5.029e-05 5.447e-05 5.614e-05 0.00% Castro::construct_new_hybrid_source() 9 4.407e-05 4.855e-05 5.275e-05 0.00% MLPoisson::prepareForSolve() 11 4.693e-05 4.92e-05 5.215e-05 0.00% Amr::initSubcycle() 1 4.893e-05 5.063e-05 5.14e-05 0.00% Castro::init(old) 3 4.403e-05 4.635e-05 5.139e-05 0.00% Gravity::set_mass_offset() 6 3.235e-05 4.267e-05 4.979e-05 0.00% Castro::retry_advance_ctu() 5 2.282e-05 3.024e-05 4.413e-05 0.00% Castro::construct_old_hybrid_source() 9 3.499e-05 3.71e-05 4.135e-05 0.00% Castro::avgDown() 4 3.372e-05 3.543e-05 4.036e-05 0.00% AmrMesh-cluster 1 0 9.976e-06 3.99e-05 0.00% DistributionMapping::KnapSackDoIt() 3 3.085e-05 3.4e-05 3.781e-05 0.00% BoxList::complementIn 12 2.46e-05 2.851e-05 3.715e-05 0.00% MLMG::MLRhsNormInf() 10 3.42e-05 3.48e-05 3.601e-05 0.00% Castro::post_init() 1 2.886e-05 3.164e-05 3.493e-05 0.00% BoxList::parallelComplementIn() 12 3.002e-05 3.16e-05 3.354e-05 0.00% Castro::construct_old_rotation_source() 9 2.449e-05 2.741e-05 3.184e-05 0.00% MLMG::compResidual() 1 2.699e-05 2.893e-05 3.12e-05 0.00% Castro::construct_new_rotation_source() 9 2.076e-05 2.316e-05 2.669e-05 0.00% Amr::bldFineLevels() 1 1.787e-05 1.993e-05 2.386e-05 0.00% Castro::computeInitialDt() 2 1.593e-05 1.93e-05 2.366e-05 0.00% Castro::derive() 6 1.822e-05 1.934e-05 2.038e-05 0.00% Castro::init() 1 1.636e-05 1.663e-05 1.679e-05 0.00% Amr::InitializeInit() 1 8.097e-06 1.074e-05 1.417e-05 0.00% Amr::init() 1 7.25e-06 9.494e-06 1.376e-05 0.00% ClusterList::chop() 1 0 2.453e-06 9.811e-06 0.00% Amr::initialInit() 1 4.194e-06 5.974e-06 8.952e-06 0.00% knapsack()swap 1 3.787e-06 3.877e-06 4.003e-06 0.00% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 244.5 244.5 244.5 100.00% Amr::coarseTimeStep() 5 231.9 231.9 231.9 94.85% Amr::timeStep() 9 230.6 230.6 230.6 94.30% Castro::advance() 5 202.8 202.8 202.8 82.96% Castro::subcycle_advance_ctu() 5 194.1 194.5 195.4 79.90% Castro::do_advance_ctu() 5 194.1 194.5 195.4 79.90% Castro::construct_ctu_hydro_source() 9 117.4 117.6 117.7 48.15% Castro::clean_state() 87 71.74 81.31 85.61 35.02% Castro::computeTemp() 88 45.51 51.53 54.47 22.28% Castro::reset_internal_energy(MultiFab) 88 30.81 34.91 37.11 15.18% Castro::evaluate_source_change() 193 17.73 19.8 20.78 8.50% Castro::enforce_min_density() 87 16.94 19.27 20.68 8.46% Castro::post_regrid() 10 15.94 17.53 18.17 7.43% Amr::regrid() 5 15.84 17.25 17.83 7.29% Castro::initialize_do_advance() 9 14.93 15.3 15.48 6.33% check_for_negative_density() 9 0.3131 4.529 15.17 6.21% amrex::Copy() 736 12.86 14.31 15.06 6.16% AmrLevel::FillPatch() 62 7.932 10.11 14.64 5.99% FillPatchIterator::Initialize 62 6.194 8.439 13.13 5.37% FabArray::ParallelCopy() 830 3.699 6.801 13.07 5.34% Castro::normalize_species() 87 10.08 11.41 12.72 5.20% Gravity::solve_phi_with_mlmg() 11 11.23 11.61 12.39 5.07% MultiFab::AddProduct() 4205 10.46 11.67 12.23 5.00% Gravity::multilevel_solve_for_new_phi() 10 10.96 11.34 12.11 4.95% Gravity::actual_multilevel_solve() 10 10.92 11.29 12.07 4.94% Castro::do_old_sources() 9 11.77 11.83 11.93 4.88% Castro::do_new_sources() 9 9.646 10.95 11.54 4.72% FillPatchSingleLevel 92 3.866 6.154 10.97 4.49% FabArray::ParallelCopy_finish() 830 1.481 4.412 10.78 4.41% Amr::init() 1 10.77 10.77 10.77 4.40% Castro::post_timestep() 9 10.4 10.41 10.41 4.26% FillPatchIterator::FillFromLevel0() 31 3.477 5.686 10.35 4.23% Amr::initialInit() 1 10.1 10.11 10.16 4.16% Castro::reset_internal_energy(Fab) 1476 8.284 9.247 9.601 3.93% Castro::estTimeStep() 26 4.022 5.691 9.534 3.90% Castro::reflux() 4 8.626 9.176 9.427 3.86% FabArray::setVal() 5078 7.347 8.215 8.742 3.58% Castro::initialize_advance() 9 7.465 8.354 8.737 3.57% MultiFab::Subtract() 175 6.646 7.433 7.77 3.18% Castro::construct_new_gravity() 9 6.44 6.794 7.551 3.09% Gravity::actual_solve_with_mlmg() 11 6.858 6.871 6.885 2.82% Castro::expand_state() 13 5.191 5.769 6.641 2.72% MLMG::solve() 10 6.35 6.362 6.38 2.61% Amr::FinalizeInit() 1 5.852 6.021 6.354 2.60% Gravity::fill_multipole_BCs() 11 4.351 4.721 5.49 2.25% MLMG::oneIter() 31 5.131 5.148 5.164 2.11% FluxRegister::Reflux() 24 3.883 4.398 4.694 1.92% Castro::finalize_do_advance() 9 1.346 2.343 4.575 1.87% MLMG::mgVcycle() 65 4.343 4.399 4.426 1.81% Castro::construct_old_gravity() 9 0.7224 1.751 4.324 1.77% Gravity::get_old_grav_vector() 9 0.7086 1.746 4.323 1.77% Amr::InitializeInit() 1 3.745 4.093 4.247 1.74% Amr::defBaseLevel() 1 3.745 4.093 4.247 1.74% Castro::problem_post_timestep() 10 2.768 3.213 4.094 1.67% FabArray::sum() 4205 3.446 3.857 4.06 1.66% Castro::computeNewDt() 8 2.183 2.699 3.96 1.62% MLCellLinOp::smooth() 688 3.296 3.355 3.407 1.39% Castro::FluxRegCrseInit 9 1.88 2.331 3.11 1.27% Castro::post_init() 1 2.453 2.598 2.916 1.19% Castro::initData() 1 2.308 2.636 2.777 1.14% FabArray::ParallelCopy_nowait() 830 2.215 2.385 2.77 1.13% Castro::wd_update() 6 1.175 1.667 2.564 1.05% Castro::construct_new_source() 63 2.09 2.344 2.491 1.02% MLPoisson::Fsmooth() 1376 2.067 2.328 2.429 0.99% MLMG::mgVcycle_down::0 31 2.002 2.044 2.083 0.85% MLCellLinOp::applyBC() 1750 1.233 1.412 1.847 0.76% FillPatchIterator::FillFromTwoLevels() 31 1.536 1.641 1.775 0.73% FillPatchTwoLevels 31 1.535 1.641 1.775 0.73% FabArray::Saxpy() 125 1.471 1.649 1.728 0.71% Amr::bldFineLevels() 1 1.428 1.584 1.647 0.67% FabArray::FillBoundary() 1561 0.8005 0.9926 1.491 0.61% Castro::Castro() 5 1.389 1.412 1.448 0.59% MLMG::mgVcycle_up::0 31 1.382 1.406 1.429 0.58% FillPatchInterp(MF) 30 1.217 1.328 1.428 0.58% Castro::construct_old_source() 63 1.095 1.229 1.409 0.58% MLCellLinOp::apply() 326 1.285 1.303 1.334 0.55% amrex::Add() 143 1.063 1.194 1.3 0.53% Amr::writePlotFile() 2 1.193 1.239 1.255 0.51% Amr::InitAmr() 1 1.23 1.23 1.23 0.50% Castro::construct_old_gravity_source() 18 0.9225 1.045 1.209 0.49% Castro::initMFs() 5 1.106 1.153 1.177 0.48% FillBoundary_finish() 1561 0.4048 0.6259 1.092 0.45% MLCellLinOp::solutionResidual() 72 0.9985 1.022 1.041 0.43% MLPoisson::Fapply() 326 0.885 0.9854 1.024 0.42% Castro::computeInitialDt() 2 0.5059 0.6595 1.008 0.41% Castro::construct_new_gravity_source() 9 0.7489 0.8431 0.913 0.37% Gravity::get_new_grav_vector() 10 0.7939 0.8088 0.8335 0.34% VisMF::Write(FabArray) 3 0.6387 0.7443 0.8239 0.34% MLCellLinOp::correctionResidual() 196 0.8129 0.8151 0.8179 0.33% Castro::volInBoundary() 42 0.5747 0.6512 0.7695 0.31% MLMG:computeResOfCorrection() 155 0.6542 0.6819 0.694 0.28% Castro::apply_source_to_state() 18 0.5897 0.6605 0.6938 0.28% StateDataPhysBCFunct::() 95 0.6149 0.6303 0.6442 0.26% StateData::FillBoundary(geom) 50 0.6105 0.6224 0.639 0.26% MLMG::computeResWithCrseSolFineCor() 17 0.4646 0.4973 0.5792 0.24% Gravity::update_max_rhs() 10 0.08111 0.2163 0.5528 0.23% MLMG::computeMLResidual() 15 0.4177 0.4546 0.505 0.21% MultiFab::max() 18 0.01806 0.1569 0.4997 0.20% Castro::init() 1 0.4757 0.4776 0.4797 0.20% AmrLevel::FillCoarsePatch() 5 0.4757 0.4776 0.4797 0.20% Castro::check_for_nan() 18 0.3752 0.42 0.4402 0.18% MultiFab::contains_nan() 18 0.375 0.4199 0.44 0.18% Gravity::test_composite_phi() 1 0.4023 0.4032 0.4049 0.17% FabArray::Xpay() 291 0.3428 0.3844 0.4001 0.16% FillBoundary_nowait() 1561 0.312 0.3621 0.3962 0.16% MLMG::mgVcycle_down::1 31 0.2985 0.3339 0.3866 0.16% Castro::init(old) 3 0.3643 0.3736 0.3778 0.15% MLCellLinOp::reflux() 30 0.187 0.2366 0.3672 0.15% MLMG::computeResidual() 31 0.3445 0.3508 0.3544 0.14% Castro::react_state() 18 0.2338 0.2698 0.3154 0.13% MLPoisson::define() 11 0.2922 0.304 0.3139 0.13% amrex::average_down 233 0.2692 0.2777 0.285 0.12% Amr::grid_places() 6 0.2686 0.2687 0.2688 0.11% AmrMesh::MakeNewGrids() 6 0.268 0.2682 0.2683 0.11% MLMG::prepareForSolve() 10 0.2018 0.2237 0.2403 0.10% FabArray::setDomainBndry() 92 0.2052 0.2199 0.2305 0.09% MLMG::mgVcycle_up::1 31 0.2067 0.2113 0.2137 0.09% MLCellLinOp::defineAuxData() 11 0.1827 0.1985 0.2089 0.09% Castro::fill_hybrid_hydro_source() 27 0.1699 0.1908 0.2069 0.08% MLMG::mgVcycle_bottom 65 0.1533 0.162 0.1874 0.08% MLMG::addInterpCorrection() 155 0.127 0.1677 0.1864 0.08% Castro::errorEst() 6 0.1213 0.14 0.182 0.07% MLMG::getGradSolution() 10 0.1688 0.1717 0.1734 0.07% MLCellLinOp::compGrad() 18 0.1687 0.1717 0.1734 0.07% MLPoisson::FFlux() 150 0.004671 0.04248 0.1558 0.06% MLMG::miniCycle() 34 0.1533 0.1537 0.154 0.06% Castro::avgDown() 4 0.1258 0.1361 0.1464 0.06% Castro::avgDown(state_indx) 20 0.1258 0.1361 0.1464 0.06% amrex::average_down_w_geom 20 0.1257 0.136 0.1463 0.06% FabArray::mult() 56 0.1167 0.1314 0.1383 0.06% AmrLevel::AmrLevel(dm) 5 0.1235 0.1327 0.1369 0.06% StateData::define() 25 0.1232 0.1325 0.1366 0.06% MLCellLinOp::updateSolBC() 25 0.07771 0.09885 0.1354 0.06% Castro::construct_new_hybrid_source() 9 0.1126 0.1265 0.1349 0.06% Castro::swap_state_time_levels() 9 0.1015 0.1147 0.1304 0.05% MLMG::ResNormInf() 54 0.07315 0.08823 0.1071 0.04% Castro::apply_tagging_restrictions() 6 0.05505 0.06461 0.09138 0.04% MLCellLinOp::defineBC() 11 0.08731 0.08884 0.09075 0.04% MLMG::computeResWithCrseCorFineCor() 17 0.08768 0.08939 0.09038 0.04% TagBoxArray::mapPRD 6 0.009295 0.06121 0.08395 0.03% MLMG::mgVcycle_up::3 31 0.02521 0.06908 0.08382 0.03% MLMG::mgVcycle_down::2 31 0.0696 0.07429 0.08205 0.03% MLMG::MLResNormInf() 15 0.04541 0.06848 0.07683 0.03% Castro::construct_old_hybrid_source() 9 0.05741 0.06446 0.07204 0.03% FabArrayBase::TheFPinfo() 25 0.005189 0.02678 0.06654 0.03% FPinfo::FPinfo() 3 0.005004 0.0266 0.06635 0.03% DistributionMapping::KnapSackDoIt() 3 0.004876 0.02646 0.0662 0.03% DistributionMapping::LeastUsedCPUs() 5 0.004813 0.02641 0.06615 0.03% FabArray::norminf(mask) 21 0.05607 0.06312 0.06584 0.03% MLMG::interpCorrection_1 17 0.02638 0.04335 0.06547 0.03% BndryData::define() 27 0.05943 0.06102 0.06307 0.03% Castro::derive() 6 0.04731 0.05399 0.06295 0.03% AmrLevel::derive() 6 0.04729 0.05397 0.06294 0.03% MLCellLinOp::updateCorBC() 17 0.05084 0.05179 0.05237 0.02% MLMG::mgVcycle_up::2 31 0.05037 0.05105 0.05166 0.02% Castro::buildMetrics() 5 0.03958 0.04241 0.04539 0.02% Gravity::swapTimeLevels() 9 0.0362 0.03835 0.04084 0.02% FabArray::norminf() 117 0.03652 0.03856 0.03973 0.02% MLLinOp::define() 11 0.002337 0.01653 0.03606 0.01% MLLinOp::defineGrids() 11 0.002257 0.01645 0.03596 0.01% MLMG::mgVcycle_down::3 31 0.03068 0.03175 0.03469 0.01% MLLinOp::makeSubCommunicator() 11 0.0003215 0.01445 0.03367 0.01% MLMG::actualBottomSolve() 31 4.15e-05 0.008393 0.03344 0.01% MLCGSolver::bicgstab 7 0 0.008182 0.03273 0.01% FabArrayBase::getCPC() 1440 0.02467 0.02694 0.03052 0.01% MLMG::MLRhsNormInf() 10 0.02549 0.02855 0.02967 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 622 0.01912 0.02253 0.02868 0.01% MultiFab::Multiply() 8 0.02233 0.02509 0.02613 0.01% Castro::FluxRegFineAdd() 9 0.0255 0.02573 0.02594 0.01% Castro::build_fine_mask() 44 0.01747 0.01974 0.02248 0.01% Gravity::average_fine_ec_onto_crse_ec() 8 0.01498 0.01705 0.02175 0.01% MLMG::compResidual() 1 0.01749 0.01909 0.02001 0.01% FabArrayBase::CPC::define() 416 0.01442 0.01632 0.01976 0.01% MLMG::mgVcycle_down::4 31 0.005197 0.008526 0.0184 0.01% AMRErrorTag::operator() 6 0.00999 0.01261 0.01825 0.01% Castro::linear_to_hybrid_momentum() 1 0.01371 0.01573 0.01688 0.01% MLMG::mgVcycle_up::4 31 0.004268 0.006726 0.01397 0.01% MLPoisson::prepareForSolve() 11 0.01171 0.0118 0.01202 0.00% MLCellLinOp::prepareForSolve() 11 0.01165 0.01176 0.01198 0.00% Castro::enforce_consistent_e() 1 0.009071 0.01024 0.01075 0.00% average_down_faces 48 0.009079 0.00986 0.01039 0.00% Castro::apply_problem_tags() 6 0.007394 0.008625 0.01038 0.00% FabArrayBase::getFB() 1561 0.005902 0.006398 0.006902 0.00% MLCellLinOp::setLevelBC() 19 0.003679 0.004576 0.005499 0.00% FabArrayBase::FB::FB() 19 0.002806 0.003157 0.003514 0.00% amrex::Dot() 115 0 0.0008461 0.003384 0.00% MLCGSolver::ParallelAllReduce 153 0 0.0008156 0.003262 0.00% FabArray::LinComb() 57 0 0.0004704 0.001882 0.00% TagBoxArray::collate() 6 0.00059 0.001005 0.00147 0.00% OwnerMask() 6 0.0008155 0.000943 0.001014 0.00% Castro::finalize_advance() 9 0.0005355 0.0006233 0.0006926 0.00% MLLinOp::makeAgglomeratedDMap 11 0.0003954 0.0004226 0.000448 0.00% makeSFC 33 0.000281 0.0003039 0.0003225 0.00% Castro::enforce_speed_limit() 87 0.0002386 0.0002492 0.0002627 0.00% AmrMesh-cluster 1 0 3.531e-05 0.0001413 0.00% DistributionMapping::Distribute() 35 0.0001148 0.0001214 0.0001319 0.00% DistributionMapping::SFCProcessorMapDoIt() 2 0.0001121 0.0001216 0.0001261 0.00% FabArray::RecvLayoutMask() 8 9.756e-05 0.0001055 0.0001218 0.00% Amr::writeSmallPlotFile() 1 9.556e-05 0.0001072 0.0001207 0.00% ClusterList::intersect() 1 0 2.288e-05 9.154e-05 0.00% knapsack() 3 7.08e-05 7.407e-05 7.912e-05 0.00% BoxList::parallelComplementIn() 12 5.462e-05 6.011e-05 6.967e-05 0.00% Castro::create_source_corrector() 9 5.596e-05 6.23e-05 6.684e-05 0.00% Amr::initSubcycle() 1 4.893e-05 5.063e-05 5.14e-05 0.00% Gravity::set_mass_offset() 6 3.235e-05 4.267e-05 4.979e-05 0.00% Castro::retry_advance_ctu() 5 2.282e-05 3.024e-05 4.413e-05 0.00% BoxList::complementIn 12 2.46e-05 2.851e-05 3.715e-05 0.00% Castro::construct_old_rotation_source() 9 2.449e-05 2.741e-05 3.184e-05 0.00% Castro::construct_new_rotation_source() 9 2.076e-05 2.316e-05 2.669e-05 0.00% ClusterList::chop() 1 0 2.453e-06 9.811e-06 0.00% knapsack()swap 1 3.787e-06 3.877e-06 4.003e-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.07-13-gb41f85ab36ac) finalized