Initializing AMReX (26.10-97-gf5d545d2f0b0)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.10-97-gf5d545d2f0b0) initialized Starting run at 07:50:40 UTC on 2026-10-22. Successfully read inputs file ... Castro git describe: 26.10 AMReX git describe: 26.10-97-gf5d545d2f Microphysics git describe: 26.10-9-gffa96fc3 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.177929586 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.002313561 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010152298 Gravity::multilevel_solve_for_new_phi() time = 0.116976329 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010176544 TIME= 0 MASS = 2.549996568e+33 TIME= 0 XMOM = 0 TIME= 0 YMOM = 0 TIME= 0 ZMOM = 0 TIME= 0 ANG MOM X = 0 TIME= 0 ANG MOM Y = 0 TIME= 0 ANG MOM Z = 0 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 = 1.213668339e-07 TIME= 0 CENTER OF MASS Y-VEL = 0 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.04705001 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004197113 seconds [Level 0 step 1] ADVANCE at time 0 with dt = 0.25  Beginning subcycle 1 starting at time 0 with dt = 0.25 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.7670713549 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.7670713549 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071148 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.769409 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477044 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035502 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -0.01295059363 at index 21, 83, 0 Invalid density = -0.01295059363 at index 21, 172, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035502 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477044 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.769409 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071148 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Castro::construct_old_gravity() time = 0.003549339 on level 0 Castro::do_old_sources() time = 0.008862018 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.505048496e+28 zmom added: 0 eint added: 0 ener added: 6.931398274e+47 Castro::do_old_sources() time = 0.057529058 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Castro::construct_old_gravity() time = 0.005748176 on level 0 Castro::do_old_sources() time = 0.007829241 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.252524248e+28 zmom added: 0 eint added: 0 ener added: 3.465699137e+47 Castro::do_old_sources() time = 0.056594422 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.125 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.0625   Beginning subcycle 1 starting at time 0 with dt = 0.0625 Estimated number of subcycles remaining: 4 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.1060744788 at index 20, 78, 0 Invalid density = -0.1060744788 at index 20, 177, 0 Castro::construct_old_gravity() time = 0.004283439 on level 0 Castro::do_old_sources() time = 0.008595037 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -2.93534348e+40 ymom added: 1.626262124e+28 zmom added: 0 eint added: 0 ener added: 1.732849568e+47 Castro::do_old_sources() time = 0.053886913 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.004567397 on level 0 Castro::do_old_sources() time = 0.010125969 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.131310621e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.052549933 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.302115204 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.013608326 Gravity::solve_for_phi() time = 0.096028283 on level 0 Castro::construct_new_gravity() time = 0.098109424 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.588996166e+27 zmom added: 0 eint added: 0 ener added: 6.222804276e+43 Castro::do_new_sources() time = 0.05337312 on level 0 Castro::construct_old_gravity() time = 0.004345101 on level 0 Castro::do_old_sources() time = 0.007427109 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.131310621e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.052900265 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.304959237 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.013356754 Gravity::solve_for_phi() time = 0.096953801 on level 0 Castro::construct_new_gravity() time = 0.098993314 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.590666703e+27 zmom added: 0 eint added: 0 ener added: 4.583540443e+43 Castro::do_new_sources() time = 0.053015758 on level 0 Castro::construct_old_gravity() time = 0.003006333 on level 0 Castro::do_old_sources() time = 0.007449474 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.131310621e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.052219382 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.308590329 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.013791053 Gravity::solve_for_phi() time = 0.097751067 on level 0 Castro::construct_new_gravity() time = 0.099802417 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.59568894e+27 zmom added: 0 eint added: 0 ener added: 4.58255499e+43 Castro::do_new_sources() time = 0.052843213 on level 0 Castro::construct_old_gravity() time = 0.002718691 on level 0 Castro::do_old_sources() time = 0.007363559 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.400909405e+26 zmom added: 0 eint added: 0 ener added: 8.721813599e+46 Castro::do_old_sources() time = 0.052308611 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.316744602 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.01381975 Gravity::solve_for_phi() time = 0.100017047 on level 0 Castro::construct_new_gravity() time = 0.102115763 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.9313312e+24 zmom added: 0 eint added: 0 ener added: 6.012622799e+43 Castro::do_new_sources() time = 0.053005629 on level 0 Castro::construct_old_gravity() time = 0.009577293 on level 0 Castro::do_old_sources() time = 0.013368594 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.400909405e+26 zmom added: 0 eint added: 0 ener added: 8.721813599e+46 Castro::do_old_sources() time = 0.052437851 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.330036271 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.013588446 Gravity::solve_for_phi() time = 0.098834046 on level 0 Castro::construct_new_gravity() time = 0.100905507 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.547561128e+24 zmom added: 0 eint added: 0 ener added: 5.281428219e+43 Castro::do_new_sources() time = 0.054550401 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.002631431 on level 0 Castro::do_old_sources() time = 0.007473392 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.526335459e+26 zmom added: 0 eint added: 0 ener added: 8.780735363e+46 Castro::do_old_sources() time = 0.054102925 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.318186357 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.013819121 Gravity::solve_for_phi() time = 0.097710565 on level 0 Castro::construct_new_gravity() time = 0.099798578 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: 4.828619729e+23 zmom added: 0 eint added: 0 ener added: 5.994705329e+43 Castro::do_new_sources() time = 0.052488292 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.004268421 on level 0 Castro::do_old_sources() time = 0.007734745 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.526335459e+26 zmom added: 0 eint added: 0 ener added: 8.780735363e+46 Castro::do_old_sources() time = 0.052854421 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.321118442 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.013567421 Gravity::solve_for_phi() time = 0.097785057 on level 0 Castro::construct_new_gravity() time = 0.099832092 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.918422486e+24 zmom added: 0 eint added: 0 ener added: 5.260150257e+43 Castro::do_new_sources() time = 0.052854213 on level 0 Castro::construct_old_gravity() time = 0.004094451 on level 0 Castro::do_old_sources() time = 0.00859761 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.588292907e+26 zmom added: 0 eint added: 0 ener added: 8.84101861e+46 Castro::do_old_sources() time = 0.052341229 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.309436541 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.013667418 Gravity::solve_for_phi() time = 0.09642704 on level 0 Castro::construct_new_gravity() time = 0.098467169 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.376492035e+23 zmom added: 0 eint added: 0 ener added: 4.637706889e+43 Castro::do_new_sources() time = 0.052231686 on level 0 Castro::construct_old_gravity() time = 0.002971563 on level 0 Castro::do_old_sources() time = 0.007460365 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.588292907e+26 zmom added: 0 eint added: 0 ener added: 8.84101861e+46 Castro::do_old_sources() time = 0.052482246 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.320826339 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.013735992 Gravity::solve_for_phi() time = 0.095345939 on level 0 Castro::construct_new_gravity() time = 0.09850467 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.04399717e+25 zmom added: 0 eint added: 0 ener added: 3.901837814e+43 Castro::do_new_sources() time = 0.052282285 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.002786148 on level 0 Castro::do_old_sources() time = 0.008764577 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.382775517e+26 zmom added: 0 eint added: 0 ener added: 8.902573541e+46 Castro::do_old_sources() time = 0.052239911 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.315716529 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.014038211 Gravity::solve_for_phi() time = 0.101904372 on level 0 Castro::construct_new_gravity() time = 0.104145727 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.27503895e+24 zmom added: 0 eint added: 0 ener added: 2.081560288e+43 Castro::do_new_sources() time = 0.052875176 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.002660328 on level 0 Castro::do_old_sources() time = 0.007340505 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.382775517e+26 zmom added: 0 eint added: 0 ener added: 8.902573541e+46 Castro::do_old_sources() time = 0.052183089 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.318965363 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.013707235 Gravity::solve_for_phi() time = 0.097226546 on level 0 Castro::construct_new_gravity() time = 0.099256881 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.760223209e+24 zmom added: 0 eint added: 0 ener added: 1.337821578e+43 Castro::do_new_sources() time = 0.052373832 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.006263118 on level 0 Castro::do_old_sources() time = 0.010746425 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.437932758e+26 zmom added: 0 eint added: 0 ener added: 8.964911027e+46 Castro::do_old_sources() time = 0.052336404 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.311654153 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.014014136 Gravity::solve_for_phi() time = 0.095407636 on level 0 Castro::construct_new_gravity() time = 0.097435165 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.725986052e+24 zmom added: 0 eint added: 0 ener added: -1.419860559e+43 Castro::do_new_sources() time = 0.053570013 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.002564564 on level 0 Castro::do_old_sources() time = 0.007285003 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.437932758e+26 zmom added: 0 eint added: 0 ener added: 8.964911027e+46 Castro::do_old_sources() time = 0.052403912 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.316795232 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.013861894 Gravity::solve_for_phi() time = 0.095353345 on level 0 Castro::construct_new_gravity() time = 0.097394771 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.640866764e+24 zmom added: 0 eint added: 0 ener added: -2.166675325e+43 Castro::do_new_sources() time = 0.052335344 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.002644075 on level 0 Castro::do_old_sources() time = 0.007604899 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.555803025e+26 zmom added: 0 eint added: 0 ener added: 9.027797292e+46 Castro::do_old_sources() time = 0.053666711 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.312708486 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.013875906 Gravity::solve_for_phi() time = 0.097596939 on level 0 Castro::construct_new_gravity() time = 0.099733215 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.400103868e+23 zmom added: 0 eint added: 0 ener added: -4.056999371e+43 Castro::do_new_sources() time = 0.055983348 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.002813689 on level 0 Castro::do_old_sources() time = 0.007440526 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.555803025e+26 zmom added: 0 eint added: 0 ener added: 9.027797292e+46 Castro::do_old_sources() time = 0.052159726 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.329479777 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.014067555 Gravity::solve_for_phi() time = 0.096112456 on level 0 Castro::construct_new_gravity() time = 0.098140147 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.386922874e+24 zmom added: 0 eint added: 0 ener added: -4.811917635e+43 Castro::do_new_sources() time = 0.052264817 on level 0 Castro::construct_old_gravity() time = 0.002816018 on level 0 Castro::do_old_sources() time = 0.007397564 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.377486467e+26 zmom added: 0 eint added: 0 ener added: 9.091359292e+46 Castro::do_old_sources() time = 0.052617355 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.311226038 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.013941192 Gravity::solve_for_phi() time = 0.09564237 on level 0 Castro::construct_new_gravity() time = 0.097675919 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.208848025e+24 zmom added: 0 eint added: 0 ener added: -1.081736549e+43 Castro::do_new_sources() time = 0.053341325 on level 0 Castro::construct_old_gravity() time = 0.002643184 on level 0 Castro::do_old_sources() time = 0.007404363 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.377486467e+26 zmom added: 0 eint added: 0 ener added: 9.091359292e+46 Castro::do_old_sources() time = 0.052254842 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.315696376 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002801249544 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.013935793 Gravity::solve_for_phi() time = 0.095272894 on level 0 Castro::construct_new_gravity() time = 0.097315697 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.854709436e+24 zmom added: 0 eint added: 0 ener added: -1.838318421e+43 Castro::do_new_sources() time = 0.052329969 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.118014598e+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.467862322e+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 = 4.854673355e-07 TIME= 0.25 CENTER OF MASS Y-VEL = 4.384376873e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 57600050.72 TIME= 0.25 MAXIMUM DENSITY = 4458760.788 TIME= 0.25 MAXIMUM T_S / T_E = 0 STEP = 1 TIME = 0.25 DT = 0.25 [STEP 1] Coarse TimeStep time: 11.75119021 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.047306208 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004166191 seconds Ending run at 07:50:54 UTC on 2026-10-22. Run time = 14.06591966 Run time without initialization = 11.80358306 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.0575 ... 14.0594 ... 14.0660 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.1472 5.1839 5.3075 37.73% Amr::InitAmr() 1 1.315 1.317 1.324 9.41% Castro::computeTemp() 80 0.7327 0.8177 0.9518 6.77% amrex::Copy() 766 0.7939 0.799 0.8206 5.83% FabArray::ParallelCopy_finish() 1469 0.2579 0.5949 0.79 5.62% MultiFab::AddProduct() 3598 0.5377 0.5409 0.5667 4.03% MLPoisson::Fsmooth() 6144 0.3391 0.365 0.4388 3.12% Castro::initMFs() 1 0.3944 0.4001 0.4204 2.99% FabArray::sum() 3598 0.4022 0.4089 0.4176 2.97% FabArray::setVal() 5232 0.3575 0.3608 0.375 2.67% Castro::reset_internal_energy(Fab) 160 0.3584 0.3615 0.3675 2.61% MultiFab::Subtract() 159 0.3337 0.3359 0.3512 2.50% Castro::do_advance_ctu() 20 0.2622 0.2654 0.2742 1.95% check_for_negative_density() 20 0.08649 0.2217 0.2722 1.94% FabArray::ParallelCopy_nowait() 1469 0.2223 0.2414 0.2624 1.87% Gravity::fill_multipole_BCs() 19 0.1851 0.2329 0.2498 1.78% FillBoundary_finish() 6380 0.1342 0.154 0.19 1.35% Castro::initData() 1 0.003256 0.02835 0.1669 1.19% FillBoundary_nowait() 6380 0.09784 0.1278 0.1599 1.14% TagBoxArray::collate() 2 0.0003061 0.03275 0.1486 1.06% MLCellLinOp::applyBC() 7148 0.08059 0.09569 0.1444 1.03% Castro::normalize_species() 79 0.1177 0.1184 0.1201 0.85% MLPoisson::Fapply() 986 0.08759 0.09206 0.1071 0.76% Castro::construct_old_gravity_source() 40 0.09341 0.0943 0.09847 0.70% MultiFab::contains_nan() 37 0.09071 0.09145 0.0934 0.66% FabArray::Saxpy() 111 0.07953 0.08022 0.08365 0.59% Castro::construct_new_gravity_source() 17 0.06672 0.06755 0.0692 0.49% FabArray::Xpay() 915 0.04992 0.05189 0.05708 0.41% MLMG::addInterpCorrection() 768 0.03324 0.03576 0.04345 0.31% AmrLevel::derive() 18 0.03921 0.03967 0.04069 0.29% amrex::average_down 768 0.02618 0.02838 0.03552 0.25% StateData::FillBoundary(geom) 193 0 0.01847 0.03546 0.25% FabArray::FillBoundary() 6380 0.02562 0.02622 0.03087 0.22% FabArray::norminf() 244 0.02296 0.02371 0.02963 0.21% Castro::reset_internal_energy(MultiFab) 80 0.02539 0.02572 0.02665 0.19% FabArrayBase::getFB() 6380 0.01976 0.02024 0.02411 0.17% VisMF::Write(FabArray) 4 0.01264 0.01369 0.02362 0.17% MLCellLinOp::defineAuxData() 19 0.01609 0.01707 0.02088 0.15% Castro::enforce_min_density() 79 0.01984 0.02019 0.02055 0.15% Castro::construct_old_source() 240 0.01949 0.0196 0.01995 0.14% MLCGSolver::bicgstab 8 0 0.001226 0.01961 0.14% FabArray::setDomainBndry() 187 0.0009061 0.008851 0.01707 0.12% Castro::evaluate_source_change() 196 0.01607 0.01633 0.01705 0.12% FillPatchIterator::Initialize 187 0.01618 0.0166 0.01687 0.12% Amr::writePlotFile() 2 0.006941 0.0155 0.01666 0.12% MLCellLinOp::smooth() 1544 0.01458 0.01492 0.01581 0.11% Castro::construct_new_source() 102 0.01512 0.01525 0.01545 0.11% amrex::Add() 128 0.01367 0.01377 0.01416 0.10% MLMG::ResNormInf() 146 0.004834 0.009026 0.01245 0.09% amrex::Dot() 142 0 0.0007384 0.01181 0.08% FabArrayBase::getCPC() 2045 0.01125 0.01139 0.01169 0.08% MLMG::prepareForSolve() 18 0.01024 0.01056 0.01164 0.08% Castro::wd_update() 2 0.00414 0.008312 0.01069 0.08% Gravity::actual_solve_with_mlmg() 19 0.007485 0.007941 0.009633 0.07% MLCGSolver::ParallelAllReduce 189 0 0.000602 0.009632 0.07% MLMG::solve() 18 0.001418 0.00405 0.009544 0.07% MLCellLinOp::apply() 986 0.003767 0.004117 0.008848 0.06% FabArray::ParallelCopy() 1469 0.006641 0.006775 0.006988 0.05% Gravity::get_old_grav_vector() 20 0.00626 0.006304 0.006402 0.05% Gravity::get_new_grav_vector() 18 0.005811 0.005923 0.006045 0.04% Castro::initialize_do_advance() 20 0.00467 0.00489 0.005058 0.04% Castro::sum_integrated_quantities() 2 0.002651 0.003378 0.003765 0.03% MLCellLinOp::correctionResidual() 768 0.003622 0.003671 0.003762 0.03% AmrLevel::FillPatch() 187 0.003365 0.003471 0.003671 0.03% MLLinOp::defineCoarseMGLevels() 19 0.002642 0.002778 0.003103 0.02% MLMG:computeResOfCorrection() 768 0.002792 0.002833 0.002906 0.02% MLMG::mgVcycle() 128 0.002579 0.002642 0.002782 0.02% StateDataPhysBCFunct::() 187 0.001195 0.001392 0.001543 0.01% Castro::save_data_for_retry() 3 0.001328 0.001373 0.001541 0.01% MLMG::actualBottomSolve() 128 0.0003225 0.0004023 0.001526 0.01% FillPatchIterator::FillFromLevel0() 187 0.001119 0.001146 0.001176 0.01% Castro::subcycle_advance_ctu() 1 0.0002585 0.0003364 0.001049 0.01% FillPatchSingleLevel 187 0.000922 0.0009421 0.0009866 0.01% main() 1 0.0007978 0.0008506 0.0008943 0.01% Castro::Castro() 1 0.0004815 0.0005683 0.0008882 0.01% MLMG::mgVcycle_down::0 128 0.0008424 0.0008514 0.0008687 0.01% MLMG::mgVcycle_down::2 128 0.0007836 0.0008133 0.0008613 0.01% MLMG::mgVcycle_down::5 128 0.0007948 0.0008057 0.0008399 0.01% MLMG::mgVcycle_down::3 128 0.0007465 0.0007786 0.0008216 0.01% MLMG::mgVcycle_down::1 128 0.0007861 0.0008015 0.0008191 0.01% MLMG::mgVcycle_down::4 128 0.0007398 0.0007568 0.0008143 0.01% Gravity::solve_for_phi() 17 0.0007294 0.0007547 0.0008008 0.01% MLCellLinOp::solutionResidual() 147 0.0007154 0.0007234 0.0007436 0.01% Castro::clean_state() 79 0.0006303 0.0006497 0.0006822 0.00% MLMG::mgVcycle_up::1 128 0.0006134 0.0006243 0.0006352 0.00% MLMG::mgVcycle_up::0 128 0.000611 0.0006197 0.0006327 0.00% AmrMesh::MakeNewGrids() 2 0.0005736 0.0005971 0.0006292 0.00% MLMG::mgVcycle_up::2 128 0.0005801 0.0006021 0.0006269 0.00% MLMG::oneIter() 128 0.0005915 0.0006029 0.0006238 0.00% MLMG::mgVcycle_up::3 128 0.0005627 0.0005756 0.0006196 0.00% MLMG::mgVcycle_up::4 128 0.0005607 0.0005712 0.0005979 0.00% Castro::do_old_sources() 40 0.0005424 0.0005566 0.0005785 0.00% MLMG::mgVcycle_up::5 128 0.0005241 0.000537 0.0005484 0.00% MLLinOp::defineGrids() 19 0.0003344 0.0003855 0.0005028 0.00% MLMG::mgVcycle_bottom 128 0.0004343 0.0004452 0.0004633 0.00% MLMG::computeResidual() 128 0.0004379 0.0004468 0.0004615 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002457 0.0002683 0.0003515 0.00% Castro::construct_new_gravity() 17 0.0002079 0.0002315 0.00033 0.00% Castro::do_new_sources() 17 0.0002504 0.000258 0.0002663 0.00% Castro::problem_post_timestep() 2 1.47e-05 3.07e-05 0.0002535 0.00% Castro::check_for_nan() 37 0.0001684 0.0001793 0.0001987 0.00% MLPoisson::define() 19 0.0001656 0.000176 0.0001837 0.00% Castro::apply_source_to_state() 37 0.0001595 0.0001668 0.0001767 0.00% Castro::construct_old_gravity() 20 0.000151 0.0001622 0.0001741 0.00% Amr::coarseTimeStep() 1 0.0001321 0.0001574 0.0001629 0.00% Castro::retry_advance_ctu() 11 8.349e-05 9.842e-05 0.0001564 0.00% Amr::defBaseLevel() 1 0.0001038 0.0001308 0.0001543 0.00% Castro::initialize_advance() 1 0.00012 0.0001268 0.0001518 0.00% TagBoxArray::mapPRD 2 0.0001227 0.0001262 0.0001327 0.00% Gravity::test_composite_phi() 1 0.0001093 0.0001146 0.0001285 0.00% MLLinOp::define() 19 0.0001004 0.0001063 0.0001136 0.00% Amr::grid_places() 2 6.402e-05 7.865e-05 0.0001035 0.00% Castro::create_source_corrector() 17 7.713e-05 8.243e-05 9.321e-05 0.00% Castro::expand_state() 20 8.117e-05 8.397e-05 8.889e-05 0.00% Amr::FinalizeInit() 1 1.52e-05 2.455e-05 8.866e-05 0.00% MLMG::computeMLResidual() 18 7.6e-05 7.94e-05 8.519e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.085e-05 6.424e-05 8.018e-05 0.00% Castro::derive() 18 6.949e-05 7.279e-05 7.723e-05 0.00% Amr::regrid() 1 1.694e-05 2.349e-05 5.808e-05 0.00% Gravity::actual_multilevel_solve() 1 3.338e-05 3.538e-05 3.699e-05 0.00% Amr::timeStep() 1 1.68e-05 2.006e-05 2.616e-05 0.00% Castro::advance() 1 1.185e-05 1.347e-05 2.201e-05 0.00% Amr::initialInit() 1 4.277e-06 5.596e-06 1.731e-05 0.00% Castro::post_init() 1 1.46e-05 1.55e-05 1.719e-05 0.00% Amr::InitializeInit() 1 9.249e-06 1.08e-05 1.638e-05 0.00% Castro::post_regrid() 2 1.389e-05 1.457e-05 1.544e-05 0.00% Amr::bldFineLevels() 1 8.921e-06 1.074e-05 1.314e-05 0.00% Amr::init() 1 8.128e-06 9.106e-06 1.086e-05 0.00% Castro::post_timestep() 1 7.558e-06 7.966e-06 9.297e-06 0.00% Other 2199 0.04999 0.06383 0.09338 0.66% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.06 14.06 14.07 100.00% Amr::coarseTimeStep() 1 11.8 11.8 11.8 83.91% Amr::timeStep() 1 11.75 11.75 11.75 83.54% Castro::advance() 1 11.69 11.7 11.7 83.16% Castro::subcycle_advance_ctu() 1 11.65 11.65 11.65 82.85% Castro::do_advance_ctu() 20 11.61 11.62 11.62 82.59% Castro::construct_ctu_hydro_source() 20 6.087 6.107 6.145 43.69% Castro::clean_state() 79 3.142 3.236 3.365 23.92% Castro::computeTemp() 80 2.162 2.253 2.384 16.95% Gravity::solve_phi_with_mlmg() 19 1.709 1.758 1.777 12.63% Castro::construct_new_gravity() 17 1.617 1.666 1.686 11.98% Gravity::solve_for_phi() 17 1.581 1.63 1.649 11.73% Gravity::actual_solve_with_mlmg() 19 1.517 1.522 1.525 10.84% Castro::reset_internal_energy(MultiFab) 80 1.427 1.435 1.479 10.51% MLMG::solve() 18 1.457 1.463 1.465 10.41% Amr::InitAmr() 1 1.315 1.317 1.324 9.41% MLMG::oneIter() 128 1.307 1.312 1.314 9.34% MLMG::mgVcycle() 128 1.293 1.297 1.299 9.24% Castro::initialize_do_advance() 20 1.195 1.23 1.258 8.94% Castro::do_old_sources() 40 1.162 1.187 1.216 8.64% Castro::evaluate_source_change() 196 1.153 1.164 1.209 8.59% FabArray::ParallelCopy() 1469 0.5269 0.8656 1.06 7.54% Amr::init() 1 0.9383 0.9384 0.9384 6.67% Castro::enforce_min_density() 79 0.883 0.8896 0.9196 6.54% Castro::do_new_sources() 17 0.8294 0.8629 0.8952 6.36% Amr::initialInit() 1 0.8866 0.8866 0.8867 6.30% MLCellLinOp::smooth() 1544 0.6722 0.7263 0.8611 6.12% amrex::Copy() 766 0.7939 0.799 0.8206 5.83% FabArray::ParallelCopy_finish() 1469 0.2579 0.5949 0.79 5.62% AmrLevel::FillPatch() 187 0.5851 0.6545 0.7132 5.07% Amr::InitializeInit() 1 0.4751 0.5046 0.6516 4.63% Amr::defBaseLevel() 1 0.4751 0.5046 0.6516 4.63% MultiFab::AddProduct() 3598 0.5377 0.5409 0.5667 4.03% FillPatchIterator::Initialize 187 0.3868 0.4574 0.5173 3.68% MLCellLinOp::applyBC() 7148 0.3869 0.4247 0.5169 3.67% MLMG::mgVcycle_up::1 128 0.1134 0.3064 0.499 3.55% FillPatchIterator::FillFromLevel0() 187 0.3589 0.4319 0.4896 3.48% FillPatchSingleLevel 187 0.3578 0.4308 0.4885 3.47% MLPoisson::Fsmooth() 6144 0.3391 0.365 0.4388 3.12% MLMG::addInterpCorrection() 768 0.07015 0.3569 0.434 3.09% Castro::Castro() 1 0.4001 0.406 0.426 3.03% Castro::initMFs() 1 0.3945 0.4003 0.4206 2.99% FabArray::sum() 3598 0.4022 0.4089 0.4176 2.97% Amr::FinalizeInit() 1 0.235 0.382 0.4114 2.93% FabArray::setVal() 5232 0.3575 0.3608 0.375 2.67% FabArray::FillBoundary() 6380 0.2959 0.329 0.3725 2.65% Castro::reset_internal_energy(Fab) 160 0.3584 0.3615 0.3675 2.61% MultiFab::Subtract() 159 0.3337 0.3359 0.3512 2.50% MLMG::mgVcycle_up::2 128 0.0209 0.1455 0.3377 2.40% FabArray::ParallelCopy_nowait() 1469 0.2443 0.2639 0.2858 2.03% check_for_negative_density() 20 0.08649 0.2217 0.2722 1.94% MLMG::mgVcycle_down::0 128 0.2667 0.2685 0.2721 1.93% Gravity::fill_multipole_BCs() 19 0.1871 0.2349 0.2518 1.79% Castro::expand_state() 20 0.1862 0.2149 0.2418 1.72% Castro::initData() 1 0.07062 0.09687 0.2355 1.67% MLMG::mgVcycle_up::3 128 0.02086 0.03653 0.2325 1.65% MLCellLinOp::apply() 986 0.1594 0.1723 0.2245 1.60% Amr::bldFineLevels() 1 0.04237 0.1892 0.2181 1.55% MLMG::mgVcycle_up::0 128 0.2089 0.2123 0.2152 1.53% FillBoundary_finish() 6380 0.1342 0.154 0.19 1.35% FillBoundary_nowait() 6380 0.1182 0.1488 0.1814 1.29% Amr::grid_places() 2 0.003778 0.1508 0.1802 1.28% AmrMesh::MakeNewGrids() 2 0.0037 0.1507 0.1802 1.28% TagBoxArray::mapPRD 2 0.0008645 0.1165 0.1783 1.27% Castro::construct_old_source() 240 0.156 0.1574 0.1631 1.16% MLMG:computeResOfCorrection() 768 0.1255 0.1352 0.159 1.13% MLCellLinOp::correctionResidual() 768 0.1227 0.1324 0.1561 1.11% Castro::post_init() 1 0.1515 0.1548 0.156 1.11% TagBoxArray::collate() 2 0.0003061 0.03275 0.1486 1.06% MLMG::mgVcycle_down::1 128 0.1348 0.1383 0.1434 1.02% Castro::construct_new_source() 102 0.1186 0.1198 0.1219 0.87% Castro::normalize_species() 79 0.1177 0.1184 0.1201 0.85% Gravity::multilevel_solve_for_new_phi() 1 0.1169 0.1169 0.117 0.83% Gravity::actual_multilevel_solve() 1 0.1168 0.1168 0.1169 0.83% MLMG::mgVcycle_bottom 128 0.0007581 0.007618 0.1103 0.78% MLMG::actualBottomSolve() 128 0.0003225 0.007173 0.1099 0.78% MLCGSolver::bicgstab 8 0 0.006704 0.1073 0.76% MLPoisson::Fapply() 986 0.08759 0.09206 0.1071 0.76% MLCellLinOp::solutionResidual() 147 0.08824 0.09361 0.09875 0.70% Castro::construct_old_gravity_source() 40 0.09341 0.0943 0.09847 0.70% Amr::writePlotFile() 2 0.09476 0.09486 0.09487 0.67% Castro::check_for_nan() 37 0.0909 0.09163 0.09359 0.67% MultiFab::contains_nan() 37 0.09071 0.09145 0.0934 0.66% MLMG::mgVcycle_down::2 128 0.02458 0.0545 0.08927 0.63% FabArray::Saxpy() 111 0.07953 0.08022 0.08365 0.59% Castro::create_source_corrector() 17 0.04997 0.06102 0.08138 0.58% Castro::post_regrid() 2 0.07511 0.07629 0.07967 0.57% MLMG::computeResidual() 128 0.07517 0.0766 0.07931 0.56% MLMG::mgVcycle_down::3 128 0.02391 0.03039 0.07436 0.53% Castro::derive() 18 0.06785 0.06836 0.06926 0.49% Castro::construct_new_gravity_source() 17 0.06672 0.06755 0.0692 0.49% AmrLevel::derive() 18 0.06778 0.06829 0.06918 0.49% Castro::construct_old_gravity() 20 0.04567 0.06134 0.06882 0.49% Gravity::get_old_grav_vector() 20 0.04551 0.06118 0.06866 0.49% amrex::average_down 768 0.03941 0.04677 0.06506 0.46% Castro::post_timestep() 1 0.05333 0.05515 0.05711 0.41% FabArray::Xpay() 915 0.04992 0.05189 0.05708 0.41% MLMG::mgVcycle_down::4 128 0.02393 0.02564 0.04773 0.34% MLPoisson::define() 19 0.03588 0.04102 0.04663 0.33% Castro::initialize_advance() 1 0.04348 0.04393 0.04622 0.33% Amr::regrid() 1 0.04009 0.04065 0.04232 0.30% MLMG::mgVcycle_up::4 128 0.0209 0.0221 0.03837 0.27% Castro::retry_advance_ctu() 11 0.03421 0.03497 0.03823 0.27% Gravity::get_new_grav_vector() 18 0.03474 0.03559 0.03739 0.27% StateDataPhysBCFunct::() 187 0.001195 0.01986 0.03691 0.26% MLMG::mgVcycle_down::5 128 0.02405 0.02498 0.03686 0.26% StateData::FillBoundary(geom) 193 0 0.01847 0.03546 0.25% MLMG::ResNormInf() 146 0.02525 0.02954 0.03292 0.23% MLMG::mgVcycle_up::5 128 0.02089 0.02171 0.03146 0.22% MLCellLinOp::defineAuxData() 19 0.02236 0.02418 0.03009 0.21% FabArray::norminf() 244 0.02296 0.02371 0.02963 0.21% Castro::apply_source_to_state() 37 0.02663 0.02684 0.0275 0.20% FabArrayBase::getFB() 6380 0.02018 0.02102 0.02572 0.18% MLMG::prepareForSolve() 18 0.02081 0.02171 0.02454 0.17% VisMF::Write(FabArray) 4 0.01264 0.01369 0.02362 0.17% Castro::sum_integrated_quantities() 2 0.023 0.02315 0.02317 0.16% MLMG::computeMLResidual() 18 0.01047 0.01674 0.02189 0.16% Castro::save_data_for_retry() 3 0.01817 0.01853 0.02118 0.15% FabArrayBase::getCPC() 2045 0.01423 0.01556 0.01919 0.14% FabArray::setDomainBndry() 187 0.0009061 0.008851 0.01707 0.12% MLCGSolver::ParallelAllReduce 189 0 0.001004 0.01607 0.11% Gravity::test_composite_phi() 1 0.01416 0.01424 0.01427 0.10% amrex::Add() 128 0.01367 0.01377 0.01416 0.10% Castro::problem_post_timestep() 2 0.01014 0.01198 0.01402 0.10% Castro::wd_update() 2 0.01013 0.01195 0.014 0.10% MLLinOp::define() 19 0.00505 0.009278 0.01278 0.09% MLLinOp::defineGrids() 19 0.004939 0.009171 0.01267 0.09% MLLinOp::defineCoarseMGLevels() 19 0.004584 0.008786 0.01222 0.09% amrex::Dot() 142 0 0.0007384 0.01181 0.08% Other 2199 0.08489 0.1017 0.1343 0.95% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Async Memory Usage: -------------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max -------------------------------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 6400 387 KiB 390 KiB 400 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 2139 B 2178 B 2197 B 89 KiB 89 KiB 89 KiB MLMG::addInterpCorrection() 5632 1 B 18 B 28 B 1024 B 1024 B 1024 B amrex::average_down 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B MLCGSolver::bicgstab 640 0 B 0 B 4 B 0 B 36 B 576 B -------------------------------------------------------------------------------------------------------------------- Pinned Memory Usage: -------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max -------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 2 B 2 B 3 B 8192 KiB 8192 KiB 8192 KiB -------------------------------------------------------------------------------------------------------------- Cpu Memory Usage: ----------------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max ----------------------------------------------------------------------------------------------------------------------- Castro::save_data_for_retry() 416 1585 KiB 1587 KiB 1587 KiB 2332 KiB 2332 KiB 2332 KiB Castro::construct_ctu_hydro_source() 6400 388 KiB 391 KiB 401 KiB 1178 KiB 1178 KiB 1178 KiB StateData::define() 192 1054 KiB 1055 KiB 1055 KiB 1166 KiB 1166 KiB 1166 KiB Castro::initMFs() 256 1033 KiB 1034 KiB 1034 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 673 KiB 674 KiB 674 KiB 804 KiB 804 KiB 804 KiB Amr::writePlotFile() 64 4352 B 4416 B 4925 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 55 KiB 55 KiB 57 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 12 KiB 13 KiB 15 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 434 KiB 434 KiB 434 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2234 B 4039 B 5156 B 226 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1394 B 1405 B 1428 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 20 KiB 21 KiB 21 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 131 KiB 131 KiB 131 KiB 158 KiB 158 KiB 158 KiB Castro::buildMetrics() 192 129 KiB 129 KiB 129 KiB 143 KiB 143 KiB 143 KiB MLMG::prepareForSolve() 6102 10 KiB 11 KiB 13 KiB 104 KiB 110 KiB 130 KiB Castro::construct_old_source() 1280 491 B 495 B 512 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 413 B 416 B 428 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 2145 B 2183 B 2202 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 18 B 20 B 37 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 34 KiB 34 KiB 34 KiB 76 KiB 76 KiB 76 KiB Gravity::test_composite_phi() 96 53 B 53 B 54 B 53 KiB 53 KiB 53 KiB Gravity::update_max_rhs() 192 3 B 7 B 9 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 112 B 114 B 121 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 148 B 200 B 225 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 32 KiB 32 KiB 32 KiB 38 KiB 38 KiB 38 KiB Gravity::fill_multipole_BCs() 2432 238 B 246 B 250 B 18 KiB 18 KiB 18 KiB MLCellLinOp::defineAuxData() 11704 873 B 1050 B 1959 B 8192 B 9856 B 18 KiB Castro::evaluate_source_change() 6272 1336 B 1349 B 1401 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 135 B 135 B 136 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 1840 B 1896 B 1919 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 0 B 0 B 1 B 16 KiB 16 KiB 16 KiB BndryData::define() 4864 860 B 863 B 865 B 8192 B 8192 B 8192 B FillBoundary_nowait() 91296 31 B 43 B 55 B 3072 B 3840 B 4608 B AmrMesh::MakeNewGrids() 128 0 B 25 B 39 B 3584 B 3584 B 3584 B OwnerMask() 64 0 B 0 B 0 B 3072 B 3072 B 3072 B MLCGSolver::bicgstab 640 0 B 0 B 9 B 0 B 80 B 1280 B MLMG::addInterpCorrection() 5632 1 B 18 B 28 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 0 B 10 B 13 B 1024 B 1024 B 1024 B amrex::average_down 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B MLCellLinOp::applyBC() 114375 0 B 0 B 1 B 256 B 256 B 256 B ----------------------------------------------------------------------------------------------------------------------- AMReX (26.10-97-gf5d545d2f0b0) finalized