Initializing AMReX (26.03-47-g384b99a08d34)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.03-47-g384b99a08d34) initialized Starting run at 06:55:58 UTC on 2026-03-09. Successfully read inputs file ... Castro git describe: 26.03-13-g58c083c59 AMReX git describe: 26.03-47-g384b99a08 Microphysics git describe: 26.03-4-g0b9339e8 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.197331288 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001218526 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010227465 Gravity::multilevel_solve_for_new_phi() time = 0.107322829 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010253077 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.048186222 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004222441 seconds [Level 0 step 1] ADVANCE at time 0 with dt = 0.25  Beginning subcycle 1 starting at time 0 with dt = 0.25 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.7670713549 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.7670713549 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071148 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.769409 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477044 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035502 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -0.01295059363 at index 21, 83, 0 Invalid density = -0.01295059363 at index 21, 172, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035502 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477044 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.769409 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071148 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Castro::construct_old_gravity() time = 0.003523524 on level 0 Castro::do_old_sources() time = 0.007498974 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.054593125 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.005162039 on level 0 Castro::do_old_sources() time = 0.007687863 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.053735695 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.003186708 on level 0 Castro::do_old_sources() time = 0.007518257 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.050804353 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.005172543 on level 0 Castro::do_old_sources() time = 0.008061417 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.050134657 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.298184099 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.014258796 Gravity::solve_for_phi() time = 0.090594098 on level 0 Castro::construct_new_gravity() time = 0.0926996 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.050098146 on level 0 Castro::construct_old_gravity() time = 0.010725947 on level 0 Castro::do_old_sources() time = 0.013178459 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.049492257 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.34728298 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.0134009 Gravity::solve_for_phi() time = 0.098448854 on level 0 Castro::construct_new_gravity() time = 0.100632838 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.057767828 on level 0 Castro::construct_old_gravity() time = 0.025683683 on level 0 Castro::do_old_sources() time = 0.008659252 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.051296757 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.30367784 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.013859821 Gravity::solve_for_phi() time = 0.090344214 on level 0 Castro::construct_new_gravity() time = 0.092399402 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.049465484 on level 0 Castro::construct_old_gravity() time = 0.002665045 on level 0 Castro::do_old_sources() time = 0.007487731 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.04875866 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.306271951 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.014635768 Gravity::solve_for_phi() time = 0.090719621 on level 0 Castro::construct_new_gravity() time = 0.092790505 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.049576004 on level 0 Castro::construct_old_gravity() time = 0.002743964 on level 0 Castro::do_old_sources() time = 0.007534728 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.048869175 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.31319371 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.014035456 Gravity::solve_for_phi() time = 0.090169598 on level 0 Castro::construct_new_gravity() time = 0.092269162 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.049815271 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.002742244 on level 0 Castro::do_old_sources() time = 0.007505435 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.049017286 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.306426666 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.01467613 Gravity::solve_for_phi() time = 0.090694569 on level 0 Castro::construct_new_gravity() time = 0.092726412 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.049587881 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.002660658 on level 0 Castro::do_old_sources() time = 0.007485253 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.048886664 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.311406243 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.014397898 Gravity::solve_for_phi() time = 0.091387106 on level 0 Castro::construct_new_gravity() time = 0.093496441 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.050222992 on level 0 Castro::construct_old_gravity() time = 0.002871659 on level 0 Castro::do_old_sources() time = 0.007508118 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.04904976 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.306806674 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.014284166 Gravity::solve_for_phi() time = 0.090701408 on level 0 Castro::construct_new_gravity() time = 0.092746009 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.049536569 on level 0 Castro::construct_old_gravity() time = 0.002740409 on level 0 Castro::do_old_sources() time = 0.00751106 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.048836351 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.311674671 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.014710651 Gravity::solve_for_phi() time = 0.09068138 on level 0 Castro::construct_new_gravity() time = 0.092723205 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.049618054 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.002761836 on level 0 Castro::do_old_sources() time = 0.007479172 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.0488412 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.306686721 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.014324166 Gravity::solve_for_phi() time = 0.090162132 on level 0 Castro::construct_new_gravity() time = 0.092207361 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.049500878 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.002691372 on level 0 Castro::do_old_sources() time = 0.007552477 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.048999874 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.312656464 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.015528962 Gravity::solve_for_phi() time = 0.091508593 on level 0 Castro::construct_new_gravity() time = 0.093610982 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.049509767 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.003042689 on level 0 Castro::do_old_sources() time = 0.007488187 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.048997205 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.306663799 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.014742141 Gravity::solve_for_phi() time = 0.090597058 on level 0 Castro::construct_new_gravity() time = 0.092633075 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.049381921 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.00269222 on level 0 Castro::do_old_sources() time = 0.007435151 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.04852685 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.311930017 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.014675814 Gravity::solve_for_phi() time = 0.090717435 on level 0 Castro::construct_new_gravity() time = 0.092759211 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.049591214 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.00275555 on level 0 Castro::do_old_sources() time = 0.007502371 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.048651547 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.307490622 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.014722116 Gravity::solve_for_phi() time = 0.091397296 on level 0 Castro::construct_new_gravity() time = 0.093512605 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.049733801 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.003279373 on level 0 Castro::do_old_sources() time = 0.007983713 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.048658659 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.315209623 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.015022679 Gravity::solve_for_phi() time = 0.091984522 on level 0 Castro::construct_new_gravity() time = 0.094073236 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.049841844 on level 0 Castro::construct_old_gravity() time = 0.002726201 on level 0 Castro::do_old_sources() time = 0.007512272 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.049382409 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.308576813 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.014619362 Gravity::solve_for_phi() time = 0.090597108 on level 0 Castro::construct_new_gravity() time = 0.092716625 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.049697859 on level 0 Castro::construct_old_gravity() time = 0.002751093 on level 0 Castro::do_old_sources() time = 0.007489489 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.048740789 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.312526289 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002889497469 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.014579898 Gravity::solve_for_phi() time = 0.090590833 on level 0 Castro::construct_new_gravity() time = 0.092631881 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.049614681 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.39124193 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.047860365 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004262638 seconds Ending run at 06:56:11 UTC on 2026-03-09. Run time = 13.66093668 Run time without initialization = 11.44431417 Average number of zones advanced per microsecond: 0.003 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 13.6534 ... 13.6557 ... 13.6610 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.1173 5.2090 5.2903 38.73% Amr::InitAmr() 1 1.284 1.286 1.292 9.45% Castro::computeTemp() 80 0.7157 0.8051 0.923 6.76% amrex::Copy() 766 0.8087 0.8217 0.8393 6.14% FabArray::ParallelCopy_finish() 1469 0.2377 0.5883 0.8247 6.04% MultiFab::AddProduct() 3598 0.5199 0.5275 0.5365 3.93% MLPoisson::Fsmooth() 6144 0.334 0.3642 0.448 3.28% Castro::initMFs() 1 0.3882 0.3929 0.4111 3.01% FabArray::setVal() 5232 0.3758 0.3832 0.3964 2.90% Castro::reset_internal_energy(Fab) 160 0.3431 0.3492 0.3558 2.60% MultiFab::Subtract() 159 0.3322 0.3369 0.3432 2.51% check_for_negative_density() 20 0.08156 0.1814 0.2796 2.05% Castro::do_advance_ctu() 20 0.261 0.2662 0.2731 2.00% FabArray::ParallelCopy_nowait() 1469 0.2255 0.2458 0.2715 1.99% Gravity::fill_multipole_BCs() 19 0.1899 0.2411 0.2643 1.93% TagBoxArray::collate() 2 0.0002813 0.03973 0.1881 1.38% FabArray::sum() 3598 0.1798 0.1826 0.1855 1.36% Castro::initData() 1 0.003161 0.02918 0.1726 1.26% FillBoundary_finish() 6380 0.1104 0.1324 0.1556 1.14% FillBoundary_nowait() 6380 0.08797 0.1183 0.1539 1.13% MLCellLinOp::applyBC() 7148 0.07045 0.08518 0.1345 0.98% Castro::normalize_species() 79 0.1145 0.1163 0.1184 0.87% MLPoisson::Fapply() 986 0.08702 0.09193 0.1091 0.80% Castro::construct_old_gravity_source() 40 0.09413 0.09564 0.09766 0.71% FabArray::Saxpy() 111 0.07937 0.08061 0.08204 0.60% Castro::construct_new_gravity_source() 17 0.0668 0.06847 0.07062 0.52% FabArray::Xpay() 915 0.04863 0.05178 0.06274 0.46% MultiFab::contains_nan() 37 0.04616 0.04712 0.05048 0.37% MLMG::addInterpCorrection() 768 0.03195 0.03473 0.04292 0.31% AmrLevel::derive() 18 0.03881 0.03944 0.04004 0.29% StateData::FillBoundary(geom) 193 0 0.01878 0.03674 0.27% amrex::average_down 768 0.0245 0.02695 0.03436 0.25% Castro::reset_internal_energy(MultiFab) 80 0.02416 0.02477 0.02544 0.19% VisMF::Write(FabArray) 4 0.01337 0.01453 0.0243 0.18% MLCellLinOp::defineAuxData() 19 0.01579 0.0169 0.02162 0.16% FabArray::FillBoundary() 6380 0.01658 0.01734 0.02091 0.15% Castro::construct_old_source() 240 0.01889 0.01924 0.01957 0.14% Castro::enforce_min_density() 79 0.01846 0.01878 0.01908 0.14% Amr::writePlotFile() 2 0.006981 0.01587 0.01762 0.13% FabArrayBase::getFB() 6380 0.01341 0.01404 0.01731 0.13% MLCGSolver::bicgstab 8 0 0.001073 0.01718 0.13% Castro::evaluate_source_change() 196 0.01433 0.01487 0.01636 0.12% FabArray::setDomainBndry() 187 0.0006844 0.008568 0.01619 0.12% FillPatchIterator::Initialize 187 0.01517 0.01554 0.01605 0.12% Castro::construct_new_source() 102 0.01476 0.01502 0.01542 0.11% FabArray::norminf() 244 0.01037 0.01083 0.0146 0.11% amrex::Add() 128 0.01341 0.01361 0.01386 0.10% MLCellLinOp::smooth() 1544 0.01072 0.01127 0.01299 0.10% Castro::wd_update() 2 0.004176 0.008733 0.01231 0.09% MLMG::prepareForSolve() 18 0.01001 0.01045 0.01176 0.09% MLMG::ResNormInf() 146 0.003003 0.007978 0.01057 0.08% MLMG::solve() 18 0.0008037 0.005235 0.01036 0.08% Gravity::actual_solve_with_mlmg() 19 0.007431 0.007939 0.009877 0.07% FabArrayBase::getCPC() 2045 0.008632 0.008932 0.009482 0.07% Gravity::get_old_grav_vector() 20 0.006091 0.006184 0.006298 0.05% Gravity::get_new_grav_vector() 18 0.005604 0.005735 0.005999 0.04% MLCGSolver::ParallelAllReduce 189 0 0.0003598 0.005756 0.04% MLCellLinOp::apply() 986 0.002388 0.002681 0.00573 0.04% Castro::initialize_do_advance() 20 0.004622 0.004855 0.005148 0.04% FabArray::ParallelCopy() 1469 0.00424 0.00443 0.004692 0.03% AmrLevel::FillPatch() 187 0.003083 0.003194 0.003604 0.03% Castro::sum_integrated_quantities() 2 0.001578 0.002159 0.002538 0.02% MLCellLinOp::correctionResidual() 768 0.002273 0.002353 0.002436 0.02% MLMG::mgVcycle() 128 0.002141 0.002221 0.002397 0.02% MLLinOp::defineGrids() 19 0.001896 0.002065 0.002347 0.02% MLMG:computeResOfCorrection() 768 0.001833 0.001896 0.002041 0.01% StateDataPhysBCFunct::() 187 0.0009869 0.001176 0.001413 0.01% Castro::save_data_for_retry() 3 0.001291 0.001321 0.00138 0.01% MLMG::actualBottomSolve() 128 0.0001831 0.0002605 0.001316 0.01% FillPatchIterator::FillFromLevel0() 187 0.0008715 0.000901 0.0009708 0.01% Castro::subcycle_advance_ctu() 1 0.0002452 0.0003169 0.0009566 0.01% main() 1 0.0007309 0.0008008 0.0008785 0.01% Gravity::solve_for_phi() 17 0.0006563 0.0006952 0.000767 0.01% MLMG::mgVcycle_down::0 128 0.0006293 0.0006487 0.0006928 0.01% FillPatchSingleLevel 187 0.0006322 0.000656 0.0006823 0.00% MLMG::mgVcycle_down::1 128 0.0006053 0.0006322 0.0006746 0.00% MLMG::mgVcycle_down::2 128 0.0005896 0.0006317 0.0006642 0.00% MLMG::mgVcycle_down::5 128 0.0005885 0.0006146 0.0006612 0.00% MLMG::mgVcycle_down::3 128 0.0005695 0.0006037 0.0006585 0.00% MLMG::mgVcycle_down::4 128 0.0005601 0.0005843 0.0006491 0.00% AmrMesh::MakeNewGrids() 2 0.0005314 0.0005479 0.000583 0.00% Castro::Castro() 1 0.0004553 0.0005095 0.000576 0.00% MLCellLinOp::solutionResidual() 147 0.0005156 0.000532 0.0005495 0.00% Castro::clean_state() 79 0.0004662 0.000495 0.0005238 0.00% Castro::do_old_sources() 40 0.0004515 0.0004771 0.0005089 0.00% MLMG::oneIter() 128 0.000431 0.0004479 0.0005043 0.00% MLMG::mgVcycle_up::1 128 0.0004197 0.0004366 0.0004776 0.00% MLMG::mgVcycle_up::0 128 0.0004186 0.0004342 0.0004693 0.00% MLMG::mgVcycle_up::2 128 0.0003978 0.0004224 0.0004565 0.00% MLMG::mgVcycle_up::4 128 0.0003926 0.0004058 0.0004439 0.00% MLMG::mgVcycle_up::3 128 0.0003858 0.0003998 0.0004372 0.00% MLMG::mgVcycle_up::5 128 0.000374 0.0003866 0.0004018 0.00% MLMG::computeResidual() 128 0.0002919 0.0003049 0.0003145 0.00% MLMG::mgVcycle_bottom 128 0.0002751 0.0002849 0.0003079 0.00% Castro::construct_new_gravity() 17 0.0001677 0.0001918 0.0003059 0.00% Castro::problem_post_timestep() 2 1.203e-05 2.995e-05 0.0002813 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002155 0.0002325 0.0002558 0.00% Castro::do_new_sources() 17 0.0002119 0.0002275 0.0002488 0.00% MLPoisson::define() 19 0.00014 0.0001654 0.0001877 0.00% Amr::coarseTimeStep() 1 0.0001424 0.0001681 0.0001727 0.00% Castro::check_for_nan() 37 0.0001199 0.0001282 0.0001483 0.00% Castro::initialize_advance() 1 0.0001083 0.0001233 0.0001455 0.00% Amr::defBaseLevel() 1 9.328e-05 0.0001055 0.0001391 0.00% Castro::construct_old_gravity() 20 0.0001178 0.0001263 0.0001368 0.00% Castro::apply_source_to_state() 37 0.0001178 0.0001234 0.0001347 0.00% Gravity::test_composite_phi() 1 0.0001084 0.000118 0.0001334 0.00% Castro::retry_advance_ctu() 11 6.809e-05 7.842e-05 0.0001301 0.00% TagBoxArray::mapPRD 2 0.0001164 0.0001203 0.0001259 0.00% MLLinOp::define() 19 8.374e-05 8.751e-05 9.242e-05 0.00% Amr::grid_places() 2 5.37e-05 5.967e-05 7.892e-05 0.00% Castro::derive() 18 4.645e-05 5.078e-05 6.634e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.93e-05 6.102e-05 6.633e-05 0.00% Castro::expand_state() 20 5.549e-05 6.005e-05 6.447e-05 0.00% Amr::FinalizeInit() 1 1.191e-05 2.069e-05 6.274e-05 0.00% Castro::create_source_corrector() 17 5.385e-05 5.741e-05 6.081e-05 0.00% MLMG::computeMLResidual() 18 5.101e-05 5.442e-05 5.71e-05 0.00% Amr::regrid() 1 2.007e-05 2.85e-05 5.101e-05 0.00% Gravity::actual_multilevel_solve() 1 3.154e-05 3.519e-05 4.978e-05 0.00% Amr::timeStep() 1 1.52e-05 1.78e-05 2.323e-05 0.00% Amr::InitializeInit() 1 8.901e-06 1.106e-05 2.169e-05 0.00% Castro::post_regrid() 2 1.15e-05 1.303e-05 1.937e-05 0.00% Castro::post_init() 1 1.229e-05 1.365e-05 1.874e-05 0.00% Castro::advance() 1 8.102e-06 1.019e-05 1.632e-05 0.00% Amr::bldFineLevels() 1 7.459e-06 8.678e-06 1.039e-05 0.00% Amr::init() 1 6.749e-06 7.707e-06 9.533e-06 0.00% Castro::post_timestep() 1 6.088e-06 6.441e-06 7.575e-06 0.00% Amr::initialInit() 1 3.067e-06 3.655e-06 5.173e-06 0.00% Other 2341 0.04808 0.05956 0.09236 0.68% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 13.65 13.66 13.66 100.00% Amr::coarseTimeStep() 1 11.44 11.44 11.44 83.77% Amr::timeStep() 1 11.39 11.39 11.39 83.39% Castro::advance() 1 11.34 11.34 11.34 83.01% Castro::subcycle_advance_ctu() 1 11.29 11.3 11.3 82.70% Castro::do_advance_ctu() 20 11.26 11.26 11.26 82.43% Castro::construct_ctu_hydro_source() 20 6.004 6.023 6.057 44.34% Castro::clean_state() 79 2.887 3.009 3.133 22.94% Castro::computeTemp() 80 2.02 2.129 2.241 16.40% Gravity::solve_phi_with_mlmg() 19 1.593 1.645 1.669 12.22% Castro::construct_new_gravity() 17 1.511 1.563 1.586 11.61% Gravity::solve_for_phi() 17 1.474 1.527 1.551 11.35% Gravity::actual_solve_with_mlmg() 19 1.397 1.401 1.404 10.28% Castro::reset_internal_energy(MultiFab) 80 1.303 1.324 1.348 9.87% MLMG::solve() 18 1.339 1.346 1.348 9.87% Amr::InitAmr() 1 1.284 1.286 1.292 9.45% MLMG::oneIter() 128 1.207 1.21 1.213 8.88% Castro::initialize_do_advance() 20 1.142 1.171 1.201 8.79% MLMG::mgVcycle() 128 1.193 1.196 1.199 8.77% Castro::do_old_sources() 40 1.106 1.126 1.142 8.36% FabArray::ParallelCopy() 1469 0.5131 0.8595 1.092 8.00% Castro::evaluate_source_change() 196 0.9207 0.9347 0.9505 6.96% Amr::init() 1 0.9249 0.9249 0.9249 6.77% Amr::initialInit() 1 0.8719 0.8719 0.872 6.38% Castro::do_new_sources() 17 0.7643 0.8096 0.8519 6.24% amrex::Copy() 766 0.8087 0.8217 0.8393 6.14% FabArray::ParallelCopy_finish() 1469 0.2377 0.5883 0.8247 6.04% MLCellLinOp::smooth() 1544 0.6206 0.6727 0.8114 5.94% Castro::enforce_min_density() 79 0.7757 0.7874 0.8019 5.87% AmrLevel::FillPatch() 187 0.5913 0.676 0.7642 5.59% Amr::InitializeInit() 1 0.4616 0.4933 0.6573 4.81% Amr::defBaseLevel() 1 0.4615 0.4933 0.6573 4.81% FillPatchIterator::Initialize 187 0.3852 0.4725 0.5598 4.10% MultiFab::AddProduct() 3598 0.5199 0.5275 0.5365 3.93% FillPatchIterator::FillFromLevel0() 187 0.3534 0.4484 0.5335 3.91% FillPatchSingleLevel 187 0.3525 0.4475 0.5326 3.90% MLMG::mgVcycle_up::1 128 0.1047 0.2832 0.4638 3.39% MLPoisson::Fsmooth() 6144 0.334 0.3642 0.448 3.28% MLCellLinOp::applyBC() 7148 0.3326 0.368 0.443 3.24% Castro::Castro() 1 0.3941 0.3989 0.4186 3.06% Castro::initMFs() 1 0.3884 0.3931 0.4113 3.01% Amr::FinalizeInit() 1 0.2146 0.3786 0.4103 3.00% MLMG::addInterpCorrection() 768 0.06536 0.3274 0.4026 2.95% FabArray::setVal() 5232 0.3758 0.3832 0.3964 2.90% Castro::reset_internal_energy(Fab) 160 0.3431 0.3492 0.3558 2.60% MultiFab::Subtract() 159 0.3322 0.3369 0.3432 2.51% FabArray::FillBoundary() 6380 0.2503 0.2828 0.3149 2.31% MLMG::mgVcycle_up::2 128 0.01553 0.131 0.3085 2.26% FabArray::ParallelCopy_nowait() 1469 0.2459 0.2668 0.2936 2.15% check_for_negative_density() 20 0.08156 0.1814 0.2796 2.05% MLMG::mgVcycle_down::0 128 0.2633 0.2662 0.2715 1.99% Gravity::fill_multipole_BCs() 19 0.1919 0.2431 0.2664 1.95% Castro::expand_state() 20 0.1954 0.2258 0.258 1.89% Castro::initData() 1 0.06558 0.09261 0.2369 1.73% Amr::bldFineLevels() 1 0.03834 0.202 0.2334 1.71% MLCellLinOp::apply() 986 0.1518 0.1633 0.2047 1.50% MLMG::mgVcycle_up::0 128 0.2019 0.2027 0.2034 1.49% MLMG::mgVcycle_up::3 128 0.01547 0.03001 0.2032 1.49% Amr::grid_places() 2 0.002862 0.1668 0.1985 1.45% AmrMesh::MakeNewGrids() 2 0.002791 0.1668 0.1985 1.45% TagBoxArray::mapPRD 2 0.0008684 0.1257 0.1966 1.44% TagBoxArray::collate() 2 0.0002813 0.03973 0.1881 1.38% FabArray::sum() 3598 0.1798 0.1826 0.1855 1.36% FillBoundary_nowait() 6380 0.1019 0.1331 0.1692 1.24% Castro::construct_old_source() 240 0.1573 0.1597 0.1624 1.19% FillBoundary_finish() 6380 0.1104 0.1324 0.1556 1.14% MLMG:computeResOfCorrection() 768 0.1162 0.1262 0.15 1.10% MLCellLinOp::correctionResidual() 768 0.1144 0.1243 0.148 1.08% Castro::post_init() 1 0.1403 0.1418 0.1428 1.05% MLMG::mgVcycle_down::1 128 0.1266 0.1299 0.1344 0.98% Castro::construct_new_source() 102 0.119 0.1215 0.1238 0.91% Castro::normalize_species() 79 0.1145 0.1163 0.1184 0.87% MLPoisson::Fapply() 986 0.08702 0.09193 0.1091 0.80% Gravity::multilevel_solve_for_new_phi() 1 0.1072 0.1072 0.1073 0.79% Gravity::actual_multilevel_solve() 1 0.1071 0.1072 0.1073 0.79% MLCellLinOp::solutionResidual() 147 0.08506 0.09171 0.09882 0.72% Castro::create_source_corrector() 17 0.05538 0.07188 0.09864 0.72% Castro::construct_old_gravity_source() 40 0.09413 0.09564 0.09766 0.71% Amr::writePlotFile() 2 0.09646 0.09656 0.09657 0.71% Castro::construct_old_gravity() 20 0.0444 0.06012 0.08704 0.64% Gravity::get_old_grav_vector() 20 0.04427 0.06 0.0869 0.64% MLMG::mgVcycle_bottom 128 0.0004625 0.005701 0.08412 0.62% MLMG::actualBottomSolve() 128 0.0001831 0.005417 0.08381 0.61% FabArray::Saxpy() 111 0.07937 0.08061 0.08204 0.60% MLCGSolver::bicgstab 8 0 0.005108 0.08173 0.60% MLMG::mgVcycle_down::2 128 0.01815 0.0479 0.08004 0.59% MLMG::computeResidual() 128 0.07374 0.0744 0.07806 0.57% Castro::post_regrid() 2 0.06897 0.06997 0.07145 0.52% Castro::construct_new_gravity_source() 17 0.0668 0.06847 0.07062 0.52% Castro::derive() 18 0.06748 0.06878 0.0699 0.51% AmrLevel::derive() 18 0.06743 0.06873 0.06985 0.51% MLMG::mgVcycle_down::3 128 0.01767 0.02437 0.06862 0.50% FabArray::Xpay() 915 0.04863 0.05178 0.06274 0.46% amrex::average_down 768 0.03431 0.04156 0.05989 0.44% Castro::post_timestep() 1 0.05081 0.05312 0.0559 0.41% Castro::check_for_nan() 37 0.04628 0.04725 0.05061 0.37% MultiFab::contains_nan() 37 0.04616 0.04712 0.05048 0.37% Castro::initialize_advance() 1 0.04128 0.04178 0.04334 0.32% MLPoisson::define() 19 0.03201 0.03706 0.04223 0.31% MLMG::mgVcycle_down::4 128 0.01759 0.01958 0.04223 0.31% Castro::retry_advance_ctu() 11 0.03535 0.03624 0.03904 0.29% StateDataPhysBCFunct::() 187 0.0009869 0.01996 0.038 0.28% Amr::regrid() 1 0.03602 0.0364 0.03684 0.27% StateData::FillBoundary(geom) 193 0 0.01878 0.03674 0.27% Gravity::get_new_grav_vector() 18 0.03484 0.03545 0.0362 0.27% MLMG::mgVcycle_up::4 128 0.01559 0.01711 0.03372 0.25% MLMG::mgVcycle_down::5 128 0.01754 0.01887 0.03088 0.23% MLCellLinOp::defineAuxData() 19 0.02046 0.02229 0.02911 0.21% Castro::apply_source_to_state() 37 0.02661 0.02701 0.02742 0.20% MLMG::mgVcycle_up::5 128 0.01567 0.01677 0.02654 0.19% MLMG::prepareForSolve() 18 0.02048 0.02159 0.02491 0.18% VisMF::Write(FabArray) 4 0.01337 0.01453 0.0243 0.18% MLMG::computeMLResidual() 18 0.01052 0.0169 0.02361 0.17% Castro::sum_integrated_quantities() 2 0.02179 0.02196 0.02198 0.16% Castro::save_data_for_retry() 3 0.01841 0.01891 0.02128 0.16% MLMG::ResNormInf() 146 0.01227 0.0172 0.01974 0.14% FabArrayBase::getFB() 6380 0.01375 0.01472 0.01876 0.14% FabArrayBase::getCPC() 2045 0.01105 0.01241 0.01624 0.12% FabArray::setDomainBndry() 187 0.0006844 0.008568 0.01619 0.12% Castro::problem_post_timestep() 2 0.009895 0.01243 0.01559 0.11% Castro::wd_update() 2 0.009613 0.0124 0.01558 0.11% FabArray::norminf() 244 0.01037 0.01083 0.0146 0.11% Gravity::test_composite_phi() 1 0.01399 0.01406 0.01411 0.10% amrex::Add() 128 0.01341 0.01361 0.01386 0.10% MLLinOp::define() 19 0.004469 0.008002 0.01127 0.08% MLLinOp::defineGrids() 19 0.004379 0.007915 0.01119 0.08% MLCGSolver::ParallelAllReduce 189 0 0.0006299 0.01008 0.07% Other 2341 0.07933 0.09157 0.1267 0.93% -------------------------------------------------------------------------------------------- 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 397 KiB 403 KiB 410 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 2203 B 2230 B 2266 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 3 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 1581 KiB 1583 KiB 1583 KiB 2329 KiB 2329 KiB 2329 KiB Castro::construct_ctu_hydro_source() 6400 397 KiB 403 KiB 410 KiB 1176 KiB 1176 KiB 1176 KiB StateData::define() 192 1053 KiB 1053 KiB 1054 KiB 1164 KiB 1164 KiB 1164 KiB Castro::initMFs() 256 1032 KiB 1033 KiB 1033 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 672 KiB 672 KiB 672 KiB 803 KiB 803 KiB 803 KiB Amr::writePlotFile() 64 4523 B 4622 B 5165 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 50 KiB 51 KiB 52 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 13 KiB 14 KiB 16 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 433 KiB 433 KiB 433 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2369 B 4422 B 5742 B 225 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1435 B 1462 B 1485 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 18 KiB 19 KiB 19 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 130 KiB 130 KiB 131 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 128 KiB 128 KiB 128 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 10 KiB 10 KiB 12 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 1280 514 B 521 B 531 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 431 B 437 B 447 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 2204 B 2230 B 2266 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 19 B 21 B 39 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 33 KiB 33 KiB 33 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 53 B 53 B 54 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 115 B 116 B 119 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 147 B 200 B 292 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 31 KiB 31 KiB 31 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 238 B 243 B 249 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 1098 B 1114 B 1132 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 128 B 128 B 128 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 1766 B 1829 B 1858 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 0 B 0 B 1 B 16 KiB 16 KiB 16 KiB MLCellLinOp::defineAuxData() 11704 466 B 520 B 722 B 4608 B 5151 B 7152 B BndryData::define() 4864 484 B 485 B 486 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 25 B 35 B 45 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 0 B 23 B 37 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 27 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 7 B 9 B 648 B 648 B 648 B MLCGSolver::bicgstab 640 0 B 0 B 3 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 0 B 0 B 0 B 8 B 8 B 8 B ----------------------------------------------------------------------------------------------------------------------- AMReX (26.03-47-g384b99a08d34) finalized