Initializing AMReX (26.01-5-g005066b9e5ef)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.01-5-g005066b9e5ef) initialized Starting run at 07:51:52 UTC on 2026-01-08. Successfully read inputs file ... Castro git describe: 26.01-1-ge6f237912 AMReX git describe: 26.01-5-g005066b9e Microphysics git describe: 26.01-4-g935d36be 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.186733835 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001191168 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010989236 Gravity::multilevel_solve_for_new_phi() time = 0.108073663 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010730056 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.054281476 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.007538048 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 = -0.01295059363 at index 21, 172, 0 Invalid density = -8.076810125 at index 20, 174, 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, Invalid density = -15.15828486 at index 20, 175, 0 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 = -108.2035502 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477044 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.769409 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071148 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071148 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.769409 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477044 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035502 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -0.01295059363 at index 21, 83, 0 Castro::construct_old_gravity() time = 0.002765792 on level 0 Castro::do_old_sources() time = 0.007888103 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.0623919 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.002971096 on level 0 Castro::do_old_sources() time = 0.007920178 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.061789892 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.003647918 on level 0 Castro::do_old_sources() time = 0.007881463 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.057927332 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.00315747 on level 0 Castro::do_old_sources() time = 0.008735991 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.05690146 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.33427894 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.015162305 Gravity::solve_for_phi() time = 0.091387624 on level 0 Castro::construct_new_gravity() time = 0.093533675 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.057540646 on level 0 Castro::construct_old_gravity() time = 0.00281108 on level 0 Castro::do_old_sources() time = 0.007848639 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.056828536 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.335623833 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.015088796 Gravity::solve_for_phi() time = 0.091809004 on level 0 Castro::construct_new_gravity() time = 0.093956579 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.07151753 on level 0 Castro::construct_old_gravity() time = 0.022848865 on level 0 Castro::do_old_sources() time = 0.008843339 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.058829689 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.348443971 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.014893755 Gravity::solve_for_phi() time = 0.09126425 on level 0 Castro::construct_new_gravity() time = 0.093430092 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.056447161 on level 0 Castro::construct_old_gravity() time = 0.002949885 on level 0 Castro::do_old_sources() time = 0.008181011 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.056601912 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.345032046 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.015287943 Gravity::solve_for_phi() time = 0.091941295 on level 0 Castro::construct_new_gravity() time = 0.094086433 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.056661625 on level 0 Castro::construct_old_gravity() time = 0.003058924 on level 0 Castro::do_old_sources() time = 0.008649761 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.056517466 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.34781002 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.014939161 Gravity::solve_for_phi() time = 0.09168934 on level 0 Castro::construct_new_gravity() time = 0.093837718 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.057242702 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.00303238 on level 0 Castro::do_old_sources() time = 0.007939369 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.090172149 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.36467826 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.015888736 Gravity::solve_for_phi() time = 0.100342247 on level 0 Castro::construct_new_gravity() time = 0.102515703 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.05746375 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.003054793 on level 0 Castro::do_old_sources() time = 0.007984303 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.056731782 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.355087221 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.015622133 Gravity::solve_for_phi() time = 0.09333552 on level 0 Castro::construct_new_gravity() time = 0.096225175 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.057194113 on level 0 Castro::construct_old_gravity() time = 0.002903942 on level 0 Castro::do_old_sources() time = 0.008081349 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.05674885 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.345563921 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.017556334 Gravity::solve_for_phi() time = 0.093723744 on level 0 Castro::construct_new_gravity() time = 0.095861258 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.057381348 on level 0 Castro::construct_old_gravity() time = 0.002873346 on level 0 Castro::do_old_sources() time = 0.007933764 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.056491239 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.34818732 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.015780011 Gravity::solve_for_phi() time = 0.092097371 on level 0 Castro::construct_new_gravity() time = 0.094239922 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.056758154 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.002796049 on level 0 Castro::do_old_sources() time = 0.007819087 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.056730754 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.344398746 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.015826295 Gravity::solve_for_phi() time = 0.092154862 on level 0 Castro::construct_new_gravity() time = 0.094291219 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.056995129 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.002789195 on level 0 Castro::do_old_sources() time = 0.007935815 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.056726157 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.349985833 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.01638791 Gravity::solve_for_phi() time = 0.092627095 on level 0 Castro::construct_new_gravity() time = 0.094800267 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.057179517 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.002906601 on level 0 Castro::do_old_sources() time = 0.00811849 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.056533349 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.345556379 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.016574755 Gravity::solve_for_phi() time = 0.09289206 on level 0 Castro::construct_new_gravity() time = 0.095060499 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.056778516 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.00281166 on level 0 Castro::do_old_sources() time = 0.00790413 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.056759504 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.352185603 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.016110989 Gravity::solve_for_phi() time = 0.092953906 on level 0 Castro::construct_new_gravity() time = 0.095176767 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.057022994 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.002889882 on level 0 Castro::do_old_sources() time = 0.007989243 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.056611637 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.344232101 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.015879821 Gravity::solve_for_phi() time = 0.092516977 on level 0 Castro::construct_new_gravity() time = 0.094663566 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.056890411 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.002991033 on level 0 Castro::do_old_sources() time = 0.008075091 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.056808414 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.350820919 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.017596166 Gravity::solve_for_phi() time = 0.094196835 on level 0 Castro::construct_new_gravity() time = 0.096350449 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.05817896 on level 0 Castro::construct_old_gravity() time = 0.003828458 on level 0 Castro::do_old_sources() time = 0.007915629 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.056775636 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.345712218 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.016548933 Gravity::solve_for_phi() time = 0.093151766 on level 0 Castro::construct_new_gravity() time = 0.095372285 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.057459663 on level 0 Castro::construct_old_gravity() time = 0.010374854 on level 0 Castro::do_old_sources() time = 0.008476716 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.056713888 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.350347596 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002603737404 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.016557401 Gravity::solve_for_phi() time = 0.095738107 on level 0 Castro::construct_new_gravity() time = 0.097950735 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.057278708 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.64136585 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.053715606 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.005498801 seconds Ending run at 07:52:07 UTC on 2026-01-08. Run time = 15.09458513 Run time without initialization = 12.70156495 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]: 15.0841 ... 15.0864 ... 15.0947 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.8095 5.8457 5.8872 39.00% Amr::InitAmr() 1 1.364 1.366 1.374 9.10% Castro::computeTemp() 80 0.9188 1.05 1.245 8.25% amrex::Copy() 766 0.8564 0.8624 0.8716 5.77% FabArray::ParallelCopy_finish() 1469 0.3208 0.6508 0.8662 5.74% MultiFab::AddProduct() 3598 0.5878 0.5895 0.5933 3.93% Castro::reset_internal_energy(Fab) 160 0.4462 0.453 0.4622 3.06% MLPoisson::Fsmooth() 6144 0.3531 0.3818 0.4557 3.02% FabArray::setVal() 5232 0.3963 0.4004 0.4062 2.69% Castro::initMFs() 1 0.3834 0.392 0.4055 2.69% MultiFab::Subtract() 159 0.3584 0.3597 0.3638 2.41% Castro::do_advance_ctu() 20 0.2938 0.3002 0.3341 2.21% Gravity::fill_multipole_BCs() 19 0.2027 0.2648 0.2899 1.92% FabArray::ParallelCopy_nowait() 1469 0.2431 0.2575 0.2762 1.83% FabArray::sum() 3598 0.1982 0.1991 0.2006 1.33% check_for_negative_density() 20 0.08409 0.132 0.1721 1.14% Castro::initData() 1 0.003379 0.02867 0.1702 1.13% TagBoxArray::collate() 2 0.0002156 0.036 0.1652 1.09% FillBoundary_finish() 6380 0.0983 0.1234 0.1567 1.04% FillBoundary_nowait() 6380 0.09044 0.1161 0.148 0.98% MLCellLinOp::applyBC() 7148 0.06663 0.0826 0.1289 0.85% Castro::normalize_species() 79 0.1151 0.1155 0.1171 0.78% MLPoisson::Fapply() 986 0.09537 0.09956 0.1151 0.76% Castro::construct_old_gravity_source() 40 0.1015 0.1023 0.1035 0.69% FabArray::Saxpy() 111 0.08525 0.08594 0.08957 0.59% Castro::sum_integrated_quantities() 2 0.06373 0.07806 0.07929 0.53% Castro::construct_new_gravity_source() 17 0.07293 0.07346 0.07391 0.49% FabArray::Xpay() 915 0.05126 0.05329 0.0588 0.39% MultiFab::contains_nan() 37 0.04902 0.04936 0.04988 0.33% AmrLevel::derive() 18 0.04462 0.04553 0.04833 0.32% MLMG::addInterpCorrection() 768 0.03385 0.03643 0.04386 0.29% StateData::FillBoundary(geom) 193 0 0.02006 0.03951 0.26% amrex::average_down 768 0.02516 0.02726 0.0337 0.22% VisMF::Write(FabArray) 4 0.01572 0.01742 0.02979 0.20% Castro::reset_internal_energy(MultiFab) 80 0.0243 0.02499 0.02603 0.17% Amr::writePlotFile() 2 0.008191 0.01884 0.02097 0.14% Castro::construct_old_source() 240 0.02049 0.02066 0.0208 0.14% Castro::enforce_min_density() 79 0.01937 0.01977 0.02015 0.13% MLCellLinOp::defineAuxData() 19 0.0156 0.01671 0.02003 0.13% FabArray::FillBoundary() 6380 0.01647 0.01692 0.01984 0.13% FabArray::setDomainBndry() 187 0.0006548 0.00912 0.01749 0.12% MLCGSolver::bicgstab 8 0 0.001058 0.01693 0.11% FabArrayBase::getFB() 6380 0.01353 0.01397 0.01684 0.11% FillPatchIterator::Initialize 187 0.01543 0.0159 0.01659 0.11% Castro::construct_new_source() 102 0.01616 0.01627 0.01643 0.11% Castro::evaluate_source_change() 196 0.01481 0.01515 0.01592 0.11% FabArray::norminf() 244 0.01106 0.01138 0.01479 0.10% amrex::Add() 128 0.01458 0.01465 0.01478 0.10% Castro::wd_update() 2 0.004693 0.009371 0.01281 0.08% MLCellLinOp::smooth() 1544 0.01064 0.01097 0.01182 0.08% MLMG::prepareForSolve() 18 0.01029 0.01072 0.01182 0.08% Gravity::actual_solve_with_mlmg() 19 0.007459 0.007936 0.009584 0.06% FabArrayBase::getCPC() 2045 0.008486 0.008738 0.009084 0.06% Amr::writeSmallPlotFile() 2 0.00626 0.008324 0.008751 0.06% Gravity::get_old_grav_vector() 20 0.006624 0.006714 0.007028 0.05% MLMG::ResNormInf() 146 0.004256 0.005576 0.006944 0.05% MLMG::solve() 18 0.0009423 0.003704 0.006395 0.04% Gravity::get_new_grav_vector() 18 0.006055 0.006157 0.006338 0.04% Castro::initialize_do_advance() 20 0.004809 0.005373 0.005955 0.04% MLCellLinOp::apply() 986 0.00239 0.002629 0.005575 0.04% FabArray::ParallelCopy() 1469 0.004209 0.004329 0.00456 0.03% AmrLevel::FillPatch() 187 0.0031 0.003232 0.003485 0.02% MLMG::mgVcycle() 128 0.002236 0.002284 0.002367 0.02% MLCellLinOp::correctionResidual() 768 0.002224 0.002269 0.002317 0.02% MLMG:computeResOfCorrection() 768 0.001767 0.001797 0.001865 0.01% Castro::save_data_for_retry() 3 0.001331 0.001442 0.001522 0.01% StateDataPhysBCFunct::() 187 0.0009855 0.001167 0.001326 0.01% MLMG::actualBottomSolve() 128 0.0001801 0.0002505 0.001242 0.01% Castro::subcycle_advance_ctu() 1 0.0002454 0.0003444 0.001115 0.01% FillPatchIterator::FillFromLevel0() 187 0.0008889 0.0009321 0.001047 0.01% main() 1 0.000725 0.0008321 0.0009772 0.01% Gravity::solve_for_phi() 17 0.000664 0.0007109 0.0007845 0.01% FillPatchSingleLevel 187 0.0006643 0.0006828 0.0007185 0.00% MLMG::mgVcycle_down::0 128 0.0006251 0.0006441 0.0006705 0.00% MLMG::mgVcycle_down::2 128 0.0005796 0.0006132 0.000652 0.00% MLMG::mgVcycle_down::1 128 0.0006063 0.0006229 0.0006494 0.00% MLMG::mgVcycle_down::3 128 0.0005526 0.0005837 0.0006474 0.00% MLMG::mgVcycle_down::4 128 0.0005439 0.000563 0.0006261 0.00% MLMG::mgVcycle_down::5 128 0.0005748 0.0005886 0.0006191 0.00% AmrMesh::MakeNewGrids() 2 0.0005434 0.0005638 0.0006113 0.00% Castro::Castro() 1 0.0004095 0.0004795 0.0005623 0.00% Castro::clean_state() 79 0.0004671 0.0005032 0.0005558 0.00% MLCellLinOp::solutionResidual() 147 0.0005064 0.0005176 0.0005366 0.00% Castro::do_old_sources() 40 0.0004586 0.0004753 0.000494 0.00% MLMG::oneIter() 128 0.0004283 0.0004413 0.0004648 0.00% MLMG::mgVcycle_up::2 128 0.0003841 0.0004063 0.0004398 0.00% MLMG::mgVcycle_up::0 128 0.0004041 0.0004139 0.0004388 0.00% MLMG::mgVcycle_up::1 128 0.0004015 0.0004124 0.0004298 0.00% MLMG::mgVcycle_up::4 128 0.0003787 0.000387 0.0003996 0.00% MLMG::mgVcycle_up::3 128 0.0003715 0.0003811 0.0003975 0.00% MLMG::mgVcycle_up::5 128 0.0003662 0.000372 0.0003873 0.00% Castro::problem_post_timestep() 2 1.153e-05 3.641e-05 0.0003867 0.00% Castro::construct_new_gravity() 17 0.000163 0.0001924 0.0003647 0.00% MLMG::computeResidual() 128 0.0002885 0.000296 0.0003102 0.00% MLMG::mgVcycle_bottom 128 0.0002715 0.0002789 0.0002895 0.00% Castro::do_new_sources() 17 0.0002132 0.0002267 0.0002613 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002174 0.0002339 0.0002527 0.00% Castro::initialize_advance() 1 0.0001477 0.0001652 0.0001775 0.00% Amr::defBaseLevel() 1 0.0001089 0.0001494 0.0001721 0.00% Castro::check_for_nan() 37 0.0001177 0.0001292 0.0001527 0.00% Amr::coarseTimeStep() 1 0.0001244 0.0001474 0.0001517 0.00% MLPoisson::define() 19 0.0001211 0.0001313 0.0001472 0.00% Castro::retry_advance_ctu() 11 6.802e-05 8.224e-05 0.0001465 0.00% Castro::construct_old_gravity() 20 0.0001207 0.0001277 0.0001452 0.00% TagBoxArray::mapPRD 2 0.0001162 0.0001237 0.0001368 0.00% Castro::apply_source_to_state() 37 0.0001182 0.0001274 0.0001361 0.00% Gravity::test_composite_phi() 1 0.0001041 0.0001136 0.000123 0.00% Amr::grid_places() 2 5.274e-05 5.947e-05 0.0001046 0.00% MLLinOp::define() 19 8.17e-05 9.046e-05 0.0001041 0.00% Castro::expand_state() 20 5.946e-05 6.652e-05 7.245e-05 0.00% Castro::create_source_corrector() 17 5.517e-05 6.167e-05 6.879e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.928e-05 6.119e-05 6.752e-05 0.00% MLMG::computeMLResidual() 18 5.184e-05 5.605e-05 6.592e-05 0.00% Castro::derive() 18 4.605e-05 4.993e-05 5.793e-05 0.00% Amr::FinalizeInit() 1 9.019e-06 1.387e-05 5.716e-05 0.00% Amr::regrid() 1 1.713e-05 2.144e-05 4.684e-05 0.00% Gravity::actual_multilevel_solve() 1 3.149e-05 3.318e-05 3.686e-05 0.00% Amr::init() 1 7.069e-06 1.545e-05 2.853e-05 0.00% Amr::timeStep() 1 1.55e-05 1.823e-05 2.424e-05 0.00% Castro::advance() 1 8.635e-06 1.123e-05 1.953e-05 0.00% Amr::InitializeInit() 1 8.444e-06 1.178e-05 1.94e-05 0.00% Castro::post_regrid() 2 1.092e-05 1.205e-05 1.728e-05 0.00% Castro::post_init() 1 1.146e-05 1.261e-05 1.435e-05 0.00% Amr::bldFineLevels() 1 7.302e-06 7.996e-06 8.883e-06 0.00% Castro::post_timestep() 1 6.017e-06 6.686e-06 7.463e-06 0.00% Amr::initialInit() 1 3.274e-06 3.947e-06 6.304e-06 0.00% Other 2547 0.0456 0.05666 0.09692 0.64% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 15.08 15.09 15.09 100.00% Amr::coarseTimeStep() 1 12.7 12.7 12.7 84.14% Amr::timeStep() 1 12.64 12.64 12.64 83.75% Castro::advance() 1 12.58 12.58 12.59 83.38% Castro::subcycle_advance_ctu() 1 12.53 12.54 12.54 83.07% Castro::do_advance_ctu() 20 12.49 12.5 12.5 82.80% Castro::construct_ctu_hydro_source() 20 6.666 6.703 6.772 44.86% Castro::clean_state() 79 3.35 3.482 3.692 24.46% Castro::computeTemp() 80 2.414 2.546 2.749 18.21% Gravity::solve_phi_with_mlmg() 19 1.609 1.674 1.702 11.28% Castro::construct_new_gravity() 17 1.529 1.593 1.621 10.74% Gravity::solve_for_phi() 17 1.49 1.555 1.583 10.49% Castro::reset_internal_energy(MultiFab) 80 1.487 1.496 1.508 9.99% Gravity::actual_solve_with_mlmg() 19 1.402 1.406 1.41 9.34% Castro::initialize_do_advance() 20 1.295 1.347 1.396 9.25% Amr::InitAmr() 1 1.364 1.366 1.374 9.10% MLMG::solve() 18 1.347 1.352 1.354 8.97% Castro::do_old_sources() 40 1.26 1.294 1.312 8.69% MLMG::oneIter() 128 1.215 1.216 1.218 8.07% MLMG::mgVcycle() 128 1.2 1.201 1.203 7.97% FabArray::ParallelCopy() 1469 0.6076 0.9345 1.152 7.63% Castro::evaluate_source_change() 196 1.018 1.022 1.029 6.82% Amr::init() 1 1.019 1.019 1.019 6.75% Castro::do_new_sources() 17 0.8659 0.9217 0.9853 6.53% Amr::initialInit() 1 0.9561 0.9562 0.9563 6.34% amrex::Copy() 766 0.8564 0.8624 0.8716 5.77% FabArray::ParallelCopy_finish() 1469 0.3208 0.6508 0.8662 5.74% AmrLevel::FillPatch() 187 0.6282 0.7772 0.8595 5.69% Castro::enforce_min_density() 79 0.8442 0.8485 0.8562 5.67% MLCellLinOp::smooth() 1544 0.6262 0.6793 0.8149 5.40% Amr::InitializeInit() 1 0.4695 0.502 0.6549 4.34% Amr::defBaseLevel() 1 0.4695 0.502 0.6549 4.34% FillPatchIterator::Initialize 187 0.4103 0.5636 0.6465 4.28% FillPatchIterator::FillFromLevel0() 187 0.3763 0.5386 0.6158 4.08% FillPatchSingleLevel 187 0.3753 0.5377 0.6149 4.07% MultiFab::AddProduct() 3598 0.5878 0.5895 0.5933 3.93% Amr::FinalizeInit() 1 0.3012 0.4542 0.4867 3.22% Castro::reset_internal_energy(Fab) 160 0.4462 0.453 0.4622 3.06% MLMG::mgVcycle_up::1 128 0.1065 0.2815 0.4574 3.03% MLPoisson::Fsmooth() 6144 0.3531 0.3818 0.4557 3.02% MLCellLinOp::applyBC() 7148 0.3094 0.3536 0.4372 2.90% Castro::Castro() 1 0.3895 0.3981 0.4118 2.73% FabArray::setVal() 5232 0.3963 0.4004 0.4062 2.69% Castro::initMFs() 1 0.3836 0.3922 0.4057 2.69% MLMG::addInterpCorrection() 768 0.06732 0.3231 0.3975 2.63% MultiFab::Subtract() 159 0.3584 0.3597 0.3638 2.41% FabArray::FillBoundary() 6380 0.2368 0.271 0.3129 2.07% MLMG::mgVcycle_up::2 128 0.01542 0.1279 0.3004 1.99% FabArray::ParallelCopy_nowait() 1469 0.2644 0.2794 0.2986 1.98% Gravity::fill_multipole_BCs() 19 0.2049 0.267 0.292 1.93% Castro::expand_state() 20 0.2108 0.2481 0.2912 1.93% MLMG::mgVcycle_down::0 128 0.2747 0.2762 0.278 1.84% Castro::initData() 1 0.07593 0.1016 0.2436 1.61% Amr::bldFineLevels() 1 0.0426 0.1955 0.2276 1.51% Castro::post_init() 1 0.2175 0.219 0.2198 1.46% MLMG::mgVcycle_up::0 128 0.2091 0.2105 0.2126 1.41% MLCellLinOp::apply() 986 0.1553 0.1673 0.2087 1.38% FabArray::sum() 3598 0.1982 0.1991 0.2006 1.33% MLMG::mgVcycle_up::3 128 0.01533 0.02924 0.1964 1.30% Amr::grid_places() 2 0.002432 0.1554 0.1879 1.24% AmrMesh::MakeNewGrids() 2 0.002379 0.1554 0.1879 1.24% TagBoxArray::mapPRD 2 0.0007846 0.118 0.1862 1.23% check_for_negative_density() 20 0.08409 0.132 0.1721 1.14% Castro::construct_old_source() 240 0.1689 0.1701 0.1712 1.13% TagBoxArray::collate() 2 0.0002156 0.036 0.1652 1.09% FillBoundary_nowait() 6380 0.1049 0.1307 0.1631 1.08% FillBoundary_finish() 6380 0.0983 0.1234 0.1567 1.04% MLMG:computeResOfCorrection() 768 0.1183 0.1291 0.1535 1.02% MLCellLinOp::correctionResidual() 768 0.1165 0.1273 0.1516 1.00% Castro::construct_new_source() 102 0.1292 0.1299 0.1304 0.86% MLMG::mgVcycle_down::1 128 0.1221 0.1258 0.1295 0.86% Castro::normalize_species() 79 0.1151 0.1155 0.1171 0.78% MLPoisson::Fapply() 986 0.09537 0.09956 0.1151 0.76% Castro::create_source_corrector() 17 0.05922 0.08532 0.1128 0.75% Amr::writePlotFile() 2 0.1085 0.1086 0.1086 0.72% Gravity::multilevel_solve_for_new_phi() 1 0.108 0.108 0.1081 0.72% Gravity::actual_multilevel_solve() 1 0.1079 0.108 0.108 0.72% Castro::construct_old_gravity_source() 40 0.1015 0.1023 0.1035 0.69% MLCellLinOp::solutionResidual() 147 0.08828 0.09412 0.09956 0.66% Castro::sum_integrated_quantities() 2 0.09767 0.09795 0.09797 0.65% FabArray::Saxpy() 111 0.08525 0.08594 0.08957 0.59% Castro::derive() 18 0.07523 0.07687 0.0882 0.58% AmrLevel::derive() 18 0.07518 0.07682 0.08815 0.58% Castro::construct_old_gravity() 20 0.04772 0.06163 0.08196 0.54% Gravity::get_old_grav_vector() 20 0.04759 0.0615 0.08183 0.54% MLMG::mgVcycle_down::2 128 0.01782 0.04801 0.08128 0.54% Castro::post_regrid() 2 0.07889 0.07974 0.08126 0.54% MLMG::mgVcycle_bottom 128 0.0004551 0.005404 0.07952 0.53% MLMG::actualBottomSolve() 128 0.0001801 0.005125 0.07923 0.52% MLMG::computeResidual() 128 0.07564 0.07725 0.0783 0.52% MLCGSolver::bicgstab 8 0 0.004827 0.07724 0.51% Castro::construct_new_gravity_source() 17 0.07293 0.07346 0.07391 0.49% MLMG::mgVcycle_down::3 128 0.01723 0.02383 0.06832 0.45% Castro::post_timestep() 1 0.05539 0.05863 0.06185 0.41% amrex::average_down 768 0.03493 0.04171 0.05921 0.39% FabArray::Xpay() 915 0.05126 0.05329 0.0588 0.39% Castro::check_for_nan() 37 0.04915 0.04949 0.05002 0.33% MultiFab::contains_nan() 37 0.04902 0.04936 0.04988 0.33% Castro::initialize_advance() 1 0.04599 0.04682 0.04755 0.32% Amr::regrid() 1 0.04083 0.04122 0.04254 0.28% MLPoisson::define() 19 0.03261 0.03615 0.04118 0.27% MLMG::mgVcycle_down::4 128 0.01726 0.01901 0.04116 0.27% StateDataPhysBCFunct::() 187 0.0009855 0.02123 0.04083 0.27% StateData::FillBoundary(geom) 193 0 0.02006 0.03951 0.26% Castro::retry_advance_ctu() 11 0.03752 0.03834 0.0394 0.26% Gravity::get_new_grav_vector() 18 0.03698 0.03753 0.03799 0.25% Castro::apply_source_to_state() 37 0.02856 0.02898 0.03264 0.22% MLMG::mgVcycle_up::4 128 0.0154 0.01668 0.03262 0.22% MLMG::mgVcycle_down::5 128 0.01722 0.01826 0.02995 0.20% VisMF::Write(FabArray) 4 0.01572 0.01742 0.02979 0.20% MLCellLinOp::defineAuxData() 19 0.02012 0.02198 0.02723 0.18% MLMG::mgVcycle_up::5 128 0.01544 0.01624 0.02516 0.17% MLMG::prepareForSolve() 18 0.021 0.02203 0.02493 0.17% MLMG::computeMLResidual() 18 0.01087 0.01644 0.02161 0.14% Castro::save_data_for_retry() 3 0.01956 0.01984 0.02053 0.14% FabArrayBase::getFB() 6380 0.01385 0.01465 0.0183 0.12% FabArray::setDomainBndry() 187 0.0006548 0.00912 0.01749 0.12% MLMG::ResNormInf() 146 0.01406 0.01532 0.01675 0.11% Castro::problem_post_timestep() 2 0.009254 0.01287 0.016 0.11% Castro::wd_update() 2 0.009241 0.01284 0.01598 0.11% FabArrayBase::getCPC() 2045 0.01096 0.01219 0.01561 0.10% FabArray::norminf() 244 0.01106 0.01138 0.01479 0.10% amrex::Add() 128 0.01458 0.01465 0.01478 0.10% Gravity::test_composite_phi() 1 0.01436 0.01441 0.01445 0.10% Amr::writeSmallPlotFile() 2 0.01328 0.01328 0.01329 0.09% MLLinOp::define() 19 0.003664 0.007596 0.01126 0.07% Other 2547 0.07372 0.09192 0.1438 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 409 KiB 412 KiB 415 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 2202 B 2349 B 2468 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 4 B 4 B 8 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 1588 KiB 1589 KiB 1589 KiB 2329 KiB 2329 KiB 2329 KiB Castro::construct_ctu_hydro_source() 6400 409 KiB 412 KiB 415 KiB 1176 KiB 1176 KiB 1176 KiB StateData::define() 192 1057 KiB 1058 KiB 1058 KiB 1164 KiB 1164 KiB 1164 KiB Castro::initMFs() 256 1036 KiB 1037 KiB 1037 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 675 KiB 675 KiB 675 KiB 803 KiB 803 KiB 803 KiB Amr::writePlotFile() 64 4532 B 4646 B 5230 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 51 KiB 52 KiB 52 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 12 KiB 15 KiB 16 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 435 KiB 435 KiB 435 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2272 B 4625 B 6117 B 225 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1464 B 1484 B 1552 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 131 KiB 131 KiB 131 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 129 KiB 129 KiB 129 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 9355 B 9911 B 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 1280 497 B 500 B 502 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 420 B 422 B 424 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 2202 B 2349 B 2468 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 25 B 27 B 44 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 50 B 50 B 50 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 3 B 4 B 6 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 110 B 111 B 113 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 143 B 185 B 248 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 224 B 232 B 238 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 1099 B 1102 B 1110 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 117 B 117 B 117 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 1615 B 1686 B 1716 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 423 B 473 B 656 B 4608 B 5151 B 7152 B BndryData::define() 4864 440 B 441 B 442 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 21 B 30 B 41 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 0 B 20 B 32 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 6 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 (26.01-5-g005066b9e5ef) finalized