Initializing AMReX (25.06-22-g1d307c2c661a-dirty)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (25.06-22-g1d307c2c661a-dirty) initialized Starting run at 07:14:43 UTC on 2025-08-04. Successfully read inputs file ... Castro git describe: 25.07-15-gb7dbecf76 AMReX git describe: 25.06-22-g1d307c2c6-dirty Microphysics git describe: joss_submit-1-ga350de6c 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.163378713 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001181565 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010693001 Gravity::multilevel_solve_for_new_phi() time = 0.10945574 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.01075037 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.049657458 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004151314 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 = -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 = -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, Invalid density = -415.2911996 at index 16, 183, 0 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index Invalid density = -764.1508347 at index 15, 184, 0 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = Invalid density = -859.0333271 at index 13, 186, 0 -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 1880 Invalid density = -6.239071148 at index 10, 66, 0 , 0 Invalid density = Invalid density = -907.0826581 at index 10, 67, 0 -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -6.239071148 at index 10, 189, 0 Invalid density = -606.7538937 at index 12, 68, 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 = Invalid density = -859.0333271 at index 13, 69, 0 -11.08517908 at index 8, 190, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -764.1508347 at index 15, 71, 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 Invalid density = -415.2911996 at index 16, 72, 0 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 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 Castro::construct_old_gravity() time = 0.00289958 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.064654719 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Castro::construct_old_gravity() time = 0.003776212 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.063871168 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.012034666 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.064259128 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.004493131 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.060861948 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.338566885 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.015922496 Gravity::solve_for_phi() time = 0.092906365 on level 0 Castro::construct_new_gravity() time = 0.095174964 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.061034174 on level 0 Castro::construct_old_gravity() time = 0.004177793 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.066702092 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.338190761 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.014372711 Gravity::solve_for_phi() time = 0.091205082 on level 0 Castro::construct_new_gravity() time = 0.093339092 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.060434915 on level 0 Castro::construct_old_gravity() time = 0.003565675 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.060106421 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.351422955 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.02214619 Gravity::solve_for_phi() time = 0.101898739 on level 0 Castro::construct_new_gravity() time = 0.104042899 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.060672158 on level 0 Castro::construct_old_gravity() time = 0.005677431 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.059514743 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.363688405 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.019818869 Gravity::solve_for_phi() time = 0.112626535 on level 0 Castro::construct_new_gravity() time = 0.114851181 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.072480705 on level 0 Castro::construct_old_gravity() time = 0.012683724 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.065051464 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.3575342 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.014933125 Gravity::solve_for_phi() time = 0.092271776 on level 0 Castro::construct_new_gravity() time = 0.094400317 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.060373214 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.003797413 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.059257126 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.346586495 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.014680284 Gravity::solve_for_phi() time = 0.097181009 on level 0 Castro::construct_new_gravity() time = 0.099324945 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.060425239 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.004191852 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.059782879 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.352085578 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.01490135 Gravity::solve_for_phi() time = 0.092980744 on level 0 Castro::construct_new_gravity() time = 0.095115837 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.059953312 on level 0 Castro::construct_old_gravity() time = 0.004283377 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.05941094 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.346232257 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.015103208 Gravity::solve_for_phi() time = 0.09374259 on level 0 Castro::construct_new_gravity() time = 0.095865791 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.060657422 on level 0 Castro::construct_old_gravity() time = 0.006184483 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.062460115 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.429390848 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.026514086 Gravity::solve_for_phi() time = 0.118769083 on level 0 Castro::construct_new_gravity() time = 0.121008232 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.065989533 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.006520007 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.060128458 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.351556406 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.016079682 Gravity::solve_for_phi() time = 0.093446496 on level 0 Castro::construct_new_gravity() time = 0.095586886 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.060410695 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.002850937 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.059655596 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.35347127 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.016379838 Gravity::solve_for_phi() time = 0.092515966 on level 0 Castro::construct_new_gravity() time = 0.0946883 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.061100246 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.006590371 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.059807422 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.346416728 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.015233692 Gravity::solve_for_phi() time = 0.094811226 on level 0 Castro::construct_new_gravity() time = 0.096951594 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.060604547 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.00481512 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.060823559 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.358890647 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.015466923 Gravity::solve_for_phi() time = 0.092234915 on level 0 Castro::construct_new_gravity() time = 0.094361611 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.059844397 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.00466075 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.062326656 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.354331005 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.015588403 Gravity::solve_for_phi() time = 0.092852256 on level 0 Castro::construct_new_gravity() time = 0.095051908 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.063074607 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.006153515 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.059263413 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.35208574 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.014974847 Gravity::solve_for_phi() time = 0.093209255 on level 0 Castro::construct_new_gravity() time = 0.095418637 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.060155051 on level 0 Castro::construct_old_gravity() time = 0.004447831 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.060034704 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.349086257 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.015788461 Gravity::solve_for_phi() time = 0.093751558 on level 0 Castro::construct_new_gravity() time = 0.095888849 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.062514413 on level 0 Castro::construct_old_gravity() time = 0.006466686 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.060068406 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.35466931 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002610153967 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.015930075 Gravity::solve_for_phi() time = 0.092272145 on level 0 Castro::construct_new_gravity() time = 0.094391455 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.060126496 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.581215356e+29 ener added: 1.581215356e+29 TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -4.341809003e+38 TIME= 0.25 YMOM = 1.164437349e+28 TIME= 0.25 ZMOM = 0 TIME= 0.25 ANG MOM X = 0 TIME= 0.25 ANG MOM Y = 0 TIME= 0.25 ANG MOM Z = -1.469212319e+38 TIME= 0.25 RHO*e = 1.325894512e+50 TIME= 0.25 RHO*K = 3.170941987e+49 TIME= 0.25 RHO*E = 1.642988711e+50 TIME= 0.25 RHO*PHI = -5.449011378e+50 TIME= 0.25 TOTAL ENERGY= -1.081516979e+50 TIME= 0.25 CENTER OF MASS X-LOC = 316424388.8 TIME= 0.25 CENTER OF MASS X-VEL = -170267.2489 TIME= 0.25 CENTER OF MASS Y-LOC = 2.427336677e-06 TIME= 0.25 CENTER OF MASS Y-VEL = 4.566427124e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 57600050.72 TIME= 0.25 MAXIMUM DENSITY = 4458760.788 TIME= 0.25 MAXIMUM T_S / T_E = 0 STEP = 1 TIME = 0.25 DT = 0.25 [STEP 1] Coarse TimeStep time: 12.61346661 [STEP 1] FAB kilobyte spread across MPI nodes: [6998 ... 6998] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.050332291 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004167147 seconds Ending run at 07:14:58 UTC on 2025-08-04. Run time = 14.87942596 Run time without initialization = 12.6690056 Average number of zones advanced per microsecond: 0.003 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 14.8712 ... 14.8737 ... 14.8796 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.7556 5.7903 5.9781 40.18% Amr::InitAmr() 1 1.301 1.304 1.309 8.80% Castro::computeTemp() 80 0.8475 0.9616 1.127 7.57% amrex::Copy() 724 0.8019 0.8082 0.8435 5.67% FabArray::ParallelCopy_finish() 1429 0.1913 0.6122 0.8082 5.43% MultiFab::AddProduct() 3439 0.5547 0.5571 0.5767 3.88% Castro::normalize_species() 79 0.5506 0.5548 0.5707 3.84% Castro::reset_internal_energy(Fab) 160 0.5236 0.5281 0.5398 3.63% MLPoisson::Fsmooth() 6144 0.3556 0.3858 0.4896 3.29% check_for_negative_density() 20 0.09296 0.3139 0.3981 2.68% Castro::initMFs() 1 0.3636 0.3686 0.3944 2.65% FabArray::setVal() 5013 0.35 0.3536 0.3713 2.50% MultiFab::Subtract() 159 0.3388 0.3406 0.3541 2.38% Gravity::fill_multipole_BCs() 19 0.2237 0.2767 0.3017 2.03% FabArray::ParallelCopy_nowait() 1429 0.2008 0.2174 0.2371 1.59% FabArray::sum() 3439 0.1985 0.1997 0.2059 1.38% Castro::initData() 1 0.003358 0.02801 0.1625 1.09% FillBoundary_finish() 6380 0.115 0.1426 0.1594 1.07% FillBoundary_nowait() 6380 0.09112 0.1212 0.1558 1.05% MLCellLinOp::applyBC() 7148 0.06407 0.08129 0.1337 0.90% TagBoxArray::collate() 2 0.0002677 0.02363 0.1269 0.85% MLPoisson::Fapply() 986 0.09542 0.09995 0.1221 0.82% Castro::construct_new_gravity_source() 17 0.07324 0.07428 0.07781 0.52% FabArray::Saxpy() 91 0.07204 0.0725 0.07435 0.50% FabArray::Xpay() 915 0.05184 0.05433 0.06159 0.41% Castro::construct_old_gravity_source() 20 0.04948 0.04973 0.05052 0.34% MultiFab::contains_nan() 37 0.04632 0.04683 0.04879 0.33% MLMG::addInterpCorrection() 768 0.03415 0.03673 0.04592 0.31% AmrLevel::derive() 18 0.04103 0.04156 0.0423 0.28% amrex::average_down 768 0.02533 0.02753 0.03529 0.24% StateData::FillBoundary(geom) 143 0 0.01767 0.03418 0.23% Castro::reset_internal_energy(MultiFab) 80 0.02348 0.02395 0.02502 0.17% VisMF::Write(FabArray) 4 0.01405 0.01508 0.02411 0.16% MLCellLinOp::defineAuxData() 19 0.01576 0.01682 0.02115 0.14% Castro::enforce_min_density() 79 0.01869 0.01892 0.01924 0.13% Amr::writePlotFile() 2 0.006765 0.01692 0.01855 0.12% FabArray::FillBoundary() 6380 0.014 0.0146 0.01776 0.12% MLCGSolver::bicgstab 8 0 0.001055 0.01688 0.11% Castro::construct_new_source() 102 0.0162 0.01629 0.01652 0.11% FabArray::norminf() 244 0.01079 0.0113 0.01597 0.11% FabArray::setDomainBndry() 147 0.0005007 0.007693 0.01574 0.11% Castro::evaluate_source_change() 196 0.01415 0.01447 0.01539 0.10% amrex::Add() 128 0.01457 0.01466 0.01499 0.10% Castro::wd_update() 2 0.005968 0.01093 0.0142 0.10% FabArrayBase::getFB() 6380 0.01067 0.0111 0.01373 0.09% FillPatchIterator::Initialize 147 0.01291 0.01307 0.01324 0.09% MLMG::ResNormInf() 146 0.006497 0.009806 0.01265 0.09% MLCellLinOp::smooth() 1544 0.01035 0.01086 0.0123 0.08% MLMG::prepareForSolve() 18 0.01025 0.01065 0.01202 0.08% Castro::construct_old_source() 120 0.0103 0.0104 0.01109 0.07% MLMG::solve() 18 0.0007977 0.006042 0.01105 0.07% Gravity::actual_solve_with_mlmg() 19 0.007372 0.007777 0.009512 0.06% FabArrayBase::getCPC() 2005 0.007784 0.007991 0.008572 0.06% Gravity::get_old_grav_vector() 20 0.006475 0.006621 0.007495 0.05% Gravity::get_new_grav_vector() 18 0.006022 0.006105 0.006359 0.04% MLCellLinOp::apply() 986 0.002073 0.002302 0.00501 0.03% Castro::initialize_do_advance() 20 0.00444 0.00469 0.004905 0.03% FabArray::ParallelCopy() 1429 0.003661 0.003815 0.00424 0.03% AmrLevel::FillPatch() 147 0.002383 0.002498 0.002811 0.02% Castro::sum_integrated_quantities() 2 0.00153 0.002066 0.002365 0.02% MLCellLinOp::correctionResidual() 768 0.001914 0.001982 0.002163 0.01% MLMG::mgVcycle() 128 0.001663 0.001744 0.001938 0.01% MLMG:computeResOfCorrection() 768 0.001519 0.00158 0.001716 0.01% Castro::save_data_for_retry() 3 0.001334 0.001366 0.001437 0.01% MLMG::actualBottomSolve() 128 0.0001325 0.000208 0.001258 0.01% StateDataPhysBCFunct::() 147 0.0007537 0.0008904 0.001039 0.01% Castro::subcycle_advance_ctu() 1 0.0002295 0.0002997 0.0009736 0.01% main() 1 0.0006979 0.0007745 0.0009384 0.01% FillPatchIterator::FillFromLevel0() 147 0.0007069 0.0007263 0.0007913 0.01% Gravity::solve_for_phi() 17 0.0006599 0.0006983 0.000756 0.01% Castro::Castro() 1 0.0004754 0.0005575 0.000652 0.00% MLMG::mgVcycle_down::2 128 0.0005418 0.0005919 0.000639 0.00% MLMG::mgVcycle_down::3 128 0.0005066 0.0005535 0.0006331 0.00% MLMG::mgVcycle_down::0 128 0.0005705 0.0005832 0.0006156 0.00% AmrMesh::MakeNewGrids() 2 0.0005389 0.0005637 0.0006152 0.00% MLMG::mgVcycle_down::1 128 0.0005605 0.0005749 0.0006119 0.00% MLMG::mgVcycle_down::4 128 0.0004886 0.0005082 0.0006081 0.00% MLMG::mgVcycle_down::5 128 0.0005129 0.0005327 0.0006025 0.00% Castro::clean_state() 79 0.0004322 0.0004668 0.0005366 0.00% FillPatchSingleLevel 147 0.0004759 0.0004908 0.0005248 0.00% MLCellLinOp::solutionResidual() 147 0.0004575 0.0004696 0.0004892 0.00% MLMG::oneIter() 128 0.0004027 0.0004204 0.0004674 0.00% Castro::do_advance_ctu() 20 0.0003994 0.0004234 0.0004562 0.00% MLMG::mgVcycle_up::0 128 0.0003765 0.0003915 0.0004473 0.00% MLMG::mgVcycle_up::1 128 0.0003717 0.0003904 0.0004286 0.00% MLMG::mgVcycle_up::2 128 0.0003262 0.0003575 0.0003942 0.00% MLMG::mgVcycle_up::3 128 0.0003223 0.0003379 0.0003883 0.00% MLMG::mgVcycle_up::4 128 0.0003261 0.0003377 0.0003747 0.00% MLMG::mgVcycle_up::5 128 0.0003135 0.0003243 0.0003559 0.00% Castro::do_old_sources() 20 0.0002868 0.0002972 0.000317 0.00% Castro::construct_new_gravity() 17 0.0001562 0.0001791 0.0003036 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002271 0.0002446 0.0002756 0.00% MLMG::computeResidual() 128 0.0002405 0.0002475 0.0002729 0.00% MLMG::mgVcycle_bottom 128 0.0002236 0.0002319 0.0002579 0.00% Castro::problem_post_timestep() 2 1.062e-05 2.795e-05 0.0002534 0.00% Castro::do_new_sources() 17 0.000214 0.0002289 0.0002494 0.00% TagBoxArray::mapPRD 2 0.0001241 0.0001428 0.0002143 0.00% MLPoisson::define() 19 0.0001349 0.0001586 0.0001908 0.00% Amr::coarseTimeStep() 1 0.0001377 0.0001655 0.0001703 0.00% Castro::check_for_nan() 37 0.0001106 0.0001194 0.0001381 0.00% Gravity::test_composite_phi() 1 0.0001055 0.0001149 0.0001375 0.00% Castro::retry_advance_ctu() 11 6.523e-05 7.879e-05 0.0001343 0.00% Castro::initialize_advance() 1 0.0001051 0.0001123 0.0001338 0.00% Castro::construct_old_gravity() 20 0.0001138 0.0001209 0.0001331 0.00% Amr::defBaseLevel() 1 8.792e-05 0.0001051 0.000132 0.00% Castro::apply_source_to_state() 37 9.827e-05 0.0001035 0.0001128 0.00% Amr::grid_places() 2 5.442e-05 6.177e-05 8.893e-05 0.00% MLLinOp::define() 19 7.355e-05 8.116e-05 8.707e-05 0.00% Castro::create_source_corrector() 17 5.805e-05 6.579e-05 8.356e-05 0.00% Amr::FinalizeInit() 1 1.531e-05 2.669e-05 8.159e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.882e-05 6.139e-05 6.525e-05 0.00% Castro::expand_state() 20 5.05e-05 5.415e-05 5.837e-05 0.00% MLMG::computeMLResidual() 18 4.047e-05 4.428e-05 5.705e-05 0.00% Castro::derive() 18 4.187e-05 4.39e-05 5.152e-05 0.00% Amr::regrid() 1 1.596e-05 2e-05 4.291e-05 0.00% Gravity::actual_multilevel_solve() 1 3.204e-05 3.317e-05 3.551e-05 0.00% Castro::advance() 1 7.047e-06 1.039e-05 3.178e-05 0.00% Amr::timeStep() 1 1.404e-05 2.02e-05 3.164e-05 0.00% Castro::post_regrid() 2 1.186e-05 1.651e-05 3.09e-05 0.00% Amr::init() 1 6.213e-06 8.461e-06 2.674e-05 0.00% Amr::bldFineLevels() 1 6.902e-06 9.587e-06 2.588e-05 0.00% Amr::InitializeInit() 1 7.886e-06 1.088e-05 2.556e-05 0.00% Castro::post_init() 1 1.212e-05 1.297e-05 1.593e-05 0.00% Castro::post_timestep() 1 5.857e-06 6.44e-06 7.132e-06 0.00% Amr::initialInit() 1 2.862e-06 3.258e-06 5.106e-06 0.00% Other 2521 0.05054 0.06281 0.09834 0.66% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.87 14.87 14.88 100.00% Amr::coarseTimeStep() 1 12.67 12.67 12.67 85.14% Amr::timeStep() 1 12.61 12.61 12.61 84.77% Castro::advance() 1 12.55 12.55 12.55 84.37% Castro::subcycle_advance_ctu() 1 12.5 12.5 12.5 84.03% Castro::do_advance_ctu() 20 12.45 12.46 12.46 83.76% Castro::construct_ctu_hydro_source() 20 6.854 6.891 6.929 46.57% Castro::clean_state() 79 3.688 3.815 4.012 26.96% Castro::computeTemp() 80 2.361 2.482 2.657 17.86% Gravity::solve_phi_with_mlmg() 19 1.674 1.731 1.757 11.81% Castro::construct_new_gravity() 17 1.591 1.648 1.673 11.25% Gravity::solve_for_phi() 17 1.554 1.611 1.637 11.00% Castro::reset_internal_energy(MultiFab) 80 1.512 1.521 1.567 10.53% Castro::initialize_do_advance() 20 1.416 1.457 1.493 10.04% Gravity::actual_solve_with_mlmg() 19 1.447 1.451 1.455 9.78% MLMG::solve() 18 1.389 1.396 1.398 9.39% Amr::InitAmr() 1 1.301 1.304 1.309 8.80% MLMG::oneIter() 128 1.248 1.25 1.253 8.42% MLMG::mgVcycle() 128 1.233 1.235 1.238 8.32% Castro::do_old_sources() 20 1.163 1.188 1.209 8.12% Castro::do_new_sources() 17 0.9262 0.9799 1.048 7.05% FabArray::ParallelCopy() 1429 0.4415 0.854 1.048 7.04% Castro::evaluate_source_change() 196 0.9664 0.9707 1.006 6.76% Amr::init() 1 0.901 0.901 0.901 6.06% MLCellLinOp::smooth() 1544 0.6375 0.6963 0.8471 5.69% Amr::initialInit() 1 0.8467 0.8467 0.8468 5.69% amrex::Copy() 724 0.8019 0.8082 0.8435 5.67% Castro::enforce_min_density() 79 0.8009 0.8054 0.8374 5.63% FabArray::ParallelCopy_finish() 1429 0.1913 0.6122 0.8082 5.43% AmrLevel::FillPatch() 147 0.5314 0.6524 0.711 4.78% Amr::InitializeInit() 1 0.4494 0.4805 0.6114 4.11% Amr::defBaseLevel() 1 0.4494 0.4805 0.6113 4.11% MultiFab::AddProduct() 3439 0.5547 0.5571 0.5767 3.88% Castro::normalize_species() 79 0.5506 0.5548 0.5707 3.84% Castro::reset_internal_energy(Fab) 160 0.5236 0.5281 0.5398 3.63% FillPatchIterator::Initialize 147 0.3388 0.468 0.5273 3.54% FillPatchIterator::FillFromLevel0() 147 0.3099 0.4473 0.5042 3.39% FillPatchSingleLevel 147 0.3092 0.4466 0.5035 3.38% MLPoisson::Fsmooth() 6144 0.3556 0.3858 0.4896 3.29% MLMG::mgVcycle_up::1 128 0.1099 0.2972 0.4852 3.26% MLCellLinOp::applyBC() 7148 0.3281 0.3715 0.4338 2.92% MLMG::addInterpCorrection() 768 0.06896 0.3409 0.4209 2.83% Castro::Castro() 1 0.3697 0.3748 0.402 2.70% check_for_negative_density() 20 0.09296 0.3139 0.3981 2.68% Amr::FinalizeInit() 1 0.2354 0.3662 0.3973 2.67% Castro::initMFs() 1 0.3638 0.3687 0.3946 2.65% FabArray::setVal() 5013 0.35 0.3536 0.3713 2.50% MultiFab::Subtract() 159 0.3388 0.3406 0.3541 2.38% FabArray::FillBoundary() 6380 0.2564 0.2902 0.3405 2.29% MLMG::mgVcycle_up::2 128 0.01388 0.1341 0.3167 2.13% Gravity::fill_multipole_BCs() 19 0.2258 0.2788 0.3038 2.04% MLMG::mgVcycle_down::0 128 0.2795 0.2823 0.2853 1.92% FabArray::ParallelCopy_nowait() 1429 0.2207 0.238 0.2585 1.74% Castro::expand_state() 20 0.1929 0.2338 0.2571 1.73% Castro::initData() 1 0.07789 0.104 0.2382 1.60% MLMG::mgVcycle_up::0 128 0.2151 0.2161 0.2172 1.46% MLCellLinOp::apply() 986 0.1597 0.172 0.2138 1.44% Amr::bldFineLevels() 1 0.046 0.1774 0.2079 1.40% FabArray::sum() 3439 0.1985 0.1997 0.2059 1.38% MLMG::mgVcycle_up::3 128 0.01383 0.02859 0.2045 1.37% FillBoundary_nowait() 6380 0.1022 0.133 0.1682 1.13% Amr::grid_places() 2 0.002609 0.1334 0.1645 1.11% AmrMesh::MakeNewGrids() 2 0.002553 0.1334 0.1645 1.11% TagBoxArray::mapPRD 2 0.0009347 0.1084 0.1627 1.09% MLMG:computeResOfCorrection() 768 0.1204 0.1308 0.1596 1.07% FillBoundary_finish() 6380 0.115 0.1426 0.1594 1.07% MLCellLinOp::correctionResidual() 768 0.1189 0.1293 0.1579 1.06% Castro::post_init() 1 0.1433 0.1451 0.1464 0.98% MLMG::mgVcycle_down::1 128 0.1302 0.1349 0.1397 0.94% Castro::construct_new_source() 102 0.1301 0.131 0.1367 0.92% TagBoxArray::collate() 2 0.0002677 0.02363 0.1269 0.85% MLPoisson::Fapply() 986 0.09542 0.09995 0.1221 0.82% Castro::create_source_corrector() 17 0.05317 0.08573 0.1121 0.75% Gravity::multilevel_solve_for_new_phi() 1 0.1094 0.1094 0.1095 0.74% Gravity::actual_multilevel_solve() 1 0.1093 0.1093 0.1094 0.74% MLCellLinOp::solutionResidual() 147 0.09157 0.09776 0.1047 0.70% Amr::writePlotFile() 2 0.1003 0.1004 0.1004 0.67% Castro::construct_old_gravity() 20 0.04501 0.0712 0.09077 0.61% Gravity::get_old_grav_vector() 20 0.0449 0.07108 0.09065 0.61% Castro::post_regrid() 2 0.08629 0.08756 0.08941 0.60% Castro::construct_old_source() 120 0.08353 0.08389 0.08495 0.57% MLMG::mgVcycle_down::2 128 0.01619 0.04851 0.08434 0.57% MLMG::computeResidual() 128 0.07824 0.07937 0.08066 0.54% Castro::construct_new_gravity_source() 17 0.07324 0.07428 0.07781 0.52% MLMG::mgVcycle_bottom 128 0.0003571 0.004993 0.07436 0.50% FabArray::Saxpy() 91 0.07204 0.0725 0.07435 0.50% MLMG::actualBottomSolve() 128 0.0001325 0.004761 0.07411 0.50% Castro::derive() 18 0.07028 0.07131 0.07394 0.50% AmrLevel::derive() 18 0.07024 0.07127 0.0739 0.50% MLCGSolver::bicgstab 8 0 0.004512 0.07219 0.49% MLMG::mgVcycle_down::3 128 0.01553 0.02267 0.07056 0.47% Castro::post_timestep() 1 0.05937 0.06333 0.06627 0.45% FabArray::Xpay() 915 0.05184 0.05433 0.06159 0.41% amrex::average_down 768 0.03456 0.04179 0.06128 0.41% Castro::initialize_advance() 1 0.04925 0.05019 0.05256 0.35% Castro::construct_old_gravity_source() 20 0.04948 0.04973 0.05052 0.34% Castro::check_for_nan() 37 0.04643 0.04695 0.04893 0.33% MultiFab::contains_nan() 37 0.04632 0.04683 0.04879 0.33% Amr::regrid() 1 0.04448 0.04509 0.04602 0.31% MLMG::mgVcycle_down::4 128 0.01538 0.01741 0.04352 0.29% MLPoisson::define() 19 0.03215 0.03743 0.04295 0.29% Castro::retry_advance_ctu() 11 0.03665 0.03741 0.03888 0.26% Gravity::get_new_grav_vector() 18 0.03625 0.03678 0.03884 0.26% StateDataPhysBCFunct::() 147 0.0007537 0.01856 0.03515 0.24% StateData::FillBoundary(geom) 143 0 0.01767 0.03418 0.23% MLMG::mgVcycle_up::4 128 0.01394 0.01538 0.03386 0.23% Castro::apply_source_to_state() 37 0.02937 0.02955 0.03002 0.20% MLMG::mgVcycle_down::5 128 0.01532 0.01648 0.02979 0.20% MLCellLinOp::defineAuxData() 19 0.02004 0.02185 0.02853 0.19% MLMG::prepareForSolve() 18 0.02067 0.02171 0.0257 0.17% MLMG::mgVcycle_up::5 128 0.01399 0.01489 0.02526 0.17% MLMG::computeMLResidual() 18 0.01141 0.01789 0.0252 0.17% VisMF::Write(FabArray) 4 0.01405 0.01508 0.02411 0.16% MLMG::ResNormInf() 146 0.01645 0.01939 0.02208 0.15% Castro::sum_integrated_quantities() 2 0.02159 0.02174 0.02176 0.15% Castro::save_data_for_retry() 3 0.01925 0.01975 0.02071 0.14% Castro::problem_post_timestep() 2 0.01028 0.01469 0.01763 0.12% Castro::wd_update() 2 0.01003 0.01466 0.01762 0.12% FabArray::norminf() 244 0.01079 0.0113 0.01597 0.11% FabArrayBase::getCPC() 2005 0.01019 0.01155 0.0159 0.11% FabArray::setDomainBndry() 147 0.0005007 0.007693 0.01574 0.11% FabArrayBase::getFB() 6380 0.01104 0.01185 0.0154 0.10% amrex::Add() 128 0.01457 0.01466 0.01499 0.10% Gravity::test_composite_phi() 1 0.01449 0.01456 0.01459 0.10% MLLinOp::define() 19 0.00522 0.008915 0.01266 0.09% Other 2521 0.08508 0.1025 0.1476 0.99% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Pinned Memory Usage: --------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max --------------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 16 9 B 14 B 28 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- Cpu Memory Usage: ------------------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max ------------------------------------------------------------------------------------------------------------------------------- Castro::save_data_for_retry() 416 416 1584 KiB 1585 KiB 1585 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 406 KiB 408 KiB 422 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1047 KiB 1047 KiB 1048 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 996 KiB 997 KiB 997 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 670 KiB 670 KiB 670 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 4209 B 4297 B 4862 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 49 KiB 49 KiB 51 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 11 KiB 13 KiB 14 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 419 KiB 419 KiB 420 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 1973 B 4073 B 5070 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1310 B 1328 B 1380 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 17 KiB 17 KiB 17 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 133 KiB 133 KiB 133 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 192 129 KiB 129 KiB 129 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 6102 9810 B 10 KiB 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 254 B 255 B 262 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 429 B 431 B 447 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 18 B 20 B 36 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 33 KiB 33 KiB 33 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 51 B 51 B 51 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 3 B 5 B 7 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 110 B 111 B 117 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 137 B 218 B 279 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 64 32 KiB 32 KiB 32 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 2432 229 B 236 B 241 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1059 B 1063 B 1101 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 120 B 120 B 120 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1711 B 1773 B 1801 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 32 0 B 0 B 1 B 16 KiB 16 KiB 16 KiB MLCellLinOp::defineAuxData() 11704 11704 443 B 495 B 687 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 461 B 462 B 463 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 22 B 34 B 47 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 18 B 28 B 2916 B 2916 B 2916 B OwnerMask() 64 64 0 B 0 B 0 B 2592 B 2592 B 2592 B MLMG::addInterpCorrection() 5632 5632 1 B 16 B 26 B 1024 B 1024 B 1024 B amrex::average_down 5632 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 64 0 B 5 B 7 B 648 B 648 B 648 B MLCGSolver::bicgstab 640 640 0 B 0 B 2 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 114375 0 B 0 B 0 B 8 B 8 B 8 B ------------------------------------------------------------------------------------------------------------------------------- AMReX (25.06-22-g1d307c2c661a-dirty) finalized