Initializing AMReX (26.01-17-g909028cb95b3)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.01-17-g909028cb95b3) initialized Starting run at 07:48:06 UTC on 2026-01-16. Successfully read inputs file ... Castro git describe: 26.01-10-g409d8e9cc AMReX git describe: 26.01-17-g909028cb9 Microphysics git describe: 26.01-12-g9148d8d1 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.184062405 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.002840801 new finest: 0 Gravity::fill_multipole_BCs() time = 0.034439329 Gravity::multilevel_solve_for_new_phi() time = 0.360563277 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.036824894 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.112736597 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.00753213 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 = -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 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 Castro::construct_old_gravity() time = 0.011342624 on level 0 Castro::do_old_sources() time = 0.028681561 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.157523671 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.013711232 on level 0 Castro::do_old_sources() time = 0.028588733 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.16457727 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.045040307 on level 0 Castro::do_old_sources() time = 0.048516921 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.151851694 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.009377406 on level 0 Castro::do_old_sources() time = 0.026765934 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.144421743 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.816833926 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.048250527 Gravity::solve_for_phi() time = 0.294934123 on level 0 Castro::construct_new_gravity() time = 0.302499775 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.147755658 on level 0 Castro::construct_old_gravity() time = 0.009685915 on level 0 Castro::do_old_sources() time = 0.026792512 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.145730292 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.806544264 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.03630206 Gravity::solve_for_phi() time = 0.26933749 on level 0 Castro::construct_new_gravity() time = 0.2768381 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.153977533 on level 0 Castro::construct_old_gravity() time = 0.009817191 on level 0 Castro::do_old_sources() time = 0.026846895 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.147978245 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.844082204 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.039195132 Gravity::solve_for_phi() time = 0.274611902 on level 0 Castro::construct_new_gravity() time = 0.282156311 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.151731235 on level 0 Castro::construct_old_gravity() time = 0.010059044 on level 0 Castro::do_old_sources() time = 0.02686263 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.147953038 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.836782846 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.040192009 Gravity::solve_for_phi() time = 0.27557868 on level 0 Castro::construct_new_gravity() time = 0.283083514 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.151741098 on level 0 Castro::construct_old_gravity() time = 0.009197699 on level 0 Castro::do_old_sources() time = 0.026383232 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.144029775 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.827502727 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.03683969 Gravity::solve_for_phi() time = 0.26955753 on level 0 Castro::construct_new_gravity() time = 0.277201944 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.147331588 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.007670109 on level 0 Castro::do_old_sources() time = 0.026420233 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.144007235 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.819318519 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.039544203 Gravity::solve_for_phi() time = 0.303489782 on level 0 Castro::construct_new_gravity() time = 0.311150621 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.147136121 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.008021639 on level 0 Castro::do_old_sources() time = 0.026450422 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.143940519 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.842993427 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.039667986 Gravity::solve_for_phi() time = 0.27483239 on level 0 Castro::construct_new_gravity() time = 0.282399695 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.150977401 on level 0 Castro::construct_old_gravity() time = 0.00974924 on level 0 Castro::do_old_sources() time = 0.027363843 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.147341181 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.828833838 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.040027308 Gravity::solve_for_phi() time = 0.2751604 on level 0 Castro::construct_new_gravity() time = 0.282706341 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.151123765 on level 0 Castro::construct_old_gravity() time = 0.009841126 on level 0 Castro::do_old_sources() time = 0.026856695 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.148270858 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.848198321 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.040009553 Gravity::solve_for_phi() time = 0.275186966 on level 0 Castro::construct_new_gravity() time = 0.28277174 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.151012313 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.009798705 on level 0 Castro::do_old_sources() time = 0.026788286 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.147338772 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.828573721 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.039437323 Gravity::solve_for_phi() time = 0.274226075 on level 0 Castro::construct_new_gravity() time = 0.281984943 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.150957439 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.010157392 on level 0 Castro::do_old_sources() time = 0.026904559 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.148287504 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.846879658 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.040148433 Gravity::solve_for_phi() time = 0.274622422 on level 0 Castro::construct_new_gravity() time = 0.282230627 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.150912098 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.00963415 on level 0 Castro::do_old_sources() time = 0.026760353 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.147363417 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.827241898 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.04042275 Gravity::solve_for_phi() time = 0.275258773 on level 0 Castro::construct_new_gravity() time = 0.282724571 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.151875756 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.009780225 on level 0 Castro::do_old_sources() time = 0.026894223 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.14718834 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.84621735 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.039777104 Gravity::solve_for_phi() time = 0.276402276 on level 0 Castro::construct_new_gravity() time = 0.28387064 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.150964377 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.063603692 on level 0 Castro::do_old_sources() time = 0.028076164 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.147224553 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.838547624 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.039992408 Gravity::solve_for_phi() time = 0.274870269 on level 0 Castro::construct_new_gravity() time = 0.282382289 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.150752214 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.00964026 on level 0 Castro::do_old_sources() time = 0.026860401 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.147369346 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.848297226 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.039511021 Gravity::solve_for_phi() time = 0.275672205 on level 0 Castro::construct_new_gravity() time = 0.283205925 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.150728153 on level 0 Castro::construct_old_gravity() time = 0.009839908 on level 0 Castro::do_old_sources() time = 0.026792264 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.147397881 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.836140832 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.045519627 Gravity::solve_for_phi() time = 0.29275447 on level 0 Castro::construct_new_gravity() time = 0.300320424 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.156493367 on level 0 Castro::construct_old_gravity() time = 0.015759205 on level 0 Castro::do_old_sources() time = 0.028730498 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.147421503 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.867406012 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.001016829165 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.039836866 Gravity::solve_for_phi() time = 0.274968604 on level 0 Castro::construct_new_gravity() time = 0.282594469 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.150745075 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.137357411e+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.469316165e+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.460231145e-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: 32.38259524 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.114216087 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.008358251 seconds Ending run at 07:48:41 UTC on 2026-01-16. Run time = 35.71823593 Run time without initialization = 32.50621124 Average number of zones advanced per microsecond: 0.001 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 35.7098 ... 35.7129 ... 35.7184 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 13.8427 13.9967 14.1814 39.70% amrex::Copy() 766 3.27 3.306 3.336 9.34% MultiFab::AddProduct() 3598 2.22 2.243 2.264 6.34% MLPoisson::Fsmooth() 6144 1.523 1.647 1.987 5.56% FabArray::ParallelCopy_finish() 1469 0.4467 1.291 1.707 4.78% Amr::InitAmr() 1 1.611 1.615 1.62 4.54% FabArray::setVal() 5232 1.366 1.386 1.409 3.95% MultiFab::Subtract() 159 1.363 1.378 1.393 3.90% Castro::computeTemp() 80 0.9819 1.071 1.211 3.39% FabArray::ParallelCopy_nowait() 1469 0.7525 0.7996 0.87 2.44% FabArray::sum() 3598 0.7516 0.76 0.7668 2.15% Gravity::fill_multipole_BCs() 19 0.658 0.7078 0.7453 2.09% MLPoisson::Fapply() 986 0.5573 0.5846 0.6591 1.85% check_for_negative_density() 20 0.166 0.4129 0.6113 1.71% Castro::do_advance_ctu() 20 0.5173 0.5242 0.5331 1.49% Castro::initMFs() 1 0.3814 0.3894 0.4002 1.12% Castro::reset_internal_energy(Fab) 160 0.3703 0.3751 0.3796 1.06% Castro::construct_old_gravity_source() 40 0.3265 0.3305 0.3337 0.93% FabArray::Saxpy() 111 0.3206 0.3235 0.3266 0.91% FillBoundary_finish() 6380 0.236 0.2806 0.3229 0.90% MLCellLinOp::applyBC() 7148 0.1283 0.1767 0.2994 0.84% Castro::construct_new_gravity_source() 17 0.2642 0.2668 0.2693 0.75% FillBoundary_nowait() 6380 0.1406 0.1859 0.2407 0.67% FabArray::Xpay() 915 0.192 0.2007 0.2221 0.62% Castro::initData() 1 0.0102 0.03537 0.1751 0.49% MultiFab::contains_nan() 37 0.1658 0.1681 0.172 0.48% TagBoxArray::collate() 2 0.0005643 0.03413 0.1563 0.44% Castro::normalize_species() 79 0.1391 0.1409 0.1423 0.40% MLMG::addInterpCorrection() 768 0.08305 0.08976 0.1077 0.30% StateData::FillBoundary(geom) 193 0 0.05098 0.1017 0.28% amrex::average_down 768 0.05939 0.06451 0.07938 0.22% Castro::construct_old_source() 240 0.07057 0.0714 0.07276 0.20% AmrLevel::derive() 18 0.05942 0.06058 0.06275 0.18% Castro::construct_new_source() 102 0.05855 0.05924 0.06003 0.17% amrex::Add() 128 0.05293 0.05366 0.05495 0.15% VisMF::Write(FabArray) 4 0.03788 0.03935 0.05082 0.14% FabArray::setDomainBndry() 187 0.0006543 0.02525 0.0498 0.14% FabArray::norminf() 244 0.03602 0.03796 0.04269 0.12% Castro::enforce_min_density() 79 0.0334 0.03388 0.03454 0.10% Gravity::get_old_grav_vector() 20 0.03002 0.03059 0.0315 0.09% MLMG::ResNormInf() 146 0.01741 0.02371 0.03048 0.09% Castro::volWgtSum() 46 0.02867 0.02913 0.02988 0.08% Gravity::get_new_grav_vector() 18 0.02735 0.02768 0.02857 0.08% Castro::reset_internal_energy(MultiFab) 80 0.02469 0.02522 0.02675 0.07% Castro::wd_update() 2 0.01253 0.01943 0.02472 0.07% Amr::writePlotFile() 2 0.007483 0.01936 0.02299 0.06% MLCellLinOp::defineAuxData() 19 0.01699 0.0182 0.02233 0.06% FabArray::FillBoundary() 6380 0.01702 0.01764 0.02103 0.06% MLMG::solve() 18 0.005626 0.01101 0.0189 0.05% MLCGSolver::bicgstab 8 0 0.001073 0.01717 0.05% Castro::evaluate_source_change() 196 0.01514 0.01553 0.01709 0.05% FillPatchIterator::Initialize 187 0.01556 0.01587 0.01646 0.05% MLMG::prepareForSolve() 18 0.01285 0.01333 0.01445 0.04% MLCellLinOp::smooth() 1544 0.0108 0.0113 0.0125 0.04% Gravity::actual_solve_with_mlmg() 19 0.00742 0.007896 0.009506 0.03% Castro::sum_integrated_quantities() 2 0.004893 0.005783 0.00668 0.02% MLCellLinOp::apply() 986 0.002479 0.002759 0.005839 0.02% Castro::initialize_do_advance() 20 0.00469 0.004953 0.005808 0.02% FabArray::ParallelCopy() 1469 0.004219 0.004388 0.004722 0.01% AmrLevel::FillPatch() 187 0.003205 0.003358 0.003921 0.01% MLCellLinOp::correctionResidual() 768 0.002353 0.002409 0.002539 0.01% MLMG::mgVcycle() 128 0.002194 0.002283 0.00246 0.01% MLLinOp::defineGrids() 19 0.001935 0.002088 0.002373 0.01% Castro::save_data_for_retry() 3 0.001485 0.001727 0.002153 0.01% MLMG:computeResOfCorrection() 768 0.001865 0.001908 0.002004 0.01% AmrMesh::MakeNewGrids() 2 0.001702 0.001749 0.001811 0.01% StateDataPhysBCFunct::() 187 0.001022 0.001222 0.001384 0.00% MLMG::actualBottomSolve() 128 0.0001858 0.0002596 0.001285 0.00% Castro::subcycle_advance_ctu() 1 0.0002465 0.0003472 0.001036 0.00% main() 1 0.0007416 0.0008524 0.001011 0.00% FillPatchIterator::FillFromLevel0() 187 0.0008885 0.0009251 0.0009965 0.00% Gravity::solve_for_phi() 17 0.0006735 0.0007112 0.0007703 0.00% FillPatchSingleLevel 187 0.0006775 0.0006937 0.0007345 0.00% MLMG::mgVcycle_down::0 128 0.0006399 0.0006658 0.0007127 0.00% MLMG::mgVcycle_down::1 128 0.0006206 0.0006365 0.0006936 0.00% MLMG::mgVcycle_down::2 128 0.0005974 0.0006362 0.000665 0.00% MLMG::mgVcycle_down::4 128 0.0005534 0.0005739 0.0006578 0.00% MLMG::mgVcycle_down::3 128 0.0005635 0.0005946 0.0006561 0.00% MLMG::mgVcycle_down::5 128 0.0005848 0.0006057 0.0006494 0.00% MLCellLinOp::solutionResidual() 147 0.0005341 0.0005534 0.0005846 0.00% Castro::do_old_sources() 40 0.0004694 0.0004926 0.0005735 0.00% Castro::clean_state() 79 0.0004923 0.0005116 0.0005386 0.00% MLMG::oneIter() 128 0.0004434 0.0004572 0.0005041 0.00% MLMG::mgVcycle_up::0 128 0.0004212 0.0004326 0.0004854 0.00% MLMG::mgVcycle_up::1 128 0.0004164 0.0004352 0.0004803 0.00% Castro::Castro() 1 0.0003855 0.0004139 0.0004644 0.00% MLMG::mgVcycle_up::2 128 0.0003892 0.0004123 0.0004421 0.00% MLMG::mgVcycle_up::3 128 0.0003789 0.000395 0.0004354 0.00% MLMG::mgVcycle_up::4 128 0.0003831 0.0003948 0.0004296 0.00% MLMG::mgVcycle_up::5 128 0.0003705 0.0003809 0.000397 0.00% MLMG::computeResidual() 128 0.0003046 0.0003134 0.000333 0.00% TagBoxArray::mapPRD 2 0.0002712 0.0002924 0.0003224 0.00% MLMG::mgVcycle_bottom 128 0.0002762 0.0002869 0.000309 0.00% Castro::construct_new_gravity() 17 0.0001762 0.0001954 0.0003059 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002218 0.0002407 0.0002681 0.00% Castro::do_new_sources() 17 0.0002212 0.0002332 0.0002638 0.00% Castro::problem_post_timestep() 2 1.315e-05 2.97e-05 0.000254 0.00% Gravity::test_composite_phi() 1 0.000109 0.0001204 0.0001958 0.00% Amr::coarseTimeStep() 1 0.0001439 0.0001659 0.0001722 0.00% Castro::construct_old_gravity() 20 0.0001189 0.000136 0.0001695 0.00% Castro::check_for_nan() 37 0.0001248 0.000135 0.0001658 0.00% MLPoisson::define() 19 0.0001326 0.0001459 0.000161 0.00% Castro::apply_source_to_state() 37 0.0001212 0.000129 0.0001511 0.00% Castro::initialize_advance() 1 0.0001092 0.000121 0.0001357 0.00% Castro::retry_advance_ctu() 11 7.553e-05 8.906e-05 0.0001318 0.00% Amr::defBaseLevel() 1 9.653e-05 0.0001029 0.0001133 0.00% MLLinOp::define() 19 8.732e-05 9.158e-05 0.0001029 0.00% Castro::create_source_corrector() 17 5.477e-05 6.251e-05 7.691e-05 0.00% Amr::grid_places() 2 5.17e-05 5.981e-05 7.542e-05 0.00% Castro::expand_state() 20 5.704e-05 6.081e-05 7.473e-05 0.00% MLMG::computeMLResidual() 18 5.268e-05 5.68e-05 6.422e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.601e-05 5.932e-05 6.292e-05 0.00% Castro::derive() 18 4.764e-05 5.212e-05 6.229e-05 0.00% Amr::FinalizeInit() 1 9.025e-06 1.269e-05 5.346e-05 0.00% Amr::regrid() 1 1.567e-05 2.341e-05 4.606e-05 0.00% Gravity::actual_multilevel_solve() 1 3.078e-05 3.184e-05 3.344e-05 0.00% Amr::timeStep() 1 1.478e-05 1.768e-05 3.336e-05 0.00% Castro::post_init() 1 1.762e-05 1.972e-05 2.59e-05 0.00% Amr::InitializeInit() 1 8.4e-06 9.896e-06 1.726e-05 0.00% Castro::post_regrid() 2 1.084e-05 1.164e-05 1.288e-05 0.00% Castro::advance() 1 7.949e-06 9.581e-06 1.257e-05 0.00% Amr::bldFineLevels() 1 7.914e-06 8.442e-06 9.293e-06 0.00% Castro::post_timestep() 1 6.187e-06 6.855e-06 8.858e-06 0.00% Amr::init() 1 6.385e-06 7.336e-06 8.417e-06 0.00% Amr::initialInit() 1 3.353e-06 3.631e-06 4.174e-06 0.00% Other 10909 0.1256 0.1585 0.2245 0.63% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 35.71 35.71 35.72 100.00% Amr::coarseTimeStep() 1 32.51 32.51 32.51 91.01% Amr::timeStep() 1 32.38 32.38 32.38 90.66% Castro::advance() 1 32.22 32.22 32.23 90.22% Castro::subcycle_advance_ctu() 1 32.09 32.1 32.1 89.87% Castro::do_advance_ctu() 20 31.94 31.95 31.96 89.47% Castro::construct_ctu_hydro_source() 20 16.25 16.28 16.32 45.70% Castro::clean_state() 79 8.395 8.541 8.747 24.49% Castro::computeTemp() 80 5.199 5.317 5.499 15.39% Gravity::solve_phi_with_mlmg() 19 5.03 5.092 5.125 14.35% Castro::construct_new_gravity() 17 4.766 4.826 4.859 13.60% Gravity::solve_for_phi() 17 4.635 4.696 4.73 13.24% Gravity::actual_solve_with_mlmg() 19 4.36 4.372 4.384 12.27% Castro::reset_internal_energy(MultiFab) 80 4.199 4.245 4.29 12.01% MLMG::solve() 18 4.274 4.28 4.282 11.99% Castro::evaluate_source_change() 196 3.733 3.774 3.808 10.66% MLMG::oneIter() 128 3.735 3.74 3.746 10.49% MLMG::mgVcycle() 128 3.682 3.686 3.691 10.33% Castro::do_old_sources() 40 3.391 3.441 3.489 9.77% Castro::initialize_do_advance() 20 3.254 3.328 3.388 9.48% amrex::Copy() 766 3.27 3.306 3.336 9.34% Castro::enforce_min_density() 79 3.113 3.147 3.173 8.88% MLCellLinOp::smooth() 1544 2.039 2.197 2.652 7.43% FabArray::ParallelCopy() 1469 1.298 2.156 2.57 7.20% Castro::do_new_sources() 17 2.444 2.5 2.555 7.15% AmrLevel::FillPatch() 187 2.073 2.171 2.278 6.38% MultiFab::AddProduct() 3598 2.22 2.243 2.264 6.34% MLPoisson::Fsmooth() 6144 1.523 1.647 1.987 5.56% FabArray::ParallelCopy_finish() 1469 0.4467 1.291 1.707 4.78% Amr::InitAmr() 1 1.611 1.615 1.62 4.54% Amr::init() 1 1.592 1.592 1.592 4.46% FillPatchIterator::Initialize 187 1.255 1.363 1.479 4.14% Amr::initialInit() 1 1.471 1.471 1.471 4.12% FillPatchIterator::FillFromLevel0() 187 1.189 1.322 1.446 4.05% FillPatchSingleLevel 187 1.188 1.321 1.445 4.04% FabArray::setVal() 5232 1.366 1.386 1.409 3.95% MultiFab::Subtract() 159 1.363 1.378 1.393 3.90% MLMG::mgVcycle_up::1 128 0.2806 0.7574 1.235 3.46% MLMG::mgVcycle_down::0 128 1.133 1.139 1.143 3.20% MLMG::addInterpCorrection() 768 0.1459 0.8437 1.057 2.96% FabArray::ParallelCopy_nowait() 1469 0.8136 0.8608 0.9325 2.61% Amr::FinalizeInit() 1 0.6958 0.8409 0.8771 2.46% MLCellLinOp::apply() 986 0.6866 0.721 0.8223 2.30% MLCellLinOp::applyBC() 7148 0.5998 0.676 0.8143 2.28% MLMG::mgVcycle_up::0 128 0.8018 0.8056 0.8093 2.27% Castro::expand_state() 20 0.6815 0.7366 0.7919 2.22% Amr::InitializeInit() 1 0.5941 0.6302 0.7754 2.17% Amr::defBaseLevel() 1 0.5941 0.6302 0.7753 2.17% MLMG::mgVcycle_up::2 128 0.01543 0.317 0.7742 2.17% FabArray::sum() 3598 0.7516 0.76 0.7668 2.15% Gravity::fill_multipole_BCs() 19 0.666 0.7158 0.7532 2.11% MLPoisson::Fapply() 986 0.5573 0.5846 0.6591 1.85% check_for_negative_density() 20 0.166 0.4129 0.6113 1.71% MLMG:computeResOfCorrection() 768 0.464 0.4998 0.5964 1.67% MLCellLinOp::correctionResidual() 768 0.4621 0.4979 0.5944 1.66% Castro::construct_old_source() 240 0.5701 0.5765 0.5817 1.63% FabArray::FillBoundary() 6380 0.4237 0.4993 0.5663 1.59% Castro::post_init() 1 0.4708 0.4805 0.4843 1.36% Castro::construct_new_source() 102 0.4697 0.4742 0.4783 1.34% MLMG::mgVcycle_up::3 128 0.01535 0.05075 0.4554 1.27% MLCellLinOp::solutionResidual() 147 0.4084 0.424 0.4336 1.21% Castro::Castro() 1 0.4 0.4082 0.4186 1.17% Castro::initMFs() 1 0.382 0.39 0.4008 1.12% MLMG::mgVcycle_down::1 128 0.3699 0.3787 0.3902 1.09% Castro::reset_internal_energy(Fab) 160 0.3703 0.3751 0.3796 1.06% MLMG::computeResidual() 128 0.3534 0.3578 0.3635 1.02% Castro::initData() 1 0.1921 0.2201 0.3619 1.01% Gravity::multilevel_solve_for_new_phi() 1 0.3605 0.3605 0.3606 1.01% Gravity::actual_multilevel_solve() 1 0.3604 0.3604 0.3605 1.01% Castro::construct_old_gravity_source() 40 0.3265 0.3305 0.3337 0.93% FabArray::Saxpy() 111 0.3206 0.3235 0.3266 0.91% FillBoundary_finish() 6380 0.236 0.2806 0.3229 0.90% Amr::bldFineLevels() 1 0.1115 0.2559 0.2902 0.81% Castro::construct_new_gravity_source() 17 0.2642 0.2668 0.2693 0.75% FillBoundary_nowait() 6380 0.1551 0.2011 0.2565 0.72% Castro::construct_old_gravity() 20 0.1768 0.2105 0.2502 0.70% Gravity::get_old_grav_vector() 20 0.1767 0.2104 0.2501 0.70% Castro::create_source_corrector() 17 0.1704 0.216 0.2466 0.69% Amr::writePlotFile() 2 0.2274 0.2274 0.2275 0.64% FabArray::Xpay() 915 0.192 0.2007 0.2221 0.62% Castro::post_regrid() 2 0.2059 0.2096 0.2194 0.61% MLMG::mgVcycle_down::2 128 0.018 0.1091 0.21 0.59% MLMG::mgVcycle_down::3 128 0.01739 0.03911 0.1892 0.53% Amr::grid_places() 2 0.005596 0.1507 0.1869 0.52% AmrMesh::MakeNewGrids() 2 0.005539 0.1507 0.1868 0.52% TagBoxArray::mapPRD 2 0.0017 0.1133 0.1829 0.51% Castro::derive() 18 0.1677 0.1703 0.1734 0.49% AmrLevel::derive() 18 0.1676 0.1703 0.1734 0.49% Castro::check_for_nan() 37 0.166 0.1682 0.1721 0.48% MultiFab::contains_nan() 37 0.1658 0.1681 0.172 0.48% Castro::post_timestep() 1 0.1563 0.1602 0.1638 0.46% TagBoxArray::collate() 2 0.0005643 0.03413 0.1563 0.44% Castro::retry_advance_ctu() 11 0.1358 0.138 0.1483 0.42% Castro::normalize_species() 79 0.1391 0.1409 0.1423 0.40% Gravity::get_new_grav_vector() 18 0.1277 0.1289 0.1302 0.36% Castro::initialize_advance() 1 0.1225 0.1249 0.1293 0.36% amrex::average_down 768 0.0719 0.08653 0.121 0.34% Amr::regrid() 1 0.1061 0.108 0.1142 0.32% Castro::apply_source_to_state() 37 0.1071 0.108 0.109 0.31% MLMG::mgVcycle_bottom 128 0.0004631 0.007149 0.1072 0.30% MLMG::actualBottomSolve() 128 0.0001858 0.006862 0.1069 0.30% MLMG::mgVcycle_down::4 128 0.01739 0.02328 0.1063 0.30% MLCGSolver::bicgstab 8 0 0.006539 0.1046 0.29% StateDataPhysBCFunct::() 187 0.001022 0.05221 0.103 0.29% StateData::FillBoundary(geom) 193 0 0.05098 0.1017 0.28% MLMG::mgVcycle_up::4 128 0.01541 0.01963 0.07823 0.22% Castro::sum_integrated_quantities() 2 0.07678 0.07693 0.07695 0.22% Castro::save_data_for_retry() 3 0.0707 0.07213 0.07576 0.21% MLMG::computeMLResidual() 18 0.04914 0.06368 0.07323 0.21% MLPoisson::define() 19 0.03757 0.05142 0.06615 0.19% MLMG::ResNormInf() 146 0.04984 0.05587 0.0624 0.17% amrex::Add() 128 0.05293 0.05366 0.05495 0.15% MLMG::prepareForSolve() 18 0.04452 0.04754 0.05422 0.15% MLMG::mgVcycle_down::5 128 0.01738 0.01982 0.051 0.14% VisMF::Write(FabArray) 4 0.03788 0.03935 0.05082 0.14% FabArray::setDomainBndry() 187 0.0006543 0.02525 0.0498 0.14% Gravity::test_composite_phi() 1 0.04593 0.046 0.04608 0.13% FabArray::norminf() 244 0.03602 0.03796 0.04269 0.12% MLMG::mgVcycle_up::5 128 0.01551 0.01736 0.04076 0.11% Castro::problem_post_timestep() 2 0.02819 0.03319 0.03773 0.11% Castro::wd_update() 2 0.02818 0.03316 0.03772 0.11% MLCellLinOp::defineAuxData() 19 0.02203 0.02411 0.03016 0.08% Castro::volWgtSum() 46 0.02867 0.02913 0.02988 0.08% MLLinOp::define() 19 0.006007 0.01782 0.02973 0.08% MLLinOp::defineGrids() 19 0.005917 0.01773 0.02964 0.08% Other 10909 0.2067 0.2514 0.3339 0.93% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Async Memory Usage: -------------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max -------------------------------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 6400 426 KiB 430 KiB 436 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 1930 B 1959 B 1978 B 89 KiB 89 KiB 89 KiB -------------------------------------------------------------------------------------------------------------------- Pinned Memory Usage: -------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max -------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 0 B 0 B 1 B 8192 KiB 8192 KiB 8192 KiB -------------------------------------------------------------------------------------------------------------- Cpu Memory Usage: ----------------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max ----------------------------------------------------------------------------------------------------------------------- Castro::save_data_for_retry() 416 1720 KiB 1720 KiB 1720 KiB 2329 KiB 2329 KiB 2329 KiB Castro::construct_ctu_hydro_source() 6400 426 KiB 430 KiB 436 KiB 1176 KiB 1176 KiB 1176 KiB StateData::define() 192 1111 KiB 1111 KiB 1111 KiB 1164 KiB 1164 KiB 1164 KiB Castro::initMFs() 256 1089 KiB 1090 KiB 1090 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 730 KiB 731 KiB 731 KiB 803 KiB 803 KiB 803 KiB Amr::writePlotFile() 64 4434 B 4513 B 4787 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 68 KiB 69 KiB 70 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 18 KiB 18 KiB 19 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 471 KiB 471 KiB 471 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 3081 B 5013 B 6412 B 225 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1236 B 1257 B 1284 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 29 KiB 29 KiB 29 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 142 KiB 142 KiB 142 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 135 KiB 135 KiB 135 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 12 KiB 12 KiB 15 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 1280 774 B 783 B 792 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 656 B 662 B 668 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 1930 B 1959 B 1978 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 19 B 20 B 30 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 35 KiB 35 KiB 35 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 67 B 67 B 68 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 4 B 8 B 10 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 164 B 165 B 167 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 227 B 270 B 322 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 34 KiB 34 KiB 34 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 311 B 313 B 318 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 1710 B 1728 B 1744 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 165 B 165 B 165 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 2125 B 2153 B 2168 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 1 B 1 B 1 B 16 KiB 16 KiB 16 KiB MLCellLinOp::defineAuxData() 11704 560 B 625 B 869 B 4608 B 5151 B 7152 B BndryData::define() 4864 588 B 588 B 589 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 20 B 30 B 38 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 0 B 8 B 13 B 2916 B 2916 B 2916 B OwnerMask() 64 0 B 0 B 0 B 2592 B 2592 B 2592 B MLMG::addInterpCorrection() 5632 1 B 18 B 28 B 1024 B 1024 B 1024 B amrex::average_down 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 0 B 2 B 3 B 648 B 648 B 648 B MLCGSolver::bicgstab 640 0 B 0 B 1 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 0 B 0 B 0 B 8 B 8 B 8 B ----------------------------------------------------------------------------------------------------------------------- AMReX (26.01-17-g909028cb95b3) finalized