Initializing AMReX (25.06-18-gdf1e7d01d26d)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (25.06-18-gdf1e7d01d26d) initialized Starting run at 07:38:15 UTC on 2025-06-20. Successfully read inputs file ... Castro git describe: 25.06-4-g95261d7fd AMReX git describe: 25.06-18-gdf1e7d01d Microphysics git describe: 25.06-14-g3b767a4e 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.175763138 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001220556 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010580526 Gravity::multilevel_solve_for_new_phi() time = 0.112856601 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010597342 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.048383147 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004077403 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.006003083 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.066299779 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Castro::construct_old_gravity() time = 0.005751705 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.064036914 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.003463858 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.060870625 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.006052121 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.058997282 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.334364291 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.015762903 Gravity::solve_for_phi() time = 0.095507224 on level 0 Castro::construct_new_gravity() time = 0.097658915 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.059485816 on level 0 Castro::construct_old_gravity() time = 0.002910373 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.059012005 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.342365265 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.015440774 Gravity::solve_for_phi() time = 0.094057553 on level 0 Castro::construct_new_gravity() time = 0.096286311 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.060543043 on level 0 Castro::construct_old_gravity() time = 0.005925469 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.059030094 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.334439889 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.013950298 Gravity::solve_for_phi() time = 0.090907816 on level 0 Castro::construct_new_gravity() time = 0.093028923 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.059020633 on level 0 Castro::construct_old_gravity() time = 0.003119733 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.059351513 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.342102077 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.014618211 Gravity::solve_for_phi() time = 0.092298969 on level 0 Castro::construct_new_gravity() time = 0.094452069 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.059041147 on level 0 Castro::construct_old_gravity() time = 0.00297696 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.059876199 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.343853076 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.014521562 Gravity::solve_for_phi() time = 0.091000972 on level 0 Castro::construct_new_gravity() time = 0.0931847 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.058939524 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.003137496 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.058674454 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.341231116 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.014175301 Gravity::solve_for_phi() time = 0.090253247 on level 0 Castro::construct_new_gravity() time = 0.092362646 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.058913172 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.002833263 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.058795299 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.344880742 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.01431231 Gravity::solve_for_phi() time = 0.091454668 on level 0 Castro::construct_new_gravity() time = 0.093564443 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.059469761 on level 0 Castro::construct_old_gravity() time = 0.003229084 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.058706386 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.339593337 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.014362596 Gravity::solve_for_phi() time = 0.090653563 on level 0 Castro::construct_new_gravity() time = 0.092773343 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.058898688 on level 0 Castro::construct_old_gravity() time = 0.002999246 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.058708916 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.348040914 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.014919702 Gravity::solve_for_phi() time = 0.091163297 on level 0 Castro::construct_new_gravity() time = 0.093317127 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.059908454 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.003166756 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.058773628 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.345275055 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.015058827 Gravity::solve_for_phi() time = 0.094907735 on level 0 Castro::construct_new_gravity() time = 0.097031867 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.061422556 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.004546629 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.060110969 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.345717926 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.01454422 Gravity::solve_for_phi() time = 0.090454444 on level 0 Castro::construct_new_gravity() time = 0.092589437 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.059356987 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.003113259 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.059494729 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.342017923 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.014960814 Gravity::solve_for_phi() time = 0.090807083 on level 0 Castro::construct_new_gravity() time = 0.093006177 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.058835276 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.003173752 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.058608434 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.347289904 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.015094609 Gravity::solve_for_phi() time = 0.091310493 on level 0 Castro::construct_new_gravity() time = 0.093444514 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.05921474 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.003434628 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.060319635 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.342831821 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.014720001 Gravity::solve_for_phi() time = 0.090856495 on level 0 Castro::construct_new_gravity() time = 0.092971748 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.059178632 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.002762936 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.058799361 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.346152624 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.014937594 Gravity::solve_for_phi() time = 0.092266546 on level 0 Castro::construct_new_gravity() time = 0.094379516 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.059297069 on level 0 Castro::construct_old_gravity() time = 0.002779271 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.059712345 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.341397738 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.015050515 Gravity::solve_for_phi() time = 0.091241817 on level 0 Castro::construct_new_gravity() time = 0.093361496 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.059098709 on level 0 Castro::construct_old_gravity() time = 0.002980997 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.058813245 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.347149837 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002697586058 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.014887437 Gravity::solve_for_phi() time = 0.090912886 on level 0 Castro::construct_new_gravity() time = 0.093027155 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.059043183 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.20654781 [STEP 1] FAB kilobyte spread across MPI nodes: [6998 ... 6998] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.096278647 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004114182 seconds Ending run at 07:38:30 UTC on 2025-06-20. Run time = 14.57672524 Run time without initialization = 12.30786536 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.5671 ... 14.5701 ... 14.5768 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.6758 5.7157 5.7711 39.59% Amr::InitAmr() 1 1.352 1.355 1.362 9.34% Castro::computeTemp() 80 0.8265 0.9351 1.094 7.51% amrex::Copy() 724 0.795 0.8 0.8119 5.57% FabArray::ParallelCopy_finish() 1429 0.2354 0.5556 0.7646 5.25% Castro::normalize_species() 79 0.5499 0.5559 0.5845 4.01% MultiFab::AddProduct() 3439 0.5469 0.5511 0.5745 3.94% Castro::reset_internal_energy(Fab) 160 0.5119 0.5159 0.5219 3.58% MLPoisson::Fsmooth() 6144 0.3557 0.3851 0.4629 3.18% MultiFab::Subtract() 159 0.3509 0.3561 0.4059 2.78% Castro::initMFs() 1 0.3558 0.3587 0.3739 2.56% FabArray::setVal() 5013 0.3505 0.3542 0.3598 2.47% Gravity::fill_multipole_BCs() 19 0.1969 0.247 0.2673 1.83% FabArray::ParallelCopy_nowait() 1429 0.2019 0.2188 0.2451 1.68% check_for_negative_density() 20 0.08232 0.1543 0.2021 1.39% FabArray::sum() 3439 0.1917 0.194 0.1964 1.35% Castro::initData() 1 0.003426 0.02783 0.1622 1.11% FillBoundary_finish() 6380 0.1118 0.1312 0.1598 1.10% FillBoundary_nowait() 6380 0.09078 0.1201 0.1569 1.08% TagBoxArray::collate() 2 0.0003051 0.03101 0.1488 1.02% MLCellLinOp::applyBC() 7148 0.06242 0.07883 0.1247 0.86% MLPoisson::Fapply() 986 0.09369 0.09836 0.1218 0.84% Castro::construct_new_gravity_source() 17 0.0714 0.07273 0.0746 0.51% VisMF::Write(FabArray) 4 0.01406 0.01817 0.07308 0.50% FabArray::Saxpy() 91 0.06834 0.06864 0.06952 0.48% Amr::writePlotFile() 2 0.006805 0.05993 0.06374 0.44% FabArray::Xpay() 915 0.05087 0.05306 0.05874 0.40% Castro::construct_old_gravity_source() 20 0.04876 0.04893 0.0496 0.34% MultiFab::contains_nan() 37 0.04486 0.04526 0.04598 0.32% MLMG::addInterpCorrection() 768 0.03356 0.0363 0.044 0.30% AmrLevel::derive() 18 0.04031 0.04045 0.04077 0.28% amrex::average_down 768 0.02631 0.02871 0.03544 0.24% StateData::FillBoundary(geom) 143 0 0.01755 0.03379 0.23% Castro::reset_internal_energy(MultiFab) 80 0.02336 0.0238 0.02432 0.17% MLCellLinOp::defineAuxData() 19 0.01553 0.01661 0.01996 0.14% Castro::enforce_min_density() 79 0.01861 0.01885 0.0192 0.13% FabArray::FillBoundary() 6380 0.0136 0.01419 0.01705 0.12% amrex::Add() 128 0.01459 0.01486 0.01669 0.11% MLCGSolver::bicgstab 8 0 0.001029 0.01647 0.11% Castro::construct_new_source() 102 0.01595 0.01605 0.01635 0.11% Castro::evaluate_source_change() 196 0.01367 0.01409 0.01501 0.10% FabArray::setDomainBndry() 147 0.0004659 0.007794 0.01481 0.10% FillPatchIterator::Initialize 147 0.01295 0.01318 0.01425 0.10% FabArray::norminf() 244 0.01066 0.01103 0.01389 0.10% FabArrayBase::getFB() 6380 0.01042 0.0109 0.01332 0.09% Castro::wd_update() 2 0.005544 0.009 0.01265 0.09% MLCellLinOp::smooth() 1544 0.01019 0.01074 0.01206 0.08% MLMG::prepareForSolve() 18 0.01024 0.01058 0.0116 0.08% Castro::construct_old_source() 120 0.01014 0.01022 0.01042 0.07% Gravity::actual_solve_with_mlmg() 19 0.007295 0.007762 0.009177 0.06% MLMG::ResNormInf() 146 0.004766 0.006708 0.008708 0.06% FabArrayBase::getCPC() 2005 0.007589 0.007855 0.00821 0.06% MLMG::solve() 18 0.0007912 0.003945 0.008033 0.06% Gravity::get_old_grav_vector() 20 0.006509 0.006561 0.006646 0.05% Gravity::get_new_grav_vector() 18 0.00599 0.006082 0.006221 0.04% Castro::initialize_do_advance() 20 0.004565 0.00488 0.005801 0.04% MLCellLinOp::apply() 986 0.001993 0.00221 0.00465 0.03% FabArray::ParallelCopy() 1429 0.003436 0.003598 0.003856 0.03% AmrLevel::FillPatch() 147 0.00239 0.00249 0.002705 0.02% MLCellLinOp::correctionResidual() 768 0.001909 0.001963 0.002068 0.01% Castro::sum_integrated_quantities() 2 0.001458 0.001798 0.001972 0.01% MLMG::mgVcycle() 128 0.00158 0.001677 0.001836 0.01% MLMG:computeResOfCorrection() 768 0.001433 0.00148 0.001587 0.01% Castro::save_data_for_retry() 3 0.001316 0.001367 0.001526 0.01% MLMG::actualBottomSolve() 128 0.0001259 0.0001981 0.001231 0.01% Castro::subcycle_advance_ctu() 1 0.0002236 0.0003142 0.001043 0.01% StateDataPhysBCFunct::() 147 0.0007391 0.0008884 0.0009902 0.01% FillPatchIterator::FillFromLevel0() 147 0.0006873 0.0007171 0.0007785 0.01% main() 1 0.0006635 0.0007141 0.0007721 0.01% Gravity::solve_for_phi() 17 0.0006496 0.0006863 0.0007422 0.01% MLMG::mgVcycle_down::1 128 0.0005613 0.000577 0.0006273 0.00% MLMG::mgVcycle_down::2 128 0.0005432 0.0005821 0.0006256 0.00% AmrMesh::MakeNewGrids() 2 0.0005488 0.0005719 0.0006229 0.00% Castro::Castro() 1 0.0004567 0.0005484 0.0006055 0.00% MLMG::mgVcycle_down::3 128 0.000513 0.0005498 0.0006051 0.00% MLMG::mgVcycle_down::0 128 0.0005554 0.0005646 0.0005979 0.00% MLMG::mgVcycle_down::4 128 0.0004946 0.0005155 0.0005886 0.00% MLMG::mgVcycle_down::5 128 0.0005171 0.0005337 0.0005847 0.00% Castro::do_advance_ctu() 20 0.0003999 0.0004312 0.0005077 0.00% FillPatchSingleLevel 147 0.000467 0.0004847 0.0005076 0.00% Castro::clean_state() 79 0.0004266 0.0004552 0.0004947 0.00% MLCellLinOp::solutionResidual() 147 0.0004627 0.0004752 0.0004876 0.00% MLMG::oneIter() 128 0.0003868 0.0004021 0.0004384 0.00% MLMG::mgVcycle_up::0 128 0.0003588 0.0003782 0.0004358 0.00% MLMG::mgVcycle_up::1 128 0.0003584 0.0003765 0.0004216 0.00% MLMG::mgVcycle_up::4 128 0.0003253 0.0003365 0.0003875 0.00% MLMG::mgVcycle_up::2 128 0.0003241 0.0003475 0.0003691 0.00% MLMG::mgVcycle_up::3 128 0.0003208 0.0003328 0.0003663 0.00% MLMG::mgVcycle_up::5 128 0.0003123 0.0003213 0.0003396 0.00% Castro::do_old_sources() 20 0.000281 0.0002923 0.0003141 0.00% Castro::construct_new_gravity() 17 0.0001588 0.0001772 0.0002796 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002173 0.0002347 0.0002693 0.00% MLMG::computeResidual() 128 0.0002389 0.0002445 0.0002557 0.00% Castro::problem_post_timestep() 2 1.13e-05 2.777e-05 0.0002547 0.00% MLMG::mgVcycle_bottom 128 0.000218 0.0002262 0.0002412 0.00% Castro::do_new_sources() 17 0.0002124 0.0002241 0.0002403 0.00% Amr::coarseTimeStep() 1 0.0001433 0.0001707 0.0001812 0.00% Amr::defBaseLevel() 1 0.0001363 0.0001507 0.0001727 0.00% Gravity::test_composite_phi() 1 0.0001049 0.0001114 0.0001424 0.00% Castro::check_for_nan() 37 0.0001046 0.0001133 0.00014 0.00% TagBoxArray::mapPRD 2 0.0001267 0.0001323 0.0001384 0.00% MLPoisson::define() 19 0.0001199 0.0001277 0.0001365 0.00% Castro::construct_old_gravity() 20 0.0001068 0.0001149 0.0001339 0.00% Castro::retry_advance_ctu() 11 6.196e-05 7.943e-05 0.0001231 0.00% Castro::apply_source_to_state() 37 9.85e-05 0.0001039 0.0001209 0.00% Castro::initialize_advance() 1 0.0001045 0.0001104 0.0001177 0.00% Amr::FinalizeInit() 1 1.489e-05 2.89e-05 0.000105 0.00% Amr::grid_places() 2 5.591e-05 6.921e-05 8.811e-05 0.00% Castro::create_source_corrector() 17 5.855e-05 6.706e-05 8.766e-05 0.00% MLLinOp::define() 19 7.705e-05 8.1e-05 8.711e-05 0.00% Castro::expand_state() 20 4.921e-05 5.368e-05 6.894e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.993e-05 6.092e-05 6.49e-05 0.00% Castro::derive() 18 3.965e-05 4.158e-05 4.768e-05 0.00% MLMG::computeMLResidual() 18 4.037e-05 4.312e-05 4.564e-05 0.00% Amr::regrid() 1 1.547e-05 2.105e-05 4.295e-05 0.00% Gravity::actual_multilevel_solve() 1 3.078e-05 3.178e-05 3.487e-05 0.00% Amr::timeStep() 1 1.475e-05 1.945e-05 2.732e-05 0.00% Amr::InitializeInit() 1 7.976e-06 9.981e-06 2.493e-05 0.00% Castro::post_regrid() 2 1.075e-05 1.477e-05 2.155e-05 0.00% Castro::advance() 1 7.194e-06 1.264e-05 1.504e-05 0.00% Castro::post_init() 1 1.19e-05 1.248e-05 1.322e-05 0.00% Amr::bldFineLevels() 1 6.688e-06 7.588e-06 9.059e-06 0.00% Castro::post_timestep() 1 5.801e-06 6.354e-06 7.643e-06 0.00% Amr::init() 1 6.549e-06 7.125e-06 7.56e-06 0.00% Amr::initialInit() 1 2.858e-06 3.377e-06 6.618e-06 0.00% Other 2521 0.04863 0.06126 0.09549 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.31 12.31 12.31 84.43% Amr::timeStep() 1 12.21 12.21 12.21 83.74% Castro::advance() 1 12.14 12.15 12.15 83.34% Castro::subcycle_advance_ctu() 1 12.09 12.1 12.1 82.99% Castro::do_advance_ctu() 20 12.05 12.06 12.06 82.73% Castro::construct_ctu_hydro_source() 20 6.621 6.647 6.708 46.02% Castro::clean_state() 79 3.673 3.776 3.929 26.96% Castro::computeTemp() 80 2.335 2.442 2.599 17.83% Gravity::solve_phi_with_mlmg() 19 1.607 1.658 1.68 11.53% Castro::construct_new_gravity() 17 1.52 1.572 1.593 10.93% Gravity::solve_for_phi() 17 1.483 1.535 1.556 10.68% Castro::reset_internal_energy(MultiFab) 80 1.496 1.507 1.547 10.62% Castro::initialize_do_advance() 20 1.38 1.424 1.474 10.11% Gravity::actual_solve_with_mlmg() 19 1.404 1.408 1.412 9.69% Amr::InitAmr() 1 1.352 1.355 1.362 9.34% MLMG::solve() 18 1.349 1.355 1.357 9.31% MLMG::oneIter() 128 1.214 1.217 1.219 8.36% MLMG::mgVcycle() 128 1.197 1.202 1.204 8.26% Castro::do_old_sources() 20 1.154 1.169 1.192 8.18% Castro::do_new_sources() 17 0.9182 0.9664 1.008 6.92% FabArray::ParallelCopy() 1429 0.4804 0.7983 1.004 6.89% Castro::evaluate_source_change() 196 0.9513 0.9581 0.9828 6.74% Amr::init() 1 0.9066 0.9066 0.9066 6.22% Amr::initialInit() 1 0.8536 0.8537 0.8538 5.86% Castro::enforce_min_density() 79 0.7982 0.8054 0.8391 5.76% MLCellLinOp::smooth() 1544 0.6269 0.6832 0.8201 5.63% amrex::Copy() 724 0.795 0.8 0.8119 5.57% FabArray::ParallelCopy_finish() 1429 0.2354 0.5556 0.7646 5.25% AmrLevel::FillPatch() 147 0.5294 0.6043 0.6896 4.73% Amr::InitializeInit() 1 0.4415 0.4697 0.6157 4.22% Amr::defBaseLevel() 1 0.4415 0.4697 0.6156 4.22% Castro::normalize_species() 79 0.5499 0.5559 0.5845 4.01% MultiFab::AddProduct() 3439 0.5469 0.5511 0.5745 3.94% Castro::reset_internal_energy(Fab) 160 0.5119 0.5159 0.5219 3.58% FillPatchIterator::Initialize 147 0.3435 0.4217 0.5061 3.47% FillPatchIterator::FillFromLevel0() 147 0.3195 0.4007 0.4785 3.28% FillPatchSingleLevel 147 0.3188 0.4 0.4778 3.28% MLPoisson::Fsmooth() 6144 0.3557 0.3851 0.4629 3.18% MLMG::mgVcycle_up::1 128 0.1072 0.2837 0.462 3.17% MLCellLinOp::applyBC() 7148 0.3196 0.356 0.4346 2.98% Amr::FinalizeInit() 1 0.238 0.384 0.4122 2.83% MultiFab::Subtract() 159 0.3509 0.3561 0.4059 2.78% MLMG::addInterpCorrection() 768 0.0663 0.3239 0.4007 2.75% Castro::Castro() 1 0.3618 0.365 0.3814 2.62% Castro::initMFs() 1 0.3559 0.3589 0.374 2.57% FabArray::setVal() 5013 0.3505 0.3542 0.3598 2.47% FabArray::FillBoundary() 6380 0.2442 0.2771 0.3139 2.15% MLMG::mgVcycle_up::2 128 0.01341 0.127 0.3 2.06% MLMG::mgVcycle_down::0 128 0.2786 0.2803 0.2821 1.94% Castro::expand_state() 20 0.1944 0.2283 0.2754 1.89% Gravity::fill_multipole_BCs() 19 0.199 0.2492 0.2694 1.85% FabArray::ParallelCopy_nowait() 1429 0.2217 0.2391 0.2666 1.83% Castro::initData() 1 0.07782 0.103 0.2367 1.62% Amr::bldFineLevels() 1 0.04595 0.1921 0.22 1.51% MLMG::mgVcycle_up::0 128 0.2135 0.217 0.2184 1.50% MLCellLinOp::apply() 986 0.1549 0.1671 0.2138 1.47% check_for_negative_density() 20 0.08232 0.1543 0.2021 1.39% FabArray::sum() 3439 0.1917 0.194 0.1964 1.35% MLMG::mgVcycle_up::3 128 0.01335 0.02721 0.1907 1.31% Amr::grid_places() 2 0.002735 0.1487 0.1769 1.21% AmrMesh::MakeNewGrids() 2 0.002667 0.1486 0.1768 1.21% TagBoxArray::mapPRD 2 0.0008994 0.1162 0.1751 1.20% FillBoundary_nowait() 6380 0.1017 0.1318 0.1691 1.16% FillBoundary_finish() 6380 0.1118 0.1312 0.1598 1.10% MLMG:computeResOfCorrection() 768 0.1169 0.1268 0.1557 1.07% MLCellLinOp::correctionResidual() 768 0.1154 0.1253 0.1541 1.06% Castro::post_init() 1 0.1463 0.1488 0.1497 1.03% TagBoxArray::collate() 2 0.0003051 0.03101 0.1488 1.02% Amr::writePlotFile() 2 0.145 0.1451 0.1451 1.00% MLMG::mgVcycle_down::1 128 0.1265 0.1294 0.1339 0.92% Castro::construct_new_source() 102 0.1264 0.1279 0.1305 0.90% MLPoisson::Fapply() 986 0.09369 0.09836 0.1218 0.84% Gravity::multilevel_solve_for_new_phi() 1 0.1127 0.1128 0.1129 0.77% Gravity::actual_multilevel_solve() 1 0.1127 0.1127 0.1128 0.77% MLCellLinOp::solutionResidual() 147 0.09048 0.09557 0.09994 0.69% Castro::post_regrid() 2 0.08552 0.08639 0.08896 0.61% Castro::create_source_corrector() 17 0.05333 0.06842 0.08687 0.60% Castro::construct_old_source() 120 0.08186 0.08217 0.08323 0.57% MLMG::computeResidual() 128 0.07714 0.07849 0.0823 0.56% MLMG::mgVcycle_down::2 128 0.01569 0.0467 0.08144 0.56% Castro::construct_new_gravity_source() 17 0.0714 0.07273 0.0746 0.51% VisMF::Write(FabArray) 4 0.01406 0.01817 0.07308 0.50% MLMG::mgVcycle_bottom 128 0.0003463 0.004795 0.07141 0.49% MLMG::actualBottomSolve() 128 0.0001259 0.004569 0.07117 0.49% Castro::derive() 18 0.06957 0.06978 0.0704 0.48% AmrLevel::derive() 18 0.06953 0.06973 0.07036 0.48% FabArray::Saxpy() 91 0.06834 0.06864 0.06952 0.48% MLCGSolver::bicgstab 8 0 0.004332 0.06931 0.48% MLMG::mgVcycle_down::3 128 0.01489 0.02188 0.06776 0.46% Castro::construct_old_gravity() 20 0.04449 0.05953 0.0668 0.46% Gravity::get_old_grav_vector() 20 0.04437 0.05941 0.06669 0.46% Castro::post_timestep() 1 0.05869 0.06096 0.06388 0.44% amrex::average_down 768 0.03519 0.0425 0.05997 0.41% FabArray::Xpay() 915 0.05087 0.05306 0.05874 0.40% Castro::initialize_advance() 1 0.04909 0.04959 0.05081 0.35% Castro::construct_old_gravity_source() 20 0.04876 0.04893 0.0496 0.34% Castro::check_for_nan() 37 0.04498 0.04538 0.04612 0.32% MultiFab::contains_nan() 37 0.04486 0.04526 0.04598 0.32% Amr::regrid() 1 0.0442 0.04459 0.04586 0.31% MLPoisson::define() 19 0.03025 0.03552 0.04141 0.28% MLMG::mgVcycle_down::4 128 0.01482 0.0167 0.03981 0.27% Castro::retry_advance_ctu() 11 0.03633 0.03699 0.03863 0.26% Gravity::get_new_grav_vector() 18 0.03591 0.03642 0.03706 0.25% StateDataPhysBCFunct::() 147 0.0007391 0.01844 0.03472 0.24% StateData::FillBoundary(geom) 143 0 0.01755 0.03379 0.23% MLMG::mgVcycle_up::4 128 0.01344 0.01483 0.03147 0.22% Castro::apply_source_to_state() 37 0.02787 0.02798 0.02826 0.19% MLMG::mgVcycle_down::5 128 0.01487 0.01595 0.02767 0.19% MLCellLinOp::defineAuxData() 19 0.01967 0.02159 0.02702 0.19% MLMG::prepareForSolve() 18 0.02064 0.02153 0.02431 0.17% MLMG::mgVcycle_up::5 128 0.01353 0.01441 0.02367 0.16% MLMG::computeMLResidual() 18 0.01124 0.01655 0.02164 0.15% Castro::sum_integrated_quantities() 2 0.02101 0.02115 0.02117 0.15% Castro::save_data_for_retry() 3 0.01909 0.01936 0.02047 0.14% MLMG::ResNormInf() 146 0.01427 0.01611 0.01809 0.12% amrex::Add() 128 0.01459 0.01486 0.01669 0.11% Castro::problem_post_timestep() 2 0.009857 0.01282 0.01609 0.11% Castro::wd_update() 2 0.009845 0.01279 0.01608 0.11% FabArrayBase::getCPC() 2005 0.00994 0.01141 0.01523 0.10% FabArrayBase::getFB() 6380 0.01076 0.01163 0.01495 0.10% FabArray::setDomainBndry() 147 0.0004659 0.007794 0.01481 0.10% Gravity::test_composite_phi() 1 0.01435 0.01439 0.01443 0.10% FabArray::norminf() 244 0.01066 0.01103 0.01389 0.10% MLLinOp::define() 19 0.003727 0.007358 0.01103 0.08% Other 2521 0.08164 0.09994 0.1428 0.98% -------------------------------------------------------------------------------------------- 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 10 B 13 B 24 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 1570 KiB 1570 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 409 KiB 411 KiB 415 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 4273 B 4474 B 7269 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 50 KiB 50 KiB 52 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 11 KiB 13 KiB 15 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 414 KiB 414 KiB 415 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 2012 B 3865 B 5151 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1321 B 1326 B 1337 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 17 KiB 17 KiB 18 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 131 KiB 131 KiB 131 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 9715 B 10 KiB 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 253 B 253 B 257 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 426 B 427 B 433 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 18 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 6 B 7 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 111 B 112 B 115 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 139 B 186 B 209 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 1064 B 1071 B 1099 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 126 B 126 B 126 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1666 B 1724 B 1749 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 491 B 681 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 457 B 458 B 459 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 25 B 34 B 45 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 20 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 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-18-gdf1e7d01d26d) finalized