Initializing AMReX (26.08-21-g3b31af42a8b8)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.08-21-g3b31af42a8b8) initialized Starting run at 06:44:53 UTC on 2026-08-18. Successfully read inputs file ... Castro git describe: 26.08-1-gf0b335c8a AMReX git describe: 26.08-21-g3b31af42a Microphysics git describe: 26.08-10-g21581f4b 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.183717469 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001171585 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010242835 Gravity::multilevel_solve_for_new_phi() time = 0.10947677 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010190676 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.053190779 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004146092 seconds [Level 0 step 1] ADVANCE at time 0 with dt = 0.25  Beginning subcycle 1 starting at time 0 with dt = 0.25 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071148 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.769409 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477044 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035502 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -0.01295059363 at index 21, 83, 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.7670713549 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.7670713549 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Invalid density = -0.01295059363 at index 21, 172, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035502 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477044 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.769409 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071148 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Castro::construct_old_gravity() time = 0.004171263 on level 0 Castro::do_old_sources() time = 0.007451631 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.054260163 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.007543641 on level 0 Castro::do_old_sources() time = 0.007479778 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.05380757 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.004843193 on level 0 Castro::do_old_sources() time = 0.009264615 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.051075921 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.003577786 on level 0 Castro::do_old_sources() time = 0.007494464 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.048901931 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.298225152 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.013797082 Gravity::solve_for_phi() time = 0.089339002 on level 0 Castro::construct_new_gravity() time = 0.091386563 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.05076906 on level 0 Castro::construct_old_gravity() time = 0.002908074 on level 0 Castro::do_old_sources() time = 0.00744351 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.050044793 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.298328966 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.014060975 Gravity::solve_for_phi() time = 0.089176174 on level 0 Castro::construct_new_gravity() time = 0.091209419 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.049506587 on level 0 Castro::construct_old_gravity() time = 0.002662204 on level 0 Castro::do_old_sources() time = 0.007443636 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.048966933 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.303900924 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.013980948 Gravity::solve_for_phi() time = 0.089082451 on level 0 Castro::construct_new_gravity() time = 0.091134118 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.049555657 on level 0 Castro::construct_old_gravity() time = 0.003027607 on level 0 Castro::do_old_sources() time = 0.007458591 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.049042422 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.308265249 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.01437964 Gravity::solve_for_phi() time = 0.090165955 on level 0 Castro::construct_new_gravity() time = 0.092247559 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.049909291 on level 0 Castro::construct_old_gravity() time = 0.002870807 on level 0 Castro::do_old_sources() time = 0.007513295 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.049197658 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.312633551 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.014259113 Gravity::solve_for_phi() time = 0.08942659 on level 0 Castro::construct_new_gravity() time = 0.091478002 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.049480816 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.002755678 on level 0 Castro::do_old_sources() time = 0.00743759 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.049083145 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.308842417 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.014271988 Gravity::solve_for_phi() time = 0.089332376 on level 0 Castro::construct_new_gravity() time = 0.091352894 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.04959586 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.004501077 on level 0 Castro::do_old_sources() time = 0.007971351 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.050373923 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.314335884 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.01418659 Gravity::solve_for_phi() time = 0.089614823 on level 0 Castro::construct_new_gravity() time = 0.091627621 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.049597122 on level 0 Castro::construct_old_gravity() time = 0.002913506 on level 0 Castro::do_old_sources() time = 0.007457827 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.048812452 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.30802102 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.014347187 Gravity::solve_for_phi() time = 0.089644471 on level 0 Castro::construct_new_gravity() time = 0.091669197 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.049493201 on level 0 Castro::construct_old_gravity() time = 0.003362574 on level 0 Castro::do_old_sources() time = 0.007430562 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.0492788 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.313547012 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.014202574 Gravity::solve_for_phi() time = 0.090072766 on level 0 Castro::construct_new_gravity() time = 0.092153648 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.049472123 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.002888733 on level 0 Castro::do_old_sources() time = 0.007426771 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.049482492 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.308440045 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.014020733 Gravity::solve_for_phi() time = 0.089201984 on level 0 Castro::construct_new_gravity() time = 0.091286637 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.049500123 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.002763456 on level 0 Castro::do_old_sources() time = 0.007460071 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.048960565 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.313178347 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.014168449 Gravity::solve_for_phi() time = 0.089444641 on level 0 Castro::construct_new_gravity() time = 0.091465215 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.049814249 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.003467672 on level 0 Castro::do_old_sources() time = 0.007427739 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.048814333 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.308695026 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.014222069 Gravity::solve_for_phi() time = 0.089478466 on level 0 Castro::construct_new_gravity() time = 0.091493145 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.049527186 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.002677345 on level 0 Castro::do_old_sources() time = 0.007447527 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.048833227 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.314274851 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.014560537 Gravity::solve_for_phi() time = 0.08984997 on level 0 Castro::construct_new_gravity() time = 0.091867884 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.049705948 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.003384978 on level 0 Castro::do_old_sources() time = 0.007660354 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.049243785 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.309628508 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.014411976 Gravity::solve_for_phi() time = 0.08990471 on level 0 Castro::construct_new_gravity() time = 0.09199823 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.0495818 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.003089794 on level 0 Castro::do_old_sources() time = 0.00762596 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.048917928 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.314273358 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.014052851 Gravity::solve_for_phi() time = 0.089458298 on level 0 Castro::construct_new_gravity() time = 0.091524489 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.049595726 on level 0 Castro::construct_old_gravity() time = 0.002705604 on level 0 Castro::do_old_sources() time = 0.007462596 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.049184829 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.30923639 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.014518821 Gravity::solve_for_phi() time = 0.089777009 on level 0 Castro::construct_new_gravity() time = 0.091803138 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.049614481 on level 0 Castro::construct_old_gravity() time = 0.002711704 on level 0 Castro::do_old_sources() time = 0.007452199 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.049066052 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.313519977 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002908291199 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.014590403 Gravity::solve_for_phi() time = 0.089930945 on level 0 Castro::construct_new_gravity() time = 0.092006238 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.04983843 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.31914876 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.049237209 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004242109 seconds Ending run at 06:45:07 UTC on 2026-08-18. Run time = 14.10016028 Run time without initialization = 11.37354502 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.0784 ... 14.0912 ... 14.1002 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.1571 5.2381 5.3167 37.71% Amr::InitAmr() 1 1.79 1.804 1.813 12.86% Castro::computeTemp() 80 0.7029 0.8031 0.9185 6.51% amrex::Copy() 766 0.7797 0.793 0.8304 5.89% FabArray::ParallelCopy_finish() 1469 0.3108 0.5452 0.7682 5.45% MultiFab::AddProduct() 3598 0.5483 0.556 0.5639 4.00% MLPoisson::Fsmooth() 6144 0.3297 0.3589 0.4259 3.02% Castro::initMFs() 1 0.3806 0.3877 0.3943 2.80% FabArray::setVal() 5232 0.3562 0.363 0.3674 2.61% Castro::reset_internal_energy(Fab) 160 0.3474 0.353 0.3583 2.54% MultiFab::Subtract() 159 0.3313 0.3358 0.3407 2.42% FabArray::ParallelCopy_nowait() 1469 0.2311 0.2523 0.2713 1.92% Castro::do_advance_ctu() 20 0.2581 0.2642 0.2703 1.92% Gravity::fill_multipole_BCs() 19 0.19 0.2366 0.2597 1.84% check_for_negative_density() 20 0.03737 0.1303 0.2416 1.71% FabArray::sum() 3598 0.191 0.194 0.1966 1.39% Castro::initData() 1 0.003157 0.0288 0.1718 1.22% FillBoundary_finish() 6380 0.1047 0.1278 0.1713 1.22% TagBoxArray::collate() 2 0.0002059 0.03368 0.1684 1.19% FillBoundary_nowait() 6380 0.09432 0.1202 0.1511 1.07% MLCellLinOp::applyBC() 7148 0.07153 0.08657 0.1318 0.93% Castro::normalize_species() 79 0.1199 0.1215 0.123 0.87% MLPoisson::Fapply() 986 0.08554 0.08998 0.1035 0.73% Castro::construct_old_gravity_source() 40 0.09472 0.09628 0.09887 0.70% FabArray::Saxpy() 111 0.07587 0.07688 0.07786 0.55% Castro::construct_new_gravity_source() 17 0.06709 0.06826 0.06958 0.49% FabArray::Xpay() 915 0.04921 0.05164 0.0563 0.40% MultiFab::contains_nan() 37 0.04373 0.04435 0.04485 0.32% AmrLevel::derive() 18 0.03898 0.03989 0.04333 0.31% MLMG::addInterpCorrection() 768 0.03283 0.03581 0.04271 0.30% StateData::FillBoundary(geom) 193 0 0.01831 0.03514 0.25% VisMF::Write(FabArray) 4 0.0124 0.01426 0.03338 0.24% amrex::average_down 768 0.02423 0.02659 0.03295 0.23% Castro::reset_internal_energy(MultiFab) 80 0.02526 0.02573 0.0263 0.19% Amr::writePlotFile() 2 0.006788 0.02282 0.02577 0.18% MLCellLinOp::defineAuxData() 19 0.01626 0.01719 0.02066 0.15% Castro::construct_old_source() 240 0.01931 0.01973 0.02016 0.14% FabArray::FillBoundary() 6380 0.01665 0.01725 0.02012 0.14% Castro::enforce_min_density() 79 0.01926 0.01962 0.02001 0.14% FillPatchIterator::Initialize 187 0.01603 0.01647 0.01691 0.12% FabArrayBase::getFB() 6380 0.01352 0.01411 0.01684 0.12% FabArray::setDomainBndry() 187 0.0006655 0.009027 0.01673 0.12% MLCGSolver::bicgstab 8 0 0.001044 0.01671 0.12% Castro::evaluate_source_change() 196 0.01446 0.01484 0.01592 0.11% Castro::construct_new_source() 102 0.01511 0.01537 0.01556 0.11% FabArray::norminf() 244 0.01067 0.01118 0.01453 0.10% amrex::Add() 128 0.01337 0.01359 0.01399 0.10% MLCellLinOp::smooth() 1544 0.01063 0.01118 0.01269 0.09% MLMG::prepareForSolve() 18 0.01014 0.01062 0.0115 0.08% Castro::wd_update() 2 0.004251 0.008073 0.01148 0.08% MLMG::ResNormInf() 146 0.00324 0.006776 0.01033 0.07% FabArrayBase::getCPC() 2045 0.008684 0.008887 0.009496 0.07% Gravity::actual_solve_with_mlmg() 19 0.007432 0.007872 0.009431 0.07% MLMG::solve() 18 0.000797 0.003697 0.006963 0.05% Gravity::get_old_grav_vector() 20 0.006144 0.006263 0.006397 0.05% Gravity::get_new_grav_vector() 18 0.005755 0.005853 0.005965 0.04% MLCellLinOp::apply() 986 0.002356 0.002613 0.005476 0.04% Castro::initialize_do_advance() 20 0.004727 0.004872 0.005164 0.04% Castro::sum_integrated_quantities() 2 0.00234 0.004656 0.005142 0.04% FabArray::ParallelCopy() 1469 0.004111 0.004284 0.004595 0.03% AmrLevel::FillPatch() 187 0.003129 0.003256 0.003636 0.03% MLLinOp::defineGrids() 19 0.002612 0.00275 0.003122 0.02% MLMG::mgVcycle() 128 0.002256 0.002304 0.002396 0.02% MLCellLinOp::correctionResidual() 768 0.002262 0.002323 0.002358 0.02% MLMG:computeResOfCorrection() 768 0.00186 0.001894 0.002006 0.01% Castro::save_data_for_retry() 3 0.001333 0.001366 0.001442 0.01% StateDataPhysBCFunct::() 187 0.0009949 0.001177 0.001375 0.01% MLMG::actualBottomSolve() 128 0.0001872 0.0002577 0.00125 0.01% main() 1 0.0007922 0.0008849 0.001129 0.01% FillPatchIterator::FillFromLevel0() 187 0.0008938 0.0009292 0.000993 0.01% Castro::subcycle_advance_ctu() 1 0.0002396 0.00032 0.0009247 0.01% Gravity::solve_for_phi() 17 0.0006892 0.0007247 0.0008266 0.01% FillPatchSingleLevel 187 0.0006663 0.0006819 0.000708 0.01% MLMG::mgVcycle_down::0 128 0.0006298 0.0006445 0.0006783 0.00% MLMG::mgVcycle_down::1 128 0.0006107 0.000629 0.0006698 0.00% MLMG::mgVcycle_down::2 128 0.000593 0.0006338 0.0006616 0.00% MLMG::mgVcycle_down::5 128 0.0005883 0.0006133 0.0006587 0.00% MLMG::mgVcycle_down::4 128 0.0005628 0.000588 0.000648 0.00% MLMG::mgVcycle_down::3 128 0.0005703 0.0006035 0.0006424 0.00% AmrMesh::MakeNewGrids() 2 0.0005458 0.0005826 0.0006139 0.00% Castro::clean_state() 79 0.000482 0.0005052 0.0005456 0.00% MLCellLinOp::solutionResidual() 147 0.0005069 0.000525 0.0005403 0.00% Castro::do_old_sources() 40 0.0004526 0.0004723 0.0005029 0.00% Castro::Castro() 1 0.0004187 0.0004392 0.000482 0.00% MLMG::mgVcycle_up::2 128 0.0003942 0.0004228 0.0004814 0.00% MLMG::oneIter() 128 0.00043 0.0004474 0.0004777 0.00% MLMG::mgVcycle_up::1 128 0.0004253 0.000435 0.0004635 0.00% MLMG::mgVcycle_up::0 128 0.0004144 0.0004255 0.000458 0.00% MLMG::mgVcycle_up::4 128 0.0003859 0.0004009 0.0004535 0.00% MLMG::mgVcycle_up::3 128 0.0003834 0.0003987 0.0004444 0.00% MLMG::mgVcycle_up::5 128 0.0003721 0.0003829 0.000398 0.00% MLMG::computeResidual() 128 0.0002897 0.0003016 0.0003223 0.00% Castro::construct_new_gravity() 17 0.0001851 0.0002072 0.0003117 0.00% MLMG::mgVcycle_bottom 128 0.0002706 0.0002808 0.0003063 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002196 0.0002364 0.0002638 0.00% Castro::do_new_sources() 17 0.0002143 0.000227 0.0002385 0.00% Castro::problem_post_timestep() 2 1.207e-05 2.793e-05 0.0002308 0.00% Amr::defBaseLevel() 1 0.0001084 0.0001407 0.0001586 0.00% Castro::construct_old_gravity() 20 0.0001289 0.0001381 0.0001538 0.00% MLPoisson::define() 19 0.0001285 0.0001353 0.0001528 0.00% Amr::coarseTimeStep() 1 0.0001158 0.0001426 0.0001502 0.00% Castro::retry_advance_ctu() 11 7.147e-05 8.569e-05 0.0001467 0.00% Castro::apply_source_to_state() 37 0.0001153 0.0001243 0.0001454 0.00% TagBoxArray::mapPRD 2 0.0001228 0.0001289 0.0001413 0.00% Castro::check_for_nan() 37 0.000121 0.0001293 0.000141 0.00% Gravity::test_composite_phi() 1 0.0001067 0.0001137 0.0001249 0.00% Castro::initialize_advance() 1 0.0001053 0.0001098 0.0001237 0.00% MLLinOp::define() 19 8.39e-05 8.983e-05 9.73e-05 0.00% Amr::grid_places() 2 5.615e-05 6.053e-05 8.688e-05 0.00% Castro::expand_state() 20 5.806e-05 6.259e-05 7.552e-05 0.00% Castro::create_source_corrector() 17 5.821e-05 6.334e-05 7.139e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.741e-05 6.006e-05 6.435e-05 0.00% Castro::derive() 18 4.844e-05 5.177e-05 5.894e-05 0.00% MLMG::computeMLResidual() 18 5.052e-05 5.361e-05 5.755e-05 0.00% Amr::FinalizeInit() 1 8.951e-06 1.246e-05 4.971e-05 0.00% Amr::regrid() 1 1.564e-05 1.942e-05 4.034e-05 0.00% Gravity::actual_multilevel_solve() 1 3.223e-05 3.333e-05 3.641e-05 0.00% Castro::post_init() 1 1.598e-05 1.893e-05 2.943e-05 0.00% Amr::timeStep() 1 1.57e-05 1.892e-05 2.288e-05 0.00% Castro::post_regrid() 2 1.133e-05 1.255e-05 1.759e-05 0.00% Castro::advance() 1 8.169e-06 1.051e-05 1.74e-05 0.00% Amr::InitializeInit() 1 8.43e-06 9.427e-06 1.401e-05 0.00% Amr::bldFineLevels() 1 7.727e-06 8.807e-06 1.234e-05 0.00% Castro::post_timestep() 1 6.201e-06 7.123e-06 8.577e-06 0.00% Amr::init() 1 6.43e-06 7.445e-06 8.453e-06 0.00% Amr::initialInit() 1 3.295e-06 3.504e-06 4.181e-06 0.00% Other 2530 0.04778 0.06096 0.09996 0.71% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.08 14.09 14.1 100.00% Amr::coarseTimeStep() 1 11.37 11.37 11.37 80.66% Amr::timeStep() 1 11.32 11.32 11.32 80.28% Castro::advance() 1 11.26 11.27 11.27 79.91% Castro::subcycle_advance_ctu() 1 11.22 11.22 11.23 79.62% Castro::do_advance_ctu() 20 11.18 11.19 11.19 79.37% Castro::construct_ctu_hydro_source() 20 5.983 6.004 6.037 42.82% Castro::clean_state() 79 2.899 3.027 3.137 22.25% Castro::computeTemp() 80 2.016 2.133 2.234 15.85% Amr::InitAmr() 1 1.791 1.804 1.813 12.86% Gravity::solve_phi_with_mlmg() 19 1.57 1.617 1.643 11.65% Castro::construct_new_gravity() 17 1.485 1.533 1.558 11.05% Gravity::solve_for_phi() 17 1.45 1.497 1.523 10.80% Gravity::actual_solve_with_mlmg() 19 1.373 1.377 1.381 9.80% Castro::reset_internal_energy(MultiFab) 80 1.312 1.33 1.349 9.57% MLMG::solve() 18 1.317 1.322 1.324 9.39% MLMG::oneIter() 128 1.187 1.19 1.193 8.46% Castro::initialize_do_advance() 20 1.137 1.164 1.191 8.45% MLMG::mgVcycle() 128 1.173 1.175 1.179 8.36% Castro::do_old_sources() 40 1.088 1.123 1.137 8.07% FabArray::ParallelCopy() 1469 0.5829 0.8223 1.041 7.39% Castro::evaluate_source_change() 196 0.9501 0.9639 0.9774 6.93% Amr::init() 1 0.9138 0.9138 0.9138 6.48% Amr::initialInit() 1 0.8559 0.856 0.856 6.07% Castro::do_new_sources() 17 0.7615 0.8096 0.844 5.99% amrex::Copy() 766 0.7797 0.793 0.8304 5.89% Castro::enforce_min_density() 79 0.7837 0.796 0.8157 5.78% MLCellLinOp::smooth() 1544 0.6113 0.666 0.7985 5.66% FabArray::ParallelCopy_finish() 1469 0.3108 0.5452 0.7682 5.45% AmrLevel::FillPatch() 187 0.5901 0.654 0.7353 5.21% Amr::InitializeInit() 1 0.4538 0.4875 0.6361 4.51% Amr::defBaseLevel() 1 0.4537 0.4875 0.6361 4.51% MultiFab::AddProduct() 3598 0.5483 0.556 0.5639 4.00% FillPatchIterator::Initialize 187 0.393 0.4584 0.5429 3.85% FillPatchIterator::FillFromLevel0() 187 0.3649 0.4329 0.5096 3.61% FillPatchSingleLevel 187 0.364 0.432 0.5087 3.61% MLCellLinOp::applyBC() 7148 0.3308 0.3666 0.4564 3.24% MLMG::mgVcycle_up::1 128 0.1038 0.2752 0.4495 3.19% MLPoisson::Fsmooth() 6144 0.3297 0.3589 0.4259 3.02% Amr::FinalizeInit() 1 0.2199 0.3685 0.4022 2.85% Castro::Castro() 1 0.3858 0.393 0.4 2.84% Castro::initMFs() 1 0.3808 0.3878 0.3945 2.80% MLMG::addInterpCorrection() 768 0.06584 0.3163 0.3908 2.77% FabArray::setVal() 5232 0.3562 0.363 0.3674 2.61% Castro::reset_internal_energy(Fab) 160 0.3474 0.353 0.3583 2.54% MultiFab::Subtract() 159 0.3313 0.3358 0.3407 2.42% FabArray::FillBoundary() 6380 0.2441 0.28 0.3246 2.30% MLMG::mgVcycle_up::2 128 0.01563 0.1268 0.2965 2.10% FabArray::ParallelCopy_nowait() 1469 0.2513 0.2729 0.2926 2.07% MLMG::mgVcycle_down::0 128 0.2609 0.2639 0.2674 1.90% Gravity::fill_multipole_BCs() 19 0.192 0.2386 0.2617 1.86% Castro::expand_state() 20 0.2014 0.223 0.2479 1.76% check_for_negative_density() 20 0.03737 0.1303 0.2416 1.71% Castro::initData() 1 0.06588 0.09277 0.2363 1.68% Amr::bldFineLevels() 1 0.03826 0.1868 0.2199 1.56% MLCellLinOp::apply() 986 0.1479 0.1611 0.2027 1.44% MLMG::mgVcycle_up::0 128 0.2006 0.2011 0.2015 1.43% FabArray::sum() 3598 0.191 0.194 0.1966 1.39% MLMG::mgVcycle_up::3 128 0.01556 0.0295 0.1959 1.39% Amr::grid_places() 2 0.002526 0.1511 0.1849 1.31% AmrMesh::MakeNewGrids() 2 0.002469 0.1511 0.1848 1.31% TagBoxArray::mapPRD 2 0.0007931 0.116 0.1825 1.29% FillBoundary_finish() 6380 0.1047 0.1278 0.1713 1.22% TagBoxArray::collate() 2 0.0002059 0.03368 0.1684 1.19% FillBoundary_nowait() 6380 0.1095 0.135 0.1662 1.18% Castro::construct_old_source() 240 0.1562 0.1585 0.1619 1.15% Castro::post_init() 1 0.1442 0.1468 0.1481 1.05% MLMG:computeResOfCorrection() 768 0.1166 0.1254 0.1474 1.05% MLCellLinOp::correctionResidual() 768 0.1147 0.1235 0.1456 1.03% MLMG::mgVcycle_down::1 128 0.1238 0.1275 0.1314 0.93% Castro::normalize_species() 79 0.1199 0.1215 0.123 0.87% Castro::construct_new_source() 102 0.118 0.1197 0.1217 0.86% Gravity::multilevel_solve_for_new_phi() 1 0.1094 0.1094 0.1095 0.78% Gravity::actual_multilevel_solve() 1 0.1094 0.1094 0.1094 0.78% MLPoisson::Fapply() 986 0.08554 0.08998 0.1035 0.73% Amr::writePlotFile() 2 0.1028 0.1029 0.1029 0.73% Castro::construct_old_gravity_source() 40 0.09472 0.09628 0.09887 0.70% MLCellLinOp::solutionResidual() 147 0.08502 0.09011 0.0963 0.68% Castro::create_source_corrector() 17 0.05272 0.06432 0.08397 0.60% MLMG::mgVcycle_bottom 128 0.0004578 0.005533 0.08146 0.58% MLMG::actualBottomSolve() 128 0.0001872 0.005252 0.08116 0.58% MLMG::mgVcycle_down::2 128 0.01819 0.0475 0.07956 0.56% MLCGSolver::bicgstab 8 0 0.004947 0.07916 0.56% FabArray::Saxpy() 111 0.07587 0.07688 0.07786 0.55% MLMG::computeResidual() 128 0.07334 0.07439 0.07499 0.53% Castro::derive() 18 0.06666 0.06837 0.07349 0.52% AmrLevel::derive() 18 0.06662 0.06832 0.07344 0.52% Castro::post_regrid() 2 0.06917 0.07053 0.07309 0.52% Castro::construct_new_gravity_source() 17 0.06709 0.06826 0.06958 0.49% Castro::construct_old_gravity() 20 0.04341 0.05411 0.06614 0.47% Gravity::get_old_grav_vector() 20 0.04328 0.05397 0.06599 0.47% MLMG::mgVcycle_down::3 128 0.0177 0.02406 0.06555 0.46% amrex::average_down 768 0.03398 0.04098 0.057 0.40% FabArray::Xpay() 915 0.04921 0.05164 0.0563 0.40% Castro::post_timestep() 1 0.0513 0.05317 0.05586 0.40% Castro::initialize_advance() 1 0.04021 0.04131 0.04804 0.34% Castro::check_for_nan() 37 0.04386 0.04448 0.04498 0.32% MultiFab::contains_nan() 37 0.04373 0.04435 0.04485 0.32% MLPoisson::define() 19 0.03396 0.03804 0.04331 0.31% MLMG::mgVcycle_down::4 128 0.01767 0.01954 0.04094 0.29% Castro::retry_advance_ctu() 11 0.03445 0.03554 0.04019 0.29% Amr::regrid() 1 0.03601 0.03677 0.03849 0.27% StateDataPhysBCFunct::() 187 0.0009949 0.01949 0.03652 0.26% Gravity::get_new_grav_vector() 18 0.03443 0.035 0.0355 0.25% StateData::FillBoundary(geom) 193 0 0.01831 0.03514 0.25% VisMF::Write(FabArray) 4 0.0124 0.01426 0.03338 0.24% MLMG::mgVcycle_up::4 128 0.01562 0.01705 0.03245 0.23% MLMG::mgVcycle_down::5 128 0.01767 0.01886 0.03003 0.21% MLCellLinOp::defineAuxData() 19 0.02097 0.02261 0.02791 0.20% Castro::apply_source_to_state() 37 0.02546 0.02581 0.02612 0.19% MLMG::mgVcycle_up::5 128 0.01566 0.01667 0.02574 0.18% Castro::sum_integrated_quantities() 2 0.02418 0.02431 0.02432 0.17% MLMG::prepareForSolve() 18 0.02025 0.02133 0.02383 0.17% Castro::save_data_for_retry() 3 0.01778 0.01862 0.02302 0.16% MLMG::computeMLResidual() 18 0.01023 0.01534 0.02071 0.15% MLMG::ResNormInf() 146 0.01288 0.01631 0.01978 0.14% FabArrayBase::getFB() 6380 0.01386 0.0148 0.01831 0.13% FabArray::setDomainBndry() 187 0.0006655 0.009027 0.01673 0.12% FabArrayBase::getCPC() 2045 0.01115 0.01243 0.01569 0.11% Castro::problem_post_timestep() 2 0.009883 0.01183 0.01486 0.11% Castro::wd_update() 2 0.009867 0.0118 0.01484 0.11% FabArray::norminf() 244 0.01067 0.01118 0.01453 0.10% amrex::Add() 128 0.01337 0.01359 0.01399 0.10% Gravity::test_composite_phi() 1 0.01386 0.0139 0.01392 0.10% MLLinOp::define() 19 0.004297 0.008166 0.01159 0.08% MLLinOp::defineGrids() 19 0.004206 0.008076 0.01149 0.08% Other 2530 0.07885 0.09468 0.1404 1.00% -------------------------------------------------------------------------------------------- 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 387 KiB 393 KiB 400 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 2101 B 2137 B 2174 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 1 B 2 B 4 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 1523 KiB 1524 KiB 1525 KiB 2332 KiB 2332 KiB 2332 KiB Castro::construct_ctu_hydro_source() 6400 388 KiB 394 KiB 400 KiB 1178 KiB 1178 KiB 1178 KiB StateData::define() 192 1014 KiB 1015 KiB 1016 KiB 1166 KiB 1166 KiB 1166 KiB Castro::initMFs() 256 994 KiB 994 KiB 995 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 647 KiB 648 KiB 648 KiB 804 KiB 804 KiB 804 KiB Amr::writePlotFile() 64 4264 B 4443 B 5357 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 49 KiB 50 KiB 51 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 417 KiB 417 KiB 417 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2556 B 4216 B 5821 B 226 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1376 B 1410 B 1527 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 18 KiB 18 KiB 19 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 126 KiB 126 KiB 126 KiB 158 KiB 158 KiB 158 KiB Castro::buildMetrics() 192 124 KiB 124 KiB 124 KiB 143 KiB 143 KiB 143 KiB MLMG::prepareForSolve() 6102 9952 B 10 KiB 12 KiB 104 KiB 110 KiB 130 KiB Castro::construct_old_source() 1280 482 B 489 B 497 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 405 B 410 B 417 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 2106 B 2142 B 2180 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 19 B 20 B 38 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 32 KiB 32 KiB 32 KiB 76 KiB 76 KiB 76 KiB Gravity::test_composite_phi() 96 52 B 52 B 52 B 53 KiB 53 KiB 53 KiB Gravity::update_max_rhs() 192 3 B 6 B 8 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 110 B 112 B 114 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 141 B 176 B 216 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 31 KiB 31 KiB 31 KiB 38 KiB 38 KiB 38 KiB Gravity::fill_multipole_BCs() 2432 244 B 251 B 258 B 18 KiB 18 KiB 18 KiB MLCellLinOp::defineAuxData() 11704 788 B 948 B 1768 B 8192 B 9856 B 18 KiB Castro::evaluate_source_change() 6272 1097 B 1113 B 1129 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 126 B 126 B 127 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 1683 B 1737 B 1768 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 776 B 779 B 780 B 8192 B 8192 B 8192 B FillBoundary_nowait() 91296 27 B 39 B 51 B 3072 B 3840 B 4608 B AmrMesh::MakeNewGrids() 128 0 B 25 B 40 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 6 B 0 B 80 B 1280 B MLMG::addInterpCorrection() 5632 1 B 16 B 25 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 0 B 10 B 13 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.08-21-g3b31af42a8b8) finalized