Initializing AMReX (25.12-15-g4fc753db6b0b)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (25.12-15-g4fc753db6b0b) initialized Starting run at 22:17:00 UTC on 2025-12-13. Successfully read inputs file ... Castro git describe: 25.12-3-gde831c614 AMReX git describe: 25.12-15-g4fc753db6 Microphysics git describe: 25.12-6-g18b6ec58 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.166988678 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001215002 new finest: 0 Gravity::fill_multipole_BCs() time = 0.01191993 Gravity::multilevel_solve_for_new_phi() time = 0.113505662 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010593045 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.05002871 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004181246 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 = -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 = -80.16477044 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 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 = -552.769409 at index 17, 181, 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 10Invalid 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 , 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, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = 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, -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 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 = -552.769409 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 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 = -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.002870353 on level 0 Castro::do_old_sources() time = 0.007910626 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.061959817 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.005120539 on level 0 Castro::do_old_sources() time = 0.00808129 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.061970845 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.006033795 on level 0 Castro::do_old_sources() time = 0.010111904 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.058363845 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.004843279 on level 0 Castro::do_old_sources() time = 0.007834634 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.057346973 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.331469607 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.014649 Gravity::solve_for_phi() time = 0.090655275 on level 0 Castro::construct_new_gravity() time = 0.092840135 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.057719042 on level 0 Castro::construct_old_gravity() time = 0.00274724 on level 0 Castro::do_old_sources() time = 0.007717026 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.057070611 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.330802786 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.014203779 Gravity::solve_for_phi() time = 0.09160245 on level 0 Castro::construct_new_gravity() time = 0.094371915 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.05743447 on level 0 Castro::construct_old_gravity() time = 0.003041539 on level 0 Castro::do_old_sources() time = 0.007914597 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.056984634 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.336899687 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.013816606 Gravity::solve_for_phi() time = 0.091416456 on level 0 Castro::construct_new_gravity() time = 0.093619085 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.05731617 on level 0 Castro::construct_old_gravity() time = 0.00286508 on level 0 Castro::do_old_sources() time = 0.00775082 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.057184288 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.341098584 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.015661247 Gravity::solve_for_phi() time = 0.09297336 on level 0 Castro::construct_new_gravity() time = 0.095701833 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.057892358 on level 0 Castro::construct_old_gravity() time = 0.002934702 on level 0 Castro::do_old_sources() time = 0.007774663 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.057327704 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.346085312 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.01409798 Gravity::solve_for_phi() time = 0.090782177 on level 0 Castro::construct_new_gravity() time = 0.092924205 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.05745768 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.002787961 on level 0 Castro::do_old_sources() time = 0.007856899 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.05726649 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.342583949 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.015880801 Gravity::solve_for_phi() time = 0.092927076 on level 0 Castro::construct_new_gravity() time = 0.09511891 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.058039706 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.002846524 on level 0 Castro::do_old_sources() time = 0.007759688 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.057528517 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.345947321 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.01443781 Gravity::solve_for_phi() time = 0.0912086 on level 0 Castro::construct_new_gravity() time = 0.093360848 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.057497215 on level 0 Castro::construct_old_gravity() time = 0.003077977 on level 0 Castro::do_old_sources() time = 0.007949339 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.057217068 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.341394114 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.015605622 Gravity::solve_for_phi() time = 0.092724292 on level 0 Castro::construct_new_gravity() time = 0.094855341 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.057820837 on level 0 Castro::construct_old_gravity() time = 0.002996255 on level 0 Castro::do_old_sources() time = 0.007874194 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.057478795 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.34649818 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.015288598 Gravity::solve_for_phi() time = 0.091027932 on level 0 Castro::construct_new_gravity() time = 0.093152848 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.057762918 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.003135677 on level 0 Castro::do_old_sources() time = 0.0078244 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.05733496 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.342340935 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.015836826 Gravity::solve_for_phi() time = 0.092262117 on level 0 Castro::construct_new_gravity() time = 0.094399009 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.057522772 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.003540368 on level 0 Castro::do_old_sources() time = 0.008460097 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.057278732 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.346961398 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.015270142 Gravity::solve_for_phi() time = 0.090968757 on level 0 Castro::construct_new_gravity() time = 0.093115963 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.057659959 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.003251639 on level 0 Castro::do_old_sources() time = 0.007825412 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.057068867 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.341362692 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.0157581 Gravity::solve_for_phi() time = 0.093201103 on level 0 Castro::construct_new_gravity() time = 0.095402337 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.057671184 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.004194272 on level 0 Castro::do_old_sources() time = 0.007806077 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.057623213 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.351404669 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.015557694 Gravity::solve_for_phi() time = 0.091805457 on level 0 Castro::construct_new_gravity() time = 0.094005474 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.058195277 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.00709034 on level 0 Castro::do_old_sources() time = 0.00812741 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.057995781 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.34305133 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.015862753 Gravity::solve_for_phi() time = 0.091541262 on level 0 Castro::construct_new_gravity() time = 0.093710417 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.058044886 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.002783453 on level 0 Castro::do_old_sources() time = 0.007669052 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.057298726 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.349980519 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.016175375 Gravity::solve_for_phi() time = 0.092544408 on level 0 Castro::construct_new_gravity() time = 0.094674034 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.057986491 on level 0 Castro::construct_old_gravity() time = 0.002750892 on level 0 Castro::do_old_sources() time = 0.007689123 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.057015424 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.342128621 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.016379293 Gravity::solve_for_phi() time = 0.092229401 on level 0 Castro::construct_new_gravity() time = 0.094362794 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.058388332 on level 0 Castro::construct_old_gravity() time = 0.003299361 on level 0 Castro::do_old_sources() time = 0.007736285 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.057366423 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.347381647 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002610418616 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.015187361 Gravity::solve_for_phi() time = 0.091485 on level 0 Castro::construct_new_gravity() time = 0.093634122 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.057565075 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.581215356e+29 ener added: 1.581215356e+29 TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -4.341809003e+38 TIME= 0.25 YMOM = 1.137357411e+28 TIME= 0.25 ZMOM = 0 TIME= 0.25 ANG MOM X = 0 TIME= 0.25 ANG MOM Y = 0 TIME= 0.25 ANG MOM Z = -1.469316165e+38 TIME= 0.25 RHO*e = 1.325894512e+50 TIME= 0.25 RHO*K = 3.170941987e+49 TIME= 0.25 RHO*E = 1.642988711e+50 TIME= 0.25 RHO*PHI = -5.449011378e+50 TIME= 0.25 TOTAL ENERGY= -1.081516979e+50 TIME= 0.25 CENTER OF MASS X-LOC = 316424388.8 TIME= 0.25 CENTER OF MASS X-VEL = -170267.2489 TIME= 0.25 CENTER OF MASS Y-LOC = 2.427336677e-06 TIME= 0.25 CENTER OF MASS Y-VEL = 4.460231145e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 57600050.72 TIME= 0.25 MAXIMUM DENSITY = 4458760.788 TIME= 0.25 MAXIMUM T_S / T_E = 0 STEP = 1 TIME = 0.25 DT = 0.25 [STEP 1] Coarse TimeStep time: 12.60988702 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.049321575 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004237512 seconds Ending run at 22:17:15 UTC on 2025-12-13. Run time = 14.96077097 Run time without initialization = 12.66444643 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.9538 ... 14.9558 ... 14.9608 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.6946 5.7493 5.8386 39.03% Amr::InitAmr() 1 1.389 1.391 1.396 9.33% Castro::computeTemp() 80 0.8189 0.9271 1.08 7.22% amrex::Copy() 766 0.8645 0.8682 0.8781 5.87% FabArray::ParallelCopy_finish() 1469 0.243 0.5469 0.7275 4.86% MultiFab::AddProduct() 3598 0.5633 0.5655 0.5706 3.81% Castro::reset_internal_energy(Fab) 160 0.5078 0.5125 0.5222 3.49% MLPoisson::Fsmooth() 6144 0.3511 0.3787 0.4558 3.05% FabArray::setVal() 5232 0.3919 0.3955 0.4005 2.68% Castro::initMFs() 1 0.3529 0.3602 0.3648 2.44% MultiFab::Subtract() 159 0.3503 0.352 0.3593 2.40% Castro::normalize_species() 79 0.3446 0.3477 0.3575 2.39% check_for_negative_density() 20 0.1519 0.228 0.3121 2.09% Gravity::fill_multipole_BCs() 19 0.1986 0.2564 0.2773 1.85% Castro::do_advance_ctu() 20 0.2704 0.2721 0.275 1.84% FabArray::ParallelCopy_nowait() 1469 0.2312 0.2499 0.2718 1.82% FabArray::sum() 3598 0.1892 0.1903 0.1921 1.28% Castro::initData() 1 0.003362 0.02837 0.1661 1.11% FillBoundary_finish() 6380 0.09769 0.1242 0.1539 1.03% TagBoxArray::collate() 2 0.0002331 0.02728 0.1514 1.01% FillBoundary_nowait() 6380 0.08865 0.1161 0.1476 0.99% MLCellLinOp::applyBC() 7148 0.06713 0.08183 0.1279 0.85% MLPoisson::Fapply() 986 0.09536 0.09965 0.1148 0.77% Castro::construct_old_gravity_source() 40 0.09835 0.0988 0.09968 0.67% FabArray::Saxpy() 111 0.08534 0.08573 0.08656 0.58% Castro::construct_new_gravity_source() 17 0.07353 0.07436 0.07545 0.50% FabArray::Xpay() 915 0.05133 0.05333 0.05893 0.39% MultiFab::contains_nan() 37 0.04813 0.04849 0.0492 0.33% MLMG::addInterpCorrection() 768 0.03339 0.03574 0.04311 0.29% AmrLevel::derive() 18 0.0401 0.04045 0.04086 0.27% StateData::FillBoundary(geom) 193 0 0.01964 0.03774 0.25% amrex::average_down 768 0.025 0.02704 0.03376 0.23% VisMF::Write(FabArray) 4 0.01374 0.01503 0.02583 0.17% Castro::reset_internal_energy(MultiFab) 80 0.02423 0.02466 0.02524 0.17% Castro::construct_old_source() 240 0.02028 0.02042 0.02076 0.14% FabArray::FillBoundary() 6380 0.01668 0.0171 0.02019 0.13% Castro::enforce_min_density() 79 0.01906 0.01939 0.02 0.13% MLCellLinOp::defineAuxData() 19 0.01531 0.01619 0.01949 0.13% FabArray::setDomainBndry() 187 0.0006464 0.009356 0.01801 0.12% Amr::writePlotFile() 2 0.007288 0.01615 0.01751 0.12% Castro::construct_new_source() 102 0.01601 0.01616 0.01687 0.11% MLCGSolver::bicgstab 8 0 0.001049 0.01679 0.11% FabArrayBase::getFB() 6380 0.01369 0.01404 0.01674 0.11% FillPatchIterator::Initialize 187 0.01523 0.0154 0.01595 0.11% amrex::Add() 128 0.01465 0.0148 0.01551 0.10% Castro::evaluate_source_change() 196 0.01478 0.01502 0.01533 0.10% FabArray::norminf() 244 0.01103 0.01137 0.01472 0.10% MLCellLinOp::smooth() 1544 0.01084 0.01109 0.01167 0.08% MLMG::prepareForSolve() 18 0.01011 0.01042 0.01152 0.08% Castro::wd_update() 2 0.004736 0.008133 0.01123 0.08% Gravity::actual_solve_with_mlmg() 19 0.00735 0.007787 0.009387 0.06% MLMG::ResNormInf() 146 0.003804 0.006418 0.009005 0.06% FabArrayBase::getCPC() 2045 0.008505 0.00869 0.008997 0.06% MLMG::solve() 18 0.0009772 0.003825 0.00776 0.05% Gravity::get_old_grav_vector() 20 0.006517 0.006647 0.00702 0.05% Gravity::get_new_grav_vector() 18 0.00605 0.006199 0.00662 0.04% Castro::initialize_do_advance() 20 0.004739 0.005096 0.006402 0.04% MLCellLinOp::apply() 986 0.002427 0.002672 0.005767 0.04% FabArray::ParallelCopy() 1469 0.004229 0.004332 0.004552 0.03% AmrLevel::FillPatch() 187 0.003014 0.003138 0.003258 0.02% Castro::sum_integrated_quantities() 2 0.001789 0.002301 0.002461 0.02% MLMG::mgVcycle() 128 0.002154 0.002221 0.002375 0.02% MLCellLinOp::correctionResidual() 768 0.002254 0.002297 0.00236 0.02% MLMG:computeResOfCorrection() 768 0.001791 0.00184 0.001911 0.01% Castro::save_data_for_retry() 3 0.001287 0.001336 0.001464 0.01% MLMG::actualBottomSolve() 128 0.0001794 0.0002516 0.001265 0.01% StateDataPhysBCFunct::() 187 0.0009687 0.001128 0.001222 0.01% Castro::subcycle_advance_ctu() 1 0.0002612 0.0003353 0.001072 0.01% main() 1 0.0006992 0.0007807 0.0009547 0.01% FillPatchIterator::FillFromLevel0() 187 0.0008945 0.0009203 0.0009477 0.01% Gravity::solve_for_phi() 17 0.0006631 0.0006888 0.0007283 0.00% FillPatchSingleLevel 187 0.0006487 0.0006604 0.0006841 0.00% MLMG::mgVcycle_down::0 128 0.0006315 0.0006462 0.0006645 0.00% MLMG::mgVcycle_down::2 128 0.0005834 0.0006145 0.0006494 0.00% MLMG::mgVcycle_down::3 128 0.0005624 0.0005875 0.0006446 0.00% MLMG::mgVcycle_down::1 128 0.0006004 0.0006119 0.0006335 0.00% MLMG::mgVcycle_down::5 128 0.0005829 0.0005928 0.0006227 0.00% MLMG::mgVcycle_down::4 128 0.0005562 0.0005682 0.0006117 0.00% AmrMesh::MakeNewGrids() 2 0.000571 0.0005922 0.0006033 0.00% Castro::Castro() 1 0.0004799 0.000529 0.0005922 0.00% MLCellLinOp::solutionResidual() 147 0.0004991 0.0005102 0.0005311 0.00% Castro::clean_state() 79 0.0004658 0.0004855 0.0005047 0.00% Castro::do_old_sources() 40 0.000457 0.0004664 0.0004925 0.00% MLMG::oneIter() 128 0.0004304 0.00044 0.0004526 0.00% MLMG::mgVcycle_up::0 128 0.0004098 0.0004198 0.0004382 0.00% MLMG::mgVcycle_up::1 128 0.0004007 0.0004085 0.0004212 0.00% MLMG::mgVcycle_up::2 128 0.0003854 0.0004002 0.0004211 0.00% MLMG::mgVcycle_up::4 128 0.0003768 0.0003878 0.0004093 0.00% MLMG::mgVcycle_up::3 128 0.0003737 0.0003798 0.0003923 0.00% MLMG::mgVcycle_up::5 128 0.0003628 0.0003717 0.0003831 0.00% Castro::construct_new_gravity() 17 0.0001679 0.0001863 0.0003143 0.00% MLMG::computeResidual() 128 0.0002924 0.0002968 0.0003014 0.00% MLMG::mgVcycle_bottom 128 0.0002738 0.0002819 0.0002972 0.00% Castro::problem_post_timestep() 2 1.207e-05 3.007e-05 0.0002659 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002126 0.0002247 0.0002522 0.00% Castro::do_new_sources() 17 0.0002106 0.0002191 0.0002305 0.00% Gravity::update_max_rhs() 2 0.0001928 0.000203 0.0002297 0.00% MLPoisson::define() 19 0.0001363 0.0001531 0.0001714 0.00% Amr::defBaseLevel() 1 0.0001269 0.0001384 0.0001694 0.00% Amr::coarseTimeStep() 1 0.0001254 0.0001507 0.0001561 0.00% Castro::retry_advance_ctu() 11 6.673e-05 8.076e-05 0.0001414 0.00% Castro::check_for_nan() 37 0.0001194 0.0001262 0.0001406 0.00% Castro::construct_old_gravity() 20 0.0001147 0.0001232 0.0001404 0.00% Castro::initialize_advance() 1 0.0001132 0.000122 0.0001364 0.00% TagBoxArray::mapPRD 2 0.0001192 0.0001243 0.0001338 0.00% Castro::apply_source_to_state() 37 0.0001133 0.0001186 0.000125 0.00% Gravity::test_composite_phi() 1 0.000106 0.000111 0.0001237 0.00% MLLinOp::define() 19 7.995e-05 8.352e-05 8.817e-05 0.00% Amr::grid_places() 2 5.439e-05 6.209e-05 8.664e-05 0.00% Castro::expand_state() 20 5.656e-05 6.094e-05 6.993e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.804e-05 6.104e-05 6.686e-05 0.00% Amr::FinalizeInit() 1 1.225e-05 2.122e-05 6.39e-05 0.00% Amr::regrid() 1 1.76e-05 2.258e-05 6.16e-05 0.00% Castro::create_source_corrector() 17 5.333e-05 5.521e-05 5.722e-05 0.00% MLMG::computeMLResidual() 18 4.86e-05 5.234e-05 5.687e-05 0.00% Castro::post_init() 1 1.726e-05 2.182e-05 5.191e-05 0.00% Castro::derive() 18 4.69e-05 4.8e-05 5.058e-05 0.00% Gravity::actual_multilevel_solve() 1 3.145e-05 3.247e-05 3.677e-05 0.00% Amr::InitializeInit() 1 8.641e-06 1.303e-05 3.349e-05 0.00% Castro::advance() 1 7.956e-06 1.051e-05 2.922e-05 0.00% Amr::timeStep() 1 1.49e-05 1.649e-05 2.033e-05 0.00% Castro::post_regrid() 2 1.126e-05 1.261e-05 1.625e-05 0.00% Amr::bldFineLevels() 1 7.482e-06 8.19e-06 9.664e-06 0.00% Amr::init() 1 6.951e-06 7.873e-06 9.461e-06 0.00% Amr::initialInit() 1 3.352e-06 4.129e-06 7.45e-06 0.00% Castro::post_timestep() 1 5.887e-06 6.371e-06 6.86e-06 0.00% Other 2547 0.0479 0.06558 0.1059 0.71% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.95 14.96 14.96 100.00% Amr::coarseTimeStep() 1 12.66 12.66 12.66 84.65% Amr::timeStep() 1 12.61 12.61 12.61 84.29% Castro::advance() 1 12.55 12.55 12.55 83.91% Castro::subcycle_advance_ctu() 1 12.5 12.5 12.51 83.59% Castro::do_advance_ctu() 20 12.46 12.46 12.47 83.32% Castro::construct_ctu_hydro_source() 20 6.699 6.727 6.773 45.27% Castro::clean_state() 79 3.495 3.617 3.773 25.22% Castro::computeTemp() 80 2.35 2.467 2.623 17.53% Gravity::solve_phi_with_mlmg() 19 1.601 1.66 1.684 11.25% Castro::construct_new_gravity() 17 1.515 1.574 1.598 10.68% Gravity::solve_for_phi() 17 1.476 1.535 1.559 10.42% Castro::reset_internal_energy(MultiFab) 80 1.531 1.54 1.554 10.39% Gravity::actual_solve_with_mlmg() 19 1.395 1.4 1.404 9.38% Castro::initialize_do_advance() 20 1.327 1.364 1.399 9.35% Amr::InitAmr() 1 1.389 1.391 1.396 9.33% MLMG::solve() 18 1.342 1.347 1.349 9.01% Castro::do_old_sources() 40 1.244 1.283 1.311 8.76% MLMG::oneIter() 128 1.207 1.208 1.211 8.09% MLMG::mgVcycle() 128 1.192 1.193 1.195 7.99% FabArray::ParallelCopy() 1469 0.5174 0.8228 1.002 6.70% Castro::evaluate_source_change() 196 0.9829 0.9874 0.9961 6.66% Castro::do_new_sources() 17 0.8905 0.9391 0.981 6.56% Amr::init() 1 0.9001 0.9002 0.9002 6.02% amrex::Copy() 766 0.8645 0.8682 0.8781 5.87% Amr::initialInit() 1 0.8453 0.8454 0.8455 5.65% Castro::enforce_min_density() 79 0.8253 0.8294 0.8402 5.62% MLCellLinOp::smooth() 1544 0.6244 0.6749 0.8112 5.42% AmrLevel::FillPatch() 187 0.6033 0.6797 0.751 5.02% FabArray::ParallelCopy_finish() 1469 0.243 0.5469 0.7275 4.86% Amr::InitializeInit() 1 0.4379 0.4715 0.6034 4.03% Amr::defBaseLevel() 1 0.4379 0.4715 0.6034 4.03% MultiFab::AddProduct() 3598 0.5633 0.5655 0.5706 3.81% FillPatchIterator::Initialize 187 0.3885 0.4656 0.5378 3.59% Castro::reset_internal_energy(Fab) 160 0.5078 0.5125 0.5222 3.49% FillPatchIterator::FillFromLevel0() 187 0.3549 0.4408 0.505 3.38% FillPatchSingleLevel 187 0.3539 0.4399 0.504 3.37% MLPoisson::Fsmooth() 6144 0.3511 0.3787 0.4558 3.05% MLMG::mgVcycle_up::1 128 0.1051 0.2789 0.4528 3.03% MLCellLinOp::applyBC() 7148 0.3133 0.354 0.4335 2.90% Amr::FinalizeInit() 1 0.242 0.3739 0.4075 2.72% FabArray::setVal() 5232 0.3919 0.3955 0.4005 2.68% MLMG::addInterpCorrection() 768 0.06514 0.3203 0.3923 2.62% Castro::Castro() 1 0.3592 0.3667 0.3712 2.48% Castro::initMFs() 1 0.3531 0.3604 0.365 2.44% MultiFab::Subtract() 159 0.3503 0.352 0.3593 2.40% Castro::normalize_species() 79 0.3446 0.3477 0.3575 2.39% FabArray::FillBoundary() 6380 0.24 0.2722 0.3135 2.10% check_for_negative_density() 20 0.1519 0.228 0.3121 2.09% MLMG::mgVcycle_up::2 128 0.01533 0.1272 0.2978 1.99% FabArray::ParallelCopy_nowait() 1469 0.2522 0.2715 0.2943 1.97% Gravity::fill_multipole_BCs() 19 0.2007 0.2586 0.2794 1.87% MLMG::mgVcycle_down::0 128 0.2726 0.2752 0.277 1.85% Castro::expand_state() 20 0.2015 0.2313 0.2602 1.74% Castro::initData() 1 0.07638 0.103 0.2395 1.60% Amr::bldFineLevels() 1 0.04469 0.1777 0.2139 1.43% MLMG::mgVcycle_up::0 128 0.2075 0.2083 0.2095 1.40% MLCellLinOp::apply() 986 0.1584 0.1692 0.2085 1.39% MLMG::mgVcycle_up::3 128 0.01525 0.02907 0.1961 1.31% FabArray::sum() 3598 0.1892 0.1903 0.1921 1.28% Amr::grid_places() 2 0.002666 0.1345 0.1682 1.12% AmrMesh::MakeNewGrids() 2 0.002579 0.1345 0.1681 1.12% Castro::construct_old_source() 240 0.1653 0.1662 0.1676 1.12% TagBoxArray::mapPRD 2 0.0009461 0.1058 0.1654 1.11% FillBoundary_nowait() 6380 0.1028 0.1309 0.1628 1.09% Castro::post_init() 1 0.1465 0.154 0.1562 1.04% FillBoundary_finish() 6380 0.09769 0.1242 0.1539 1.03% MLMG:computeResOfCorrection() 768 0.1202 0.1293 0.1526 1.02% TagBoxArray::collate() 2 0.0002331 0.02728 0.1514 1.01% MLCellLinOp::correctionResidual() 768 0.1183 0.1275 0.1507 1.01% Castro::construct_new_source() 102 0.1293 0.1304 0.1317 0.88% MLMG::mgVcycle_down::1 128 0.1228 0.1268 0.1304 0.87% MLPoisson::Fapply() 986 0.09536 0.09965 0.1148 0.77% Gravity::multilevel_solve_for_new_phi() 1 0.1134 0.1135 0.1135 0.76% Gravity::actual_multilevel_solve() 1 0.1134 0.1134 0.1134 0.76% MLCellLinOp::solutionResidual() 147 0.08911 0.09599 0.1012 0.68% Amr::writePlotFile() 2 0.09978 0.09988 0.09988 0.67% Castro::construct_old_gravity_source() 40 0.09835 0.0988 0.09968 0.67% Castro::post_regrid() 2 0.08306 0.08529 0.09276 0.62% FabArray::Saxpy() 111 0.08534 0.08573 0.08656 0.58% MLMG::computeResidual() 128 0.07797 0.07919 0.07997 0.53% MLMG::mgVcycle_bottom 128 0.0004532 0.005399 0.0794 0.53% MLMG::actualBottomSolve() 128 0.0001794 0.005117 0.07911 0.53% MLMG::mgVcycle_down::2 128 0.01775 0.04704 0.07903 0.53% Castro::create_source_corrector() 17 0.05335 0.06419 0.07808 0.52% MLCGSolver::bicgstab 8 0 0.00482 0.07712 0.52% Castro::construct_new_gravity_source() 17 0.07353 0.07436 0.07545 0.50% Castro::derive() 18 0.07079 0.07133 0.07205 0.48% AmrLevel::derive() 18 0.07074 0.07128 0.072 0.48% MLMG::mgVcycle_down::3 128 0.01721 0.02359 0.06744 0.45% Castro::construct_old_gravity() 20 0.05261 0.05756 0.06358 0.43% Gravity::get_old_grav_vector() 20 0.05249 0.05743 0.06346 0.42% Castro::post_timestep() 1 0.05648 0.05877 0.06137 0.41% FabArray::Xpay() 915 0.05133 0.05333 0.05893 0.39% amrex::average_down 768 0.03486 0.04132 0.05854 0.39% Amr::regrid() 1 0.04311 0.04437 0.05276 0.35% Castro::initialize_advance() 1 0.04778 0.04869 0.04945 0.33% Castro::check_for_nan() 37 0.04825 0.04861 0.04934 0.33% MultiFab::contains_nan() 37 0.04813 0.04849 0.0492 0.33% MLMG::mgVcycle_down::4 128 0.01722 0.01885 0.04132 0.28% MLPoisson::define() 19 0.03027 0.03521 0.04057 0.27% Castro::retry_advance_ctu() 11 0.03711 0.03824 0.0405 0.27% StateDataPhysBCFunct::() 187 0.0009687 0.02077 0.03894 0.26% Gravity::get_new_grav_vector() 18 0.03675 0.03789 0.0389 0.26% StateData::FillBoundary(geom) 193 0 0.01964 0.03774 0.25% MLMG::mgVcycle_up::4 128 0.01528 0.01646 0.03256 0.22% MLMG::mgVcycle_down::5 128 0.01727 0.01815 0.0294 0.20% Castro::apply_source_to_state() 37 0.02856 0.02871 0.02893 0.19% MLCellLinOp::defineAuxData() 19 0.01984 0.02145 0.02666 0.18% VisMF::Write(FabArray) 4 0.01374 0.01503 0.02583 0.17% MLMG::mgVcycle_up::5 128 0.01536 0.01608 0.02506 0.17% MLMG::prepareForSolve() 18 0.02074 0.02167 0.02451 0.16% MLMG::computeMLResidual() 18 0.0109 0.01637 0.02157 0.14% Castro::sum_integrated_quantities() 2 0.02124 0.02139 0.0214 0.14% Castro::save_data_for_retry() 3 0.01944 0.01992 0.02136 0.14% MLMG::ResNormInf() 146 0.01366 0.01617 0.01874 0.13% FabArrayBase::getFB() 6380 0.01402 0.01474 0.01821 0.12% FabArray::setDomainBndry() 187 0.0006464 0.009356 0.01801 0.12% FabArrayBase::getCPC() 2045 0.01084 0.01213 0.01553 0.10% amrex::Add() 128 0.01465 0.0148 0.01551 0.10% FabArray::norminf() 244 0.01103 0.01137 0.01472 0.10% Castro::problem_post_timestep() 2 0.00924 0.01161 0.01438 0.10% Castro::wd_update() 2 0.009228 0.01158 0.01437 0.10% Gravity::test_composite_phi() 1 0.01419 0.01423 0.01427 0.10% Gravity::update_max_rhs() 2 0.001601 0.009065 0.01127 0.08% MLLinOp::define() 19 0.003686 0.007273 0.01087 0.07% Other 2547 0.07946 0.1004 0.1494 1.00% -------------------------------------------------------------------------------------------- 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 405 KiB 410 KiB 416 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 2061 B 2078 B 2098 B 89 KiB 89 KiB 89 KiB -------------------------------------------------------------------------------------------------------------------- Pinned Memory Usage: -------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max -------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 1 B 1 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 1590 KiB 1592 KiB 1592 KiB 2329 KiB 2329 KiB 2329 KiB Castro::construct_ctu_hydro_source() 6400 405 KiB 410 KiB 416 KiB 1176 KiB 1176 KiB 1176 KiB StateData::define() 192 1054 KiB 1055 KiB 1055 KiB 1164 KiB 1164 KiB 1164 KiB Castro::initMFs() 256 1034 KiB 1034 KiB 1035 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 679 KiB 679 KiB 679 KiB 803 KiB 803 KiB 803 KiB Amr::writePlotFile() 64 4303 B 4375 B 4868 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 52 KiB 53 KiB 53 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 438 KiB 438 KiB 438 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2159 B 4008 B 5161 B 225 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1367 B 1379 B 1392 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 18 KiB 18 KiB 18 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 132 KiB 132 KiB 132 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 128 KiB 128 KiB 128 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 9414 B 9964 B 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 1280 499 B 501 B 507 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 420 B 422 B 432 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 2061 B 2078 B 2098 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 17 B 20 B 36 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 33 KiB 33 KiB 33 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 49 B 49 B 50 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 3 B 11 B 14 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 110 B 114 B 117 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 160 B 175 B 194 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 31 KiB 31 KiB 31 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 223 B 229 B 233 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 1071 B 1075 B 1085 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 124 B 124 B 124 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 1615 B 1679 B 1705 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 0 B 0 B 1 B 16 KiB 16 KiB 16 KiB MLCellLinOp::defineAuxData() 11704 425 B 475 B 659 B 4608 B 5151 B 7152 B BndryData::define() 4864 442 B 443 B 445 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 22 B 31 B 41 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 0 B 17 B 28 B 2916 B 2916 B 2916 B OwnerMask() 64 0 B 0 B 0 B 2592 B 2592 B 2592 B MLMG::addInterpCorrection() 5632 1 B 15 B 24 B 1024 B 1024 B 1024 B amrex::average_down 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 0 B 5 B 7 B 648 B 648 B 648 B MLCGSolver::bicgstab 640 0 B 0 B 2 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 0 B 0 B 0 B 8 B 8 B 8 B ----------------------------------------------------------------------------------------------------------------------- AMReX (25.12-15-g4fc753db6b0b) finalized