Initializing AMReX (25.06-22-g1d307c2c661a-dirty)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (25.06-22-g1d307c2c661a-dirty) initialized Starting run at 06:48:25 UTC on 2025-07-24. Successfully read inputs file ... Castro git describe: 25.07-12-g2c8f2fc1e AMReX git describe: 25.06-22-g1d307c2c6-dirty Microphysics git describe: joss_submit 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.178874506 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.00117933 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010681883 Gravity::multilevel_solve_for_new_phi() time = 0.109038804 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.01064859 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.04952769 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004487188 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.002880651 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.06563287 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.003575274 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.063434753 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.003370192 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.060327277 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.003483658 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.059451595 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.331963214 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.014068015 Gravity::solve_for_phi() time = 0.09164342 on level 0 Castro::construct_new_gravity() time = 0.093768142 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.060074659 on level 0 Castro::construct_old_gravity() time = 0.002701537 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.058880626 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.331679713 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.013816389 Gravity::solve_for_phi() time = 0.090392594 on level 0 Castro::construct_new_gravity() time = 0.092529622 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.05950243 on level 0 Castro::construct_old_gravity() time = 0.002944557 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.05944584 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.338397171 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.013980845 Gravity::solve_for_phi() time = 0.090554032 on level 0 Castro::construct_new_gravity() time = 0.09267892 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.059699736 on level 0 Castro::construct_old_gravity() time = 0.004245712 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.059483409 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.342284277 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.014684566 Gravity::solve_for_phi() time = 0.092001053 on level 0 Castro::construct_new_gravity() time = 0.094112182 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.059591494 on level 0 Castro::construct_old_gravity() time = 0.003220159 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.060075756 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.349983577 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.014598592 Gravity::solve_for_phi() time = 0.092239175 on level 0 Castro::construct_new_gravity() time = 0.094405173 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.05975551 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.003408958 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.059617636 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.346144025 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.014211231 Gravity::solve_for_phi() time = 0.09165241 on level 0 Castro::construct_new_gravity() time = 0.093815418 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.059638907 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.003423337 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.059488199 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.345125252 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.014410958 Gravity::solve_for_phi() time = 0.091643233 on level 0 Castro::construct_new_gravity() time = 0.093849477 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.059834929 on level 0 Castro::construct_old_gravity() time = 0.003115473 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.05944553 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.344302676 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.014413809 Gravity::solve_for_phi() time = 0.091645489 on level 0 Castro::construct_new_gravity() time = 0.093805515 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.060648235 on level 0 Castro::construct_old_gravity() time = 0.002790867 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.059426892 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.347346344 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.014754696 Gravity::solve_for_phi() time = 0.091548949 on level 0 Castro::construct_new_gravity() time = 0.093668674 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.05973245 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.002728203 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.059454961 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.34321781 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.014886668 Gravity::solve_for_phi() time = 0.091552831 on level 0 Castro::construct_new_gravity() time = 0.093671749 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.060695814 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.002674728 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.059508967 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.348600674 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.014878425 Gravity::solve_for_phi() time = 0.092424958 on level 0 Castro::construct_new_gravity() time = 0.094554757 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.059412592 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.003642089 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.058822217 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.345501299 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.014585475 Gravity::solve_for_phi() time = 0.092706668 on level 0 Castro::construct_new_gravity() time = 0.09483649 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.060926245 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.002790392 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.059275809 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.348788525 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.01482054 Gravity::solve_for_phi() time = 0.094137394 on level 0 Castro::construct_new_gravity() time = 0.096261377 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.061107152 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.004673145 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.061221156 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.356678221 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.015383281 Gravity::solve_for_phi() time = 0.093554872 on level 0 Castro::construct_new_gravity() time = 0.095704922 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.059738979 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.002688771 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.059330382 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.349160016 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.015056939 Gravity::solve_for_phi() time = 0.092953141 on level 0 Castro::construct_new_gravity() time = 0.095086289 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.060250685 on level 0 Castro::construct_old_gravity() time = 0.002734358 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.059719198 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.345591719 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.015410359 Gravity::solve_for_phi() time = 0.09399517 on level 0 Castro::construct_new_gravity() time = 0.096126333 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.062566715 on level 0 Castro::construct_old_gravity() time = 0.002790132 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.059585648 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.350302244 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002684001322 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.015584276 Gravity::solve_for_phi() time = 0.09236565 on level 0 Castro::construct_new_gravity() time = 0.094486708 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.059744517 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.164437349e+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.469212319e+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.566427124e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 57600050.72 TIME= 0.25 MAXIMUM DENSITY = 4458760.788 TIME= 0.25 MAXIMUM T_S / T_E = 0 STEP = 1 TIME = 0.25 DT = 0.25 [STEP 1] Coarse TimeStep time: 12.26736323 [STEP 1] FAB kilobyte spread across MPI nodes: [6998 ... 6998] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.048946468 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004064352 seconds Ending run at 06:48:39 UTC on 2025-07-24. Run time = 14.54094497 Run time without initialization = 12.32129285 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.5318 ... 14.5344 ... 14.5411 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.7227 5.7470 5.7965 39.86% Amr::InitAmr() 1 1.3 1.303 1.309 9.01% Castro::computeTemp() 80 0.8423 0.9567 1.112 7.64% amrex::Copy() 724 0.7983 0.8036 0.8173 5.62% FabArray::ParallelCopy_finish() 1429 0.2226 0.5596 0.7654 5.26% MultiFab::AddProduct() 3439 0.5514 0.5545 0.564 3.88% Castro::normalize_species() 79 0.5475 0.5502 0.5555 3.82% Castro::reset_internal_energy(Fab) 160 0.5203 0.5253 0.5309 3.65% MLPoisson::Fsmooth() 6144 0.3544 0.3855 0.4757 3.27% Castro::initMFs() 1 0.3579 0.3652 0.373 2.57% FabArray::setVal() 5013 0.348 0.3515 0.359 2.47% MultiFab::Subtract() 159 0.3368 0.3388 0.3438 2.36% Gravity::fill_multipole_BCs() 19 0.1963 0.2459 0.2661 1.83% FabArray::ParallelCopy_nowait() 1429 0.1992 0.2167 0.2389 1.64% check_for_negative_density() 20 0.08256 0.1477 0.2062 1.42% FabArray::sum() 3439 0.197 0.1985 0.2011 1.38% Castro::initData() 1 0.003361 0.02823 0.1669 1.15% TagBoxArray::collate() 2 0.0002818 0.03395 0.1594 1.10% FillBoundary_nowait() 6380 0.09339 0.1204 0.1552 1.07% FillBoundary_finish() 6380 0.09684 0.133 0.1487 1.02% MLCellLinOp::applyBC() 7148 0.0635 0.08028 0.127 0.87% MLPoisson::Fapply() 986 0.09491 0.09921 0.1153 0.79% Castro::construct_new_gravity_source() 17 0.07243 0.07319 0.07458 0.51% FabArray::Saxpy() 91 0.0722 0.07235 0.07277 0.50% FabArray::Xpay() 915 0.0515 0.05476 0.06356 0.44% Castro::construct_old_gravity_source() 20 0.04914 0.04965 0.05087 0.35% MultiFab::contains_nan() 37 0.04581 0.04621 0.04653 0.32% MLMG::addInterpCorrection() 768 0.03377 0.03628 0.04398 0.30% AmrLevel::derive() 18 0.04117 0.04159 0.04217 0.29% amrex::average_down 768 0.02507 0.02732 0.03413 0.23% StateData::FillBoundary(geom) 143 0 0.01754 0.03361 0.23% VisMF::Write(FabArray) 4 0.01402 0.01517 0.02508 0.17% Castro::reset_internal_energy(MultiFab) 80 0.02306 0.02343 0.0241 0.17% MLCellLinOp::defineAuxData() 19 0.0155 0.01654 0.02009 0.14% Castro::enforce_min_density() 79 0.01855 0.0187 0.019 0.13% FabArray::FillBoundary() 6380 0.01395 0.0145 0.01746 0.12% Amr::writePlotFile() 2 0.006816 0.01547 0.01695 0.12% MLCGSolver::bicgstab 8 0 0.001032 0.01651 0.11% Castro::construct_new_source() 102 0.01612 0.0162 0.01642 0.11% Castro::evaluate_source_change() 196 0.01389 0.01428 0.01512 0.10% amrex::Add() 128 0.01444 0.01454 0.01476 0.10% FabArray::setDomainBndry() 147 0.0004935 0.007597 0.0145 0.10% FabArray::norminf() 244 0.0107 0.01111 0.01406 0.10% FillPatchIterator::Initialize 147 0.01276 0.0131 0.01347 0.09% FabArrayBase::getFB() 6380 0.01057 0.01108 0.01343 0.09% Castro::wd_update() 2 0.005628 0.008916 0.01275 0.09% MLCellLinOp::smooth() 1544 0.01021 0.01074 0.01218 0.08% MLMG::prepareForSolve() 18 0.01019 0.01053 0.01152 0.08% MLMG::ResNormInf() 146 0.003212 0.006842 0.01043 0.07% Castro::construct_old_source() 120 0.01026 0.01032 0.01042 0.07% Gravity::actual_solve_with_mlmg() 19 0.007305 0.007668 0.009088 0.06% MLMG::solve() 18 0.0007495 0.003991 0.00864 0.06% FabArrayBase::getCPC() 2005 0.007717 0.007904 0.008225 0.06% Gravity::get_old_grav_vector() 20 0.006467 0.006552 0.006755 0.05% Gravity::get_new_grav_vector() 18 0.00603 0.006114 0.006264 0.04% Castro::initialize_do_advance() 20 0.004382 0.004675 0.005016 0.03% MLCellLinOp::apply() 986 0.002038 0.002279 0.00488 0.03% FabArray::ParallelCopy() 1429 0.003628 0.003778 0.004051 0.03% AmrLevel::FillPatch() 147 0.002338 0.002477 0.002709 0.02% MLCellLinOp::correctionResidual() 768 0.001889 0.001958 0.002073 0.01% Castro::sum_integrated_quantities() 2 0.001569 0.001842 0.002043 0.01% MLMG::mgVcycle() 128 0.001649 0.001723 0.001874 0.01% MLMG:computeResOfCorrection() 768 0.001501 0.001561 0.001664 0.01% Castro::save_data_for_retry() 3 0.001293 0.001323 0.001349 0.01% MLMG::actualBottomSolve() 128 0.000132 0.0002021 0.001209 0.01% Castro::subcycle_advance_ctu() 1 0.0002285 0.0003061 0.001121 0.01% StateDataPhysBCFunct::() 147 0.0007293 0.0008816 0.00102 0.01% main() 1 0.0006622 0.0007273 0.0008666 0.01% FillPatchIterator::FillFromLevel0() 147 0.000693 0.0007198 0.0007784 0.01% Gravity::solve_for_phi() 17 0.000651 0.0006867 0.0007545 0.01% MLMG::mgVcycle_down::1 128 0.0005591 0.0005745 0.0006172 0.00% MLMG::mgVcycle_down::3 128 0.0005044 0.0005425 0.0006167 0.00% Castro::Castro() 1 0.0004647 0.0005445 0.0006141 0.00% MLMG::mgVcycle_down::2 128 0.0005416 0.0005788 0.00061 0.00% MLMG::mgVcycle_down::0 128 0.0005655 0.000577 0.0006025 0.00% MLMG::mgVcycle_down::5 128 0.0005104 0.0005263 0.0005952 0.00% MLMG::mgVcycle_down::4 128 0.0004869 0.000509 0.0005883 0.00% AmrMesh::MakeNewGrids() 2 0.0005428 0.0005634 0.0005792 0.00% FillPatchSingleLevel 147 0.0004692 0.0004833 0.0005122 0.00% Castro::clean_state() 79 0.0004256 0.0004554 0.0005055 0.00% MLCellLinOp::solutionResidual() 147 0.0004439 0.0004605 0.0004925 0.00% Castro::do_advance_ctu() 20 0.000379 0.0004172 0.0004545 0.00% MLMG::oneIter() 128 0.0003909 0.0004032 0.0004417 0.00% MLMG::mgVcycle_up::0 128 0.0003761 0.0003861 0.0004249 0.00% MLMG::mgVcycle_up::1 128 0.0003667 0.0003857 0.000406 0.00% MLMG::mgVcycle_up::2 128 0.0003243 0.0003523 0.0003862 0.00% MLMG::mgVcycle_up::4 128 0.0003233 0.0003335 0.0003722 0.00% MLMG::mgVcycle_up::3 128 0.0003219 0.0003335 0.0003642 0.00% MLMG::mgVcycle_up::5 128 0.0003116 0.0003188 0.0003412 0.00% Castro::do_old_sources() 20 0.0002777 0.0002918 0.0003166 0.00% Castro::construct_new_gravity() 17 0.0001628 0.0001775 0.0002982 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002174 0.0002331 0.0002674 0.00% Castro::problem_post_timestep() 2 1.072e-05 2.716e-05 0.0002603 0.00% MLMG::computeResidual() 128 0.0002375 0.0002458 0.0002584 0.00% Castro::do_new_sources() 17 0.0002131 0.0002284 0.0002504 0.00% MLMG::mgVcycle_bottom 128 0.0002264 0.0002315 0.0002452 0.00% MLPoisson::define() 19 0.000132 0.000146 0.0001858 0.00% Amr::coarseTimeStep() 1 0.0001374 0.0001653 0.0001729 0.00% Castro::construct_old_gravity() 20 0.0001088 0.000122 0.0001426 0.00% TagBoxArray::mapPRD 2 0.000124 0.0001335 0.0001419 0.00% Castro::check_for_nan() 37 0.0001066 0.0001158 0.0001404 0.00% Gravity::test_composite_phi() 1 0.0001052 0.0001102 0.0001266 0.00% Castro::initialize_advance() 1 0.0001037 0.0001124 0.0001255 0.00% Castro::retry_advance_ctu() 11 6.223e-05 7.224e-05 0.0001215 0.00% Amr::defBaseLevel() 1 9.739e-05 0.0001077 0.0001179 0.00% Castro::apply_source_to_state() 37 9.593e-05 0.0001019 0.0001124 0.00% MLLinOp::define() 19 7.406e-05 8.169e-05 9.512e-05 0.00% Castro::create_source_corrector() 17 5.464e-05 6.412e-05 8.744e-05 0.00% Amr::grid_places() 2 5.334e-05 5.972e-05 7.282e-05 0.00% Amr::FinalizeInit() 1 1.533e-05 2.708e-05 7.265e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.833e-05 6.091e-05 6.425e-05 0.00% Castro::expand_state() 20 4.837e-05 5.307e-05 5.811e-05 0.00% MLMG::computeMLResidual() 18 4.122e-05 4.352e-05 4.615e-05 0.00% Castro::derive() 18 4.216e-05 4.342e-05 4.596e-05 0.00% Amr::regrid() 1 1.567e-05 2.013e-05 4.347e-05 0.00% Gravity::actual_multilevel_solve() 1 3.172e-05 3.27e-05 3.466e-05 0.00% Castro::post_regrid() 2 1.139e-05 1.43e-05 2.248e-05 0.00% Amr::timeStep() 1 1.413e-05 1.546e-05 1.999e-05 0.00% Amr::InitializeInit() 1 8.02e-06 1.094e-05 1.965e-05 0.00% Castro::advance() 1 6.698e-06 8.019e-06 1.514e-05 0.00% Castro::post_init() 1 1.207e-05 1.271e-05 1.337e-05 0.00% Amr::bldFineLevels() 1 7.028e-06 8.758e-06 1.026e-05 0.00% Amr::init() 1 6.146e-06 6.922e-06 7.754e-06 0.00% Amr::initialInit() 1 2.883e-06 3.524e-06 7.455e-06 0.00% Castro::post_timestep() 1 5.673e-06 6.177e-06 7.172e-06 0.00% Other 2521 0.04988 0.06141 0.09604 0.66% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.53 14.53 14.54 100.00% Amr::coarseTimeStep() 1 12.32 12.32 12.32 84.73% Amr::timeStep() 1 12.27 12.27 12.27 84.36% Castro::advance() 1 12.2 12.21 12.21 83.96% Castro::subcycle_advance_ctu() 1 12.13 12.15 12.16 83.62% Castro::do_advance_ctu() 20 12.09 12.11 12.12 83.35% Castro::construct_ctu_hydro_source() 20 6.643 6.673 6.721 46.22% Castro::clean_state() 79 3.668 3.793 3.947 27.14% Castro::computeTemp() 80 2.35 2.469 2.626 18.06% Gravity::solve_phi_with_mlmg() 19 1.608 1.662 1.685 11.59% Castro::construct_new_gravity() 17 1.526 1.579 1.602 11.02% Gravity::solve_for_phi() 17 1.489 1.542 1.566 10.77% Castro::reset_internal_energy(MultiFab) 80 1.502 1.513 1.531 10.53% Castro::initialize_do_advance() 20 1.407 1.443 1.478 10.17% Gravity::actual_solve_with_mlmg() 19 1.409 1.412 1.416 9.74% MLMG::solve() 18 1.353 1.359 1.361 9.36% Amr::InitAmr() 1 1.3 1.303 1.309 9.01% MLMG::oneIter() 128 1.217 1.22 1.225 8.43% MLMG::mgVcycle() 128 1.202 1.205 1.21 8.32% Castro::do_old_sources() 20 1.138 1.172 1.2 8.25% Castro::do_new_sources() 17 0.9242 0.9732 1.019 7.00% FabArray::ParallelCopy() 1429 0.466 0.8005 1.004 6.91% Castro::evaluate_source_change() 196 0.9598 0.9657 0.9813 6.75% Amr::init() 1 0.91 0.91 0.91 6.26% Amr::initialInit() 1 0.8555 0.8555 0.8556 5.88% MLCellLinOp::smooth() 1544 0.6271 0.6859 0.8276 5.69% amrex::Copy() 724 0.7983 0.8036 0.8173 5.62% Castro::enforce_min_density() 79 0.7954 0.8007 0.8151 5.61% FabArray::ParallelCopy_finish() 1429 0.2226 0.5596 0.7654 5.26% AmrLevel::FillPatch() 147 0.54 0.6125 0.6819 4.69% Amr::InitializeInit() 1 0.4435 0.4771 0.6211 4.27% Amr::defBaseLevel() 1 0.4435 0.4771 0.621 4.27% MultiFab::AddProduct() 3439 0.5514 0.5545 0.564 3.88% Castro::normalize_species() 79 0.5475 0.5502 0.5555 3.82% Castro::reset_internal_energy(Fab) 160 0.5203 0.5253 0.5309 3.65% FillPatchIterator::Initialize 147 0.3565 0.4295 0.4996 3.44% MLPoisson::Fsmooth() 6144 0.3544 0.3855 0.4757 3.27% FillPatchIterator::FillFromLevel0() 147 0.3289 0.4088 0.4725 3.25% FillPatchSingleLevel 147 0.3281 0.4081 0.4717 3.24% MLMG::mgVcycle_up::1 128 0.1068 0.2824 0.4603 3.17% MLCellLinOp::applyBC() 7148 0.3091 0.3599 0.4272 2.94% Amr::FinalizeInit() 1 0.2344 0.3784 0.4119 2.83% MLMG::addInterpCorrection() 768 0.06635 0.3229 0.3992 2.75% Castro::Castro() 1 0.3639 0.3714 0.3794 2.61% Castro::initMFs() 1 0.3581 0.3654 0.3732 2.57% FabArray::setVal() 5013 0.348 0.3515 0.359 2.47% MultiFab::Subtract() 159 0.3368 0.3388 0.3438 2.36% FabArray::FillBoundary() 6380 0.239 0.2796 0.3244 2.23% MLMG::mgVcycle_up::2 128 0.01379 0.127 0.2994 2.06% MLMG::mgVcycle_down::0 128 0.2806 0.2841 0.2882 1.98% Castro::expand_state() 20 0.2078 0.2429 0.2796 1.92% Gravity::fill_multipole_BCs() 19 0.1984 0.2481 0.2682 1.84% FabArray::ParallelCopy_nowait() 1429 0.2192 0.2372 0.2605 1.79% Castro::initData() 1 0.07775 0.1041 0.2417 1.66% Amr::bldFineLevels() 1 0.04601 0.1904 0.2233 1.54% MLMG::mgVcycle_up::0 128 0.2108 0.2126 0.2157 1.48% check_for_negative_density() 20 0.08256 0.1477 0.2062 1.42% MLCellLinOp::apply() 986 0.1558 0.1697 0.2055 1.41% FabArray::sum() 3439 0.197 0.1985 0.2011 1.38% MLMG::mgVcycle_up::3 128 0.01371 0.02773 0.1933 1.33% Amr::grid_places() 2 0.002541 0.1465 0.18 1.24% AmrMesh::MakeNewGrids() 2 0.002485 0.1464 0.18 1.24% TagBoxArray::mapPRD 2 0.0008731 0.1111 0.1783 1.23% FillBoundary_nowait() 6380 0.1056 0.1322 0.1675 1.15% TagBoxArray::collate() 2 0.0002818 0.03395 0.1594 1.10% MLMG:computeResOfCorrection() 768 0.1205 0.1304 0.1547 1.06% MLCellLinOp::correctionResidual() 768 0.1189 0.1288 0.153 1.05% FillBoundary_finish() 6380 0.09684 0.133 0.1487 1.02% Castro::post_init() 1 0.1428 0.1442 0.1455 1.00% MLMG::mgVcycle_down::1 128 0.1276 0.1325 0.1376 0.95% Castro::construct_new_source() 102 0.1287 0.1296 0.1311 0.90% MLPoisson::Fapply() 986 0.09491 0.09921 0.1153 0.79% Gravity::multilevel_solve_for_new_phi() 1 0.1089 0.109 0.109 0.75% Gravity::actual_multilevel_solve() 1 0.1089 0.1089 0.109 0.75% MLCellLinOp::solutionResidual() 147 0.09034 0.09642 0.1015 0.70% Amr::writePlotFile() 2 0.0988 0.09889 0.0989 0.68% Castro::create_source_corrector() 17 0.05733 0.07075 0.09278 0.64% Castro::post_regrid() 2 0.08637 0.08764 0.08906 0.61% Castro::construct_old_source() 120 0.08305 0.08364 0.08521 0.59% MLMG::computeResidual() 128 0.07806 0.07918 0.0813 0.56% Castro::initialize_advance() 1 0.04959 0.05427 0.08049 0.55% MLMG::mgVcycle_down::2 128 0.01604 0.04679 0.08023 0.55% Castro::construct_new_gravity_source() 17 0.07243 0.07319 0.07458 0.51% FabArray::Saxpy() 91 0.0722 0.07235 0.07277 0.50% MLMG::mgVcycle_bottom 128 0.00036 0.004871 0.07246 0.50% Castro::derive() 18 0.07046 0.0712 0.07239 0.50% AmrLevel::derive() 18 0.07041 0.07116 0.07235 0.50% MLMG::actualBottomSolve() 128 0.000132 0.00464 0.07221 0.50% MLCGSolver::bicgstab 8 0 0.004397 0.07036 0.48% MLMG::mgVcycle_down::3 128 0.01541 0.02224 0.06792 0.47% Castro::post_timestep() 1 0.05845 0.06091 0.06396 0.44% FabArray::Xpay() 915 0.0515 0.05476 0.06356 0.44% amrex::average_down 768 0.03414 0.04133 0.05911 0.41% Castro::construct_old_gravity() 20 0.04777 0.05384 0.05772 0.40% Gravity::get_old_grav_vector() 20 0.04764 0.05372 0.05759 0.40% Castro::construct_old_gravity_source() 20 0.04914 0.04965 0.05087 0.35% Castro::check_for_nan() 37 0.04592 0.04633 0.04667 0.32% MultiFab::contains_nan() 37 0.04581 0.04621 0.04653 0.32% Amr::regrid() 1 0.04443 0.04506 0.04567 0.31% MLPoisson::define() 19 0.03213 0.03553 0.04151 0.29% MLMG::mgVcycle_down::4 128 0.01532 0.01719 0.04069 0.28% Castro::retry_advance_ctu() 11 0.03611 0.03669 0.03826 0.26% Gravity::get_new_grav_vector() 18 0.03607 0.03641 0.03675 0.25% StateDataPhysBCFunct::() 147 0.0007293 0.01842 0.0346 0.24% StateData::FillBoundary(geom) 143 0 0.01754 0.03361 0.23% MLMG::mgVcycle_up::4 128 0.01379 0.01516 0.03207 0.22% Castro::apply_source_to_state() 37 0.02949 0.02957 0.02974 0.20% MLMG::mgVcycle_down::5 128 0.01534 0.01641 0.02823 0.19% MLCellLinOp::defineAuxData() 19 0.01969 0.02154 0.02717 0.19% VisMF::Write(FabArray) 4 0.01402 0.01517 0.02508 0.17% MLMG::prepareForSolve() 18 0.02057 0.02148 0.02429 0.17% MLMG::mgVcycle_up::5 128 0.01388 0.01478 0.02423 0.17% MLMG::computeMLResidual() 18 0.01126 0.01671 0.02198 0.15% Castro::sum_integrated_quantities() 2 0.0214 0.02155 0.02157 0.15% Castro::save_data_for_retry() 3 0.01927 0.01943 0.02035 0.14% MLMG::ResNormInf() 146 0.01274 0.01632 0.01984 0.14% Castro::problem_post_timestep() 2 0.01022 0.01266 0.01612 0.11% Castro::wd_update() 2 0.01021 0.01264 0.01611 0.11% FabArrayBase::getCPC() 2005 0.01001 0.01147 0.01522 0.10% FabArrayBase::getFB() 6380 0.01089 0.01183 0.01514 0.10% amrex::Add() 128 0.01444 0.01454 0.01476 0.10% FabArray::setDomainBndry() 147 0.0004935 0.007597 0.0145 0.10% Gravity::test_composite_phi() 1 0.01436 0.01442 0.01446 0.10% FabArray::norminf() 244 0.0107 0.01111 0.01406 0.10% MLLinOp::define() 19 0.003917 0.007406 0.01102 0.08% Other 2521 0.08388 0.09958 0.1429 0.98% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Pinned Memory Usage: --------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max --------------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 16 9 B 12 B 20 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- Cpu Memory Usage: ------------------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max ------------------------------------------------------------------------------------------------------------------------------- Castro::save_data_for_retry() 416 416 1570 KiB 1571 KiB 1572 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 413 KiB 415 KiB 418 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1044 KiB 1045 KiB 1045 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 993 KiB 994 KiB 995 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 666 KiB 667 KiB 667 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 4312 B 4392 B 4888 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 50 KiB 50 KiB 51 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 12 KiB 13 KiB 15 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 416 KiB 417 KiB 418 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 2367 B 4159 B 5453 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1343 B 1358 B 1382 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 17 KiB 17 KiB 17 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 132 KiB 132 KiB 132 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 192 129 KiB 129 KiB 129 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 6102 9770 B 10 KiB 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 260 B 260 B 263 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 437 B 438 B 441 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 18 B 21 B 38 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 33 KiB 33 KiB 33 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 51 B 51 B 52 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 4 B 5 B 8 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 112 B 112 B 114 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 149 B 168 B 181 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 64 31 KiB 32 KiB 32 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 2432 230 B 236 B 240 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1077 B 1082 B 1100 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 122 B 122 B 122 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1676 B 1736 B 1764 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 32 0 B 0 B 1 B 16 KiB 16 KiB 16 KiB MLCellLinOp::defineAuxData() 11704 11704 442 B 494 B 685 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 459 B 461 B 462 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 23 B 34 B 46 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 19 B 32 B 2916 B 2916 B 2916 B OwnerMask() 64 64 0 B 0 B 0 B 2592 B 2592 B 2592 B MLMG::addInterpCorrection() 5632 5632 1 B 16 B 25 B 1024 B 1024 B 1024 B amrex::average_down 5632 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 64 0 B 5 B 7 B 648 B 648 B 648 B MLCGSolver::bicgstab 640 640 0 B 0 B 2 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 114375 0 B 0 B 0 B 8 B 8 B 8 B ------------------------------------------------------------------------------------------------------------------------------- AMReX (25.06-22-g1d307c2c661a-dirty) finalized