Initializing AMReX (25.06-22-g1d307c2c661a-dirty)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (25.06-22-g1d307c2c661a-dirty) initialized Starting run at 06:47:45 UTC on 2025-07-25. Successfully read inputs file ... Castro git describe: 25.07-12-g2c8f2fc1e AMReX git describe: 25.06-22-g1d307c2c6-dirty Microphysics git describe: joss_submit reading extern runtime parameters ... 13 Species: He4 C12 O16 Ne20 Mg24 Si28 S32 Ar36 Ca40 Ti44 Cr48 Fe52 Ni56 Successfully read inputs file ... Creating primary with CO core with mass fractions C = 0.5and O = 0.5 and a He shell of solar mass =0. Creating secondary with CO core with mass fractions C = 0.5and O = 0.5 and a He shell of solar mass =0. Generated initial model for primary WD of mass 0.64 solar masses, central density 4.500e+06 g cm**-3, and radius 8.734e+08 cm. Generated initial model for secondary WD of mass 0.64 solar masses, central density 4.500e+06 g cm**-3, and radius 8.734e+08 cm. Castro::numpts_1d at level 0 is 294 Setting the gravity type to PoissonGrav Initializing the data at level 0 Done initializing the level 0 data Now regridding at level lbase = 0 grid_places() time: 0.17950971 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001190335 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010785599 Gravity::multilevel_solve_for_new_phi() time = 0.110097093 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010710352 TIME= 0 MASS = 2.549996568e+33 TIME= 0 XMOM = 0 TIME= 0 YMOM = 3.48170636e+26 TIME= 0 ZMOM = 0 TIME= 0 ANG MOM X = 0 TIME= 0 ANG MOM Y = 0 TIME= 0 ANG MOM Z = 2.076918743e+34 TIME= 0 RHO*e = 1.325915365e+50 TIME= 0 RHO*K = 3.09960185e+49 TIME= 0 RHO*E = 1.63587555e+50 TIME= 0 RHO*PHI = -5.434813569e+50 TIME= 0 TOTAL ENERGY= -1.081531234e+50 TIME= 0 CENTER OF MASS X-LOC = 316439254.4 TIME= 0 CENTER OF MASS X-VEL = 0 TIME= 0 CENTER OF MASS Y-LOC = -3.641005016e-07 TIME= 0 CENTER OF MASS Y-VEL = 1.365376881e-07 TIME= 0 CENTER OF MASS Z-LOC = 0 TIME= 0 CENTER OF MASS Z-VEL = 0 TIME= 0 MAXIMUM TEMPERATURE = 10000000 TIME= 0 MAXIMUM DENSITY = 4476474 TIME= 0 MAXIMUM T_S / T_E = 0 INITIAL GRIDS Level 0 32 grids 32768 cells 100 % of domain smallest grid: 32 x 32 biggest grid: 32 x 32 PLOTFILE: file = wdmerger_retry-SDC_plt00000 Write plotfile time = 0.049609897 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004141701 seconds [Level 0 step 1] ADVANCE at time 0 with dt = 0.25  Beginning subcycle 1 starting at time 0 with dt = 0.25 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.01295059363 at index 21, 172, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035502 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477044 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.769409 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071148 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.7670713549 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.7670713549 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071148 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.769409 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477044 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035502 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -0.01295059363 at index 21, 83, 0 Castro::construct_old_gravity() time = 0.002900928 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: -5.111789233e+47 ener added: 3.959219145e+48 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.174137392e+41 ymom added: 6.503960463e+28 zmom added: 0 eint added: 0 ener added: 6.931398274e+47 Castro::do_old_sources() time = 0.065963144 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Castro::construct_old_gravity() time = 0.004324409 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: -5.024842118e+47 ener added: 6.151153052e+47 Contributions to the state from old-time sources: mass added: 0 xmom added: -5.870686959e+40 ymom added: 3.251980232e+28 zmom added: 0 eint added: 0 ener added: 3.465699137e+47 Castro::do_old_sources() time = 0.065117038 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.125 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.0625   Beginning subcycle 1 starting at time 0 with dt = 0.0625 Estimated number of subcycles remaining: 4 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.1060744788 at index 20, 177, 0 Invalid density = -0.1060744788 at index 20, 78, 0 Castro::construct_old_gravity() time = 0.004800734 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -2.93534348e+40 ymom added: 1.625990116e+28 zmom added: 0 eint added: 0 ener added: 1.732849568e+47 Castro::do_old_sources() time = 0.061779736 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.0625 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.03125   Beginning subcycle 1 starting at time 0 with dt = 0.03125 Estimated number of subcycles remaining: 8 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Castro::construct_old_gravity() time = 0.005128908 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.129950579e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.060825108 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at level 0 Castro::construct_ctu_hydro_source() time = 0.340683181 on level 0 Ending SDC iteration 1 of 3. Beginning SDC iteration 2 of 3. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.015759261 Gravity::solve_for_phi() time = 0.094057606 on level 0 Castro::construct_new_gravity() time = 0.096229471 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.673145426e+44 ener added: 3.81935025e+42 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.752993387e+35 ymom added: -3.588950123e+27 zmom added: 0 eint added: 0 ener added: 6.222804276e+43 Castro::do_new_sources() time = 0.061437126 on level 0 Castro::construct_old_gravity() time = 0.00376066 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.129950579e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.060812491 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at level 0 Castro::construct_ctu_hydro_source() time = 0.338065811 on level 0 Ending SDC iteration 2 of 3. Beginning SDC iteration 3 of 3. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.01584121 Gravity::solve_for_phi() time = 0.094105359 on level 0 Castro::construct_new_gravity() time = 0.096265923 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.837104563e+44 ener added: 3.792000243e+42 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.742618664e+35 ymom added: -3.590636008e+27 zmom added: 0 eint added: 0 ener added: 4.583540443e+43 Castro::do_new_sources() time = 0.061286098 on level 0 Castro::construct_old_gravity() time = 0.003959279 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.129950579e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.060886328 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at level 0 Castro::construct_ctu_hydro_source() time = 0.343824025 on level 0 Ending SDC iteration 3 of 3.  Subcycle completed  Beginning subcycle 2 starting at time 0.03125 with dt = 0.03125 Estimated number of subcycles remaining: 7 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.015869324 Gravity::solve_for_phi() time = 0.100363127 on level 0 Castro::construct_new_gravity() time = 0.102461822 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.837201172e+44 ener added: 3.792181427e+42 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.742608901e+35 ymom added: -3.595677134e+27 zmom added: 0 eint added: 0 ener added: 4.58255499e+43 Castro::do_new_sources() time = 0.061398214 on level 0 Castro::construct_old_gravity() time = 0.002974779 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.985021797e+46 ener added: 5.427586462e+40 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467706592e+40 ymom added: 9.380508782e+26 zmom added: 0 eint added: 0 ener added: 8.721813599e+46 Castro::do_old_sources() time = 0.060843985 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.985569545e+46 ener added: 2.746016564e+42 Castro::construct_ctu_hydro_source() time = 0.349199061 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016205626 Gravity::solve_for_phi() time = 0.095661705 on level 0 Castro::construct_new_gravity() time = 0.097835038 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.269601218e+44 ener added: 4.892522312e+43 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.763301424e+35 ymom added: 4.919525284e+24 zmom added: 0 eint added: 0 ener added: 6.012622799e+43 Castro::do_new_sources() time = 0.061136673 on level 0 Castro::construct_old_gravity() time = 0.00353732 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.985021797e+46 ener added: 5.427586462e+40 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467706592e+40 ymom added: 9.380508782e+26 zmom added: 0 eint added: 0 ener added: 8.721813599e+46 Castro::do_old_sources() time = 0.060859741 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.986269155e+46 ener added: 2.7526843e+42 Castro::construct_ctu_hydro_source() time = 0.356709685 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 3 starting at time 0.0625 with dt = 0.03125 Estimated number of subcycles remaining: 6 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.01625096 Gravity::solve_for_phi() time = 0.094845127 on level 0 Castro::construct_new_gravity() time = 0.096931214 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.343102187e+44 ener added: 4.888810292e+43 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.757196357e+35 ymom added: 6.54874172e+24 zmom added: 0 eint added: 0 ener added: 5.281428219e+43 Castro::do_new_sources() time = 0.061196944 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.910635222e+29 ener added: 1.910635222e+29 Castro::construct_old_gravity() time = 0.003656619 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.984353231e+46 ener added: 6.990279267e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467741736e+40 ymom added: 9.545224924e+26 zmom added: 0 eint added: 0 ener added: 8.780735363e+46 Castro::do_old_sources() time = 0.060846257 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.987321945e+46 ener added: 2.18047792e+43 Castro::construct_ctu_hydro_source() time = 0.349225645 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016219147 Gravity::solve_for_phi() time = 0.095315577 on level 0 Castro::construct_new_gravity() time = 0.097501657 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.031202547e+43 ener added: 1.184458535e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.781734557e+35 ymom added: 5.277244545e+23 zmom added: 0 eint added: 0 ener added: 5.994705329e+43 Castro::do_new_sources() time = 0.062448687 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.910635222e+29 ener added: 1.910635222e+29 Castro::construct_old_gravity() time = 0.004009991 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.984353231e+46 ener added: 6.990279267e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467741736e+40 ymom added: 9.545224924e+26 zmom added: 0 eint added: 0 ener added: 8.780735363e+46 Castro::do_old_sources() time = 0.060836045 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.988028457e+46 ener added: 2.185423907e+43 Castro::construct_ctu_hydro_source() time = 0.354245142 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 4 starting at time 0.09375 with dt = 0.03125 Estimated number of subcycles remaining: 5 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016130019 Gravity::solve_for_phi() time = 0.093847924 on level 0 Castro::construct_new_gravity() time = 0.095973783 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.773882032e+43 ener added: 1.183651189e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.77560788e+35 ymom added: 2.905435979e+24 zmom added: 0 eint added: 0 ener added: 5.260150257e+43 Castro::do_new_sources() time = 0.061330558 on level 0 Castro::construct_old_gravity() time = 0.003003747 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.983917963e+46 ener added: 1.154261437e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467777248e+40 ymom added: 9.587537328e+26 zmom added: 0 eint added: 0 ener added: 8.84101861e+46 Castro::do_old_sources() time = 0.061054435 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.984973203e+46 ener added: 2.291971287e+43 Castro::construct_ctu_hydro_source() time = 0.348614473 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.01627449 Gravity::solve_for_phi() time = 0.093129891 on level 0 Castro::construct_new_gravity() time = 0.095235607 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: -4.612861023e+43 ener added: 1.488715317e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.801491314e+35 ymom added: 3.612610359e+23 zmom added: 0 eint added: 0 ener added: 4.637706889e+43 Castro::do_new_sources() time = 0.061268154 on level 0 Castro::construct_old_gravity() time = 0.003532434 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.983917963e+46 ener added: 1.154261437e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467777248e+40 ymom added: 9.587537328e+26 zmom added: 0 eint added: 0 ener added: 8.84101861e+46 Castro::do_old_sources() time = 0.060892519 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.985683204e+46 ener added: 2.297906432e+43 Castro::construct_ctu_hydro_source() time = 0.354053213 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 5 starting at time 0.125 with dt = 0.03125 Estimated number of subcycles remaining: 4 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016264772 Gravity::solve_for_phi() time = 0.093148586 on level 0 Castro::construct_new_gravity() time = 0.095282403 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: -5.357824616e+43 ener added: 1.487804348e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.795347886e+35 ymom added: -1.042462401e+25 zmom added: 0 eint added: 0 ener added: 3.901837814e+43 Castro::do_new_sources() time = 0.061273085 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.371823034e+29 ener added: 2.371823034e+29 Castro::construct_old_gravity() time = 0.003025603 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.983564871e+46 ener added: 1.521590942e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467813155e+40 ymom added: 9.391842461e+26 zmom added: 0 eint added: 0 ener added: 8.902573541e+46 Castro::do_old_sources() time = 0.060704969 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.992394606e+46 ener added: 2.321200641e+43 Castro::construct_ctu_hydro_source() time = 0.350315404 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016346217 Gravity::solve_for_phi() time = 0.093399825 on level 0 Castro::construct_new_gravity() time = 0.095528301 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.122723364e+43 ener added: 1.619563371e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.826309838e+35 ymom added: 1.288025458e+24 zmom added: 0 eint added: 0 ener added: 2.081560288e+43 Castro::do_new_sources() time = 0.061387413 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.371823034e+29 ener added: 2.371823034e+29 Castro::construct_old_gravity() time = 0.003623555 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.983564871e+46 ener added: 1.521590942e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467813155e+40 ymom added: 9.391842461e+26 zmom added: 0 eint added: 0 ener added: 8.902573541e+46 Castro::do_old_sources() time = 0.061771838 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.993114585e+46 ener added: 2.326635407e+43 Castro::construct_ctu_hydro_source() time = 0.355354941 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 6 starting at time 0.15625 with dt = 0.03125 Estimated number of subcycles remaining: 3 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016434992 Gravity::solve_for_phi() time = 0.09360208 on level 0 Castro::construct_new_gravity() time = 0.095729806 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.874849314e+43 ener added: 1.618719464e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.820136286e+35 ymom added: 2.815711015e+24 zmom added: 0 eint added: 0 ener added: 1.337821578e+43 Castro::do_new_sources() time = 0.061259528 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: 5.138949907e+29 ener added: 5.138949907e+29 Castro::construct_old_gravity() time = 0.003674235 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.983077706e+46 ener added: 2.015109106e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467849558e+40 ymom added: 9.409976348e+26 zmom added: 0 eint added: 0 ener added: 8.964911027e+46 Castro::do_old_sources() time = 0.060742401 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.000666253e+46 ener added: 2.405830666e+43 Castro::construct_ctu_hydro_source() time = 0.349052224 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016585867 Gravity::solve_for_phi() time = 0.09378193 on level 0 Castro::construct_new_gravity() time = 0.095908511 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: -8.911535206e+43 ener added: 1.713555602e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.844610083e+35 ymom added: -2.677581796e+24 zmom added: 0 eint added: 0 ener added: -1.419860559e+43 Castro::do_new_sources() time = 0.061551148 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: 5.138949907e+29 ener added: 5.138949907e+29 Castro::construct_old_gravity() time = 0.003871176 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.983077706e+46 ener added: 2.015109106e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467849558e+40 ymom added: 9.409976348e+26 zmom added: 0 eint added: 0 ener added: 8.964911027e+46 Castro::do_old_sources() time = 0.060798085 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.001390031e+46 ener added: 2.412276394e+43 Castro::construct_ctu_hydro_source() time = 0.354563025 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 7 starting at time 0.1875 with dt = 0.03125 Estimated number of subcycles remaining: 2 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016751234 Gravity::solve_for_phi() time = 0.093725716 on level 0 Castro::construct_new_gravity() time = 0.095839312 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.667367417e+43 ener added: 1.712648092e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.83840701e+35 ymom added: 5.636144397e+24 zmom added: 0 eint added: 0 ener added: -2.166675325e+43 Castro::do_new_sources() time = 0.061275938 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.038898262e+29 ener added: 3.038898262e+29 Castro::construct_old_gravity() time = 0.003764269 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.982914015e+46 ener added: 2.174744554e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467886326e+40 ymom added: 9.585270592e+26 zmom added: 0 eint added: 0 ener added: 9.027797292e+46 Castro::do_old_sources() time = 0.060711699 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.9938659e+46 ener added: 2.675763102e+43 Castro::construct_ctu_hydro_source() time = 0.349160705 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016578975 Gravity::solve_for_phi() time = 0.095942068 on level 0 Castro::construct_new_gravity() time = 0.098435849 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.153618601e+44 ener added: 1.743191878e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.87433057e+35 ymom added: -3.282044706e+23 zmom added: 0 eint added: 0 ener added: -4.056999371e+43 Castro::do_new_sources() time = 0.066448789 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.038898262e+29 ener added: 3.038898262e+29 Castro::construct_old_gravity() time = 0.007036931 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.982914015e+46 ener added: 2.174744554e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467886326e+40 ymom added: 9.585270592e+26 zmom added: 0 eint added: 0 ener added: 9.027797292e+46 Castro::do_old_sources() time = 0.060812942 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.99459825e+46 ener added: 2.682276214e+43 Castro::construct_ctu_hydro_source() time = 0.354655725 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 8 starting at time 0.21875 with dt = 0.03125 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016697343 Gravity::solve_for_phi() time = 0.093689276 on level 0 Castro::construct_new_gravity() time = 0.095827498 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.230005451e+44 ener added: 1.742291039e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.86809364e+35 ymom added: -8.38456169e+24 zmom added: 0 eint added: 0 ener added: -4.811917635e+43 Castro::do_new_sources() time = 0.061337675 on level 0 Castro::construct_old_gravity() time = 0.002864877 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.982547556e+46 ener added: 2.529310204e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467923688e+40 ymom added: 9.382019939e+26 zmom added: 0 eint added: 0 ener added: 9.091359292e+46 Castro::do_old_sources() time = 0.06067386 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.993424638e+46 ener added: 2.936823048e+43 Castro::construct_ctu_hydro_source() time = 0.349201859 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.017021701 Gravity::solve_for_phi() time = 0.094074729 on level 0 Castro::construct_new_gravity() time = 0.096240871 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: -8.584986492e+43 ener added: 1.766412994e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.939839844e+35 ymom added: 3.243085182e+24 zmom added: 0 eint added: 0 ener added: -1.081736549e+43 Castro::do_new_sources() time = 0.061443871 on level 0 Castro::construct_old_gravity() time = 0.003577559 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.982547556e+46 ener added: 2.529310204e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467923688e+40 ymom added: 9.382019939e+26 zmom added: 0 eint added: 0 ener added: 9.091359292e+46 Castro::do_old_sources() time = 0.060821441 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.994157717e+46 ener added: 2.94327526e+43 Castro::construct_ctu_hydro_source() time = 0.358438646 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002646285019 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.015908591 Gravity::solve_for_phi() time = 0.094671164 on level 0 Castro::construct_new_gravity() time = 0.096864791 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.350350804e+43 ener added: 1.765529791e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.933604547e+35 ymom added: 1.827555829e+24 zmom added: 0 eint added: 0 ener added: -1.838318421e+43 Castro::do_new_sources() time = 0.061563861 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.581215356e+29 ener added: 1.581215356e+29 TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -4.341809003e+38 TIME= 0.25 YMOM = 1.164437349e+28 TIME= 0.25 ZMOM = 0 TIME= 0.25 ANG MOM X = 0 TIME= 0.25 ANG MOM Y = 0 TIME= 0.25 ANG MOM Z = -1.469212319e+38 TIME= 0.25 RHO*e = 1.325894512e+50 TIME= 0.25 RHO*K = 3.170941987e+49 TIME= 0.25 RHO*E = 1.642988711e+50 TIME= 0.25 RHO*PHI = -5.449011378e+50 TIME= 0.25 TOTAL ENERGY= -1.081516979e+50 TIME= 0.25 CENTER OF MASS X-LOC = 316424388.8 TIME= 0.25 CENTER OF MASS X-VEL = -170267.2489 TIME= 0.25 CENTER OF MASS Y-LOC = 2.427336677e-06 TIME= 0.25 CENTER OF MASS Y-VEL = 4.566427124e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 57600050.72 TIME= 0.25 MAXIMUM DENSITY = 4458760.788 TIME= 0.25 MAXIMUM T_S / T_E = 0 STEP = 1 TIME = 0.25 DT = 0.25 [STEP 1] Coarse TimeStep time: 12.44487221 [STEP 1] FAB kilobyte spread across MPI nodes: [6998 ... 6998] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.049795497 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004169053 seconds Ending run at 06:48:00 UTC on 2025-07-25. Run time = 14.90260136 Run time without initialization = 12.49980629 Average number of zones advanced per microsecond: 0.003 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 14.8940 ... 14.8963 ... 14.9026 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.6944 5.8084 5.9190 39.72% Amr::InitAmr() 1 1.485 1.487 1.493 10.02% Castro::computeTemp() 80 0.8395 0.9658 1.141 7.66% FabArray::ParallelCopy_finish() 1429 0.3084 0.5771 0.8391 5.63% amrex::Copy() 724 0.7961 0.8097 0.825 5.54% MultiFab::AddProduct() 3439 0.5494 0.5626 0.6087 4.08% Castro::normalize_species() 79 0.5439 0.5541 0.5651 3.79% Castro::reset_internal_energy(Fab) 160 0.5166 0.53 0.5474 3.67% MLPoisson::Fsmooth() 6144 0.3525 0.3873 0.4702 3.16% Castro::initMFs() 1 0.355 0.3654 0.3817 2.56% FabArray::setVal() 5013 0.3478 0.354 0.3604 2.42% MultiFab::Subtract() 159 0.3355 0.3417 0.3484 2.34% check_for_negative_density() 20 0.0379 0.1774 0.3278 2.20% Gravity::fill_multipole_BCs() 19 0.1997 0.2637 0.2951 1.98% FabArray::ParallelCopy_nowait() 1429 0.2039 0.2182 0.2359 1.58% FabArray::sum() 3439 0.1963 0.2001 0.2042 1.37% Castro::initData() 1 0.003395 0.02831 0.1672 1.12% FillBoundary_finish() 6380 0.09994 0.1316 0.1642 1.10% TagBoxArray::collate() 2 0.0002224 0.03228 0.1583 1.06% FillBoundary_nowait() 6380 0.09272 0.1205 0.1531 1.03% MLCellLinOp::applyBC() 7148 0.06338 0.08067 0.1299 0.87% MLPoisson::Fapply() 986 0.09433 0.09984 0.114 0.76% Castro::construct_new_gravity_source() 17 0.07289 0.07486 0.08094 0.54% FabArray::Saxpy() 91 0.07172 0.07311 0.07485 0.50% FabArray::Xpay() 915 0.05137 0.05543 0.06374 0.43% Castro::construct_old_gravity_source() 20 0.04922 0.05012 0.05095 0.34% MultiFab::contains_nan() 37 0.04586 0.04675 0.04773 0.32% MLMG::addInterpCorrection() 768 0.03363 0.03671 0.04729 0.32% AmrLevel::derive() 18 0.04131 0.04198 0.04293 0.29% StateData::FillBoundary(geom) 143 0 0.01768 0.03438 0.23% amrex::average_down 768 0.02495 0.02743 0.03423 0.23% VisMF::Write(FabArray) 4 0.01404 0.0152 0.02616 0.18% Castro::reset_internal_energy(MultiFab) 80 0.02311 0.02376 0.02428 0.16% MLCellLinOp::defineAuxData() 19 0.01589 0.01669 0.02095 0.14% Castro::enforce_min_density() 79 0.01847 0.01891 0.01929 0.13% FabArray::FillBoundary() 6380 0.01408 0.01459 0.01772 0.12% Amr::writePlotFile() 2 0.006809 0.0157 0.01746 0.12% MLCGSolver::bicgstab 8 0 0.001067 0.01707 0.11% Castro::construct_new_source() 102 0.01607 0.01639 0.01681 0.11% amrex::Add() 128 0.01441 0.01466 0.01497 0.10% Castro::evaluate_source_change() 196 0.0136 0.01429 0.01475 0.10% FabArray::setDomainBndry() 147 0.0004969 0.007677 0.01455 0.10% Castro::wd_update() 2 0.005061 0.01057 0.01405 0.09% FabArray::norminf() 244 0.01063 0.01129 0.01399 0.09% FabArrayBase::getFB() 6380 0.01061 0.01107 0.01371 0.09% FillPatchIterator::Initialize 147 0.01272 0.01309 0.01348 0.09% MLCellLinOp::smooth() 1544 0.01024 0.01077 0.01262 0.08% MLMG::prepareForSolve() 18 0.01024 0.0106 0.01177 0.08% Castro::construct_old_source() 120 0.01019 0.01041 0.01059 0.07% MLMG::ResNormInf() 146 0.003991 0.008177 0.01046 0.07% Gravity::actual_solve_with_mlmg() 19 0.00733 0.007765 0.009378 0.06% MLMG::solve() 18 0.001044 0.004915 0.008569 0.06% FabArrayBase::getCPC() 2005 0.007686 0.007947 0.008301 0.06% Gravity::get_old_grav_vector() 20 0.006422 0.006593 0.00674 0.05% Gravity::get_new_grav_vector() 18 0.006024 0.006192 0.00658 0.04% MLCellLinOp::apply() 986 0.00206 0.002295 0.005037 0.03% Castro::initialize_do_advance() 20 0.004513 0.004721 0.004986 0.03% FabArray::ParallelCopy() 1429 0.003621 0.003803 0.004172 0.03% AmrLevel::FillPatch() 147 0.002287 0.002482 0.002592 0.02% Castro::sum_integrated_quantities() 2 0.001794 0.001965 0.002174 0.01% MLCellLinOp::correctionResidual() 768 0.001899 0.001963 0.002057 0.01% MLMG::mgVcycle() 128 0.001652 0.00173 0.001877 0.01% MLMG:computeResOfCorrection() 768 0.001505 0.001572 0.001714 0.01% Castro::save_data_for_retry() 3 0.001313 0.001343 0.001356 0.01% MLMG::actualBottomSolve() 128 0.000132 0.0002054 0.001247 0.01% Castro::subcycle_advance_ctu() 1 0.0002196 0.0003069 0.00112 0.01% StateDataPhysBCFunct::() 147 0.0007322 0.0008916 0.0009821 0.01% main() 1 0.0006606 0.0007569 0.0009789 0.01% FillPatchIterator::FillFromLevel0() 147 0.0007017 0.0007239 0.0007674 0.01% Gravity::solve_for_phi() 17 0.0006634 0.0006975 0.0007536 0.01% MLMG::mgVcycle_down::3 128 0.0005037 0.0005465 0.0006282 0.00% MLMG::mgVcycle_down::2 128 0.0005416 0.0005796 0.0006229 0.00% MLMG::mgVcycle_down::4 128 0.0004833 0.000509 0.0006088 0.00% MLMG::mgVcycle_down::5 128 0.0005095 0.0005285 0.0005966 0.00% MLMG::mgVcycle_down::0 128 0.0005574 0.0005732 0.0005958 0.00% AmrMesh::MakeNewGrids() 2 0.0005487 0.000561 0.0005809 0.00% MLMG::mgVcycle_down::1 128 0.0005571 0.000571 0.0005786 0.00% Castro::Castro() 1 0.0004141 0.0004658 0.0005621 0.00% FillPatchSingleLevel 147 0.0004665 0.0004849 0.0004959 0.00% MLCellLinOp::solutionResidual() 147 0.0004473 0.0004628 0.0004778 0.00% Castro::clean_state() 79 0.0004272 0.0004498 0.0004672 0.00% Castro::do_advance_ctu() 20 0.0003985 0.0004174 0.0004333 0.00% MLMG::oneIter() 128 0.0003944 0.0004077 0.000427 0.00% MLMG::mgVcycle_up::1 128 0.0003605 0.000382 0.0004092 0.00% MLMG::mgVcycle_up::0 128 0.0003724 0.000387 0.0004072 0.00% MLMG::mgVcycle_up::2 128 0.0003293 0.0003538 0.0003873 0.00% MLMG::mgVcycle_up::3 128 0.0003197 0.0003354 0.0003849 0.00% MLMG::mgVcycle_up::4 128 0.0003235 0.0003365 0.0003834 0.00% MLMG::mgVcycle_up::5 128 0.0003102 0.0003251 0.0003597 0.00% Castro::do_old_sources() 20 0.0002784 0.0002952 0.0003153 0.00% Castro::construct_new_gravity() 17 0.0001617 0.0001812 0.0003144 0.00% MLMG::computeResidual() 128 0.0002379 0.0002467 0.0002622 0.00% Castro::problem_post_timestep() 2 1.129e-05 2.766e-05 0.0002602 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002221 0.000238 0.0002545 0.00% MLMG::mgVcycle_bottom 128 0.0002289 0.0002347 0.0002506 0.00% Castro::do_new_sources() 17 0.0002037 0.0002249 0.0002411 0.00% MLPoisson::define() 19 0.0001401 0.0001692 0.0002049 0.00% TagBoxArray::mapPRD 2 0.0001246 0.0001367 0.0001884 0.00% Amr::coarseTimeStep() 1 0.0001283 0.0001518 0.0001585 0.00% Gravity::test_composite_phi() 1 0.0001055 0.0001105 0.0001499 0.00% Castro::retry_advance_ctu() 11 6.502e-05 7.386e-05 0.0001393 0.00% Castro::construct_old_gravity() 20 0.0001084 0.0001197 0.0001361 0.00% Castro::check_for_nan() 37 0.000106 0.0001155 0.0001316 0.00% Castro::initialize_advance() 1 0.0001068 0.0001105 0.0001176 0.00% Amr::defBaseLevel() 1 8.599e-05 0.0001011 0.0001138 0.00% Castro::apply_source_to_state() 37 9.723e-05 0.0001025 0.0001101 0.00% Amr::FinalizeInit() 1 8.23e-06 2.129e-05 0.0001047 0.00% Amr::grid_places() 2 5.364e-05 6.607e-05 8.794e-05 0.00% MLLinOp::define() 19 7.567e-05 8.046e-05 8.459e-05 0.00% Castro::create_source_corrector() 17 5.942e-05 6.576e-05 7.588e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.948e-05 6.148e-05 6.401e-05 0.00% Castro::expand_state() 20 5.089e-05 5.332e-05 6.188e-05 0.00% MLMG::computeMLResidual() 18 4.047e-05 4.432e-05 5.003e-05 0.00% Amr::regrid() 1 1.593e-05 2.078e-05 4.705e-05 0.00% Castro::derive() 18 4.288e-05 4.416e-05 4.654e-05 0.00% Gravity::actual_multilevel_solve() 1 3.143e-05 3.281e-05 3.55e-05 0.00% Castro::advance() 1 7.167e-06 9.916e-06 3.092e-05 0.00% Amr::init() 1 5.782e-06 8.022e-06 2.64e-05 0.00% Amr::timeStep() 1 1.452e-05 1.761e-05 2.313e-05 0.00% Amr::InitializeInit() 1 7.735e-06 1.346e-05 1.959e-05 0.00% Castro::post_regrid() 2 1.138e-05 1.254e-05 1.476e-05 0.00% Amr::initialInit() 1 2.917e-06 3.892e-06 1.474e-05 0.00% Castro::post_init() 1 1.197e-05 1.299e-05 1.411e-05 0.00% Amr::bldFineLevels() 1 7.237e-06 8.091e-06 1.081e-05 0.00% Castro::post_timestep() 1 6.065e-06 6.485e-06 6.814e-06 0.00% Other 2521 0.04947 0.06121 0.09912 0.67% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.89 14.9 14.9 100.00% Amr::coarseTimeStep() 1 12.5 12.5 12.5 83.87% Amr::timeStep() 1 12.44 12.44 12.44 83.51% Castro::advance() 1 12.38 12.38 12.38 83.10% Castro::subcycle_advance_ctu() 1 12.33 12.33 12.33 82.76% Castro::do_advance_ctu() 20 12.29 12.29 12.3 82.51% Castro::construct_ctu_hydro_source() 20 6.73 6.773 6.827 45.81% Castro::clean_state() 79 3.653 3.828 4.047 27.15% Castro::computeTemp() 80 2.34 2.492 2.69 18.05% Gravity::solve_phi_with_mlmg() 19 1.625 1.693 1.728 11.60% Castro::construct_new_gravity() 17 1.542 1.61 1.644 11.03% Gravity::solve_for_phi() 17 1.504 1.572 1.607 10.78% Castro::reset_internal_energy(MultiFab) 80 1.494 1.526 1.559 10.46% Amr::InitAmr() 1 1.485 1.487 1.493 10.02% Castro::initialize_do_advance() 20 1.41 1.45 1.485 9.97% Gravity::actual_solve_with_mlmg() 19 1.422 1.426 1.43 9.59% MLMG::solve() 18 1.366 1.372 1.374 9.22% MLMG::oneIter() 128 1.229 1.231 1.233 8.27% Castro::do_old_sources() 20 1.143 1.194 1.228 8.24% MLMG::mgVcycle() 128 1.214 1.216 1.218 8.17% FabArray::ParallelCopy() 1429 0.5491 0.8197 1.079 7.24% Castro::do_new_sources() 17 0.9215 0.9845 1.044 7.00% Castro::evaluate_source_change() 196 0.9557 0.9769 1.016 6.82% Amr::init() 1 0.9093 0.9094 0.9095 6.10% Amr::initialInit() 1 0.8551 0.8552 0.8553 5.74% FabArray::ParallelCopy_finish() 1429 0.3084 0.5771 0.8391 5.63% MLCellLinOp::smooth() 1544 0.6346 0.6887 0.8323 5.59% amrex::Copy() 724 0.7961 0.8097 0.825 5.54% Castro::enforce_min_density() 79 0.7939 0.8091 0.8239 5.53% AmrLevel::FillPatch() 147 0.51 0.627 0.7431 4.99% Amr::InitializeInit() 1 0.4414 0.4781 0.6195 4.16% Amr::defBaseLevel() 1 0.4414 0.4781 0.6195 4.16% MultiFab::AddProduct() 3439 0.5494 0.5626 0.6087 4.08% Castro::normalize_species() 79 0.5439 0.5541 0.5651 3.79% FillPatchIterator::Initialize 147 0.3228 0.4423 0.5615 3.77% Castro::reset_internal_energy(Fab) 160 0.5166 0.53 0.5474 3.67% FillPatchIterator::FillFromLevel0() 147 0.295 0.4215 0.5417 3.63% FillPatchSingleLevel 147 0.2943 0.4208 0.541 3.63% MLPoisson::Fsmooth() 6144 0.3525 0.3873 0.4702 3.16% MLMG::mgVcycle_up::1 128 0.1064 0.2872 0.4699 3.15% MLCellLinOp::applyBC() 7148 0.3118 0.3591 0.4389 2.95% Amr::FinalizeInit() 1 0.2357 0.3771 0.4138 2.78% MLMG::addInterpCorrection() 768 0.07 0.3308 0.4082 2.74% Castro::Castro() 1 0.3604 0.3713 0.3885 2.61% Castro::initMFs() 1 0.3552 0.3655 0.3818 2.56% FabArray::setVal() 5013 0.3478 0.354 0.3604 2.42% MultiFab::Subtract() 159 0.3355 0.3417 0.3484 2.34% check_for_negative_density() 20 0.0379 0.1774 0.3278 2.20% FabArray::FillBoundary() 6380 0.2401 0.2785 0.3222 2.16% MLMG::mgVcycle_up::2 128 0.01374 0.1299 0.3072 2.06% Gravity::fill_multipole_BCs() 19 0.2018 0.2659 0.2972 1.99% Castro::expand_state() 20 0.1995 0.2323 0.2879 1.93% MLMG::mgVcycle_down::0 128 0.2806 0.2828 0.2868 1.92% FabArray::ParallelCopy_nowait() 1429 0.2241 0.2388 0.2573 1.73% Castro::initData() 1 0.07886 0.1048 0.2443 1.64% Amr::bldFineLevels() 1 0.04696 0.188 0.2242 1.50% MLMG::mgVcycle_up::0 128 0.2153 0.2167 0.2189 1.47% MLCellLinOp::apply() 986 0.1569 0.1686 0.2066 1.39% FabArray::sum() 3439 0.1963 0.2001 0.2042 1.37% MLMG::mgVcycle_up::3 128 0.01366 0.02811 0.1993 1.34% Amr::grid_places() 2 0.002527 0.144 0.1807 1.21% AmrMesh::MakeNewGrids() 2 0.002473 0.1439 0.1806 1.21% TagBoxArray::mapPRD 2 0.0008397 0.1103 0.1772 1.19% FillBoundary_nowait() 6380 0.1039 0.1323 0.1653 1.11% FillBoundary_finish() 6380 0.09994 0.1316 0.1642 1.10% TagBoxArray::collate() 2 0.0002224 0.03228 0.1583 1.06% MLMG:computeResOfCorrection() 768 0.1192 0.1294 0.1546 1.04% MLCellLinOp::correctionResidual() 768 0.1177 0.1279 0.1529 1.03% Castro::post_init() 1 0.1441 0.1455 0.147 0.99% Castro::construct_new_source() 102 0.1291 0.132 0.1394 0.94% MLMG::mgVcycle_down::1 128 0.1254 0.1308 0.135 0.91% Castro::create_source_corrector() 17 0.05183 0.07753 0.1256 0.84% MLPoisson::Fapply() 986 0.09433 0.09984 0.114 0.76% Gravity::multilevel_solve_for_new_phi() 1 0.11 0.11 0.1101 0.74% Gravity::actual_multilevel_solve() 1 0.1099 0.11 0.11 0.74% MLCellLinOp::solutionResidual() 147 0.09169 0.09685 0.1023 0.69% Amr::writePlotFile() 2 0.09975 0.09985 0.09987 0.67% Castro::post_regrid() 2 0.08614 0.08762 0.08902 0.60% Castro::construct_old_source() 120 0.08294 0.08445 0.08581 0.58% MLMG::mgVcycle_down::2 128 0.01601 0.04752 0.08164 0.55% Castro::construct_new_gravity_source() 17 0.07289 0.07486 0.08094 0.54% MLMG::computeResidual() 128 0.07844 0.07934 0.08087 0.54% FabArray::Saxpy() 91 0.07172 0.07311 0.07485 0.50% MLMG::mgVcycle_bottom 128 0.000363 0.004995 0.07437 0.50% MLMG::actualBottomSolve() 128 0.000132 0.00476 0.07412 0.50% Castro::construct_old_gravity() 20 0.04161 0.0561 0.0738 0.50% Gravity::get_old_grav_vector() 20 0.04149 0.05598 0.07367 0.49% Castro::derive() 18 0.07077 0.07183 0.07336 0.49% AmrLevel::derive() 18 0.07072 0.07179 0.07331 0.49% MLCGSolver::bicgstab 8 0 0.004513 0.07221 0.48% MLMG::mgVcycle_down::3 128 0.01532 0.02233 0.06946 0.47% Castro::post_timestep() 1 0.06128 0.06404 0.06695 0.45% FabArray::Xpay() 915 0.05137 0.05543 0.06374 0.43% amrex::average_down 768 0.03399 0.0414 0.05955 0.40% Castro::construct_old_gravity_source() 20 0.04922 0.05012 0.05095 0.34% Castro::initialize_advance() 1 0.04911 0.05003 0.05084 0.34% Castro::check_for_nan() 37 0.04597 0.04687 0.04785 0.32% MultiFab::contains_nan() 37 0.04586 0.04675 0.04773 0.32% Amr::regrid() 1 0.04446 0.04522 0.04602 0.31% MLPoisson::define() 19 0.03066 0.03569 0.04131 0.28% MLMG::mgVcycle_down::4 128 0.01524 0.01717 0.04108 0.28% Gravity::get_new_grav_vector() 18 0.03623 0.03711 0.03784 0.25% Castro::retry_advance_ctu() 11 0.03612 0.03706 0.0377 0.25% StateDataPhysBCFunct::() 147 0.0007322 0.01857 0.03532 0.24% StateData::FillBoundary(geom) 143 0 0.01768 0.03438 0.23% MLMG::mgVcycle_up::4 128 0.01372 0.01519 0.03302 0.22% Castro::apply_source_to_state() 37 0.02927 0.02982 0.03072 0.21% MLMG::mgVcycle_down::5 128 0.01523 0.0164 0.02895 0.19% MLCellLinOp::defineAuxData() 19 0.0202 0.02173 0.02816 0.19% VisMF::Write(FabArray) 4 0.01404 0.0152 0.02616 0.18% MLMG::mgVcycle_up::5 128 0.01384 0.01479 0.02483 0.17% MLMG::prepareForSolve() 18 0.02061 0.02165 0.02456 0.16% MLMG::computeMLResidual() 18 0.01105 0.01694 0.02266 0.15% Castro::sum_integrated_quantities() 2 0.02171 0.02187 0.02189 0.15% MLMG::ResNormInf() 146 0.01371 0.01771 0.01985 0.13% Castro::save_data_for_retry() 3 0.01923 0.01953 0.01979 0.13% Castro::problem_post_timestep() 2 0.01086 0.01435 0.01743 0.12% Castro::wd_update() 2 0.01084 0.01432 0.01742 0.12% FabArrayBase::getCPC() 2005 0.01014 0.01151 0.01539 0.10% FabArrayBase::getFB() 6380 0.01095 0.01181 0.01532 0.10% amrex::Add() 128 0.01441 0.01466 0.01497 0.10% Gravity::test_composite_phi() 1 0.01458 0.01463 0.01466 0.10% FabArray::setDomainBndry() 147 0.0004969 0.007677 0.01455 0.10% FabArray::norminf() 244 0.01063 0.01129 0.01399 0.09% MLLinOp::define() 19 0.003709 0.007341 0.01104 0.07% Other 2521 0.08211 0.09907 0.1483 0.99% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Pinned Memory Usage: --------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max --------------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 16 9 B 12 B 22 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- Cpu Memory Usage: ------------------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max ------------------------------------------------------------------------------------------------------------------------------- Castro::save_data_for_retry() 416 416 1559 KiB 1560 KiB 1560 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 401 KiB 409 KiB 416 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1033 KiB 1033 KiB 1034 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 983 KiB 983 KiB 984 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 660 KiB 660 KiB 660 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 4232 B 4324 B 4820 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 49 KiB 50 KiB 51 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 11 KiB 13 KiB 15 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 413 KiB 413 KiB 413 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 1973 B 3960 B 5892 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1317 B 1337 B 1365 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 16 KiB 17 KiB 17 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 131 KiB 131 KiB 131 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 192 127 KiB 127 KiB 128 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 6102 9639 B 9 KiB 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 252 B 256 B 261 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 426 B 433 B 443 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 18 B 20 B 36 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 33 KiB 33 KiB 33 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 51 B 51 B 51 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 3 B 4 B 6 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 110 B 112 B 115 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 127 B 171 B 226 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 64 31 KiB 31 KiB 31 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 2432 225 B 234 B 240 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1046 B 1068 B 1112 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 120 B 120 B 120 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1653 B 1727 B 1766 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 32 0 B 0 B 1 B 16 KiB 16 KiB 16 KiB MLCellLinOp::defineAuxData() 11704 11704 435 B 486 B 675 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 452 B 454 B 455 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 23 B 33 B 46 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 18 B 31 B 2916 B 2916 B 2916 B OwnerMask() 64 64 0 B 0 B 0 B 2592 B 2592 B 2592 B MLMG::addInterpCorrection() 5632 5632 1 B 16 B 25 B 1024 B 1024 B 1024 B amrex::average_down 5632 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 64 0 B 5 B 7 B 648 B 648 B 648 B MLCGSolver::bicgstab 640 640 0 B 0 B 2 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 114375 0 B 0 B 0 B 8 B 8 B 8 B ------------------------------------------------------------------------------------------------------------------------------- AMReX (25.06-22-g1d307c2c661a-dirty) finalized