Initializing AMReX (26.01-21-ge5851400bcf7)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.01-21-ge5851400bcf7) initialized Starting run at 09:16:04 UTC on 2026-01-20. Successfully read inputs file ... Castro git describe: 26.01-11-g3d7eded2f AMReX git describe: 26.01-21-ge5851400b Microphysics git describe: 26.01-17-gfd03c5b2 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.16455602 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.003214192 new finest: 0 Gravity::fill_multipole_BCs() time = 0.033641868 Gravity::multilevel_solve_for_new_phi() time = 0.330422987 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.033411958 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.109712383 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.007501296 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 = -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 = -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 = -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.009520441 on level 0 Castro::do_old_sources() time = 0.026337658 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.147621549 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.009133741 on level 0 Castro::do_old_sources() time = 0.026450229 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.146758656 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.020699847 on level 0 Castro::do_old_sources() time = 0.031191161 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.162083982 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.011940555 on level 0 Castro::do_old_sources() time = 0.026838738 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.146561047 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.817117892 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.039793553 Gravity::solve_for_phi() time = 0.275125814 on level 0 Castro::construct_new_gravity() time = 0.282725809 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.150659416 on level 0 Castro::construct_old_gravity() time = 0.009845273 on level 0 Castro::do_old_sources() time = 0.026962363 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.146653083 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.816879311 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.039890946 Gravity::solve_for_phi() time = 0.277914238 on level 0 Castro::construct_new_gravity() time = 0.285487087 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.150500592 on level 0 Castro::construct_old_gravity() time = 0.009751286 on level 0 Castro::do_old_sources() time = 0.027008435 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.146662216 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.838149337 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.039833981 Gravity::solve_for_phi() time = 0.273792076 on level 0 Castro::construct_new_gravity() time = 0.28138547 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.150676168 on level 0 Castro::construct_old_gravity() time = 0.00990123 on level 0 Castro::do_old_sources() time = 0.026780857 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.146523003 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.828133887 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.039989503 Gravity::solve_for_phi() time = 0.273892961 on level 0 Castro::construct_new_gravity() time = 0.281410818 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.150555102 on level 0 Castro::construct_old_gravity() time = 0.009826033 on level 0 Castro::do_old_sources() time = 0.026951998 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.146873077 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.84617333 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.040305616 Gravity::solve_for_phi() time = 0.277185343 on level 0 Castro::construct_new_gravity() time = 0.284756133 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.150463347 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.009822779 on level 0 Castro::do_old_sources() time = 0.027005316 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.146586658 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.827809414 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.040143488 Gravity::solve_for_phi() time = 0.273875458 on level 0 Castro::construct_new_gravity() time = 0.281378483 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.150243236 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.009735746 on level 0 Castro::do_old_sources() time = 0.026935177 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.146430021 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.845322356 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.040130476 Gravity::solve_for_phi() time = 0.275233924 on level 0 Castro::construct_new_gravity() time = 0.28287221 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.150460924 on level 0 Castro::construct_old_gravity() time = 0.009788875 on level 0 Castro::do_old_sources() time = 0.026898589 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.146433358 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.827894685 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.040807015 Gravity::solve_for_phi() time = 0.276356549 on level 0 Castro::construct_new_gravity() time = 0.283927518 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.150441752 on level 0 Castro::construct_old_gravity() time = 0.009926488 on level 0 Castro::do_old_sources() time = 0.026971916 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.1465858 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.905730095 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.039664946 Gravity::solve_for_phi() time = 0.274219647 on level 0 Castro::construct_new_gravity() time = 0.281860392 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.150833006 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.00985662 on level 0 Castro::do_old_sources() time = 0.027003386 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.146566219 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.825628175 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.040458318 Gravity::solve_for_phi() time = 0.287510068 on level 0 Castro::construct_new_gravity() time = 0.29499922 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.150230901 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.00979717 on level 0 Castro::do_old_sources() time = 0.026996003 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.147439552 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.868779343 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.044990692 Gravity::solve_for_phi() time = 0.278971982 on level 0 Castro::construct_new_gravity() time = 0.286598805 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.150399033 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.009854577 on level 0 Castro::do_old_sources() time = 0.026875477 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.14632862 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.854123193 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.038431131 Gravity::solve_for_phi() time = 0.272865076 on level 0 Castro::construct_new_gravity() time = 0.280357158 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.151036549 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.009787907 on level 0 Castro::do_old_sources() time = 0.026783732 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.146572391 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.860481361 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.038678918 Gravity::solve_for_phi() time = 0.269444033 on level 0 Castro::construct_new_gravity() time = 0.276994771 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.148160614 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.014644571 on level 0 Castro::do_old_sources() time = 0.027245973 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.145441683 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.831864478 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.051305601 Gravity::solve_for_phi() time = 0.285272865 on level 0 Castro::construct_new_gravity() time = 0.292885728 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.148267004 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.009490635 on level 0 Castro::do_old_sources() time = 0.026343906 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.143741099 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.892673332 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.039478724 Gravity::solve_for_phi() time = 0.278691796 on level 0 Castro::construct_new_gravity() time = 0.286190634 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.15005096 on level 0 Castro::construct_old_gravity() time = 0.009597665 on level 0 Castro::do_old_sources() time = 0.026632357 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.145822112 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.821445871 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.039033495 Gravity::solve_for_phi() time = 0.269580944 on level 0 Castro::construct_new_gravity() time = 0.277205564 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.147212034 on level 0 Castro::construct_old_gravity() time = 0.01078103 on level 0 Castro::do_old_sources() time = 0.028559736 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.144100679 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.839485114 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.001016887975 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.040166229 Gravity::solve_for_phi() time = 0.274213812 on level 0 Castro::construct_new_gravity() time = 0.281707848 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.148728352 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.581215356e+29 ener added: 1.581215356e+29 TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -4.341809003e+38 TIME= 0.25 YMOM = 1.137357411e+28 TIME= 0.25 ZMOM = 0 TIME= 0.25 ANG MOM X = 0 TIME= 0.25 ANG MOM Y = 0 TIME= 0.25 ANG MOM Z = -1.469316165e+38 TIME= 0.25 RHO*e = 1.325894512e+50 TIME= 0.25 RHO*K = 3.170941987e+49 TIME= 0.25 RHO*E = 1.642988711e+50 TIME= 0.25 RHO*PHI = -5.449011378e+50 TIME= 0.25 TOTAL ENERGY= -1.081516979e+50 TIME= 0.25 CENTER OF MASS X-LOC = 316424388.8 TIME= 0.25 CENTER OF MASS X-VEL = -170267.2489 TIME= 0.25 CENTER OF MASS Y-LOC = 2.427336677e-06 TIME= 0.25 CENTER OF MASS Y-VEL = 4.460231145e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 57600050.72 TIME= 0.25 MAXIMUM DENSITY = 4458760.788 TIME= 0.25 MAXIMUM T_S / T_E = 0 STEP = 1 TIME = 0.25 DT = 0.25 [STEP 1] Coarse TimeStep time: 32.37924383 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.109996228 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.006815378 seconds Ending run at 09:16:39 UTC on 2026-01-20. Run time = 35.36685185 Run time without initialization = 32.49697986 Average number of zones advanced per microsecond: 0.001 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 35.3588 ... 35.3616 ... 35.3669 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 13.8310 14.0187 14.2313 40.24% amrex::Copy() 766 3.255 3.288 3.324 9.40% MultiFab::AddProduct() 3598 2.219 2.243 2.267 6.41% MLPoisson::Fsmooth() 6144 1.509 1.634 1.938 5.48% FabArray::ParallelCopy_finish() 1469 0.5373 1.203 1.6 4.52% FabArray::setVal() 5232 1.362 1.381 1.402 3.96% MultiFab::Subtract() 159 1.351 1.368 1.394 3.94% Amr::InitAmr() 1 1.336 1.338 1.344 3.80% Castro::computeTemp() 80 0.9621 1.064 1.18 3.34% FabArray::ParallelCopy_nowait() 1469 0.7336 0.7968 0.8684 2.46% check_for_negative_density() 20 0.2969 0.5525 0.7743 2.19% FabArray::sum() 3598 0.745 0.7529 0.7612 2.15% Gravity::fill_multipole_BCs() 19 0.6488 0.7063 0.7448 2.11% MLPoisson::Fapply() 986 0.5583 0.5845 0.6483 1.83% Castro::do_advance_ctu() 20 0.5169 0.5248 0.5382 1.52% Castro::initMFs() 1 0.3886 0.3917 0.418 1.18% Castro::reset_internal_energy(Fab) 160 0.368 0.3724 0.3775 1.07% FabArray::Saxpy() 111 0.3212 0.3267 0.3519 0.99% Castro::construct_old_gravity_source() 40 0.3214 0.3258 0.333 0.94% FillBoundary_finish() 6380 0.214 0.2588 0.3315 0.94% MLCellLinOp::applyBC() 7148 0.1256 0.175 0.2936 0.83% Castro::construct_new_gravity_source() 17 0.2623 0.2655 0.269 0.76% FillBoundary_nowait() 6380 0.1353 0.1844 0.2398 0.68% FabArray::Xpay() 915 0.1911 0.1998 0.218 0.62% Castro::initData() 1 0.0102 0.03471 0.1698 0.48% MultiFab::contains_nan() 37 0.1643 0.1664 0.1686 0.48% Castro::normalize_species() 79 0.1384 0.14 0.1418 0.40% TagBoxArray::collate() 2 0.0006933 0.02473 0.1257 0.36% MLMG::addInterpCorrection() 768 0.08301 0.08965 0.106 0.30% StateData::FillBoundary(geom) 193 0 0.05085 0.09965 0.28% amrex::average_down 768 0.05923 0.06443 0.07726 0.22% Castro::construct_old_source() 240 0.07088 0.07194 0.07327 0.21% AmrLevel::derive() 18 0.05916 0.05951 0.06078 0.17% Castro::construct_new_source() 102 0.05905 0.05982 0.06072 0.17% amrex::Add() 128 0.05283 0.05334 0.05386 0.15% VisMF::Write(FabArray) 4 0.03661 0.03791 0.05053 0.14% FabArray::setDomainBndry() 187 0.0006585 0.02508 0.04844 0.14% FabArray::norminf() 244 0.03613 0.03678 0.04096 0.12% Castro::enforce_min_density() 79 0.03336 0.03368 0.03423 0.10% Gravity::get_old_grav_vector() 20 0.03025 0.03092 0.0335 0.09% Castro::volWgtSum() 46 0.0283 0.0284 0.02859 0.08% Gravity::get_new_grav_vector() 18 0.02736 0.02769 0.02801 0.08% FabArray::mult() 59 0.0262 0.0267 0.02756 0.08% Castro::reset_internal_energy(MultiFab) 80 0.02445 0.02488 0.02576 0.07% MLCellLinOp::defineAuxData() 19 0.01725 0.01833 0.02221 0.06% Castro::wd_update() 2 0.01028 0.01658 0.02077 0.06% FabArray::FillBoundary() 6380 0.01688 0.01752 0.0204 0.06% MLMG::ResNormInf() 146 0.004983 0.01066 0.01802 0.05% Amr::writePlotFile() 2 0.006878 0.01668 0.01788 0.05% Castro::evaluate_source_change() 196 0.01492 0.01542 0.01676 0.05% MLCGSolver::bicgstab 8 0 0.001044 0.01671 0.05% FillPatchIterator::Initialize 187 0.0154 0.01572 0.01605 0.05% MLMG::prepareForSolve() 18 0.01294 0.01336 0.01426 0.04% MLCellLinOp::smooth() 1544 0.01068 0.01126 0.01273 0.04% MLMG::solve() 18 0.002039 0.00506 0.01175 0.03% Gravity::actual_solve_with_mlmg() 19 0.007406 0.007856 0.009429 0.03% MLCellLinOp::apply() 986 0.002508 0.002779 0.005859 0.02% Castro::initialize_do_advance() 20 0.004748 0.005035 0.00551 0.02% FabArray::ParallelCopy() 1469 0.004156 0.004335 0.004593 0.01% Castro::sum_integrated_quantities() 2 0.003304 0.003788 0.004057 0.01% AmrLevel::FillPatch() 187 0.003188 0.003316 0.003668 0.01% MLCellLinOp::correctionResidual() 768 0.002287 0.002358 0.002435 0.01% MLMG::mgVcycle() 128 0.002146 0.002251 0.002347 0.01% AmrMesh::MakeNewGrids() 2 0.001707 0.001758 0.002056 0.01% MLMG:computeResOfCorrection() 768 0.001863 0.001918 0.002018 0.01% Castro::save_data_for_retry() 3 0.001687 0.001723 0.001854 0.01% StateDataPhysBCFunct::() 187 0.001031 0.001218 0.001438 0.00% MLMG::actualBottomSolve() 128 0.0001826 0.0002577 0.001277 0.00% Castro::subcycle_advance_ctu() 1 0.0002457 0.0003374 0.00114 0.00% FillPatchIterator::FillFromLevel0() 187 0.0009017 0.0009323 0.0009848 0.00% main() 1 0.0006855 0.0007714 0.0008916 0.00% Castro::Castro() 1 0.0004597 0.0005319 0.0008772 0.00% Gravity::solve_for_phi() 17 0.0006922 0.0007151 0.0007726 0.00% FillPatchSingleLevel 187 0.0006568 0.0006858 0.0007228 0.00% MLMG::mgVcycle_down::0 128 0.0006538 0.0006712 0.0007174 0.00% MLMG::mgVcycle_down::1 128 0.00062 0.0006364 0.0007115 0.00% MLMG::mgVcycle_down::2 128 0.0005867 0.0006334 0.0006702 0.00% MLMG::mgVcycle_down::5 128 0.0005801 0.0006001 0.0006617 0.00% MLMG::mgVcycle_down::3 128 0.0005576 0.0005926 0.0006518 0.00% MLMG::mgVcycle_down::4 128 0.0005515 0.0005725 0.0006224 0.00% MLCellLinOp::solutionResidual() 147 0.0005355 0.0005501 0.0005771 0.00% Castro::clean_state() 79 0.0004974 0.0005132 0.0005406 0.00% MLMG::oneIter() 128 0.0004508 0.0004705 0.0005243 0.00% Castro::do_old_sources() 40 0.0004562 0.0004788 0.0005202 0.00% MLMG::mgVcycle_up::0 128 0.0004141 0.0004318 0.0004688 0.00% MLMG::mgVcycle_up::1 128 0.0004138 0.0004322 0.0004607 0.00% MLMG::mgVcycle_up::2 128 0.0003832 0.0004125 0.0004566 0.00% MLMG::mgVcycle_up::3 128 0.0003725 0.0003866 0.0004423 0.00% MLMG::mgVcycle_up::4 128 0.0003782 0.0003919 0.0004374 0.00% MLMG::mgVcycle_up::5 128 0.0003667 0.0003796 0.000401 0.00% MLMG::computeResidual() 128 0.0002987 0.0003107 0.0003312 0.00% Castro::construct_new_gravity() 17 0.0001785 0.0001981 0.0003117 0.00% MLMG::mgVcycle_bottom 128 0.0002774 0.0002838 0.0002993 0.00% TagBoxArray::mapPRD 2 0.000267 0.0002852 0.0002936 0.00% Gravity::solve_phi_with_mlmg() 19 0.00022 0.0002383 0.0002717 0.00% Castro::problem_post_timestep() 2 1.261e-05 2.963e-05 0.0002553 0.00% Castro::do_new_sources() 17 0.0002182 0.0002268 0.0002407 0.00% Amr::coarseTimeStep() 1 0.0001292 0.0001507 0.000157 0.00% MLPoisson::define() 19 0.0001279 0.0001402 0.0001533 0.00% Castro::construct_old_gravity() 20 0.0001206 0.0001316 0.0001495 0.00% Castro::check_for_nan() 37 0.0001259 0.0001324 0.0001484 0.00% Castro::initialize_advance() 1 0.0001117 0.0001247 0.0001469 0.00% Castro::apply_source_to_state() 37 0.0001209 0.000129 0.0001423 0.00% Gravity::test_composite_phi() 1 0.0001053 0.0001104 0.0001412 0.00% Castro::retry_advance_ctu() 11 8.315e-05 9.916e-05 0.0001353 0.00% Amr::defBaseLevel() 1 9.432e-05 0.0001042 0.0001169 0.00% Amr::FinalizeInit() 1 1.275e-05 2.229e-05 8.812e-05 0.00% Amr::grid_places() 2 5.245e-05 6.359e-05 7.832e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.759e-05 6.087e-05 7.146e-05 0.00% MLMG::computeMLResidual() 18 5.26e-05 5.718e-05 6.579e-05 0.00% Castro::expand_state() 20 5.576e-05 6.035e-05 6.541e-05 0.00% Castro::create_source_corrector() 17 5.436e-05 5.819e-05 6.246e-05 0.00% Castro::derive() 18 4.801e-05 4.918e-05 5.328e-05 0.00% Amr::regrid() 1 1.673e-05 2.088e-05 4.68e-05 0.00% Gravity::actual_multilevel_solve() 1 3.213e-05 3.687e-05 4.323e-05 0.00% Castro::post_init() 1 1.696e-05 1.993e-05 2.864e-05 0.00% Amr::timeStep() 1 1.554e-05 1.855e-05 2.55e-05 0.00% Amr::InitializeInit() 1 8.868e-06 1.087e-05 1.638e-05 0.00% Castro::post_regrid() 2 1.103e-05 1.233e-05 1.56e-05 0.00% Castro::advance() 1 8.208e-06 9.775e-06 1.368e-05 0.00% Amr::bldFineLevels() 1 7.661e-06 8.409e-06 1.013e-05 0.00% Amr::init() 1 6.609e-06 7.371e-06 7.933e-06 0.00% Castro::post_timestep() 1 6.304e-06 6.744e-06 7.382e-06 0.00% Amr::initialInit() 1 3.353e-06 3.552e-06 3.87e-06 0.00% Other 10888 0.0974 0.1226 0.1834 0.52% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 35.36 35.36 35.37 100.00% Amr::coarseTimeStep() 1 32.5 32.5 32.5 91.88% Amr::timeStep() 1 32.38 32.38 32.38 91.55% Castro::advance() 1 32.22 32.22 32.23 91.12% Castro::subcycle_advance_ctu() 1 32.09 32.1 32.1 90.77% Castro::do_advance_ctu() 20 31.95 31.96 31.96 90.37% Castro::construct_ctu_hydro_source() 20 16.39 16.44 16.47 46.57% Castro::clean_state() 79 8.329 8.502 8.67 24.52% Castro::computeTemp() 80 5.143 5.289 5.424 15.34% Gravity::solve_phi_with_mlmg() 19 4.957 5.019 5.049 14.28% Castro::construct_new_gravity() 17 4.726 4.786 4.815 13.62% Gravity::solve_for_phi() 17 4.595 4.657 4.687 13.25% Gravity::actual_solve_with_mlmg() 19 4.288 4.3 4.313 12.19% Castro::reset_internal_energy(MultiFab) 80 4.181 4.225 4.269 12.07% MLMG::solve() 18 4.205 4.211 4.213 11.91% Castro::evaluate_source_change() 196 3.724 3.764 3.804 10.75% MLMG::oneIter() 128 3.699 3.703 3.708 10.48% MLMG::mgVcycle() 128 3.645 3.649 3.654 10.33% Castro::do_old_sources() 40 3.375 3.422 3.469 9.81% amrex::Copy() 766 3.255 3.288 3.324 9.40% Castro::initialize_do_advance() 20 3.245 3.272 3.306 9.35% Castro::enforce_min_density() 79 3.101 3.135 3.181 8.99% MLCellLinOp::smooth() 1544 2.013 2.168 2.619 7.41% Castro::do_new_sources() 17 2.428 2.493 2.545 7.20% FabArray::ParallelCopy() 1469 1.37 2.064 2.46 6.96% MultiFab::AddProduct() 3598 2.219 2.243 2.267 6.41% AmrLevel::FillPatch() 187 1.989 2.085 2.209 6.25% MLPoisson::Fsmooth() 6144 1.509 1.634 1.938 5.48% FabArray::ParallelCopy_finish() 1469 0.5373 1.203 1.6 4.52% Amr::init() 1 1.526 1.526 1.526 4.32% FillPatchIterator::Initialize 187 1.179 1.282 1.414 4.00% Amr::initialInit() 1 1.408 1.408 1.408 3.98% FabArray::setVal() 5232 1.362 1.381 1.402 3.96% MultiFab::Subtract() 159 1.351 1.368 1.394 3.94% FillPatchIterator::FillFromLevel0() 187 1.131 1.241 1.366 3.86% FillPatchSingleLevel 187 1.13 1.24 1.365 3.86% Amr::InitAmr() 1 1.336 1.338 1.344 3.80% MLMG::mgVcycle_up::1 128 0.2801 0.7504 1.225 3.46% MLMG::mgVcycle_down::0 128 1.12 1.127 1.132 3.20% MLMG::addInterpCorrection() 768 0.1432 0.8332 1.047 2.96% FabArray::ParallelCopy_nowait() 1469 0.7929 0.8569 0.9296 2.63% MLCellLinOp::applyBC() 7148 0.5808 0.6511 0.8347 2.36% MLCellLinOp::apply() 986 0.6788 0.7113 0.8157 2.31% Amr::FinalizeInit() 1 0.6464 0.7797 0.808 2.28% MLMG::mgVcycle_up::0 128 0.7854 0.7883 0.7908 2.24% check_for_negative_density() 20 0.2969 0.5525 0.7743 2.19% Castro::expand_state() 20 0.6847 0.7288 0.7725 2.18% MLMG::mgVcycle_up::2 128 0.01551 0.3132 0.7638 2.16% Amr::InitializeInit() 1 0.6003 0.6287 0.762 2.15% Amr::defBaseLevel() 1 0.6003 0.6287 0.762 2.15% FabArray::sum() 3598 0.745 0.7529 0.7612 2.15% Gravity::fill_multipole_BCs() 19 0.6568 0.7143 0.7527 2.13% MLPoisson::Fapply() 986 0.5583 0.5845 0.6483 1.83% MLMG:computeResOfCorrection() 768 0.4673 0.5022 0.5953 1.68% MLCellLinOp::correctionResidual() 768 0.4654 0.5003 0.5933 1.68% Castro::construct_old_source() 240 0.5652 0.5733 0.5879 1.66% FabArray::FillBoundary() 6380 0.4081 0.4761 0.541 1.53% Castro::construct_new_source() 102 0.4682 0.4742 0.4847 1.37% MLMG::mgVcycle_up::3 128 0.01536 0.05025 0.4473 1.26% Castro::Castro() 1 0.4078 0.4111 0.4403 1.25% Castro::post_init() 1 0.4315 0.4348 0.4359 1.23% Castro::initMFs() 1 0.3893 0.3924 0.4186 1.18% MLCellLinOp::solutionResidual() 147 0.3994 0.411 0.4185 1.18% MLMG::mgVcycle_down::1 128 0.3735 0.3826 0.3917 1.11% Castro::reset_internal_energy(Fab) 160 0.368 0.3724 0.3775 1.07% MLMG::computeResidual() 128 0.3481 0.351 0.3539 1.00% FabArray::Saxpy() 111 0.3212 0.3267 0.3519 0.99% Castro::initData() 1 0.1905 0.2159 0.3512 0.99% Castro::construct_old_gravity_source() 40 0.3214 0.3258 0.333 0.94% FillBoundary_finish() 6380 0.214 0.2588 0.3315 0.94% Gravity::multilevel_solve_for_new_phi() 1 0.3303 0.3304 0.3304 0.93% Gravity::actual_multilevel_solve() 1 0.3303 0.3303 0.3304 0.93% Amr::bldFineLevels() 1 0.1089 0.2421 0.2701 0.76% Castro::construct_new_gravity_source() 17 0.2623 0.2655 0.269 0.76% FillBoundary_nowait() 6380 0.1498 0.1998 0.2558 0.72% Amr::writePlotFile() 2 0.2201 0.2202 0.2202 0.62% FabArray::Xpay() 915 0.1911 0.1998 0.218 0.62% Castro::post_regrid() 2 0.2043 0.2054 0.2087 0.59% MLMG::mgVcycle_down::2 128 0.01807 0.1094 0.2079 0.59% Castro::construct_old_gravity() 20 0.1558 0.1805 0.2077 0.59% Gravity::get_old_grav_vector() 20 0.1557 0.1804 0.2076 0.59% Castro::create_source_corrector() 17 0.172 0.1829 0.2039 0.58% MLMG::mgVcycle_down::3 128 0.0175 0.03902 0.1885 0.53% Castro::check_for_nan() 37 0.1645 0.1665 0.1688 0.48% MultiFab::contains_nan() 37 0.1643 0.1664 0.1686 0.48% Castro::derive() 18 0.1663 0.1669 0.1683 0.48% AmrLevel::derive() 18 0.1662 0.1668 0.1682 0.48% Amr::grid_places() 2 0.006171 0.1394 0.1678 0.47% AmrMesh::MakeNewGrids() 2 0.006114 0.1393 0.1677 0.47% TagBoxArray::mapPRD 2 0.001942 0.1113 0.1628 0.46% Castro::post_timestep() 1 0.1521 0.156 0.1592 0.45% Castro::normalize_species() 79 0.1384 0.14 0.1418 0.40% Castro::retry_advance_ctu() 11 0.1348 0.1367 0.1414 0.40% Gravity::get_new_grav_vector() 18 0.1269 0.1285 0.1301 0.37% Castro::initialize_advance() 1 0.1232 0.124 0.1281 0.36% TagBoxArray::collate() 2 0.0006933 0.02473 0.1257 0.36% amrex::average_down 768 0.0717 0.08689 0.1224 0.35% Castro::apply_source_to_state() 37 0.1074 0.1095 0.1182 0.33% Amr::regrid() 1 0.1054 0.1059 0.1076 0.30% MLMG::mgVcycle_bottom 128 0.0004641 0.007051 0.1057 0.30% MLMG::actualBottomSolve() 128 0.0001826 0.006767 0.1054 0.30% MLMG::mgVcycle_down::4 128 0.01746 0.02326 0.1039 0.29% MLCGSolver::bicgstab 8 0 0.006446 0.1031 0.29% StateDataPhysBCFunct::() 187 0.001031 0.05207 0.1009 0.29% StateData::FillBoundary(geom) 193 0 0.05085 0.09965 0.28% MLMG::mgVcycle_up::4 128 0.01546 0.01963 0.07651 0.22% Castro::save_data_for_retry() 3 0.07077 0.07154 0.07538 0.21% Castro::sum_integrated_quantities() 2 0.07316 0.07333 0.07335 0.21% MLMG::computeMLResidual() 18 0.04816 0.05752 0.06396 0.18% MLPoisson::define() 19 0.03513 0.04915 0.06394 0.18% amrex::Add() 128 0.05283 0.05334 0.05386 0.15% MLMG::prepareForSolve() 18 0.04471 0.04729 0.05318 0.15% VisMF::Write(FabArray) 4 0.03661 0.03791 0.05053 0.14% MLMG::mgVcycle_down::5 128 0.01742 0.01985 0.05012 0.14% MLMG::ResNormInf() 146 0.03739 0.04285 0.05007 0.14% FabArray::setDomainBndry() 187 0.0006585 0.02508 0.04844 0.14% Gravity::test_composite_phi() 1 0.04145 0.04149 0.04151 0.12% FabArray::norminf() 244 0.03613 0.03678 0.04096 0.12% MLMG::mgVcycle_up::5 128 0.01551 0.01742 0.04014 0.11% Castro::problem_post_timestep() 2 0.02529 0.02983 0.03342 0.09% Castro::wd_update() 2 0.02528 0.0298 0.03341 0.09% MLCellLinOp::defineAuxData() 19 0.02215 0.02415 0.02988 0.08% Castro::volWgtSum() 46 0.0283 0.0284 0.02859 0.08% FabArray::mult() 59 0.0262 0.0267 0.02756 0.08% Other 10888 0.1832 0.2376 0.3344 0.95% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Async Memory Usage: -------------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max -------------------------------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 6400 430 KiB 436 KiB 442 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 1943 B 1969 B 1988 B 89 KiB 89 KiB 89 KiB -------------------------------------------------------------------------------------------------------------------- Pinned Memory Usage: -------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max -------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 0 B 0 B 1 B 8192 KiB 8192 KiB 8192 KiB -------------------------------------------------------------------------------------------------------------- Cpu Memory Usage: ----------------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max ----------------------------------------------------------------------------------------------------------------------- Castro::save_data_for_retry() 416 1738 KiB 1739 KiB 1739 KiB 2329 KiB 2329 KiB 2329 KiB Castro::construct_ctu_hydro_source() 6400 430 KiB 436 KiB 442 KiB 1176 KiB 1176 KiB 1176 KiB StateData::define() 192 1119 KiB 1120 KiB 1120 KiB 1164 KiB 1164 KiB 1164 KiB Castro::initMFs() 256 1098 KiB 1098 KiB 1098 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 737 KiB 738 KiB 738 KiB 803 KiB 803 KiB 803 KiB Amr::writePlotFile() 64 4432 B 4459 B 4690 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 69 KiB 69 KiB 70 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 17 KiB 18 KiB 19 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 476 KiB 476 KiB 476 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 3464 B 4889 B 6281 B 225 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1238 B 1242 B 1253 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 29 KiB 29 KiB 30 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 143 KiB 143 KiB 143 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 136 KiB 136 KiB 136 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 12 KiB 12 KiB 15 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 1280 783 B 795 B 831 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 664 B 673 B 701 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 1943 B 1969 B 1988 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 19 B 20 B 29 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 35 KiB 35 KiB 35 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 61 B 61 B 61 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 4 B 5 B 6 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 164 B 166 B 168 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 202 B 234 B 270 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 34 KiB 34 KiB 34 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 309 B 313 B 318 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 1722 B 1740 B 1759 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 152 B 152 B 152 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 2128 B 2156 B 2170 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 1 B 1 B 1 B 16 KiB 16 KiB 16 KiB MLCellLinOp::defineAuxData() 11704 556 B 621 B 863 B 4608 B 5151 B 7152 B BndryData::define() 4864 584 B 584 B 585 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 21 B 28 B 36 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 0 B 8 B 12 B 2916 B 2916 B 2916 B OwnerMask() 64 0 B 0 B 0 B 2592 B 2592 B 2592 B MLMG::addInterpCorrection() 5632 1 B 17 B 28 B 1024 B 1024 B 1024 B amrex::average_down 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 0 B 1 B 3 B 648 B 648 B 648 B MLCGSolver::bicgstab 640 0 B 0 B 1 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 0 B 0 B 0 B 8 B 8 B 8 B ----------------------------------------------------------------------------------------------------------------------- AMReX (26.01-21-ge5851400bcf7) finalized