Initializing AMReX (26.05-41-gba3cb0c68650)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.05-41-gba3cb0c68650) initialized Starting run at 06:44:04 UTC on 2026-05-14. Successfully read inputs file ... Castro git describe: 26.05-3-g308427cb3 AMReX git describe: 26.05-41-gba3cb0c68 Microphysics git describe: 26.05-1-g106dcbb2 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.190974096 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001133532 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010481575 Gravity::multilevel_solve_for_new_phi() time = 0.105867042 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.011162867 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.047126816 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004118124 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.004029836 on level 0 Castro::do_old_sources() time = 0.008797986 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.054686407 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 = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Castro::construct_old_gravity() time = 0.00361211 on level 0 Castro::do_old_sources() time = 0.007603864 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.053767559 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.007792036 on level 0 Castro::do_old_sources() time = 0.008060279 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.051418961 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.004386379 on level 0 Castro::do_old_sources() time = 0.00816712 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.050403909 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.300011551 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.014186129 Gravity::solve_for_phi() time = 0.093514836 on level 0 Castro::construct_new_gravity() time = 0.095574146 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.05062523 on level 0 Castro::construct_old_gravity() time = 0.002703313 on level 0 Castro::do_old_sources() time = 0.007545816 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.050436694 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.298239547 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.014481331 Gravity::solve_for_phi() time = 0.095306443 on level 0 Castro::construct_new_gravity() time = 0.097751679 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.050365039 on level 0 Castro::construct_old_gravity() time = 0.003116326 on level 0 Castro::do_old_sources() time = 0.007780386 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.050406537 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.302856093 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.01410974 Gravity::solve_for_phi() time = 0.090599271 on level 0 Castro::construct_new_gravity() time = 0.092649801 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.050378166 on level 0 Castro::construct_old_gravity() time = 0.002753564 on level 0 Castro::do_old_sources() time = 0.007533712 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.05029048 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.307182058 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.014622467 Gravity::solve_for_phi() time = 0.092110678 on level 0 Castro::construct_new_gravity() time = 0.094140248 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.050614141 on level 0 Castro::construct_old_gravity() time = 0.003824876 on level 0 Castro::do_old_sources() time = 0.008341787 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.05031959 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.313830993 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.015120169 Gravity::solve_for_phi() time = 0.095106736 on level 0 Castro::construct_new_gravity() time = 0.09716749 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.050369302 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.00274262 on level 0 Castro::do_old_sources() time = 0.007560812 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.050316926 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.30799692 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.014452739 Gravity::solve_for_phi() time = 0.095117735 on level 0 Castro::construct_new_gravity() time = 0.097134022 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.053752943 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.002757052 on level 0 Castro::do_old_sources() time = 0.007565916 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.050311275 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.31608899 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.015452383 Gravity::solve_for_phi() time = 0.090807381 on level 0 Castro::construct_new_gravity() time = 0.092900571 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.050893254 on level 0 Castro::construct_old_gravity() time = 0.002846404 on level 0 Castro::do_old_sources() time = 0.008140353 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.050190459 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.31071396 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.014990972 Gravity::solve_for_phi() time = 0.092784009 on level 0 Castro::construct_new_gravity() time = 0.094820478 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.050983926 on level 0 Castro::construct_old_gravity() time = 0.002914449 on level 0 Castro::do_old_sources() time = 0.007603943 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.050230532 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.314523066 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.015288225 Gravity::solve_for_phi() time = 0.091025126 on level 0 Castro::construct_new_gravity() time = 0.093077224 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.050935017 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.002710081 on level 0 Castro::do_old_sources() time = 0.007585138 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.050236033 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.309838232 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.015334894 Gravity::solve_for_phi() time = 0.091385011 on level 0 Castro::construct_new_gravity() time = 0.093457773 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.05058074 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.003249004 on level 0 Castro::do_old_sources() time = 0.008130813 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.050639933 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.313463627 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.014172551 Gravity::solve_for_phi() time = 0.09075677 on level 0 Castro::construct_new_gravity() time = 0.092799331 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.049084156 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.00274487 on level 0 Castro::do_old_sources() time = 0.007459021 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.048680365 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.306977256 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.014729478 Gravity::solve_for_phi() time = 0.092244071 on level 0 Castro::construct_new_gravity() time = 0.094522325 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.050555168 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.002784882 on level 0 Castro::do_old_sources() time = 0.008675573 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.050235919 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.310148097 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.013589443 Gravity::solve_for_phi() time = 0.090282093 on level 0 Castro::construct_new_gravity() time = 0.092409144 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.050495663 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.003960346 on level 0 Castro::do_old_sources() time = 0.008875658 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.048916013 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.310684506 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.013836066 Gravity::solve_for_phi() time = 0.091675409 on level 0 Castro::construct_new_gravity() time = 0.093718568 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.049291417 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.004296412 on level 0 Castro::do_old_sources() time = 0.008898876 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.048592032 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.326171008 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.014099772 Gravity::solve_for_phi() time = 0.093282867 on level 0 Castro::construct_new_gravity() time = 0.09532845 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.054258033 on level 0 Castro::construct_old_gravity() time = 0.008693744 on level 0 Castro::do_old_sources() time = 0.008320708 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.050687547 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.319267204 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.014334716 Gravity::solve_for_phi() time = 0.092985147 on level 0 Castro::construct_new_gravity() time = 0.095041407 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.049205371 on level 0 Castro::construct_old_gravity() time = 0.003526592 on level 0 Castro::do_old_sources() time = 0.007479858 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.050633977 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.309094234 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002882557239 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.013942515 Gravity::solve_for_phi() time = 0.090102492 on level 0 Castro::construct_new_gravity() time = 0.093186633 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.049267091 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: 11.42855452 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.046964092 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004036489 seconds Ending run at 06:44:18 UTC on 2026-05-14. Run time = 14.05752105 Run time without initialization = 11.48047907 Average number of zones advanced per microsecond: 0.003 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 14.0479 ... 14.0517 ... 14.0576 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.1037 5.1947 5.2960 37.67% Amr::InitAmr() 1 1.654 1.658 1.664 11.84% Castro::computeTemp() 80 0.7024 0.7985 0.9374 6.67% amrex::Copy() 766 0.81 0.8223 0.8486 6.04% FabArray::ParallelCopy_finish() 1469 0.2706 0.5923 0.8222 5.85% MultiFab::AddProduct() 3598 0.5199 0.5279 0.5522 3.93% MLPoisson::Fsmooth() 6144 0.3293 0.3574 0.4323 3.08% Castro::initMFs() 1 0.3816 0.3888 0.4068 2.89% FabArray::setVal() 5232 0.3656 0.3727 0.3803 2.71% Castro::reset_internal_energy(Fab) 160 0.345 0.3503 0.3575 2.54% MultiFab::Subtract() 159 0.3422 0.3474 0.357 2.54% check_for_negative_density() 20 0.05463 0.1851 0.2883 2.05% Castro::do_advance_ctu() 20 0.2587 0.2644 0.2748 1.95% FabArray::ParallelCopy_nowait() 1469 0.234 0.2509 0.2685 1.91% Gravity::fill_multipole_BCs() 19 0.1933 0.241 0.2641 1.88% FabArray::sum() 3598 0.185 0.1878 0.1949 1.39% FillBoundary_finish() 6380 0.1223 0.1592 0.1894 1.35% Castro::initData() 1 0.003168 0.0303 0.1797 1.28% TagBoxArray::collate() 2 0.0002066 0.02984 0.1625 1.16% FillBoundary_nowait() 6380 0.09163 0.1184 0.1493 1.06% MLCellLinOp::applyBC() 7148 0.07369 0.08476 0.1309 0.93% Castro::normalize_species() 79 0.1176 0.1193 0.1221 0.87% MLPoisson::Fapply() 986 0.08538 0.09025 0.1051 0.75% Castro::construct_old_gravity_source() 40 0.09525 0.09681 0.09859 0.70% FabArray::Saxpy() 111 0.07927 0.08062 0.0868 0.62% Castro::construct_new_gravity_source() 17 0.0672 0.06875 0.07357 0.52% FabArray::Xpay() 915 0.04854 0.05166 0.05747 0.41% MultiFab::contains_nan() 37 0.0457 0.04646 0.04813 0.34% MLMG::addInterpCorrection() 768 0.0321 0.03482 0.04229 0.30% AmrLevel::derive() 18 0.03851 0.03904 0.0399 0.28% StateData::FillBoundary(geom) 193 0 0.01842 0.0361 0.26% amrex::average_down 768 0.02444 0.0268 0.03368 0.24% Castro::reset_internal_energy(MultiFab) 80 0.02496 0.02571 0.02616 0.19% VisMF::Write(FabArray) 4 0.01236 0.01356 0.02372 0.17% MLCellLinOp::defineAuxData() 19 0.01631 0.01738 0.02107 0.15% Castro::construct_old_source() 240 0.01935 0.01979 0.0209 0.15% FabArray::FillBoundary() 6380 0.01664 0.01718 0.02032 0.14% Castro::enforce_min_density() 79 0.01957 0.01991 0.02031 0.14% FillPatchIterator::Initialize 187 0.01568 0.01659 0.01819 0.13% FabArrayBase::getFB() 6380 0.01351 0.0141 0.01723 0.12% Amr::writePlotFile() 2 0.006395 0.01546 0.01698 0.12% MLCGSolver::bicgstab 8 0 0.001056 0.01689 0.12% MLMG::ResNormInf() 146 0.006038 0.01302 0.01684 0.12% Castro::wd_update() 2 0.006303 0.01282 0.0167 0.12% FabArray::setDomainBndry() 187 0.0006697 0.008708 0.01636 0.12% Castro::evaluate_source_change() 196 0.01441 0.0149 0.01606 0.11% Castro::construct_new_source() 102 0.01505 0.01532 0.01588 0.11% FabArray::norminf() 244 0.01083 0.01127 0.01466 0.10% amrex::Add() 128 0.01355 0.01373 0.01387 0.10% MLCellLinOp::smooth() 1544 0.01038 0.01083 0.01219 0.09% MLMG::prepareForSolve() 18 0.01031 0.01073 0.01194 0.08% Gravity::actual_solve_with_mlmg() 19 0.007452 0.008008 0.00975 0.07% FabArrayBase::getCPC() 2045 0.008803 0.00904 0.009519 0.07% Castro::sum_integrated_quantities() 2 0.002738 0.007924 0.008725 0.06% MLMG::solve() 18 0.000814 0.004167 0.007049 0.05% Gravity::get_old_grav_vector() 20 0.00611 0.006192 0.006388 0.05% Gravity::get_new_grav_vector() 18 0.005631 0.005775 0.005939 0.04% Castro::initialize_do_advance() 20 0.004659 0.004957 0.005565 0.04% MLCellLinOp::apply() 986 0.00235 0.002609 0.005564 0.04% FabArray::ParallelCopy() 1469 0.004185 0.004312 0.004606 0.03% AmrLevel::FillPatch() 187 0.003091 0.003265 0.003584 0.03% MLLinOp::defineGrids() 19 0.002566 0.002778 0.003278 0.02% MLCellLinOp::correctionResidual() 768 0.002284 0.002356 0.00244 0.02% MLMG::mgVcycle() 128 0.002154 0.002233 0.002391 0.02% MLMG:computeResOfCorrection() 768 0.001806 0.00186 0.002008 0.01% Castro::save_data_for_retry() 3 0.001317 0.001364 0.001552 0.01% StateDataPhysBCFunct::() 187 0.001023 0.001212 0.001345 0.01% MLMG::actualBottomSolve() 128 0.0001907 0.0002625 0.00128 0.01% FillPatchIterator::FillFromLevel0() 187 0.00089 0.0009281 0.001011 0.01% Castro::subcycle_advance_ctu() 1 0.000241 0.000313 0.0009399 0.01% main() 1 0.0007827 0.0008348 0.00092 0.01% Gravity::solve_for_phi() 17 0.000687 0.0007318 0.0007774 0.01% FillPatchSingleLevel 187 0.0006413 0.0006636 0.0006828 0.00% MLMG::mgVcycle_down::0 128 0.0006311 0.0006419 0.0006762 0.00% MLMG::mgVcycle_down::5 128 0.0005795 0.0005975 0.0006626 0.00% MLMG::mgVcycle_down::1 128 0.0006055 0.0006219 0.0006611 0.00% MLMG::mgVcycle_down::3 128 0.0005539 0.0005879 0.0006553 0.00% MLMG::mgVcycle_down::2 128 0.0005778 0.0006137 0.0006449 0.00% MLMG::mgVcycle_down::4 128 0.0005466 0.0005689 0.000633 0.00% AmrMesh::MakeNewGrids() 2 0.0005472 0.0005753 0.0006078 0.00% MLCellLinOp::solutionResidual() 147 0.0005101 0.0005253 0.0005508 0.00% Castro::clean_state() 79 0.0004756 0.0004989 0.0005298 0.00% Castro::Castro() 1 0.0003885 0.000421 0.0004929 0.00% Castro::do_old_sources() 40 0.0004392 0.0004668 0.000492 0.00% MLMG::oneIter() 128 0.0004274 0.0004413 0.0004859 0.00% MLMG::mgVcycle_up::1 128 0.0004071 0.0004212 0.0004713 0.00% MLMG::mgVcycle_up::0 128 0.0004048 0.0004185 0.0004679 0.00% MLMG::mgVcycle_up::4 128 0.0003756 0.0003893 0.0004429 0.00% MLMG::mgVcycle_up::2 128 0.0003839 0.0004095 0.0004368 0.00% MLMG::mgVcycle_up::3 128 0.0003701 0.0003865 0.0004251 0.00% MLMG::mgVcycle_up::5 128 0.0003634 0.0003763 0.0003987 0.00% Castro::construct_new_gravity() 17 0.0001937 0.0002155 0.0003299 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002211 0.0002481 0.0003173 0.00% MLMG::computeResidual() 128 0.000289 0.0002963 0.0003062 0.00% MLMG::mgVcycle_bottom 128 0.0002696 0.0002787 0.0003019 0.00% Castro::problem_post_timestep() 2 1.386e-05 2.99e-05 0.0002427 0.00% Castro::do_new_sources() 17 0.0002104 0.000225 0.0002392 0.00% Amr::coarseTimeStep() 1 0.0001296 0.0001609 0.0001661 0.00% MLPoisson::define() 19 0.000127 0.0001376 0.0001566 0.00% Castro::construct_old_gravity() 20 0.0001237 0.0001375 0.0001506 0.00% Castro::apply_source_to_state() 37 0.0001182 0.0001229 0.0001447 0.00% Castro::check_for_nan() 37 0.0001214 0.0001304 0.0001434 0.00% Castro::retry_advance_ctu() 11 7.41e-05 8.874e-05 0.0001371 0.00% TagBoxArray::mapPRD 2 0.0001184 0.0001241 0.0001339 0.00% Gravity::test_composite_phi() 1 0.0001075 0.0001133 0.0001264 0.00% Castro::initialize_advance() 1 0.0001077 0.0001106 0.0001137 0.00% Amr::defBaseLevel() 1 9.179e-05 0.0001005 0.0001064 0.00% MLLinOp::define() 19 8.797e-05 9.17e-05 9.956e-05 0.00% Castro::create_source_corrector() 17 5.727e-05 6.275e-05 7.531e-05 0.00% Castro::expand_state() 20 5.695e-05 6.173e-05 7.208e-05 0.00% MLMG::computeMLResidual() 18 4.808e-05 5.545e-05 7.095e-05 0.00% Amr::grid_places() 2 5.371e-05 5.824e-05 6.964e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.824e-05 6.06e-05 6.56e-05 0.00% Castro::derive() 18 4.659e-05 4.917e-05 5.877e-05 0.00% Amr::FinalizeInit() 1 8.812e-06 1.211e-05 5.188e-05 0.00% Amr::regrid() 1 1.403e-05 1.89e-05 4.261e-05 0.00% Gravity::actual_multilevel_solve() 1 3.133e-05 3.243e-05 3.6e-05 0.00% Amr::timeStep() 1 1.497e-05 1.655e-05 2.248e-05 0.00% Amr::InitializeInit() 1 1.023e-05 1.397e-05 1.642e-05 0.00% Castro::advance() 1 8.103e-06 1.06e-05 1.598e-05 0.00% Castro::post_init() 1 1.286e-05 1.365e-05 1.488e-05 0.00% Castro::post_timestep() 1 5.843e-06 7.045e-06 1.391e-05 0.00% Castro::post_regrid() 2 1.078e-05 1.179e-05 1.305e-05 0.00% Amr::bldFineLevels() 1 7.256e-06 8.494e-06 1.041e-05 0.00% Amr::init() 1 6.803e-06 7.51e-06 8.614e-06 0.00% Amr::initialInit() 1 3.27e-06 3.434e-06 3.858e-06 0.00% Other 2530 0.04771 0.06061 0.09904 0.70% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.05 14.05 14.06 100.00% Amr::coarseTimeStep() 1 11.48 11.48 11.48 81.66% Amr::timeStep() 1 11.43 11.43 11.43 81.30% Castro::advance() 1 11.36 11.37 11.37 80.87% Castro::subcycle_advance_ctu() 1 11.32 11.32 11.33 80.57% Castro::do_advance_ctu() 20 11.28 11.29 11.29 80.31% Castro::construct_ctu_hydro_source() 20 5.99 6.014 6.054 43.06% Castro::clean_state() 79 2.9 3.02 3.178 22.61% Castro::computeTemp() 80 2.018 2.131 2.282 16.23% Gravity::solve_phi_with_mlmg() 19 1.607 1.659 1.685 11.99% Amr::InitAmr() 1 1.654 1.658 1.664 11.84% Castro::construct_new_gravity() 17 1.528 1.579 1.604 11.41% Gravity::solve_for_phi() 17 1.49 1.542 1.568 11.15% Gravity::actual_solve_with_mlmg() 19 1.411 1.415 1.419 10.09% Castro::reset_internal_energy(MultiFab) 80 1.312 1.332 1.367 9.72% MLMG::solve() 18 1.353 1.359 1.361 9.68% MLMG::oneIter() 128 1.214 1.218 1.222 8.69% Castro::initialize_do_advance() 20 1.153 1.184 1.21 8.61% MLMG::mgVcycle() 128 1.2 1.204 1.208 8.59% Castro::do_old_sources() 40 1.107 1.134 1.159 8.24% FabArray::ParallelCopy() 1469 0.547 0.8688 1.094 7.78% Castro::evaluate_source_change() 196 0.9202 0.9349 0.9747 6.93% Amr::init() 1 0.9128 0.9128 0.9129 6.49% Amr::initialInit() 1 0.8611 0.8611 0.8612 6.13% Castro::do_new_sources() 17 0.7652 0.8112 0.8538 6.07% amrex::Copy() 766 0.81 0.8223 0.8486 6.04% Castro::enforce_min_density() 79 0.7823 0.7941 0.8248 5.87% FabArray::ParallelCopy_finish() 1469 0.2706 0.5923 0.8222 5.85% MLCellLinOp::smooth() 1544 0.6343 0.688 0.8218 5.85% AmrLevel::FillPatch() 187 0.6105 0.694 0.7707 5.48% Amr::InitializeInit() 1 0.4548 0.4903 0.6445 4.58% Amr::defBaseLevel() 1 0.4548 0.4903 0.6444 4.58% FillPatchIterator::Initialize 187 0.4055 0.4912 0.5701 4.06% MultiFab::AddProduct() 3598 0.5199 0.5279 0.5522 3.93% FillPatchIterator::FillFromLevel0() 187 0.3721 0.4659 0.5468 3.89% FillPatchSingleLevel 187 0.3712 0.465 0.5459 3.88% MLCellLinOp::applyBC() 7148 0.3515 0.3943 0.477 3.39% MLMG::mgVcycle_up::1 128 0.1081 0.283 0.459 3.26% MLPoisson::Fsmooth() 6144 0.3293 0.3574 0.4323 3.08% Castro::Castro() 1 0.387 0.3941 0.4122 2.93% Castro::initMFs() 1 0.3818 0.389 0.4069 2.89% Amr::FinalizeInit() 1 0.2166 0.3708 0.4063 2.89% MLMG::addInterpCorrection() 768 0.06504 0.3193 0.3949 2.81% FabArray::setVal() 5232 0.3656 0.3727 0.3803 2.71% Castro::reset_internal_energy(Fab) 160 0.345 0.3503 0.3575 2.54% MultiFab::Subtract() 159 0.3422 0.3474 0.357 2.54% FabArray::FillBoundary() 6380 0.2753 0.3096 0.3567 2.54% MLMG::mgVcycle_up::2 128 0.01533 0.1272 0.2985 2.12% FabArray::ParallelCopy_nowait() 1469 0.2546 0.2722 0.2902 2.06% check_for_negative_density() 20 0.05463 0.1851 0.2883 2.05% MLMG::mgVcycle_down::0 128 0.2647 0.2705 0.2771 1.97% Gravity::fill_multipole_BCs() 19 0.1954 0.243 0.2661 1.89% Castro::expand_state() 20 0.2028 0.2336 0.2607 1.85% Castro::initData() 1 0.06599 0.09449 0.2447 1.74% Amr::bldFineLevels() 1 0.03833 0.1922 0.2271 1.62% MLMG::mgVcycle_up::0 128 0.2073 0.208 0.2086 1.48% MLCellLinOp::apply() 986 0.1551 0.1653 0.2077 1.48% MLMG::mgVcycle_up::3 128 0.01534 0.02948 0.1988 1.41% FabArray::sum() 3598 0.185 0.1878 0.1949 1.39% Amr::grid_places() 2 0.00243 0.1566 0.1921 1.37% AmrMesh::MakeNewGrids() 2 0.002375 0.1565 0.192 1.37% TagBoxArray::mapPRD 2 0.0007817 0.1253 0.1904 1.35% FillBoundary_finish() 6380 0.1223 0.1592 0.1894 1.35% FillBoundary_nowait() 6380 0.1057 0.1332 0.1644 1.17% Castro::construct_old_source() 240 0.1581 0.1609 0.1641 1.17% TagBoxArray::collate() 2 0.0002066 0.02984 0.1625 1.16% MLMG:computeResOfCorrection() 768 0.1173 0.1285 0.1524 1.08% MLCellLinOp::correctionResidual() 768 0.1155 0.1266 0.1505 1.07% Castro::post_init() 1 0.1407 0.1435 0.145 1.03% MLMG::mgVcycle_down::1 128 0.1281 0.136 0.1428 1.02% Castro::construct_new_source() 102 0.1195 0.1217 0.126 0.90% Castro::normalize_species() 79 0.1176 0.1193 0.1221 0.87% Gravity::multilevel_solve_for_new_phi() 1 0.1058 0.1058 0.1059 0.75% Gravity::actual_multilevel_solve() 1 0.1057 0.1058 0.1058 0.75% MLPoisson::Fapply() 986 0.08538 0.09025 0.1051 0.75% MLCellLinOp::solutionResidual() 147 0.08548 0.09123 0.1009 0.72% Castro::construct_old_gravity_source() 40 0.09525 0.09681 0.09859 0.70% Castro::create_source_corrector() 17 0.05801 0.07421 0.09793 0.70% Amr::writePlotFile() 2 0.0945 0.0946 0.0946 0.67% FabArray::Saxpy() 111 0.07927 0.08062 0.0868 0.62% MLMG::mgVcycle_bottom 128 0.0004646 0.005538 0.08153 0.58% MLMG::actualBottomSolve() 128 0.0001907 0.005259 0.08123 0.58% MLMG::mgVcycle_down::2 128 0.01788 0.04756 0.08035 0.57% MLCGSolver::bicgstab 8 0 0.004949 0.07918 0.56% MLMG::computeResidual() 128 0.07368 0.07502 0.07858 0.56% Castro::construct_new_gravity_source() 17 0.0672 0.06875 0.07357 0.52% Castro::post_regrid() 2 0.0692 0.07065 0.07348 0.52% Castro::derive() 18 0.06724 0.06855 0.07325 0.52% AmrLevel::derive() 18 0.0672 0.0685 0.07319 0.52% Castro::construct_old_gravity() 20 0.0464 0.05725 0.06983 0.50% Gravity::get_old_grav_vector() 20 0.04626 0.05711 0.06969 0.50% MLMG::mgVcycle_down::3 128 0.01731 0.02376 0.06647 0.47% Castro::post_timestep() 1 0.06014 0.06246 0.06565 0.47% amrex::average_down 768 0.03415 0.04143 0.05832 0.41% FabArray::Xpay() 915 0.04854 0.05166 0.05747 0.41% Castro::check_for_nan() 37 0.04582 0.04659 0.04825 0.34% MultiFab::contains_nan() 37 0.0457 0.04646 0.04813 0.34% MLPoisson::define() 19 0.03309 0.03796 0.04323 0.31% Castro::initialize_advance() 1 0.04029 0.041 0.04192 0.30% MLMG::mgVcycle_down::4 128 0.01728 0.01918 0.04105 0.29% Amr::regrid() 1 0.03599 0.03679 0.03887 0.28% Castro::retry_advance_ctu() 11 0.03505 0.03644 0.03836 0.27% StateDataPhysBCFunct::() 187 0.001023 0.01963 0.03739 0.27% Gravity::get_new_grav_vector() 18 0.03512 0.0366 0.03731 0.27% StateData::FillBoundary(geom) 193 0 0.01842 0.0361 0.26% MLMG::mgVcycle_up::4 128 0.01535 0.01674 0.03252 0.23% Castro::apply_source_to_state() 37 0.02665 0.0271 0.03032 0.22% MLMG::mgVcycle_down::5 128 0.0173 0.01851 0.03013 0.21% MLCellLinOp::defineAuxData() 19 0.02101 0.02282 0.02841 0.20% Castro::sum_integrated_quantities() 2 0.02796 0.02808 0.0281 0.20% MLMG::ResNormInf() 146 0.01662 0.02263 0.02631 0.19% MLMG::mgVcycle_up::5 128 0.01547 0.01644 0.02569 0.18% MLMG::prepareForSolve() 18 0.02062 0.02165 0.0247 0.18% VisMF::Write(FabArray) 4 0.01236 0.01356 0.02372 0.17% MLMG::computeMLResidual() 18 0.01032 0.0158 0.02171 0.15% Castro::save_data_for_retry() 3 0.01841 0.01896 0.02127 0.15% Castro::problem_post_timestep() 2 0.01154 0.0167 0.02 0.14% Castro::wd_update() 2 0.01152 0.01667 0.01998 0.14% FabArrayBase::getFB() 6380 0.01384 0.01479 0.01872 0.13% FabArray::setDomainBndry() 187 0.0006697 0.008708 0.01636 0.12% FabArrayBase::getCPC() 2045 0.01118 0.01256 0.01604 0.11% Gravity::test_composite_phi() 1 0.01493 0.01497 0.015 0.11% FabArray::norminf() 244 0.01083 0.01127 0.01466 0.10% amrex::Add() 128 0.01355 0.01373 0.01387 0.10% MLLinOp::define() 19 0.00503 0.008311 0.01178 0.08% MLLinOp::defineGrids() 19 0.004941 0.008219 0.01169 0.08% Other 2530 0.07875 0.09437 0.1403 1.00% -------------------------------------------------------------------------------------------- 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 384 KiB 391 KiB 399 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 2151 B 2173 B 2191 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 1 B 1 B 3 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 1543 KiB 1544 KiB 1544 KiB 2332 KiB 2332 KiB 2332 KiB Castro::construct_ctu_hydro_source() 6400 385 KiB 392 KiB 399 KiB 1178 KiB 1178 KiB 1178 KiB StateData::define() 192 1026 KiB 1027 KiB 1027 KiB 1166 KiB 1166 KiB 1166 KiB Castro::initMFs() 256 1005 KiB 1006 KiB 1006 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 656 KiB 656 KiB 656 KiB 804 KiB 804 KiB 804 KiB Amr::writePlotFile() 64 4305 B 4405 B 4939 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 49 KiB 50 KiB 51 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 13 KiB 14 KiB 16 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 422 KiB 422 KiB 422 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2472 B 4508 B 6081 B 226 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1382 B 1402 B 1434 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 18 KiB 18 KiB 19 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 127 KiB 127 KiB 127 KiB 158 KiB 158 KiB 158 KiB Castro::buildMetrics() 192 125 KiB 125 KiB 125 KiB 143 KiB 143 KiB 143 KiB MLMG::prepareForSolve() 6102 10 KiB 10 KiB 12 KiB 104 KiB 110 KiB 130 KiB Castro::construct_old_source() 1280 496 B 506 B 533 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 418 B 423 B 447 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 2157 B 2179 B 2197 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 18 B 20 B 37 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 33 KiB 33 KiB 33 KiB 76 KiB 76 KiB 76 KiB Gravity::test_composite_phi() 96 56 B 56 B 57 B 53 KiB 53 KiB 53 KiB Gravity::update_max_rhs() 192 3 B 6 B 8 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 113 B 118 B 120 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 151 B 186 B 229 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 31 KiB 31 KiB 31 KiB 38 KiB 38 KiB 38 KiB Gravity::fill_multipole_BCs() 2432 246 B 252 B 257 B 18 KiB 18 KiB 18 KiB MLCellLinOp::defineAuxData() 11704 812 B 976 B 1822 B 8192 B 9856 B 18 KiB Castro::evaluate_source_change() 6272 1066 B 1083 B 1129 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 123 B 123 B 123 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 1735 B 1794 B 1825 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 0 B 0 B 1 B 16 KiB 16 KiB 16 KiB BndryData::define() 4864 800 B 803 B 804 B 8192 B 8192 B 8192 B FillBoundary_nowait() 91296 34 B 45 B 59 B 3072 B 3840 B 4608 B AmrMesh::MakeNewGrids() 128 0 B 27 B 41 B 3584 B 3584 B 3584 B OwnerMask() 64 0 B 0 B 0 B 3072 B 3072 B 3072 B MLCGSolver::bicgstab 640 0 B 0 B 6 B 0 B 80 B 1280 B MLMG::addInterpCorrection() 5632 1 B 16 B 26 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 0 B 10 B 13 B 1024 B 1024 B 1024 B amrex::average_down 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B MLCellLinOp::applyBC() 114375 0 B 0 B 1 B 256 B 256 B 256 B ----------------------------------------------------------------------------------------------------------------------- AMReX (26.05-41-gba3cb0c68650) finalized