Initializing AMReX (26.03-91-gf546b598fb75)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.03-91-gf546b598fb75) initialized Starting run at 07:43:31 UTC on 2026-03-18. Successfully read inputs file ... Castro git describe: 26.03-22-g44e4ce9a5 AMReX git describe: 26.03-91-gf546b598f Microphysics git describe: 26.03-7-g736a456b 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.187779068 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001179276 new finest: 0 Gravity::fill_multipole_BCs() time = 0.01054103 Gravity::multilevel_solve_for_new_phi() time = 0.105807691 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010245313 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.046615653 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004097904 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.002552992 on level 0 Castro::do_old_sources() time = 0.007213053 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.052844834 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.007436805 on level 0 Castro::do_old_sources() time = 0.007385886 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.052374244 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.043569208 on level 0 Castro::do_old_sources() time = 0.024937206 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.072206402 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.005227851 on level 0 Castro::do_old_sources() time = 0.008120982 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.049072542 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.298320107 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.014384469 Gravity::solve_for_phi() time = 0.091899978 on level 0 Castro::construct_new_gravity() time = 0.093916971 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.049899243 on level 0 Castro::construct_old_gravity() time = 0.002631222 on level 0 Castro::do_old_sources() time = 0.00717383 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.048397557 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.302226901 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.012889234 Gravity::solve_for_phi() time = 0.090854733 on level 0 Castro::construct_new_gravity() time = 0.092855757 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.048911187 on level 0 Castro::construct_old_gravity() time = 0.004728103 on level 0 Castro::do_old_sources() time = 0.008710755 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.048592238 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.304279568 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.012934084 Gravity::solve_for_phi() time = 0.090942006 on level 0 Castro::construct_new_gravity() time = 0.092904402 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.048906729 on level 0 Castro::construct_old_gravity() time = 0.003223324 on level 0 Castro::do_old_sources() time = 0.00722813 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.048793229 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.301859146 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.013231712 Gravity::solve_for_phi() time = 0.087182981 on level 0 Castro::construct_new_gravity() time = 0.089155672 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.049121461 on level 0 Castro::construct_old_gravity() time = 0.002746014 on level 0 Castro::do_old_sources() time = 0.007367273 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.048303802 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.309236049 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.013488065 Gravity::solve_for_phi() time = 0.087757704 on level 0 Castro::construct_new_gravity() time = 0.089786729 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.048796077 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.002565928 on level 0 Castro::do_old_sources() time = 0.007238655 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.048269939 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.301241798 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.013498302 Gravity::solve_for_phi() time = 0.088060715 on level 0 Castro::construct_new_gravity() time = 0.090087765 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.048706258 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.002608166 on level 0 Castro::do_old_sources() time = 0.007302782 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.048530709 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.305663564 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.013282013 Gravity::solve_for_phi() time = 0.087964108 on level 0 Castro::construct_new_gravity() time = 0.089949753 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.048713065 on level 0 Castro::construct_old_gravity() time = 0.003001425 on level 0 Castro::do_old_sources() time = 0.007553571 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.048516875 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.301116432 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.01343222 Gravity::solve_for_phi() time = 0.08817792 on level 0 Castro::construct_new_gravity() time = 0.090158796 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.049467525 on level 0 Castro::construct_old_gravity() time = 0.007107489 on level 0 Castro::do_old_sources() time = 0.010898973 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.084642627 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.365671519 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.013454625 Gravity::solve_for_phi() time = 0.087837944 on level 0 Castro::construct_new_gravity() time = 0.089803452 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.049128397 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.003020098 on level 0 Castro::do_old_sources() time = 0.00723819 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.048683679 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.30265211 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.013381484 Gravity::solve_for_phi() time = 0.088513192 on level 0 Castro::construct_new_gravity() time = 0.090522388 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.049256977 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.002675343 on level 0 Castro::do_old_sources() time = 0.007244238 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.054074351 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.309219394 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.013619916 Gravity::solve_for_phi() time = 0.088058045 on level 0 Castro::construct_new_gravity() time = 0.090084614 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.04878536 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.003150601 on level 0 Castro::do_old_sources() time = 0.008307736 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.048268854 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.30248315 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.014016086 Gravity::solve_for_phi() time = 0.088208636 on level 0 Castro::construct_new_gravity() time = 0.090232322 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.052457645 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.002710992 on level 0 Castro::do_old_sources() time = 0.007187124 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.048240891 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.306571075 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.013651942 Gravity::solve_for_phi() time = 0.087572391 on level 0 Castro::construct_new_gravity() time = 0.089539674 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.049530959 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.002709927 on level 0 Castro::do_old_sources() time = 0.00726621 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.049838802 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.305281557 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.013905758 Gravity::solve_for_phi() time = 0.08852193 on level 0 Castro::construct_new_gravity() time = 0.090480438 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.04902286 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.002974809 on level 0 Castro::do_old_sources() time = 0.007436538 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.048346647 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.310521423 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.013784614 Gravity::solve_for_phi() time = 0.087904472 on level 0 Castro::construct_new_gravity() time = 0.089930612 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.048724313 on level 0 Castro::construct_old_gravity() time = 0.006566092 on level 0 Castro::do_old_sources() time = 0.008035462 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.050683987 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.308606578 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.014962915 Gravity::solve_for_phi() time = 0.090269889 on level 0 Castro::construct_new_gravity() time = 0.092300735 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.050023484 on level 0 Castro::construct_old_gravity() time = 0.002853333 on level 0 Castro::do_old_sources() time = 0.007510863 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.049655833 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.326165737 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002889518778 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.027390635 Gravity::solve_for_phi() time = 0.15179836 on level 0 Castro::construct_new_gravity() time = 0.153837771 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.051944741 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.39561072 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.046711156 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004130632 seconds Ending run at 07:43:45 UTC on 2026-03-18. Run time = 13.74930733 Run time without initialization = 11.44737324 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.7411 ... 13.7426 ... 13.7494 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.1301 5.1703 5.2212 37.97% Amr::InitAmr() 1 1.383 1.385 1.395 10.14% Castro::computeTemp() 80 0.704 0.7951 0.9238 6.72% FabArray::ParallelCopy_finish() 1469 0.2995 0.6778 0.8858 6.44% amrex::Copy() 766 0.7654 0.7722 0.7966 5.79% MultiFab::AddProduct() 3598 0.5502 0.5546 0.5677 4.13% MLPoisson::Fsmooth() 6144 0.3287 0.3556 0.4438 3.23% Castro::initMFs() 1 0.3865 0.3913 0.4087 2.97% FabArray::setVal() 5232 0.3622 0.3667 0.3789 2.76% Castro::reset_internal_energy(Fab) 160 0.3433 0.3464 0.3515 2.56% MultiFab::Subtract() 159 0.3288 0.3302 0.3333 2.42% Castro::do_advance_ctu() 20 0.2594 0.264 0.2849 2.07% FabArray::ParallelCopy_nowait() 1469 0.2253 0.2418 0.2627 1.91% Gravity::fill_multipole_BCs() 19 0.2046 0.244 0.2619 1.90% check_for_negative_density() 20 0.0929 0.165 0.258 1.88% FabArray::sum() 3598 0.1855 0.1871 0.1917 1.39% FillBoundary_finish() 6380 0.1357 0.1595 0.178 1.29% TagBoxArray::collate() 2 0.0002911 0.03637 0.1757 1.28% Castro::initData() 1 0.003214 0.02816 0.1658 1.21% FillBoundary_nowait() 6380 0.0908 0.1194 0.1498 1.09% MLCellLinOp::applyBC() 7148 0.06939 0.08501 0.1337 0.97% Castro::normalize_species() 79 0.1195 0.1205 0.1224 0.89% MLPoisson::Fapply() 986 0.08541 0.08992 0.1093 0.79% Castro::construct_old_gravity_source() 40 0.09252 0.09333 0.09469 0.69% FabArray::Saxpy() 111 0.07664 0.07712 0.07797 0.57% Castro::construct_new_gravity_source() 17 0.06794 0.06866 0.06962 0.51% FabArray::Xpay() 915 0.04889 0.05106 0.05859 0.43% MultiFab::contains_nan() 37 0.0439 0.04428 0.04492 0.33% MLMG::addInterpCorrection() 768 0.03259 0.03505 0.04362 0.32% AmrLevel::derive() 18 0.03867 0.03918 0.0397 0.29% StateData::FillBoundary(geom) 193 0 0.01861 0.03655 0.27% amrex::average_down 768 0.02422 0.0264 0.03389 0.25% Castro::reset_internal_energy(MultiFab) 80 0.02432 0.02463 0.02522 0.18% VisMF::Write(FabArray) 4 0.01355 0.01448 0.02434 0.18% FabArray::FillBoundary() 6380 0.0168 0.01736 0.02066 0.15% Castro::construct_old_source() 240 0.01916 0.01947 0.02025 0.15% MLCellLinOp::defineAuxData() 19 0.01584 0.01667 0.0199 0.14% Castro::enforce_min_density() 79 0.0184 0.01862 0.01885 0.14% Castro::wd_update() 2 0.005223 0.01444 0.0176 0.13% MLCGSolver::bicgstab 8 0 0.001066 0.01706 0.12% FabArrayBase::getFB() 6380 0.0133 0.01381 0.01687 0.12% Amr::writePlotFile() 2 0.00678 0.01531 0.01667 0.12% FabArray::setDomainBndry() 187 0.0006764 0.008684 0.01664 0.12% Castro::evaluate_source_change() 196 0.01438 0.01466 0.01606 0.12% FillPatchIterator::Initialize 187 0.0152 0.0155 0.01592 0.12% Castro::construct_new_source() 102 0.01497 0.01515 0.01591 0.12% FabArray::norminf() 244 0.01056 0.01095 0.01481 0.11% amrex::Add() 128 0.0136 0.01367 0.0142 0.10% MLCellLinOp::smooth() 1544 0.01084 0.01134 0.01366 0.10% MLMG::ResNormInf() 146 0.006541 0.009559 0.01307 0.10% MLMG::prepareForSolve() 18 0.01012 0.01046 0.01145 0.08% Gravity::actual_solve_with_mlmg() 19 0.007434 0.007888 0.009526 0.07% FabArrayBase::getCPC() 2045 0.008657 0.008854 0.009451 0.07% MLMG::solve() 18 0.0008436 0.003143 0.007542 0.05% Gravity::get_old_grav_vector() 20 0.006065 0.006141 0.006271 0.05% MLCellLinOp::apply() 986 0.002428 0.002692 0.005848 0.04% Gravity::get_new_grav_vector() 18 0.0056 0.005711 0.00583 0.04% MLCGSolver::ParallelAllReduce 189 0 0.0003579 0.005726 0.04% Castro::initialize_do_advance() 20 0.004698 0.004859 0.005051 0.04% FabArray::ParallelCopy() 1469 0.004246 0.004405 0.004872 0.04% AmrLevel::FillPatch() 187 0.003071 0.003239 0.003624 0.03% Castro::sum_integrated_quantities() 2 0.001529 0.002466 0.002623 0.02% MLCellLinOp::correctionResidual() 768 0.002288 0.002325 0.002433 0.02% MLMG::mgVcycle() 128 0.002179 0.002247 0.002405 0.02% MLLinOp::defineGrids() 19 0.001926 0.002043 0.002331 0.02% MLMG:computeResOfCorrection() 768 0.001802 0.001836 0.002018 0.01% Castro::save_data_for_retry() 3 0.001291 0.001339 0.001499 0.01% StateDataPhysBCFunct::() 187 0.0009686 0.001175 0.001393 0.01% MLMG::actualBottomSolve() 128 0.0001898 0.0002624 0.001293 0.01% Castro::subcycle_advance_ctu() 1 0.0002541 0.0003305 0.001054 0.01% FillPatchIterator::FillFromLevel0() 187 0.0008693 0.0009086 0.0009732 0.01% main() 1 0.0008144 0.0008519 0.0009595 0.01% Gravity::solve_for_phi() 17 0.0006608 0.0006885 0.0007861 0.01% MLMG::mgVcycle_down::2 128 0.0005932 0.0006256 0.0007111 0.01% MLMG::mgVcycle_down::1 128 0.000611 0.0006311 0.0007001 0.01% MLMG::mgVcycle_down::0 128 0.0006333 0.0006477 0.0006998 0.01% FillPatchSingleLevel 187 0.0006439 0.0006637 0.0006932 0.01% MLMG::mgVcycle_down::5 128 0.0005881 0.0006068 0.0006862 0.00% MLMG::mgVcycle_down::3 128 0.0005642 0.0005927 0.0006759 0.00% MLMG::mgVcycle_down::4 128 0.0005599 0.0005771 0.0006404 0.00% Castro::Castro() 1 0.0004424 0.000512 0.0006019 0.00% AmrMesh::MakeNewGrids() 2 0.0005196 0.0005358 0.0005599 0.00% MLCellLinOp::solutionResidual() 147 0.00051 0.0005252 0.0005525 0.00% Castro::clean_state() 79 0.0004782 0.0005034 0.0005519 0.00% Castro::do_old_sources() 40 0.0004565 0.0004769 0.0005156 0.00% MLMG::oneIter() 128 0.0004324 0.0004449 0.0005079 0.00% MLMG::mgVcycle_up::0 128 0.0004115 0.0004272 0.0004727 0.00% MLMG::mgVcycle_up::1 128 0.0004071 0.00042 0.0004654 0.00% MLMG::mgVcycle_up::2 128 0.0003913 0.0004091 0.0004611 0.00% MLMG::mgVcycle_up::3 128 0.0003773 0.00039 0.0004498 0.00% MLMG::mgVcycle_up::4 128 0.0003795 0.0003911 0.0004398 0.00% MLMG::mgVcycle_up::5 128 0.0003722 0.0003799 0.0004138 0.00% MLMG::computeResidual() 128 0.0002915 0.0003016 0.0003406 0.00% MLMG::mgVcycle_bottom 128 0.0002687 0.0002795 0.0003056 0.00% Castro::construct_new_gravity() 17 0.0001703 0.0001866 0.0002929 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002159 0.00023 0.0002712 0.00% Castro::do_new_sources() 17 0.0002147 0.0002245 0.0002534 0.00% Castro::problem_post_timestep() 2 1.173e-05 2.837e-05 0.0002489 0.00% Amr::defBaseLevel() 1 0.0001123 0.0001323 0.0001632 0.00% Amr::coarseTimeStep() 1 0.0001241 0.0001516 0.0001581 0.00% MLPoisson::define() 19 0.0001243 0.0001335 0.0001454 0.00% Castro::construct_old_gravity() 20 0.0001189 0.0001281 0.0001425 0.00% Castro::check_for_nan() 37 0.0001212 0.0001268 0.0001421 0.00% Castro::retry_advance_ctu() 11 7.253e-05 8.939e-05 0.0001407 0.00% Castro::initialize_advance() 1 0.0001085 0.0001218 0.0001391 0.00% Castro::apply_source_to_state() 37 0.0001135 0.0001188 0.000134 0.00% Gravity::test_composite_phi() 1 0.0001052 0.0001093 0.0001245 0.00% TagBoxArray::mapPRD 2 0.0001151 0.0001188 0.0001226 0.00% Amr::grid_places() 2 5.363e-05 5.968e-05 9.943e-05 0.00% MLLinOp::define() 19 8.487e-05 9.081e-05 9.909e-05 0.00% Amr::FinalizeInit() 1 9.249e-06 2.02e-05 7.259e-05 0.00% MLMG::computeMLResidual() 18 5.64e-05 6.228e-05 6.825e-05 0.00% Castro::create_source_corrector() 17 5.337e-05 5.828e-05 6.596e-05 0.00% Castro::expand_state() 20 5.875e-05 6.194e-05 6.583e-05 0.00% Amr::regrid() 1 2.009e-05 3.038e-05 6.218e-05 0.00% Castro::derive() 18 4.789e-05 5.029e-05 6.051e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.739e-05 5.859e-05 6.043e-05 0.00% MLMG::getGradSolution() 18 5.066e-05 5.369e-05 5.645e-05 0.00% Gravity::actual_multilevel_solve() 1 3.076e-05 3.14e-05 3.226e-05 0.00% Amr::timeStep() 1 1.473e-05 1.878e-05 2.573e-05 0.00% Amr::InitializeInit() 1 9.644e-06 1.213e-05 1.987e-05 0.00% Castro::post_regrid() 2 1.147e-05 1.329e-05 1.854e-05 0.00% Castro::post_init() 1 1.235e-05 1.323e-05 1.427e-05 0.00% Castro::advance() 1 8.065e-06 9.279e-06 1.303e-05 0.00% Amr::bldFineLevels() 1 7.456e-06 9.169e-06 1.068e-05 0.00% Amr::init() 1 6.785e-06 7.744e-06 9.213e-06 0.00% Castro::post_timestep() 1 6.318e-06 7.126e-06 8.949e-06 0.00% Amr::initialInit() 1 3.344e-06 3.637e-06 4.29e-06 0.00% Other 2323 0.04683 0.0606 0.1008 0.73% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 13.74 13.74 13.75 100.00% Amr::coarseTimeStep() 1 11.45 11.45 11.45 83.25% Amr::timeStep() 1 11.4 11.4 11.4 82.88% Castro::advance() 1 11.33 11.34 11.34 82.48% Castro::subcycle_advance_ctu() 1 11.29 11.29 11.3 82.18% Castro::do_advance_ctu() 20 11.25 11.26 11.26 81.91% Castro::construct_ctu_hydro_source() 20 5.936 5.962 6.023 43.80% Castro::clean_state() 79 2.883 2.983 3.125 22.73% Castro::computeTemp() 80 2.007 2.102 2.239 16.29% Gravity::solve_phi_with_mlmg() 19 1.622 1.665 1.686 12.27% Castro::construct_new_gravity() 17 1.541 1.584 1.605 11.67% Gravity::solve_for_phi() 17 1.506 1.549 1.57 11.42% Gravity::actual_solve_with_mlmg() 19 1.414 1.418 1.422 10.34% Amr::InitAmr() 1 1.383 1.385 1.395 10.14% MLMG::solve() 18 1.359 1.365 1.366 9.94% Castro::reset_internal_energy(MultiFab) 80 1.296 1.307 1.329 9.67% MLMG::oneIter() 128 1.227 1.23 1.233 8.97% MLMG::mgVcycle() 128 1.213 1.216 1.219 8.87% Castro::do_old_sources() 40 1.086 1.148 1.182 8.60% Castro::initialize_do_advance() 20 1.131 1.15 1.181 8.59% FabArray::ParallelCopy() 1469 0.5687 0.9441 1.151 8.37% Castro::evaluate_source_change() 196 0.9483 0.9563 0.9795 7.12% Amr::init() 1 0.9071 0.91 0.9102 6.62% FabArray::ParallelCopy_finish() 1469 0.2995 0.6778 0.8858 6.44% Amr::initialInit() 1 0.8559 0.8588 0.8591 6.25% Castro::do_new_sources() 17 0.7644 0.8011 0.8363 6.08% MLCellLinOp::smooth() 1544 0.6312 0.6869 0.8264 6.01% AmrLevel::FillPatch() 187 0.6486 0.7487 0.8088 5.88% amrex::Copy() 766 0.7654 0.7722 0.7966 5.79% Castro::enforce_min_density() 79 0.7784 0.7839 0.7956 5.79% Amr::InitializeInit() 1 0.4595 0.4927 0.6488 4.72% Amr::defBaseLevel() 1 0.4594 0.4927 0.6488 4.72% FillPatchIterator::Initialize 187 0.458 0.5578 0.6198 4.51% FillPatchIterator::FillFromLevel0() 187 0.4258 0.5336 0.5979 4.35% FillPatchSingleLevel 187 0.4249 0.5327 0.597 4.34% MultiFab::AddProduct() 3598 0.5502 0.5546 0.5677 4.13% MLMG::mgVcycle_up::1 128 0.1075 0.2908 0.4725 3.44% MLCellLinOp::applyBC() 7148 0.3585 0.3957 0.4679 3.40% MLPoisson::Fsmooth() 6144 0.3287 0.3556 0.4438 3.23% Castro::Castro() 1 0.3922 0.3974 0.4163 3.03% Castro::initMFs() 1 0.3867 0.3915 0.4089 2.97% MLMG::addInterpCorrection() 768 0.06582 0.3316 0.4087 2.97% Amr::FinalizeInit() 1 0.2101 0.3661 0.3964 2.88% FabArray::setVal() 5232 0.3622 0.3667 0.3789 2.76% FabArray::FillBoundary() 6380 0.2794 0.3107 0.3565 2.59% Castro::reset_internal_energy(Fab) 160 0.3433 0.3464 0.3515 2.56% MultiFab::Subtract() 159 0.3288 0.3302 0.3333 2.42% MLMG::mgVcycle_up::2 128 0.01552 0.131 0.3092 2.25% FabArray::ParallelCopy_nowait() 1469 0.2449 0.2619 0.2837 2.06% MLMG::mgVcycle_down::0 128 0.2654 0.2703 0.2736 1.99% Gravity::fill_multipole_BCs() 19 0.2065 0.2459 0.2638 1.92% check_for_negative_density() 20 0.0929 0.165 0.258 1.88% Castro::expand_state() 20 0.2011 0.221 0.2521 1.83% Castro::initData() 1 0.06502 0.09022 0.2283 1.66% Amr::bldFineLevels() 1 0.03687 0.193 0.2235 1.63% MLCellLinOp::apply() 986 0.154 0.1663 0.2108 1.53% MLMG::mgVcycle_up::0 128 0.202 0.2036 0.2047 1.49% MLMG::mgVcycle_up::3 128 0.01552 0.02993 0.2033 1.48% FabArray::sum() 3598 0.1855 0.1871 0.1917 1.39% Amr::grid_places() 2 0.002626 0.1586 0.1889 1.37% AmrMesh::MakeNewGrids() 2 0.002567 0.1586 0.1889 1.37% TagBoxArray::mapPRD 2 0.000942 0.1209 0.1872 1.36% FillBoundary_finish() 6380 0.1357 0.1595 0.178 1.29% TagBoxArray::collate() 2 0.0002911 0.03637 0.1757 1.28% FillBoundary_nowait() 6380 0.1044 0.1339 0.1645 1.20% Castro::construct_old_source() 240 0.1545 0.1559 0.1582 1.15% MLMG:computeResOfCorrection() 768 0.1192 0.1292 0.1554 1.13% MLCellLinOp::correctionResidual() 768 0.1173 0.1274 0.1535 1.12% MLMG::mgVcycle_down::1 128 0.132 0.1368 0.1436 1.04% Castro::post_init() 1 0.138 0.1386 0.139 1.01% Castro::normalize_species() 79 0.1195 0.1205 0.1224 0.89% Castro::construct_new_source() 102 0.119 0.12 0.1209 0.88% MLPoisson::Fapply() 986 0.08541 0.08992 0.1093 0.79% Castro::construct_old_gravity() 20 0.05089 0.08176 0.1076 0.78% Gravity::get_old_grav_vector() 20 0.05076 0.08163 0.1075 0.78% Gravity::multilevel_solve_for_new_phi() 1 0.1057 0.1058 0.1058 0.77% Gravity::actual_multilevel_solve() 1 0.1057 0.1057 0.1058 0.77% MLCellLinOp::solutionResidual() 147 0.08669 0.09077 0.09563 0.70% Castro::construct_old_gravity_source() 40 0.09252 0.09333 0.09469 0.69% Amr::writePlotFile() 2 0.09373 0.09381 0.09382 0.68% MLMG::mgVcycle_down::2 128 0.01811 0.04972 0.08593 0.62% MLMG::mgVcycle_bottom 128 0.0004591 0.005677 0.08374 0.61% MLMG::actualBottomSolve() 128 0.0001898 0.005397 0.08345 0.61% MLCGSolver::bicgstab 8 0 0.005087 0.08139 0.59% FabArray::Saxpy() 111 0.07664 0.07712 0.07797 0.57% MLMG::computeResidual() 128 0.0732 0.07474 0.07643 0.56% Castro::create_source_corrector() 17 0.05336 0.06415 0.07329 0.53% Castro::construct_new_gravity_source() 17 0.06794 0.06866 0.06962 0.51% Castro::post_regrid() 2 0.06842 0.06874 0.06944 0.51% MLMG::mgVcycle_down::3 128 0.01763 0.02427 0.06767 0.49% Castro::derive() 18 0.06608 0.06676 0.06742 0.49% AmrLevel::derive() 18 0.06603 0.06671 0.06737 0.49% Castro::post_timestep() 1 0.05526 0.05991 0.06231 0.45% amrex::average_down 768 0.034 0.04116 0.05958 0.43% FabArray::Xpay() 915 0.04889 0.05106 0.05859 0.43% Castro::check_for_nan() 37 0.04402 0.0444 0.04507 0.33% MultiFab::contains_nan() 37 0.0439 0.04428 0.04492 0.33% MLMG::mgVcycle_down::4 128 0.01762 0.01947 0.04242 0.31% Castro::initialize_advance() 1 0.04032 0.04063 0.04197 0.31% MLPoisson::define() 19 0.03074 0.03567 0.04106 0.30% Castro::retry_advance_ctu() 11 0.03387 0.03561 0.03962 0.29% StateDataPhysBCFunct::() 187 0.0009686 0.01979 0.03781 0.28% StateData::FillBoundary(geom) 193 0 0.01861 0.03655 0.27% Amr::regrid() 1 0.03525 0.03555 0.03582 0.26% Gravity::get_new_grav_vector() 18 0.03415 0.03488 0.03571 0.26% MLMG::mgVcycle_up::4 128 0.01551 0.01687 0.03319 0.24% MLMG::mgVcycle_down::5 128 0.01757 0.01866 0.03071 0.22% MLCellLinOp::defineAuxData() 19 0.02039 0.02199 0.02714 0.20% Castro::apply_source_to_state() 37 0.02565 0.02578 0.02645 0.19% MLMG::mgVcycle_up::5 128 0.01563 0.0165 0.02603 0.19% VisMF::Write(FabArray) 4 0.01355 0.01448 0.02434 0.18% MLMG::prepareForSolve() 18 0.02016 0.02112 0.02372 0.17% MLMG::ResNormInf() 146 0.01584 0.01887 0.02233 0.16% Castro::sum_integrated_quantities() 2 0.02159 0.02174 0.02175 0.16% Castro::save_data_for_retry() 3 0.01748 0.01822 0.02118 0.15% Castro::problem_post_timestep() 2 0.009438 0.01823 0.02108 0.15% Castro::wd_update() 2 0.009189 0.0182 0.02106 0.15% MLMG::computeMLResidual() 18 0.01013 0.01565 0.02014 0.15% FabArrayBase::getFB() 6380 0.01362 0.0145 0.01834 0.13% FabArray::setDomainBndry() 187 0.0006764 0.008684 0.01664 0.12% FabArrayBase::getCPC() 2045 0.01101 0.0123 0.01581 0.12% FabArray::norminf() 244 0.01056 0.01095 0.01481 0.11% amrex::Add() 128 0.0136 0.01367 0.0142 0.10% Gravity::test_composite_phi() 1 0.01386 0.01391 0.01393 0.10% MLLinOp::define() 19 0.003762 0.007037 0.01044 0.08% MLLinOp::defineGrids() 19 0.003673 0.006946 0.01035 0.08% MLCGSolver::ParallelAllReduce 189 0 0.0006284 0.01005 0.07% MLMG::getGradSolution() 18 0.007188 0.007654 0.008903 0.06% Other 2323 0.07025 0.08728 0.1336 0.97% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Async Memory Usage: -------------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max -------------------------------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 6400 395 KiB 398 KiB 402 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 2157 B 2239 B 2380 B 89 KiB 89 KiB 89 KiB -------------------------------------------------------------------------------------------------------------------- Pinned Memory Usage: -------------------------------------------------------------------------------------------------------------- Name Nalloc AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max -------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 1 B 1 B 2 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 1573 KiB 1574 KiB 1574 KiB 2329 KiB 2329 KiB 2329 KiB Castro::construct_ctu_hydro_source() 6400 395 KiB 398 KiB 402 KiB 1176 KiB 1176 KiB 1176 KiB StateData::define() 192 1045 KiB 1046 KiB 1046 KiB 1164 KiB 1164 KiB 1164 KiB Castro::initMFs() 256 1024 KiB 1025 KiB 1025 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 668 KiB 668 KiB 668 KiB 803 KiB 803 KiB 803 KiB Amr::writePlotFile() 64 4391 B 4471 B 4990 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 50 KiB 50 KiB 51 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 13 KiB 15 KiB 17 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 430 KiB 430 KiB 431 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2779 B 4698 B 5838 B 225 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1399 B 1413 B 1429 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 18 KiB 19 KiB 19 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 130 KiB 130 KiB 130 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 127 KiB 127 KiB 127 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 10 KiB 10 KiB 12 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 1280 499 B 505 B 526 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 418 B 420 B 427 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 2157 B 2239 B 2380 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 20 B 21 B 38 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 33 KiB 33 KiB 33 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 53 B 53 B 53 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 3 B 3 B 4 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 111 B 114 B 117 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 168 B 272 B 359 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 31 KiB 31 KiB 31 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 237 B 242 B 245 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 1124 B 1133 B 1160 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 125 B 125 B 125 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 1792 B 1843 B 1869 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 0 B 0 B 1 B 16 KiB 16 KiB 16 KiB MLCellLinOp::defineAuxData() 11704 469 B 524 B 727 B 4608 B 5151 B 7152 B BndryData::define() 4864 488 B 489 B 490 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 27 B 40 B 54 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 0 B 22 B 35 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 27 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 6 B 8 B 648 B 648 B 648 B MLCGSolver::bicgstab 640 0 B 0 B 3 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 0 B 0 B 0 B 8 B 8 B 8 B ----------------------------------------------------------------------------------------------------------------------- AMReX (26.03-91-gf546b598fb75) finalized