Initializing AMReX (26.01-17-g909028cb95b3)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.01-17-g909028cb95b3) initialized Starting run at 03:43:58 UTC on 2026-01-16. Successfully read inputs file ... Castro git describe: 26.01-10-g409d8e9cc AMReX git describe: 26.01-17-g909028cb9 Microphysics git describe: 26.01-21-gbace8d29 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.166976124 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.002804672 new finest: 0 Gravity::fill_multipole_BCs() time = 0.033697329 Gravity::multilevel_solve_for_new_phi() time = 0.326608215 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.033629496 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.111140739 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.00685837 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 = -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 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 Castro::construct_old_gravity() time = 0.009981295 on level 0 Castro::do_old_sources() time = 0.026583412 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.148831255 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.015565681 on level 0 Castro::do_old_sources() time = 0.027589274 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.147094075 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.125 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.0625   Beginning subcycle 1 starting at time 0 with dt = 0.0625 Estimated number of subcycles remaining: 4 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.1060744788 at index 20, 177, 0 Invalid density = -0.1060744788 at index 20, 78, 0 Castro::construct_old_gravity() time = 0.044667554 on level 0 Castro::do_old_sources() time = 0.041934295 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.149456177 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.008979146 on level 0 Castro::do_old_sources() time = 0.026237967 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.143619067 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.809444228 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.036012935 Gravity::solve_for_phi() time = 0.270270583 on level 0 Castro::construct_new_gravity() time = 0.277849826 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.149021641 on level 0 Castro::construct_old_gravity() time = 0.009684932 on level 0 Castro::do_old_sources() time = 0.026699724 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.144671837 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.811222015 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.035830341 Gravity::solve_for_phi() time = 0.267287544 on level 0 Castro::construct_new_gravity() time = 0.274789191 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.147893312 on level 0 Castro::construct_old_gravity() time = 0.010406163 on level 0 Castro::do_old_sources() time = 0.026680056 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.144583932 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.842849856 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.035521262 Gravity::solve_for_phi() time = 0.26499338 on level 0 Castro::construct_new_gravity() time = 0.272541331 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.148197957 on level 0 Castro::construct_old_gravity() time = 0.009460731 on level 0 Castro::do_old_sources() time = 0.026932468 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.144914058 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.819161331 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.035854488 Gravity::solve_for_phi() time = 0.275088031 on level 0 Castro::construct_new_gravity() time = 0.286445609 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.155548716 on level 0 Castro::construct_old_gravity() time = 0.009593573 on level 0 Castro::do_old_sources() time = 0.026982879 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.143715412 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.836554817 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.035459411 Gravity::solve_for_phi() time = 0.268492462 on level 0 Castro::construct_new_gravity() time = 0.276012327 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.147090185 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.009842402 on level 0 Castro::do_old_sources() time = 0.026290051 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.143554747 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.875854157 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.036642637 Gravity::solve_for_phi() time = 0.268116723 on level 0 Castro::construct_new_gravity() time = 0.275646978 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.147775369 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.009767652 on level 0 Castro::do_old_sources() time = 0.027343203 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.145145589 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.842145579 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.035535508 Gravity::solve_for_phi() time = 0.26766599 on level 0 Castro::construct_new_gravity() time = 0.275186313 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.148174533 on level 0 Castro::construct_old_gravity() time = 0.009793903 on level 0 Castro::do_old_sources() time = 0.026293745 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.14481021 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.824684938 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.035869627 Gravity::solve_for_phi() time = 0.267271853 on level 0 Castro::construct_new_gravity() time = 0.275060465 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.148907749 on level 0 Castro::construct_old_gravity() time = 0.009930869 on level 0 Castro::do_old_sources() time = 0.026513457 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.145476194 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.864916729 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.040494373 Gravity::solve_for_phi() time = 0.273863255 on level 0 Castro::construct_new_gravity() time = 0.281435707 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.147541579 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.009848357 on level 0 Castro::do_old_sources() time = 0.026612385 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.144641812 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.820317999 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.036142776 Gravity::solve_for_phi() time = 0.268925389 on level 0 Castro::construct_new_gravity() time = 0.276424491 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.148755076 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.0096992 on level 0 Castro::do_old_sources() time = 0.026738399 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.144570861 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.849094786 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.046155352 Gravity::solve_for_phi() time = 0.294529207 on level 0 Castro::construct_new_gravity() time = 0.302761612 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.15824507 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.016005641 on level 0 Castro::do_old_sources() time = 0.02808001 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.154294537 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.877163343 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.045615569 Gravity::solve_for_phi() time = 0.281626546 on level 0 Castro::construct_new_gravity() time = 0.289119152 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.15000955 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.009933986 on level 0 Castro::do_old_sources() time = 0.030417249 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.148310675 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.852191137 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.039757558 Gravity::solve_for_phi() time = 0.273308947 on level 0 Castro::construct_new_gravity() time = 0.280834858 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.14973033 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.00987529 on level 0 Castro::do_old_sources() time = 0.026838216 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.146286573 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.828938922 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.039233684 Gravity::solve_for_phi() time = 0.271789454 on level 0 Castro::construct_new_gravity() time = 0.279308441 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.1507706 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.009794916 on level 0 Castro::do_old_sources() time = 0.02682445 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.146249597 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.847995909 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.039723407 Gravity::solve_for_phi() time = 0.279525914 on level 0 Castro::construct_new_gravity() time = 0.288465253 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.184665215 on level 0 Castro::construct_old_gravity() time = 0.017426144 on level 0 Castro::do_old_sources() time = 0.061616565 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.149317862 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.844731721 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.039491237 Gravity::solve_for_phi() time = 0.271974919 on level 0 Castro::construct_new_gravity() time = 0.279464719 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.150277418 on level 0 Castro::construct_old_gravity() time = 0.009793444 on level 0 Castro::do_old_sources() time = 0.026817078 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.14640033 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.849640552 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.001017328664 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.040379173 Gravity::solve_for_phi() time = 0.276067677 on level 0 Castro::construct_new_gravity() time = 0.283624421 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.151396766 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: 32.36936486 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.112266869 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.006910297 seconds Ending run at 03:44:34 UTC on 2026-01-16. Run time = 35.47634482 Run time without initialization = 32.4894538 Average number of zones advanced per microsecond: 0.001 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 35.4703 ... 35.4721 ... 35.4764 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 13.8230 13.9912 14.1542 39.90% amrex::Copy() 766 3.253 3.289 3.329 9.38% MultiFab::AddProduct() 3598 2.201 2.222 2.251 6.34% MLPoisson::Fsmooth() 6144 1.469 1.593 1.902 5.36% FabArray::ParallelCopy_finish() 1469 0.5578 1.335 1.75 4.93% Amr::InitAmr() 1 1.462 1.464 1.468 4.14% FabArray::setVal() 5232 1.36 1.382 1.409 3.97% MultiFab::Subtract() 159 1.37 1.383 1.399 3.94% Castro::computeTemp() 80 0.9805 1.062 1.178 3.32% FabArray::ParallelCopy_nowait() 1469 0.7366 0.8023 0.8734 2.46% FabArray::sum() 3598 0.7475 0.7548 0.7637 2.15% Gravity::fill_multipole_BCs() 19 0.6521 0.68 0.7108 2.00% check_for_negative_density() 20 0.2975 0.4889 0.6706 1.89% MLPoisson::Fapply() 986 0.5587 0.5865 0.6541 1.84% Castro::do_advance_ctu() 20 0.5174 0.5242 0.5309 1.50% Castro::initMFs() 1 0.38 0.3874 0.4284 1.21% Castro::reset_internal_energy(Fab) 160 0.3682 0.372 0.3759 1.06% Castro::construct_old_gravity_source() 40 0.3228 0.3283 0.348 0.98% FabArray::Saxpy() 111 0.3174 0.3203 0.324 0.91% FillBoundary_finish() 6380 0.2243 0.263 0.3204 0.90% MLCellLinOp::applyBC() 7148 0.1249 0.1738 0.2946 0.83% Castro::construct_new_gravity_source() 17 0.2646 0.2679 0.2713 0.76% FillBoundary_nowait() 6380 0.1391 0.1861 0.2421 0.68% FabArray::Xpay() 915 0.1903 0.1994 0.2185 0.62% MultiFab::contains_nan() 37 0.1648 0.1672 0.17 0.48% Castro::initData() 1 0.01027 0.03484 0.1698 0.48% Castro::normalize_species() 79 0.1388 0.1406 0.1455 0.41% TagBoxArray::collate() 2 0.000546 0.02327 0.113 0.32% MLMG::addInterpCorrection() 768 0.08261 0.08921 0.1058 0.30% StateData::FillBoundary(geom) 193 0 0.05114 0.1009 0.28% amrex::average_down 768 0.05895 0.06421 0.07738 0.22% Castro::construct_old_source() 240 0.07072 0.07144 0.07273 0.21% Castro::construct_new_source() 102 0.05882 0.05942 0.061 0.17% AmrLevel::derive() 18 0.0588 0.05973 0.06063 0.17% amrex::Add() 128 0.0524 0.05314 0.05399 0.15% VisMF::Write(FabArray) 4 0.03735 0.03852 0.05269 0.15% FabArray::setDomainBndry() 187 0.0006601 0.02525 0.04927 0.14% FabArray::norminf() 244 0.03587 0.03681 0.04106 0.12% Castro::enforce_min_density() 79 0.03312 0.03354 0.034 0.10% Gravity::get_old_grav_vector() 20 0.03008 0.03056 0.03118 0.09% Castro::volWgtSum() 46 0.02835 0.02865 0.02929 0.08% Gravity::get_new_grav_vector() 18 0.02723 0.02773 0.0285 0.08% MLMG::ResNormInf() 146 0.0107 0.01856 0.02681 0.08% Castro::reset_internal_energy(MultiFab) 80 0.02439 0.02489 0.02538 0.07% Castro::wd_update() 2 0.01243 0.018 0.02404 0.07% MLCellLinOp::defineAuxData() 19 0.01697 0.01804 0.02188 0.06% Amr::writePlotFile() 2 0.00688 0.01829 0.02164 0.06% FabArray::FillBoundary() 6380 0.01698 0.01768 0.02074 0.06% MLCGSolver::bicgstab 8 0 0.001053 0.01685 0.05% Castro::evaluate_source_change() 196 0.01472 0.01532 0.0164 0.05% FillPatchIterator::Initialize 187 0.01558 0.01584 0.01631 0.05% MLMG::prepareForSolve() 18 0.01281 0.01323 0.01419 0.04% MLCellLinOp::smooth() 1544 0.01058 0.01127 0.01277 0.04% MLMG::solve() 18 0.001502 0.004394 0.01207 0.03% Gravity::actual_solve_with_mlmg() 19 0.007467 0.007888 0.009486 0.03% MLCellLinOp::apply() 986 0.002444 0.002729 0.005789 0.02% Castro::initialize_do_advance() 20 0.004703 0.004922 0.005375 0.02% Castro::sum_integrated_quantities() 2 0.003313 0.004359 0.00516 0.01% FabArray::ParallelCopy() 1469 0.004187 0.004364 0.004636 0.01% AmrLevel::FillPatch() 187 0.003263 0.003389 0.00373 0.01% MLCellLinOp::correctionResidual() 768 0.002301 0.002379 0.00249 0.01% MLMG::mgVcycle() 128 0.00216 0.002247 0.002404 0.01% MLMG:computeResOfCorrection() 768 0.001839 0.001906 0.002012 0.01% AmrMesh::MakeNewGrids() 2 0.001753 0.001787 0.001828 0.01% Castro::save_data_for_retry() 3 0.001317 0.001681 0.001769 0.00% StateDataPhysBCFunct::() 187 0.001019 0.001211 0.001454 0.00% MLMG::actualBottomSolve() 128 0.000186 0.0002614 0.001301 0.00% Castro::subcycle_advance_ctu() 1 0.0002515 0.0003456 0.001174 0.00% FillPatchIterator::FillFromLevel0() 187 0.000895 0.0009331 0.0009789 0.00% main() 1 0.0007114 0.0007606 0.0008186 0.00% Gravity::solve_for_phi() 17 0.0006855 0.0007077 0.0007842 0.00% MLMG::mgVcycle_down::0 128 0.0006413 0.0006667 0.0007045 0.00% FillPatchSingleLevel 187 0.0006735 0.0006854 0.0007042 0.00% MLMG::mgVcycle_down::1 128 0.0006122 0.0006344 0.0006733 0.00% MLMG::mgVcycle_down::2 128 0.0005806 0.0006254 0.0006634 0.00% MLMG::mgVcycle_down::3 128 0.0005665 0.0005991 0.0006548 0.00% MLMG::mgVcycle_down::5 128 0.0005926 0.0006097 0.0006472 0.00% Castro::Castro() 1 0.0003826 0.000459 0.0006472 0.00% MLMG::mgVcycle_down::4 128 0.0005595 0.0005786 0.0006274 0.00% Castro::clean_state() 79 0.0004893 0.0005212 0.0005892 0.00% MLCellLinOp::solutionResidual() 147 0.0005187 0.0005414 0.0005707 0.00% Castro::do_old_sources() 40 0.0004748 0.000492 0.0005478 0.00% MLMG::oneIter() 128 0.0004357 0.0004538 0.0004833 0.00% MLMG::mgVcycle_up::0 128 0.0004102 0.0004248 0.0004604 0.00% MLMG::mgVcycle_up::1 128 0.0004079 0.0004238 0.0004563 0.00% MLMG::mgVcycle_up::3 128 0.0003739 0.0003928 0.0004407 0.00% MLMG::mgVcycle_up::2 128 0.0003859 0.000409 0.0004337 0.00% MLMG::mgVcycle_up::4 128 0.0003772 0.0003903 0.000422 0.00% MLMG::mgVcycle_up::5 128 0.0003687 0.0003825 0.0004004 0.00% MLMG::computeResidual() 128 0.0002989 0.0003107 0.0003261 0.00% Castro::construct_new_gravity() 17 0.0001789 0.0001969 0.0003135 0.00% TagBoxArray::mapPRD 2 0.000261 0.0002907 0.0003091 0.00% MLMG::mgVcycle_bottom 128 0.0002728 0.000284 0.0002993 0.00% Castro::initialize_advance() 1 0.0001108 0.0001313 0.0002802 0.00% Castro::problem_post_timestep() 2 1.257e-05 3.005e-05 0.000272 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002194 0.0002422 0.0002672 0.00% Castro::do_new_sources() 17 0.0002234 0.0002327 0.000248 0.00% Amr::coarseTimeStep() 1 0.0001329 0.0001556 0.0001601 0.00% MLPoisson::define() 19 0.0001265 0.0001364 0.0001486 0.00% Castro::check_for_nan() 37 0.0001251 0.0001308 0.0001442 0.00% Amr::defBaseLevel() 1 9.582e-05 0.0001054 0.0001431 0.00% Castro::retry_advance_ctu() 11 7.76e-05 8.875e-05 0.0001412 0.00% Castro::apply_source_to_state() 37 0.0001245 0.00013 0.0001407 0.00% Castro::construct_old_gravity() 20 0.0001236 0.0001296 0.0001402 0.00% Gravity::test_composite_phi() 1 0.0001055 0.0001157 0.0001267 0.00% Amr::grid_places() 2 5.254e-05 5.964e-05 8.118e-05 0.00% Castro::create_source_corrector() 17 5.593e-05 6.043e-05 7.243e-05 0.00% Castro::expand_state() 20 5.542e-05 6.044e-05 6.587e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.731e-05 5.857e-05 6.348e-05 0.00% MLMG::computeMLResidual() 18 5.216e-05 5.482e-05 5.905e-05 0.00% MLMG::getGradSolution() 18 5.183e-05 5.428e-05 5.817e-05 0.00% Amr::FinalizeInit() 1 8.856e-06 1.572e-05 5.587e-05 0.00% Castro::derive() 18 4.767e-05 4.994e-05 5.404e-05 0.00% Amr::regrid() 1 1.702e-05 2.108e-05 5.073e-05 0.00% Gravity::actual_multilevel_solve() 1 3.041e-05 3.157e-05 3.563e-05 0.00% Amr::InitializeInit() 1 8.23e-06 1.147e-05 3.359e-05 0.00% Castro::post_regrid() 2 1.122e-05 1.374e-05 2.749e-05 0.00% Amr::timeStep() 1 1.511e-05 1.677e-05 2.233e-05 0.00% Castro::post_init() 1 1.638e-05 1.818e-05 1.963e-05 0.00% Castro::advance() 1 8.141e-06 9.625e-06 1.369e-05 0.00% Amr::bldFineLevels() 1 7.477e-06 8.403e-06 9.909e-06 0.00% Amr::init() 1 6.643e-06 7.79e-06 9.196e-06 0.00% Castro::post_timestep() 1 6.257e-06 6.697e-06 7.648e-06 0.00% Amr::initialInit() 1 3.124e-06 3.447e-06 4.086e-06 0.00% Other 10929 0.1232 0.1482 0.2099 0.59% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 35.47 35.47 35.48 100.00% Amr::coarseTimeStep() 1 32.49 32.49 32.49 91.58% Amr::timeStep() 1 32.37 32.37 32.37 91.24% Castro::advance() 1 32.21 32.21 32.21 90.80% Castro::subcycle_advance_ctu() 1 32.08 32.09 32.09 90.45% Castro::do_advance_ctu() 20 31.93 31.94 31.95 90.06% Castro::construct_ctu_hydro_source() 20 16.3 16.34 16.39 46.20% Castro::clean_state() 79 8.362 8.497 8.595 24.23% Castro::computeTemp() 80 5.18 5.285 5.368 15.13% Gravity::solve_phi_with_mlmg() 19 4.933 4.965 4.999 14.09% Castro::construct_new_gravity() 17 4.707 4.738 4.77 13.45% Gravity::solve_for_phi() 17 4.575 4.607 4.64 13.08% Gravity::actual_solve_with_mlmg() 19 4.26 4.272 4.285 12.08% Castro::reset_internal_energy(MultiFab) 80 4.181 4.223 4.271 12.04% MLMG::solve() 18 4.178 4.184 4.185 11.80% Castro::evaluate_source_change() 196 3.707 3.744 3.789 10.68% MLMG::oneIter() 128 3.659 3.665 3.67 10.35% MLMG::mgVcycle() 128 3.606 3.611 3.616 10.19% Castro::do_old_sources() 40 3.378 3.421 3.463 9.76% Castro::initialize_do_advance() 20 3.318 3.368 3.405 9.60% amrex::Copy() 766 3.253 3.289 3.329 9.38% Castro::enforce_min_density() 79 3.102 3.134 3.178 8.96% FabArray::ParallelCopy() 1469 1.439 2.203 2.61 7.36% MLCellLinOp::smooth() 1544 1.975 2.13 2.577 7.26% Castro::do_new_sources() 17 2.445 2.491 2.523 7.11% AmrLevel::FillPatch() 187 2.165 2.225 2.335 6.58% MultiFab::AddProduct() 3598 2.201 2.222 2.251 6.34% MLPoisson::Fsmooth() 6144 1.469 1.593 1.902 5.36% FabArray::ParallelCopy_finish() 1469 0.5578 1.335 1.75 4.93% FillPatchIterator::Initialize 187 1.355 1.422 1.54 4.34% Amr::init() 1 1.518 1.518 1.518 4.28% FillPatchIterator::FillFromLevel0() 187 1.307 1.381 1.507 4.25% FillPatchSingleLevel 187 1.306 1.38 1.506 4.25% Amr::InitAmr() 1 1.462 1.464 1.469 4.14% FabArray::setVal() 5232 1.36 1.382 1.409 3.97% Amr::initialInit() 1 1.4 1.4 1.4 3.95% MultiFab::Subtract() 159 1.37 1.383 1.399 3.94% MLMG::mgVcycle_up::1 128 0.2751 0.7454 1.218 3.43% MLMG::mgVcycle_down::0 128 1.105 1.111 1.118 3.15% MLMG::addInterpCorrection() 768 0.1455 0.834 1.045 2.95% FabArray::ParallelCopy_nowait() 1469 0.7963 0.8629 0.9359 2.64% Castro::expand_state() 20 0.7294 0.7897 0.8475 2.39% MLCellLinOp::applyBC() 7148 0.5805 0.6556 0.8259 2.33% MLCellLinOp::apply() 986 0.6804 0.7152 0.8191 2.31% Amr::FinalizeInit() 1 0.6433 0.774 0.8075 2.28% MLMG::mgVcycle_up::0 128 0.771 0.7759 0.7798 2.20% MLMG::mgVcycle_up::2 128 0.01535 0.3136 0.766 2.16% FabArray::sum() 3598 0.7475 0.7548 0.7637 2.15% Amr::InitializeInit() 1 0.5922 0.6257 0.7564 2.13% Amr::defBaseLevel() 1 0.5922 0.6257 0.7564 2.13% Gravity::fill_multipole_BCs() 19 0.6601 0.688 0.7187 2.03% check_for_negative_density() 20 0.2975 0.4889 0.6706 1.89% MLPoisson::Fapply() 986 0.5587 0.5865 0.6541 1.84% MLMG:computeResOfCorrection() 768 0.4644 0.5016 0.5945 1.68% Castro::construct_old_source() 240 0.5658 0.5731 0.593 1.67% MLCellLinOp::correctionResidual() 768 0.4626 0.4997 0.5925 1.67% FabArray::FillBoundary() 6380 0.4195 0.4818 0.539 1.52% Castro::construct_new_source() 102 0.4693 0.4743 0.4793 1.35% Castro::Castro() 1 0.3986 0.4063 0.4477 1.26% MLMG::mgVcycle_up::3 128 0.01527 0.05009 0.4468 1.26% Castro::post_init() 1 0.428 0.4311 0.4326 1.22% Castro::initMFs() 1 0.3806 0.388 0.4291 1.21% MLCellLinOp::solutionResidual() 147 0.4045 0.415 0.4211 1.19% MLMG::mgVcycle_down::1 128 0.3693 0.3792 0.3921 1.11% Castro::reset_internal_energy(Fab) 160 0.3682 0.372 0.3759 1.06% MLMG::computeResidual() 128 0.3503 0.3542 0.3568 1.01% Castro::initData() 1 0.1917 0.2175 0.3512 0.99% Castro::construct_old_gravity_source() 40 0.3228 0.3283 0.348 0.98% Gravity::multilevel_solve_for_new_phi() 1 0.3265 0.3265 0.3266 0.92% Gravity::actual_multilevel_solve() 1 0.3264 0.3265 0.3266 0.92% FabArray::Saxpy() 111 0.3174 0.3203 0.324 0.91% FillBoundary_finish() 6380 0.2243 0.263 0.3204 0.90% Castro::create_source_corrector() 17 0.1768 0.2177 0.2769 0.78% Amr::bldFineLevels() 1 0.1084 0.2396 0.2726 0.77% Castro::construct_new_gravity_source() 17 0.2646 0.2679 0.2713 0.76% FillBoundary_nowait() 6380 0.1543 0.2011 0.258 0.73% Castro::construct_old_gravity() 20 0.1742 0.2043 0.2395 0.68% Gravity::get_old_grav_vector() 20 0.174 0.2041 0.2394 0.67% Amr::writePlotFile() 2 0.2238 0.2239 0.2239 0.63% FabArray::Xpay() 915 0.1903 0.1994 0.2185 0.62% Castro::post_regrid() 2 0.2051 0.2066 0.2097 0.59% MLMG::mgVcycle_down::2 128 0.01797 0.108 0.2076 0.59% MLMG::mgVcycle_down::3 128 0.01732 0.03904 0.1899 0.54% Castro::derive() 18 0.166 0.1682 0.171 0.48% AmrLevel::derive() 18 0.166 0.1682 0.1709 0.48% Castro::check_for_nan() 37 0.1649 0.1673 0.1701 0.48% MultiFab::contains_nan() 37 0.1648 0.1672 0.17 0.48% Amr::grid_places() 2 0.005615 0.1363 0.1698 0.48% AmrMesh::MakeNewGrids() 2 0.00556 0.1362 0.1697 0.48% TagBoxArray::mapPRD 2 0.001702 0.1096 0.1657 0.47% Castro::post_timestep() 1 0.1551 0.1593 0.1631 0.46% Castro::retry_advance_ctu() 11 0.1353 0.1375 0.156 0.44% Castro::normalize_species() 79 0.1388 0.1406 0.1455 0.41% Gravity::get_new_grav_vector() 18 0.1272 0.1302 0.1374 0.39% Castro::initialize_advance() 1 0.1226 0.1238 0.1273 0.36% amrex::average_down 768 0.07147 0.08649 0.1234 0.35% TagBoxArray::collate() 2 0.000546 0.02327 0.113 0.32% Castro::apply_source_to_state() 37 0.1059 0.107 0.1084 0.31% Amr::regrid() 1 0.1055 0.106 0.1075 0.30% MLMG::mgVcycle_bottom 128 0.0004602 0.007062 0.1059 0.30% MLMG::actualBottomSolve() 128 0.000186 0.006778 0.1056 0.30% MLMG::mgVcycle_down::4 128 0.01726 0.02312 0.1045 0.29% MLCGSolver::bicgstab 8 0 0.006454 0.1033 0.29% StateDataPhysBCFunct::() 187 0.001019 0.05235 0.1023 0.29% StateData::FillBoundary(geom) 193 0 0.05114 0.1009 0.28% Castro::save_data_for_retry() 3 0.07066 0.07228 0.08835 0.25% MLMG::mgVcycle_up::4 128 0.01535 0.01947 0.07624 0.21% Castro::sum_integrated_quantities() 2 0.07416 0.07434 0.07436 0.21% MLMG::computeMLResidual() 18 0.04869 0.05827 0.06448 0.18% MLPoisson::define() 19 0.03741 0.04885 0.0639 0.18% MLMG::ResNormInf() 146 0.04387 0.05073 0.05852 0.16% amrex::Add() 128 0.0524 0.05314 0.05399 0.15% MLMG::prepareForSolve() 18 0.04459 0.04701 0.05334 0.15% VisMF::Write(FabArray) 4 0.03735 0.03852 0.05269 0.15% MLMG::mgVcycle_down::5 128 0.01731 0.01973 0.04992 0.14% FabArray::setDomainBndry() 187 0.0006601 0.02525 0.04927 0.14% Gravity::test_composite_phi() 1 0.04138 0.04145 0.04149 0.12% FabArray::norminf() 244 0.03587 0.03681 0.04106 0.12% MLMG::mgVcycle_up::5 128 0.01539 0.0173 0.03967 0.11% Castro::problem_post_timestep() 2 0.02514 0.03114 0.03651 0.10% Castro::wd_update() 2 0.02512 0.03111 0.0365 0.10% MLCellLinOp::defineAuxData() 19 0.02212 0.0239 0.02955 0.08% Castro::volWgtSum() 46 0.02835 0.02865 0.02929 0.08% MLMG::getGradSolution() 18 0.02672 0.02724 0.02774 0.08% Other 10929 0.1825 0.2361 0.3365 0.95% -------------------------------------------------------------------------------------------- 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 428 KiB 433 KiB 438 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 1953 B 1999 B 2050 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 0 B 0 B 1 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 1732 KiB 1732 KiB 1732 KiB 2329 KiB 2329 KiB 2329 KiB Castro::construct_ctu_hydro_source() 6400 428 KiB 433 KiB 438 KiB 1176 KiB 1176 KiB 1176 KiB StateData::define() 192 1116 KiB 1116 KiB 1116 KiB 1164 KiB 1164 KiB 1164 KiB Castro::initMFs() 256 1094 KiB 1094 KiB 1094 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 735 KiB 735 KiB 735 KiB 803 KiB 803 KiB 803 KiB Amr::writePlotFile() 64 4416 B 4490 B 4754 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 69 KiB 69 KiB 70 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 18 KiB 19 KiB 20 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 474 KiB 474 KiB 474 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 3565 B 5534 B 7264 B 225 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1233 B 1249 B 1269 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 29 KiB 29 KiB 30 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 143 KiB 143 KiB 143 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 136 KiB 136 KiB 136 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 12 KiB 12 KiB 14 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 1280 778 B 784 B 794 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 658 B 663 B 672 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 1953 B 1999 B 2050 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 20 B 20 B 27 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 35 KiB 35 KiB 35 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 61 B 61 B 61 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 4 B 5 B 6 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 164 B 168 B 177 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 225 B 264 B 311 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 34 KiB 34 KiB 34 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 307 B 311 B 315 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 1709 B 1726 B 1747 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 150 B 150 B 150 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 2111 B 2126 B 2142 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 1 B 1 B 1 B 16 KiB 16 KiB 16 KiB MLCellLinOp::defineAuxData() 11704 551 B 616 B 855 B 4608 B 5151 B 7152 B BndryData::define() 4864 578 B 578 B 579 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 18 B 28 B 37 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 0 B 7 B 12 B 2916 B 2916 B 2916 B OwnerMask() 64 0 B 0 B 0 B 2592 B 2592 B 2592 B MLMG::addInterpCorrection() 5632 1 B 17 B 28 B 1024 B 1024 B 1024 B amrex::average_down 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 0 B 2 B 3 B 648 B 648 B 648 B MLCGSolver::bicgstab 640 0 B 0 B 1 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 0 B 0 B 0 B 8 B 8 B 8 B ----------------------------------------------------------------------------------------------------------------------- AMReX (26.01-17-g909028cb95b3) finalized