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:55:18 UTC on 2025-07-08. Successfully read inputs file ... Castro git describe: 25.07-5-g5d746e0a1 AMReX git describe: 25.06-22-g1d307c2c6-dirty Microphysics git describe: 25.07-3-g3881c8fa 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.164370284 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001205474 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010660416 Gravity::multilevel_solve_for_new_phi() time = 0.109562572 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010629157 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.048664914 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004081678 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.002724198 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.064571623 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.008002172 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.063812216 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.002664093 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.060016597 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.003760389 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.059468302 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.332885024 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.014508566 Gravity::solve_for_phi() time = 0.090696198 on level 0 Castro::construct_new_gravity() time = 0.092867282 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.059777684 on level 0 Castro::construct_old_gravity() time = 0.002785678 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.059187755 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.335140305 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.013980537 Gravity::solve_for_phi() time = 0.090483094 on level 0 Castro::construct_new_gravity() time = 0.092595334 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.059519993 on level 0 Castro::construct_old_gravity() time = 0.002698608 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.059237497 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.340935495 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.013866216 Gravity::solve_for_phi() time = 0.089243474 on level 0 Castro::construct_new_gravity() time = 0.09136918 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.059960909 on level 0 Castro::construct_old_gravity() time = 0.00466881 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.059029835 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.343711858 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.014285254 Gravity::solve_for_phi() time = 0.09109884 on level 0 Castro::construct_new_gravity() time = 0.093254241 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.060366372 on level 0 Castro::construct_old_gravity() time = 0.004497975 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.05961757 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.348470177 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.014855459 Gravity::solve_for_phi() time = 0.091556898 on level 0 Castro::construct_new_gravity() time = 0.093725328 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.059449523 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.002786977 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.059077185 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.341814229 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.014211088 Gravity::solve_for_phi() time = 0.090918011 on level 0 Castro::construct_new_gravity() time = 0.09302795 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.059331411 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.002712858 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.059180206 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.348346298 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.015363729 Gravity::solve_for_phi() time = 0.093459129 on level 0 Castro::construct_new_gravity() time = 0.095604851 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.05940913 on level 0 Castro::construct_old_gravity() time = 0.002732841 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.059060416 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.342706826 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.014457155 Gravity::solve_for_phi() time = 0.091315163 on level 0 Castro::construct_new_gravity() time = 0.09350119 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.05969837 on level 0 Castro::construct_old_gravity() time = 0.002845114 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.059313709 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.34791439 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.014693529 Gravity::solve_for_phi() time = 0.090883236 on level 0 Castro::construct_new_gravity() time = 0.093003382 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.059496046 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.002799135 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.059110906 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.342535148 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.014539497 Gravity::solve_for_phi() time = 0.091082772 on level 0 Castro::construct_new_gravity() time = 0.093206992 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.059458323 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.002724072 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.059243558 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.347904477 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.014747039 Gravity::solve_for_phi() time = 0.091036168 on level 0 Castro::construct_new_gravity() time = 0.093156471 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.059521047 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.002808777 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.059175153 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.342322726 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.014994329 Gravity::solve_for_phi() time = 0.092401255 on level 0 Castro::construct_new_gravity() time = 0.094521752 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.059425274 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.002661891 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.059135998 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.348530218 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.01534528 Gravity::solve_for_phi() time = 0.090925968 on level 0 Castro::construct_new_gravity() time = 0.093047083 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.059661418 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.003680259 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.058998957 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.346864769 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.015630855 Gravity::solve_for_phi() time = 0.092370475 on level 0 Castro::construct_new_gravity() time = 0.094544254 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.059830094 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.002904697 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.059193645 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.3488947 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.01492348 Gravity::solve_for_phi() time = 0.091823784 on level 0 Castro::construct_new_gravity() time = 0.093937882 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.059667249 on level 0 Castro::construct_old_gravity() time = 0.002791294 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.05908275 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.344276719 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.015088887 Gravity::solve_for_phi() time = 0.091299544 on level 0 Castro::construct_new_gravity() time = 0.093424886 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.059790854 on level 0 Castro::construct_old_gravity() time = 0.003069174 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.059260296 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.349678186 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002699805793 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.015273429 Gravity::solve_for_phi() time = 0.090908704 on level 0 Castro::construct_new_gravity() time = 0.093033802 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.059406829 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.19662874 [STEP 1] FAB kilobyte spread across MPI nodes: [6998 ... 6998] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.048669422 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004087252 seconds Ending run at 06:55:33 UTC on 2025-07-08. Run time = 14.50247443 Run time without initialization = 12.25054676 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.4919 ... 14.4953 ... 14.5028 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.7099 5.7423 5.8473 40.32% Amr::InitAmr() 1 1.345 1.348 1.356 9.35% Castro::computeTemp() 80 0.8433 0.9539 1.116 7.69% amrex::Copy() 724 0.7985 0.8016 0.8092 5.58% FabArray::ParallelCopy_finish() 1429 0.2208 0.5186 0.7046 4.86% Castro::normalize_species() 79 0.5465 0.5493 0.5631 3.88% MultiFab::AddProduct() 3439 0.551 0.5523 0.5582 3.85% Castro::reset_internal_energy(Fab) 160 0.5191 0.523 0.5416 3.73% MLPoisson::Fsmooth() 6144 0.3541 0.3822 0.4594 3.17% Castro::initMFs() 1 0.3609 0.3625 0.3651 2.52% FabArray::setVal() 5013 0.3474 0.3502 0.3544 2.44% MultiFab::Subtract() 159 0.3369 0.3377 0.3403 2.35% Gravity::fill_multipole_BCs() 19 0.195 0.2473 0.2675 1.84% FabArray::ParallelCopy_nowait() 1429 0.1984 0.2159 0.2371 1.63% check_for_negative_density() 20 0.03081 0.1446 0.212 1.46% FabArray::sum() 3439 0.1966 0.1977 0.1996 1.38% Castro::initData() 1 0.003359 0.02782 0.1623 1.12% TagBoxArray::collate() 2 0.0002566 0.02897 0.1611 1.11% FillBoundary_finish() 6380 0.106 0.1267 0.1609 1.11% FillBoundary_nowait() 6380 0.09067 0.12 0.153 1.05% MLCellLinOp::applyBC() 7148 0.06378 0.07992 0.1252 0.86% MLPoisson::Fapply() 986 0.09474 0.09931 0.1134 0.78% Castro::construct_new_gravity_source() 17 0.07226 0.07305 0.0745 0.51% FabArray::Saxpy() 91 0.07165 0.07198 0.07225 0.50% FabArray::Xpay() 915 0.05154 0.05361 0.05911 0.41% Castro::construct_old_gravity_source() 20 0.04926 0.04963 0.05089 0.35% MultiFab::contains_nan() 37 0.04593 0.04624 0.04679 0.32% MLMG::addInterpCorrection() 768 0.03371 0.03619 0.04351 0.30% AmrLevel::derive() 18 0.04101 0.04138 0.04169 0.29% amrex::average_down 768 0.02515 0.02738 0.03383 0.23% StateData::FillBoundary(geom) 143 0 0.0175 0.03351 0.23% VisMF::Write(FabArray) 4 0.01415 0.01506 0.02431 0.17% Castro::reset_internal_energy(MultiFab) 80 0.0231 0.02342 0.02381 0.16% MLCellLinOp::defineAuxData() 19 0.01556 0.01656 0.01989 0.14% Castro::enforce_min_density() 79 0.01854 0.01869 0.01907 0.13% FabArray::FillBoundary() 6380 0.01387 0.01443 0.01734 0.12% MLCGSolver::bicgstab 8 0 0.001016 0.01626 0.11% Castro::construct_new_source() 102 0.01599 0.01609 0.01621 0.11% Amr::writePlotFile() 2 0.006564 0.01454 0.01575 0.11% Castro::evaluate_source_change() 196 0.01386 0.01424 0.015 0.10% FabArray::setDomainBndry() 147 0.0004898 0.007612 0.0146 0.10% amrex::Add() 128 0.01444 0.01448 0.01458 0.10% FabArray::norminf() 244 0.01066 0.01106 0.01391 0.10% FillPatchIterator::Initialize 147 0.0128 0.01304 0.01332 0.09% FabArrayBase::getFB() 6380 0.01058 0.01099 0.01328 0.09% MLCellLinOp::smooth() 1544 0.01018 0.01071 0.01219 0.08% Castro::wd_update() 2 0.005027 0.009053 0.01216 0.08% MLMG::prepareForSolve() 18 0.01021 0.01054 0.01139 0.08% Castro::construct_old_source() 120 0.01026 0.0103 0.01038 0.07% MLMG::ResNormInf() 146 0.003336 0.007477 0.009436 0.07% Gravity::actual_solve_with_mlmg() 19 0.007269 0.007613 0.008912 0.06% FabArrayBase::getCPC() 2005 0.007745 0.007891 0.008206 0.06% MLMG::solve() 18 0.0008224 0.003469 0.007031 0.05% Gravity::get_old_grav_vector() 20 0.006419 0.006531 0.006774 0.05% Gravity::get_new_grav_vector() 18 0.006011 0.006064 0.006213 0.04% Castro::initialize_do_advance() 20 0.004485 0.004673 0.004911 0.03% MLCellLinOp::apply() 986 0.002048 0.002272 0.004805 0.03% FabArray::ParallelCopy() 1429 0.00361 0.003764 0.004 0.03% AmrLevel::FillPatch() 147 0.002352 0.002485 0.002697 0.02% MLCellLinOp::correctionResidual() 768 0.001906 0.00195 0.002038 0.01% Castro::sum_integrated_quantities() 2 0.001528 0.001798 0.001947 0.01% MLMG::mgVcycle() 128 0.001629 0.001709 0.001825 0.01% MLMG:computeResOfCorrection() 768 0.001511 0.001555 0.001664 0.01% Castro::save_data_for_retry() 3 0.001295 0.001358 0.001535 0.01% MLMG::actualBottomSolve() 128 0.000131 0.0001994 0.001178 0.01% Castro::subcycle_advance_ctu() 1 0.0002109 0.0002982 0.001085 0.01% main() 1 0.0006872 0.00082 0.001045 0.01% StateDataPhysBCFunct::() 147 0.0007264 0.0008782 0.001004 0.01% FillPatchIterator::FillFromLevel0() 147 0.0006929 0.000724 0.0007986 0.01% Gravity::solve_for_phi() 17 0.0006498 0.0006891 0.0007702 0.01% MLMG::mgVcycle_down::2 128 0.0005389 0.0005802 0.0006244 0.00% MLMG::mgVcycle_down::1 128 0.0005516 0.0005732 0.0006195 0.00% MLMG::mgVcycle_down::0 128 0.0005606 0.0005754 0.0006151 0.00% MLMG::mgVcycle_down::3 128 0.0005046 0.0005421 0.0006063 0.00% Castro::Castro() 1 0.000487 0.0005374 0.0006043 0.00% AmrMesh::MakeNewGrids() 2 0.0005419 0.0005646 0.0005871 0.00% MLMG::mgVcycle_down::4 128 0.000485 0.0005072 0.000586 0.00% MLMG::mgVcycle_down::5 128 0.0005066 0.0005228 0.0005737 0.00% FillPatchSingleLevel 147 0.0004632 0.0004817 0.0005008 0.00% Castro::clean_state() 79 0.000433 0.0004561 0.0004869 0.00% MLCellLinOp::solutionResidual() 147 0.0004445 0.0004568 0.0004686 0.00% MLMG::oneIter() 128 0.0003931 0.0004044 0.0004497 0.00% Castro::do_advance_ctu() 20 0.0003957 0.0004144 0.0004415 0.00% MLMG::mgVcycle_up::0 128 0.0003745 0.0003872 0.0004269 0.00% MLMG::mgVcycle_up::1 128 0.0003655 0.0003818 0.0004137 0.00% MLMG::mgVcycle_up::2 128 0.0003242 0.0003488 0.0003797 0.00% MLMG::mgVcycle_up::4 128 0.0003213 0.0003347 0.0003727 0.00% MLMG::mgVcycle_up::3 128 0.0003194 0.0003306 0.0003688 0.00% MLMG::mgVcycle_up::5 128 0.0003099 0.0003204 0.0003444 0.00% Castro::do_old_sources() 20 0.0002801 0.0002917 0.0003094 0.00% Castro::construct_new_gravity() 17 0.0001608 0.0001772 0.0003007 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002243 0.0002339 0.0002603 0.00% MLMG::mgVcycle_bottom 128 0.0002227 0.0002337 0.0002588 0.00% Castro::problem_post_timestep() 2 1.083e-05 2.712e-05 0.0002558 0.00% MLMG::computeResidual() 128 0.0002364 0.0002418 0.0002495 0.00% Castro::do_new_sources() 17 0.0002099 0.0002232 0.0002338 0.00% MLPoisson::define() 19 0.0001275 0.0001462 0.000166 0.00% Amr::coarseTimeStep() 1 0.0001349 0.00016 0.0001654 0.00% TagBoxArray::mapPRD 2 0.0001263 0.0001322 0.0001382 0.00% Amr::defBaseLevel() 1 0.0001107 0.0001192 0.0001324 0.00% Castro::check_for_nan() 37 0.0001081 0.0001176 0.0001323 0.00% Castro::retry_advance_ctu() 11 6.347e-05 7.602e-05 0.000131 0.00% Castro::construct_old_gravity() 20 0.0001093 0.0001194 0.000131 0.00% Gravity::test_composite_phi() 1 0.000105 0.0001093 0.000122 0.00% Castro::initialize_advance() 1 0.0001059 0.0001123 0.0001182 0.00% Castro::apply_source_to_state() 37 9.651e-05 0.000102 0.000117 0.00% Amr::FinalizeInit() 1 1.508e-05 2.983e-05 9.402e-05 0.00% MLLinOp::define() 19 7.322e-05 7.857e-05 8.416e-05 0.00% Amr::grid_places() 2 5.35e-05 6.375e-05 7.68e-05 0.00% Castro::create_source_corrector() 17 5.516e-05 6.259e-05 7.582e-05 0.00% Castro::expand_state() 20 5.177e-05 5.543e-05 6.631e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.884e-05 6.041e-05 6.493e-05 0.00% MLMG::computeMLResidual() 18 3.928e-05 4.459e-05 6.349e-05 0.00% Castro::derive() 18 4.228e-05 4.431e-05 4.908e-05 0.00% Amr::regrid() 1 1.553e-05 2.265e-05 4.283e-05 0.00% Gravity::actual_multilevel_solve() 1 3.133e-05 3.275e-05 3.517e-05 0.00% Castro::post_init() 1 1.181e-05 1.367e-05 3.008e-05 0.00% Amr::timeStep() 1 1.396e-05 1.643e-05 2.347e-05 0.00% Castro::post_regrid() 2 1.142e-05 1.39e-05 1.969e-05 0.00% Castro::advance() 1 6.65e-06 8.887e-06 1.528e-05 0.00% Amr::InitializeInit() 1 7.772e-06 9.135e-06 1.434e-05 0.00% Amr::bldFineLevels() 1 7.003e-06 8.369e-06 1.328e-05 0.00% Amr::init() 1 5.89e-06 7.581e-06 1.192e-05 0.00% Castro::post_timestep() 1 5.744e-06 6.428e-06 7.45e-06 0.00% Amr::initialInit() 1 2.703e-06 3.214e-06 6.516e-06 0.00% Other 2521 0.04884 0.0601 0.09246 0.64% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.49 14.5 14.5 100.00% Amr::coarseTimeStep() 1 12.25 12.25 12.25 84.46% Amr::timeStep() 1 12.2 12.2 12.2 84.10% Castro::advance() 1 12.13 12.14 12.14 83.69% Castro::subcycle_advance_ctu() 1 12.08 12.09 12.09 83.35% Castro::do_advance_ctu() 20 12.04 12.05 12.05 83.08% Castro::construct_ctu_hydro_source() 20 6.635 6.665 6.709 46.26% Castro::clean_state() 79 3.666 3.78 3.939 27.16% Castro::computeTemp() 80 2.345 2.459 2.62 18.07% Gravity::solve_phi_with_mlmg() 19 1.593 1.647 1.671 11.52% Castro::construct_new_gravity() 17 1.509 1.563 1.587 10.94% Gravity::solve_for_phi() 17 1.472 1.526 1.551 10.69% Castro::reset_internal_energy(MultiFab) 80 1.499 1.505 1.525 10.51% Castro::initialize_do_advance() 20 1.379 1.408 1.435 9.89% Gravity::actual_solve_with_mlmg() 19 1.392 1.396 1.4 9.65% Amr::InitAmr() 1 1.345 1.348 1.356 9.35% MLMG::solve() 18 1.337 1.343 1.345 9.27% MLMG::oneIter() 128 1.203 1.205 1.208 8.33% Castro::do_old_sources() 20 1.139 1.168 1.194 8.23% MLMG::mgVcycle() 128 1.189 1.19 1.193 8.22% Castro::do_new_sources() 17 0.9221 0.9699 1.013 6.99% Castro::evaluate_source_change() 196 0.9593 0.9618 0.972 6.70% FabArray::ParallelCopy() 1429 0.4608 0.7587 0.9436 6.51% Amr::init() 1 0.8964 0.8964 0.8964 6.18% Amr::initialInit() 1 0.8431 0.8431 0.8432 5.81% MLCellLinOp::smooth() 1544 0.6242 0.6788 0.8198 5.65% amrex::Copy() 724 0.7985 0.8016 0.8092 5.58% Castro::enforce_min_density() 79 0.7957 0.7982 0.8064 5.56% FabArray::ParallelCopy_finish() 1429 0.2208 0.5186 0.7046 4.86% AmrLevel::FillPatch() 147 0.5053 0.5777 0.6418 4.43% Amr::InitializeInit() 1 0.4473 0.4733 0.6102 4.21% Amr::defBaseLevel() 1 0.4473 0.4733 0.6102 4.21% Castro::normalize_species() 79 0.5465 0.5493 0.5631 3.88% MultiFab::AddProduct() 3439 0.551 0.5523 0.5582 3.85% Castro::reset_internal_energy(Fab) 160 0.5191 0.523 0.5416 3.73% MLPoisson::Fsmooth() 6144 0.3541 0.3822 0.4594 3.17% FillPatchIterator::Initialize 147 0.3224 0.3946 0.4594 3.17% MLMG::mgVcycle_up::1 128 0.1082 0.282 0.4572 3.15% MLCellLinOp::applyBC() 7148 0.3192 0.3528 0.4353 3.00% FillPatchIterator::FillFromLevel0() 147 0.2951 0.374 0.432 2.98% FillPatchSingleLevel 147 0.2944 0.3732 0.4312 2.97% Amr::FinalizeInit() 1 0.2329 0.3699 0.3959 2.73% MLMG::addInterpCorrection() 768 0.06614 0.3179 0.3932 2.71% Castro::Castro() 1 0.367 0.3686 0.3712 2.56% Castro::initMFs() 1 0.3611 0.3626 0.3653 2.52% FabArray::setVal() 5013 0.3474 0.3502 0.3544 2.44% MultiFab::Subtract() 159 0.3369 0.3377 0.3403 2.35% FabArray::FillBoundary() 6380 0.2437 0.2729 0.3146 2.17% MLMG::mgVcycle_up::2 128 0.01374 0.1254 0.2945 2.03% MLMG::mgVcycle_down::0 128 0.275 0.2768 0.2798 1.93% Gravity::fill_multipole_BCs() 19 0.1971 0.2494 0.2697 1.86% FabArray::ParallelCopy_nowait() 1429 0.2181 0.2363 0.2586 1.78% Castro::expand_state() 20 0.1955 0.2195 0.2498 1.72% Castro::initData() 1 0.07828 0.103 0.2374 1.64% check_for_negative_density() 20 0.03081 0.1446 0.212 1.46% MLMG::mgVcycle_up::0 128 0.2102 0.2109 0.2117 1.46% Amr::bldFineLevels() 1 0.04596 0.1829 0.2088 1.44% MLCellLinOp::apply() 986 0.1559 0.1665 0.2031 1.40% FabArray::sum() 3439 0.1966 0.1977 0.1996 1.38% MLMG::mgVcycle_up::3 128 0.01368 0.02726 0.1883 1.30% Amr::grid_places() 2 0.002604 0.1396 0.1656 1.14% AmrMesh::MakeNewGrids() 2 0.002547 0.1395 0.1655 1.14% FillBoundary_nowait() 6380 0.1015 0.1318 0.1652 1.14% TagBoxArray::mapPRD 2 0.0009209 0.1092 0.1637 1.13% TagBoxArray::collate() 2 0.0002566 0.02897 0.1611 1.11% FillBoundary_finish() 6380 0.106 0.1267 0.1609 1.11% MLMG:computeResOfCorrection() 768 0.1171 0.1276 0.1506 1.04% MLCellLinOp::correctionResidual() 768 0.1156 0.126 0.149 1.03% Castro::post_init() 1 0.143 0.1438 0.144 0.99% MLMG::mgVcycle_down::1 128 0.1266 0.1301 0.1335 0.92% Castro::construct_new_source() 102 0.1283 0.1292 0.1307 0.90% MLPoisson::Fapply() 986 0.09474 0.09931 0.1134 0.78% Gravity::multilevel_solve_for_new_phi() 1 0.1095 0.1095 0.1096 0.76% Gravity::actual_multilevel_solve() 1 0.1094 0.1095 0.1095 0.76% MLCellLinOp::solutionResidual() 147 0.08879 0.0948 0.09927 0.68% Amr::writePlotFile() 2 0.0977 0.09779 0.09781 0.67% Castro::post_regrid() 2 0.08621 0.08647 0.08726 0.60% Castro::construct_old_source() 120 0.08306 0.08353 0.0849 0.59% MLMG::computeResidual() 128 0.0771 0.07815 0.08296 0.57% MLMG::mgVcycle_down::2 128 0.01602 0.04655 0.07914 0.55% Castro::create_source_corrector() 17 0.05049 0.06022 0.07912 0.55% Castro::construct_new_gravity_source() 17 0.07226 0.07305 0.0745 0.51% FabArray::Saxpy() 91 0.07165 0.07198 0.07225 0.50% Castro::derive() 18 0.07046 0.07103 0.07163 0.49% AmrLevel::derive() 18 0.07042 0.07099 0.07158 0.49% MLMG::mgVcycle_bottom 128 0.0003554 0.004797 0.07126 0.49% MLMG::actualBottomSolve() 128 0.000131 0.004563 0.071 0.49% MLCGSolver::bicgstab 8 0 0.004324 0.06918 0.48% MLMG::mgVcycle_down::3 128 0.01527 0.02199 0.06668 0.46% Castro::post_timestep() 1 0.05909 0.06133 0.06416 0.44% Castro::construct_old_gravity() 20 0.04074 0.05422 0.06222 0.43% Gravity::get_old_grav_vector() 20 0.04062 0.0541 0.06209 0.43% FabArray::Xpay() 915 0.05154 0.05361 0.05911 0.41% amrex::average_down 768 0.03426 0.04118 0.05873 0.40% Castro::construct_old_gravity_source() 20 0.04926 0.04963 0.05089 0.35% Castro::initialize_advance() 1 0.04921 0.04939 0.05007 0.35% Castro::check_for_nan() 37 0.04604 0.04635 0.04691 0.32% MultiFab::contains_nan() 37 0.04593 0.04624 0.04679 0.32% Amr::regrid() 1 0.04436 0.04454 0.04499 0.31% Castro::retry_advance_ctu() 11 0.03624 0.03732 0.04343 0.30% MLPoisson::define() 19 0.03065 0.03601 0.04175 0.29% MLMG::mgVcycle_down::4 128 0.01527 0.01704 0.03967 0.27% Gravity::get_new_grav_vector() 18 0.03595 0.03633 0.03673 0.25% StateDataPhysBCFunct::() 147 0.0007264 0.01838 0.03446 0.24% StateData::FillBoundary(geom) 143 0 0.0175 0.03351 0.23% MLMG::mgVcycle_up::4 128 0.01374 0.01504 0.03131 0.22% Castro::apply_source_to_state() 37 0.02917 0.02937 0.02956 0.20% MLMG::mgVcycle_down::5 128 0.01522 0.01627 0.02771 0.19% MLCellLinOp::defineAuxData() 19 0.0197 0.02157 0.02689 0.19% Castro::save_data_for_retry() 3 0.01928 0.01996 0.02597 0.18% VisMF::Write(FabArray) 4 0.01415 0.01506 0.02431 0.17% MLMG::prepareForSolve() 18 0.02055 0.02151 0.02406 0.17% MLMG::mgVcycle_up::5 128 0.01386 0.01469 0.02375 0.16% Castro::sum_integrated_quantities() 2 0.02129 0.02146 0.02147 0.15% MLMG::computeMLResidual() 18 0.0108 0.01616 0.02105 0.15% MLMG::ResNormInf() 146 0.01283 0.01693 0.01887 0.13% Castro::problem_post_timestep() 2 0.01023 0.01279 0.01558 0.11% Castro::wd_update() 2 0.01022 0.01276 0.01557 0.11% FabArrayBase::getCPC() 2005 0.01002 0.01145 0.0151 0.10% FabArrayBase::getFB() 6380 0.01087 0.01172 0.0149 0.10% FabArray::setDomainBndry() 147 0.0004898 0.007612 0.0146 0.10% amrex::Add() 128 0.01444 0.01448 0.01458 0.10% Gravity::test_composite_phi() 1 0.01429 0.01433 0.01436 0.10% FabArray::norminf() 244 0.01066 0.01106 0.01391 0.10% MLLinOp::define() 19 0.003625 0.007872 0.01153 0.08% Other 2521 0.08119 0.09731 0.1375 0.95% -------------------------------------------------------------------------------------------- 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 20 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 1571 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 413 KiB 415 KiB 423 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1040 KiB 1041 KiB 1041 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 4332 B 4395 B 4852 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 12 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 2001 B 3639 B 4808 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1347 B 1358 B 1369 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 128 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 6102 9677 B 10 KiB 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 260 B 260 B 261 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 436 B 437 B 440 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 19 B 21 B 37 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 51 B 51 B 52 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 3 B 3 B 4 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 111 B 112 B 114 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 127 B 170 B 196 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 231 B 237 B 241 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1078 B 1081 B 1093 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 123 B 123 B 123 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1663 B 1723 B 1751 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 438 B 489 B 678 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 455 B 456 B 457 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 24 B 32 B 43 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