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 06:48:30 UTC on 2025-07-22. Successfully read inputs file ... Castro git describe: 25.07-12-g2c8f2fc1e AMReX git describe: 25.06-22-g1d307c2c6-dirty Microphysics git describe: joss_submit 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.169654506 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.00118224 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010697653 Gravity::multilevel_solve_for_new_phi() time = 0.107593775 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010496345 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.04897429 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004145436 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 = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071148 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Invalid density = -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 Castro::construct_old_gravity() time = 0.002715465 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.065188362 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.003501466 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.06414597 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.004002923 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.06251544 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.003187643 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.060133098 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.336290918 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.016335561 Gravity::solve_for_phi() time = 0.096160191 on level 0 Castro::construct_new_gravity() time = 0.098356274 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.060780414 on level 0 Castro::construct_old_gravity() time = 0.002737399 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.060226177 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.336581372 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.01503933 Gravity::solve_for_phi() time = 0.091357047 on level 0 Castro::construct_new_gravity() time = 0.093948274 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.060543091 on level 0 Castro::construct_old_gravity() time = 0.002802914 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.060035884 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.342463506 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.015106467 Gravity::solve_for_phi() time = 0.090929306 on level 0 Castro::construct_new_gravity() time = 0.09305753 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.06058471 on level 0 Castro::construct_old_gravity() time = 0.002768156 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.059727397 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.344397458 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.015186566 Gravity::solve_for_phi() time = 0.091816359 on level 0 Castro::construct_new_gravity() time = 0.093952716 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.059726716 on level 0 Castro::construct_old_gravity() time = 0.004224047 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.059200272 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.347892187 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.014743086 Gravity::solve_for_phi() time = 0.090744922 on level 0 Castro::construct_new_gravity() time = 0.092866989 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.059570125 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.002763951 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.059459346 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.342276876 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.014434602 Gravity::solve_for_phi() time = 0.092386633 on level 0 Castro::construct_new_gravity() time = 0.094506385 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.060485436 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.002763277 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.059505624 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.347672781 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.015248857 Gravity::solve_for_phi() time = 0.09202207 on level 0 Castro::construct_new_gravity() time = 0.094157045 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.059859056 on level 0 Castro::construct_old_gravity() time = 0.003530559 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.059562941 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.34304911 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.014751379 Gravity::solve_for_phi() time = 0.092096795 on level 0 Castro::construct_new_gravity() time = 0.094226047 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.059758705 on level 0 Castro::construct_old_gravity() time = 0.002815694 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.05921344 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.347484711 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.014845044 Gravity::solve_for_phi() time = 0.096858804 on level 0 Castro::construct_new_gravity() time = 0.098986654 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.061358108 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.002790494 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.060831528 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.349400545 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.015625947 Gravity::solve_for_phi() time = 0.093142983 on level 0 Castro::construct_new_gravity() time = 0.095279041 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.061167018 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.002825992 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.060859553 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.354637866 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.016110463 Gravity::solve_for_phi() time = 0.093887175 on level 0 Castro::construct_new_gravity() time = 0.0960085 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.061181074 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.002873211 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.060664392 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.349082182 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.016146284 Gravity::solve_for_phi() time = 0.093418826 on level 0 Castro::construct_new_gravity() time = 0.09562609 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.061318821 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.002843578 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.060727816 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.353752766 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.016639951 Gravity::solve_for_phi() time = 0.093757675 on level 0 Castro::construct_new_gravity() time = 0.095932842 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.061323138 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.003590983 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.061089738 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.349548364 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.016254501 Gravity::solve_for_phi() time = 0.093208475 on level 0 Castro::construct_new_gravity() time = 0.095361981 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.061394811 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.00283632 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.060765885 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.371256517 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.016573992 Gravity::solve_for_phi() time = 0.093812177 on level 0 Castro::construct_new_gravity() time = 0.095933966 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.061373836 on level 0 Castro::construct_old_gravity() time = 0.00282947 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.060770086 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.350126518 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.016952316 Gravity::solve_for_phi() time = 0.094236429 on level 0 Castro::construct_new_gravity() time = 0.096404398 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.061324503 on level 0 Castro::construct_old_gravity() time = 0.003135113 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.060745076 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.354853672 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002669808887 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.016147303 Gravity::solve_for_phi() time = 0.095073015 on level 0 Castro::construct_new_gravity() time = 0.09720912 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.061302468 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.33508398 [STEP 1] FAB kilobyte spread across MPI nodes: [6998 ... 6998] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.049512427 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004173399 seconds Ending run at 06:48:45 UTC on 2025-07-22. Run time = 14.96903587 Run time without initialization = 12.38965383 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.9611 ... 14.9634 ... 14.9691 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.7039 5.7717 5.8558 39.12% Amr::InitAmr() 1 1.674 1.676 1.682 11.24% Castro::computeTemp() 80 0.8525 0.9601 1.131 7.55% amrex::Copy() 724 0.7976 0.8048 0.813 5.43% FabArray::ParallelCopy_finish() 1429 0.2075 0.5559 0.7842 5.24% MultiFab::AddProduct() 3439 0.551 0.5606 0.6267 4.19% Castro::normalize_species() 79 0.5456 0.5514 0.5593 3.74% Castro::reset_internal_energy(Fab) 160 0.5185 0.5272 0.5379 3.59% MLPoisson::Fsmooth() 6144 0.3544 0.3846 0.4653 3.11% Castro::initMFs() 1 0.3572 0.3595 0.3676 2.46% FabArray::setVal() 5013 0.3468 0.3516 0.3597 2.40% MultiFab::Subtract() 159 0.3371 0.3399 0.3441 2.30% Gravity::fill_multipole_BCs() 19 0.1986 0.2573 0.2828 1.89% check_for_negative_density() 20 0.05748 0.1642 0.2602 1.74% FabArray::ParallelCopy_nowait() 1429 0.2022 0.2164 0.2351 1.57% FabArray::sum() 3439 0.1971 0.1991 0.2009 1.34% Castro::initData() 1 0.003354 0.02804 0.1651 1.10% TagBoxArray::collate() 2 0.0001981 0.03105 0.1649 1.10% FillBoundary_finish() 6380 0.1042 0.129 0.1566 1.05% FillBoundary_nowait() 6380 0.09329 0.1207 0.1531 1.02% MLCellLinOp::applyBC() 7148 0.06455 0.08069 0.1282 0.86% MLPoisson::Fapply() 986 0.09463 0.09927 0.115 0.77% Castro::construct_new_gravity_source() 17 0.07329 0.07411 0.07489 0.50% FabArray::Saxpy() 91 0.07199 0.07268 0.07325 0.49% FabArray::Xpay() 915 0.05142 0.05414 0.05993 0.40% Castro::construct_old_gravity_source() 20 0.04936 0.04987 0.05055 0.34% MultiFab::contains_nan() 37 0.0459 0.04651 0.04749 0.32% MLMG::addInterpCorrection() 768 0.03381 0.03656 0.04446 0.30% AmrLevel::derive() 18 0.0417 0.04195 0.04243 0.28% amrex::average_down 768 0.02502 0.02731 0.03423 0.23% StateData::FillBoundary(geom) 143 0 0.01757 0.0341 0.23% VisMF::Write(FabArray) 4 0.01421 0.01521 0.02461 0.16% Castro::reset_internal_energy(MultiFab) 80 0.02311 0.02357 0.02411 0.16% MLCellLinOp::defineAuxData() 19 0.01556 0.01665 0.02058 0.14% Castro::enforce_min_density() 79 0.01855 0.01883 0.01908 0.13% FabArray::FillBoundary() 6380 0.01387 0.01461 0.01801 0.12% Castro::construct_new_source() 102 0.01609 0.01634 0.01736 0.12% MLCGSolver::bicgstab 8 0 0.001078 0.01725 0.12% Amr::writePlotFile() 2 0.006792 0.01472 0.01549 0.10% amrex::Add() 128 0.01449 0.01476 0.01526 0.10% Castro::evaluate_source_change() 196 0.01378 0.01422 0.01492 0.10% FabArray::setDomainBndry() 147 0.0005004 0.007639 0.01457 0.10% FabArray::norminf() 244 0.01068 0.01109 0.01418 0.09% FabArrayBase::getFB() 6380 0.01057 0.01107 0.01358 0.09% Castro::wd_update() 2 0.005481 0.009821 0.01338 0.09% FillPatchIterator::Initialize 147 0.01274 0.013 0.01327 0.09% MLCellLinOp::smooth() 1544 0.01025 0.01086 0.01228 0.08% MLMG::prepareForSolve() 18 0.01014 0.01066 0.01185 0.08% Castro::construct_old_source() 120 0.01026 0.01041 0.01115 0.07% MLMG::solve() 18 0.0008916 0.004322 0.009267 0.06% Gravity::actual_solve_with_mlmg() 19 0.007297 0.007737 0.009141 0.06% MLMG::ResNormInf() 146 0.003924 0.006932 0.008936 0.06% FabArrayBase::getCPC() 2005 0.007664 0.007906 0.008301 0.06% Gravity::get_old_grav_vector() 20 0.006351 0.006567 0.00684 0.05% Gravity::get_new_grav_vector() 18 0.006018 0.006162 0.006535 0.04% MLCellLinOp::apply() 986 0.002044 0.002295 0.005017 0.03% Castro::initialize_do_advance() 20 0.004574 0.004681 0.004935 0.03% FabArray::ParallelCopy() 1429 0.003605 0.0038 0.004212 0.03% AmrLevel::FillPatch() 147 0.002381 0.002466 0.002694 0.02% Castro::sum_integrated_quantities() 2 0.001894 0.002383 0.002569 0.02% MLCellLinOp::correctionResidual() 768 0.001896 0.001962 0.002111 0.01% MLMG::mgVcycle() 128 0.001642 0.001711 0.00187 0.01% MLMG:computeResOfCorrection() 768 0.001515 0.001571 0.001668 0.01% Castro::save_data_for_retry() 3 0.00129 0.001351 0.001466 0.01% MLMG::actualBottomSolve() 128 0.0001306 0.0002014 0.001214 0.01% Castro::subcycle_advance_ctu() 1 0.0002134 0.000289 0.001013 0.01% StateDataPhysBCFunct::() 147 0.0007187 0.0008719 0.0009701 0.01% FillPatchIterator::FillFromLevel0() 147 0.000697 0.0007216 0.0007728 0.01% Gravity::solve_for_phi() 17 0.0006527 0.0006867 0.0007463 0.00% main() 1 0.0006608 0.0007048 0.0007446 0.00% MLMG::mgVcycle_down::3 128 0.0005076 0.000549 0.0006528 0.00% MLMG::mgVcycle_down::2 128 0.0005406 0.0005828 0.00062 0.00% AmrMesh::MakeNewGrids() 2 0.0005374 0.0005689 0.0006135 0.00% MLMG::mgVcycle_down::4 128 0.0004854 0.0005102 0.0006081 0.00% MLMG::mgVcycle_down::1 128 0.0005532 0.0005709 0.0006067 0.00% MLMG::mgVcycle_down::5 128 0.0005119 0.0005275 0.0006038 0.00% MLMG::mgVcycle_down::0 128 0.0005577 0.0005733 0.0005992 0.00% FillPatchSingleLevel 147 0.0004697 0.0004845 0.0005217 0.00% Castro::clean_state() 79 0.000428 0.0004524 0.0004936 0.00% Castro::Castro() 1 0.0003981 0.0004307 0.0004885 0.00% MLCellLinOp::solutionResidual() 147 0.0004529 0.000467 0.0004857 0.00% MLMG::oneIter() 128 0.0003901 0.0004058 0.0004528 0.00% Castro::do_advance_ctu() 20 0.0003817 0.0004079 0.0004392 0.00% MLMG::mgVcycle_up::0 128 0.000376 0.0003913 0.0004326 0.00% MLMG::mgVcycle_up::1 128 0.000368 0.0003872 0.0004269 0.00% MLMG::mgVcycle_up::2 128 0.0003263 0.0003544 0.0003985 0.00% MLMG::mgVcycle_up::4 128 0.0003234 0.0003355 0.0003805 0.00% MLMG::mgVcycle_up::3 128 0.0003215 0.0003364 0.000378 0.00% MLMG::mgVcycle_up::5 128 0.0003122 0.0003224 0.0003564 0.00% Castro::problem_post_timestep() 2 1.119e-05 3.334e-05 0.0003391 0.00% Castro::do_old_sources() 20 0.0002782 0.000294 0.0003101 0.00% Castro::construct_new_gravity() 17 0.0001624 0.0001797 0.0002892 0.00% MLMG::computeResidual() 128 0.0002387 0.0002471 0.0002746 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002195 0.0002359 0.0002718 0.00% MLMG::mgVcycle_bottom 128 0.0002215 0.0002295 0.0002455 0.00% Castro::do_new_sources() 17 0.0002073 0.0002221 0.0002405 0.00% MLPoisson::define() 19 0.0001339 0.0001495 0.0001815 0.00% Amr::coarseTimeStep() 1 0.0001323 0.0001614 0.0001677 0.00% Castro::retry_advance_ctu() 11 6.243e-05 7.817e-05 0.0001676 0.00% Castro::check_for_nan() 37 0.0001029 0.0001165 0.000145 0.00% Gravity::test_composite_phi() 1 0.0001051 0.0001098 0.0001401 0.00% TagBoxArray::mapPRD 2 0.0001207 0.0001257 0.0001334 0.00% Castro::construct_old_gravity() 20 0.000109 0.0001153 0.000127 0.00% Amr::defBaseLevel() 1 9.601e-05 0.0001022 0.000116 0.00% Castro::initialize_advance() 1 0.0001031 0.0001083 0.0001141 0.00% Amr::grid_places() 2 5.252e-05 5.884e-05 0.0001127 0.00% Castro::apply_source_to_state() 37 9.544e-05 0.000102 0.0001109 0.00% Amr::FinalizeInit() 1 8.075e-06 1.39e-05 8.95e-05 0.00% MLLinOp::define() 19 7.539e-05 7.945e-05 8.688e-05 0.00% Castro::create_source_corrector() 17 5.119e-05 6.05e-05 6.77e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.013e-05 6.185e-05 6.547e-05 0.00% Castro::expand_state() 20 5.062e-05 5.357e-05 5.937e-05 0.00% Amr::regrid() 1 1.543e-05 2.018e-05 5.306e-05 0.00% Castro::derive() 18 4.34e-05 4.55e-05 5.058e-05 0.00% MLMG::computeMLResidual() 18 3.978e-05 4.33e-05 4.675e-05 0.00% Gravity::actual_multilevel_solve() 1 3.265e-05 3.344e-05 3.673e-05 0.00% Amr::timeStep() 1 1.492e-05 1.76e-05 2.308e-05 0.00% Castro::post_regrid() 2 1.053e-05 1.2e-05 2.143e-05 0.00% Amr::InitializeInit() 1 8.935e-06 1.223e-05 1.666e-05 0.00% Castro::advance() 1 6.958e-06 8.742e-06 1.634e-05 0.00% Castro::post_init() 1 1.189e-05 1.275e-05 1.566e-05 0.00% Amr::bldFineLevels() 1 6.713e-06 7.782e-06 1.084e-05 0.00% Amr::init() 1 6.111e-06 6.984e-06 8.364e-06 0.00% Castro::post_timestep() 1 6.345e-06 7.124e-06 7.978e-06 0.00% Amr::initialInit() 1 2.863e-06 3.27e-06 4.925e-06 0.00% Other 2521 0.04972 0.06171 0.09538 0.64% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.96 14.96 14.97 100.00% Amr::coarseTimeStep() 1 12.39 12.39 12.39 82.76% Amr::timeStep() 1 12.34 12.34 12.34 82.40% Castro::advance() 1 12.27 12.27 12.27 82.00% Castro::subcycle_advance_ctu() 1 12.22 12.22 12.22 81.66% Castro::do_advance_ctu() 20 12.18 12.18 12.19 81.41% Castro::construct_ctu_hydro_source() 20 6.683 6.72 6.781 45.30% Castro::clean_state() 79 3.67 3.808 4.006 26.76% Castro::computeTemp() 80 2.354 2.478 2.662 17.79% Gravity::solve_phi_with_mlmg() 19 1.612 1.674 1.7 11.36% Amr::InitAmr() 1 1.674 1.676 1.682 11.24% Castro::construct_new_gravity() 17 1.531 1.593 1.618 10.81% Gravity::solve_for_phi() 17 1.493 1.556 1.582 10.57% Castro::reset_internal_energy(MultiFab) 80 1.499 1.518 1.537 10.27% Castro::initialize_do_advance() 20 1.394 1.43 1.466 9.79% Gravity::actual_solve_with_mlmg() 19 1.409 1.413 1.418 9.47% MLMG::solve() 18 1.353 1.36 1.361 9.10% MLMG::oneIter() 128 1.218 1.22 1.222 8.17% Castro::do_old_sources() 20 1.148 1.186 1.215 8.12% MLMG::mgVcycle() 128 1.203 1.205 1.207 8.06% Castro::evaluate_source_change() 196 0.9596 0.9728 1.036 6.92% Castro::do_new_sources() 17 0.93 0.9801 1.033 6.90% FabArray::ParallelCopy() 1429 0.4518 0.7967 1.022 6.83% Amr::init() 1 0.8972 0.8972 0.8973 5.99% Amr::initialInit() 1 0.8436 0.8436 0.8437 5.64% Castro::enforce_min_density() 79 0.7961 0.8054 0.8342 5.57% MLCellLinOp::smooth() 1544 0.6256 0.683 0.8275 5.53% amrex::Copy() 724 0.7976 0.8048 0.813 5.43% FabArray::ParallelCopy_finish() 1429 0.2075 0.5559 0.7842 5.24% AmrLevel::FillPatch() 147 0.5054 0.6054 0.6925 4.63% MultiFab::AddProduct() 3439 0.551 0.5606 0.6267 4.19% Amr::InitializeInit() 1 0.443 0.4701 0.6114 4.08% Amr::defBaseLevel() 1 0.443 0.4701 0.6113 4.08% Castro::normalize_species() 79 0.5456 0.5514 0.5593 3.74% Castro::reset_internal_energy(Fab) 160 0.5185 0.5272 0.5379 3.59% FillPatchIterator::Initialize 147 0.3196 0.422 0.511 3.41% FillPatchIterator::FillFromLevel0() 147 0.2917 0.4014 0.4886 3.26% FillPatchSingleLevel 147 0.291 0.4006 0.4879 3.26% MLMG::mgVcycle_up::1 128 0.108 0.2863 0.4666 3.12% MLPoisson::Fsmooth() 6144 0.3544 0.3846 0.4653 3.11% MLCellLinOp::applyBC() 7148 0.3169 0.3568 0.4389 2.93% MLMG::addInterpCorrection() 768 0.06692 0.3268 0.4051 2.71% Amr::FinalizeInit() 1 0.2322 0.3735 0.4006 2.68% Castro::Castro() 1 0.3628 0.365 0.373 2.49% Castro::initMFs() 1 0.3574 0.3596 0.3677 2.46% FabArray::setVal() 5013 0.3468 0.3516 0.3597 2.40% MultiFab::Subtract() 159 0.3371 0.3399 0.3441 2.30% FabArray::FillBoundary() 6380 0.2383 0.2762 0.3246 2.17% MLMG::mgVcycle_up::2 128 0.01373 0.1288 0.3038 2.03% Gravity::fill_multipole_BCs() 19 0.2008 0.2594 0.2849 1.90% MLMG::mgVcycle_down::0 128 0.2763 0.2785 0.2817 1.88% check_for_negative_density() 20 0.05748 0.1642 0.2602 1.74% FabArray::ParallelCopy_nowait() 1429 0.2224 0.237 0.2567 1.71% Castro::expand_state() 20 0.1935 0.225 0.2531 1.69% Castro::initData() 1 0.07839 0.1034 0.2412 1.61% MLMG::mgVcycle_up::0 128 0.213 0.2137 0.2142 1.43% Amr::bldFineLevels() 1 0.04652 0.1873 0.2141 1.43% MLCellLinOp::apply() 986 0.1567 0.1689 0.2071 1.38% FabArray::sum() 3439 0.1971 0.1991 0.2009 1.34% MLMG::mgVcycle_up::3 128 0.01371 0.028 0.196 1.31% Amr::grid_places() 2 0.002458 0.1437 0.1708 1.14% AmrMesh::MakeNewGrids() 2 0.002406 0.1437 0.1708 1.14% TagBoxArray::mapPRD 2 0.0008357 0.1113 0.1686 1.13% FillBoundary_nowait() 6380 0.1044 0.1325 0.1654 1.10% TagBoxArray::collate() 2 0.0001981 0.03105 0.1649 1.10% FillBoundary_finish() 6380 0.1042 0.129 0.1566 1.05% MLMG:computeResOfCorrection() 768 0.1188 0.1291 0.1538 1.03% MLCellLinOp::correctionResidual() 768 0.1173 0.1275 0.1521 1.02% Castro::post_init() 1 0.1413 0.1429 0.1434 0.96% MLMG::mgVcycle_down::1 128 0.1261 0.1301 0.1347 0.90% Castro::construct_new_source() 102 0.1296 0.131 0.1329 0.89% MLPoisson::Fapply() 986 0.09463 0.09927 0.115 0.77% Gravity::multilevel_solve_for_new_phi() 1 0.1075 0.1075 0.1076 0.72% Gravity::actual_multilevel_solve() 1 0.1075 0.1075 0.1075 0.72% MLCellLinOp::solutionResidual() 147 0.08954 0.09617 0.1041 0.70% Amr::writePlotFile() 2 0.09887 0.09897 0.09899 0.66% Castro::create_source_corrector() 17 0.05199 0.0699 0.09223 0.62% Castro::post_regrid() 2 0.08637 0.08681 0.08843 0.59% Castro::construct_old_source() 120 0.08334 0.0841 0.08543 0.57% MLMG::mgVcycle_down::2 128 0.01604 0.04708 0.08093 0.54% MLMG::computeResidual() 128 0.07807 0.07846 0.0798 0.53% Castro::construct_new_gravity_source() 17 0.07329 0.07411 0.07489 0.50% MLMG::mgVcycle_bottom 128 0.000352 0.004974 0.07415 0.50% MLMG::actualBottomSolve() 128 0.0001306 0.004745 0.07391 0.49% FabArray::Saxpy() 91 0.07199 0.07268 0.07325 0.49% Castro::derive() 18 0.07159 0.07187 0.07253 0.48% AmrLevel::derive() 18 0.07154 0.07183 0.07248 0.48% MLCGSolver::bicgstab 8 0 0.004502 0.07204 0.48% MLMG::mgVcycle_down::3 128 0.01536 0.02226 0.06816 0.46% Castro::post_timestep() 1 0.06077 0.06377 0.06687 0.45% FabArray::Xpay() 915 0.05142 0.05414 0.05993 0.40% amrex::average_down 768 0.03403 0.04111 0.05873 0.39% Castro::construct_old_gravity() 20 0.04484 0.05331 0.05798 0.39% Gravity::get_old_grav_vector() 20 0.04473 0.0532 0.05786 0.39% Castro::construct_old_gravity_source() 20 0.04936 0.04987 0.05055 0.34% Castro::initialize_advance() 1 0.04878 0.04905 0.04983 0.33% Castro::check_for_nan() 37 0.04601 0.04663 0.04764 0.32% MultiFab::contains_nan() 37 0.0459 0.04651 0.04749 0.32% Amr::regrid() 1 0.04441 0.04471 0.04578 0.31% MLPoisson::define() 19 0.03104 0.03636 0.04176 0.28% MLMG::mgVcycle_down::4 128 0.0153 0.01727 0.04075 0.27% Castro::retry_advance_ctu() 11 0.03604 0.03691 0.03823 0.26% Gravity::get_new_grav_vector() 18 0.03637 0.03708 0.03773 0.25% StateDataPhysBCFunct::() 147 0.0007187 0.01844 0.03504 0.23% StateData::FillBoundary(geom) 143 0 0.01757 0.0341 0.23% MLMG::mgVcycle_up::4 128 0.01377 0.01523 0.03208 0.21% Castro::apply_source_to_state() 37 0.0293 0.02959 0.02984 0.20% MLMG::mgVcycle_down::5 128 0.01531 0.01644 0.02844 0.19% MLCellLinOp::defineAuxData() 19 0.01976 0.02168 0.02789 0.19% MLMG::prepareForSolve() 18 0.02053 0.02166 0.02476 0.17% VisMF::Write(FabArray) 4 0.01421 0.01521 0.02461 0.16% MLMG::mgVcycle_up::5 128 0.01384 0.0148 0.02445 0.16% MLMG::computeMLResidual() 18 0.01096 0.01718 0.02354 0.16% Castro::sum_integrated_quantities() 2 0.02202 0.02226 0.02228 0.15% Castro::save_data_for_retry() 3 0.01912 0.0194 0.02042 0.14% MLMG::ResNormInf() 146 0.01335 0.01641 0.0183 0.12% Castro::problem_post_timestep() 2 0.01054 0.01361 0.01682 0.11% Castro::wd_update() 2 0.01053 0.01358 0.01681 0.11% FabArrayBase::getCPC() 2005 0.01003 0.01146 0.01549 0.10% amrex::Add() 128 0.01449 0.01476 0.01526 0.10% FabArrayBase::getFB() 6380 0.01092 0.01181 0.01519 0.10% FabArray::setDomainBndry() 147 0.0005004 0.007639 0.01457 0.10% Gravity::test_composite_phi() 1 0.01422 0.01425 0.01427 0.10% FabArray::norminf() 244 0.01068 0.01109 0.01418 0.09% MLLinOp::define() 19 0.004065 0.008074 0.01166 0.08% Other 2521 0.0829 0.1001 0.142 0.95% -------------------------------------------------------------------------------------------- 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 12 B 21 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 1539 KiB 1540 KiB 1540 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 400 KiB 405 KiB 411 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1019 KiB 1019 KiB 1020 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 970 KiB 970 KiB 970 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 651 KiB 651 KiB 652 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 4270 B 4312 B 4755 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 48 KiB 49 KiB 50 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 11 KiB 12 KiB 14 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 407 KiB 408 KiB 408 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 1895 B 3755 B 4803 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1328 B 1333 B 1347 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 16 KiB 17 KiB 17 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 129 KiB 129 KiB 129 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 192 126 KiB 126 KiB 126 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 6102 9494 B 9 KiB 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 252 B 254 B 262 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 425 B 429 B 440 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 19 B 20 B 36 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 32 KiB 32 KiB 32 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 49 B 49 B 50 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 3 B 4 B 5 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 109 B 111 B 113 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 136 B 162 B 177 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 64 31 KiB 31 KiB 31 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 2432 227 B 231 B 238 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1045 B 1059 B 1128 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 117 B 117 B 117 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1634 B 1702 B 1730 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 429 B 480 B 666 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 446 B 448 B 449 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 23 B 32 B 42 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 19 B 29 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 25 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