Initializing AMReX (26.01-22-gf9456be2b517)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.01-22-gf9456be2b517) initialized Starting run at 08:17:23 UTC on 2026-01-21. Successfully read inputs file ... Castro git describe: 26.01-11-g3d7eded2f AMReX git describe: 26.01-22-gf9456be2b Microphysics git describe: 26.01-18-gd7041c94 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.167670985 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.002819826 new finest: 0 Gravity::fill_multipole_BCs() time = 0.033819299 Gravity::multilevel_solve_for_new_phi() time = 0.329185783 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.034135588 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.111626309 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.007055653 seconds [Level 0 step 1] ADVANCE at time 0 with dt = 0.25  Beginning subcycle 1 starting at time 0 with dt = 0.25 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.01295059363 at index 21, 172, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035502 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477044 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.769409 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071148 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.7670713549 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.7670713549 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071148 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.769409 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477044 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035502 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -0.01295059363 at index 21, 83, 0 Castro::construct_old_gravity() time = 0.010551838 on level 0 Castro::do_old_sources() time = 0.028763755 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.157207776 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.018940299 on level 0 Castro::do_old_sources() time = 0.032253203 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.1574895 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.017099816 on level 0 Castro::do_old_sources() time = 0.033784817 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.154056842 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.014527916 on level 0 Castro::do_old_sources() time = 0.029070388 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.144437955 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.810405233 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.036536146 Gravity::solve_for_phi() time = 0.269328849 on level 0 Castro::construct_new_gravity() time = 0.276932651 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.148265271 on level 0 Castro::construct_old_gravity() time = 0.008095922 on level 0 Castro::do_old_sources() time = 0.026314751 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.144659815 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.81059553 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.036533353 Gravity::solve_for_phi() time = 0.26784139 on level 0 Castro::construct_new_gravity() time = 0.275455961 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.148232749 on level 0 Castro::construct_old_gravity() time = 0.009534911 on level 0 Castro::do_old_sources() time = 0.02637287 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.144630549 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.824605421 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.037071509 Gravity::solve_for_phi() time = 0.26941329 on level 0 Castro::construct_new_gravity() time = 0.277081604 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.148243055 on level 0 Castro::construct_old_gravity() time = 0.014618469 on level 0 Castro::do_old_sources() time = 0.03035645 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.144266574 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.817113383 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.037360099 Gravity::solve_for_phi() time = 0.270403842 on level 0 Castro::construct_new_gravity() time = 0.277908816 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.148064334 on level 0 Castro::construct_old_gravity() time = 0.009462032 on level 0 Castro::do_old_sources() time = 0.026546956 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.144527878 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.834296248 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.036908082 Gravity::solve_for_phi() time = 0.267651926 on level 0 Castro::construct_new_gravity() time = 0.275248959 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.147814472 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.008698924 on level 0 Castro::do_old_sources() time = 0.026395392 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.144677898 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.816365345 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.037167187 Gravity::solve_for_phi() time = 0.269105922 on level 0 Castro::construct_new_gravity() time = 0.276647385 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.14778814 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.008754104 on level 0 Castro::do_old_sources() time = 0.02655338 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.144608479 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.838224384 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.03741883 Gravity::solve_for_phi() time = 0.270613577 on level 0 Castro::construct_new_gravity() time = 0.278232607 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.148314356 on level 0 Castro::construct_old_gravity() time = 0.009527336 on level 0 Castro::do_old_sources() time = 0.026468236 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.144211363 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.811092494 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.037610804 Gravity::solve_for_phi() time = 0.274789907 on level 0 Castro::construct_new_gravity() time = 0.282538955 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.148463578 on level 0 Castro::construct_old_gravity() time = 0.009457419 on level 0 Castro::do_old_sources() time = 0.026472174 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.144655197 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.837158254 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.037353944 Gravity::solve_for_phi() time = 0.270169323 on level 0 Castro::construct_new_gravity() time = 0.277761114 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.148284003 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.009469955 on level 0 Castro::do_old_sources() time = 0.026291812 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.144547939 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.817344373 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.037208509 Gravity::solve_for_phi() time = 0.269277242 on level 0 Castro::construct_new_gravity() time = 0.276813537 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.14794651 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.009426326 on level 0 Castro::do_old_sources() time = 0.026460352 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.145062109 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.831119394 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.037535581 Gravity::solve_for_phi() time = 0.268775918 on level 0 Castro::construct_new_gravity() time = 0.276430151 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.148148904 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.00912952 on level 0 Castro::do_old_sources() time = 0.026748633 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.144322348 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.830680411 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.045498886 Gravity::solve_for_phi() time = 0.288508731 on level 0 Castro::construct_new_gravity() time = 0.296111719 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.148533163 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.009433335 on level 0 Castro::do_old_sources() time = 0.026688178 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.144201264 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.85522849 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.037181939 Gravity::solve_for_phi() time = 0.286242159 on level 0 Castro::construct_new_gravity() time = 0.293757096 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.1493149 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.009740901 on level 0 Castro::do_old_sources() time = 0.026346262 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.146605289 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.82913226 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.038165501 Gravity::solve_for_phi() time = 0.278973772 on level 0 Castro::construct_new_gravity() time = 0.286492155 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.14858158 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.010823069 on level 0 Castro::do_old_sources() time = 0.028391048 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.146082562 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.866379061 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.038133807 Gravity::solve_for_phi() time = 0.276211358 on level 0 Castro::construct_new_gravity() time = 0.283729406 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.156796091 on level 0 Castro::construct_old_gravity() time = 0.010793916 on level 0 Castro::do_old_sources() time = 0.02693978 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.146555006 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.856464671 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.037675636 Gravity::solve_for_phi() time = 0.271821356 on level 0 Castro::construct_new_gravity() time = 0.279319303 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.150918358 on level 0 Castro::construct_old_gravity() time = 0.010159931 on level 0 Castro::do_old_sources() time = 0.026989794 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.146955439 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.843337899 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.001024863084 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.037547774 Gravity::solve_for_phi() time = 0.269125243 on level 0 Castro::construct_new_gravity() time = 0.276661909 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.14800705 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.12628373 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.110220165 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.006878125 seconds Ending run at 08:17:58 UTC on 2026-01-21. Run time = 35.12591769 Run time without initialization = 32.2443322 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.1180 ... 35.1205 ... 35.1260 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 13.8862 13.9514 14.1680 40.33% amrex::Copy() 766 3.272 3.279 3.317 9.44% MultiFab::AddProduct() 3598 2.23 2.235 2.262 6.44% MLPoisson::Fsmooth() 6144 1.513 1.626 1.95 5.55% FabArray::ParallelCopy_finish() 1469 0.4321 1.194 1.537 4.38% FabArray::setVal() 5232 1.369 1.377 1.392 3.96% MultiFab::Subtract() 159 1.358 1.362 1.373 3.91% Amr::InitAmr() 1 1.342 1.344 1.35 3.84% Castro::computeTemp() 80 0.9688 1.062 1.216 3.46% FabArray::ParallelCopy_nowait() 1469 0.7441 0.7942 0.8578 2.44% FabArray::sum() 3598 0.7476 0.7497 0.7586 2.16% Gravity::fill_multipole_BCs() 19 0.6325 0.6724 0.6991 1.99% MLPoisson::Fapply() 986 0.5596 0.5823 0.6534 1.86% check_for_negative_density() 20 0.226 0.482 0.5679 1.62% Castro::do_advance_ctu() 20 0.5197 0.522 0.5298 1.51% Castro::initMFs() 1 0.3882 0.3927 0.4239 1.21% Castro::reset_internal_energy(Fab) 160 0.3699 0.3722 0.3804 1.08% Castro::construct_old_gravity_source() 40 0.3235 0.3247 0.3304 0.94% FabArray::Saxpy() 111 0.3229 0.3245 0.3286 0.94% MLCellLinOp::applyBC() 7148 0.1241 0.1744 0.2957 0.84% FillBoundary_finish() 6380 0.2205 0.2488 0.2918 0.83% Castro::construct_new_gravity_source() 17 0.2631 0.2644 0.2688 0.77% FillBoundary_nowait() 6380 0.1352 0.1841 0.2388 0.68% FabArray::Xpay() 915 0.1917 0.199 0.2191 0.62% Castro::initData() 1 0.01035 0.03518 0.1713 0.49% MultiFab::contains_nan() 37 0.1649 0.1659 0.1681 0.48% Castro::normalize_species() 79 0.1394 0.1401 0.1428 0.41% TagBoxArray::collate() 2 0.0005163 0.02553 0.1264 0.36% MLMG::addInterpCorrection() 768 0.08335 0.0893 0.1066 0.30% StateData::FillBoundary(geom) 193 0 0.0508 0.1025 0.29% amrex::average_down 768 0.05924 0.06404 0.0776 0.22% Castro::construct_old_source() 240 0.07114 0.0716 0.07257 0.21% AmrLevel::derive() 18 0.05981 0.06012 0.06061 0.17% Castro::construct_new_source() 102 0.05914 0.05937 0.05992 0.17% amrex::Add() 128 0.05289 0.05312 0.05347 0.15% VisMF::Write(FabArray) 4 0.03677 0.03812 0.04993 0.14% FabArray::setDomainBndry() 187 0.0006546 0.02515 0.04983 0.14% FabArray::norminf() 244 0.03616 0.03663 0.04114 0.12% Castro::enforce_min_density() 79 0.03324 0.0335 0.03391 0.10% Gravity::get_old_grav_vector() 20 0.03023 0.03045 0.03077 0.09% Castro::volWgtSum() 46 0.02839 0.02854 0.02888 0.08% Gravity::get_new_grav_vector() 18 0.02754 0.02768 0.02814 0.08% Castro::reset_internal_energy(MultiFab) 80 0.0244 0.02478 0.02544 0.07% MLCellLinOp::defineAuxData() 19 0.01729 0.01842 0.02244 0.06% FabArray::FillBoundary() 6380 0.01692 0.01746 0.02066 0.06% Castro::wd_update() 2 0.01028 0.01514 0.01829 0.05% Amr::writePlotFile() 2 0.006692 0.01649 0.01819 0.05% Castro::evaluate_source_change() 196 0.01486 0.01548 0.01683 0.05% MLCGSolver::bicgstab 8 0 0.001049 0.01679 0.05% FillPatchIterator::Initialize 187 0.01527 0.01573 0.01624 0.05% MLMG::ResNormInf() 146 0.009252 0.01393 0.01604 0.05% MLMG::prepareForSolve() 18 0.01298 0.01335 0.01476 0.04% MLCellLinOp::smooth() 1544 0.01072 0.01124 0.0132 0.04% Gravity::actual_solve_with_mlmg() 19 0.007387 0.007862 0.009481 0.03% MLMG::solve() 18 0.001414 0.004555 0.008323 0.02% MLCellLinOp::apply() 986 0.002514 0.002777 0.005865 0.02% Castro::initialize_do_advance() 20 0.004673 0.004996 0.005497 0.02% FabArray::ParallelCopy() 1469 0.004172 0.004321 0.004789 0.01% Castro::sum_integrated_quantities() 2 0.002995 0.003774 0.004071 0.01% AmrLevel::FillPatch() 187 0.003205 0.003308 0.003642 0.01% MLCellLinOp::correctionResidual() 768 0.002284 0.002345 0.002464 0.01% MLMG::mgVcycle() 128 0.002168 0.002236 0.002338 0.01% MLMG:computeResOfCorrection() 768 0.001856 0.001904 0.00204 0.01% AmrMesh::MakeNewGrids() 2 0.001698 0.001744 0.001798 0.01% Castro::save_data_for_retry() 3 0.001368 0.001706 0.00178 0.01% StateDataPhysBCFunct::() 187 0.001013 0.001208 0.001342 0.00% MLMG::actualBottomSolve() 128 0.0001844 0.0002567 0.001268 0.00% Castro::subcycle_advance_ctu() 1 0.0002542 0.000345 0.001178 0.00% FillPatchIterator::FillFromLevel0() 187 0.0009041 0.0009334 0.001003 0.00% main() 1 0.0007006 0.0007889 0.0009216 0.00% Gravity::solve_for_phi() 17 0.0006726 0.0007152 0.0007898 0.00% MLMG::mgVcycle_down::2 128 0.00059 0.0006249 0.0007294 0.00% MLMG::mgVcycle_down::0 128 0.0006536 0.000669 0.000718 0.00% FillPatchSingleLevel 187 0.0006603 0.0006792 0.0007073 0.00% MLMG::mgVcycle_down::1 128 0.0006194 0.0006334 0.000691 0.00% MLMG::mgVcycle_down::3 128 0.0005607 0.0005893 0.0006569 0.00% MLMG::mgVcycle_down::5 128 0.0005832 0.0006001 0.0006481 0.00% MLMG::mgVcycle_down::4 128 0.0005531 0.0005734 0.0006466 0.00% Castro::Castro() 1 0.0004458 0.000516 0.0005996 0.00% MLCellLinOp::solutionResidual() 147 0.0005409 0.0005539 0.0005796 0.00% Castro::clean_state() 79 0.0004875 0.0005128 0.0005518 0.00% Castro::do_old_sources() 40 0.0004552 0.0004879 0.0005311 0.00% MLMG::oneIter() 128 0.0004477 0.0004595 0.0005006 0.00% MLMG::mgVcycle_up::1 128 0.0004134 0.0004323 0.0004876 0.00% MLMG::mgVcycle_up::0 128 0.0004171 0.0004286 0.0004702 0.00% MLMG::mgVcycle_up::2 128 0.0003859 0.0004065 0.0004605 0.00% MLMG::mgVcycle_up::3 128 0.0003732 0.0003876 0.000434 0.00% MLMG::mgVcycle_up::4 128 0.0003794 0.0003893 0.0004289 0.00% MLMG::mgVcycle_up::5 128 0.0003669 0.0003749 0.0003894 0.00% MLMG::computeResidual() 128 0.0002976 0.0003087 0.0003253 0.00% Castro::construct_new_gravity() 17 0.0001764 0.0001937 0.0003189 0.00% MLMG::mgVcycle_bottom 128 0.0002738 0.0002846 0.0003051 0.00% TagBoxArray::mapPRD 2 0.0002621 0.0002857 0.0003042 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002211 0.000238 0.0002779 0.00% Castro::problem_post_timestep() 2 1.295e-05 3.06e-05 0.0002739 0.00% Castro::do_new_sources() 17 0.0002235 0.0002324 0.0002524 0.00% Castro::retry_advance_ctu() 11 7.98e-05 9.618e-05 0.0002095 0.00% Amr::coarseTimeStep() 1 0.0001356 0.0001607 0.0001655 0.00% Castro::check_for_nan() 37 0.0001236 0.000133 0.000153 0.00% Castro::construct_old_gravity() 20 0.0001202 0.0001316 0.0001489 0.00% MLPoisson::define() 19 0.000126 0.0001339 0.0001457 0.00% Castro::initialize_advance() 1 0.0001135 0.000124 0.000143 0.00% Castro::apply_source_to_state() 37 0.0001216 0.000129 0.0001416 0.00% Gravity::test_composite_phi() 1 0.0001054 0.0001106 0.000139 0.00% Amr::defBaseLevel() 1 9.714e-05 0.000106 0.0001237 0.00% Amr::FinalizeInit() 1 1.248e-05 2.309e-05 0.000114 0.00% Amr::grid_places() 2 5.28e-05 6.649e-05 0.0001136 0.00% Castro::expand_state() 20 5.549e-05 5.978e-05 6.843e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.69e-05 5.919e-05 6.832e-05 0.00% MLMG::computeMLResidual() 18 5.425e-05 5.795e-05 6.747e-05 0.00% MLMG::getGradSolution() 18 5.361e-05 5.885e-05 6.531e-05 0.00% Castro::create_source_corrector() 17 5.438e-05 5.83e-05 6.236e-05 0.00% Amr::regrid() 1 1.66e-05 2.346e-05 5.678e-05 0.00% Castro::derive() 18 4.79e-05 4.986e-05 5.454e-05 0.00% Gravity::actual_multilevel_solve() 1 3.034e-05 3.166e-05 3.487e-05 0.00% Amr::timeStep() 1 1.469e-05 1.9e-05 3.293e-05 0.00% Amr::init() 1 6.237e-06 8.465e-06 2.53e-05 0.00% Castro::post_timestep() 1 6.198e-06 7.58e-06 2.237e-05 0.00% Castro::post_regrid() 2 1.115e-05 1.245e-05 2.049e-05 0.00% Castro::post_init() 1 1.607e-05 1.7e-05 1.865e-05 0.00% Amr::InitializeInit() 1 9.026e-06 1.032e-05 1.745e-05 0.00% Castro::advance() 1 8.16e-06 9.461e-06 1.359e-05 0.00% Amr::bldFineLevels() 1 7.593e-06 8.558e-06 1.299e-05 0.00% Amr::initialInit() 1 3.141e-06 3.557e-06 4.9e-06 0.00% Other 10929 0.1237 0.1467 0.2081 0.59% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 35.12 35.12 35.13 100.00% Amr::coarseTimeStep() 1 32.24 32.24 32.24 91.79% Amr::timeStep() 1 32.13 32.13 32.13 91.46% Castro::advance() 1 31.97 31.97 31.97 91.02% Castro::subcycle_advance_ctu() 1 31.84 31.84 31.85 90.67% Castro::do_advance_ctu() 20 31.7 31.7 31.71 90.27% Castro::construct_ctu_hydro_source() 20 16.25 16.29 16.34 46.51% Castro::clean_state() 79 8.369 8.472 8.661 24.66% Castro::computeTemp() 80 5.175 5.275 5.449 15.51% Gravity::solve_phi_with_mlmg() 19 4.922 4.972 4.993 14.22% Castro::construct_new_gravity() 17 4.69 4.74 4.762 13.56% Gravity::solve_for_phi() 17 4.561 4.611 4.633 13.19% Gravity::actual_solve_with_mlmg() 19 4.275 4.287 4.299 12.24% Castro::reset_internal_energy(MultiFab) 80 4.203 4.212 4.266 12.14% MLMG::solve() 18 4.192 4.198 4.2 11.96% Castro::evaluate_source_change() 196 3.74 3.75 3.795 10.80% MLMG::oneIter() 128 3.685 3.688 3.69 10.50% MLMG::mgVcycle() 128 3.631 3.634 3.636 10.35% Castro::do_old_sources() 40 3.373 3.406 3.454 9.83% amrex::Copy() 766 3.272 3.279 3.317 9.44% Castro::initialize_do_advance() 20 3.22 3.245 3.27 9.31% Castro::enforce_min_density() 79 3.114 3.121 3.152 8.97% MLCellLinOp::smooth() 1544 1.992 2.149 2.602 7.41% Castro::do_new_sources() 17 2.438 2.48 2.518 7.17% FabArray::ParallelCopy() 1469 1.272 2.053 2.402 6.84% MultiFab::AddProduct() 3598 2.23 2.235 2.262 6.44% AmrLevel::FillPatch() 187 1.95 2.064 2.121 6.04% MLPoisson::Fsmooth() 6144 1.513 1.626 1.95 5.55% FabArray::ParallelCopy_finish() 1469 0.4321 1.194 1.537 4.38% Amr::init() 1 1.531 1.531 1.532 4.36% Amr::initialInit() 1 1.412 1.412 1.412 4.02% FabArray::setVal() 5232 1.369 1.377 1.392 3.96% MultiFab::Subtract() 159 1.358 1.362 1.373 3.91% Amr::InitAmr() 1 1.342 1.344 1.35 3.84% FillPatchIterator::Initialize 187 1.144 1.263 1.323 3.77% FillPatchIterator::FillFromLevel0() 187 1.078 1.222 1.29 3.67% FillPatchSingleLevel 187 1.077 1.221 1.289 3.67% MLMG::mgVcycle_up::1 128 0.2781 0.7544 1.23 3.50% MLMG::mgVcycle_down::0 128 1.112 1.116 1.121 3.19% MLMG::addInterpCorrection() 768 0.1453 0.8416 1.056 3.01% FabArray::ParallelCopy_nowait() 1469 0.8041 0.8549 0.9197 2.62% MLCellLinOp::apply() 986 0.6784 0.7085 0.8146 2.32% Amr::FinalizeInit() 1 0.6447 0.7792 0.8095 2.30% MLCellLinOp::applyBC() 7148 0.5776 0.6401 0.7994 2.28% MLMG::mgVcycle_up::0 128 0.7803 0.7817 0.783 2.23% MLMG::mgVcycle_up::2 128 0.01546 0.3146 0.7704 2.19% Amr::InitializeInit() 1 0.6026 0.6329 0.7675 2.18% Amr::defBaseLevel() 1 0.6026 0.6329 0.7675 2.18% FabArray::sum() 3598 0.7476 0.7497 0.7586 2.16% Castro::expand_state() 20 0.67 0.7124 0.7494 2.13% Gravity::fill_multipole_BCs() 19 0.6405 0.6804 0.707 2.01% MLPoisson::Fapply() 986 0.5596 0.5823 0.6534 1.86% MLMG:computeResOfCorrection() 768 0.4631 0.498 0.5943 1.69% MLCellLinOp::correctionResidual() 768 0.4613 0.4961 0.5923 1.69% Castro::construct_old_source() 240 0.5686 0.5707 0.5775 1.64% check_for_negative_density() 20 0.226 0.482 0.5679 1.62% FabArray::FillBoundary() 6380 0.4046 0.4657 0.5298 1.51% Castro::construct_new_source() 102 0.4696 0.4718 0.4782 1.36% MLMG::mgVcycle_up::3 128 0.01545 0.05054 0.4523 1.29% Castro::Castro() 1 0.4076 0.4122 0.4449 1.27% Castro::post_init() 1 0.431 0.4323 0.4332 1.23% Castro::initMFs() 1 0.3888 0.3933 0.4246 1.21% MLCellLinOp::solutionResidual() 147 0.4002 0.4115 0.418 1.19% MLMG::mgVcycle_down::1 128 0.3703 0.3774 0.3872 1.10% Castro::reset_internal_energy(Fab) 160 0.3699 0.3722 0.3804 1.08% Castro::initData() 1 0.1931 0.219 0.3544 1.01% MLMG::computeResidual() 128 0.3479 0.3496 0.3521 1.00% Castro::construct_old_gravity_source() 40 0.3235 0.3247 0.3304 0.94% Gravity::multilevel_solve_for_new_phi() 1 0.3284 0.3285 0.3292 0.94% Gravity::actual_multilevel_solve() 1 0.3284 0.3285 0.3291 0.94% FabArray::Saxpy() 111 0.3229 0.3245 0.3286 0.94% FillBoundary_finish() 6380 0.2205 0.2488 0.2918 0.83% Amr::bldFineLevels() 1 0.1097 0.2436 0.2742 0.78% Castro::construct_new_gravity_source() 17 0.2631 0.2644 0.2688 0.77% FillBoundary_nowait() 6380 0.1498 0.1995 0.2547 0.72% Amr::writePlotFile() 2 0.2223 0.2224 0.2224 0.63% FabArray::Xpay() 915 0.1917 0.199 0.2191 0.62% MLMG::mgVcycle_down::2 128 0.01797 0.1107 0.2108 0.60% Castro::post_regrid() 2 0.2057 0.2066 0.2079 0.59% Castro::construct_old_gravity() 20 0.1528 0.1797 0.1986 0.57% Gravity::get_old_grav_vector() 20 0.1527 0.1796 0.1985 0.57% Castro::create_source_corrector() 17 0.1711 0.1795 0.1949 0.55% MLMG::mgVcycle_down::3 128 0.01752 0.03911 0.1901 0.54% Amr::grid_places() 2 0.005659 0.1402 0.1705 0.49% AmrMesh::MakeNewGrids() 2 0.0056 0.1401 0.1704 0.49% Castro::derive() 18 0.168 0.1688 0.1702 0.48% AmrLevel::derive() 18 0.1679 0.1687 0.1701 0.48% Castro::check_for_nan() 37 0.1651 0.166 0.1682 0.48% MultiFab::contains_nan() 37 0.1649 0.1659 0.1681 0.48% TagBoxArray::mapPRD 2 0.001692 0.1113 0.1665 0.47% Castro::post_timestep() 1 0.1528 0.1552 0.1577 0.45% Castro::retry_advance_ctu() 11 0.1348 0.1376 0.1446 0.41% Castro::normalize_species() 79 0.1394 0.1401 0.1428 0.41% Gravity::get_new_grav_vector() 18 0.1277 0.1281 0.13 0.37% Castro::initialize_advance() 1 0.1241 0.1256 0.1278 0.36% TagBoxArray::collate() 2 0.0005163 0.02553 0.1264 0.36% amrex::average_down 768 0.07177 0.08618 0.123 0.35% Castro::apply_source_to_state() 37 0.1077 0.1086 0.1104 0.31% Amr::regrid() 1 0.1057 0.1062 0.1069 0.30% MLMG::mgVcycle_bottom 128 0.0004598 0.007091 0.1064 0.30% MLMG::actualBottomSolve() 128 0.0001844 0.006807 0.1061 0.30% MLMG::mgVcycle_down::4 128 0.0175 0.02327 0.1055 0.30% StateDataPhysBCFunct::() 187 0.001013 0.05201 0.1038 0.30% MLCGSolver::bicgstab 8 0 0.006487 0.1038 0.30% StateData::FillBoundary(geom) 193 0 0.0508 0.1025 0.29% MLMG::mgVcycle_up::4 128 0.01549 0.01959 0.07703 0.22% Castro::save_data_for_retry() 3 0.07057 0.07236 0.0769 0.22% Castro::sum_integrated_quantities() 2 0.0737 0.07387 0.07389 0.21% MLMG::computeMLResidual() 18 0.04833 0.05923 0.06623 0.19% MLPoisson::define() 19 0.03501 0.04915 0.06411 0.18% MLMG::prepareForSolve() 18 0.0449 0.04719 0.05357 0.15% amrex::Add() 128 0.05289 0.05312 0.05347 0.15% MLMG::mgVcycle_down::5 128 0.01746 0.01976 0.05025 0.14% VisMF::Write(FabArray) 4 0.03677 0.03812 0.04993 0.14% FabArray::setDomainBndry() 187 0.0006546 0.02515 0.04983 0.14% MLMG::ResNormInf() 146 0.0415 0.046 0.04809 0.14% Gravity::test_composite_phi() 1 0.04127 0.04196 0.04203 0.12% FabArray::norminf() 244 0.03616 0.03663 0.04114 0.12% MLMG::mgVcycle_up::5 128 0.01556 0.01733 0.04038 0.11% Castro::problem_post_timestep() 2 0.02595 0.02852 0.03122 0.09% Castro::wd_update() 2 0.02593 0.02849 0.0312 0.09% MLCellLinOp::defineAuxData() 19 0.02227 0.02428 0.03022 0.09% Castro::volWgtSum() 46 0.02839 0.02854 0.02888 0.08% MLMG::getGradSolution() 18 0.02655 0.02698 0.02771 0.08% Other 10929 0.1822 0.2322 0.3296 0.94% -------------------------------------------------------------------------------------------- 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 434 KiB 436 KiB 443 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 1948 B 1975 B 1991 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 1733 KiB 1733 KiB 1733 KiB 2329 KiB 2329 KiB 2329 KiB Castro::construct_ctu_hydro_source() 6400 434 KiB 436 KiB 443 KiB 1176 KiB 1176 KiB 1176 KiB StateData::define() 192 1119 KiB 1119 KiB 1119 KiB 1164 KiB 1164 KiB 1164 KiB Castro::initMFs() 256 1097 KiB 1098 KiB 1098 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 737 KiB 737 KiB 737 KiB 803 KiB 803 KiB 803 KiB Amr::writePlotFile() 64 4504 B 4542 B 4772 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 70 KiB 70 KiB 71 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 17 KiB 18 KiB 19 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 475 KiB 475 KiB 475 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2957 B 4800 B 6205 B 225 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1259 B 1266 B 1276 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 29 KiB 29 KiB 30 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 143 KiB 143 KiB 143 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 136 KiB 136 KiB 136 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 793 B 796 B 806 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 671 B 673 B 680 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 1948 B 1975 B 1991 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 19 B 20 B 28 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 61 B 62 B 63 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 4 B 4 B 5 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 166 B 167 B 169 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 200 B 235 B 260 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 310 B 313 B 318 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 1742 B 1745 B 1767 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 153 B 153 B 153 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 2126 B 2149 B 2160 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 558 B 624 B 866 B 4608 B 5151 B 7152 B BndryData::define() 4864 586 B 586 B 587 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 19 B 27 B 37 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 0 B 7 B 12 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 1 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-22-gf9456be2b517) finalized