Initializing AMReX (26.07-55-gd1b9e4ec5555)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.07-55-gd1b9e4ec5555) initialized Starting run at 06:51:55 UTC on 2026-07-31. Successfully read inputs file ... Castro git describe: 26.07-10-gd4842676d AMReX git describe: 26.07-55-gd1b9e4ec5 Microphysics git describe: 26.07-4-gdff9d5a4 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.206704201 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.00119415 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010272567 Gravity::multilevel_solve_for_new_phi() time = 0.108289877 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010290416 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.047093277 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004127973 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, Invalid density = -0.01295059363 at index 21, 172, 0 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 Invalid density = -8.076810125 at index 20, 174, 0 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -625.0368072 at index 5, 65, 0Invalid density = -108.2035502 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 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.4708539Invalid density = -80.16477044 at index 18, 179, 0 at index 9, 66, 0 Invalid density = -6.239071148 at index 10, 66, 0 Invalid density = -160.2179808 at index 18, 180Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = , 0 Invalid density = -552.769409 at index -241.113047 at index 11, 67, 0 17, 181, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -415.2911996 at index 16, 183Invalid density = -859.0333271 at index 13, 69, 0 , 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -415.2911996 at index 16, 72, 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 = -552.769409 at index 17, 74, 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 = -160.2179808 at index 18, 75, 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, Invalid density = -80.16477044 at index 18, 76, 0 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, Invalid density = -61.58489855 at index 19, 77, 0 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, 191Invalid density = -108.2035502 at index 19, 78, 0 , 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.002728377 on level 0 Castro::do_old_sources() time = 0.007389143 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.054501095 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.003066555 on level 0 Castro::do_old_sources() time = 0.007372115 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.052911443 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.005998644 on level 0 Castro::do_old_sources() time = 0.009265338 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.049444602 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.004549968 on level 0 Castro::do_old_sources() time = 0.007865378 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.04899505 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.292583565 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.013169721 Gravity::solve_for_phi() time = 0.089820884 on level 0 Castro::construct_new_gravity() time = 0.091912966 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.049051699 on level 0 Castro::construct_old_gravity() time = 0.002639003 on level 0 Castro::do_old_sources() time = 0.007308905 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.04867448 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.290751169 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.01277616 Gravity::solve_for_phi() time = 0.090494898 on level 0 Castro::construct_new_gravity() time = 0.09250913 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.048721927 on level 0 Castro::construct_old_gravity() time = 0.002683487 on level 0 Castro::do_old_sources() time = 0.007385361 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.048721955 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.296317678 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.012751068 Gravity::solve_for_phi() time = 0.089445004 on level 0 Castro::construct_new_gravity() time = 0.091448068 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.048605551 on level 0 Castro::construct_old_gravity() time = 0.00259029 on level 0 Castro::do_old_sources() time = 0.007224177 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.048694958 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.302901603 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.013288189 Gravity::solve_for_phi() time = 0.091162581 on level 0 Castro::construct_new_gravity() time = 0.093204854 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.048717273 on level 0 Castro::construct_old_gravity() time = 0.002648549 on level 0 Castro::do_old_sources() time = 0.007270983 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.0503253 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.345919939 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.040382771 Gravity::solve_for_phi() time = 0.120537969 on level 0 Castro::construct_new_gravity() time = 0.122617831 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.049748509 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.00359771 on level 0 Castro::do_old_sources() time = 0.007361659 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.058238818 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.340390461 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.017033018 Gravity::solve_for_phi() time = 0.095197945 on level 0 Castro::construct_new_gravity() time = 0.097231001 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.048607011 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.010448246 on level 0 Castro::do_old_sources() time = 0.007935789 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.048786238 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.307574192 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.013031436 Gravity::solve_for_phi() time = 0.091096612 on level 0 Castro::construct_new_gravity() time = 0.093088974 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.048878689 on level 0 Castro::construct_old_gravity() time = 0.00283578 on level 0 Castro::do_old_sources() time = 0.007333542 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.049956273 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.304384603 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.013137626 Gravity::solve_for_phi() time = 0.090564782 on level 0 Castro::construct_new_gravity() time = 0.092616993 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.049125831 on level 0 Castro::construct_old_gravity() time = 0.002758388 on level 0 Castro::do_old_sources() time = 0.00734019 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.048495459 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.30750108 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.013187957 Gravity::solve_for_phi() time = 0.090135252 on level 0 Castro::construct_new_gravity() time = 0.092173451 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.04891768 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.002679034 on level 0 Castro::do_old_sources() time = 0.007589075 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.048569842 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.302570139 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.013383662 Gravity::solve_for_phi() time = 0.089797911 on level 0 Castro::construct_new_gravity() time = 0.091803089 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.049120885 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.002576727 on level 0 Castro::do_old_sources() time = 0.008845458 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.048404572 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.306926329 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.013416222 Gravity::solve_for_phi() time = 0.08962533 on level 0 Castro::construct_new_gravity() time = 0.091646119 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.048975454 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.002605581 on level 0 Castro::do_old_sources() time = 0.007408807 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.048338302 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.302761246 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.013515033 Gravity::solve_for_phi() time = 0.090236264 on level 0 Castro::construct_new_gravity() time = 0.092310374 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.049158702 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.002631074 on level 0 Castro::do_old_sources() time = 0.007254235 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.048360058 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.30828565 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.013368791 Gravity::solve_for_phi() time = 0.090188575 on level 0 Castro::construct_new_gravity() time = 0.092242072 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.048828723 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.002641134 on level 0 Castro::do_old_sources() time = 0.007261172 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.0485026 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.304694776 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.013624928 Gravity::solve_for_phi() time = 0.090187491 on level 0 Castro::construct_new_gravity() time = 0.09221495 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.048905895 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.002658374 on level 0 Castro::do_old_sources() time = 0.007381707 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.048367046 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.310145928 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.01356882 Gravity::solve_for_phi() time = 0.090519714 on level 0 Castro::construct_new_gravity() time = 0.09258288 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.048806251 on level 0 Castro::construct_old_gravity() time = 0.00270285 on level 0 Castro::do_old_sources() time = 0.007338933 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.048487644 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.303379377 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.013972999 Gravity::solve_for_phi() time = 0.091346053 on level 0 Castro::construct_new_gravity() time = 0.093361674 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.049104842 on level 0 Castro::construct_old_gravity() time = 0.002635552 on level 0 Castro::do_old_sources() time = 0.00734045 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.048990602 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.310565022 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002926730302 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.013647242 Gravity::solve_for_phi() time = 0.090738685 on level 0 Castro::construct_new_gravity() time = 0.09275788 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.049075931 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.24598999 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.046366058 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004157322 seconds Ending run at 06:52:08 UTC on 2026-07-31. Run time = 13.51011772 Run time without initialization = 11.29742998 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]: 13.5018 ... 13.5042 ... 13.5102 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.1303 5.1673 5.2091 38.56% Amr::InitAmr() 1 1.267 1.27 1.276 9.45% Castro::computeTemp() 80 0.705 0.7936 0.9211 6.82% amrex::Copy() 766 0.7794 0.7836 0.8039 5.95% FabArray::ParallelCopy_finish() 1469 0.2479 0.5517 0.7454 5.52% MultiFab::AddProduct() 3598 0.5497 0.5513 0.5546 4.11% MLPoisson::Fsmooth() 6144 0.3289 0.3583 0.427 3.16% Castro::initMFs() 1 0.389 0.395 0.4225 3.13% FabArray::setVal() 5232 0.3571 0.3598 0.3644 2.70% Castro::reset_internal_energy(Fab) 160 0.3474 0.3501 0.3571 2.64% MultiFab::Subtract() 159 0.3275 0.3294 0.3392 2.51% Gravity::fill_multipole_BCs() 19 0.2157 0.2564 0.2743 2.03% FabArray::ParallelCopy_nowait() 1469 0.2306 0.249 0.2702 2.00% Castro::do_advance_ctu() 20 0.2603 0.2632 0.2675 1.98% FillBoundary_finish() 6380 0.109 0.1599 0.2006 1.48% FabArray::sum() 3598 0.186 0.1867 0.1885 1.40% TagBoxArray::collate() 2 0.0002303 0.04521 0.1865 1.38% check_for_negative_density() 20 0.0683 0.1415 0.186 1.38% Castro::initData() 1 0.003153 0.02919 0.1736 1.28% FillBoundary_nowait() 6380 0.09018 0.1192 0.1532 1.13% MLCellLinOp::applyBC() 7148 0.07045 0.08678 0.1322 0.98% Castro::normalize_species() 79 0.1197 0.1208 0.1294 0.96% MLPoisson::Fapply() 986 0.08603 0.09052 0.1042 0.77% Castro::construct_old_gravity_source() 40 0.09219 0.09287 0.09373 0.69% FabArray::Saxpy() 111 0.07654 0.07689 0.07736 0.57% Castro::construct_new_gravity_source() 17 0.0668 0.06752 0.06917 0.51% FabArray::Xpay() 915 0.04901 0.0511 0.05621 0.42% MultiFab::contains_nan() 37 0.04441 0.04465 0.0451 0.33% MLMG::addInterpCorrection() 768 0.03249 0.03488 0.04215 0.31% AmrLevel::derive() 18 0.03893 0.03932 0.04024 0.30% StateData::FillBoundary(geom) 193 0 0.01829 0.03544 0.26% amrex::average_down 768 0.02404 0.02621 0.03273 0.24% Castro::reset_internal_energy(MultiFab) 80 0.02507 0.02555 0.02593 0.19% VisMF::Write(FabArray) 4 0.01266 0.01399 0.02398 0.18% MLCellLinOp::defineAuxData() 19 0.01593 0.01696 0.02051 0.15% Castro::enforce_min_density() 79 0.01921 0.0196 0.01997 0.15% FabArray::FillBoundary() 6380 0.01646 0.01693 0.01996 0.15% Castro::construct_old_source() 240 0.01897 0.01921 0.0195 0.14% Amr::writePlotFile() 2 0.00659 0.01538 0.01683 0.12% FabArray::setDomainBndry() 187 0.0006792 0.008866 0.01672 0.12% FillPatchIterator::Initialize 187 0.01576 0.01614 0.01666 0.12% MLCGSolver::bicgstab 8 0 0.001038 0.01661 0.12% FabArrayBase::getFB() 6380 0.01319 0.01369 0.01652 0.12% Castro::evaluate_source_change() 196 0.0144 0.01491 0.01624 0.12% Castro::construct_new_source() 102 0.01482 0.01497 0.01537 0.11% FabArray::norminf() 244 0.01058 0.01102 0.01443 0.11% MLMG::ResNormInf() 146 0.003657 0.008158 0.01419 0.11% amrex::Add() 128 0.01331 0.01345 0.0139 0.10% MLCellLinOp::smooth() 1544 0.01025 0.01063 0.01234 0.09% MLMG::prepareForSolve() 18 0.01012 0.0105 0.01141 0.08% Castro::wd_update() 2 0.00468 0.007435 0.01028 0.08% Gravity::actual_solve_with_mlmg() 19 0.007391 0.007847 0.0094 0.07% FabArrayBase::getCPC() 2045 0.008631 0.008815 0.009282 0.07% MLMG::solve() 18 0.001006 0.003879 0.007937 0.06% Gravity::get_old_grav_vector() 20 0.006039 0.006096 0.006215 0.05% Gravity::get_new_grav_vector() 18 0.005603 0.005697 0.005832 0.04% MLCellLinOp::apply() 986 0.002361 0.002602 0.005539 0.04% Castro::initialize_do_advance() 20 0.004708 0.004859 0.005061 0.04% FabArray::ParallelCopy() 1469 0.004131 0.004256 0.004577 0.03% AmrLevel::FillPatch() 187 0.00312 0.003237 0.003612 0.03% MLLinOp::defineGrids() 19 0.002611 0.002756 0.0031 0.02% MLCellLinOp::correctionResidual() 768 0.002247 0.002292 0.002362 0.02% MLMG::mgVcycle() 128 0.002152 0.002219 0.002308 0.02% Castro::sum_integrated_quantities() 2 0.001587 0.001826 0.002115 0.02% MLMG:computeResOfCorrection() 768 0.001796 0.001839 0.001985 0.01% Castro::save_data_for_retry() 3 0.001319 0.001343 0.001394 0.01% StateDataPhysBCFunct::() 187 0.001031 0.001202 0.001335 0.01% MLMG::actualBottomSolve() 128 0.0001844 0.0002523 0.001233 0.01% main() 1 0.0007863 0.0008611 0.001212 0.01% FillPatchIterator::FillFromLevel0() 187 0.0008988 0.0009303 0.001014 0.01% Castro::subcycle_advance_ctu() 1 0.0002382 0.0003034 0.0009228 0.01% Gravity::solve_for_phi() 17 0.0006941 0.0007308 0.0007808 0.01% FillPatchSingleLevel 187 0.0006455 0.000671 0.0007152 0.01% MLMG::mgVcycle_down::1 128 0.0006139 0.0006268 0.0006845 0.01% MLMG::mgVcycle_down::0 128 0.0006339 0.0006481 0.0006816 0.01% MLMG::mgVcycle_down::2 128 0.0005915 0.0006245 0.0006687 0.00% MLMG::mgVcycle_down::5 128 0.0005759 0.0005957 0.0006513 0.00% MLMG::mgVcycle_down::3 128 0.00056 0.00059 0.0006419 0.00% Castro::Castro() 1 0.0004794 0.0005459 0.0006412 0.00% MLMG::mgVcycle_down::4 128 0.0005498 0.0005745 0.0006359 0.00% AmrMesh::MakeNewGrids() 2 0.0005589 0.0005914 0.0006321 0.00% Castro::clean_state() 79 0.0004573 0.0004985 0.0005426 0.00% MLCellLinOp::solutionResidual() 147 0.0005184 0.0005301 0.0005423 0.00% Castro::do_old_sources() 40 0.0004447 0.0004634 0.0005221 0.00% MLMG::oneIter() 128 0.0004279 0.0004396 0.0004765 0.00% MLMG::mgVcycle_up::1 128 0.0004065 0.0004184 0.000453 0.00% MLMG::mgVcycle_up::0 128 0.0004013 0.0004116 0.0004477 0.00% MLMG::mgVcycle_up::2 128 0.0003913 0.0004093 0.000442 0.00% MLMG::mgVcycle_up::3 128 0.0003782 0.0003912 0.0004336 0.00% MLMG::mgVcycle_up::4 128 0.0003808 0.0003917 0.0004311 0.00% MLMG::mgVcycle_up::5 128 0.000365 0.0003773 0.0003948 0.00% MLMG::computeResidual() 128 0.0002927 0.000299 0.0003161 0.00% MLMG::mgVcycle_bottom 128 0.0002644 0.0002755 0.0003047 0.00% Castro::construct_new_gravity() 17 0.0001899 0.0002079 0.0002982 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002214 0.0002414 0.0002673 0.00% Castro::problem_post_timestep() 2 1.136e-05 2.837e-05 0.0002516 0.00% Castro::do_new_sources() 17 0.0002093 0.0002187 0.0002433 0.00% Amr::defBaseLevel() 1 0.0001184 0.0001428 0.0001669 0.00% Amr::coarseTimeStep() 1 0.0001373 0.0001626 0.0001666 0.00% MLPoisson::define() 19 0.0001272 0.0001362 0.0001568 0.00% Castro::construct_old_gravity() 20 0.0001224 0.0001346 0.0001524 0.00% Castro::check_for_nan() 37 0.0001197 0.0001294 0.0001486 0.00% Castro::initialize_advance() 1 0.0001094 0.0001182 0.0001452 0.00% TagBoxArray::mapPRD 2 0.0001194 0.0001267 0.0001381 0.00% Castro::apply_source_to_state() 37 0.0001167 0.0001236 0.0001343 0.00% Gravity::test_composite_phi() 1 0.0001071 0.0001138 0.0001337 0.00% Castro::retry_advance_ctu() 11 6.825e-05 7.536e-05 0.0001274 0.00% Amr::grid_places() 2 5.996e-05 7.347e-05 0.0001195 0.00% MLLinOp::define() 19 8.565e-05 8.961e-05 9.881e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.878e-05 6.256e-05 7.38e-05 0.00% Castro::expand_state() 20 5.825e-05 6.202e-05 7.246e-05 0.00% Castro::create_source_corrector() 17 5.587e-05 6.018e-05 6.942e-05 0.00% Amr::FinalizeInit() 1 1.4e-05 2.336e-05 6.587e-05 0.00% MLMG::computeMLResidual() 18 5.186e-05 5.475e-05 6.468e-05 0.00% Castro::derive() 18 4.652e-05 4.928e-05 5.696e-05 0.00% Amr::regrid() 1 1.589e-05 1.968e-05 4.296e-05 0.00% Gravity::actual_multilevel_solve() 1 3.237e-05 3.365e-05 3.874e-05 0.00% Amr::timeStep() 1 1.546e-05 1.828e-05 2.778e-05 0.00% Castro::post_init() 1 1.528e-05 1.747e-05 2.326e-05 0.00% Castro::advance() 1 8.303e-06 1.074e-05 1.565e-05 0.00% Castro::post_regrid() 2 1.095e-05 1.243e-05 1.475e-05 0.00% Amr::InitializeInit() 1 8.912e-06 1.023e-05 1.433e-05 0.00% Amr::bldFineLevels() 1 7.93e-06 9.259e-06 1.092e-05 0.00% Amr::init() 1 6.429e-06 7.236e-06 9.08e-06 0.00% Castro::post_timestep() 1 5.906e-06 6.371e-06 7.21e-06 0.00% Amr::initialInit() 1 3.216e-06 3.556e-06 4.102e-06 0.00% Other 2530 0.04744 0.05845 0.09517 0.70% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 13.5 13.5 13.51 100.00% Amr::coarseTimeStep() 1 11.3 11.3 11.3 83.62% Amr::timeStep() 1 11.25 11.25 11.25 83.24% Castro::advance() 1 11.19 11.19 11.2 82.88% Castro::subcycle_advance_ctu() 1 11.15 11.15 11.16 82.57% Castro::do_advance_ctu() 20 11.11 11.12 11.12 82.31% Castro::construct_ctu_hydro_source() 20 5.911 5.937 5.977 44.24% Castro::clean_state() 79 2.901 2.988 3.118 23.08% Castro::computeTemp() 80 2.022 2.106 2.236 16.55% Gravity::solve_phi_with_mlmg() 19 1.624 1.669 1.689 12.50% Castro::construct_new_gravity() 17 1.54 1.585 1.604 11.88% Gravity::solve_for_phi() 17 1.505 1.55 1.57 11.62% Gravity::actual_solve_with_mlmg() 19 1.405 1.409 1.413 10.46% MLMG::solve() 18 1.349 1.355 1.356 10.04% Castro::reset_internal_energy(MultiFab) 80 1.307 1.312 1.32 9.77% Amr::InitAmr() 1 1.267 1.27 1.276 9.45% MLMG::oneIter() 128 1.212 1.217 1.225 9.06% MLMG::mgVcycle() 128 1.198 1.203 1.211 8.96% Castro::initialize_do_advance() 20 1.11 1.137 1.159 8.58% Castro::do_old_sources() 40 1.082 1.111 1.138 8.42% FabArray::ParallelCopy() 1469 0.5208 0.8253 1.017 7.53% Castro::evaluate_source_change() 196 0.9462 0.9496 0.9561 7.08% Amr::init() 1 0.9363 0.9363 0.9363 6.93% Amr::initialInit() 1 0.8846 0.8846 0.8847 6.55% Castro::do_new_sources() 17 0.762 0.7997 0.8324 6.16% MLCellLinOp::smooth() 1544 0.6308 0.6879 0.8215 6.08% Castro::enforce_min_density() 79 0.7805 0.785 0.8118 6.01% amrex::Copy() 766 0.7794 0.7836 0.8039 5.95% FabArray::ParallelCopy_finish() 1469 0.2479 0.5517 0.7454 5.52% AmrLevel::FillPatch() 187 0.5802 0.6365 0.6862 5.08% Amr::InitializeInit() 1 0.4644 0.4958 0.6697 4.96% Amr::defBaseLevel() 1 0.4644 0.4958 0.6697 4.96% MultiFab::AddProduct() 3598 0.5497 0.5513 0.5546 4.11% FillPatchIterator::Initialize 187 0.3823 0.4432 0.4929 3.65% MLCellLinOp::applyBC() 7148 0.3369 0.3972 0.4866 3.60% FillPatchIterator::FillFromLevel0() 187 0.3495 0.4182 0.4663 3.45% FillPatchSingleLevel 187 0.3486 0.4173 0.4654 3.44% MLMG::mgVcycle_up::1 128 0.1046 0.2785 0.4538 3.36% Castro::Castro() 1 0.3949 0.401 0.4297 3.18% MLPoisson::Fsmooth() 6144 0.3289 0.3583 0.427 3.16% Castro::initMFs() 1 0.3892 0.3952 0.4226 3.13% Amr::FinalizeInit() 1 0.2149 0.3888 0.4201 3.11% MLMG::addInterpCorrection() 768 0.06481 0.3191 0.3936 2.91% FabArray::setVal() 5232 0.3571 0.3598 0.3644 2.70% Castro::reset_internal_energy(Fab) 160 0.3474 0.3501 0.3571 2.64% FabArray::FillBoundary() 6380 0.2596 0.3104 0.3545 2.62% MultiFab::Subtract() 159 0.3275 0.3294 0.3392 2.51% MLMG::mgVcycle_up::2 128 0.0154 0.1279 0.2989 2.21% FabArray::ParallelCopy_nowait() 1469 0.2505 0.2694 0.2913 2.16% MLMG::mgVcycle_down::0 128 0.271 0.2766 0.2808 2.08% Gravity::fill_multipole_BCs() 19 0.2177 0.2584 0.2762 2.04% Amr::bldFineLevels() 1 0.03833 0.2119 0.2437 1.80% Castro::initData() 1 0.06601 0.09301 0.2382 1.76% Castro::expand_state() 20 0.1874 0.212 0.2327 1.72% MLMG::mgVcycle_up::0 128 0.2079 0.2112 0.215 1.59% MLCellLinOp::apply() 986 0.1569 0.1692 0.2113 1.56% Amr::grid_places() 2 0.002676 0.1766 0.2079 1.54% AmrMesh::MakeNewGrids() 2 0.002556 0.1765 0.2078 1.54% TagBoxArray::mapPRD 2 0.0009145 0.1299 0.2061 1.53% FillBoundary_finish() 6380 0.109 0.1599 0.2006 1.48% MLMG::mgVcycle_up::3 128 0.01536 0.02936 0.1975 1.46% FabArray::sum() 3598 0.186 0.1867 0.1885 1.40% TagBoxArray::collate() 2 0.0002303 0.04521 0.1865 1.38% check_for_negative_density() 20 0.0683 0.1415 0.186 1.38% FillBoundary_nowait() 6380 0.1037 0.1336 0.168 1.24% Castro::construct_old_source() 240 0.1535 0.1545 0.156 1.15% MLMG:computeResOfCorrection() 768 0.1155 0.129 0.1522 1.13% MLCellLinOp::correctionResidual() 768 0.1136 0.1272 0.1503 1.11% Castro::post_init() 1 0.141 0.142 0.143 1.06% MLMG::mgVcycle_down::1 128 0.125 0.13 0.1369 1.01% Castro::normalize_species() 79 0.1197 0.1208 0.1294 0.96% Castro::construct_new_source() 102 0.1175 0.1184 0.12 0.89% Gravity::multilevel_solve_for_new_phi() 1 0.1082 0.1082 0.1083 0.80% Gravity::actual_multilevel_solve() 1 0.1081 0.1082 0.1082 0.80% MLPoisson::Fapply() 986 0.08603 0.09052 0.1042 0.77% MLCellLinOp::solutionResidual() 147 0.08863 0.09393 0.09921 0.73% Amr::writePlotFile() 2 0.09388 0.09397 0.09398 0.70% Castro::construct_old_gravity_source() 40 0.09219 0.09287 0.09373 0.69% MLMG::mgVcycle_bottom 128 0.0004489 0.005539 0.08172 0.60% MLMG::actualBottomSolve() 128 0.0001844 0.005264 0.08141 0.60% MLMG::computeResidual() 128 0.07412 0.07782 0.08064 0.60% MLCGSolver::bicgstab 8 0 0.004965 0.07943 0.59% MLMG::mgVcycle_down::2 128 0.0181 0.04716 0.07902 0.58% FabArray::Saxpy() 111 0.07654 0.07689 0.07736 0.57% Castro::post_regrid() 2 0.06913 0.07017 0.07117 0.53% Castro::construct_new_gravity_source() 17 0.0668 0.06752 0.06917 0.51% Castro::create_source_corrector() 17 0.05145 0.05869 0.06913 0.51% Castro::derive() 18 0.06652 0.06726 0.06833 0.51% AmrLevel::derive() 18 0.06647 0.06721 0.06828 0.51% MLMG::mgVcycle_down::3 128 0.01749 0.02402 0.06814 0.50% Castro::construct_old_gravity() 20 0.04403 0.05188 0.06106 0.45% Gravity::get_old_grav_vector() 20 0.0439 0.05174 0.06091 0.45% amrex::average_down 768 0.03377 0.04061 0.05869 0.43% FabArray::Xpay() 915 0.04901 0.0511 0.05621 0.42% Castro::post_timestep() 1 0.04933 0.05124 0.05366 0.40% Castro::check_for_nan() 37 0.04453 0.04478 0.04525 0.33% MultiFab::contains_nan() 37 0.04441 0.04465 0.0451 0.33% Castro::initialize_advance() 1 0.04035 0.04116 0.04298 0.32% MLPoisson::define() 19 0.03249 0.03748 0.04277 0.32% MLMG::mgVcycle_down::4 128 0.01744 0.01927 0.04115 0.30% Castro::retry_advance_ctu() 11 0.03394 0.03457 0.03833 0.28% Amr::regrid() 1 0.03602 0.03653 0.03702 0.27% StateDataPhysBCFunct::() 187 0.001031 0.01949 0.03671 0.27% StateData::FillBoundary(geom) 193 0 0.01829 0.03544 0.26% Gravity::get_new_grav_vector() 18 0.03412 0.03454 0.03479 0.26% MLMG::mgVcycle_up::4 128 0.01544 0.01674 0.03252 0.24% MLMG::mgVcycle_down::5 128 0.01749 0.01852 0.02986 0.22% MLCellLinOp::defineAuxData() 19 0.02061 0.02235 0.02779 0.21% Castro::apply_source_to_state() 37 0.02562 0.02581 0.02595 0.19% MLMG::mgVcycle_up::5 128 0.01545 0.01637 0.02557 0.19% VisMF::Write(FabArray) 4 0.01266 0.01399 0.02398 0.18% MLMG::prepareForSolve() 18 0.02019 0.02116 0.02377 0.18% MLMG::ResNormInf() 146 0.01308 0.01753 0.02355 0.17% Castro::save_data_for_retry() 3 0.01787 0.01834 0.0219 0.16% Castro::sum_integrated_quantities() 2 0.02114 0.02128 0.02129 0.16% MLMG::computeMLResidual() 18 0.01029 0.01572 0.02073 0.15% FabArrayBase::getFB() 6380 0.01351 0.01438 0.01796 0.13% FabArray::setDomainBndry() 187 0.0006792 0.008866 0.01672 0.12% FabArrayBase::getCPC() 2045 0.01111 0.01231 0.01574 0.12% FabArray::norminf() 244 0.01058 0.01102 0.01443 0.11% Gravity::test_composite_phi() 1 0.01403 0.01409 0.01412 0.10% amrex::Add() 128 0.01331 0.01345 0.0139 0.10% Castro::problem_post_timestep() 2 0.009412 0.01109 0.01363 0.10% Castro::wd_update() 2 0.009399 0.01106 0.01361 0.10% MLLinOp::define() 19 0.00435 0.00792 0.0113 0.08% MLLinOp::defineGrids() 19 0.004261 0.007831 0.0112 0.08% Other 2530 0.07815 0.08997 0.1328 0.98% -------------------------------------------------------------------------------------------- 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 402 KiB 405 KiB 408 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 2178 B 2206 B 2270 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 2 B 3 B 8192 KiB 8192 KiB 8192 KiB -------------------------------------------------------------------------------------------------------------- Cpu Memory Usage: ----------------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max ----------------------------------------------------------------------------------------------------------------------- Castro::save_data_for_retry() 416 1580 KiB 1580 KiB 1581 KiB 2332 KiB 2332 KiB 2332 KiB Castro::construct_ctu_hydro_source() 6400 403 KiB 405 KiB 408 KiB 1178 KiB 1178 KiB 1178 KiB StateData::define() 192 1054 KiB 1055 KiB 1055 KiB 1166 KiB 1166 KiB 1166 KiB Castro::initMFs() 256 1033 KiB 1033 KiB 1033 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 671 KiB 672 KiB 672 KiB 804 KiB 804 KiB 804 KiB Amr::writePlotFile() 64 4471 B 4553 B 5098 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 51 KiB 51 KiB 52 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 13 KiB 14 KiB 15 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 432 KiB 433 KiB 433 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2355 B 4165 B 5248 B 226 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1428 B 1442 B 1457 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 19 KiB 19 KiB 20 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 130 KiB 130 KiB 130 KiB 158 KiB 158 KiB 158 KiB Castro::buildMetrics() 192 129 KiB 129 KiB 129 KiB 143 KiB 143 KiB 143 KiB MLMG::prepareForSolve() 6102 10 KiB 11 KiB 12 KiB 104 KiB 110 KiB 130 KiB Castro::construct_old_source() 1280 503 B 506 B 511 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 421 B 423 B 427 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 2183 B 2211 B 2276 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 19 B 21 B 39 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 34 KiB 34 KiB 34 KiB 76 KiB 76 KiB 76 KiB Gravity::test_composite_phi() 96 55 B 55 B 56 B 53 KiB 53 KiB 53 KiB Gravity::update_max_rhs() 192 3 B 4 B 6 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 115 B 115 B 117 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 149 B 176 B 208 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 32 KiB 32 KiB 32 KiB 38 KiB 38 KiB 38 KiB Gravity::fill_multipole_BCs() 2432 258 B 264 B 269 B 18 KiB 18 KiB 18 KiB MLCellLinOp::defineAuxData() 11704 842 B 1012 B 1890 B 8192 B 9856 B 18 KiB Castro::evaluate_source_change() 6272 1141 B 1144 B 1153 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 131 B 131 B 131 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 1824 B 1877 B 1901 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 0 B 0 B 1 B 16 KiB 16 KiB 16 KiB BndryData::define() 4864 830 B 833 B 834 B 8192 B 8192 B 8192 B FillBoundary_nowait() 91296 37 B 49 B 63 B 3072 B 3840 B 4608 B AmrMesh::MakeNewGrids() 128 0 B 29 B 47 B 3584 B 3584 B 3584 B OwnerMask() 64 0 B 0 B 0 B 3072 B 3072 B 3072 B MLCGSolver::bicgstab 640 0 B 0 B 6 B 0 B 80 B 1280 B MLMG::addInterpCorrection() 5632 1 B 17 B 27 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 0 B 12 B 15 B 1024 B 1024 B 1024 B amrex::average_down 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B MLCellLinOp::applyBC() 114375 0 B 0 B 1 B 256 B 256 B 256 B ----------------------------------------------------------------------------------------------------------------------- AMReX (26.07-55-gd1b9e4ec5555) finalized