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 17:50:45 UTC on 2025-08-04. Successfully read inputs file ... Castro git describe: 25.07-101-g9a9d15b8d AMReX git describe: 25.06-22-g1d307c2c6-dirty Microphysics git describe: joss_submit-1-ga350de6c 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.172877101 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001180598 new finest: 0 Gravity::fill_multipole_BCs() time = 0.011575711 Gravity::multilevel_solve_for_new_phi() time = 0.160954399 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010448783 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.050304249 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004062739 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 = -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 = -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 = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 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 = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 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 = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188Invalid 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 , 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 = Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 -6.239071148 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 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, 69, 0 Invalid density = -908.1731819 at index 14, 70, 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, 0Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 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.005711224 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.099769454 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.012346993 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.072568355 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.005217882 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.06091226 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.007931073 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.062017345 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.351060685 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.014371497 Gravity::solve_for_phi() time = 0.091288894 on level 0 Castro::construct_new_gravity() time = 0.093663857 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.059851749 on level 0 Castro::construct_old_gravity() time = 0.002960994 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.059621671 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.355754434 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.014790735 Gravity::solve_for_phi() time = 0.092329841 on level 0 Castro::construct_new_gravity() time = 0.094711996 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.059675724 on level 0 Castro::construct_old_gravity() time = 0.004757347 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.060159835 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.357041001 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.016275756 Gravity::solve_for_phi() time = 0.094511685 on level 0 Castro::construct_new_gravity() time = 0.096828754 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.060751508 on level 0 Castro::construct_old_gravity() time = 0.00401163 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.060114068 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.358024059 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.014861153 Gravity::solve_for_phi() time = 0.092008734 on level 0 Castro::construct_new_gravity() time = 0.09433874 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.059935804 on level 0 Castro::construct_old_gravity() time = 0.003781472 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.060224269 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.37099105 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.015296833 Gravity::solve_for_phi() time = 0.09035856 on level 0 Castro::construct_new_gravity() time = 0.092717239 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.059308982 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.005320006 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.060176841 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.385862748 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.014716716 Gravity::solve_for_phi() time = 0.089874175 on level 0 Castro::construct_new_gravity() time = 0.092248176 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.060916633 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.004683859 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.061537579 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.366975566 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.015566907 Gravity::solve_for_phi() time = 0.092350618 on level 0 Castro::construct_new_gravity() time = 0.094688212 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.061304461 on level 0 Castro::construct_old_gravity() time = 0.003138457 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.061133029 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.363911095 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.016809018 Gravity::solve_for_phi() time = 0.092915428 on level 0 Castro::construct_new_gravity() time = 0.09524384 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.060970519 on level 0 Castro::construct_old_gravity() time = 0.003254649 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.061460366 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.368635933 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.016196344 Gravity::solve_for_phi() time = 0.092204754 on level 0 Castro::construct_new_gravity() time = 0.094528931 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.061013956 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.003240983 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.061472362 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.363130298 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.016887922 Gravity::solve_for_phi() time = 0.093494358 on level 0 Castro::construct_new_gravity() time = 0.095831144 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.061367971 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.003205143 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.061401884 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.370282602 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.016827043 Gravity::solve_for_phi() time = 0.092979486 on level 0 Castro::construct_new_gravity() time = 0.095288096 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.061285534 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.00312362 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.060998167 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.363428104 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.016806507 Gravity::solve_for_phi() time = 0.093117249 on level 0 Castro::construct_new_gravity() time = 0.095518423 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.061150721 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.004829924 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.061417043 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.369614247 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.016596324 Gravity::solve_for_phi() time = 0.092399567 on level 0 Castro::construct_new_gravity() time = 0.094777499 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.061551072 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.003147084 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.061165463 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.363734332 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.015993845 Gravity::solve_for_phi() time = 0.091877856 on level 0 Castro::construct_new_gravity() time = 0.094227938 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.061152938 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.003290514 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.061050827 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.368262578 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.016274415 Gravity::solve_for_phi() time = 0.092434131 on level 0 Castro::construct_new_gravity() time = 0.094761605 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.061496805 on level 0 Castro::construct_old_gravity() time = 0.003088951 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.06148587 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.363675935 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.017267705 Gravity::solve_for_phi() time = 0.093511782 on level 0 Castro::construct_new_gravity() time = 0.095878665 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.061102135 on level 0 Castro::construct_old_gravity() time = 0.003240748 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.060991362 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.369335798 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002575422407 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.016499058 Gravity::solve_for_phi() time = 0.092222297 on level 0 Castro::construct_new_gravity() time = 0.094557071 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.061354085 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.78656176 [STEP 1] FAB kilobyte spread across MPI nodes: [7120 ... 7120] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.049078645 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.00407594 seconds Ending run at 17:51:01 UTC on 2025-08-04. Run time = 15.6035423 Run time without initialization = 12.84081847 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]: 15.5959 ... 15.5984 ... 15.6036 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.9751 6.0592 6.1739 39.57% Amr::InitAmr() 1 1.799 1.802 1.807 11.58% Castro::computeTemp() 80 0.8376 0.9602 1.127 7.23% amrex::Copy() 724 0.7946 0.8043 0.8178 5.24% FabArray::ParallelCopy_finish() 1429 0.301 0.6013 0.8054 5.16% MultiFab::AddProduct() 3439 0.5638 0.5792 0.6922 4.44% Castro::normalize_species() 79 0.5528 0.5746 0.5907 3.79% Castro::reset_internal_energy(Fab) 160 0.5202 0.5282 0.544 3.49% MLPoisson::Fsmooth() 6144 0.3503 0.3829 0.4545 2.91% Castro::initMFs() 1 0.3567 0.3663 0.3773 2.42% MultiFab::Subtract() 159 0.3517 0.3563 0.3613 2.32% FabArray::setVal() 5013 0.3418 0.3487 0.3586 2.30% check_for_negative_density() 20 0.1105 0.2435 0.3419 2.19% Gravity::fill_multipole_BCs() 19 0.2037 0.2605 0.2906 1.86% FabArray::ParallelCopy_nowait() 1429 0.2153 0.2332 0.258 1.65% FillBoundary_finish() 6380 0.1283 0.1627 0.2082 1.33% FabArray::sum() 3439 0.175 0.1775 0.1834 1.18% Castro::initData() 1 0.003405 0.02777 0.1611 1.03% FillBoundary_nowait() 6380 0.09283 0.1208 0.1577 1.01% TagBoxArray::collate() 2 0.0002302 0.03033 0.1447 0.93% MLCellLinOp::applyBC() 7148 0.06301 0.0797 0.1251 0.80% MLPoisson::Fapply() 986 0.09305 0.09877 0.1118 0.72% Castro::construct_new_gravity_source() 17 0.0723 0.07339 0.07501 0.48% FabArray::Saxpy() 91 0.06794 0.06876 0.07012 0.45% FabArray::Xpay() 915 0.05032 0.05308 0.05835 0.37% Castro::construct_old_gravity_source() 20 0.04853 0.04902 0.04983 0.32% MultiFab::contains_nan() 37 0.04501 0.0457 0.04648 0.30% MLMG::addInterpCorrection() 768 0.03381 0.03657 0.0438 0.28% AmrLevel::derive() 18 0.04078 0.04133 0.04216 0.27% StateData::FillBoundary(geom) 143 0 0.01941 0.03779 0.24% amrex::average_down 768 0.02518 0.02756 0.0338 0.22% VisMF::Write(FabArray) 4 0.01411 0.01539 0.02757 0.18% Castro::reset_internal_energy(MultiFab) 80 0.02338 0.02376 0.02433 0.16% MLCellLinOp::defineAuxData() 19 0.01578 0.01653 0.02009 0.13% MLMG::ResNormInf() 146 0.007234 0.01463 0.0195 0.12% Castro::enforce_min_density() 79 0.01829 0.01852 0.01876 0.12% Amr::writePlotFile() 2 0.006633 0.01671 0.01852 0.12% FabArray::FillBoundary() 6380 0.01365 0.01434 0.01725 0.11% FabArray::setDomainBndry() 147 0.0004898 0.008548 0.01629 0.10% Castro::construct_new_source() 102 0.01588 0.01608 0.01628 0.10% Castro::wd_update() 2 0.005668 0.01088 0.01627 0.10% MLCGSolver::bicgstab 8 0 0.001009 0.01614 0.10% amrex::Add() 128 0.01446 0.01472 0.01589 0.10% Castro::evaluate_source_change() 196 0.01377 0.01415 0.01489 0.10% FabArray::norminf() 244 0.0109 0.01142 0.01411 0.09% FabArrayBase::getFB() 6380 0.01057 0.0111 0.01352 0.09% FillPatchIterator::Initialize 147 0.01295 0.01319 0.01349 0.09% MLMG::prepareForSolve() 18 0.01011 0.01052 0.01149 0.07% MLCellLinOp::smooth() 1544 0.01012 0.01058 0.01139 0.07% Castro::construct_old_source() 120 0.01013 0.01025 0.01042 0.07% Gravity::actual_solve_with_mlmg() 19 0.007341 0.007715 0.009 0.06% FabArrayBase::getCPC() 2005 0.007553 0.007784 0.008057 0.05% Gravity::get_old_grav_vector() 20 0.006587 0.006748 0.007115 0.05% MLMG::solve() 18 0.0008836 0.003416 0.007075 0.05% Gravity::get_new_grav_vector() 18 0.006183 0.006298 0.006734 0.04% Castro::initialize_do_advance() 20 0.00448 0.004691 0.004899 0.03% MLCellLinOp::apply() 986 0.001985 0.002225 0.004689 0.03% FabArray::ParallelCopy() 1429 0.003557 0.003732 0.004078 0.03% AmrLevel::FillPatch() 147 0.002443 0.002521 0.002582 0.02% Castro::sum_integrated_quantities() 2 0.001801 0.001965 0.002103 0.01% MLCellLinOp::correctionResidual() 768 0.001941 0.001994 0.002067 0.01% MLMG::mgVcycle() 128 0.001647 0.001713 0.001834 0.01% MLMG:computeResOfCorrection() 768 0.001489 0.001544 0.001624 0.01% Castro::save_data_for_retry() 3 0.001342 0.001388 0.001418 0.01% MLMG::actualBottomSolve() 128 0.000129 0.0002019 0.001215 0.01% Castro::subcycle_advance_ctu() 1 0.0002185 0.0002971 0.001115 0.01% StateDataPhysBCFunct::() 147 0.0007336 0.0008925 0.0009867 0.01% main() 1 0.000679 0.000756 0.0009497 0.01% Gravity::solve_for_phi() 17 0.0006473 0.0006857 0.0007803 0.01% FillPatchIterator::FillFromLevel0() 147 0.0006826 0.0007122 0.0007433 0.00% AmrMesh::MakeNewGrids() 2 0.0005533 0.0005814 0.0006277 0.00% MLMG::mgVcycle_down::2 128 0.0005401 0.0005817 0.0006264 0.00% MLMG::mgVcycle_down::3 128 0.0005011 0.0005427 0.0006167 0.00% MLMG::mgVcycle_down::1 128 0.0005604 0.0005736 0.0006044 0.00% MLMG::mgVcycle_down::4 128 0.0004806 0.0005065 0.0005937 0.00% MLMG::mgVcycle_down::0 128 0.0005549 0.0005684 0.0005862 0.00% MLMG::mgVcycle_down::5 128 0.0005043 0.0005223 0.0005824 0.00% TagBoxArray::mapPRD 2 0.0001217 0.0001527 0.0005211 0.00% Castro::Castro() 1 0.0004057 0.0004305 0.0005205 0.00% FillPatchSingleLevel 147 0.0004642 0.0004768 0.00049 0.00% Castro::clean_state() 79 0.0004299 0.0004602 0.0004839 0.00% MLCellLinOp::solutionResidual() 147 0.000452 0.0004657 0.0004819 0.00% Castro::do_advance_ctu() 20 0.0003987 0.000414 0.0004507 0.00% MLMG::oneIter() 128 0.000389 0.0004168 0.0004473 0.00% MLMG::mgVcycle_up::0 128 0.0003705 0.0003864 0.0004087 0.00% MLMG::mgVcycle_up::2 128 0.0003239 0.0003507 0.0004001 0.00% MLMG::mgVcycle_up::1 128 0.0003629 0.0003789 0.0003972 0.00% MLMG::mgVcycle_up::3 128 0.0003182 0.0003338 0.0003771 0.00% MLMG::mgVcycle_up::4 128 0.0003226 0.0003346 0.000363 0.00% Amr::timeStep() 1 1.431e-05 3.718e-05 0.0003428 0.00% MLMG::mgVcycle_up::5 128 0.0003105 0.0003208 0.0003406 0.00% Castro::construct_new_gravity() 17 0.0001587 0.0001778 0.0003386 0.00% Castro::do_old_sources() 20 0.0002817 0.0002912 0.000303 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002221 0.0002341 0.0002665 0.00% MLMG::computeResidual() 128 0.0002384 0.0002444 0.0002546 0.00% Castro::do_new_sources() 17 0.0002079 0.0002206 0.0002472 0.00% MLMG::mgVcycle_bottom 128 0.0002186 0.0002277 0.0002469 0.00% Castro::problem_post_timestep() 2 1.181e-05 2.737e-05 0.000241 0.00% Amr::coarseTimeStep() 1 0.0001408 0.0001646 0.0001695 0.00% Amr::defBaseLevel() 1 0.0001033 0.00012 0.0001683 0.00% MLPoisson::define() 19 0.0001284 0.0001362 0.0001518 0.00% Gravity::test_composite_phi() 1 0.0001056 0.0001123 0.0001399 0.00% Castro::construct_old_gravity() 20 0.0001102 0.0001236 0.0001376 0.00% Castro::retry_advance_ctu() 11 6.888e-05 7.698e-05 0.0001348 0.00% Castro::check_for_nan() 37 0.0001109 0.0001186 0.0001305 0.00% Castro::initialize_advance() 1 0.0001098 0.0001154 0.0001252 0.00% Castro::apply_source_to_state() 37 9.693e-05 0.0001032 0.0001199 0.00% Castro::create_source_corrector() 17 5.839e-05 6.831e-05 9.288e-05 0.00% Amr::FinalizeInit() 1 8.666e-06 1.443e-05 8.88e-05 0.00% MLLinOp::define() 19 7.513e-05 7.971e-05 8.432e-05 0.00% Amr::grid_places() 2 5.25e-05 6.667e-05 8.097e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.815e-05 6.043e-05 6.33e-05 0.00% Castro::expand_state() 20 5.195e-05 5.545e-05 6.025e-05 0.00% Castro::derive() 18 4.235e-05 4.351e-05 4.607e-05 0.00% MLMG::computeMLResidual() 18 4.114e-05 4.306e-05 4.527e-05 0.00% Amr::regrid() 1 1.64e-05 2.077e-05 4.266e-05 0.00% Gravity::actual_multilevel_solve() 1 3.144e-05 3.275e-05 3.476e-05 0.00% Castro::post_init() 1 1.672e-05 1.881e-05 2.269e-05 0.00% Amr::InitializeInit() 1 9.294e-06 1.253e-05 1.817e-05 0.00% Castro::advance() 1 6.517e-06 8.157e-06 1.598e-05 0.00% Castro::post_regrid() 2 1.086e-05 1.157e-05 1.27e-05 0.00% Amr::bldFineLevels() 1 6.74e-06 8.106e-06 1.155e-05 0.00% Castro::post_timestep() 1 5.735e-06 6.489e-06 1.144e-05 0.00% Amr::init() 1 6.088e-06 7.107e-06 8.365e-06 0.00% Amr::initialInit() 1 2.964e-06 3.354e-06 3.952e-06 0.00% Other 2520 0.05062 0.06399 0.09879 0.63% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 15.6 15.6 15.6 100.00% Amr::coarseTimeStep() 1 12.84 12.84 12.84 82.29% Amr::timeStep() 1 12.79 12.79 12.79 81.95% Castro::advance() 1 12.72 12.72 12.73 81.56% Castro::subcycle_advance_ctu() 1 12.67 12.67 12.68 81.24% Castro::do_advance_ctu() 20 12.63 12.63 12.64 80.99% Castro::construct_ctu_hydro_source() 20 7.056 7.092 7.152 45.84% Castro::clean_state() 79 3.661 3.838 4.019 25.76% Castro::computeTemp() 80 2.339 2.482 2.657 17.03% Amr::InitAmr() 1 1.799 1.802 1.807 11.58% Gravity::solve_phi_with_mlmg() 19 1.652 1.713 1.739 11.15% Castro::construct_new_gravity() 17 1.522 1.582 1.607 10.30% Castro::reset_internal_energy(MultiFab) 80 1.5 1.522 1.568 10.05% Gravity::solve_for_phi() 17 1.481 1.541 1.568 10.05% Castro::initialize_do_advance() 20 1.414 1.456 1.504 9.64% Gravity::actual_solve_with_mlmg() 19 1.445 1.449 1.454 9.32% MLMG::solve() 18 1.39 1.395 1.397 8.95% Castro::do_old_sources() 20 1.191 1.222 1.258 8.06% MLMG::oneIter() 128 1.245 1.248 1.252 8.02% MLMG::mgVcycle() 128 1.23 1.233 1.236 7.92% Castro::evaluate_source_change() 196 0.9457 0.9662 1.08 6.92% FabArray::ParallelCopy() 1429 0.5555 0.8584 1.059 6.79% Castro::do_new_sources() 17 0.9188 0.981 1.033 6.62% Amr::init() 1 0.9553 0.9554 0.9555 6.12% Amr::initialInit() 1 0.9005 0.9006 0.9007 5.77% Castro::enforce_min_density() 79 0.7944 0.8086 0.8664 5.55% MLCellLinOp::smooth() 1544 0.6529 0.7092 0.8537 5.47% amrex::Copy() 724 0.7946 0.8043 0.8178 5.24% AmrLevel::FillPatch() 147 0.586 0.6791 0.8056 5.16% FabArray::ParallelCopy_finish() 1429 0.301 0.6013 0.8054 5.16% MultiFab::AddProduct() 3439 0.5638 0.5792 0.6922 4.44% FillPatchIterator::Initialize 147 0.3981 0.4917 0.6205 3.98% Amr::InitializeInit() 1 0.4423 0.4783 0.6138 3.93% Amr::defBaseLevel() 1 0.4423 0.4783 0.6137 3.93% FillPatchIterator::FillFromLevel0() 147 0.3684 0.4699 0.5915 3.79% FillPatchSingleLevel 147 0.3677 0.4692 0.5908 3.79% Castro::normalize_species() 79 0.5528 0.5746 0.5907 3.79% Castro::reset_internal_energy(Fab) 160 0.5202 0.5282 0.544 3.49% MLCellLinOp::applyBC() 7148 0.3493 0.3894 0.482 3.09% MLMG::mgVcycle_up::1 128 0.1101 0.2901 0.4699 3.01% Amr::FinalizeInit() 1 0.2868 0.4223 0.4583 2.94% MLPoisson::Fsmooth() 6144 0.3503 0.3829 0.4545 2.91% MLMG::addInterpCorrection() 768 0.06625 0.3263 0.4054 2.60% Castro::Castro() 1 0.3622 0.372 0.3834 2.46% Castro::initMFs() 1 0.3569 0.3665 0.3775 2.42% MultiFab::Subtract() 159 0.3517 0.3563 0.3613 2.32% FabArray::setVal() 5013 0.3418 0.3487 0.3586 2.30% FabArray::FillBoundary() 6380 0.2729 0.3097 0.3569 2.29% check_for_negative_density() 20 0.1105 0.2435 0.3419 2.19% MLMG::mgVcycle_up::2 128 0.01355 0.1278 0.302 1.94% Gravity::fill_multipole_BCs() 19 0.2058 0.2626 0.2927 1.88% MLMG::mgVcycle_down::0 128 0.283 0.2863 0.2925 1.87% FabArray::ParallelCopy_nowait() 1429 0.2349 0.2534 0.2793 1.79% Castro::expand_state() 20 0.1963 0.2304 0.2758 1.77% Castro::initData() 1 0.07804 0.1045 0.2359 1.51% MLMG::mgVcycle_up::0 128 0.22 0.2217 0.2238 1.43% Amr::bldFineLevels() 1 0.04555 0.1815 0.2175 1.39% MLCellLinOp::apply() 986 0.1611 0.1727 0.2089 1.34% FillBoundary_finish() 6380 0.1283 0.1627 0.2082 1.33% Castro::post_init() 1 0.1944 0.1976 0.1984 1.27% MLMG::mgVcycle_up::3 128 0.0135 0.02716 0.1881 1.21% FabArray::sum() 3439 0.175 0.1775 0.1834 1.18% Amr::grid_places() 2 0.002486 0.138 0.174 1.12% AmrMesh::MakeNewGrids() 2 0.002434 0.1379 0.174 1.11% TagBoxArray::mapPRD 2 0.0008209 0.1062 0.1723 1.10% FillBoundary_nowait() 6380 0.104 0.1327 0.1705 1.09% Gravity::multilevel_solve_for_new_phi() 1 0.1609 0.1609 0.161 1.03% Gravity::actual_multilevel_solve() 1 0.1608 0.1609 0.1609 1.03% MLMG:computeResOfCorrection() 768 0.1194 0.1304 0.1546 0.99% MLCellLinOp::correctionResidual() 768 0.1179 0.1289 0.153 0.98% MLMG::mgVcycle_down::1 128 0.1364 0.1403 0.1452 0.93% TagBoxArray::collate() 2 0.0002302 0.03033 0.1447 0.93% Castro::construct_new_source() 102 0.1266 0.1282 0.1302 0.83% MLPoisson::Fapply() 986 0.09305 0.09877 0.1118 0.72% Castro::create_source_corrector() 17 0.06365 0.08151 0.105 0.67% MLCellLinOp::solutionResidual() 147 0.08912 0.09766 0.1042 0.67% Amr::writePlotFile() 2 0.09973 0.09982 0.09983 0.64% Castro::post_regrid() 2 0.08582 0.08668 0.0899 0.58% MLMG::mgVcycle_down::2 128 0.01582 0.04814 0.08404 0.54% Castro::construct_old_source() 120 0.08132 0.08219 0.08379 0.54% MLMG::computeResidual() 128 0.07801 0.08048 0.08368 0.54% Castro::construct_old_gravity() 20 0.05404 0.06903 0.0813 0.52% Gravity::get_old_grav_vector() 20 0.05392 0.0689 0.08117 0.52% Castro::construct_new_gravity_source() 17 0.0723 0.07339 0.07501 0.48% Castro::derive() 18 0.06967 0.07058 0.0724 0.46% AmrLevel::derive() 18 0.06962 0.07054 0.07235 0.46% MLMG::mgVcycle_bottom 128 0.00035 0.004776 0.07101 0.46% MLMG::actualBottomSolve() 128 0.000129 0.004549 0.07076 0.45% FabArray::Saxpy() 91 0.06794 0.06876 0.07012 0.45% MLMG::mgVcycle_down::3 128 0.01514 0.02212 0.06894 0.44% MLCGSolver::bicgstab 8 0 0.004308 0.06892 0.44% Castro::post_timestep() 1 0.06054 0.06415 0.0684 0.44% amrex::average_down 768 0.0341 0.04189 0.06163 0.39% FabArray::Xpay() 915 0.05032 0.05308 0.05835 0.37% Castro::initialize_advance() 1 0.0495 0.0502 0.05322 0.34% Castro::construct_old_gravity_source() 20 0.04853 0.04902 0.04983 0.32% Castro::check_for_nan() 37 0.04512 0.04582 0.04659 0.30% MultiFab::contains_nan() 37 0.04501 0.0457 0.04648 0.30% Amr::regrid() 1 0.04421 0.04464 0.04617 0.30% MLPoisson::define() 19 0.03069 0.03595 0.04206 0.27% Gravity::get_new_grav_vector() 18 0.03954 0.04018 0.0407 0.26% MLMG::mgVcycle_down::4 128 0.01509 0.01692 0.03952 0.25% Castro::retry_advance_ctu() 11 0.03642 0.03706 0.03909 0.25% StateDataPhysBCFunct::() 147 0.0007336 0.0203 0.03877 0.25% StateData::FillBoundary(geom) 143 0 0.01941 0.03779 0.24% MLMG::mgVcycle_up::4 128 0.01356 0.01489 0.03131 0.20% MLMG::ResNormInf() 146 0.01768 0.02442 0.02913 0.19% Castro::apply_source_to_state() 37 0.0276 0.02793 0.02852 0.18% MLMG::mgVcycle_down::5 128 0.0151 0.01615 0.02764 0.18% VisMF::Write(FabArray) 4 0.01411 0.01539 0.02757 0.18% MLCellLinOp::defineAuxData() 19 0.01993 0.0215 0.02712 0.17% MLMG::prepareForSolve() 18 0.0203 0.02134 0.02391 0.15% MLMG::mgVcycle_up::5 128 0.01363 0.01451 0.0237 0.15% Castro::sum_integrated_quantities() 2 0.02146 0.02161 0.02164 0.14% MLMG::computeMLResidual() 18 0.01066 0.01668 0.02157 0.14% Castro::save_data_for_retry() 3 0.01937 0.01969 0.02064 0.13% Castro::problem_post_timestep() 2 0.01062 0.01465 0.01965 0.13% Castro::wd_update() 2 0.01061 0.01462 0.01964 0.13% FabArray::setDomainBndry() 147 0.0004898 0.008548 0.01629 0.10% amrex::Add() 128 0.01446 0.01472 0.01589 0.10% FabArrayBase::getFB() 6380 0.01088 0.01184 0.01511 0.10% FabArrayBase::getCPC() 2005 0.009785 0.01124 0.01494 0.10% Gravity::test_composite_phi() 1 0.01411 0.01415 0.01418 0.09% FabArray::norminf() 244 0.0109 0.01142 0.01411 0.09% MLLinOp::define() 19 0.004055 0.007886 0.01162 0.07% Other 2520 0.08471 0.1051 0.1492 0.96% -------------------------------------------------------------------------------------------- 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 10 B 12 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 1561 KiB 1562 KiB 1562 KiB 2346 KiB 2346 KiB 2346 KiB StateData::define() 192 192 1036 KiB 1037 KiB 1037 KiB 1173 KiB 1173 KiB 1173 KiB Castro::construct_ctu_hydro_source() 6400 6400 403 KiB 409 KiB 417 KiB 1164 KiB 1164 KiB 1164 KiB Castro::initMFs() 256 256 966 KiB 966 KiB 966 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 667 KiB 668 KiB 668 KiB 812 KiB 812 KiB 812 KiB Amr::writePlotFile() 64 64 3989 B 4079 B 4610 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 46 KiB 47 KiB 49 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 405 KiB 405 KiB 406 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 2225 B 4240 B 5827 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1240 B 1256 B 1290 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 16 KiB 16 KiB 17 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 142 KiB 142 KiB 142 KiB 175 KiB 175 KiB 175 KiB Castro::buildMetrics() 192 192 125 KiB 125 KiB 125 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 6102 9361 B 9906 B 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 233 B 236 B 244 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 392 B 396 B 401 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 17 B 19 B 34 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 32 KiB 32 KiB 32 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 47 B 47 B 47 B 52 KiB 52 KiB 52 KiB Gravity::get_new_grav_vector() 576 576 127 B 128 B 131 B 50 KiB 50 KiB 50 KiB Gravity::get_old_grav_vector() 640 640 175 B 223 B 264 B 50 KiB 50 KiB 50 KiB Gravity::update_max_rhs() 192 192 4 B 7 B 8 B 49 KiB 49 KiB 49 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 218 B 223 B 228 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 988 B 1009 B 1129 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 168 B 168 B 168 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1554 B 1617 B 1645 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 422 B 472 B 655 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 439 B 440 B 442 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 27 B 36 B 49 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 17 B 28 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 15 B 24 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