Initializing AMReX (26.08-22-g13221b4ab1c6)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.08-22-g13221b4ab1c6) initialized Starting run at 06:45:30 UTC on 2026-08-20. Successfully read inputs file ... Castro git describe: 26.08-1-gf0b335c8a AMReX git describe: 26.08-22-g13221b4ab Microphysics git describe: 26.08-13-g5ff9065d 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.189227603 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001164889 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010285964 Gravity::multilevel_solve_for_new_phi() time = 0.107267643 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.01041675 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.046831283 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004112935 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.002690383 on level 0 Castro::do_old_sources() time = 0.007429546 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.053152024 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.006042196 on level 0 Castro::do_old_sources() time = 0.00760823 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.052700197 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.002860728 on level 0 Castro::do_old_sources() time = 0.007573788 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.049900061 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.003382827 on level 0 Castro::do_old_sources() time = 0.007559091 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.049108885 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.293005707 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.013053622 Gravity::solve_for_phi() time = 0.088044023 on level 0 Castro::construct_new_gravity() time = 0.090082721 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.049250703 on level 0 Castro::construct_old_gravity() time = 0.002581546 on level 0 Castro::do_old_sources() time = 0.007346176 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.049313408 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.292939667 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.012999988 Gravity::solve_for_phi() time = 0.089846719 on level 0 Castro::construct_new_gravity() time = 0.091858965 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.051030352 on level 0 Castro::construct_old_gravity() time = 0.002821935 on level 0 Castro::do_old_sources() time = 0.007577553 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.049117052 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.300057623 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.013109397 Gravity::solve_for_phi() time = 0.089077195 on level 0 Castro::construct_new_gravity() time = 0.091112449 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.049155499 on level 0 Castro::construct_old_gravity() time = 0.002580006 on level 0 Castro::do_old_sources() time = 0.007385532 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.048978833 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.303948933 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.013552784 Gravity::solve_for_phi() time = 0.088743584 on level 0 Castro::construct_new_gravity() time = 0.090808494 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.049732282 on level 0 Castro::construct_old_gravity() time = 0.002598888 on level 0 Castro::do_old_sources() time = 0.007336234 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.049306444 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.30956484 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.013681355 Gravity::solve_for_phi() time = 0.08859746 on level 0 Castro::construct_new_gravity() time = 0.0906695 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.049585803 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.002759624 on level 0 Castro::do_old_sources() time = 0.007345083 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.049011383 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.301562192 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.013124107 Gravity::solve_for_phi() time = 0.088725195 on level 0 Castro::construct_new_gravity() time = 0.090749944 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.049218412 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.002619796 on level 0 Castro::do_old_sources() time = 0.00740011 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.048973708 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.307612271 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.013952742 Gravity::solve_for_phi() time = 0.090674065 on level 0 Castro::construct_new_gravity() time = 0.092743362 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.049311607 on level 0 Castro::construct_old_gravity() time = 0.002565865 on level 0 Castro::do_old_sources() time = 0.007421927 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.049227595 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.303710872 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.013798511 Gravity::solve_for_phi() time = 0.088484842 on level 0 Castro::construct_new_gravity() time = 0.090503582 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.04940096 on level 0 Castro::construct_old_gravity() time = 0.002678079 on level 0 Castro::do_old_sources() time = 0.00758379 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.049113069 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.313080567 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.013572812 Gravity::solve_for_phi() time = 0.09021695 on level 0 Castro::construct_new_gravity() time = 0.092305741 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.049544514 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.002720961 on level 0 Castro::do_old_sources() time = 0.007562503 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.049468645 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.308566817 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.013425359 Gravity::solve_for_phi() time = 0.088237427 on level 0 Castro::construct_new_gravity() time = 0.090340596 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.049284738 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.002654848 on level 0 Castro::do_old_sources() time = 0.00739453 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.049538204 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.319447314 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.014248411 Gravity::solve_for_phi() time = 0.089255355 on level 0 Castro::construct_new_gravity() time = 0.091271062 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.049238459 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.005268776 on level 0 Castro::do_old_sources() time = 0.008612589 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.051380921 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.308887335 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.013781827 Gravity::solve_for_phi() time = 0.090307373 on level 0 Castro::construct_new_gravity() time = 0.092393304 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.049584418 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.002764724 on level 0 Castro::do_old_sources() time = 0.007579831 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.049186175 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.313404217 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.013920009 Gravity::solve_for_phi() time = 0.090072259 on level 0 Castro::construct_new_gravity() time = 0.092092113 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.049588645 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.002759649 on level 0 Castro::do_old_sources() time = 0.007577564 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.049055156 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.309176136 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.0142699 Gravity::solve_for_phi() time = 0.090442399 on level 0 Castro::construct_new_gravity() time = 0.092508258 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.049566424 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.002708903 on level 0 Castro::do_old_sources() time = 0.007586202 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.049149124 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.315962992 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.014220467 Gravity::solve_for_phi() time = 0.090610998 on level 0 Castro::construct_new_gravity() time = 0.092633999 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.049778842 on level 0 Castro::construct_old_gravity() time = 0.002736591 on level 0 Castro::do_old_sources() time = 0.007579373 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.049070485 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.309439947 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.014242746 Gravity::solve_for_phi() time = 0.090437006 on level 0 Castro::construct_new_gravity() time = 0.092448073 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.049702693 on level 0 Castro::construct_old_gravity() time = 0.002737681 on level 0 Castro::do_old_sources() time = 0.007987044 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.049282238 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.337964812 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002923782298 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.013967921 Gravity::solve_for_phi() time = 0.09042191 on level 0 Castro::construct_new_gravity() time = 0.092441513 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.049730417 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.581215356e+29 ener added: 1.581215356e+29 TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -4.341809003e+38 TIME= 0.25 YMOM = 1.137357411e+28 TIME= 0.25 ZMOM = 0 TIME= 0.25 ANG MOM X = 0 TIME= 0.25 ANG MOM Y = 0 TIME= 0.25 ANG MOM Z = -1.469316165e+38 TIME= 0.25 RHO*e = 1.325894512e+50 TIME= 0.25 RHO*K = 3.170941987e+49 TIME= 0.25 RHO*E = 1.642988711e+50 TIME= 0.25 RHO*PHI = -5.449011378e+50 TIME= 0.25 TOTAL ENERGY= -1.081516979e+50 TIME= 0.25 CENTER OF MASS X-LOC = 316424388.8 TIME= 0.25 CENTER OF MASS X-VEL = -170267.2489 TIME= 0.25 CENTER OF MASS Y-LOC = 2.427336677e-06 TIME= 0.25 CENTER OF MASS Y-VEL = 4.460231145e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 57600050.72 TIME= 0.25 MAXIMUM DENSITY = 4458760.788 TIME= 0.25 MAXIMUM T_S / T_E = 0 STEP = 1 TIME = 0.25 DT = 0.25 [STEP 1] Coarse TimeStep time: 11.25762195 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.047526962 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004170365 seconds Ending run at 06:45:44 UTC on 2026-08-20. Run time = 13.66496628 Run time without initialization = 11.31018056 Average number of zones advanced per microsecond: 0.003 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 13.6536 ... 13.6585 ... 13.6651 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.1544 5.2027 5.2755 38.61% Amr::InitAmr() 1 1.435 1.44 1.447 10.59% Castro::computeTemp() 80 0.7107 0.7968 0.9237 6.76% amrex::Copy() 766 0.778 0.7846 0.7967 5.83% FabArray::ParallelCopy_finish() 1469 0.2372 0.546 0.7318 5.36% MultiFab::AddProduct() 3598 0.5487 0.5522 0.5604 4.10% MLPoisson::Fsmooth() 6144 0.3295 0.356 0.4316 3.16% Castro::initMFs() 1 0.3796 0.3917 0.416 3.04% FabArray::setVal() 5232 0.3575 0.3615 0.3686 2.70% Castro::reset_internal_energy(Fab) 160 0.3478 0.3507 0.3539 2.59% MultiFab::Subtract() 159 0.3312 0.3338 0.338 2.47% FabArray::ParallelCopy_nowait() 1469 0.2317 0.25 0.2732 2.00% Castro::do_advance_ctu() 20 0.258 0.2624 0.2686 1.97% Gravity::fill_multipole_BCs() 19 0.1901 0.2307 0.2501 1.83% check_for_negative_density() 20 0.06627 0.1483 0.2088 1.53% FabArray::sum() 3598 0.1914 0.1926 0.1953 1.43% Castro::initData() 1 0.003152 0.0287 0.1694 1.24% FillBoundary_finish() 6380 0.113 0.1347 0.164 1.20% FillBoundary_nowait() 6380 0.09107 0.1195 0.1545 1.13% TagBoxArray::collate() 2 0.0002487 0.02875 0.149 1.09% MLCellLinOp::applyBC() 7148 0.0709 0.08601 0.1336 0.98% Castro::normalize_species() 79 0.1195 0.1205 0.1226 0.90% MLPoisson::Fapply() 986 0.08508 0.08911 0.1044 0.76% Castro::construct_old_gravity_source() 40 0.09471 0.09601 0.1023 0.75% FabArray::Saxpy() 111 0.07561 0.07616 0.07693 0.56% Castro::construct_new_gravity_source() 17 0.06725 0.06838 0.07045 0.52% FabArray::Xpay() 915 0.04902 0.05121 0.05684 0.42% MultiFab::contains_nan() 37 0.04365 0.04408 0.04498 0.33% MLMG::addInterpCorrection() 768 0.03288 0.03604 0.04463 0.33% AmrLevel::derive() 18 0.03906 0.03959 0.04139 0.30% StateData::FillBoundary(geom) 193 0 0.01814 0.03523 0.26% amrex::average_down 768 0.02418 0.02649 0.03336 0.24% Castro::reset_internal_energy(MultiFab) 80 0.02505 0.02563 0.0262 0.19% VisMF::Write(FabArray) 4 0.01278 0.01374 0.0247 0.18% MLCellLinOp::defineAuxData() 19 0.01601 0.01704 0.02097 0.15% FabArray::FillBoundary() 6380 0.01653 0.01715 0.02045 0.15% Castro::construct_old_source() 240 0.01931 0.01956 0.0201 0.15% Castro::enforce_min_density() 79 0.01889 0.01933 0.0198 0.14% Amr::writePlotFile() 2 0.006862 0.01641 0.01796 0.13% MLCGSolver::bicgstab 8 0 0.001065 0.01704 0.12% FabArrayBase::getFB() 6380 0.0135 0.01399 0.01702 0.12% FabArray::setDomainBndry() 187 0.0006563 0.009056 0.01685 0.12% FillPatchIterator::Initialize 187 0.01566 0.01626 0.01671 0.12% Castro::evaluate_source_change() 196 0.01429 0.01472 0.01595 0.12% Castro::construct_new_source() 102 0.01523 0.01528 0.01539 0.11% FabArray::norminf() 244 0.0106 0.01108 0.01461 0.11% amrex::Add() 128 0.01341 0.01354 0.0137 0.10% MLCellLinOp::smooth() 1544 0.01073 0.01113 0.01298 0.09% MLMG::prepareForSolve() 18 0.01018 0.01059 0.01175 0.09% Castro::wd_update() 2 0.004936 0.008162 0.01095 0.08% Gravity::actual_solve_with_mlmg() 19 0.007415 0.007894 0.009595 0.07% MLMG::solve() 18 0.0007471 0.003652 0.009504 0.07% FabArrayBase::getCPC() 2045 0.008591 0.00886 0.009478 0.07% MLMG::ResNormInf() 146 0.002999 0.004956 0.006514 0.05% Gravity::get_old_grav_vector() 20 0.006104 0.006189 0.006277 0.05% Gravity::get_new_grav_vector() 18 0.005735 0.005831 0.006016 0.04% MLCellLinOp::apply() 986 0.002358 0.002592 0.005501 0.04% Castro::initialize_do_advance() 20 0.004764 0.00489 0.005108 0.04% FabArray::ParallelCopy() 1469 0.004102 0.004254 0.004694 0.03% AmrLevel::FillPatch() 187 0.003089 0.003221 0.003583 0.03% MLLinOp::defineGrids() 19 0.00257 0.00275 0.002987 0.02% Castro::sum_integrated_quantities() 2 0.001605 0.002553 0.002891 0.02% MLCellLinOp::correctionResidual() 768 0.002271 0.002309 0.002404 0.02% MLMG::mgVcycle() 128 0.002219 0.002277 0.002393 0.02% MLMG:computeResOfCorrection() 768 0.00183 0.001879 0.002014 0.01% Castro::save_data_for_retry() 3 0.001308 0.001349 0.001406 0.01% StateDataPhysBCFunct::() 187 0.0009853 0.001169 0.001294 0.01% MLMG::actualBottomSolve() 128 0.0001856 0.0002587 0.001278 0.01% Castro::subcycle_advance_ctu() 1 0.0002425 0.0003196 0.001038 0.01% main() 1 0.0007174 0.0008444 0.001009 0.01% FillPatchIterator::FillFromLevel0() 187 0.0008988 0.0009234 0.0009735 0.01% Gravity::solve_for_phi() 17 0.000701 0.0007359 0.0008161 0.01% FillPatchSingleLevel 187 0.0006531 0.0006759 0.0007438 0.01% MLMG::mgVcycle_down::2 128 0.0005963 0.000627 0.0007013 0.01% MLMG::mgVcycle_down::3 128 0.0005716 0.0005974 0.0006663 0.00% MLMG::mgVcycle_down::0 128 0.0006184 0.0006333 0.0006657 0.00% MLMG::mgVcycle_down::5 128 0.0005907 0.0006054 0.000664 0.00% MLMG::mgVcycle_down::1 128 0.0006082 0.0006206 0.0006617 0.00% MLMG::mgVcycle_down::4 128 0.0005615 0.000579 0.0006315 0.00% AmrMesh::MakeNewGrids() 2 0.0005451 0.0005792 0.0006309 0.00% Castro::Castro() 1 0.0004167 0.0004771 0.000619 0.00% MLCellLinOp::solutionResidual() 147 0.0005078 0.0005238 0.0005398 0.00% Castro::clean_state() 79 0.0004844 0.0004995 0.0005204 0.00% Castro::do_old_sources() 40 0.0004526 0.000466 0.0004985 0.00% MLMG::oneIter() 128 0.0004315 0.0004417 0.0004721 0.00% MLMG::mgVcycle_up::1 128 0.0004181 0.0004319 0.000467 0.00% MLMG::mgVcycle_up::0 128 0.0004153 0.0004249 0.000464 0.00% MLMG::mgVcycle_up::2 128 0.0003974 0.0004181 0.0004619 0.00% MLMG::mgVcycle_up::4 128 0.0003853 0.0003944 0.0004379 0.00% MLMG::mgVcycle_up::3 128 0.000383 0.0003963 0.0004357 0.00% MLMG::mgVcycle_up::5 128 0.0003728 0.0003805 0.0003965 0.00% MLMG::computeResidual() 128 0.0002889 0.0002991 0.000315 0.00% Castro::construct_new_gravity() 17 0.0001917 0.000213 0.0003147 0.00% MLMG::mgVcycle_bottom 128 0.000268 0.0002771 0.0003024 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002207 0.0002393 0.0002613 0.00% Castro::problem_post_timestep() 2 1.155e-05 2.828e-05 0.0002546 0.00% Castro::do_new_sources() 17 0.0002089 0.0002209 0.0002452 0.00% Amr::defBaseLevel() 1 0.0001163 0.0001329 0.0001555 0.00% Castro::construct_old_gravity() 20 0.0001264 0.0001379 0.000152 0.00% Amr::coarseTimeStep() 1 0.0001185 0.000146 0.0001515 0.00% MLPoisson::define() 19 0.0001265 0.0001365 0.0001446 0.00% Castro::check_for_nan() 37 0.0001208 0.0001283 0.000144 0.00% TagBoxArray::mapPRD 2 0.0001183 0.0001241 0.0001314 0.00% Castro::initialize_advance() 1 0.0001043 0.0001128 0.0001314 0.00% Castro::apply_source_to_state() 37 0.0001163 0.00012 0.0001305 0.00% Gravity::test_composite_phi() 1 0.0001094 0.000114 0.0001273 0.00% Castro::retry_advance_ctu() 11 7.001e-05 8.366e-05 0.0001263 0.00% Amr::FinalizeInit() 1 9.168e-06 1.857e-05 0.000105 0.00% MLLinOp::define() 19 8.38e-05 9.028e-05 9.668e-05 0.00% Amr::grid_places() 2 5.501e-05 6.383e-05 8.391e-05 0.00% Castro::create_source_corrector() 17 5.579e-05 6.131e-05 7.555e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.769e-05 6.011e-05 7.035e-05 0.00% Castro::expand_state() 20 5.697e-05 6.16e-05 6.75e-05 0.00% MLMG::computeMLResidual() 18 5.205e-05 5.438e-05 5.804e-05 0.00% Castro::derive() 18 4.734e-05 4.973e-05 5.472e-05 0.00% Amr::regrid() 1 1.501e-05 1.927e-05 4.266e-05 0.00% Castro::advance() 1 7.863e-06 1.076e-05 3.953e-05 0.00% Gravity::actual_multilevel_solve() 1 3.233e-05 3.317e-05 3.715e-05 0.00% Castro::post_init() 1 1.59e-05 1.821e-05 3.285e-05 0.00% Amr::timeStep() 1 1.492e-05 1.802e-05 2.478e-05 0.00% Amr::InitializeInit() 1 8.112e-06 9.451e-06 1.582e-05 0.00% Castro::post_regrid() 2 1.095e-05 1.201e-05 1.384e-05 0.00% Amr::bldFineLevels() 1 8.107e-06 9.137e-06 1.369e-05 0.00% Amr::init() 1 6.758e-06 7.746e-06 8.723e-06 0.00% Castro::post_timestep() 1 6.074e-06 6.411e-06 6.776e-06 0.00% Amr::initialInit() 1 3.164e-06 3.478e-06 4.178e-06 0.00% Other 2530 0.04744 0.05858 0.09586 0.70% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 13.65 13.66 13.67 100.00% Amr::coarseTimeStep() 1 11.31 11.31 11.31 82.76% Amr::timeStep() 1 11.26 11.26 11.26 82.38% Castro::advance() 1 11.2 11.21 11.21 82.02% Castro::subcycle_advance_ctu() 1 11.16 11.16 11.17 81.71% Castro::do_advance_ctu() 20 11.12 11.13 11.13 81.45% Castro::construct_ctu_hydro_source() 20 5.957 5.982 6.019 44.05% Castro::clean_state() 79 2.924 3.004 3.132 22.92% Castro::computeTemp() 80 2.035 2.117 2.246 16.43% Gravity::solve_phi_with_mlmg() 19 1.574 1.619 1.638 11.99% Castro::construct_new_gravity() 17 1.492 1.536 1.555 11.38% Gravity::solve_for_phi() 17 1.456 1.501 1.52 11.12% Amr::InitAmr() 1 1.435 1.44 1.447 10.59% Gravity::actual_solve_with_mlmg() 19 1.381 1.385 1.389 10.16% Castro::reset_internal_energy(MultiFab) 80 1.309 1.32 1.339 9.80% MLMG::solve() 18 1.324 1.33 1.332 9.75% MLMG::oneIter() 128 1.199 1.2 1.202 8.80% MLMG::mgVcycle() 128 1.185 1.186 1.188 8.69% Castro::initialize_do_advance() 20 1.117 1.146 1.17 8.56% Castro::do_old_sources() 40 1.089 1.115 1.137 8.32% FabArray::ParallelCopy() 1469 0.5125 0.8208 1.004 7.35% Castro::evaluate_source_change() 196 0.9509 0.9572 0.9717 7.11% Amr::init() 1 0.9077 0.9078 0.9078 6.64% Amr::initialInit() 1 0.8562 0.8563 0.8564 6.27% Castro::do_new_sources() 17 0.7672 0.8038 0.8409 6.15% MLCellLinOp::smooth() 1544 0.6115 0.6679 0.8029 5.88% Castro::enforce_min_density() 79 0.785 0.7901 0.8008 5.86% amrex::Copy() 766 0.778 0.7846 0.7967 5.83% FabArray::ParallelCopy_finish() 1469 0.2372 0.546 0.7318 5.36% AmrLevel::FillPatch() 187 0.5767 0.635 0.6916 5.06% Amr::InitializeInit() 1 0.4539 0.4921 0.6418 4.70% Amr::defBaseLevel() 1 0.4539 0.4921 0.6417 4.70% MultiFab::AddProduct() 3598 0.5487 0.5522 0.5604 4.10% FillPatchIterator::Initialize 187 0.383 0.4416 0.4991 3.65% FillPatchIterator::FillFromLevel0() 187 0.3501 0.4163 0.477 3.49% FillPatchSingleLevel 187 0.3492 0.4153 0.476 3.48% MLMG::mgVcycle_up::1 128 0.1054 0.2834 0.4604 3.37% MLCellLinOp::applyBC() 7148 0.3335 0.372 0.4536 3.32% MLPoisson::Fsmooth() 6144 0.3295 0.356 0.4316 3.16% Castro::Castro() 1 0.3848 0.3973 0.4233 3.10% Castro::initMFs() 1 0.3798 0.3919 0.4162 3.05% MLMG::addInterpCorrection() 768 0.06634 0.3252 0.4036 2.95% Amr::FinalizeInit() 1 0.2145 0.3642 0.4024 2.94% FabArray::setVal() 5232 0.3575 0.3615 0.3686 2.70% Castro::reset_internal_energy(Fab) 160 0.3478 0.3507 0.3539 2.59% MultiFab::Subtract() 159 0.3312 0.3338 0.338 2.47% FabArray::FillBoundary() 6380 0.2517 0.286 0.3212 2.35% MLMG::mgVcycle_up::2 128 0.01563 0.1282 0.3028 2.22% FabArray::ParallelCopy_nowait() 1469 0.2516 0.2705 0.2953 2.16% MLMG::mgVcycle_down::0 128 0.2611 0.2625 0.2649 1.94% Gravity::fill_multipole_BCs() 19 0.1921 0.2326 0.252 1.84% Castro::expand_state() 20 0.1897 0.2175 0.2434 1.78% Castro::initData() 1 0.06607 0.09308 0.2328 1.70% Amr::bldFineLevels() 1 0.03737 0.1876 0.2263 1.66% check_for_negative_density() 20 0.06627 0.1483 0.2088 1.53% MLMG::mgVcycle_up::0 128 0.2003 0.2012 0.2027 1.48% MLCellLinOp::apply() 986 0.1481 0.1609 0.2023 1.48% MLMG::mgVcycle_up::3 128 0.01553 0.02964 0.1992 1.46% FabArray::sum() 3598 0.1914 0.1926 0.1953 1.43% Amr::grid_places() 2 0.002548 0.1522 0.1904 1.39% AmrMesh::MakeNewGrids() 2 0.002488 0.1521 0.1903 1.39% TagBoxArray::mapPRD 2 0.000772 0.122 0.1886 1.38% FillBoundary_nowait() 6380 0.1049 0.1341 0.1696 1.24% Castro::construct_old_source() 240 0.1559 0.1577 0.1641 1.20% FillBoundary_finish() 6380 0.113 0.1347 0.164 1.20% TagBoxArray::collate() 2 0.0002487 0.02875 0.149 1.09% MLMG:computeResOfCorrection() 768 0.1154 0.1249 0.1488 1.09% MLCellLinOp::correctionResidual() 768 0.1135 0.1231 0.1469 1.07% Castro::post_init() 1 0.1406 0.1416 0.1426 1.04% MLMG::mgVcycle_down::1 128 0.125 0.1287 0.1325 0.97% Castro::normalize_species() 79 0.1195 0.1205 0.1226 0.90% Castro::construct_new_source() 102 0.1182 0.1194 0.1218 0.89% Gravity::multilevel_solve_for_new_phi() 1 0.1072 0.1072 0.1073 0.79% Gravity::actual_multilevel_solve() 1 0.1071 0.1071 0.1072 0.78% MLPoisson::Fapply() 986 0.08508 0.08911 0.1044 0.76% Castro::construct_old_gravity_source() 40 0.09471 0.09601 0.1023 0.75% MLCellLinOp::solutionResidual() 147 0.08339 0.08987 0.09515 0.70% Amr::writePlotFile() 2 0.09479 0.09488 0.0949 0.69% MLMG::mgVcycle_bottom 128 0.0004563 0.005605 0.08269 0.61% MLMG::actualBottomSolve() 128 0.0001856 0.005328 0.08239 0.60% MLMG::mgVcycle_down::2 128 0.01814 0.04864 0.08237 0.60% MLCGSolver::bicgstab 8 0 0.005022 0.08035 0.59% FabArray::Saxpy() 111 0.07561 0.07616 0.07693 0.56% MLMG::computeResidual() 128 0.07226 0.07268 0.07298 0.53% Castro::post_regrid() 2 0.06938 0.07038 0.0713 0.52% Castro::derive() 18 0.06677 0.06765 0.0705 0.52% Castro::construct_new_gravity_source() 17 0.06725 0.06838 0.07045 0.52% AmrLevel::derive() 18 0.06673 0.06761 0.07045 0.52% MLMG::mgVcycle_down::3 128 0.01767 0.02411 0.06787 0.50% Castro::create_source_corrector() 17 0.05027 0.0592 0.06591 0.48% amrex::average_down 768 0.03389 0.04116 0.05993 0.44% Castro::construct_old_gravity() 20 0.04465 0.05116 0.05732 0.42% Gravity::get_old_grav_vector() 20 0.04451 0.05102 0.05717 0.42% FabArray::Xpay() 915 0.04902 0.05121 0.05684 0.42% Castro::post_timestep() 1 0.05018 0.05245 0.05428 0.40% Castro::check_for_nan() 37 0.04378 0.04421 0.04512 0.33% MultiFab::contains_nan() 37 0.04365 0.04408 0.04498 0.33% Castro::initialize_advance() 1 0.04021 0.04082 0.04405 0.32% MLPoisson::define() 19 0.0323 0.03753 0.04344 0.32% MLMG::mgVcycle_down::4 128 0.01759 0.01942 0.04135 0.30% Castro::retry_advance_ctu() 11 0.03385 0.03478 0.03811 0.28% Amr::regrid() 1 0.03597 0.03657 0.03711 0.27% StateDataPhysBCFunct::() 187 0.0009853 0.01931 0.03652 0.27% Gravity::get_new_grav_vector() 18 0.0343 0.03479 0.03525 0.26% StateData::FillBoundary(geom) 193 0 0.01814 0.03523 0.26% MLMG::mgVcycle_up::4 128 0.0156 0.01693 0.0328 0.24% MLMG::mgVcycle_down::5 128 0.01764 0.01871 0.03031 0.22% MLCellLinOp::defineAuxData() 19 0.02054 0.02242 0.02839 0.21% Castro::apply_source_to_state() 37 0.02531 0.02553 0.02598 0.19% MLMG::mgVcycle_up::5 128 0.01565 0.01656 0.02597 0.19% VisMF::Write(FabArray) 4 0.01278 0.01374 0.0247 0.18% MLMG::prepareForSolve() 18 0.02029 0.02127 0.02421 0.18% Castro::sum_integrated_quantities() 2 0.02197 0.02212 0.02213 0.16% MLMG::computeMLResidual() 18 0.01042 0.01675 0.02185 0.16% Castro::save_data_for_retry() 3 0.01762 0.01814 0.02121 0.16% FabArrayBase::getFB() 6380 0.01381 0.01468 0.01848 0.14% FabArray::setDomainBndry() 187 0.0006563 0.009056 0.01685 0.12% FabArrayBase::getCPC() 2045 0.01106 0.01238 0.0161 0.12% MLMG::ResNormInf() 146 0.01243 0.01439 0.01582 0.12% FabArray::norminf() 244 0.0106 0.01108 0.01461 0.11% Castro::problem_post_timestep() 2 0.009186 0.01186 0.01425 0.10% Gravity::test_composite_phi() 1 0.01416 0.01422 0.01424 0.10% Castro::wd_update() 2 0.008932 0.01183 0.01424 0.10% amrex::Add() 128 0.01341 0.01354 0.0137 0.10% MLLinOp::define() 19 0.004321 0.007831 0.01122 0.08% MLLinOp::defineGrids() 19 0.004232 0.00774 0.01113 0.08% Other 2530 0.07852 0.09037 0.1342 0.98% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Async Memory Usage: -------------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max -------------------------------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 6400 399 KiB 403 KiB 408 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 2134 B 2171 B 2196 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 2 B 2 B 3 B 8192 KiB 8192 KiB 8192 KiB -------------------------------------------------------------------------------------------------------------- Cpu Memory Usage: ----------------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max ----------------------------------------------------------------------------------------------------------------------- Castro::save_data_for_retry() 416 1565 KiB 1565 KiB 1566 KiB 2332 KiB 2332 KiB 2332 KiB Castro::construct_ctu_hydro_source() 6400 400 KiB 404 KiB 409 KiB 1178 KiB 1178 KiB 1178 KiB StateData::define() 192 1041 KiB 1041 KiB 1041 KiB 1166 KiB 1166 KiB 1166 KiB Castro::initMFs() 256 1020 KiB 1020 KiB 1020 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 665 KiB 665 KiB 665 KiB 804 KiB 804 KiB 804 KiB Amr::writePlotFile() 64 4407 B 4495 B 5081 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 51 KiB 51 KiB 52 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 12 KiB 14 KiB 15 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 428 KiB 428 KiB 428 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2345 B 4206 B 5457 B 226 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1418 B 1436 B 1491 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 19 KiB 19 KiB 19 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 129 KiB 129 KiB 129 KiB 158 KiB 158 KiB 158 KiB Castro::buildMetrics() 192 127 KiB 127 KiB 127 KiB 143 KiB 143 KiB 143 KiB MLMG::prepareForSolve() 6102 10 KiB 10 KiB 12 KiB 104 KiB 110 KiB 130 KiB Castro::construct_old_source() 1280 497 B 501 B 508 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 419 B 420 B 422 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 2139 B 2177 B 2202 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 19 B 21 B 39 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 33 KiB 33 KiB 33 KiB 76 KiB 76 KiB 76 KiB Gravity::test_composite_phi() 96 55 B 55 B 55 B 53 KiB 53 KiB 53 KiB Gravity::update_max_rhs() 192 3 B 4 B 6 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 114 B 115 B 116 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 149 B 171 B 192 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 31 KiB 31 KiB 31 KiB 38 KiB 38 KiB 38 KiB Gravity::fill_multipole_BCs() 2432 247 B 256 B 264 B 18 KiB 18 KiB 18 KiB MLCellLinOp::defineAuxData() 11704 818 B 983 B 1835 B 8192 B 9856 B 18 KiB Castro::evaluate_source_change() 6272 1134 B 1141 B 1158 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 128 B 128 B 128 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 1744 B 1797 B 1820 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 0 B 0 B 1 B 16 KiB 16 KiB 16 KiB BndryData::define() 4864 806 B 809 B 810 B 8192 B 8192 B 8192 B FillBoundary_nowait() 91296 30 B 41 B 52 B 3072 B 3840 B 4608 B AmrMesh::MakeNewGrids() 128 0 B 27 B 42 B 3584 B 3584 B 3584 B OwnerMask() 64 0 B 0 B 0 B 3072 B 3072 B 3072 B MLCGSolver::bicgstab 640 0 B 0 B 6 B 0 B 80 B 1280 B MLMG::addInterpCorrection() 5632 1 B 17 B 27 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 0 B 10 B 14 B 1024 B 1024 B 1024 B amrex::average_down 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B MLCellLinOp::applyBC() 114375 0 B 0 B 1 B 256 B 256 B 256 B ----------------------------------------------------------------------------------------------------------------------- AMReX (26.08-22-g13221b4ab1c6) finalized