Initializing AMReX (25.01-16-g92d35c2c8163)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (25.01-16-g92d35c2c8163) initialized Starting run at 09:12:28 UTC on 2025-01-20. Successfully read inputs file ... Castro git describe: 25.01-5-g046cf076f AMReX git describe: 25.01-16-g92d35c2c8 Microphysics git describe: 25.01-29-g66a8b4a8 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: 1.393927591 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.023366229 new finest: 0 Gravity::fill_multipole_BCs() time = 0.021072214 Gravity::multilevel_solve_for_new_phi() time = 0.174834795 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.021304336 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.049856644 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004712155 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.005137851 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.066096421 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.00915726 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.065468092 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.006497154 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.063091851 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.006513498 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.061149238 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.355862749 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.014566497 Gravity::solve_for_phi() time = 0.102783642 on level 0 Castro::construct_new_gravity() time = 0.105142452 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.061723534 on level 0 Castro::construct_old_gravity() time = 0.002984787 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.060379648 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.353417498 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.015865967 Gravity::solve_for_phi() time = 0.097000971 on level 0 Castro::construct_new_gravity() time = 0.099210173 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.060477813 on level 0 Castro::construct_old_gravity() time = 0.002890776 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.060098711 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.348985796 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.015397975 Gravity::solve_for_phi() time = 0.095522397 on level 0 Castro::construct_new_gravity() time = 0.097810579 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.060543702 on level 0 Castro::construct_old_gravity() time = 0.002847927 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.060027534 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.351116521 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.014641495 Gravity::solve_for_phi() time = 0.094000112 on level 0 Castro::construct_new_gravity() time = 0.096212294 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.061519319 on level 0 Castro::construct_old_gravity() time = 0.003257968 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.06027315 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.355639962 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.014591408 Gravity::solve_for_phi() time = 0.094090733 on level 0 Castro::construct_new_gravity() time = 0.096309989 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.061131909 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.00351843 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.060150383 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.349634178 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.014576071 Gravity::solve_for_phi() time = 0.094565545 on level 0 Castro::construct_new_gravity() time = 0.096795154 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.061006102 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.002941373 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.060256067 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.357088228 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.014772849 Gravity::solve_for_phi() time = 0.094279354 on level 0 Castro::construct_new_gravity() time = 0.096497594 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.06305466 on level 0 Castro::construct_old_gravity() time = 0.003890959 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.061353517 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.35560546 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.015061884 Gravity::solve_for_phi() time = 0.095650285 on level 0 Castro::construct_new_gravity() time = 0.097886964 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.061062984 on level 0 Castro::construct_old_gravity() time = 0.003038859 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.060282658 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.3559157 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.014617727 Gravity::solve_for_phi() time = 0.093917786 on level 0 Castro::construct_new_gravity() time = 0.096168012 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.060290837 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.002963361 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.059968885 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.351771007 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.014960653 Gravity::solve_for_phi() time = 0.095804409 on level 0 Castro::construct_new_gravity() time = 0.098012462 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.060274174 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.002957092 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.060010786 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.35709236 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.015333138 Gravity::solve_for_phi() time = 0.095957614 on level 0 Castro::construct_new_gravity() time = 0.098173212 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.061246642 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.002965587 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.060058922 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.355834128 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.01534957 Gravity::solve_for_phi() time = 0.094276375 on level 0 Castro::construct_new_gravity() time = 0.096491543 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.06109468 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.002933406 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.060125507 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.357106894 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.015579028 Gravity::solve_for_phi() time = 0.09568568 on level 0 Castro::construct_new_gravity() time = 0.097895903 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.062250216 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.00422108 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.060126903 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.352037679 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.015845706 Gravity::solve_for_phi() time = 0.096269015 on level 0 Castro::construct_new_gravity() time = 0.098494011 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.061195875 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.002991946 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.060415078 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.360977775 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.015249102 Gravity::solve_for_phi() time = 0.09462968 on level 0 Castro::construct_new_gravity() time = 0.096888913 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.060343777 on level 0 Castro::construct_old_gravity() time = 0.002890327 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.060069133 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.353332171 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.016425092 Gravity::solve_for_phi() time = 0.096049699 on level 0 Castro::construct_new_gravity() time = 0.098257677 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.061263721 on level 0 Castro::construct_old_gravity() time = 0.003108932 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.060179493 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.358309377 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002610369075 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.016117872 Gravity::solve_for_phi() time = 0.097067041 on level 0 Castro::construct_new_gravity() time = 0.099315274 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.061323051 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.61420187 [STEP 1] FAB kilobyte spread across MPI nodes: [6998 ... 6998] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.050069294 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004266315 seconds Ending run at 09:12:45 UTC on 2025-01-20. Run time = 17.1834654 Run time without initialization = 12.66952946 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]: 17.1833 ... 17.1834 ... 17.1835 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.8863 5.9306 5.9937 34.88% Amr::InitAmr() 1 2.059 2.059 2.059 11.98% FabArray::ParallelCopy_finish() 1429 0.4201 1.168 1.941 11.30% TagBoxArray::collate() 2 0.003484 0.5573 1.36 7.92% Castro::initMFs() 1 0.3582 0.4388 1.352 7.87% Castro::computeTemp() 80 0.8742 0.9833 1.122 6.53% amrex::Copy() 724 0.8216 0.8435 0.932 5.42% MultiFab::AddProduct() 3439 0.5869 0.5993 0.6774 3.94% Castro::normalize_species() 79 0.56 0.5763 0.6235 3.63% Castro::reset_internal_energy(Fab) 160 0.5266 0.5437 0.5795 3.37% MLPoisson::Fsmooth() 6144 0.3675 0.4003 0.4762 2.77% FabArray::setVal() 5013 0.3543 0.3672 0.4391 2.56% MultiFab::Subtract() 159 0.3497 0.3603 0.4154 2.42% Gravity::fill_multipole_BCs() 19 0.2273 0.2678 0.2981 1.73% FabArray::ParallelCopy_nowait() 1429 0.2032 0.2255 0.2464 1.43% FillBoundary_finish() 6380 0.1273 0.1784 0.2196 1.28% FabArray::sum() 3439 0.1842 0.1891 0.2146 1.25% MultiFab::max() 2 0.0001749 0.1803 0.213 1.24% check_for_negative_density() 20 0.0794 0.14 0.1889 1.10% FillBoundary_nowait() 6380 0.09516 0.1268 0.166 0.97% Castro::initData() 1 0.003482 0.03007 0.1636 0.95% MLCellLinOp::applyBC() 7148 0.06485 0.08223 0.1283 0.75% MLPoisson::Fapply() 986 0.09717 0.1025 0.1168 0.68% Castro::enforce_min_density() 79 0.01891 0.02424 0.09762 0.57% Castro::construct_new_gravity_source() 17 0.06907 0.07189 0.08244 0.48% FabArray::Saxpy() 264 0.07031 0.07181 0.07922 0.46% FabArray::Xpay() 945 0.05327 0.0558 0.06275 0.37% MultiFab::contains_nan() 37 0.04645 0.04758 0.04905 0.29% Castro::construct_old_gravity_source() 20 0.04652 0.04742 0.04835 0.28% MLMG::addInterpCorrection() 768 0.03482 0.03771 0.04542 0.26% AmrLevel::derive() 18 0.04095 0.04167 0.04288 0.25% amrex::average_down 768 0.02614 0.0287 0.03567 0.21% StateData::FillBoundary(geom) 143 0 0.01828 0.03489 0.20% Castro::reset_internal_energy(MultiFab) 80 0.02458 0.02601 0.02905 0.17% VisMF::Write(FabArray) 4 0.01384 0.0148 0.02449 0.14% MLCellLinOp::defineAuxData() 19 0.01602 0.01719 0.0206 0.12% Castro::evaluate_source_change() 196 0.01458 0.0153 0.01781 0.10% FabArray::FillBoundary() 6380 0.01402 0.01476 0.01771 0.10% MLCGSolver::bicgstab 8 0 0.001103 0.01765 0.10% MLMG::ResNormInf() 146 0.005186 0.01209 0.01738 0.10% Amr::writePlotFile() 2 0.008217 0.01618 0.01732 0.10% amrex::Add() 128 0.01516 0.01541 0.01725 0.10% Castro::construct_new_source() 102 0.01638 0.01663 0.01686 0.10% FabArray::norminf() 244 0.01208 0.01239 0.01532 0.09% FabArray::setDomainBndry() 147 0.0004884 0.00794 0.01513 0.09% MLCellLinOp::smooth() 3072 0.01208 0.01272 0.01459 0.08% Castro::wd_update() 2 0.00853 0.01132 0.0145 0.08% FillPatchIterator::Initialize 147 0.01327 0.01361 0.01422 0.08% FabArrayBase::getFB() 6380 0.01065 0.01116 0.01367 0.08% MLMG::prepareForSolve() 18 0.01047 0.01094 0.01207 0.07% Castro::construct_old_source() 120 0.01046 0.01074 0.01129 0.07% Castro::sum_integrated_quantities() 2 0.0028 0.01003 0.01119 0.07% MLMG::solve() 18 0.00264 0.005859 0.01076 0.06% AmrMesh::MakeNewGrids() 2 0.0006013 0.001374 0.009805 0.06% Gravity::actual_solve_with_mlmg() 19 0.007688 0.008122 0.009368 0.05% FabArrayBase::getCPC() 2005 0.007679 0.007919 0.008252 0.05% Gravity::get_old_grav_vector() 20 0.006676 0.00692 0.007296 0.04% Gravity::get_new_grav_vector() 18 0.006287 0.006386 0.006606 0.04% Castro::initialize_do_advance() 20 0.004555 0.004776 0.004999 0.03% MLCellLinOp::apply() 986 0.002086 0.002327 0.004873 0.03% FabArray::ParallelCopy() 1429 0.003557 0.003751 0.004078 0.02% AmrLevel::FillPatch() 147 0.002482 0.002617 0.002896 0.02% Amr::grid_places() 2 5.877e-05 0.0002008 0.002126 0.01% MLCellLinOp::correctionResidual() 768 0.001932 0.002007 0.002103 0.01% TagBoxArray::mapPRD 2 0.0001267 0.0002829 0.001981 0.01% MLMG::mgVcycle() 128 0.001709 0.001791 0.001913 0.01% Castro::save_data_for_retry() 3 0.001325 0.00138 0.001709 0.01% MLMG:computeResOfCorrection() 768 0.001548 0.001612 0.001706 0.01% Castro::Castro() 1 0.000383 0.0006787 0.001422 0.01% StateDataPhysBCFunct::() 147 0.0007555 0.0009047 0.001127 0.01% MLMG::actualBottomSolve() 128 0.0001382 0.0002019 0.001115 0.01% Castro::subcycle_advance_ctu() 1 0.0002454 0.0003329 0.001069 0.01% main() 1 0.0007214 0.0007731 0.0008861 0.01% FillPatchIterator::FillFromLevel0() 147 0.0007086 0.0007513 0.000856 0.00% Amr::bldFineLevels() 1 8.602e-06 6.459e-05 0.0008446 0.00% MLMG::mgVcycle_down::2 128 0.0006637 0.0007185 0.0008326 0.00% Gravity::solve_for_phi() 17 0.0006958 0.0007384 0.00079 0.00% MLMG::mgVcycle_down::3 128 0.0006238 0.0006734 0.0007683 0.00% MLMG::mgVcycle_down::0 128 0.0007036 0.0007248 0.0007649 0.00% MLMG::mgVcycle_down::1 128 0.0006947 0.0007171 0.0007582 0.00% MLMG::mgVcycle_down::4 128 0.0006062 0.0006379 0.0007221 0.00% MLMG::mgVcycle_down::5 128 0.0006303 0.0006533 0.0007213 0.00% MLMG::mgVcycle_up::0 128 0.0005203 0.000539 0.0005761 0.00% MLMG::mgVcycle_up::1 128 0.0005063 0.0005305 0.0005747 0.00% Castro::clean_state() 79 0.0004659 0.0004943 0.0005595 0.00% MLMG::mgVcycle_up::2 128 0.0004395 0.0004819 0.0005484 0.00% FillPatchSingleLevel 147 0.0004769 0.0005041 0.0005374 0.00% MLMG::mgVcycle_up::3 128 0.0004364 0.0004603 0.000524 0.00% MLCellLinOp::solutionResidual() 147 0.0004624 0.0004797 0.000505 0.00% Castro::do_advance_ctu() 20 0.0004056 0.0004376 0.0004983 0.00% MLMG::mgVcycle_up::4 128 0.0004432 0.0004598 0.0004944 0.00% MLMG::mgVcycle_up::5 128 0.0004358 0.0004525 0.0004891 0.00% MLMG::oneIter() 128 0.0004178 0.0004292 0.0004646 0.00% Castro::do_old_sources() 20 0.0002967 0.0003127 0.0003525 0.00% Gravity::update_max_rhs() 2 0.0002064 0.0002361 0.0003444 0.00% Castro::construct_new_gravity() 17 0.0001761 0.0001923 0.0003116 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002297 0.0002488 0.0002812 0.00% MLMG::computeResidual() 128 0.0002474 0.0002568 0.0002764 0.00% Castro::problem_post_timestep() 2 1.139e-05 3.412e-05 0.0002725 0.00% MLMG::mgVcycle_bottom 128 0.000236 0.0002483 0.0002645 0.00% Castro::do_new_sources() 17 0.0002228 0.0002351 0.0002492 0.00% Amr::defBaseLevel() 1 9.908e-05 0.0001202 0.0001746 0.00% MLPoisson::define() 19 0.0001333 0.0001466 0.0001707 0.00% Castro::check_for_nan() 37 0.0001161 0.0001278 0.0001699 0.00% Amr::coarseTimeStep() 1 0.0001344 0.0001576 0.0001636 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.603e-05 8.785e-05 0.0001613 0.00% Gravity::test_composite_phi() 1 0.000113 0.0001233 0.0001565 0.00% Castro::construct_old_gravity() 20 0.0001072 0.0001191 0.0001481 0.00% Castro::create_source_corrector() 17 5.779e-05 6.966e-05 0.0001451 0.00% Castro::retry_advance_ctu() 11 6.251e-05 7.455e-05 0.0001412 0.00% Castro::apply_source_to_state() 37 0.0001076 0.0001131 0.000127 0.00% Castro::initialize_advance() 1 0.0001056 0.0001099 0.0001236 0.00% Gravity::actual_multilevel_solve() 1 3.391e-05 4.215e-05 6.933e-05 0.00% Castro::expand_state() 20 5.128e-05 5.641e-05 6.399e-05 0.00% Castro::derive() 18 4.179e-05 4.593e-05 5.82e-05 0.00% Amr::FinalizeInit() 1 1.045e-05 1.806e-05 5.182e-05 0.00% MLMG::computeMLResidual() 18 4.281e-05 4.571e-05 5.039e-05 0.00% Amr::regrid() 1 2.002e-05 3.501e-05 4.886e-05 0.00% Amr::timeStep() 1 1.508e-05 1.793e-05 2.557e-05 0.00% Castro::post_regrid() 2 1.244e-05 1.633e-05 2.546e-05 0.00% Castro::advance() 1 7.196e-06 9.734e-06 1.934e-05 0.00% Castro::post_init() 1 1.385e-05 1.509e-05 1.932e-05 0.00% Amr::InitializeInit() 1 7.913e-06 9.394e-06 1.558e-05 0.00% Amr::init() 1 6.693e-06 7.843e-06 1.055e-05 0.00% Amr::initialInit() 1 3.208e-06 4.143e-06 9.687e-06 0.00% Castro::post_timestep() 1 6.043e-06 6.594e-06 7.333e-06 0.00% Other 2436 0.05069 0.06285 0.09819 0.57% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 17.18 17.18 17.18 100.00% Amr::coarseTimeStep() 1 12.67 12.67 12.67 73.73% Amr::timeStep() 1 12.61 12.61 12.61 73.41% Castro::advance() 1 12.55 12.55 12.55 73.05% Castro::subcycle_advance_ctu() 1 12.5 12.5 12.5 72.76% Castro::do_advance_ctu() 20 12.46 12.46 12.46 72.53% Castro::construct_ctu_hydro_source() 20 6.85 6.876 6.92 40.27% Castro::clean_state() 79 3.801 3.955 4.308 25.07% Castro::computeTemp() 80 2.45 2.569 2.874 16.73% Amr::init() 1 2.455 2.455 2.455 14.29% Amr::initialInit() 1 2.4 2.4 2.4 13.97% FabArray::ParallelCopy() 1429 0.6925 1.419 2.191 12.75% Amr::InitAmr() 1 2.059 2.059 2.059 11.98% Amr::FinalizeInit() 1 0.5648 1.823 1.95 11.35% FabArray::ParallelCopy_finish() 1429 0.4201 1.168 1.941 11.30% Amr::InitializeInit() 1 0.4493 0.5769 1.835 10.68% Amr::defBaseLevel() 1 0.4493 0.5769 1.835 10.68% Gravity::solve_phi_with_mlmg() 19 1.75 1.793 1.824 10.61% Castro::reset_internal_energy(MultiFab) 80 1.543 1.585 1.796 10.45% Castro::construct_new_gravity() 17 1.591 1.634 1.665 9.69% Gravity::solve_for_phi() 17 1.553 1.596 1.627 9.47% Amr::bldFineLevels() 1 0.2711 1.36 1.532 8.92% Gravity::actual_solve_with_mlmg() 19 1.518 1.522 1.526 8.88% Castro::initialize_do_advance() 20 1.426 1.462 1.496 8.71% MLMG::solve() 18 1.459 1.465 1.467 8.54% Amr::grid_places() 2 0.03106 1.289 1.417 8.24% AmrMesh::MakeNewGrids() 2 0.03091 1.289 1.417 8.24% TagBoxArray::mapPRD 2 0.005884 0.7291 1.411 8.21% Castro::Castro() 1 0.3638 0.4464 1.386 8.07% TagBoxArray::collate() 2 0.003484 0.5573 1.36 7.92% Castro::initMFs() 1 0.3584 0.439 1.353 7.88% MLMG::oneIter() 128 1.304 1.309 1.314 7.65% MLMG::mgVcycle() 128 1.289 1.293 1.297 7.55% Castro::do_old_sources() 20 1.176 1.202 1.217 7.08% Castro::evaluate_source_change() 196 0.9881 1.011 1.144 6.66% Castro::do_new_sources() 17 0.9641 0.9982 1.035 6.03% Castro::enforce_min_density() 79 0.8194 0.8453 0.9484 5.52% amrex::Copy() 724 0.8216 0.8435 0.932 5.42% MLCellLinOp::smooth() 3072 0.6852 0.7452 0.8902 5.18% AmrLevel::FillPatch() 147 0.5402 0.6062 0.6918 4.03% MultiFab::AddProduct() 3439 0.5869 0.5993 0.6774 3.94% Castro::normalize_species() 79 0.56 0.5763 0.6235 3.63% Castro::reset_internal_energy(Fab) 160 0.5266 0.5437 0.5795 3.37% FillPatchIterator::Initialize 147 0.3484 0.4156 0.5045 2.94% MLCellLinOp::applyBC() 7148 0.3585 0.4141 0.5042 2.93% MLMG::mgVcycle_up::1 128 0.1159 0.3027 0.4907 2.86% FillPatchIterator::FillFromLevel0() 147 0.3198 0.394 0.4815 2.80% FillPatchSingleLevel 147 0.3191 0.3933 0.4808 2.80% MLPoisson::Fsmooth() 6144 0.3675 0.4003 0.4762 2.77% Castro::initData() 1 0.08 0.1287 0.4468 2.60% Castro::post_init() 1 0.2325 0.4127 0.4454 2.59% FabArray::setVal() 5013 0.3543 0.3672 0.4391 2.56% MLMG::addInterpCorrection() 768 0.06977 0.3395 0.4224 2.46% MultiFab::Subtract() 159 0.3497 0.3603 0.4154 2.42% FabArray::FillBoundary() 6380 0.2839 0.3319 0.3778 2.20% MLMG::mgVcycle_up::2 128 0.01413 0.1327 0.315 1.83% MLMG::mgVcycle_down::0 128 0.3 0.3035 0.311 1.81% Castro::post_regrid() 2 0.0883 0.1209 0.3011 1.75% Gravity::fill_multipole_BCs() 19 0.2295 0.27 0.3003 1.75% FabArray::ParallelCopy_nowait() 1429 0.2233 0.2466 0.2685 1.56% Amr::regrid() 1 0.06726 0.09306 0.2634 1.53% Castro::expand_state() 20 0.2028 0.2284 0.2618 1.52% MLMG::mgVcycle_up::0 128 0.2303 0.2327 0.2348 1.37% MLCellLinOp::apply() 986 0.1735 0.1846 0.2234 1.30% FillBoundary_finish() 6380 0.1273 0.1784 0.2196 1.28% FabArray::sum() 3439 0.1842 0.1891 0.2146 1.25% Gravity::update_max_rhs() 2 0.001591 0.1818 0.2144 1.25% MultiFab::max() 2 0.0001749 0.1803 0.213 1.24% MLMG::mgVcycle_up::3 128 0.01412 0.02863 0.1993 1.16% check_for_negative_density() 20 0.0794 0.14 0.1889 1.10% FillBoundary_nowait() 6380 0.1063 0.1387 0.1786 1.04% Gravity::multilevel_solve_for_new_phi() 1 0.1743 0.1744 0.1748 1.02% Gravity::actual_multilevel_solve() 1 0.1742 0.1743 0.1747 1.02% MLMG:computeResOfCorrection() 768 0.1249 0.1373 0.1624 0.95% MLCellLinOp::correctionResidual() 768 0.1233 0.1356 0.1608 0.94% MLMG::mgVcycle_down::1 128 0.141 0.1481 0.1548 0.90% Castro::construct_new_source() 102 0.1253 0.1288 0.1397 0.81% MLPoisson::Fapply() 986 0.09717 0.1025 0.1168 0.68% MLCellLinOp::solutionResidual() 147 0.09866 0.1054 0.1111 0.65% Amr::writePlotFile() 2 0.1003 0.1003 0.1004 0.58% MLMG::computeResidual() 128 0.08264 0.08405 0.08576 0.50% Castro::construct_old_source() 120 0.08037 0.08209 0.0846 0.49% MLMG::mgVcycle_down::2 128 0.01644 0.0493 0.08383 0.49% Castro::construct_new_gravity_source() 17 0.06907 0.07189 0.08244 0.48% FabArray::Saxpy() 264 0.07031 0.07181 0.07922 0.46% Castro::create_source_corrector() 17 0.05477 0.06657 0.0792 0.46% Castro::derive() 18 0.0712 0.07253 0.0773 0.45% AmrLevel::derive() 18 0.07116 0.07249 0.07724 0.45% MLMG::mgVcycle_bottom 128 0.0003742 0.005125 0.07617 0.44% MLMG::actualBottomSolve() 128 0.0001382 0.004877 0.07591 0.44% MLCGSolver::bicgstab 8 0 0.004645 0.07431 0.43% MLMG::mgVcycle_down::3 128 0.01572 0.02311 0.07182 0.42% Castro::construct_old_gravity() 20 0.04911 0.06143 0.07081 0.41% Gravity::get_old_grav_vector() 20 0.04899 0.06131 0.07069 0.41% Castro::post_timestep() 1 0.06114 0.06355 0.06549 0.38% amrex::average_down 768 0.03535 0.04358 0.06528 0.38% FabArray::Xpay() 945 0.05327 0.0558 0.06275 0.37% Castro::initialize_advance() 1 0.05008 0.05074 0.05216 0.30% Castro::check_for_nan() 37 0.04657 0.04771 0.04922 0.29% MultiFab::contains_nan() 37 0.04645 0.04758 0.04905 0.29% Castro::construct_old_gravity_source() 20 0.04652 0.04742 0.04835 0.28% MLPoisson::define() 19 0.03212 0.03799 0.04387 0.26% MLMG::mgVcycle_down::4 128 0.01567 0.01766 0.04128 0.24% Castro::retry_advance_ctu() 11 0.03697 0.03816 0.04115 0.24% Gravity::get_new_grav_vector() 18 0.03751 0.03861 0.03949 0.23% StateDataPhysBCFunct::() 147 0.0007555 0.01918 0.03583 0.21% StateData::FillBoundary(geom) 143 0 0.01828 0.03489 0.20% MLMG::mgVcycle_up::4 128 0.01413 0.01561 0.03274 0.19% Castro::sum_integrated_quantities() 2 0.03091 0.03108 0.0311 0.18% Castro::apply_source_to_state() 37 0.0287 0.02914 0.03081 0.18% MLMG::mgVcycle_down::5 128 0.01565 0.01681 0.02886 0.17% MLCellLinOp::defineAuxData() 19 0.02022 0.02235 0.02785 0.16% MLMG::ResNormInf() 146 0.01588 0.02223 0.0274 0.16% Gravity::test_composite_phi() 1 0.02616 0.02662 0.02671 0.16% MLMG::computeMLResidual() 18 0.01293 0.02008 0.02519 0.15% MLMG::prepareForSolve() 18 0.02107 0.02218 0.02484 0.14% MLMG::mgVcycle_up::5 128 0.01427 0.01522 0.02477 0.14% VisMF::Write(FabArray) 4 0.01384 0.0148 0.02449 0.14% Castro::save_data_for_retry() 3 0.01969 0.02042 0.02389 0.14% Castro::problem_post_timestep() 2 0.01273 0.01529 0.018 0.10% Castro::wd_update() 2 0.01246 0.01525 0.01797 0.10% amrex::Add() 128 0.01516 0.01541 0.01725 0.10% FabArrayBase::getCPC() 2005 0.01006 0.01166 0.01541 0.09% FabArrayBase::getFB() 6380 0.01101 0.01192 0.01534 0.09% FabArray::norminf() 244 0.01208 0.01239 0.01532 0.09% FabArray::setDomainBndry() 147 0.0004884 0.00794 0.01513 0.09% Other 2436 0.08891 0.1097 0.1591 0.93% -------------------------------------------------------------------------------------------- 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 8 B 9 B 11 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 1369 KiB 1370 KiB 1370 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 360 KiB 362 KiB 366 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1010 KiB 1010 KiB 1011 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 961 KiB 962 KiB 962 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 580 KiB 580 KiB 580 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 3696 B 3749 B 4139 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 43 KiB 45 KiB 51 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 10 KiB 11 KiB 12 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 363 KiB 363 KiB 363 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 1795 B 3221 B 4158 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1146 B 1160 B 1185 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 15 KiB 15 KiB 17 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 115 KiB 115 KiB 115 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 192 125 KiB 125 KiB 125 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 6102 8916 B 9437 B 10 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 219 B 224 B 238 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 368 B 373 B 379 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 16 B 17 B 33 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 31 KiB 35 KiB 35 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 80 B 81 B 82 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 3 B 174 B 205 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 98 B 101 B 104 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 130 B 162 B 188 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 64 28 KiB 28 KiB 28 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 2432 204 B 213 B 234 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 937 B 958 B 1085 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 166 B 166 B 166 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1479 B 1520 B 1550 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 32 1 B 1 B 1 B 16 KiB 16 KiB 16 KiB MLCellLinOp::defineAuxData() 11704 11704 403 B 450 B 624 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 419 B 420 B 421 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 25 B 36 B 49 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 1 B 109 B 212 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 14 B 23 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 48 B 53 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.01-16-g92d35c2c8163) finalized