Initializing AMReX (26.04-100-g9aaff4128aa5)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.04-100-g9aaff4128aa5) initialized Starting run at 07:27:50 UTC on 2026-04-23. Successfully read inputs file ... Castro git describe: 26.04-6-g9e0439ba0 AMReX git describe: 26.04-100-g9aaff4128 Microphysics git describe: 26.04-3-ge03e94b4 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.17966253 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001176461 new finest: 0 Gravity::fill_multipole_BCs() time = 0.011244166 Gravity::multilevel_solve_for_new_phi() time = 0.109214319 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010250436 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.048477152 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004328183 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.003514304 on level 0 Castro::do_old_sources() time = 0.007359757 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.054173342 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.003334332 on level 0 Castro::do_old_sources() time = 0.008539929 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.053477999 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.00569606 on level 0 Castro::do_old_sources() time = 0.007453181 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.050886392 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.004453007 on level 0 Castro::do_old_sources() time = 0.008566868 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.048629569 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.299977552 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.014168967 Gravity::solve_for_phi() time = 0.093255266 on level 0 Castro::construct_new_gravity() time = 0.095258404 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.049437142 on level 0 Castro::construct_old_gravity() time = 0.002691794 on level 0 Castro::do_old_sources() time = 0.008452675 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.048834775 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.301427833 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.014219218 Gravity::solve_for_phi() time = 0.090430028 on level 0 Castro::construct_new_gravity() time = 0.092403475 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.04927321 on level 0 Castro::construct_old_gravity() time = 0.004682987 on level 0 Castro::do_old_sources() time = 0.007355643 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.048621947 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.30302854 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.014044365 Gravity::solve_for_phi() time = 0.091737151 on level 0 Castro::construct_new_gravity() time = 0.093727852 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.04904874 on level 0 Castro::construct_old_gravity() time = 0.003881125 on level 0 Castro::do_old_sources() time = 0.008411542 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.048660671 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.312381605 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.014413992 Gravity::solve_for_phi() time = 0.091031844 on level 0 Castro::construct_new_gravity() time = 0.093090686 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.049616081 on level 0 Castro::construct_old_gravity() time = 0.003203349 on level 0 Castro::do_old_sources() time = 0.007334018 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.048720172 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.310685642 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.014370616 Gravity::solve_for_phi() time = 0.094598752 on level 0 Castro::construct_new_gravity() time = 0.096583767 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.049248848 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.004376241 on level 0 Castro::do_old_sources() time = 0.008431266 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.048708181 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.305355101 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.014300523 Gravity::solve_for_phi() time = 0.106140281 on level 0 Castro::construct_new_gravity() time = 0.108219007 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.054013745 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.003037056 on level 0 Castro::do_old_sources() time = 0.007495897 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.049071743 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.31346938 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.013723412 Gravity::solve_for_phi() time = 0.09083904 on level 0 Castro::construct_new_gravity() time = 0.092892135 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.053897921 on level 0 Castro::construct_old_gravity() time = 0.007693167 on level 0 Castro::do_old_sources() time = 0.007300356 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.048921244 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.302262103 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.013844831 Gravity::solve_for_phi() time = 0.089096055 on level 0 Castro::construct_new_gravity() time = 0.091105135 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.048961285 on level 0 Castro::construct_old_gravity() time = 0.002616777 on level 0 Castro::do_old_sources() time = 0.008430337 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.048569966 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.314436208 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.013927685 Gravity::solve_for_phi() time = 0.092785716 on level 0 Castro::construct_new_gravity() time = 0.094788328 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.048895212 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.002542753 on level 0 Castro::do_old_sources() time = 0.007227749 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.048697887 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.305134606 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.013678355 Gravity::solve_for_phi() time = 0.090942449 on level 0 Castro::construct_new_gravity() time = 0.092922912 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.04893822 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.002696654 on level 0 Castro::do_old_sources() time = 0.008496868 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.04868583 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.316951606 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.013790745 Gravity::solve_for_phi() time = 0.09217748 on level 0 Castro::construct_new_gravity() time = 0.094161411 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.048983707 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.002979257 on level 0 Castro::do_old_sources() time = 0.007297033 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.048576041 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.311812956 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.014159996 Gravity::solve_for_phi() time = 0.090290579 on level 0 Castro::construct_new_gravity() time = 0.092314836 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.048948428 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.002756395 on level 0 Castro::do_old_sources() time = 0.007259795 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.048794487 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.3171444 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.013956478 Gravity::solve_for_phi() time = 0.08918237 on level 0 Castro::construct_new_gravity() time = 0.09116397 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.049009052 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.002607569 on level 0 Castro::do_old_sources() time = 0.008495917 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.048697417 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.306158069 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.015013542 Gravity::solve_for_phi() time = 0.091902502 on level 0 Castro::construct_new_gravity() time = 0.093975 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.049230902 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.002781328 on level 0 Castro::do_old_sources() time = 0.008440925 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.048651549 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.31737409 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.014093713 Gravity::solve_for_phi() time = 0.088018539 on level 0 Castro::construct_new_gravity() time = 0.090021313 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.049023124 on level 0 Castro::construct_old_gravity() time = 0.002773464 on level 0 Castro::do_old_sources() time = 0.00733164 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.048735542 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.334674512 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.021223216 Gravity::solve_for_phi() time = 0.106603537 on level 0 Castro::construct_new_gravity() time = 0.108611587 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.049160236 on level 0 Castro::construct_old_gravity() time = 0.002576234 on level 0 Castro::do_old_sources() time = 0.007278207 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.048616692 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.311606795 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002901335391 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.01413772 Gravity::solve_for_phi() time = 0.092272316 on level 0 Castro::construct_new_gravity() time = 0.094340318 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.049026295 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.34452561 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.046662361 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.00408333 seconds Ending run at 07:28:04 UTC on 2026-04-23. Run time = 14.18587169 Run time without initialization = 11.39628114 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.1770 ... 14.1804 ... 14.1859 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.1285 5.1883 5.2492 37.00% Amr::InitAmr() 1 1.877 1.881 1.887 13.30% Castro::computeTemp() 80 0.7149 0.8018 0.9356 6.60% amrex::Copy() 766 0.7814 0.7903 0.7994 5.63% FabArray::ParallelCopy_finish() 1469 0.243 0.5559 0.7568 5.34% MultiFab::AddProduct() 3598 0.5482 0.554 0.5616 3.96% Castro::initMFs() 1 0.3815 0.3916 0.437 3.08% MLPoisson::Fsmooth() 6144 0.3285 0.3571 0.432 3.05% FabArray::setVal() 5232 0.3483 0.3542 0.3621 2.55% Castro::reset_internal_energy(Fab) 160 0.3463 0.3516 0.359 2.53% MultiFab::Subtract() 159 0.3393 0.3433 0.3532 2.49% Castro::do_advance_ctu() 20 0.2593 0.2627 0.2685 1.89% Gravity::fill_multipole_BCs() 19 0.1967 0.2432 0.2663 1.88% check_for_negative_density() 20 0.1275 0.192 0.2605 1.84% FabArray::ParallelCopy_nowait() 1469 0.2154 0.2369 0.2604 1.84% FillBoundary_finish() 6380 0.1448 0.1623 0.1966 1.39% FabArray::sum() 3598 0.1773 0.1792 0.1823 1.28% Castro::initData() 1 0.003206 0.02892 0.1703 1.20% FillBoundary_nowait() 6380 0.08991 0.1201 0.1544 1.09% TagBoxArray::collate() 2 0.0002022 0.02786 0.1526 1.08% MLCellLinOp::applyBC() 7148 0.06857 0.08481 0.1311 0.92% Castro::normalize_species() 79 0.1192 0.1208 0.1226 0.86% MLPoisson::Fapply() 986 0.08685 0.09162 0.1061 0.75% Castro::construct_old_gravity_source() 40 0.09221 0.09339 0.09493 0.67% FabArray::Saxpy() 111 0.07575 0.07666 0.07737 0.55% Castro::construct_new_gravity_source() 17 0.06724 0.06809 0.06973 0.49% FabArray::Xpay() 915 0.04908 0.05145 0.05676 0.40% MultiFab::contains_nan() 37 0.04382 0.04434 0.04492 0.32% MLMG::addInterpCorrection() 768 0.03194 0.0347 0.04208 0.30% AmrLevel::derive() 18 0.03859 0.03919 0.03986 0.28% StateData::FillBoundary(geom) 193 0 0.01854 0.0358 0.25% amrex::average_down 768 0.02459 0.02692 0.03362 0.24% Castro::reset_internal_energy(MultiFab) 80 0.0252 0.02581 0.02668 0.19% VisMF::Write(FabArray) 4 0.01284 0.01417 0.02538 0.18% Castro::construct_old_source() 240 0.0193 0.01957 0.02097 0.15% FabArray::FillBoundary() 6380 0.01689 0.0175 0.02075 0.15% MLCellLinOp::defineAuxData() 19 0.01601 0.01717 0.02072 0.15% Castro::enforce_min_density() 79 0.01932 0.01994 0.02046 0.14% Amr::writePlotFile() 2 0.006615 0.01668 0.01822 0.13% Castro::construct_new_source() 102 0.01511 0.01537 0.01775 0.13% FabArrayBase::getFB() 6380 0.01345 0.01402 0.01713 0.12% MLCGSolver::bicgstab 8 0 0.001069 0.0171 0.12% FabArray::setDomainBndry() 187 0.0006807 0.008836 0.01704 0.12% FillPatchIterator::Initialize 187 0.01593 0.01646 0.01689 0.12% Castro::evaluate_source_change() 196 0.01468 0.01506 0.01657 0.12% MLMG::ResNormInf() 146 0.00991 0.01226 0.01547 0.11% FabArray::norminf() 244 0.01056 0.01106 0.01454 0.10% amrex::Add() 128 0.01356 0.01377 0.01418 0.10% MLCellLinOp::smooth() 1544 0.01056 0.01102 0.01342 0.09% MLMG::prepareForSolve() 18 0.01015 0.01064 0.01191 0.08% Castro::wd_update() 2 0.004547 0.008115 0.01067 0.08% MLMG::solve() 18 0.002351 0.004444 0.009821 0.07% FabArrayBase::getCPC() 2045 0.008827 0.009084 0.009816 0.07% Gravity::actual_solve_with_mlmg() 19 0.00751 0.008056 0.009721 0.07% Gravity::get_old_grav_vector() 20 0.00616 0.006267 0.006489 0.05% Gravity::get_new_grav_vector() 18 0.005791 0.005933 0.00622 0.04% MLCellLinOp::apply() 986 0.002391 0.002659 0.00563 0.04% Castro::initialize_do_advance() 20 0.004714 0.004931 0.005212 0.04% FabArray::ParallelCopy() 1469 0.00427 0.004441 0.004948 0.03% AmrLevel::FillPatch() 187 0.003051 0.003254 0.003662 0.03% MLLinOp::defineGrids() 19 0.00256 0.002709 0.003044 0.02% MLCellLinOp::correctionResidual() 768 0.002281 0.002346 0.002549 0.02% MLMG::mgVcycle() 128 0.002146 0.00223 0.002415 0.02% Castro::sum_integrated_quantities() 2 0.001579 0.001864 0.002085 0.01% MLMG:computeResOfCorrection() 768 0.001793 0.001859 0.002076 0.01% Castro::save_data_for_retry() 3 0.001325 0.001366 0.001514 0.01% StateDataPhysBCFunct::() 187 0.0009868 0.00119 0.001458 0.01% MLMG::actualBottomSolve() 128 0.0001892 0.0002654 0.001286 0.01% Castro::subcycle_advance_ctu() 1 0.0002388 0.000336 0.001079 0.01% FillPatchIterator::FillFromLevel0() 187 0.0008849 0.0009265 0.001042 0.01% main() 1 0.0007425 0.000836 0.000996 0.01% Gravity::solve_for_phi() 17 0.0007081 0.0007308 0.0007833 0.01% MLMG::mgVcycle_down::2 128 0.0005942 0.0006303 0.0007373 0.01% Castro::Castro() 1 0.0004123 0.0004553 0.0007122 0.01% FillPatchSingleLevel 187 0.0006654 0.0006783 0.0007042 0.00% MLMG::mgVcycle_down::0 128 0.0006358 0.0006521 0.0007033 0.00% MLMG::mgVcycle_down::1 128 0.0006068 0.000626 0.0006958 0.00% MLMG::mgVcycle_down::3 128 0.0005647 0.000595 0.0006781 0.00% MLMG::mgVcycle_down::5 128 0.0005867 0.0006084 0.0006757 0.00% MLMG::mgVcycle_down::4 128 0.000559 0.0005776 0.0006565 0.00% AmrMesh::MakeNewGrids() 2 0.0005665 0.0005977 0.0006338 0.00% MLCellLinOp::solutionResidual() 147 0.0005223 0.0005391 0.000576 0.00% Castro::clean_state() 79 0.0004781 0.0005028 0.0005744 0.00% Castro::do_old_sources() 40 0.0004626 0.0004781 0.0005208 0.00% MLMG::oneIter() 128 0.0004314 0.0004494 0.0004981 0.00% MLMG::mgVcycle_up::2 128 0.0003916 0.0004126 0.0004871 0.00% MLMG::mgVcycle_up::1 128 0.0004074 0.0004232 0.0004833 0.00% MLMG::mgVcycle_up::0 128 0.0004103 0.0004261 0.0004604 0.00% MLMG::mgVcycle_up::4 128 0.0003814 0.0003928 0.0004401 0.00% MLMG::mgVcycle_up::3 128 0.000379 0.0003902 0.0004378 0.00% MLMG::mgVcycle_up::5 128 0.0003669 0.0003803 0.0004202 0.00% MLMG::computeResidual() 128 0.000294 0.0003046 0.0003244 0.00% MLMG::mgVcycle_bottom 128 0.0002745 0.0002902 0.0003218 0.00% Castro::construct_new_gravity() 17 0.0001917 0.0002081 0.000311 0.00% Gravity::solve_phi_with_mlmg() 19 0.000214 0.0002358 0.000286 0.00% Castro::do_new_sources() 17 0.0002205 0.0002346 0.0002723 0.00% Castro::problem_post_timestep() 2 1.139e-05 2.793e-05 0.0002398 0.00% Amr::coarseTimeStep() 1 0.0001397 0.000163 0.0001663 0.00% Castro::check_for_nan() 37 0.0001193 0.0001329 0.0001608 0.00% Castro::construct_old_gravity() 20 0.000128 0.000138 0.0001535 0.00% Castro::initialize_advance() 1 0.000108 0.0001274 0.0001505 0.00% MLPoisson::define() 19 0.0001268 0.0001358 0.0001479 0.00% TagBoxArray::mapPRD 2 0.0001205 0.0001255 0.0001446 0.00% Castro::apply_source_to_state() 37 0.0001164 0.0001233 0.0001422 0.00% Castro::retry_advance_ctu() 11 7.875e-05 9.105e-05 0.0001343 0.00% Gravity::test_composite_phi() 1 0.0001073 0.0001125 0.0001296 0.00% Amr::defBaseLevel() 1 9.732e-05 0.0001064 0.00012 0.00% MLLinOp::define() 19 8.672e-05 9.216e-05 9.919e-05 0.00% Castro::create_source_corrector() 17 5.897e-05 6.308e-05 7.158e-05 0.00% Amr::grid_places() 2 5.299e-05 5.851e-05 7.001e-05 0.00% MLMG::computeMLResidual() 18 5.556e-05 6.067e-05 6.569e-05 0.00% Castro::expand_state() 20 6.015e-05 6.272e-05 6.507e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.807e-05 6.008e-05 6.406e-05 0.00% Amr::FinalizeInit() 1 8.931e-06 1.262e-05 5.318e-05 0.00% Castro::derive() 18 4.611e-05 4.801e-05 5.141e-05 0.00% Amr::regrid() 1 1.322e-05 1.97e-05 4.293e-05 0.00% Gravity::actual_multilevel_solve() 1 3.184e-05 3.347e-05 3.775e-05 0.00% Castro::post_init() 1 1.3e-05 1.486e-05 3.017e-05 0.00% Amr::timeStep() 1 1.555e-05 1.736e-05 2.132e-05 0.00% Amr::InitializeInit() 1 9.803e-06 1.401e-05 1.84e-05 0.00% Amr::bldFineLevels() 1 7.645e-06 8.82e-06 1.387e-05 0.00% Castro::post_regrid() 2 1.1e-05 1.188e-05 1.381e-05 0.00% Castro::advance() 1 8.249e-06 8.909e-06 1.155e-05 0.00% Amr::init() 1 6.943e-06 7.682e-06 8.303e-06 0.00% Castro::post_timestep() 1 6.045e-06 6.393e-06 6.896e-06 0.00% Amr::initialInit() 1 3.181e-06 3.633e-06 4.178e-06 0.00% Other 2530 0.04833 0.06035 0.09861 0.70% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.18 14.18 14.19 100.00% Amr::coarseTimeStep() 1 11.4 11.4 11.4 80.33% Amr::timeStep() 1 11.34 11.34 11.34 79.97% Castro::advance() 1 11.29 11.29 11.29 79.62% Castro::subcycle_advance_ctu() 1 11.25 11.25 11.25 79.33% Castro::do_advance_ctu() 20 11.21 11.22 11.22 79.08% Castro::construct_ctu_hydro_source() 20 5.993 6.013 6.055 42.68% Castro::clean_state() 79 2.932 3.009 3.128 22.05% Castro::computeTemp() 80 2.049 2.123 2.249 15.85% Amr::InitAmr() 1 1.877 1.881 1.887 13.30% Gravity::solve_phi_with_mlmg() 19 1.631 1.678 1.699 11.98% Castro::construct_new_gravity() 17 1.546 1.593 1.614 11.38% Gravity::solve_for_phi() 17 1.511 1.558 1.58 11.14% Gravity::actual_solve_with_mlmg() 19 1.428 1.431 1.436 10.12% MLMG::solve() 18 1.369 1.375 1.377 9.71% Castro::reset_internal_energy(MultiFab) 80 1.305 1.322 1.335 9.41% MLMG::oneIter() 128 1.231 1.233 1.235 8.70% MLMG::mgVcycle() 128 1.217 1.219 1.221 8.61% Castro::initialize_do_advance() 20 1.106 1.152 1.184 8.35% Castro::do_old_sources() 40 1.087 1.106 1.125 7.93% FabArray::ParallelCopy() 1469 0.5046 0.8179 1.017 7.17% Castro::evaluate_source_change() 196 0.9316 0.9422 0.9567 6.74% Amr::init() 1 0.9029 0.9029 0.9029 6.36% Amr::initialInit() 1 0.8496 0.8496 0.8497 5.99% Castro::do_new_sources() 17 0.772 0.805 0.8428 5.94% MLCellLinOp::smooth() 1544 0.634 0.6911 0.8311 5.86% Castro::enforce_min_density() 79 0.781 0.7889 0.802 5.65% amrex::Copy() 766 0.7814 0.7903 0.7994 5.63% FabArray::ParallelCopy_finish() 1469 0.243 0.5559 0.7568 5.34% AmrLevel::FillPatch() 187 0.5747 0.6352 0.6889 4.86% Amr::InitializeInit() 1 0.455 0.4922 0.6333 4.46% Amr::defBaseLevel() 1 0.455 0.4921 0.6333 4.46% MultiFab::AddProduct() 3598 0.5482 0.554 0.5616 3.96% FillPatchIterator::Initialize 187 0.3817 0.4399 0.4961 3.50% MLCellLinOp::applyBC() 7148 0.348 0.3994 0.4845 3.42% MLMG::mgVcycle_up::1 128 0.1096 0.2928 0.4744 3.34% FillPatchIterator::FillFromLevel0() 187 0.3498 0.4146 0.4688 3.30% FillPatchSingleLevel 187 0.3489 0.4137 0.4679 3.30% Castro::Castro() 1 0.3866 0.3968 0.4427 3.12% Castro::initMFs() 1 0.3817 0.3918 0.4372 3.08% MLPoisson::Fsmooth() 6144 0.3285 0.3571 0.432 3.05% MLMG::addInterpCorrection() 768 0.06766 0.3297 0.4108 2.90% Amr::FinalizeInit() 1 0.2163 0.3575 0.3946 2.78% FabArray::setVal() 5232 0.3483 0.3542 0.3621 2.55% FabArray::FillBoundary() 6380 0.2671 0.3146 0.3611 2.55% Castro::reset_internal_energy(Fab) 160 0.3463 0.3516 0.359 2.53% MultiFab::Subtract() 159 0.3393 0.3433 0.3532 2.49% MLMG::mgVcycle_up::2 128 0.01556 0.1302 0.3088 2.18% FabArray::ParallelCopy_nowait() 1469 0.2352 0.2575 0.282 1.99% MLMG::mgVcycle_down::0 128 0.2694 0.2732 0.2757 1.94% Gravity::fill_multipole_BCs() 19 0.1987 0.2452 0.2682 1.89% check_for_negative_density() 20 0.1275 0.192 0.2605 1.84% Castro::expand_state() 20 0.1864 0.215 0.2384 1.68% Castro::initData() 1 0.06652 0.09359 0.235 1.66% Amr::bldFineLevels() 1 0.03851 0.1791 0.2157 1.52% MLCellLinOp::apply() 986 0.156 0.1687 0.2067 1.46% MLMG::mgVcycle_up::0 128 0.2019 0.2029 0.2043 1.44% MLMG::mgVcycle_up::3 128 0.01542 0.02957 0.1988 1.40% FillBoundary_finish() 6380 0.1448 0.1623 0.1966 1.39% FabArray::sum() 3598 0.1773 0.1792 0.1823 1.28% Amr::grid_places() 2 0.002552 0.1437 0.1808 1.27% AmrMesh::MakeNewGrids() 2 0.002496 0.1436 0.1808 1.27% TagBoxArray::mapPRD 2 0.0007712 0.1144 0.1789 1.26% FillBoundary_nowait() 6380 0.1038 0.1348 0.1696 1.20% Castro::construct_old_source() 240 0.1533 0.1553 0.1592 1.12% TagBoxArray::collate() 2 0.0002022 0.02786 0.1526 1.08% MLMG:computeResOfCorrection() 768 0.1189 0.1294 0.1515 1.07% MLCellLinOp::correctionResidual() 768 0.1171 0.1275 0.1496 1.05% MLMG::mgVcycle_down::1 128 0.1326 0.1377 0.1451 1.02% Castro::post_init() 1 0.1418 0.1433 0.1447 1.02% Castro::normalize_species() 79 0.1192 0.1208 0.1226 0.86% Castro::construct_new_source() 102 0.1177 0.1191 0.1209 0.85% Gravity::multilevel_solve_for_new_phi() 1 0.1091 0.1092 0.1092 0.77% Gravity::actual_multilevel_solve() 1 0.1091 0.1091 0.1092 0.77% MLPoisson::Fapply() 986 0.08685 0.09162 0.1061 0.75% MLCellLinOp::solutionResidual() 147 0.08529 0.09353 0.09921 0.70% Amr::writePlotFile() 2 0.09554 0.09563 0.09564 0.67% Castro::construct_old_gravity_source() 40 0.09221 0.09339 0.09493 0.67% MLMG::mgVcycle_down::2 128 0.01807 0.04954 0.08544 0.60% MLMG::mgVcycle_bottom 128 0.0004666 0.005615 0.08253 0.58% MLMG::actualBottomSolve() 128 0.0001892 0.005324 0.08223 0.58% MLCGSolver::bicgstab 8 0 0.005012 0.08019 0.57% Castro::create_source_corrector() 17 0.05038 0.06445 0.07831 0.55% FabArray::Saxpy() 111 0.07575 0.07666 0.07737 0.55% MLMG::computeResidual() 128 0.07427 0.07538 0.07652 0.54% Castro::post_regrid() 2 0.06903 0.07039 0.07187 0.51% Castro::construct_new_gravity_source() 17 0.06724 0.06809 0.06973 0.49% Castro::derive() 18 0.0664 0.06735 0.06912 0.49% MLMG::mgVcycle_down::3 128 0.01751 0.02411 0.06908 0.49% AmrLevel::derive() 18 0.06635 0.06731 0.06907 0.49% Castro::construct_old_gravity() 20 0.04581 0.05379 0.06103 0.43% Gravity::get_old_grav_vector() 20 0.04567 0.05365 0.0609 0.43% amrex::average_down 768 0.03456 0.04247 0.06056 0.43% FabArray::Xpay() 915 0.04908 0.05145 0.05676 0.40% Castro::post_timestep() 1 0.0497 0.05166 0.05379 0.38% Castro::check_for_nan() 37 0.04395 0.04448 0.04506 0.32% MultiFab::contains_nan() 37 0.04382 0.04434 0.04492 0.32% MLPoisson::define() 19 0.03346 0.03864 0.04427 0.31% Castro::initialize_advance() 1 0.03993 0.04056 0.04122 0.29% MLMG::mgVcycle_down::4 128 0.01752 0.01934 0.04105 0.29% Amr::regrid() 1 0.03592 0.03655 0.03719 0.26% StateDataPhysBCFunct::() 187 0.0009868 0.01973 0.03704 0.26% Castro::retry_advance_ctu() 11 0.03424 0.03492 0.03642 0.26% Gravity::get_new_grav_vector() 18 0.03381 0.03467 0.03589 0.25% StateData::FillBoundary(geom) 193 0 0.01854 0.0358 0.25% MLMG::mgVcycle_up::4 128 0.01548 0.0169 0.03276 0.23% MLMG::mgVcycle_down::5 128 0.01752 0.01867 0.0303 0.21% MLCellLinOp::defineAuxData() 19 0.02073 0.02269 0.02821 0.20% MLMG::mgVcycle_up::5 128 0.01556 0.01652 0.02596 0.18% Castro::apply_source_to_state() 37 0.02524 0.02552 0.02588 0.18% VisMF::Write(FabArray) 4 0.01284 0.01417 0.02538 0.18% MLMG::ResNormInf() 146 0.01932 0.0217 0.02488 0.18% MLMG::prepareForSolve() 18 0.02008 0.02128 0.02426 0.17% MLMG::computeMLResidual() 18 0.01082 0.01771 0.02288 0.16% Castro::sum_integrated_quantities() 2 0.02096 0.02109 0.0211 0.15% Castro::save_data_for_retry() 3 0.01793 0.0183 0.01957 0.14% FabArrayBase::getFB() 6380 0.0138 0.01473 0.01861 0.13% FabArray::setDomainBndry() 187 0.0006807 0.008836 0.01704 0.12% FabArrayBase::getCPC() 2045 0.01126 0.01265 0.01619 0.11% FabArray::norminf() 244 0.01056 0.01106 0.01454 0.10% amrex::Add() 128 0.01356 0.01377 0.01418 0.10% Gravity::test_composite_phi() 1 0.01396 0.01402 0.01406 0.10% Castro::problem_post_timestep() 2 0.009643 0.01186 0.01405 0.10% Castro::wd_update() 2 0.009632 0.01183 0.01403 0.10% MLLinOp::define() 19 0.005638 0.00916 0.01262 0.09% MLLinOp::defineGrids() 19 0.005546 0.009068 0.01253 0.09% Other 2530 0.07938 0.09247 0.1369 0.97% -------------------------------------------------------------------------------------------- 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 387 KiB 392 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 2100 B 2113 B 2125 B 89 KiB 89 KiB 89 KiB -------------------------------------------------------------------------------------------------------------------- Pinned Memory Usage: -------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max -------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 1 B 1 B 3 B 8192 KiB 8192 KiB 8192 KiB -------------------------------------------------------------------------------------------------------------- Cpu Memory Usage: ----------------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max ----------------------------------------------------------------------------------------------------------------------- Castro::save_data_for_retry() 416 1518 KiB 1519 KiB 1520 KiB 2332 KiB 2332 KiB 2332 KiB Castro::construct_ctu_hydro_source() 6400 384 KiB 388 KiB 393 KiB 1178 KiB 1178 KiB 1178 KiB StateData::define() 192 1009 KiB 1010 KiB 1010 KiB 1166 KiB 1166 KiB 1166 KiB Castro::initMFs() 256 989 KiB 989 KiB 989 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 645 KiB 645 KiB 645 KiB 804 KiB 804 KiB 804 KiB Amr::writePlotFile() 64 4274 B 4357 B 4942 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 49 KiB 49 KiB 50 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 415 KiB 416 KiB 416 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2180 B 3969 B 5140 B 226 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1356 B 1377 B 1422 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 125 KiB 125 KiB 125 KiB 158 KiB 158 KiB 158 KiB Castro::buildMetrics() 192 123 KiB 123 KiB 123 KiB 143 KiB 143 KiB 143 KiB MLMG::prepareForSolve() 6102 10 KiB 10 KiB 12 KiB 104 KiB 110 KiB 130 KiB Castro::construct_old_source() 1280 477 B 483 B 512 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 400 B 404 B 423 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 2106 B 2118 B 2131 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 17 B 20 B 38 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 32 KiB 32 KiB 32 KiB 76 KiB 76 KiB 76 KiB Gravity::test_composite_phi() 96 52 B 52 B 53 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 108 B 110 B 115 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 148 B 173 B 198 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 31 KiB 31 KiB 31 KiB 38 KiB 38 KiB 38 KiB Gravity::fill_multipole_BCs() 2432 243 B 250 B 257 B 18 KiB 18 KiB 18 KiB MLCellLinOp::defineAuxData() 11704 814 B 978 B 1826 B 8192 B 9856 B 18 KiB Castro::evaluate_source_change() 6272 1069 B 1081 B 1098 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 125 B 125 B 125 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 1743 B 1797 B 1822 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 802 B 805 B 806 B 8192 B 8192 B 8192 B FillBoundary_nowait() 91296 31 B 45 B 59 B 3072 B 3840 B 4608 B AmrMesh::MakeNewGrids() 128 0 B 24 B 39 B 3584 B 3584 B 3584 B OwnerMask() 64 0 B 0 B 0 B 3072 B 3072 B 3072 B MLCGSolver::bicgstab 640 0 B 0 B 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 9 B 12 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.04-100-g9aaff4128aa5) finalized