Initializing AMReX (26.03-110-g388b0db16adb)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.03-110-g388b0db16adb) initialized Starting run at 06:49:57 UTC on 2026-03-31. Successfully read inputs file ... Castro git describe: 26.03-32-g1f44b02a5 AMReX git describe: 26.03-110-g388b0db16 Microphysics git describe: 26.03-15-g0f912491 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.163964027 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001321382 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010366618 Gravity::multilevel_solve_for_new_phi() time = 0.104713605 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.009918135 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.04648466 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004183239 seconds [Level 0 step 1] ADVANCE at time 0 with dt = 0.25  Beginning subcycle 1 starting at time 0 with dt = 0.25 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.7670713549 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.7670713549 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071148 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.769409 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477044 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035502 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -0.01295059363 at index 21, 83, 0 Invalid density = -0.01295059363 at index 21, 172, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035502 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477044 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.769409 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071148 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Castro::construct_old_gravity() time = 0.00267028 on level 0 Castro::do_old_sources() time = 0.007320843 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.052759181 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Castro::construct_old_gravity() time = 0.003359503 on level 0 Castro::do_old_sources() time = 0.009010215 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.055007983 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.125 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.0625   Beginning subcycle 1 starting at time 0 with dt = 0.0625 Estimated number of subcycles remaining: 4 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.1060744788 at index 20, 78, 0 Invalid density = -0.1060744788 at index 20, 177, 0 Castro::construct_old_gravity() time = 0.002605587 on level 0 Castro::do_old_sources() time = 0.007313302 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.052908274 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.002914186 on level 0 Castro::do_old_sources() time = 0.007312488 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.04862642 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.290885462 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.012825403 Gravity::solve_for_phi() time = 0.087616228 on level 0 Castro::construct_new_gravity() time = 0.089633217 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.04886955 on level 0 Castro::construct_old_gravity() time = 0.002536832 on level 0 Castro::do_old_sources() time = 0.007380958 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.048364877 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.290303825 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.012622199 Gravity::solve_for_phi() time = 0.086815787 on level 0 Castro::construct_new_gravity() time = 0.088796977 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.048592089 on level 0 Castro::construct_old_gravity() time = 0.002639136 on level 0 Castro::do_old_sources() time = 0.00734736 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.048743275 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.296606925 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.012621756 Gravity::solve_for_phi() time = 0.087154164 on level 0 Castro::construct_new_gravity() time = 0.089178341 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.04852081 on level 0 Castro::construct_old_gravity() time = 0.002786701 on level 0 Castro::do_old_sources() time = 0.007314948 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.048381266 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.301580843 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.013344088 Gravity::solve_for_phi() time = 0.087501896 on level 0 Castro::construct_new_gravity() time = 0.089460318 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.048573121 on level 0 Castro::construct_old_gravity() time = 0.002655806 on level 0 Castro::do_old_sources() time = 0.007266044 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.048579214 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.305918656 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.012917059 Gravity::solve_for_phi() time = 0.087168297 on level 0 Castro::construct_new_gravity() time = 0.089162137 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.048648172 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.002603764 on level 0 Castro::do_old_sources() time = 0.007368967 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.048143851 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.301810532 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.013021367 Gravity::solve_for_phi() time = 0.088480391 on level 0 Castro::construct_new_gravity() time = 0.090503574 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.049493491 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.002786934 on level 0 Castro::do_old_sources() time = 0.007435581 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.048504508 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.305095977 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.013210539 Gravity::solve_for_phi() time = 0.08743931 on level 0 Castro::construct_new_gravity() time = 0.089454804 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.048234772 on level 0 Castro::construct_old_gravity() time = 0.002981676 on level 0 Castro::do_old_sources() time = 0.007591356 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.048736217 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.299906537 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.013017754 Gravity::solve_for_phi() time = 0.087518979 on level 0 Castro::construct_new_gravity() time = 0.089477681 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.048820721 on level 0 Castro::construct_old_gravity() time = 0.002660453 on level 0 Castro::do_old_sources() time = 0.007283108 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.048533774 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.30491348 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.012902069 Gravity::solve_for_phi() time = 0.093370602 on level 0 Castro::construct_new_gravity() time = 0.095344522 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.048565475 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.00262349 on level 0 Castro::do_old_sources() time = 0.007309963 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.048613612 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.300213508 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.013044977 Gravity::solve_for_phi() time = 0.087235369 on level 0 Castro::construct_new_gravity() time = 0.089222403 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.048591949 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.002618892 on level 0 Castro::do_old_sources() time = 0.007382255 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.048443894 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.308761315 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.013215325 Gravity::solve_for_phi() time = 0.087718257 on level 0 Castro::construct_new_gravity() time = 0.089747044 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.048710612 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.002528209 on level 0 Castro::do_old_sources() time = 0.007459107 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.048200166 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.30490773 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.013230126 Gravity::solve_for_phi() time = 0.087578988 on level 0 Castro::construct_new_gravity() time = 0.089559991 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.048855441 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.002678621 on level 0 Castro::do_old_sources() time = 0.007270806 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.048263962 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.312997058 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.013358313 Gravity::solve_for_phi() time = 0.087454826 on level 0 Castro::construct_new_gravity() time = 0.089485462 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.048787672 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.002568794 on level 0 Castro::do_old_sources() time = 0.007354636 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.048272895 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.300792039 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.013457821 Gravity::solve_for_phi() time = 0.088920842 on level 0 Castro::construct_new_gravity() time = 0.090899149 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.049475993 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.003553735 on level 0 Castro::do_old_sources() time = 0.007387974 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.048244805 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.307287574 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.01342884 Gravity::solve_for_phi() time = 0.087513282 on level 0 Castro::construct_new_gravity() time = 0.089511272 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.048899512 on level 0 Castro::construct_old_gravity() time = 0.011854574 on level 0 Castro::do_old_sources() time = 0.01597569 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.048287004 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.307253233 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.013633037 Gravity::solve_for_phi() time = 0.088460709 on level 0 Castro::construct_new_gravity() time = 0.090442324 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.04887529 on level 0 Castro::construct_old_gravity() time = 0.002673144 on level 0 Castro::do_old_sources() time = 0.007613322 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.048041034 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.30641977 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002971962068 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.01358212 Gravity::solve_for_phi() time = 0.088520249 on level 0 Castro::construct_new_gravity() time = 0.090559543 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.048881792 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: 11.07563698 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.046356552 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.005004165 seconds Ending run at 06:50:11 UTC on 2026-03-31. Run time = 14.09924429 Run time without initialization = 11.12799627 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.0918 ... 14.0940 ... 14.0993 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.1074 5.1334 5.1874 36.79% Amr::InitAmr() 1 2.041 2.043 2.048 14.53% Castro::computeTemp() 80 0.7027 0.7905 0.9138 6.48% amrex::Copy() 766 0.7591 0.7637 0.773 5.48% FabArray::ParallelCopy_finish() 1469 0.2126 0.5135 0.6874 4.88% MultiFab::AddProduct() 3598 0.5401 0.5427 0.5514 3.91% MLPoisson::Fsmooth() 6144 0.3303 0.3565 0.4283 3.04% Castro::initMFs() 1 0.3824 0.3842 0.3897 2.76% FabArray::setVal() 5232 0.3614 0.3643 0.3697 2.62% Castro::reset_internal_energy(Fab) 160 0.343 0.3445 0.3488 2.47% MultiFab::Subtract() 159 0.3385 0.3397 0.3426 2.43% Castro::do_advance_ctu() 20 0.2596 0.2626 0.2667 1.89% FabArray::ParallelCopy_nowait() 1469 0.2211 0.2412 0.2622 1.86% Gravity::fill_multipole_BCs() 19 0.1853 0.224 0.2394 1.70% FabArray::sum() 3598 0.1828 0.1844 0.1865 1.32% Castro::initData() 1 0.003233 0.02822 0.1656 1.17% FillBoundary_finish() 6380 0.1061 0.1267 0.1637 1.16% TagBoxArray::collate() 2 0.0001904 0.02818 0.1586 1.13% FillBoundary_nowait() 6380 0.08827 0.1185 0.152 1.08% check_for_negative_density() 20 0.06173 0.123 0.1464 1.04% MLCellLinOp::applyBC() 7148 0.06887 0.08425 0.1299 0.92% Castro::normalize_species() 79 0.119 0.1195 0.1205 0.85% MLPoisson::Fapply() 986 0.0867 0.09048 0.1052 0.75% Castro::construct_old_gravity_source() 40 0.09429 0.09506 0.09616 0.68% FabArray::Saxpy() 111 0.07934 0.07979 0.08066 0.57% Castro::construct_new_gravity_source() 17 0.06664 0.06724 0.06775 0.48% FabArray::Xpay() 915 0.04829 0.05011 0.05527 0.39% MultiFab::contains_nan() 37 0.04542 0.04572 0.04621 0.33% Castro::sum_integrated_quantities() 2 0.02814 0.042 0.04311 0.31% AmrLevel::derive() 18 0.03853 0.039 0.0411 0.29% MLMG::addInterpCorrection() 768 0.03135 0.03378 0.04071 0.29% StateData::FillBoundary(geom) 193 0 0.01815 0.03573 0.25% amrex::average_down 768 0.02425 0.0263 0.03284 0.23% Castro::reset_internal_energy(MultiFab) 80 0.02439 0.02472 0.02531 0.18% VisMF::Write(FabArray) 4 0.01351 0.01465 0.02495 0.18% FabArray::FillBoundary() 6380 0.0167 0.01719 0.02019 0.14% MLCellLinOp::defineAuxData() 19 0.01576 0.01676 0.02011 0.14% Castro::enforce_min_density() 79 0.0188 0.01898 0.01919 0.14% Castro::construct_old_source() 240 0.01885 0.01899 0.01918 0.14% FabArrayBase::getFB() 6380 0.01343 0.01391 0.01669 0.12% MLCGSolver::bicgstab 8 0 0.001042 0.01667 0.12% Castro::evaluate_source_change() 196 0.01488 0.01521 0.0166 0.12% Amr::writePlotFile() 2 0.006721 0.01535 0.01642 0.12% FabArray::setDomainBndry() 187 0.0006653 0.008563 0.01603 0.11% FillPatchIterator::Initialize 187 0.01517 0.01543 0.01575 0.11% Castro::construct_new_source() 102 0.01477 0.01486 0.01523 0.11% FabArray::norminf() 244 0.01019 0.01061 0.01406 0.10% amrex::Add() 128 0.01351 0.01359 0.01379 0.10% MLCellLinOp::smooth() 1544 0.01074 0.01119 0.01351 0.10% MLMG::prepareForSolve() 18 0.01008 0.01044 0.01149 0.08% Castro::wd_update() 2 0.004102 0.008038 0.01065 0.08% FabArrayBase::getCPC() 2045 0.008771 0.008929 0.009461 0.07% Gravity::actual_solve_with_mlmg() 19 0.007307 0.007765 0.00934 0.07% MLMG::solve() 18 0.0007957 0.003447 0.008604 0.06% MLMG::ResNormInf() 146 0.002732 0.004826 0.006533 0.05% Gravity::get_old_grav_vector() 20 0.006053 0.006169 0.00627 0.04% Gravity::get_new_grav_vector() 18 0.005606 0.005696 0.005829 0.04% MLCellLinOp::apply() 986 0.002403 0.002666 0.00573 0.04% Castro::initialize_do_advance() 20 0.00466 0.004889 0.005107 0.04% FabArray::ParallelCopy() 1469 0.004192 0.004327 0.004768 0.03% AmrLevel::FillPatch() 187 0.003074 0.003227 0.003705 0.03% MLLinOp::defineGrids() 19 0.002533 0.002705 0.003101 0.02% MLCellLinOp::correctionResidual() 768 0.00225 0.002316 0.002438 0.02% MLMG::mgVcycle() 128 0.002158 0.002231 0.002357 0.02% MLMG:computeResOfCorrection() 768 0.001808 0.00186 0.002055 0.01% StateDataPhysBCFunct::() 187 0.0009872 0.001174 0.001416 0.01% Castro::save_data_for_retry() 3 0.001287 0.001311 0.001345 0.01% MLMG::actualBottomSolve() 128 0.0001844 0.0002553 0.001251 0.01% Castro::subcycle_advance_ctu() 1 0.000242 0.000328 0.001096 0.01% FillPatchIterator::FillFromLevel0() 187 0.0008731 0.0009011 0.0009674 0.01% main() 1 0.0007578 0.0008787 0.0009659 0.01% Gravity::solve_for_phi() 17 0.0006569 0.0006899 0.0007709 0.01% AmrMesh::MakeNewGrids() 2 0.0005746 0.0006076 0.0007255 0.01% MLMG::mgVcycle_down::2 128 0.0005792 0.0006134 0.0007014 0.00% MLMG::mgVcycle_down::0 128 0.0006346 0.0006484 0.0006929 0.00% MLMG::mgVcycle_down::1 128 0.0006041 0.0006181 0.0006895 0.00% FillPatchSingleLevel 187 0.0006388 0.0006622 0.0006818 0.00% MLMG::mgVcycle_down::5 128 0.0005811 0.0005956 0.0006584 0.00% MLMG::mgVcycle_down::3 128 0.0005571 0.0005856 0.0006485 0.00% MLMG::mgVcycle_down::4 128 0.0005561 0.0005706 0.0006281 0.00% MLCellLinOp::solutionResidual() 147 0.0005153 0.0005265 0.0005531 0.00% Castro::clean_state() 79 0.0004813 0.0005039 0.0005474 0.00% Castro::do_old_sources() 40 0.0004539 0.0004692 0.000531 0.00% MLMG::oneIter() 128 0.0004368 0.0004487 0.0005027 0.00% MLMG::mgVcycle_up::2 128 0.0003926 0.0004099 0.0004722 0.00% Castro::Castro() 1 0.0003906 0.00043 0.0004706 0.00% MLMG::mgVcycle_up::0 128 0.0004189 0.0004304 0.0004688 0.00% MLMG::mgVcycle_up::1 128 0.0004075 0.0004227 0.0004658 0.00% MLMG::mgVcycle_up::4 128 0.0003872 0.0003947 0.0004467 0.00% MLMG::mgVcycle_up::3 128 0.0003801 0.0003914 0.0004441 0.00% MLMG::mgVcycle_up::5 128 0.0003744 0.0003824 0.000404 0.00% MLMG::computeResidual() 128 0.0002919 0.0002986 0.0003173 0.00% MLMG::mgVcycle_bottom 128 0.0002705 0.0002793 0.0003034 0.00% Castro::construct_new_gravity() 17 0.0001683 0.0001901 0.0002922 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002123 0.0002326 0.0002614 0.00% Castro::problem_post_timestep() 2 1.172e-05 2.901e-05 0.00026 0.00% Castro::do_new_sources() 17 0.0002145 0.0002228 0.0002433 0.00% Amr::coarseTimeStep() 1 0.0001426 0.0001671 0.0001704 0.00% Amr::defBaseLevel() 1 0.0001223 0.0001479 0.0001663 0.00% Castro::check_for_nan() 37 0.0001201 0.0001323 0.0001512 0.00% MLPoisson::define() 19 0.0001288 0.0001367 0.000144 0.00% Castro::initialize_advance() 1 0.0001086 0.0001256 0.0001403 0.00% Castro::construct_old_gravity() 20 0.0001132 0.0001255 0.000138 0.00% TagBoxArray::mapPRD 2 0.0001195 0.0001258 0.0001325 0.00% Castro::retry_advance_ctu() 11 7.192e-05 8.165e-05 0.0001274 0.00% Castro::apply_source_to_state() 37 0.0001139 0.0001191 0.000127 0.00% Gravity::test_composite_phi() 1 0.000106 0.0001111 0.0001233 0.00% MLLinOp::define() 19 8.64e-05 9.056e-05 9.748e-05 0.00% Amr::regrid() 1 1.609e-05 2.329e-05 7.798e-05 0.00% MLMG::computeMLResidual() 18 5.822e-05 6.157e-05 7.757e-05 0.00% Amr::FinalizeInit() 1 8.485e-06 1.431e-05 7.232e-05 0.00% Amr::grid_places() 2 5.393e-05 5.939e-05 6.879e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.8e-05 6.108e-05 6.517e-05 0.00% Castro::expand_state() 20 5.452e-05 5.925e-05 6.419e-05 0.00% Castro::create_source_corrector() 17 5.347e-05 5.789e-05 6.394e-05 0.00% Castro::derive() 18 4.708e-05 4.981e-05 5.567e-05 0.00% Gravity::actual_multilevel_solve() 1 3.178e-05 3.31e-05 3.603e-05 0.00% Castro::post_init() 1 1.236e-05 1.599e-05 3.41e-05 0.00% Amr::timeStep() 1 1.547e-05 1.782e-05 2.432e-05 0.00% Castro::post_regrid() 2 1.133e-05 1.424e-05 2.271e-05 0.00% Amr::InitializeInit() 1 1.054e-05 1.425e-05 1.729e-05 0.00% Castro::advance() 1 8.21e-06 9.337e-06 1.252e-05 0.00% Amr::bldFineLevels() 1 7.737e-06 9.076e-06 1.134e-05 0.00% Amr::init() 1 6.811e-06 7.779e-06 8.376e-06 0.00% Castro::post_timestep() 1 6.264e-06 6.702e-06 7.649e-06 0.00% Amr::initialInit() 1 3.097e-06 3.581e-06 4.147e-06 0.00% Other 2530 0.04794 0.06065 0.1 0.71% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.09 14.09 14.1 100.00% Amr::coarseTimeStep() 1 11.13 11.13 11.13 78.92% Amr::timeStep() 1 11.08 11.08 11.08 78.55% Castro::advance() 1 11.02 11.02 11.03 78.20% Castro::subcycle_advance_ctu() 1 10.98 10.98 10.99 77.92% Castro::do_advance_ctu() 20 10.95 10.95 10.95 77.66% Castro::construct_ctu_hydro_source() 20 5.861 5.88 5.91 41.92% Castro::clean_state() 79 2.885 2.968 3.085 21.88% Castro::computeTemp() 80 2.006 2.091 2.211 15.68% Amr::InitAmr() 1 2.041 2.043 2.048 14.53% Gravity::solve_phi_with_mlmg() 19 1.551 1.592 1.609 11.41% Castro::construct_new_gravity() 17 1.471 1.511 1.528 10.84% Gravity::solve_for_phi() 17 1.436 1.477 1.494 10.60% Gravity::actual_solve_with_mlmg() 19 1.361 1.364 1.368 9.71% MLMG::solve() 18 1.306 1.311 1.313 9.31% Castro::reset_internal_energy(MultiFab) 80 1.295 1.3 1.31 9.29% MLMG::oneIter() 128 1.18 1.183 1.185 8.41% MLMG::mgVcycle() 128 1.167 1.169 1.171 8.31% Castro::initialize_do_advance() 20 1.093 1.121 1.143 8.11% Castro::do_old_sources() 40 1.08 1.099 1.12 7.94% Castro::evaluate_source_change() 196 0.9385 0.9419 0.9534 6.76% FabArray::ParallelCopy() 1469 0.4753 0.7794 0.95 6.74% Amr::init() 1 0.9226 0.9226 0.9226 6.54% Amr::initialInit() 1 0.8714 0.8714 0.8715 6.18% Castro::do_new_sources() 17 0.7616 0.7976 0.829 5.88% MLCellLinOp::smooth() 1544 0.6062 0.6603 0.7945 5.64% Castro::enforce_min_density() 79 0.7774 0.7808 0.7902 5.60% amrex::Copy() 766 0.7591 0.7637 0.773 5.48% FabArray::ParallelCopy_finish() 1469 0.2126 0.5135 0.6874 4.88% AmrLevel::FillPatch() 187 0.5489 0.607 0.6603 4.68% Amr::InitializeInit() 1 0.4545 0.4818 0.6172 4.38% Amr::defBaseLevel() 1 0.4545 0.4817 0.6172 4.38% MultiFab::AddProduct() 3598 0.5401 0.5427 0.5514 3.91% FillPatchIterator::Initialize 187 0.3609 0.4182 0.4682 3.32% MLMG::mgVcycle_up::1 128 0.104 0.2787 0.4518 3.20% MLCellLinOp::applyBC() 7148 0.3256 0.3612 0.4472 3.17% FillPatchIterator::FillFromLevel0() 187 0.3296 0.3942 0.4417 3.13% FillPatchSingleLevel 187 0.3287 0.3933 0.4408 3.13% MLPoisson::Fsmooth() 6144 0.3303 0.3565 0.4283 3.04% Amr::FinalizeInit() 1 0.2542 0.3896 0.4168 2.96% Castro::Castro() 1 0.3876 0.3894 0.3951 2.80% MLMG::addInterpCorrection() 768 0.06237 0.3176 0.3909 2.77% Castro::initMFs() 1 0.3825 0.3844 0.3899 2.77% FabArray::setVal() 5232 0.3614 0.3643 0.3697 2.62% Castro::reset_internal_energy(Fab) 160 0.343 0.3445 0.3488 2.47% MultiFab::Subtract() 159 0.3385 0.3397 0.3426 2.43% FabArray::FillBoundary() 6380 0.2469 0.277 0.3173 2.25% MLMG::mgVcycle_up::2 128 0.0156 0.1265 0.2972 2.11% FabArray::ParallelCopy_nowait() 1469 0.2407 0.2615 0.2831 2.01% MLMG::mgVcycle_down::0 128 0.2585 0.2598 0.2616 1.86% Gravity::fill_multipole_BCs() 19 0.1872 0.2259 0.2414 1.71% Castro::expand_state() 20 0.1817 0.2056 0.2284 1.62% Castro::initData() 1 0.0652 0.09056 0.2277 1.62% MLCellLinOp::apply() 986 0.1496 0.1597 0.2022 1.43% Amr::bldFineLevels() 1 0.03721 0.1733 0.1998 1.42% MLMG::mgVcycle_up::0 128 0.1975 0.1988 0.1996 1.42% MLMG::mgVcycle_up::3 128 0.01546 0.02926 0.1956 1.39% FabArray::sum() 3598 0.1828 0.1844 0.1865 1.32% Castro::post_init() 1 0.1778 0.1805 0.1815 1.29% FillBoundary_nowait() 6380 0.1021 0.1331 0.167 1.18% Amr::grid_places() 2 0.00261 0.1381 0.1653 1.17% AmrMesh::MakeNewGrids() 2 0.002556 0.138 0.1652 1.17% FillBoundary_finish() 6380 0.1061 0.1267 0.1637 1.16% TagBoxArray::mapPRD 2 0.0009243 0.1084 0.1635 1.16% Castro::construct_old_source() 240 0.1564 0.1575 0.1589 1.13% TagBoxArray::collate() 2 0.0001904 0.02818 0.1586 1.13% MLMG:computeResOfCorrection() 768 0.1139 0.1235 0.1472 1.04% check_for_negative_density() 20 0.06173 0.123 0.1464 1.04% MLCellLinOp::correctionResidual() 768 0.1121 0.1217 0.1453 1.03% MLMG::mgVcycle_down::1 128 0.1212 0.1248 0.1277 0.91% Castro::normalize_species() 79 0.119 0.1195 0.1205 0.85% Castro::construct_new_source() 102 0.1187 0.1192 0.1203 0.85% MLPoisson::Fapply() 986 0.0867 0.09048 0.1052 0.75% Gravity::multilevel_solve_for_new_phi() 1 0.1046 0.1047 0.1047 0.74% Gravity::actual_multilevel_solve() 1 0.1045 0.1046 0.1047 0.74% Castro::construct_old_gravity_source() 40 0.09429 0.09506 0.09616 0.68% MLCellLinOp::solutionResidual() 147 0.08366 0.08905 0.09463 0.67% Amr::writePlotFile() 2 0.09328 0.09338 0.0934 0.66% MLMG::mgVcycle_down::2 128 0.01812 0.0481 0.08178 0.58% MLMG::mgVcycle_bottom 128 0.0004569 0.005516 0.08125 0.58% MLMG::actualBottomSolve() 128 0.0001844 0.005236 0.08095 0.57% FabArray::Saxpy() 111 0.07934 0.07979 0.08066 0.57% MLCGSolver::bicgstab 8 0 0.004935 0.07896 0.56% Castro::derive() 18 0.0658 0.06697 0.07809 0.55% AmrLevel::derive() 18 0.06575 0.06692 0.07803 0.55% Castro::post_regrid() 2 0.06998 0.07105 0.07368 0.52% MLMG::computeResidual() 128 0.07158 0.07282 0.07344 0.52% Castro::construct_new_gravity_source() 17 0.06664 0.06724 0.06775 0.48% MLMG::mgVcycle_down::3 128 0.01754 0.02377 0.06471 0.46% Castro::sum_integrated_quantities() 2 0.06198 0.06214 0.06215 0.44% Castro::create_source_corrector() 17 0.04704 0.05432 0.06193 0.44% Castro::construct_old_gravity() 20 0.04187 0.05135 0.06133 0.43% Gravity::get_old_grav_vector() 20 0.04175 0.05123 0.0612 0.43% amrex::average_down 768 0.03409 0.04053 0.05775 0.41% FabArray::Xpay() 915 0.04829 0.05011 0.05527 0.39% Castro::post_timestep() 1 0.04953 0.05132 0.05333 0.38% Castro::check_for_nan() 37 0.04554 0.04585 0.04634 0.33% MultiFab::contains_nan() 37 0.04542 0.04572 0.04621 0.33% MLPoisson::define() 19 0.03133 0.03644 0.04186 0.30% Castro::initialize_advance() 1 0.03995 0.04034 0.04113 0.29% MLMG::mgVcycle_down::4 128 0.01759 0.01932 0.04078 0.29% Amr::regrid() 1 0.03579 0.03658 0.03811 0.27% StateDataPhysBCFunct::() 187 0.0009872 0.01933 0.03701 0.26% StateData::FillBoundary(geom) 193 0 0.01815 0.03573 0.25% Castro::retry_advance_ctu() 11 0.03394 0.0344 0.03527 0.25% Gravity::get_new_grav_vector() 18 0.03351 0.03406 0.03452 0.24% MLMG::mgVcycle_up::4 128 0.01555 0.01683 0.03234 0.23% MLMG::mgVcycle_down::5 128 0.0176 0.01864 0.02982 0.21% MLCellLinOp::defineAuxData() 19 0.02036 0.02218 0.02744 0.19% Castro::apply_source_to_state() 37 0.02655 0.02676 0.02717 0.19% MLMG::mgVcycle_up::5 128 0.01563 0.01646 0.02557 0.18% VisMF::Write(FabArray) 4 0.01351 0.01465 0.02495 0.18% MLMG::prepareForSolve() 18 0.02015 0.02109 0.02385 0.17% MLMG::computeMLResidual() 18 0.01031 0.01584 0.0208 0.15% FabArrayBase::getFB() 6380 0.01378 0.0146 0.01817 0.13% Castro::save_data_for_retry() 3 0.01731 0.01749 0.01798 0.13% FabArray::setDomainBndry() 187 0.0006653 0.008563 0.01603 0.11% FabArrayBase::getCPC() 2045 0.01109 0.01241 0.01578 0.11% MLMG::ResNormInf() 146 0.01168 0.01389 0.01563 0.11% FabArray::norminf() 244 0.01019 0.01061 0.01406 0.10% Castro::problem_post_timestep() 2 0.01003 0.01179 0.01398 0.10% Castro::wd_update() 2 0.01002 0.01176 0.01396 0.10% amrex::Add() 128 0.01351 0.01359 0.01379 0.10% Gravity::test_composite_phi() 1 0.01353 0.01358 0.01361 0.10% MLLinOp::define() 19 0.00433 0.007645 0.01105 0.08% MLLinOp::defineGrids() 19 0.004242 0.007554 0.01096 0.08% Other 2530 0.07842 0.09353 0.139 0.99% -------------------------------------------------------------------------------------------- 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 383 KiB 385 KiB 390 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 2057 B 2080 B 2093 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 2 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 1489 KiB 1490 KiB 1490 KiB 2329 KiB 2329 KiB 2329 KiB Castro::construct_ctu_hydro_source() 6400 383 KiB 385 KiB 390 KiB 1176 KiB 1176 KiB 1176 KiB StateData::define() 192 994 KiB 994 KiB 995 KiB 1164 KiB 1164 KiB 1164 KiB Castro::initMFs() 256 974 KiB 975 KiB 975 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 633 KiB 633 KiB 633 KiB 803 KiB 803 KiB 803 KiB Amr::writePlotFile() 64 4271 B 4331 B 4839 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 49 KiB 49 KiB 49 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 12 KiB 13 KiB 14 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 408 KiB 408 KiB 408 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2130 B 3806 B 4842 B 225 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1357 B 1370 B 1428 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 123 KiB 123 KiB 123 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 121 KiB 121 KiB 121 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 9723 B 10 KiB 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 1280 490 B 492 B 498 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 414 B 415 B 421 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 2057 B 2080 B 2093 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 18 B 21 B 38 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 31 KiB 31 KiB 31 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 6 B 8 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 107 B 108 B 110 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 134 B 165 B 198 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 29 KiB 29 KiB 29 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 226 B 231 B 235 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 1084 B 1088 B 1101 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 121 B 121 B 121 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 1667 B 1713 B 1734 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 440 B 491 B 682 B 4608 B 5151 B 7152 B BndryData::define() 4864 457 B 458 B 459 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 24 B 33 B 44 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 0 B 19 B 30 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 16 B 26 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 3 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 0 B 0 B 0 B 8 B 8 B 8 B ----------------------------------------------------------------------------------------------------------------------- AMReX (26.03-110-g388b0db16adb) finalized