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:49:07 UTC on 2025-07-21. Successfully read inputs file ... Castro git describe: 25.07-11-gcb51e8cc9 AMReX git describe: 25.06-22-g1d307c2c6-dirty Microphysics git describe: joss_submit 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.178673899 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001167534 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010666949 Gravity::multilevel_solve_for_new_phi() time = 0.113296077 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.011504897 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.052715324 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.007662239 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.005155411 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.064652982 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.004722078 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.064728853 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.006468449 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.060622277 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.005018037 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.060601602 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.34028681 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.0144044 Gravity::solve_for_phi() time = 0.092849236 on level 0 Castro::construct_new_gravity() time = 0.095163804 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.061675839 on level 0 Castro::construct_old_gravity() time = 0.002839603 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.05991205 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.340604534 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.017482144 Gravity::solve_for_phi() time = 0.096450104 on level 0 Castro::construct_new_gravity() time = 0.098554015 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.060784161 on level 0 Castro::construct_old_gravity() time = 0.004309357 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.059272369 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.340487626 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.014174508 Gravity::solve_for_phi() time = 0.090281625 on level 0 Castro::construct_new_gravity() time = 0.092455334 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.05955717 on level 0 Castro::construct_old_gravity() time = 0.003043893 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.059339581 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.346605804 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.014414244 Gravity::solve_for_phi() time = 0.091660658 on level 0 Castro::construct_new_gravity() time = 0.093788404 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.060020258 on level 0 Castro::construct_old_gravity() time = 0.004295626 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.059548608 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.348120709 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.014809726 Gravity::solve_for_phi() time = 0.091656593 on level 0 Castro::construct_new_gravity() time = 0.093825771 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.059503324 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.003912732 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.05945725 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.342291909 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.014795925 Gravity::solve_for_phi() time = 0.091788912 on level 0 Castro::construct_new_gravity() time = 0.093902567 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.059626616 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.004069816 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.059623808 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.349521193 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.014888022 Gravity::solve_for_phi() time = 0.091078525 on level 0 Castro::construct_new_gravity() time = 0.093278538 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.060156558 on level 0 Castro::construct_old_gravity() time = 0.00274672 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.059326904 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.34418648 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.014608077 Gravity::solve_for_phi() time = 0.091550875 on level 0 Castro::construct_new_gravity() time = 0.09366204 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.059615672 on level 0 Castro::construct_old_gravity() time = 0.003009672 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.059394897 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.348652725 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.014787304 Gravity::solve_for_phi() time = 0.092350707 on level 0 Castro::construct_new_gravity() time = 0.094555706 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.060040972 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.00279477 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.059506622 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.343717911 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.014930579 Gravity::solve_for_phi() time = 0.092255866 on level 0 Castro::construct_new_gravity() time = 0.094458902 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.061317902 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.003050838 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.061333248 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.356923691 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.015522589 Gravity::solve_for_phi() time = 0.091786349 on level 0 Castro::construct_new_gravity() time = 0.093914452 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.060361787 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.00359841 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.058875925 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.34683517 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.01488708 Gravity::solve_for_phi() time = 0.092304797 on level 0 Castro::construct_new_gravity() time = 0.09443046 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.060246554 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.003677404 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.065500548 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.352112057 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.015714522 Gravity::solve_for_phi() time = 0.091454829 on level 0 Castro::construct_new_gravity() time = 0.093671465 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.060058676 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.004618965 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.059445229 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.34416631 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.015126936 Gravity::solve_for_phi() time = 0.09132582 on level 0 Castro::construct_new_gravity() time = 0.093445933 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.059849902 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.003075969 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.059233108 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.348938916 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.01513883 Gravity::solve_for_phi() time = 0.091649761 on level 0 Castro::construct_new_gravity() time = 0.093783515 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.060061279 on level 0 Castro::construct_old_gravity() time = 0.00277969 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.059292299 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.344001276 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.015306 Gravity::solve_for_phi() time = 0.09274329 on level 0 Castro::construct_new_gravity() time = 0.094848513 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.059781891 on level 0 Castro::construct_old_gravity() time = 0.002870302 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.059145477 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.354780687 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002679866097 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.021061427 Gravity::solve_for_phi() time = 0.097306678 on level 0 Castro::construct_new_gravity() time = 0.09942116 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.059980244 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.2877532 [STEP 1] FAB kilobyte spread across MPI nodes: [6998 ... 6998] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.050951451 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.006570896 seconds Ending run at 06:49:22 UTC on 2025-07-21. Run time = 14.57860387 Run time without initialization = 12.34663145 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.5700 ... 14.5727 ... 14.5787 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.7187 5.7562 5.8382 40.05% Amr::InitAmr() 1 1.296 1.299 1.305 8.95% Castro::computeTemp() 80 0.8444 0.9577 1.113 7.64% amrex::Copy() 724 0.7997 0.8048 0.8296 5.69% FabArray::ParallelCopy_finish() 1429 0.1915 0.5513 0.7299 5.01% Castro::normalize_species() 79 0.5469 0.556 0.6142 4.21% MultiFab::AddProduct() 3439 0.5528 0.5549 0.5627 3.86% Castro::reset_internal_energy(Fab) 160 0.5194 0.5251 0.5317 3.65% MLPoisson::Fsmooth() 6144 0.3542 0.3826 0.4689 3.22% Castro::initMFs() 1 0.3594 0.3676 0.3852 2.64% FabArray::setVal() 5013 0.3488 0.3523 0.3648 2.50% MultiFab::Subtract() 159 0.3376 0.3388 0.3431 2.35% Gravity::fill_multipole_BCs() 19 0.2044 0.2565 0.2784 1.91% FabArray::ParallelCopy_nowait() 1429 0.1985 0.2185 0.2386 1.64% check_for_negative_density() 20 0.05228 0.1602 0.2355 1.62% FabArray::sum() 3439 0.1978 0.1989 0.2013 1.38% Castro::initData() 1 0.003355 0.0283 0.1669 1.14% TagBoxArray::collate() 2 0.000203 0.03237 0.1576 1.08% FillBoundary_nowait() 6380 0.09069 0.1202 0.1547 1.06% FillBoundary_finish() 6380 0.1049 0.1245 0.1446 0.99% MLCellLinOp::applyBC() 7148 0.06371 0.07988 0.1292 0.89% MLPoisson::Fapply() 986 0.09475 0.09907 0.1174 0.81% Castro::construct_new_gravity_source() 17 0.07258 0.07346 0.07472 0.51% FabArray::Saxpy() 91 0.07201 0.07251 0.07402 0.51% FabArray::Xpay() 915 0.05166 0.05403 0.06026 0.41% Castro::construct_old_gravity_source() 20 0.0494 0.04974 0.0519 0.36% MultiFab::contains_nan() 37 0.04626 0.04655 0.0471 0.32% MLMG::addInterpCorrection() 768 0.03379 0.03618 0.04456 0.31% AmrLevel::derive() 18 0.04121 0.04165 0.04323 0.30% amrex::average_down 768 0.02515 0.02723 0.03445 0.24% StateData::FillBoundary(geom) 143 0 0.01761 0.03372 0.23% VisMF::Write(FabArray) 4 0.0166 0.0187 0.03099 0.21% Castro::reset_internal_energy(MultiFab) 80 0.02322 0.02386 0.02491 0.17% MLCellLinOp::defineAuxData() 19 0.01556 0.01659 0.02075 0.14% Amr::writePlotFile() 2 0.008324 0.01898 0.02034 0.14% Castro::enforce_min_density() 79 0.0188 0.01917 0.0197 0.14% FabArray::FillBoundary() 6380 0.01391 0.01446 0.01776 0.12% MLCGSolver::bicgstab 8 0 0.00105 0.0168 0.12% Castro::construct_new_source() 102 0.01616 0.01623 0.01646 0.11% Castro::evaluate_source_change() 196 0.01396 0.01431 0.01543 0.11% FabArray::setDomainBndry() 147 0.0004974 0.007621 0.01493 0.10% amrex::Add() 128 0.0145 0.01454 0.01471 0.10% Castro::wd_update() 2 0.005947 0.01047 0.01421 0.10% FabArray::norminf() 244 0.01075 0.01114 0.01412 0.10% FillPatchIterator::Initialize 147 0.01298 0.01335 0.01362 0.09% FabArrayBase::getFB() 6380 0.01062 0.01093 0.01352 0.09% MLCellLinOp::smooth() 1544 0.01019 0.01062 0.0123 0.08% MLMG::prepareForSolve() 18 0.01014 0.01056 0.01196 0.08% Castro::construct_old_source() 120 0.01026 0.01039 0.01116 0.08% MLMG::solve() 18 0.0008597 0.005114 0.009823 0.07% Amr::writeSmallPlotFile() 2 0.00674 0.009096 0.009438 0.06% Gravity::actual_solve_with_mlmg() 19 0.007215 0.007657 0.009275 0.06% FabArrayBase::getCPC() 2005 0.007717 0.007912 0.008511 0.06% MLMG::ResNormInf() 146 0.004515 0.006128 0.007827 0.05% Gravity::get_old_grav_vector() 20 0.006497 0.006673 0.007211 0.05% Gravity::get_new_grav_vector() 18 0.005987 0.006136 0.006579 0.05% Castro::initialize_do_advance() 20 0.004449 0.004964 0.005421 0.04% MLCellLinOp::apply() 986 0.002044 0.002275 0.004973 0.03% FabArray::ParallelCopy() 1429 0.003609 0.003763 0.004166 0.03% Castro::sum_integrated_quantities() 2 0.001581 0.002424 0.002693 0.02% AmrLevel::FillPatch() 147 0.002323 0.002481 0.002637 0.02% MLCellLinOp::correctionResidual() 768 0.001911 0.001957 0.002131 0.01% MLMG::mgVcycle() 128 0.00165 0.001721 0.00199 0.01% MLMG:computeResOfCorrection() 768 0.001511 0.001555 0.001664 0.01% Castro::save_data_for_retry() 3 0.001307 0.00137 0.001461 0.01% MLMG::actualBottomSolve() 128 0.0001315 0.000207 0.001269 0.01% Castro::subcycle_advance_ctu() 1 0.0002352 0.0003217 0.001139 0.01% StateDataPhysBCFunct::() 147 0.0007154 0.0008797 0.0009903 0.01% main() 1 0.0007116 0.0008109 0.0009697 0.01% Gravity::solve_for_phi() 17 0.0006588 0.0007013 0.0007733 0.01% FillPatchIterator::FillFromLevel0() 147 0.0007018 0.0007221 0.0007645 0.01% MLMG::mgVcycle_down::3 128 0.0005107 0.0005419 0.0006272 0.00% MLMG::mgVcycle_down::2 128 0.0005436 0.0005773 0.0006205 0.00% MLMG::mgVcycle_down::0 128 0.0005559 0.0005755 0.0006106 0.00% MLMG::mgVcycle_down::4 128 0.0004852 0.0005054 0.000602 0.00% MLMG::mgVcycle_down::5 128 0.0005098 0.0005234 0.0005946 0.00% MLMG::mgVcycle_down::1 128 0.0005646 0.0005714 0.0005932 0.00% AmrMesh::MakeNewGrids() 2 0.0005446 0.000559 0.0005791 0.00% Castro::Castro() 1 0.0004746 0.0005165 0.0005713 0.00% Castro::clean_state() 79 0.0004339 0.0004647 0.0005366 0.00% FillPatchSingleLevel 147 0.0004759 0.000485 0.0005014 0.00% MLCellLinOp::solutionResidual() 147 0.0004478 0.0004598 0.0004902 0.00% Castro::do_advance_ctu() 20 0.0003783 0.0004299 0.0004892 0.00% MLMG::oneIter() 128 0.000393 0.000406 0.0004321 0.00% MLMG::mgVcycle_up::0 128 0.0003715 0.0003797 0.000401 0.00% MLMG::mgVcycle_up::1 128 0.0003645 0.0003765 0.0003948 0.00% MLMG::mgVcycle_up::2 128 0.0003239 0.0003475 0.0003846 0.00% MLMG::mgVcycle_up::3 128 0.0003211 0.0003341 0.000373 0.00% Castro::construct_new_gravity() 17 0.0001585 0.0001832 0.0003728 0.00% MLMG::mgVcycle_up::4 128 0.0003226 0.0003307 0.0003705 0.00% MLMG::mgVcycle_up::5 128 0.00031 0.0003188 0.0003645 0.00% Castro::do_old_sources() 20 0.0002772 0.0002962 0.000342 0.00% Castro::problem_post_timestep() 2 1.147e-05 3.028e-05 0.0002844 0.00% MLMG::computeResidual() 128 0.000237 0.0002437 0.0002558 0.00% MLMG::mgVcycle_bottom 128 0.0002215 0.0002296 0.0002521 0.00% Gravity::solve_phi_with_mlmg() 19 0.000216 0.0002336 0.0002491 0.00% Castro::do_new_sources() 17 0.0002088 0.0002223 0.000245 0.00% Amr::coarseTimeStep() 1 0.0001467 0.0001752 0.000184 0.00% Castro::initialize_advance() 1 0.0001304 0.0001476 0.0001819 0.00% MLPoisson::define() 19 0.0001287 0.0001407 0.0001562 0.00% TagBoxArray::mapPRD 2 0.0001229 0.0001302 0.000146 0.00% Castro::construct_old_gravity() 20 0.0001128 0.0001223 0.0001338 0.00% Castro::retry_advance_ctu() 11 6.215e-05 7.549e-05 0.0001333 0.00% Castro::check_for_nan() 37 0.000104 0.0001156 0.000133 0.00% Gravity::test_composite_phi() 1 0.0001046 0.0001101 0.0001283 0.00% Amr::defBaseLevel() 1 9.08e-05 0.0001096 0.0001233 0.00% Castro::apply_source_to_state() 37 9.726e-05 0.0001061 0.0001186 0.00% MLLinOp::define() 19 8.097e-05 8.746e-05 0.0001005 0.00% Castro::create_source_corrector() 17 5.535e-05 6.509e-05 7.38e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.947e-05 6.362e-05 6.729e-05 0.00% Amr::FinalizeInit() 1 1.375e-05 1.937e-05 6.543e-05 0.00% Amr::grid_places() 2 5.284e-05 5.727e-05 6.366e-05 0.00% Castro::expand_state() 20 5.066e-05 5.597e-05 6.144e-05 0.00% MLMG::computeMLResidual() 18 3.993e-05 4.593e-05 5.224e-05 0.00% Castro::derive() 18 4.357e-05 4.616e-05 5.005e-05 0.00% Amr::regrid() 1 1.564e-05 2.052e-05 4.749e-05 0.00% Gravity::actual_multilevel_solve() 1 3.186e-05 3.391e-05 3.618e-05 0.00% Amr::timeStep() 1 1.362e-05 1.843e-05 3.259e-05 0.00% Amr::init() 1 6.216e-06 8.097e-06 2.156e-05 0.00% Amr::InitializeInit() 1 8.425e-06 9.847e-06 1.572e-05 0.00% Castro::post_regrid() 2 1.122e-05 1.387e-05 1.544e-05 0.00% Castro::post_init() 1 1.202e-05 1.274e-05 1.39e-05 0.00% Castro::advance() 1 7.314e-06 9.169e-06 1.388e-05 0.00% Castro::post_timestep() 1 5.821e-06 7.364e-06 1.248e-05 0.00% Amr::bldFineLevels() 1 6.921e-06 7.706e-06 8.591e-06 0.00% Amr::initialInit() 1 2.951e-06 3.189e-06 3.761e-06 0.00% Other 2519 0.04647 0.05964 0.09586 0.66% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.57 14.57 14.58 100.00% Amr::coarseTimeStep() 1 12.35 12.35 12.35 84.68% Amr::timeStep() 1 12.29 12.29 12.29 84.29% Castro::advance() 1 12.22 12.22 12.23 83.87% Castro::subcycle_advance_ctu() 1 12.17 12.17 12.18 83.53% Castro::do_advance_ctu() 20 12.13 12.14 12.14 83.27% Castro::construct_ctu_hydro_source() 20 6.672 6.698 6.74 46.23% Castro::clean_state() 79 3.673 3.802 3.944 27.05% Castro::computeTemp() 80 2.352 2.471 2.622 17.98% Gravity::solve_phi_with_mlmg() 19 1.615 1.668 1.69 11.59% Castro::construct_new_gravity() 17 1.527 1.58 1.602 10.99% Gravity::solve_for_phi() 17 1.489 1.543 1.565 10.73% Castro::reset_internal_energy(MultiFab) 80 1.502 1.513 1.533 10.52% Castro::initialize_do_advance() 20 1.381 1.425 1.456 9.99% Gravity::actual_solve_with_mlmg() 19 1.403 1.408 1.412 9.68% MLMG::solve() 18 1.348 1.354 1.356 9.30% Amr::InitAmr() 1 1.296 1.299 1.305 8.95% MLMG::oneIter() 128 1.213 1.216 1.217 8.35% MLMG::mgVcycle() 128 1.198 1.201 1.202 8.25% Castro::do_old_sources() 20 1.145 1.178 1.196 8.21% Castro::do_new_sources() 17 0.9255 0.9749 1.018 6.98% Castro::evaluate_source_change() 196 0.9627 0.9665 0.9811 6.73% FabArray::ParallelCopy() 1429 0.441 0.7939 0.974 6.68% Amr::init() 1 0.9274 0.9274 0.9274 6.36% Amr::initialInit() 1 0.8665 0.8665 0.8666 5.94% amrex::Copy() 724 0.7997 0.8048 0.8296 5.69% MLCellLinOp::smooth() 1544 0.6212 0.6759 0.8229 5.64% Castro::enforce_min_density() 79 0.7987 0.8022 0.8205 5.63% FabArray::ParallelCopy_finish() 1429 0.1915 0.5513 0.7299 5.01% AmrLevel::FillPatch() 147 0.519 0.5985 0.6811 4.67% Amr::InitializeInit() 1 0.4458 0.4804 0.6231 4.27% Amr::defBaseLevel() 1 0.4458 0.4804 0.6231 4.27% Castro::normalize_species() 79 0.5469 0.556 0.6142 4.21% MultiFab::AddProduct() 3439 0.5528 0.5549 0.5627 3.86% Castro::reset_internal_energy(Fab) 160 0.5194 0.5251 0.5317 3.65% FillPatchIterator::Initialize 147 0.3304 0.415 0.4986 3.42% FillPatchIterator::FillFromLevel0() 147 0.3018 0.394 0.4756 3.26% FillPatchSingleLevel 147 0.3011 0.3933 0.4749 3.26% MLPoisson::Fsmooth() 6144 0.3542 0.3826 0.4689 3.22% MLMG::mgVcycle_up::1 128 0.1053 0.2864 0.4667 3.20% MLCellLinOp::applyBC() 7148 0.3115 0.3507 0.4281 2.94% Amr::FinalizeInit() 1 0.2435 0.3861 0.4207 2.89% MLMG::addInterpCorrection() 768 0.06699 0.3307 0.4069 2.79% Castro::Castro() 1 0.3654 0.3738 0.3914 2.68% Castro::initMFs() 1 0.3596 0.3678 0.3853 2.64% FabArray::setVal() 5013 0.3488 0.3523 0.3648 2.50% MultiFab::Subtract() 159 0.3376 0.3388 0.3431 2.35% FabArray::FillBoundary() 6380 0.2386 0.2708 0.3097 2.12% MLMG::mgVcycle_up::2 128 0.01382 0.1296 0.3058 2.10% Gravity::fill_multipole_BCs() 19 0.2066 0.2586 0.2806 1.92% MLMG::mgVcycle_down::0 128 0.2763 0.2778 0.2796 1.92% FabArray::ParallelCopy_nowait() 1429 0.2182 0.2388 0.26 1.78% Castro::expand_state() 20 0.192 0.2235 0.2521 1.73% Castro::initData() 1 0.07858 0.1048 0.2443 1.68% check_for_negative_density() 20 0.05228 0.1602 0.2355 1.62% Amr::bldFineLevels() 1 0.04712 0.1894 0.2229 1.53% MLMG::mgVcycle_up::0 128 0.209 0.2098 0.2112 1.45% MLCellLinOp::apply() 986 0.1557 0.1673 0.207 1.42% FabArray::sum() 3439 0.1978 0.1989 0.2013 1.38% MLMG::mgVcycle_up::3 128 0.0138 0.02796 0.1978 1.36% Amr::grid_places() 2 0.002547 0.1452 0.1798 1.23% AmrMesh::MakeNewGrids() 2 0.002493 0.1452 0.1798 1.23% TagBoxArray::mapPRD 2 0.0008226 0.1115 0.1756 1.20% FillBoundary_nowait() 6380 0.1016 0.1319 0.167 1.15% TagBoxArray::collate() 2 0.000203 0.03237 0.1576 1.08% Castro::post_init() 1 0.1495 0.1529 0.1549 1.06% MLMG:computeResOfCorrection() 768 0.1172 0.1277 0.154 1.06% MLCellLinOp::correctionResidual() 768 0.1157 0.1262 0.1524 1.05% FillBoundary_finish() 6380 0.1049 0.1245 0.1446 0.99% MLMG::mgVcycle_down::1 128 0.1249 0.1297 0.1335 0.92% Castro::construct_new_source() 102 0.1292 0.1301 0.1318 0.90% MLPoisson::Fapply() 986 0.09475 0.09907 0.1174 0.81% Gravity::multilevel_solve_for_new_phi() 1 0.1132 0.1132 0.1133 0.78% Gravity::actual_multilevel_solve() 1 0.1131 0.1132 0.1132 0.78% Amr::writePlotFile() 2 0.104 0.1041 0.1041 0.71% MLCellLinOp::solutionResidual() 147 0.08947 0.09589 0.1028 0.71% Castro::post_regrid() 2 0.08591 0.08791 0.09134 0.63% Castro::construct_old_source() 120 0.08336 0.0839 0.0864 0.59% MLMG::mgVcycle_down::2 128 0.01607 0.0475 0.08141 0.56% Castro::create_source_corrector() 17 0.05207 0.06722 0.08128 0.56% MLMG::computeResidual() 128 0.07717 0.07829 0.07913 0.54% Castro::construct_new_gravity_source() 17 0.07258 0.07346 0.07472 0.51% Castro::derive() 18 0.07074 0.0714 0.07451 0.51% AmrLevel::derive() 18 0.07069 0.07135 0.07446 0.51% FabArray::Saxpy() 91 0.07201 0.07251 0.07402 0.51% MLMG::mgVcycle_bottom 128 0.000353 0.004938 0.07355 0.50% MLMG::actualBottomSolve() 128 0.0001315 0.004709 0.07329 0.50% MLCGSolver::bicgstab 8 0 0.00446 0.07136 0.49% Castro::construct_old_gravity() 20 0.04393 0.05856 0.06977 0.48% Gravity::get_old_grav_vector() 20 0.0438 0.05844 0.06965 0.48% MLMG::mgVcycle_down::3 128 0.01542 0.02214 0.06847 0.47% Castro::post_timestep() 1 0.06021 0.06324 0.06633 0.45% FabArray::Xpay() 915 0.05166 0.05403 0.06026 0.41% amrex::average_down 768 0.03429 0.04114 0.05888 0.40% Castro::initialize_advance() 1 0.04961 0.05027 0.05261 0.36% Castro::construct_old_gravity_source() 20 0.0494 0.04974 0.0519 0.36% Castro::check_for_nan() 37 0.04637 0.04667 0.04724 0.32% MultiFab::contains_nan() 37 0.04626 0.04655 0.0471 0.32% Amr::regrid() 1 0.04426 0.04526 0.04698 0.32% MLPoisson::define() 19 0.03037 0.03579 0.04131 0.28% MLMG::mgVcycle_down::4 128 0.0154 0.01713 0.04124 0.28% Castro::retry_advance_ctu() 11 0.03595 0.03677 0.03795 0.26% Gravity::get_new_grav_vector() 18 0.03606 0.0366 0.03775 0.26% StateDataPhysBCFunct::() 147 0.0007154 0.01849 0.03471 0.24% StateData::FillBoundary(geom) 143 0 0.01761 0.03372 0.23% MLMG::mgVcycle_up::4 128 0.01378 0.01505 0.03264 0.22% VisMF::Write(FabArray) 4 0.0166 0.0187 0.03099 0.21% Castro::apply_source_to_state() 37 0.02932 0.02958 0.02989 0.21% MLMG::mgVcycle_down::5 128 0.01532 0.0162 0.02839 0.19% MLCellLinOp::defineAuxData() 19 0.0197 0.02159 0.0279 0.19% MLMG::prepareForSolve() 18 0.02044 0.02147 0.02497 0.17% MLMG::mgVcycle_up::5 128 0.01388 0.01465 0.02489 0.17% MLMG::computeMLResidual() 18 0.01114 0.01708 0.02393 0.16% Castro::sum_integrated_quantities() 2 0.02209 0.02225 0.02227 0.15% Castro::save_data_for_retry() 3 0.01914 0.01957 0.02061 0.14% Castro::problem_post_timestep() 2 0.01078 0.01423 0.01758 0.12% Castro::wd_update() 2 0.01076 0.0142 0.01757 0.12% MLMG::ResNormInf() 146 0.01407 0.01559 0.01729 0.12% FabArrayBase::getCPC() 2005 0.01 0.01147 0.01563 0.11% Gravity::test_composite_phi() 1 0.01524 0.01529 0.01532 0.11% FabArrayBase::getFB() 6380 0.01092 0.01168 0.01519 0.10% FabArray::setDomainBndry() 147 0.0004974 0.007621 0.01493 0.10% Amr::writeSmallPlotFile() 2 0.01472 0.01472 0.01472 0.10% amrex::Add() 128 0.0145 0.01454 0.01471 0.10% FabArray::norminf() 244 0.01075 0.01114 0.01412 0.10% MLLinOp::define() 19 0.004126 0.007616 0.01129 0.08% Other 2519 0.07567 0.09712 0.1443 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 12 B 16 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 1573 KiB 1575 KiB 1575 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 412 KiB 414 KiB 420 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1045 KiB 1046 KiB 1046 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 994 KiB 995 KiB 995 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 666 KiB 667 KiB 667 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 4395 B 4476 B 5082 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 50 KiB 50 KiB 51 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 11 KiB 13 KiB 14 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 417 KiB 417 KiB 417 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 1966 B 3759 B 4846 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1342 B 1358 B 1427 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 129 KiB 129 KiB 129 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 6102 9719 B 10 KiB 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 259 B 261 B 271 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 436 B 438 B 449 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 30 B 32 B 52 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 7 B 8 B 11 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 111 B 113 B 116 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 136 B 183 B 219 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 64 32 KiB 32 KiB 32 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 2432 231 B 237 B 242 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1076 B 1080 B 1096 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 127 B 127 B 127 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1673 B 1732 B 1758 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 439 B 490 B 681 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 456 B 458 B 459 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 23 B 32 B 44 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 19 B 31 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 26 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