Initializing AMReX (26.07-51-g50c73c8ed200)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.07-51-g50c73c8ed200) initialized Starting run at 06:45:55 UTC on 2026-07-27. Successfully read inputs file ... Castro git describe: 26.07-7-g553efa2e4 AMReX git describe: 26.07-51-g50c73c8ed Microphysics git describe: 26.07-3-g6beccff2 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.163613496 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001161568 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010024162 Gravity::multilevel_solve_for_new_phi() time = 0.10471558 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.009969452 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.046364168 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004082431 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.002576251 on level 0 Castro::do_old_sources() time = 0.007393271 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.052433863 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.002969132 on level 0 Castro::do_old_sources() time = 0.007421242 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.05195523 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.003110767 on level 0 Castro::do_old_sources() time = 0.007344064 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.049745174 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.003778725 on level 0 Castro::do_old_sources() time = 0.007936806 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.049141911 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.29642878 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.013464599 Gravity::solve_for_phi() time = 0.088343941 on level 0 Castro::construct_new_gravity() time = 0.090415415 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.048635438 on level 0 Castro::construct_old_gravity() time = 0.002689574 on level 0 Castro::do_old_sources() time = 0.007348242 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.048029885 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.295027869 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.01359973 Gravity::solve_for_phi() time = 0.088066694 on level 0 Castro::construct_new_gravity() time = 0.09013428 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.048518344 on level 0 Castro::construct_old_gravity() time = 0.002605616 on level 0 Castro::do_old_sources() time = 0.007792206 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.048418596 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.298356238 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.014032857 Gravity::solve_for_phi() time = 0.088682955 on level 0 Castro::construct_new_gravity() time = 0.090661267 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.048551477 on level 0 Castro::construct_old_gravity() time = 0.002613174 on level 0 Castro::do_old_sources() time = 0.007619002 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.048070939 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.302183854 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.0137297 Gravity::solve_for_phi() time = 0.088121895 on level 0 Castro::construct_new_gravity() time = 0.090189796 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.048700191 on level 0 Castro::construct_old_gravity() time = 0.002655533 on level 0 Castro::do_old_sources() time = 0.007311428 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.048181767 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.309152892 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.014699997 Gravity::solve_for_phi() time = 0.090636304 on level 0 Castro::construct_new_gravity() time = 0.092713213 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.049070968 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.004701771 on level 0 Castro::do_old_sources() time = 0.008083822 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.048250954 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.302696637 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.014341996 Gravity::solve_for_phi() time = 0.088932436 on level 0 Castro::construct_new_gravity() time = 0.090960343 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.048624096 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.002601789 on level 0 Castro::do_old_sources() time = 0.007300764 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.048492518 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.309026292 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.01445635 Gravity::solve_for_phi() time = 0.089312338 on level 0 Castro::construct_new_gravity() time = 0.091307755 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.04858377 on level 0 Castro::construct_old_gravity() time = 0.00264732 on level 0 Castro::do_old_sources() time = 0.007411285 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.048254208 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.304443517 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.013905452 Gravity::solve_for_phi() time = 0.0889902 on level 0 Castro::construct_new_gravity() time = 0.090988408 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.048783926 on level 0 Castro::construct_old_gravity() time = 0.00302863 on level 0 Castro::do_old_sources() time = 0.007449181 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.048306073 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.307694387 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.013977098 Gravity::solve_for_phi() time = 0.089020481 on level 0 Castro::construct_new_gravity() time = 0.091033588 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.048643862 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.003597259 on level 0 Castro::do_old_sources() time = 0.008479681 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.049412332 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.305557609 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.01413441 Gravity::solve_for_phi() time = 0.088729869 on level 0 Castro::construct_new_gravity() time = 0.090760192 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.048755502 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.002628476 on level 0 Castro::do_old_sources() time = 0.007436234 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.048384505 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.310527703 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.014053306 Gravity::solve_for_phi() time = 0.090037928 on level 0 Castro::construct_new_gravity() time = 0.092081302 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.048556811 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.002725916 on level 0 Castro::do_old_sources() time = 0.007389981 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.048298873 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.3033347 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.013687917 Gravity::solve_for_phi() time = 0.087790136 on level 0 Castro::construct_new_gravity() time = 0.08984004 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.048618715 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.002582089 on level 0 Castro::do_old_sources() time = 0.007478317 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.048319041 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.309822 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.013802572 Gravity::solve_for_phi() time = 0.08948637 on level 0 Castro::construct_new_gravity() time = 0.091484815 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.048458589 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.002719363 on level 0 Castro::do_old_sources() time = 0.007801149 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.048457525 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.305478883 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.014742026 Gravity::solve_for_phi() time = 0.089425059 on level 0 Castro::construct_new_gravity() time = 0.091445238 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.048582463 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.002607143 on level 0 Castro::do_old_sources() time = 0.007383963 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.048311353 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.310403374 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.01415716 Gravity::solve_for_phi() time = 0.089069247 on level 0 Castro::construct_new_gravity() time = 0.091125408 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.050004272 on level 0 Castro::construct_old_gravity() time = 0.002700736 on level 0 Castro::do_old_sources() time = 0.007326318 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.048286117 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.306566579 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.014528051 Gravity::solve_for_phi() time = 0.089073319 on level 0 Castro::construct_new_gravity() time = 0.091143275 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.048631024 on level 0 Castro::construct_old_gravity() time = 0.003581775 on level 0 Castro::do_old_sources() time = 0.008010849 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.048219088 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.313958855 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002957365414 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.014684516 Gravity::solve_for_phi() time = 0.089306799 on level 0 Castro::construct_new_gravity() time = 0.091374298 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.048635839 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.581215356e+29 ener added: 1.581215356e+29 TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -4.341809003e+38 TIME= 0.25 YMOM = 1.137357411e+28 TIME= 0.25 ZMOM = 0 TIME= 0.25 ANG MOM X = 0 TIME= 0.25 ANG MOM Y = 0 TIME= 0.25 ANG MOM Z = -1.469316165e+38 TIME= 0.25 RHO*e = 1.325894512e+50 TIME= 0.25 RHO*K = 3.170941987e+49 TIME= 0.25 RHO*E = 1.642988711e+50 TIME= 0.25 RHO*PHI = -5.449011378e+50 TIME= 0.25 TOTAL ENERGY= -1.081516979e+50 TIME= 0.25 CENTER OF MASS X-LOC = 316424388.8 TIME= 0.25 CENTER OF MASS X-VEL = -170267.2489 TIME= 0.25 CENTER OF MASS Y-LOC = 2.427336677e-06 TIME= 0.25 CENTER OF MASS Y-VEL = 4.460231145e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 57600050.72 TIME= 0.25 MAXIMUM DENSITY = 4458760.788 TIME= 0.25 MAXIMUM T_S / T_E = 0 STEP = 1 TIME = 0.25 DT = 0.25 [STEP 1] Coarse TimeStep time: 11.13146033 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.045901394 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004134472 seconds Ending run at 06:46:08 UTC on 2026-07-27. Run time = 13.33167611 Run time without initialization = 11.18259472 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]: 13.3220 ... 13.3258 ... 13.3317 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.1314 5.1542 5.2372 39.28% Amr::InitAmr() 1 1.257 1.261 1.267 9.50% Castro::computeTemp() 80 0.7043 0.7948 0.9165 6.87% amrex::Copy() 766 0.7614 0.7642 0.7764 5.82% FabArray::ParallelCopy_finish() 1469 0.2181 0.5119 0.6753 5.07% MultiFab::AddProduct() 3598 0.5358 0.5378 0.5441 4.08% MLPoisson::Fsmooth() 6144 0.3246 0.3502 0.422 3.17% Castro::initMFs() 1 0.386 0.3898 0.4069 3.05% FabArray::setVal() 5232 0.3585 0.3613 0.3671 2.75% Castro::reset_internal_energy(Fab) 160 0.3455 0.3482 0.3518 2.64% MultiFab::Subtract() 159 0.3304 0.3321 0.3363 2.52% FabArray::ParallelCopy_nowait() 1469 0.2267 0.2455 0.2672 2.00% Castro::do_advance_ctu() 20 0.2558 0.2588 0.2633 1.98% Gravity::fill_multipole_BCs() 19 0.1887 0.2403 0.2568 1.93% check_for_negative_density() 20 0.03873 0.1463 0.1848 1.39% FabArray::sum() 3598 0.1793 0.18 0.1818 1.36% FillBoundary_finish() 6380 0.09922 0.1285 0.1653 1.24% Castro::initData() 1 0.003142 0.02805 0.1651 1.24% FillBoundary_nowait() 6380 0.08857 0.1179 0.1516 1.14% TagBoxArray::collate() 2 0.0002555 0.02519 0.1379 1.03% MLCellLinOp::applyBC() 7148 0.06902 0.08507 0.1317 0.99% Castro::normalize_species() 79 0.1192 0.1198 0.1216 0.91% MLPoisson::Fapply() 986 0.08494 0.08901 0.103 0.77% Castro::construct_old_gravity_source() 40 0.0965 0.09689 0.09849 0.74% FabArray::Saxpy() 111 0.07945 0.07967 0.08133 0.61% Castro::construct_new_gravity_source() 17 0.06659 0.06765 0.06953 0.52% FabArray::Xpay() 915 0.04872 0.05109 0.0564 0.42% MultiFab::contains_nan() 37 0.04579 0.04616 0.04722 0.35% MLMG::addInterpCorrection() 768 0.032 0.03458 0.04181 0.31% AmrLevel::derive() 18 0.03859 0.03875 0.03918 0.29% StateData::FillBoundary(geom) 193 0 0.01824 0.0352 0.26% amrex::average_down 768 0.02398 0.02621 0.03271 0.25% Castro::reset_internal_energy(MultiFab) 80 0.02519 0.02549 0.02593 0.19% VisMF::Write(FabArray) 4 0.01257 0.014 0.02381 0.18% Castro::enforce_min_density() 79 0.01907 0.01952 0.02322 0.17% MLCellLinOp::defineAuxData() 19 0.01602 0.01702 0.02076 0.16% FabArray::FillBoundary() 6380 0.01649 0.01704 0.02026 0.15% Castro::construct_old_source() 240 0.01906 0.01925 0.01965 0.15% FillPatchIterator::Initialize 187 0.01605 0.01634 0.01697 0.13% FabArray::setDomainBndry() 187 0.0006726 0.008923 0.01671 0.13% MLCGSolver::bicgstab 8 0 0.001042 0.01668 0.13% Amr::writePlotFile() 2 0.007077 0.01552 0.0165 0.12% FabArrayBase::getFB() 6380 0.01306 0.01358 0.01637 0.12% Castro::evaluate_source_change() 196 0.01421 0.01452 0.01572 0.12% Castro::construct_new_source() 102 0.01495 0.01508 0.01531 0.11% FabArray::norminf() 244 0.0102 0.01065 0.01406 0.11% amrex::Add() 128 0.0135 0.01356 0.01386 0.10% Castro::wd_update() 2 0.004152 0.009536 0.01238 0.09% MLCellLinOp::smooth() 1544 0.01048 0.01085 0.0123 0.09% MLMG::prepareForSolve() 18 0.01022 0.01058 0.0117 0.09% MLMG::ResNormInf() 146 0.002218 0.007392 0.009612 0.07% Gravity::actual_solve_with_mlmg() 19 0.007414 0.007887 0.009578 0.07% FabArrayBase::getCPC() 2045 0.008721 0.008881 0.009522 0.07% MLMG::solve() 18 0.0007627 0.004105 0.008714 0.07% Gravity::get_old_grav_vector() 20 0.006014 0.006144 0.006986 0.05% Gravity::get_new_grav_vector() 18 0.005704 0.005794 0.005991 0.04% MLCGSolver::ParallelAllReduce 189 0 0.0003453 0.005525 0.04% MLCellLinOp::apply() 986 0.002346 0.002585 0.005449 0.04% Castro::initialize_do_advance() 20 0.004681 0.004914 0.005322 0.04% FabArray::ParallelCopy() 1469 0.004129 0.004325 0.004775 0.04% AmrLevel::FillPatch() 187 0.003097 0.003205 0.003568 0.03% MLLinOp::defineGrids() 19 0.002583 0.002727 0.003117 0.02% MLMG::mgVcycle() 128 0.002154 0.002229 0.002371 0.02% MLCellLinOp::correctionResidual() 768 0.002209 0.002267 0.002337 0.02% MLMG:computeResOfCorrection() 768 0.00177 0.001812 0.00195 0.01% Castro::sum_integrated_quantities() 2 0.001558 0.001767 0.001943 0.01% Castro::save_data_for_retry() 3 0.001314 0.001335 0.001384 0.01% StateDataPhysBCFunct::() 187 0.001018 0.001181 0.00133 0.01% MLMG::actualBottomSolve() 128 0.000184 0.000256 0.001256 0.01% main() 1 0.0007954 0.0008985 0.001123 0.01% FillPatchIterator::FillFromLevel0() 187 0.0008776 0.0009214 0.0009847 0.01% Castro::subcycle_advance_ctu() 1 0.0002362 0.0003032 0.0008688 0.01% Gravity::solve_for_phi() 17 0.0006831 0.0007333 0.0008146 0.01% FillPatchSingleLevel 187 0.0006385 0.0006602 0.0006875 0.01% MLMG::mgVcycle_down::0 128 0.0006121 0.0006268 0.000674 0.01% MLMG::mgVcycle_down::5 128 0.0005781 0.0005966 0.0006673 0.01% MLMG::mgVcycle_down::1 128 0.0005979 0.0006146 0.0006669 0.01% MLMG::mgVcycle_down::2 128 0.0005783 0.0006138 0.0006661 0.00% Castro::Castro() 1 0.0004784 0.000544 0.0006469 0.00% MLMG::mgVcycle_down::3 128 0.0005531 0.0005856 0.0006357 0.00% MLMG::mgVcycle_down::4 128 0.0005508 0.000566 0.000612 0.00% AmrMesh::MakeNewGrids() 2 0.0005737 0.0005899 0.0006081 0.00% MLCellLinOp::solutionResidual() 147 0.0005043 0.0005171 0.0005597 0.00% Castro::clean_state() 79 0.0004669 0.0005004 0.0005477 0.00% MLMG::oneIter() 128 0.000429 0.0004467 0.0004985 0.00% Castro::do_old_sources() 40 0.0004433 0.0004587 0.0004977 0.00% MLMG::mgVcycle_up::1 128 0.0004069 0.0004232 0.0004792 0.00% MLMG::mgVcycle_up::0 128 0.0004003 0.0004127 0.0004522 0.00% MLMG::mgVcycle_up::2 128 0.0003871 0.0004087 0.0004432 0.00% MLMG::mgVcycle_up::3 128 0.0003732 0.0003878 0.0004339 0.00% MLMG::mgVcycle_up::4 128 0.0003783 0.000386 0.0004309 0.00% MLMG::mgVcycle_up::5 128 0.0003687 0.0003764 0.0003966 0.00% MLMG::computeResidual() 128 0.0002831 0.0002934 0.0003146 0.00% Castro::construct_new_gravity() 17 0.0001961 0.0002175 0.0003108 0.00% MLMG::mgVcycle_bottom 128 0.0002729 0.000281 0.0002965 0.00% Gravity::solve_phi_with_mlmg() 19 0.000217 0.0002431 0.0002804 0.00% Castro::problem_post_timestep() 2 1.197e-05 2.822e-05 0.0002541 0.00% Castro::do_new_sources() 17 0.0002034 0.0002167 0.0002311 0.00% Amr::coarseTimeStep() 1 0.000131 0.0001575 0.0001628 0.00% Castro::check_for_nan() 37 0.000116 0.0001292 0.000161 0.00% MLPoisson::define() 19 0.0001304 0.0001425 0.0001604 0.00% Castro::construct_old_gravity() 20 0.0001243 0.0001335 0.0001509 0.00% Amr::defBaseLevel() 1 0.0001032 0.0001279 0.0001503 0.00% Gravity::test_composite_phi() 1 0.0001124 0.000118 0.0001325 0.00% TagBoxArray::mapPRD 2 0.0001197 0.0001238 0.000132 0.00% Castro::apply_source_to_state() 37 0.0001118 0.0001189 0.0001318 0.00% Castro::initialize_advance() 1 0.0001101 0.0001164 0.0001294 0.00% Castro::retry_advance_ctu() 11 6.801e-05 7.519e-05 0.0001277 0.00% Amr::grid_places() 2 6.026e-05 7.109e-05 0.0001022 0.00% MLLinOp::define() 19 8.425e-05 9.024e-05 9.54e-05 0.00% Amr::FinalizeInit() 1 1.478e-05 2.512e-05 8.67e-05 0.00% Castro::expand_state() 20 5.761e-05 6.188e-05 7.201e-05 0.00% Castro::create_source_corrector() 17 5.629e-05 6.144e-05 7.085e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.753e-05 6.042e-05 6.584e-05 0.00% MLMG::computeMLResidual() 18 5.064e-05 5.349e-05 5.965e-05 0.00% Castro::derive() 18 4.746e-05 4.93e-05 5.567e-05 0.00% Amr::regrid() 1 1.609e-05 1.965e-05 4.38e-05 0.00% Gravity::actual_multilevel_solve() 1 3.2e-05 3.279e-05 3.678e-05 0.00% Amr::timeStep() 1 1.516e-05 1.796e-05 2.65e-05 0.00% Castro::post_init() 1 1.291e-05 1.363e-05 1.613e-05 0.00% Castro::advance() 1 8.136e-06 9.291e-06 1.557e-05 0.00% Castro::post_regrid() 2 1.091e-05 1.21e-05 1.452e-05 0.00% Amr::InitializeInit() 1 8.94e-06 9.693e-06 1.387e-05 0.00% Amr::init() 1 6.596e-06 7.775e-06 1.098e-05 0.00% Amr::bldFineLevels() 1 8.142e-06 9.34e-06 1.097e-05 0.00% Castro::post_timestep() 1 5.931e-06 6.34e-06 7.399e-06 0.00% Amr::initialInit() 1 3.111e-06 3.432e-06 3.809e-06 0.00% Other 2341 0.0484 0.05921 0.09127 0.68% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 13.32 13.33 13.33 100.00% Amr::coarseTimeStep() 1 11.18 11.18 11.18 83.87% Amr::timeStep() 1 11.13 11.13 11.13 83.50% Castro::advance() 1 11.08 11.08 11.08 83.11% Castro::subcycle_advance_ctu() 1 11.04 11.04 11.04 82.81% Castro::do_advance_ctu() 20 11 11 11 82.54% Castro::construct_ctu_hydro_source() 20 5.905 5.923 5.951 44.64% Castro::clean_state() 79 2.865 2.959 3.072 23.04% Castro::computeTemp() 80 1.998 2.09 2.206 16.55% Gravity::solve_phi_with_mlmg() 19 1.555 1.608 1.627 12.21% Castro::construct_new_gravity() 17 1.476 1.528 1.547 11.60% Gravity::solve_for_phi() 17 1.44 1.493 1.512 11.34% Gravity::actual_solve_with_mlmg() 19 1.36 1.364 1.369 10.27% Castro::reset_internal_energy(MultiFab) 80 1.29 1.296 1.314 9.86% MLMG::solve() 18 1.303 1.309 1.311 9.83% Amr::InitAmr() 1 1.257 1.261 1.267 9.50% MLMG::oneIter() 128 1.176 1.178 1.18 8.85% MLMG::mgVcycle() 128 1.162 1.164 1.166 8.75% Castro::initialize_do_advance() 20 1.093 1.123 1.15 8.63% Castro::do_old_sources() 40 1.075 1.099 1.121 8.41% FabArray::ParallelCopy() 1469 0.4855 0.7817 0.9442 7.08% Castro::evaluate_source_change() 196 0.9256 0.9295 0.9398 7.05% Amr::init() 1 0.882 0.882 0.882 6.62% Amr::initialInit() 1 0.831 0.831 0.8311 6.23% Castro::do_new_sources() 17 0.7546 0.7951 0.8273 6.21% MLCellLinOp::smooth() 1544 0.6019 0.6549 0.7913 5.94% Castro::enforce_min_density() 79 0.7692 0.7716 0.7791 5.84% amrex::Copy() 766 0.7614 0.7642 0.7764 5.82% FabArray::ParallelCopy_finish() 1469 0.2181 0.5119 0.6753 5.07% AmrLevel::FillPatch() 187 0.5584 0.6101 0.6602 4.95% Amr::InitializeInit() 1 0.4594 0.4882 0.6216 4.66% Amr::defBaseLevel() 1 0.4594 0.4882 0.6216 4.66% MultiFab::AddProduct() 3598 0.5358 0.5378 0.5441 4.08% FillPatchIterator::Initialize 187 0.368 0.421 0.4717 3.54% MLMG::mgVcycle_up::1 128 0.104 0.2768 0.452 3.39% MLCellLinOp::applyBC() 7148 0.3293 0.3628 0.4501 3.38% FillPatchIterator::FillFromLevel0() 187 0.3363 0.3957 0.4432 3.32% FillPatchSingleLevel 187 0.3354 0.3948 0.4422 3.32% MLPoisson::Fsmooth() 6144 0.3246 0.3502 0.422 3.17% Castro::Castro() 1 0.3918 0.3956 0.4126 3.09% Castro::initMFs() 1 0.3861 0.3899 0.4071 3.05% MLMG::addInterpCorrection() 768 0.06398 0.3182 0.3919 2.94% Amr::FinalizeInit() 1 0.2094 0.3428 0.3716 2.79% FabArray::setVal() 5232 0.3585 0.3613 0.3671 2.75% Castro::reset_internal_energy(Fab) 160 0.3455 0.3482 0.3518 2.64% MultiFab::Subtract() 159 0.3304 0.3321 0.3363 2.52% FabArray::FillBoundary() 6380 0.2383 0.2777 0.3184 2.39% MLMG::mgVcycle_up::2 128 0.01549 0.1275 0.2978 2.23% FabArray::ParallelCopy_nowait() 1469 0.2462 0.2655 0.2882 2.16% MLMG::mgVcycle_down::0 128 0.2557 0.258 0.2626 1.97% Gravity::fill_multipole_BCs() 19 0.1907 0.2422 0.2587 1.94% Castro::expand_state() 20 0.1873 0.2105 0.2354 1.77% Castro::initData() 1 0.06537 0.09086 0.2274 1.71% Amr::bldFineLevels() 1 0.03698 0.1706 0.1994 1.50% MLCellLinOp::apply() 986 0.1459 0.1585 0.1989 1.49% MLMG::mgVcycle_up::3 128 0.01526 0.02936 0.1978 1.48% MLMG::mgVcycle_up::0 128 0.1953 0.1962 0.197 1.48% check_for_negative_density() 20 0.03873 0.1463 0.1848 1.39% FabArray::sum() 3598 0.1793 0.18 0.1818 1.36% FillBoundary_nowait() 6380 0.1019 0.1322 0.1663 1.25% FillBoundary_finish() 6380 0.09922 0.1285 0.1653 1.24% Amr::grid_places() 2 0.002596 0.136 0.1648 1.24% AmrMesh::MakeNewGrids() 2 0.002526 0.1359 0.1647 1.24% TagBoxArray::mapPRD 2 0.000916 0.1094 0.1631 1.22% Castro::construct_old_source() 240 0.1589 0.1595 0.1625 1.22% MLMG:computeResOfCorrection() 768 0.114 0.1235 0.1458 1.09% MLCellLinOp::correctionResidual() 768 0.1122 0.1217 0.1439 1.08% Castro::post_init() 1 0.1367 0.1379 0.1384 1.04% TagBoxArray::collate() 2 0.0002555 0.02519 0.1379 1.03% MLMG::mgVcycle_down::1 128 0.1227 0.1267 0.1306 0.98% Castro::construct_new_source() 102 0.1182 0.1196 0.1227 0.92% Castro::normalize_species() 79 0.1192 0.1198 0.1216 0.91% Gravity::multilevel_solve_for_new_phi() 1 0.1046 0.1047 0.1047 0.79% Gravity::actual_multilevel_solve() 1 0.1046 0.1046 0.1047 0.79% MLPoisson::Fapply() 986 0.08494 0.08901 0.103 0.77% Castro::construct_old_gravity_source() 40 0.0965 0.09689 0.09849 0.74% MLCellLinOp::solutionResidual() 147 0.08208 0.08871 0.09745 0.73% Amr::writePlotFile() 2 0.09269 0.09278 0.0928 0.70% MLMG::mgVcycle_bottom 128 0.0004595 0.005543 0.08165 0.61% MLMG::actualBottomSolve() 128 0.000184 0.005262 0.08135 0.61% FabArray::Saxpy() 111 0.07945 0.07967 0.08133 0.61% MLCGSolver::bicgstab 8 0 0.00496 0.07936 0.60% MLMG::mgVcycle_down::2 128 0.01792 0.04697 0.07857 0.59% MLMG::computeResidual() 128 0.07178 0.07249 0.0756 0.57% Castro::post_regrid() 2 0.06833 0.0689 0.07011 0.53% Castro::construct_new_gravity_source() 17 0.06659 0.06765 0.06953 0.52% Castro::derive() 18 0.06558 0.06604 0.06682 0.50% AmrLevel::derive() 18 0.06553 0.06599 0.06677 0.50% MLMG::mgVcycle_down::3 128 0.01731 0.02373 0.06609 0.50% Castro::create_source_corrector() 17 0.04862 0.05366 0.06367 0.48% amrex::average_down 768 0.03377 0.0406 0.05697 0.43% FabArray::Xpay() 915 0.04872 0.05109 0.0564 0.42% Castro::construct_old_gravity() 20 0.04044 0.04857 0.05564 0.42% Gravity::get_old_grav_vector() 20 0.04031 0.04844 0.05551 0.42% Castro::post_timestep() 1 0.05101 0.05291 0.05532 0.41% Castro::check_for_nan() 37 0.04591 0.04629 0.04735 0.36% MultiFab::contains_nan() 37 0.04579 0.04616 0.04722 0.35% MLPoisson::define() 19 0.03263 0.03793 0.04316 0.32% Castro::initialize_advance() 1 0.04012 0.04055 0.04189 0.31% MLMG::mgVcycle_down::4 128 0.01733 0.0191 0.0406 0.30% Amr::regrid() 1 0.03552 0.03579 0.0369 0.28% StateDataPhysBCFunct::() 187 0.001018 0.01942 0.03643 0.27% StateData::FillBoundary(geom) 193 0 0.01824 0.0352 0.26% Gravity::get_new_grav_vector() 18 0.03405 0.03451 0.03512 0.26% Castro::retry_advance_ctu() 11 0.03379 0.034 0.03494 0.26% MLMG::mgVcycle_up::4 128 0.01535 0.01671 0.03269 0.25% MLMG::mgVcycle_down::5 128 0.01724 0.01842 0.02971 0.22% MLCellLinOp::defineAuxData() 19 0.02079 0.02245 0.02806 0.21% Castro::apply_source_to_state() 37 0.02652 0.02665 0.02706 0.20% MLMG::mgVcycle_up::5 128 0.01542 0.01643 0.02607 0.20% MLMG::prepareForSolve() 18 0.02018 0.02112 0.02402 0.18% VisMF::Write(FabArray) 4 0.01257 0.014 0.02381 0.18% MLMG::computeMLResidual() 18 0.01013 0.01582 0.02125 0.16% Castro::sum_integrated_quantities() 2 0.02075 0.02091 0.02093 0.16% MLMG::ResNormInf() 146 0.01125 0.01646 0.01866 0.14% Castro::save_data_for_retry() 3 0.01734 0.01753 0.018 0.13% FabArrayBase::getFB() 6380 0.01337 0.01426 0.01783 0.13% FabArray::setDomainBndry() 187 0.0006726 0.008923 0.01671 0.13% FabArrayBase::getCPC() 2045 0.01107 0.01239 0.01581 0.12% Castro::problem_post_timestep() 2 0.0105 0.01322 0.0157 0.12% Castro::wd_update() 2 0.01049 0.01319 0.01569 0.12% FabArray::norminf() 244 0.0102 0.01065 0.01406 0.11% amrex::Add() 128 0.0135 0.01356 0.01386 0.10% Gravity::test_composite_phi() 1 0.01362 0.01366 0.01369 0.10% MLLinOp::define() 19 0.004707 0.008216 0.0116 0.09% MLLinOp::defineGrids() 19 0.004623 0.008126 0.01151 0.09% MLCGSolver::ParallelAllReduce 189 0 0.0006099 0.009758 0.07% Other 2341 0.07921 0.09079 0.1255 0.94% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Async Memory Usage: -------------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max -------------------------------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 6400 407 KiB 409 KiB 417 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 2170 B 2187 B 2212 B 89 KiB 89 KiB 89 KiB -------------------------------------------------------------------------------------------------------------------- Pinned Memory Usage: -------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max -------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 4 B 5 B 7 B 8192 KiB 8192 KiB 8192 KiB -------------------------------------------------------------------------------------------------------------- Cpu Memory Usage: ----------------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max ----------------------------------------------------------------------------------------------------------------------- Castro::save_data_for_retry() 416 1585 KiB 1587 KiB 1587 KiB 2332 KiB 2332 KiB 2332 KiB Castro::construct_ctu_hydro_source() 6400 408 KiB 410 KiB 417 KiB 1178 KiB 1178 KiB 1178 KiB StateData::define() 192 1053 KiB 1054 KiB 1054 KiB 1166 KiB 1166 KiB 1166 KiB Castro::initMFs() 256 1032 KiB 1032 KiB 1033 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 673 KiB 674 KiB 674 KiB 804 KiB 804 KiB 804 KiB Amr::writePlotFile() 64 4476 B 4535 B 5067 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 51 KiB 51 KiB 52 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 13 KiB 14 KiB 15 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 434 KiB 434 KiB 434 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2343 B 4106 B 5129 B 226 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1427 B 1436 B 1455 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 19 KiB 19 KiB 19 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 131 KiB 131 KiB 131 KiB 158 KiB 158 KiB 158 KiB Castro::buildMetrics() 192 129 KiB 129 KiB 129 KiB 143 KiB 143 KiB 143 KiB MLMG::prepareForSolve() 6102 10 KiB 10 KiB 12 KiB 104 KiB 110 KiB 130 KiB Castro::construct_old_source() 1280 519 B 521 B 532 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 435 B 438 B 449 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 2175 B 2192 B 2218 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 18 B 21 B 39 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 33 KiB 34 KiB 34 KiB 76 KiB 76 KiB 76 KiB Gravity::test_composite_phi() 96 54 B 54 B 55 B 53 KiB 53 KiB 53 KiB Gravity::update_max_rhs() 192 3 B 4 B 6 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 116 B 117 B 119 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 138 B 166 B 191 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 32 KiB 32 KiB 32 KiB 38 KiB 38 KiB 38 KiB Gravity::fill_multipole_BCs() 2432 256 B 264 B 268 B 18 KiB 18 KiB 18 KiB MLCellLinOp::defineAuxData() 11704 826 B 993 B 1852 B 8192 B 9856 B 18 KiB Castro::evaluate_source_change() 6272 1131 B 1135 B 1148 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 128 B 128 B 128 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 1768 B 1832 B 1856 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 0 B 0 B 1 B 16 KiB 16 KiB 16 KiB BndryData::define() 4864 813 B 816 B 817 B 8192 B 8192 B 8192 B FillBoundary_nowait() 91296 31 B 40 B 53 B 3072 B 3840 B 4608 B AmrMesh::MakeNewGrids() 128 0 B 24 B 37 B 3584 B 3584 B 3584 B OwnerMask() 64 0 B 0 B 0 B 3072 B 3072 B 3072 B MLCGSolver::bicgstab 640 0 B 0 B 7 B 0 B 80 B 1280 B MLMG::addInterpCorrection() 5632 1 B 17 B 27 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 0 B 9 B 12 B 1024 B 1024 B 1024 B amrex::average_down 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B MLCellLinOp::applyBC() 114375 0 B 0 B 1 B 256 B 256 B 256 B ----------------------------------------------------------------------------------------------------------------------- AMReX (26.07-51-g50c73c8ed200) finalized