Initializing AMReX (25.06-22-g1d307c2c661a-dirty)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (25.06-22-g1d307c2c661a-dirty) initialized Starting run at 06:48:15 UTC on 2025-07-14. Successfully read inputs file ... Castro git describe: 25.07-5-g5d746e0a1 AMReX git describe: 25.06-22-g1d307c2c6-dirty Microphysics git describe: 25.07-5-gb321ef91 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.165042821 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001187513 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010849025 Gravity::multilevel_solve_for_new_phi() time = 0.108244528 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.011596012 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.051012811 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.007059061 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 = -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 = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.7670713549 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.7670713549 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Invalid density = -0.01295059363 at index 21, 172, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035502 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477044 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.769409 at index 17, 181, 0 Invalid density = -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.004639115 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.064835892 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.009439016 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.065255985 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, 78, 0 Invalid density = -0.1060744788 at index 20, 177, 0 Castro::construct_old_gravity() time = 0.005454351 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.061456427 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.006354036 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.061550729 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.338080506 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.014456306 Gravity::solve_for_phi() time = 0.093366871 on level 0 Castro::construct_new_gravity() time = 0.095624239 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.060327555 on level 0 Castro::construct_old_gravity() time = 0.004683856 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.059287169 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.339201549 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.013951514 Gravity::solve_for_phi() time = 0.090560748 on level 0 Castro::construct_new_gravity() time = 0.092747126 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.05959088 on level 0 Castro::construct_old_gravity() time = 0.002709967 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.059235058 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.337871812 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.013971774 Gravity::solve_for_phi() time = 0.091237582 on level 0 Castro::construct_new_gravity() time = 0.093386684 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.059850586 on level 0 Castro::construct_old_gravity() time = 0.003245918 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.060401323 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.343159253 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.014994896 Gravity::solve_for_phi() time = 0.091581541 on level 0 Castro::construct_new_gravity() time = 0.093731248 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.059727164 on level 0 Castro::construct_old_gravity() time = 0.002663499 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.059483545 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.347365203 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.014592803 Gravity::solve_for_phi() time = 0.091069047 on level 0 Castro::construct_new_gravity() time = 0.093231513 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.059773294 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.002716115 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.059347456 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.34875761 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.01452029 Gravity::solve_for_phi() time = 0.091670297 on level 0 Castro::construct_new_gravity() time = 0.09388171 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.059669862 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.003392645 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.062771955 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.349100491 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.014344992 Gravity::solve_for_phi() time = 0.090891811 on level 0 Castro::construct_new_gravity() time = 0.093026488 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.059754099 on level 0 Castro::construct_old_gravity() time = 0.00285928 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.059880895 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.343663546 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.014591786 Gravity::solve_for_phi() time = 0.091000153 on level 0 Castro::construct_new_gravity() time = 0.093145913 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.059884979 on level 0 Castro::construct_old_gravity() time = 0.002998881 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.059327777 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.348068681 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.014573098 Gravity::solve_for_phi() time = 0.091040641 on level 0 Castro::construct_new_gravity() time = 0.093180009 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.060806926 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.004151663 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.059606151 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.345236277 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.014602865 Gravity::solve_for_phi() time = 0.092300246 on level 0 Castro::construct_new_gravity() time = 0.09444492 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.059808946 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.003204022 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.059433222 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.350033853 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.015208357 Gravity::solve_for_phi() time = 0.091294156 on level 0 Castro::construct_new_gravity() time = 0.093425146 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.059725624 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.003002307 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.059108807 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.343934989 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.0147199 Gravity::solve_for_phi() time = 0.091491207 on level 0 Castro::construct_new_gravity() time = 0.09363961 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.059792002 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.003935816 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.059880511 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.353751463 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.014873851 Gravity::solve_for_phi() time = 0.091207736 on level 0 Castro::construct_new_gravity() time = 0.093394933 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.060096373 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.003791121 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.059859505 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.346812748 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.015475541 Gravity::solve_for_phi() time = 0.093912869 on level 0 Castro::construct_new_gravity() time = 0.096070674 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.060295984 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.003388272 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.065191733 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.348874493 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.015057437 Gravity::solve_for_phi() time = 0.091688051 on level 0 Castro::construct_new_gravity() time = 0.093832631 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.073249873 on level 0 Castro::construct_old_gravity() time = 0.002798266 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.060141076 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.347557179 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.016502196 Gravity::solve_for_phi() time = 0.092904905 on level 0 Castro::construct_new_gravity() time = 0.095070189 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.060448324 on level 0 Castro::construct_old_gravity() time = 0.002777522 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.059626924 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.352798978 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.0026752427 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.01556234 Gravity::solve_for_phi() time = 0.092488256 on level 0 Castro::construct_new_gravity() time = 0.094606721 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.059970199 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.581215356e+29 ener added: 1.581215356e+29 TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -4.341809003e+38 TIME= 0.25 YMOM = 1.164437349e+28 TIME= 0.25 ZMOM = 0 TIME= 0.25 ANG MOM X = 0 TIME= 0.25 ANG MOM Y = 0 TIME= 0.25 ANG MOM Z = -1.469212319e+38 TIME= 0.25 RHO*e = 1.325894512e+50 TIME= 0.25 RHO*K = 3.170941987e+49 TIME= 0.25 RHO*E = 1.642988711e+50 TIME= 0.25 RHO*PHI = -5.449011378e+50 TIME= 0.25 TOTAL ENERGY= -1.081516979e+50 TIME= 0.25 CENTER OF MASS X-LOC = 316424388.8 TIME= 0.25 CENTER OF MASS X-VEL = -170267.2489 TIME= 0.25 CENTER OF MASS Y-LOC = 2.427336677e-06 TIME= 0.25 CENTER OF MASS Y-VEL = 4.566427124e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 57600050.72 TIME= 0.25 MAXIMUM DENSITY = 4458760.788 TIME= 0.25 MAXIMUM T_S / T_E = 0 STEP = 1 TIME = 0.25 DT = 0.25 [STEP 1] Coarse TimeStep time: 12.30968548 [STEP 1] FAB kilobyte spread across MPI nodes: [6998 ... 6998] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.050192926 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.005201538 seconds Ending run at 06:48:30 UTC on 2025-07-14. Run time = 14.63313612 Run time without initialization = 12.36605807 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.6243 ... 14.6270 ... 14.6333 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.7201 5.7573 5.8081 39.69% Amr::InitAmr() 1 1.353 1.356 1.362 9.31% Castro::computeTemp() 80 0.844 0.9601 1.122 7.67% amrex::Copy() 724 0.8002 0.8046 0.8167 5.58% FabArray::ParallelCopy_finish() 1429 0.2691 0.5673 0.7599 5.19% MultiFab::AddProduct() 3439 0.552 0.554 0.5604 3.83% Castro::normalize_species() 79 0.5468 0.5498 0.5549 3.79% Castro::reset_internal_energy(Fab) 160 0.5223 0.5255 0.5325 3.64% MLPoisson::Fsmooth() 6144 0.356 0.3822 0.4599 3.14% Castro::initMFs() 1 0.359 0.3637 0.3768 2.57% FabArray::setVal() 5013 0.3488 0.3522 0.3575 2.44% MultiFab::Subtract() 159 0.3377 0.3383 0.3405 2.33% Gravity::fill_multipole_BCs() 19 0.1973 0.2475 0.2693 1.84% FabArray::ParallelCopy_nowait() 1429 0.2033 0.2189 0.2438 1.67% check_for_negative_density() 20 0.1024 0.1652 0.227 1.55% FabArray::sum() 3439 0.1974 0.1985 0.2014 1.38% Castro::initData() 1 0.003359 0.02793 0.1631 1.11% FillBoundary_finish() 6380 0.1013 0.1286 0.1627 1.11% FillBoundary_nowait() 6380 0.0919 0.1205 0.1561 1.07% TagBoxArray::collate() 2 0.0001979 0.0263 0.1457 1.00% MLCellLinOp::applyBC() 7148 0.06973 0.08023 0.1261 0.86% MLPoisson::Fapply() 986 0.09492 0.0993 0.1141 0.78% Castro::construct_new_gravity_source() 17 0.07294 0.07368 0.07445 0.51% FabArray::Saxpy() 91 0.07198 0.07239 0.07383 0.50% FabArray::Xpay() 915 0.05165 0.05366 0.05925 0.40% Castro::construct_old_gravity_source() 20 0.04892 0.04959 0.05065 0.35% MultiFab::contains_nan() 37 0.04625 0.04644 0.04664 0.32% MLMG::addInterpCorrection() 768 0.03392 0.03643 0.04386 0.30% AmrLevel::derive() 18 0.04104 0.04153 0.04209 0.29% amrex::average_down 768 0.02515 0.02731 0.0338 0.23% StateData::FillBoundary(geom) 143 0 0.01753 0.03375 0.23% VisMF::Write(FabArray) 4 0.01578 0.01735 0.02865 0.20% Castro::reset_internal_energy(MultiFab) 80 0.02311 0.02377 0.02448 0.17% MLCellLinOp::defineAuxData() 19 0.01584 0.01671 0.02004 0.14% Castro::enforce_min_density() 79 0.01894 0.01911 0.01956 0.13% Amr::writePlotFile() 2 0.007802 0.01715 0.01931 0.13% FabArray::FillBoundary() 6380 0.01398 0.01453 0.01746 0.12% MLCGSolver::bicgstab 8 0 0.001032 0.01651 0.11% Castro::construct_new_source() 102 0.01612 0.01624 0.01639 0.11% Castro::evaluate_source_change() 196 0.01405 0.01431 0.01506 0.10% amrex::Add() 128 0.01451 0.01455 0.01468 0.10% FabArray::setDomainBndry() 147 0.0005174 0.007602 0.01444 0.10% FabArray::norminf() 244 0.0107 0.01107 0.01404 0.10% FillPatchIterator::Initialize 147 0.01303 0.01333 0.01359 0.09% FabArrayBase::getFB() 6380 0.0106 0.01105 0.01342 0.09% Castro::wd_update() 2 0.005826 0.0094 0.01302 0.09% MLCellLinOp::smooth() 1544 0.01023 0.01081 0.01224 0.08% MLMG::prepareForSolve() 18 0.01022 0.01064 0.01149 0.08% Castro::construct_old_source() 120 0.01025 0.01031 0.01054 0.07% Gravity::actual_solve_with_mlmg() 19 0.00732 0.007706 0.009157 0.06% FabArrayBase::getCPC() 2005 0.007736 0.007968 0.008253 0.06% MLMG::ResNormInf() 146 0.002757 0.005319 0.00813 0.06% MLMG::solve() 18 0.001079 0.003487 0.007241 0.05% Gravity::get_old_grav_vector() 20 0.006458 0.006597 0.006824 0.05% Gravity::get_new_grav_vector() 18 0.006027 0.006116 0.006282 0.04% Castro::initialize_do_advance() 20 0.004711 0.005083 0.005411 0.04% MLCellLinOp::apply() 986 0.002057 0.002278 0.004864 0.03% FabArray::ParallelCopy() 1429 0.003633 0.003791 0.00409 0.03% Castro::sum_integrated_quantities() 2 0.002614 0.003475 0.003701 0.03% AmrLevel::FillPatch() 147 0.002436 0.002492 0.002653 0.02% MLLinOp::defineGrids() 19 0.002002 0.002123 0.00242 0.02% MLCellLinOp::correctionResidual() 768 0.001907 0.001959 0.002042 0.01% MLMG::mgVcycle() 128 0.001653 0.001717 0.001884 0.01% MLMG:computeResOfCorrection() 768 0.001524 0.001568 0.001674 0.01% Castro::save_data_for_retry() 3 0.001336 0.001403 0.001493 0.01% Castro::subcycle_advance_ctu() 1 0.0002275 0.0003164 0.001228 0.01% MLMG::actualBottomSolve() 128 0.0001337 0.0002035 0.001197 0.01% StateDataPhysBCFunct::() 147 0.0007854 0.0008913 0.0009897 0.01% Gravity::solve_for_phi() 17 0.0006765 0.0007123 0.00093 0.01% main() 1 0.0006767 0.0007933 0.000924 0.01% FillPatchIterator::FillFromLevel0() 147 0.0007044 0.0007189 0.0007568 0.01% MLMG::mgVcycle_down::3 128 0.0004995 0.0005471 0.0006216 0.00% AmrMesh::MakeNewGrids() 2 0.0005359 0.0005641 0.0006213 0.00% MLMG::mgVcycle_down::2 128 0.0005402 0.000577 0.0006181 0.00% MLMG::mgVcycle_down::0 128 0.0005578 0.0005728 0.0006051 0.00% MLMG::mgVcycle_down::1 128 0.0005558 0.0005684 0.0006029 0.00% MLMG::mgVcycle_down::4 128 0.0004799 0.0005118 0.0005929 0.00% MLMG::mgVcycle_down::5 128 0.0005105 0.0005306 0.0005842 0.00% Castro::Castro() 1 0.0004649 0.0005118 0.0005723 0.00% Castro::clean_state() 79 0.000436 0.0004702 0.000525 0.00% FillPatchSingleLevel 147 0.0004656 0.0004878 0.0005232 0.00% MLCellLinOp::solutionResidual() 147 0.0004479 0.0004602 0.0004774 0.00% Castro::do_advance_ctu() 20 0.0003945 0.0004275 0.000466 0.00% MLMG::oneIter() 128 0.0003972 0.0004098 0.0004532 0.00% MLMG::mgVcycle_up::0 128 0.0003731 0.0003852 0.0004284 0.00% MLMG::mgVcycle_up::1 128 0.0003695 0.0003834 0.000418 0.00% Castro::construct_new_gravity() 17 0.0001661 0.0001885 0.0003826 0.00% MLMG::mgVcycle_up::4 128 0.0003229 0.0003352 0.0003783 0.00% MLMG::mgVcycle_up::2 128 0.0003247 0.0003515 0.0003782 0.00% MLMG::mgVcycle_up::3 128 0.0003191 0.0003345 0.0003714 0.00% MLMG::mgVcycle_up::5 128 0.000308 0.000323 0.0003464 0.00% Castro::do_old_sources() 20 0.00028 0.0002971 0.0003082 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002223 0.0002391 0.000288 0.00% MLMG::computeResidual() 128 0.0002366 0.0002449 0.0002607 0.00% Castro::problem_post_timestep() 2 1.12e-05 2.805e-05 0.000249 0.00% MLMG::mgVcycle_bottom 128 0.0002211 0.0002326 0.0002463 0.00% Castro::do_new_sources() 17 0.0002133 0.000223 0.0002376 0.00% Castro::initialize_advance() 1 0.0001444 0.0001617 0.0001862 0.00% MLPoisson::define() 19 0.0001349 0.0001564 0.0001843 0.00% Castro::retry_advance_ctu() 11 6.348e-05 8.085e-05 0.0001761 0.00% Amr::coarseTimeStep() 1 0.0001456 0.0001529 0.000169 0.00% TagBoxArray::mapPRD 2 0.0001226 0.0001312 0.0001405 0.00% Castro::construct_old_gravity() 20 0.0001141 0.0001236 0.0001404 0.00% Gravity::test_composite_phi() 1 0.0001047 0.0001098 0.0001401 0.00% Castro::check_for_nan() 37 0.0001076 0.0001176 0.0001393 0.00% Amr::grid_places() 2 5.288e-05 6.103e-05 0.0001214 0.00% Amr::defBaseLevel() 1 0.0001075 0.0001126 0.0001183 0.00% Amr::FinalizeInit() 1 1.337e-05 2.253e-05 0.000118 0.00% Castro::apply_source_to_state() 37 9.836e-05 0.0001054 0.0001119 0.00% MLLinOp::define() 19 7.865e-05 8.387e-05 9.449e-05 0.00% Castro::expand_state() 20 5.288e-05 5.984e-05 7.649e-05 0.00% Castro::create_source_corrector() 17 5.627e-05 6.466e-05 7.144e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.992e-05 6.385e-05 6.661e-05 0.00% Castro::derive() 18 4.185e-05 4.584e-05 5.775e-05 0.00% MLMG::computeMLResidual() 18 4.05e-05 4.575e-05 5.181e-05 0.00% Amr::regrid() 1 1.566e-05 2.014e-05 5.172e-05 0.00% Gravity::actual_multilevel_solve() 1 3.179e-05 3.449e-05 3.681e-05 0.00% Amr::InitializeInit() 1 8.369e-06 1.045e-05 3.436e-05 0.00% Amr::timeStep() 1 1.407e-05 1.842e-05 2.935e-05 0.00% Castro::advance() 1 7.224e-06 1.068e-05 2.182e-05 0.00% Castro::post_regrid() 2 1.132e-05 1.402e-05 2.161e-05 0.00% Castro::post_init() 1 1.184e-05 1.307e-05 1.479e-05 0.00% Amr::bldFineLevels() 1 7.077e-06 7.863e-06 1.153e-05 0.00% Amr::init() 1 5.848e-06 7.149e-06 1.037e-05 0.00% Castro::post_timestep() 1 5.963e-06 6.37e-06 6.892e-06 0.00% Amr::initialInit() 1 2.797e-06 3.283e-06 4.995e-06 0.00% Other 2502 0.05246 0.06689 0.1017 0.69% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.62 14.63 14.63 100.00% Amr::coarseTimeStep() 1 12.37 12.37 12.37 84.50% Amr::timeStep() 1 12.31 12.31 12.31 84.12% Castro::advance() 1 12.24 12.25 12.25 83.71% Castro::subcycle_advance_ctu() 1 12.19 12.2 12.2 83.36% Castro::do_advance_ctu() 20 12.15 12.16 12.16 83.10% Castro::construct_ctu_hydro_source() 20 6.671 6.704 6.748 46.12% Castro::clean_state() 79 3.672 3.795 3.957 27.04% Castro::computeTemp() 80 2.351 2.471 2.635 18.00% Gravity::solve_phi_with_mlmg() 19 1.6 1.654 1.677 11.46% Castro::construct_new_gravity() 17 1.517 1.571 1.594 10.89% Gravity::solve_for_phi() 17 1.48 1.534 1.557 10.64% Castro::reset_internal_energy(MultiFab) 80 1.506 1.511 1.521 10.39% Castro::initialize_do_advance() 20 1.403 1.438 1.464 10.00% Gravity::actual_solve_with_mlmg() 19 1.399 1.403 1.407 9.62% Amr::InitAmr() 1 1.353 1.356 1.362 9.31% MLMG::solve() 18 1.342 1.348 1.349 9.22% MLMG::oneIter() 128 1.209 1.211 1.215 8.30% Castro::do_old_sources() 20 1.158 1.186 1.206 8.24% MLMG::mgVcycle() 128 1.194 1.196 1.2 8.20% Castro::do_new_sources() 17 0.9246 0.9753 1.018 6.96% FabArray::ParallelCopy() 1429 0.5114 0.8106 0.9998 6.83% Castro::evaluate_source_change() 196 0.9607 0.965 0.9773 6.68% Amr::init() 1 0.9048 0.9048 0.9048 6.18% Amr::initialInit() 1 0.8461 0.8462 0.8463 5.78% MLCellLinOp::smooth() 1544 0.6238 0.6797 0.8229 5.62% amrex::Copy() 724 0.8002 0.8046 0.8167 5.58% Castro::enforce_min_density() 79 0.7981 0.8012 0.8089 5.53% FabArray::ParallelCopy_finish() 1429 0.2691 0.5673 0.7599 5.19% AmrLevel::FillPatch() 147 0.5502 0.6274 0.7034 4.81% Amr::InitializeInit() 1 0.445 0.4751 0.6088 4.16% Amr::defBaseLevel() 1 0.445 0.4751 0.6088 4.16% MultiFab::AddProduct() 3439 0.552 0.554 0.5604 3.83% Castro::normalize_species() 79 0.5468 0.5498 0.5549 3.79% Castro::reset_internal_energy(Fab) 160 0.5223 0.5255 0.5325 3.64% FillPatchIterator::Initialize 147 0.3669 0.4439 0.5208 3.56% FillPatchIterator::FillFromLevel0() 147 0.3392 0.423 0.4929 3.37% FillPatchSingleLevel 147 0.3385 0.4222 0.4922 3.36% MLPoisson::Fsmooth() 6144 0.356 0.3822 0.4599 3.14% MLMG::mgVcycle_up::1 128 0.1067 0.2824 0.4586 3.13% MLCellLinOp::applyBC() 7148 0.3214 0.3556 0.44 3.01% Amr::FinalizeInit() 1 0.2374 0.3711 0.4011 2.74% MLMG::addInterpCorrection() 768 0.0665 0.3217 0.3976 2.72% Castro::Castro() 1 0.3651 0.3698 0.3829 2.62% Castro::initMFs() 1 0.3591 0.3639 0.3769 2.58% FabArray::setVal() 5013 0.3488 0.3522 0.3575 2.44% MultiFab::Subtract() 159 0.3377 0.3383 0.3405 2.33% FabArray::FillBoundary() 6380 0.2474 0.2754 0.3154 2.16% MLMG::mgVcycle_up::2 128 0.01385 0.1266 0.2978 2.04% MLMG::mgVcycle_down::0 128 0.2744 0.2769 0.2798 1.91% Gravity::fill_multipole_BCs() 19 0.1995 0.2497 0.2714 1.85% Castro::expand_state() 20 0.212 0.2359 0.2678 1.83% FabArray::ParallelCopy_nowait() 1429 0.2232 0.2395 0.2656 1.82% Castro::initData() 1 0.07817 0.1034 0.2383 1.63% check_for_negative_density() 20 0.1024 0.1652 0.227 1.55% MLMG::mgVcycle_up::0 128 0.2136 0.2143 0.2152 1.47% Amr::bldFineLevels() 1 0.04605 0.1797 0.2091 1.43% MLCellLinOp::apply() 986 0.1582 0.1684 0.2059 1.41% FabArray::sum() 3439 0.1974 0.1985 0.2014 1.38% MLMG::mgVcycle_up::3 128 0.0138 0.02764 0.1916 1.31% FillBoundary_nowait() 6380 0.103 0.1323 0.1686 1.15% Amr::grid_places() 2 0.002503 0.1362 0.1662 1.14% AmrMesh::MakeNewGrids() 2 0.00245 0.1361 0.1662 1.14% TagBoxArray::mapPRD 2 0.0008348 0.1085 0.1641 1.12% FillBoundary_finish() 6380 0.1013 0.1286 0.1627 1.11% MLMG:computeResOfCorrection() 768 0.119 0.1286 0.1529 1.04% MLCellLinOp::correctionResidual() 768 0.1175 0.1271 0.1513 1.03% Castro::post_init() 1 0.1437 0.1482 0.149 1.02% TagBoxArray::collate() 2 0.0001979 0.0263 0.1457 1.00% MLMG::mgVcycle_down::1 128 0.1266 0.1297 0.1341 0.92% Castro::construct_new_source() 102 0.1293 0.1304 0.131 0.90% MLPoisson::Fapply() 986 0.09492 0.0993 0.1141 0.78% Gravity::multilevel_solve_for_new_phi() 1 0.1082 0.1082 0.1082 0.74% Gravity::actual_multilevel_solve() 1 0.1081 0.1081 0.1082 0.74% Amr::writePlotFile() 2 0.1016 0.1017 0.1017 0.70% MLCellLinOp::solutionResidual() 147 0.08898 0.09577 0.1004 0.69% Castro::create_source_corrector() 17 0.05505 0.06978 0.09276 0.63% Castro::post_regrid() 2 0.08585 0.08667 0.09109 0.62% Castro::construct_old_source() 120 0.08293 0.08355 0.08491 0.58% MLMG::computeResidual() 128 0.0772 0.07817 0.08138 0.56% MLMG::mgVcycle_down::2 128 0.01603 0.04627 0.0794 0.54% Castro::construct_new_gravity_source() 17 0.07294 0.07368 0.07445 0.51% FabArray::Saxpy() 91 0.07198 0.07239 0.07383 0.50% MLMG::mgVcycle_bottom 128 0.0003563 0.004886 0.07263 0.50% MLMG::actualBottomSolve() 128 0.0001337 0.004653 0.07239 0.49% Castro::derive() 18 0.07056 0.07127 0.07217 0.49% AmrLevel::derive() 18 0.07051 0.07123 0.07212 0.49% MLCGSolver::bicgstab 8 0 0.004409 0.07055 0.48% Castro::construct_old_gravity() 20 0.04747 0.06012 0.07044 0.48% Gravity::get_old_grav_vector() 20 0.04735 0.06 0.07031 0.48% MLMG::mgVcycle_down::3 128 0.01536 0.02226 0.06824 0.47% Castro::post_timestep() 1 0.06019 0.06264 0.06571 0.45% FabArray::Xpay() 915 0.05165 0.05366 0.05925 0.40% amrex::average_down 768 0.03424 0.04114 0.05784 0.40% Castro::initialize_advance() 1 0.04956 0.0505 0.05407 0.37% Castro::construct_old_gravity_source() 20 0.04892 0.04959 0.05065 0.35% Amr::regrid() 1 0.04406 0.04467 0.04677 0.32% Castro::check_for_nan() 37 0.04638 0.04656 0.04676 0.32% MultiFab::contains_nan() 37 0.04625 0.04644 0.04664 0.32% MLPoisson::define() 19 0.03274 0.03777 0.04405 0.30% Castro::retry_advance_ctu() 11 0.03622 0.03753 0.0432 0.30% MLMG::mgVcycle_down::4 128 0.01529 0.01716 0.04009 0.27% Gravity::get_new_grav_vector() 18 0.03628 0.03653 0.03681 0.25% StateDataPhysBCFunct::() 147 0.0007854 0.01842 0.03471 0.24% StateData::FillBoundary(geom) 143 0 0.01753 0.03375 0.23% MLMG::mgVcycle_up::4 128 0.01378 0.01516 0.03162 0.22% Castro::apply_source_to_state() 37 0.02932 0.02954 0.03087 0.21% VisMF::Write(FabArray) 4 0.01578 0.01735 0.02865 0.20% MLMG::mgVcycle_down::5 128 0.01527 0.01635 0.02796 0.19% MLCellLinOp::defineAuxData() 19 0.02004 0.02173 0.02713 0.19% Castro::save_data_for_retry() 3 0.01933 0.02007 0.02544 0.17% MLMG::prepareForSolve() 18 0.0206 0.02164 0.0243 0.17% MLMG::mgVcycle_up::5 128 0.01392 0.01479 0.02399 0.16% Castro::sum_integrated_quantities() 2 0.02304 0.0232 0.02321 0.16% MLMG::computeMLResidual() 18 0.01119 0.01706 0.02197 0.15% MLMG::ResNormInf() 146 0.01219 0.01479 0.01757 0.12% Castro::problem_post_timestep() 2 0.0101 0.01318 0.01646 0.11% Castro::wd_update() 2 0.01008 0.01315 0.01645 0.11% Gravity::test_composite_phi() 1 0.01529 0.01532 0.01534 0.10% FabArrayBase::getCPC() 2005 0.01008 0.01153 0.01532 0.10% FabArrayBase::getFB() 6380 0.01096 0.01179 0.01503 0.10% amrex::Add() 128 0.01451 0.01455 0.01468 0.10% FabArray::setDomainBndry() 147 0.0005174 0.007602 0.01444 0.10% FabArray::norminf() 244 0.0107 0.01107 0.01404 0.10% MLLinOp::define() 19 0.005864 0.009418 0.0132 0.09% MLLinOp::defineGrids() 19 0.005777 0.009334 0.01311 0.09% Other 2502 0.08487 0.1038 0.1444 0.99% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Pinned Memory Usage: --------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max --------------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 16 9 B 11 B 19 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- Cpu Memory Usage: ------------------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max ------------------------------------------------------------------------------------------------------------------------------- Castro::save_data_for_retry() 416 416 1569 KiB 1571 KiB 1572 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 410 KiB 413 KiB 416 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1040 KiB 1042 KiB 1042 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 990 KiB 991 KiB 991 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 664 KiB 665 KiB 665 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 4312 B 4429 B 4948 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 50 KiB 50 KiB 50 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 12 KiB 13 KiB 14 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 416 KiB 416 KiB 416 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 2335 B 4072 B 5300 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1338 B 1353 B 1370 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 17 KiB 17 KiB 17 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 132 KiB 132 KiB 132 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 192 128 KiB 128 KiB 129 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 6102 9622 B 9 KiB 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 258 B 258 B 260 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 435 B 437 B 440 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 26 B 27 B 43 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 33 KiB 33 KiB 33 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 55 B 55 B 55 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 3 B 8 B 9 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 111 B 112 B 113 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 147 B 187 B 220 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 64 31 KiB 31 KiB 31 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 2432 229 B 235 B 240 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1071 B 1075 B 1088 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 121 B 121 B 121 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1657 B 1716 B 1743 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 32 0 B 0 B 1 B 16 KiB 16 KiB 16 KiB MLCellLinOp::defineAuxData() 11704 11704 435 B 486 B 675 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 452 B 454 B 455 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 24 B 33 B 42 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 19 B 29 B 2916 B 2916 B 2916 B OwnerMask() 64 64 0 B 0 B 0 B 2592 B 2592 B 2592 B MLMG::addInterpCorrection() 5632 5632 1 B 16 B 25 B 1024 B 1024 B 1024 B amrex::average_down 5632 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 64 0 B 5 B 7 B 648 B 648 B 648 B MLCGSolver::bicgstab 640 640 0 B 0 B 2 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 114375 0 B 0 B 0 B 8 B 8 B 8 B ------------------------------------------------------------------------------------------------------------------------------- AMReX (25.06-22-g1d307c2c661a-dirty) finalized