Initializing AMReX (26.01-7-gc2a1abe7911b)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (26.01-7-gc2a1abe7911b) initialized Starting run at 08:42:26 UTC on 2026-01-09. Successfully read inputs file ... Castro git describe: 26.01-1-ge6f237912 AMReX git describe: 26.01-7-gc2a1abe79 Microphysics git describe: 26.01-4-g935d36be 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.173298993 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001166002 new finest: 0 Gravity::fill_multipole_BCs() time = 0.011377367 Gravity::multilevel_solve_for_new_phi() time = 0.109298305 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.011324312 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.051906675 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004220664 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.003669752 on level 0 Castro::do_old_sources() time = 0.008283237 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.062464745 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 = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Castro::construct_old_gravity() time = 0.005641965 on level 0 Castro::do_old_sources() time = 0.008584934 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.083724602 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, 177, 0 Invalid density = -0.1060744788 at index 20, 78, 0 Castro::construct_old_gravity() time = 0.003071903 on level 0 Castro::do_old_sources() time = 0.009494265 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.057929674 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.003093212 on level 0 Castro::do_old_sources() time = 0.007889924 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.056866151 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.333715645 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.01475744 Gravity::solve_for_phi() time = 0.09268698 on level 0 Castro::construct_new_gravity() time = 0.094949495 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.058663476 on level 0 Castro::construct_old_gravity() time = 0.003366231 on level 0 Castro::do_old_sources() time = 0.008026533 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.056765692 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.33372878 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.014721545 Gravity::solve_for_phi() time = 0.091606729 on level 0 Castro::construct_new_gravity() time = 0.093752249 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.057001255 on level 0 Castro::construct_old_gravity() time = 0.002785516 on level 0 Castro::do_old_sources() time = 0.008003764 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.056593132 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.341513985 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.014760892 Gravity::solve_for_phi() time = 0.09119622 on level 0 Castro::construct_new_gravity() time = 0.093343789 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.056953425 on level 0 Castro::construct_old_gravity() time = 0.002947519 on level 0 Castro::do_old_sources() time = 0.008018884 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.057061299 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.347847668 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.015624913 Gravity::solve_for_phi() time = 0.091781545 on level 0 Castro::construct_new_gravity() time = 0.093974371 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.057661148 on level 0 Castro::construct_old_gravity() time = 0.004742839 on level 0 Castro::do_old_sources() time = 0.008038306 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.056699156 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.349029722 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.015106897 Gravity::solve_for_phi() time = 0.091649237 on level 0 Castro::construct_new_gravity() time = 0.093848502 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.056874354 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.002824961 on level 0 Castro::do_old_sources() time = 0.007996466 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.056441418 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.344317011 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.015361638 Gravity::solve_for_phi() time = 0.09213692 on level 0 Castro::construct_new_gravity() time = 0.094673674 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.056688015 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.00275441 on level 0 Castro::do_old_sources() time = 0.007922043 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.056573891 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.350155079 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.016498395 Gravity::solve_for_phi() time = 0.096225213 on level 0 Castro::construct_new_gravity() time = 0.098379627 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.05756678 on level 0 Castro::construct_old_gravity() time = 0.002817065 on level 0 Castro::do_old_sources() time = 0.00830499 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.056540548 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.348337753 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.020228102 Gravity::solve_for_phi() time = 0.097865771 on level 0 Castro::construct_new_gravity() time = 0.100067554 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.05794271 on level 0 Castro::construct_old_gravity() time = 0.002681986 on level 0 Castro::do_old_sources() time = 0.007880851 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.056588471 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.353747614 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.021480057 Gravity::solve_for_phi() time = 0.10278416 on level 0 Castro::construct_new_gravity() time = 0.104930691 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.057075909 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.0031317 on level 0 Castro::do_old_sources() time = 0.007825017 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.056378569 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.349500403 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.02080788 Gravity::solve_for_phi() time = 0.097435753 on level 0 Castro::construct_new_gravity() time = 0.099581145 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.056686934 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.002724925 on level 0 Castro::do_old_sources() time = 0.007927409 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.056318726 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.350356947 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.015226057 Gravity::solve_for_phi() time = 0.09177127 on level 0 Castro::construct_new_gravity() time = 0.093988102 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.056647141 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.002883659 on level 0 Castro::do_old_sources() time = 0.007908726 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.057327768 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.343988815 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.016276614 Gravity::solve_for_phi() time = 0.093682137 on level 0 Castro::construct_new_gravity() time = 0.095899412 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.056719248 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.003115026 on level 0 Castro::do_old_sources() time = 0.008075341 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.056494616 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.352319245 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.015985516 Gravity::solve_for_phi() time = 0.093224515 on level 0 Castro::construct_new_gravity() time = 0.095371546 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.058123219 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.002866816 on level 0 Castro::do_old_sources() time = 0.007826222 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.056724483 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.344231052 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.0157385 Gravity::solve_for_phi() time = 0.093029318 on level 0 Castro::construct_new_gravity() time = 0.095171054 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.057039081 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.002761212 on level 0 Castro::do_old_sources() time = 0.007887345 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.056279227 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.36186022 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.015926887 Gravity::solve_for_phi() time = 0.094020422 on level 0 Castro::construct_new_gravity() time = 0.096173781 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.099762692 on level 0 Castro::construct_old_gravity() time = 0.011605939 on level 0 Castro::do_old_sources() time = 0.008221333 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.058466047 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.361830066 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.015902977 Gravity::solve_for_phi() time = 0.092922593 on level 0 Castro::construct_new_gravity() time = 0.095073198 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.056642448 on level 0 Castro::construct_old_gravity() time = 0.003596442 on level 0 Castro::do_old_sources() time = 0.007954048 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.05684357 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.351058372 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002597363434 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.016204391 Gravity::solve_for_phi() time = 0.092711988 on level 0 Castro::construct_new_gravity() time = 0.094890575 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.057339617 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: 12.67197047 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.05156993 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004263384 seconds Ending run at 08:42:41 UTC on 2026-01-09. Run time = 15.13316982 Run time without initialization = 12.72877477 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]: 15.1267 ... 15.1281 ... 15.1332 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.8077 5.8439 5.9397 39.25% Amr::InitAmr() 1 1.478 1.479 1.484 9.81% Castro::computeTemp() 80 0.9233 1.048 1.233 8.15% amrex::Copy() 766 0.858 0.864 0.9001 5.95% FabArray::ParallelCopy_finish() 1469 0.2641 0.6366 0.8313 5.49% MultiFab::AddProduct() 3598 0.5862 0.5909 0.6056 4.00% Castro::reset_internal_energy(Fab) 160 0.4497 0.4549 0.4747 3.14% MLPoisson::Fsmooth() 6144 0.354 0.3799 0.4644 3.07% Castro::initMFs() 1 0.3826 0.3932 0.4146 2.74% FabArray::setVal() 5232 0.3953 0.3992 0.4133 2.73% MultiFab::Subtract() 159 0.3587 0.3596 0.3657 2.42% Castro::do_advance_ctu() 20 0.2945 0.2979 0.3079 2.03% Gravity::fill_multipole_BCs() 19 0.2165 0.2734 0.2978 1.97% FabArray::ParallelCopy_nowait() 1469 0.2356 0.2557 0.2802 1.85% check_for_negative_density() 20 0.05486 0.165 0.22 1.45% FabArray::sum() 3598 0.1976 0.1989 0.2028 1.34% Castro::initData() 1 0.003371 0.02863 0.1639 1.08% FillBoundary_nowait() 6380 0.08814 0.1161 0.1491 0.99% TagBoxArray::collate() 2 0.0002711 0.0273 0.1378 0.91% FillBoundary_finish() 6380 0.1069 0.1222 0.137 0.91% MLCellLinOp::applyBC() 7148 0.0673 0.08237 0.1337 0.88% Castro::normalize_species() 79 0.1151 0.1161 0.1223 0.81% MLPoisson::Fapply() 986 0.09555 0.09964 0.1164 0.77% Castro::construct_old_gravity_source() 40 0.1013 0.102 0.1045 0.69% FabArray::Saxpy() 111 0.08532 0.08588 0.08711 0.58% Castro::construct_new_gravity_source() 17 0.07307 0.07378 0.07644 0.51% FabArray::Xpay() 915 0.05132 0.0533 0.05971 0.39% MultiFab::contains_nan() 37 0.04881 0.0491 0.04978 0.33% AmrLevel::derive() 18 0.04438 0.04467 0.04538 0.30% MLMG::addInterpCorrection() 768 0.03393 0.03626 0.04444 0.29% StateData::FillBoundary(geom) 193 0 0.02002 0.0388 0.26% amrex::average_down 768 0.02519 0.02718 0.03453 0.23% Castro::reset_internal_energy(MultiFab) 80 0.02447 0.02469 0.02509 0.17% VisMF::Write(FabArray) 4 0.01416 0.01531 0.02502 0.17% Castro::construct_old_source() 240 0.02047 0.02061 0.02092 0.14% MLCellLinOp::defineAuxData() 19 0.01564 0.01649 0.02082 0.14% FabArray::FillBoundary() 6380 0.0165 0.01686 0.02013 0.13% Castro::enforce_min_density() 79 0.0193 0.01951 0.02004 0.13% FabArray::setDomainBndry() 187 0.0006619 0.00915 0.01861 0.12% MLCGSolver::bicgstab 8 0 0.001084 0.01735 0.11% Amr::writePlotFile() 2 0.006776 0.01592 0.01712 0.11% FabArrayBase::getFB() 6380 0.01353 0.0139 0.01711 0.11% Castro::construct_new_source() 102 0.01614 0.01622 0.01648 0.11% Castro::evaluate_source_change() 196 0.01464 0.01499 0.01575 0.10% FillPatchIterator::Initialize 187 0.01518 0.01539 0.01565 0.10% FabArray::norminf() 244 0.01104 0.01142 0.01485 0.10% amrex::Add() 128 0.01457 0.01461 0.01471 0.10% Castro::wd_update() 2 0.005271 0.01048 0.01415 0.09% MLMG::prepareForSolve() 18 0.01034 0.01066 0.01226 0.08% MLMG::solve() 18 0.00116 0.005887 0.01188 0.08% MLCellLinOp::smooth() 1544 0.01061 0.01086 0.01159 0.08% Gravity::actual_solve_with_mlmg() 19 0.007455 0.007906 0.01 0.07% FabArrayBase::getCPC() 2045 0.00848 0.008643 0.009307 0.06% MLMG::ResNormInf() 146 0.002949 0.005078 0.006935 0.05% Gravity::get_old_grav_vector() 20 0.006526 0.006621 0.006772 0.04% Gravity::get_new_grav_vector() 18 0.006136 0.006195 0.00643 0.04% MLCellLinOp::apply() 986 0.002371 0.002622 0.005682 0.04% Castro::initialize_do_advance() 20 0.004714 0.004877 0.005135 0.03% FabArray::ParallelCopy() 1469 0.004191 0.004293 0.004563 0.03% AmrLevel::FillPatch() 187 0.003064 0.003166 0.003351 0.02% Castro::sum_integrated_quantities() 2 0.001556 0.002338 0.002491 0.02% MLMG::mgVcycle() 128 0.002226 0.002265 0.002416 0.02% MLCellLinOp::correctionResidual() 768 0.002217 0.002267 0.002385 0.02% MLMG:computeResOfCorrection() 768 0.00175 0.001785 0.001856 0.01% Castro::save_data_for_retry() 3 0.001313 0.001328 0.001345 0.01% MLMG::actualBottomSolve() 128 0.0001809 0.0002552 0.001268 0.01% StateDataPhysBCFunct::() 187 0.0009449 0.001136 0.001252 0.01% Castro::subcycle_advance_ctu() 1 0.000249 0.0003318 0.001108 0.01% FillPatchIterator::FillFromLevel0() 187 0.0008869 0.0009067 0.0009329 0.01% main() 1 0.0006976 0.0007795 0.0008717 0.01% Gravity::solve_for_phi() 17 0.0006667 0.000695 0.0007898 0.01% FillPatchSingleLevel 187 0.0006607 0.0006722 0.000685 0.00% MLMG::mgVcycle_down::0 128 0.0006429 0.0006572 0.00068 0.00% MLMG::mgVcycle_down::2 128 0.0005791 0.000613 0.0006516 0.00% MLMG::mgVcycle_down::1 128 0.0006174 0.0006266 0.0006441 0.00% MLMG::mgVcycle_down::4 128 0.000545 0.0005597 0.0006428 0.00% MLMG::mgVcycle_down::3 128 0.0005481 0.0005783 0.0006411 0.00% MLMG::mgVcycle_down::5 128 0.0005713 0.0005843 0.0006284 0.00% Castro::Castro() 1 0.0004106 0.0004687 0.0005835 0.00% AmrMesh::MakeNewGrids() 2 0.0005379 0.0005505 0.0005704 0.00% Castro::clean_state() 79 0.0004625 0.0004874 0.0005414 0.00% MLCellLinOp::solutionResidual() 147 0.0005051 0.0005156 0.0005331 0.00% Castro::do_old_sources() 40 0.0004529 0.0004678 0.0004834 0.00% MLMG::oneIter() 128 0.0004228 0.0004351 0.0004571 0.00% MLMG::mgVcycle_up::1 128 0.0004021 0.0004094 0.0004217 0.00% MLMG::mgVcycle_up::0 128 0.0004058 0.0004105 0.0004207 0.00% MLMG::mgVcycle_up::2 128 0.0003833 0.0003983 0.0004145 0.00% MLMG::mgVcycle_up::3 128 0.0003718 0.0003809 0.0004017 0.00% MLMG::mgVcycle_up::4 128 0.0003761 0.0003834 0.0003977 0.00% MLMG::mgVcycle_up::5 128 0.000363 0.0003689 0.0003843 0.00% Castro::construct_new_gravity() 17 0.0001663 0.0001853 0.0003787 0.00% MLMG::computeResidual() 128 0.000286 0.0002935 0.0003097 0.00% MLMG::mgVcycle_bottom 128 0.0002696 0.0002795 0.0002972 0.00% Castro::problem_post_timestep() 2 1.165e-05 3.032e-05 0.0002757 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002129 0.0002276 0.0002681 0.00% Castro::do_new_sources() 17 0.0002134 0.0002202 0.0002306 0.00% Amr::coarseTimeStep() 1 0.0001494 0.0001767 0.0001832 0.00% Castro::check_for_nan() 37 0.0001159 0.0001256 0.0001474 0.00% MLPoisson::define() 19 0.0001238 0.0001329 0.0001469 0.00% Castro::initialize_advance() 1 0.0001093 0.0001208 0.000141 0.00% Amr::defBaseLevel() 1 9.624e-05 0.0001084 0.0001407 0.00% Castro::retry_advance_ctu() 11 7.191e-05 8.322e-05 0.0001397 0.00% Castro::apply_source_to_state() 37 0.0001152 0.0001224 0.0001363 0.00% TagBoxArray::mapPRD 2 0.0001151 0.0001212 0.0001302 0.00% Castro::construct_old_gravity() 20 0.0001135 0.0001219 0.0001276 0.00% Gravity::test_composite_phi() 1 0.0001042 0.0001099 0.0001273 0.00% MLLinOp::define() 19 8.189e-05 8.782e-05 9.555e-05 0.00% Castro::expand_state() 20 5.746e-05 6.31e-05 7.75e-05 0.00% Amr::grid_places() 2 5.241e-05 5.827e-05 7.684e-05 0.00% Amr::FinalizeInit() 1 8.407e-06 1.553e-05 6.236e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.771e-05 5.875e-05 5.973e-05 0.00% Castro::create_source_corrector() 17 5.132e-05 5.537e-05 5.896e-05 0.00% Amr::regrid() 1 1.662e-05 2.512e-05 5.753e-05 0.00% MLMG::computeMLResidual() 18 4.932e-05 5.233e-05 5.594e-05 0.00% Castro::derive() 18 4.49e-05 4.678e-05 4.9e-05 0.00% Gravity::actual_multilevel_solve() 1 3.106e-05 3.199e-05 3.279e-05 0.00% Amr::timeStep() 1 1.553e-05 1.845e-05 2.457e-05 0.00% Amr::InitializeInit() 1 8.531e-06 1.45e-05 2.336e-05 0.00% Castro::post_init() 1 1.165e-05 1.314e-05 1.909e-05 0.00% Castro::post_regrid() 2 1.069e-05 1.165e-05 1.521e-05 0.00% Castro::advance() 1 8.175e-06 9.192e-06 1.291e-05 0.00% Amr::bldFineLevels() 1 7.262e-06 8.151e-06 1.249e-05 0.00% Amr::init() 1 6.146e-06 7.298e-06 1.025e-05 0.00% Castro::post_timestep() 1 5.788e-06 6.337e-06 7.727e-06 0.00% Amr::initialInit() 1 3.29e-06 3.657e-06 5.306e-06 0.00% Other 2549 0.04972 0.06289 0.1002 0.66% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 15.13 15.13 15.13 100.00% Amr::coarseTimeStep() 1 12.73 12.73 12.73 84.11% Amr::timeStep() 1 12.67 12.67 12.67 83.74% Castro::advance() 1 12.61 12.61 12.62 83.37% Castro::subcycle_advance_ctu() 1 12.56 12.57 12.57 83.06% Castro::do_advance_ctu() 20 12.52 12.53 12.53 82.79% Castro::construct_ctu_hydro_source() 20 6.703 6.736 6.79 44.87% Castro::clean_state() 79 3.352 3.487 3.664 24.21% Castro::computeTemp() 80 2.416 2.548 2.73 18.04% Gravity::solve_phi_with_mlmg() 19 1.626 1.687 1.716 11.34% Castro::construct_new_gravity() 17 1.543 1.604 1.632 10.79% Gravity::solve_for_phi() 17 1.505 1.566 1.595 10.54% Castro::reset_internal_energy(MultiFab) 80 1.49 1.5 1.547 10.22% Amr::InitAmr() 1 1.478 1.479 1.484 9.81% Gravity::actual_solve_with_mlmg() 19 1.405 1.41 1.414 9.35% Castro::initialize_do_advance() 20 1.31 1.355 1.388 9.17% MLMG::solve() 18 1.346 1.354 1.355 8.95% Castro::do_old_sources() 40 1.255 1.291 1.314 8.68% MLMG::oneIter() 128 1.213 1.214 1.216 8.04% MLMG::mgVcycle() 128 1.198 1.199 1.201 7.94% FabArray::ParallelCopy() 1469 0.5489 0.9181 1.112 7.35% Castro::evaluate_source_change() 196 1.015 1.023 1.038 6.86% Castro::do_new_sources() 17 0.8672 0.9241 0.9749 6.44% Amr::init() 1 0.92 0.92 0.92 6.08% amrex::Copy() 766 0.858 0.864 0.9001 5.95% Castro::enforce_min_density() 79 0.844 0.8505 0.879 5.81% Amr::initialInit() 1 0.8633 0.8633 0.8634 5.71% AmrLevel::FillPatch() 187 0.6714 0.764 0.8478 5.60% FabArray::ParallelCopy_finish() 1469 0.2641 0.6366 0.8313 5.49% MLCellLinOp::smooth() 1544 0.6218 0.6732 0.8119 5.37% Amr::InitializeInit() 1 0.4658 0.5024 0.6378 4.21% Amr::defBaseLevel() 1 0.4658 0.5024 0.6377 4.21% FillPatchIterator::Initialize 187 0.4594 0.5518 0.6363 4.20% MultiFab::AddProduct() 3598 0.5862 0.5909 0.6056 4.00% FillPatchIterator::FillFromLevel0() 187 0.4255 0.5272 0.604 3.99% FillPatchSingleLevel 187 0.4246 0.5263 0.6031 3.99% Castro::reset_internal_energy(Fab) 160 0.4497 0.4549 0.4747 3.14% MLPoisson::Fsmooth() 6144 0.354 0.3799 0.4644 3.07% MLMG::mgVcycle_up::1 128 0.1041 0.2806 0.4563 3.02% MLCellLinOp::applyBC() 7148 0.3168 0.3522 0.4254 2.81% Castro::Castro() 1 0.3883 0.3991 0.4221 2.79% Castro::initMFs() 1 0.3828 0.3933 0.4147 2.74% FabArray::setVal() 5232 0.3953 0.3992 0.4133 2.73% MLMG::addInterpCorrection() 768 0.06714 0.3267 0.3987 2.63% Amr::FinalizeInit() 1 0.2256 0.3609 0.3976 2.63% MultiFab::Subtract() 159 0.3587 0.3596 0.3657 2.42% FabArray::FillBoundary() 6380 0.2382 0.2698 0.3106 2.05% MLMG::mgVcycle_up::2 128 0.01543 0.1294 0.3033 2.00% FabArray::ParallelCopy_nowait() 1469 0.2564 0.2772 0.3026 2.00% Gravity::fill_multipole_BCs() 19 0.2187 0.2756 0.3 1.98% Castro::expand_state() 20 0.2019 0.2627 0.2909 1.92% MLMG::mgVcycle_down::0 128 0.2742 0.2761 0.2793 1.85% Castro::initData() 1 0.07566 0.1015 0.2366 1.56% check_for_negative_density() 20 0.05486 0.165 0.22 1.45% Amr::bldFineLevels() 1 0.04236 0.1776 0.214 1.41% MLCellLinOp::apply() 986 0.1563 0.17 0.2097 1.39% MLMG::mgVcycle_up::0 128 0.2067 0.2074 0.2079 1.37% FabArray::sum() 3598 0.1976 0.1989 0.2028 1.34% MLMG::mgVcycle_up::3 128 0.01532 0.02952 0.2019 1.33% Amr::grid_places() 2 0.002463 0.1378 0.1745 1.15% AmrMesh::MakeNewGrids() 2 0.00241 0.1378 0.1744 1.15% Castro::construct_old_source() 240 0.169 0.1698 0.1737 1.15% TagBoxArray::mapPRD 2 0.0007792 0.1091 0.1728 1.14% FillBoundary_nowait() 6380 0.1021 0.1307 0.1641 1.08% MLMG:computeResOfCorrection() 768 0.1202 0.1301 0.1545 1.02% MLCellLinOp::correctionResidual() 768 0.1184 0.1283 0.1527 1.01% Castro::post_init() 1 0.1427 0.1437 0.144 0.95% TagBoxArray::collate() 2 0.0002711 0.0273 0.1378 0.91% FillBoundary_finish() 6380 0.1069 0.1222 0.137 0.91% Castro::construct_new_source() 102 0.1294 0.1302 0.1327 0.88% MLMG::mgVcycle_down::1 128 0.1243 0.1281 0.1323 0.87% Castro::normalize_species() 79 0.1151 0.1161 0.1223 0.81% Castro::create_source_corrector() 17 0.05626 0.0807 0.1175 0.78% MLPoisson::Fapply() 986 0.09555 0.09964 0.1164 0.77% Gravity::multilevel_solve_for_new_phi() 1 0.1092 0.1092 0.1093 0.72% Gravity::actual_multilevel_solve() 1 0.1091 0.1092 0.1092 0.72% Castro::construct_old_gravity_source() 40 0.1013 0.102 0.1045 0.69% Amr::writePlotFile() 2 0.1039 0.104 0.104 0.69% MLCellLinOp::solutionResidual() 147 0.08863 0.09577 0.1029 0.68% FabArray::Saxpy() 111 0.08532 0.08588 0.08711 0.58% MLMG::mgVcycle_bottom 128 0.0004546 0.005486 0.08077 0.53% MLMG::actualBottomSolve() 128 0.0001809 0.005206 0.08048 0.53% Castro::post_regrid() 2 0.07905 0.07938 0.08034 0.53% MLMG::mgVcycle_down::2 128 0.01788 0.04699 0.07925 0.52% MLCGSolver::bicgstab 8 0 0.004904 0.07847 0.52% MLMG::computeResidual() 128 0.07665 0.0771 0.07758 0.51% Castro::construct_new_gravity_source() 17 0.07307 0.07378 0.07644 0.51% Castro::derive() 18 0.07496 0.07535 0.07627 0.50% AmrLevel::derive() 18 0.07491 0.07531 0.07623 0.50% MLMG::mgVcycle_down::3 128 0.01725 0.02363 0.06785 0.45% Castro::construct_old_gravity() 20 0.04237 0.0558 0.06512 0.43% Gravity::get_old_grav_vector() 20 0.04225 0.05568 0.065 0.43% Castro::post_timestep() 1 0.05605 0.05963 0.06271 0.41% FabArray::Xpay() 915 0.05132 0.0533 0.05971 0.39% amrex::average_down 768 0.03501 0.0415 0.05898 0.39% Castro::check_for_nan() 37 0.04894 0.04922 0.04992 0.33% MultiFab::contains_nan() 37 0.04881 0.0491 0.04978 0.33% Castro::initialize_advance() 1 0.04582 0.04621 0.04743 0.31% MLPoisson::define() 19 0.03349 0.0385 0.04419 0.29% MLMG::mgVcycle_down::4 128 0.01726 0.0189 0.04227 0.28% Amr::regrid() 1 0.04069 0.04095 0.04181 0.28% Castro::retry_advance_ctu() 11 0.03794 0.03849 0.04037 0.27% StateDataPhysBCFunct::() 187 0.0009449 0.02115 0.04004 0.26% StateData::FillBoundary(geom) 193 0 0.02002 0.0388 0.26% Gravity::get_new_grav_vector() 18 0.03671 0.03734 0.0379 0.25% MLMG::mgVcycle_up::4 128 0.01538 0.01655 0.03343 0.22% MLMG::mgVcycle_down::5 128 0.01725 0.01812 0.02987 0.20% Castro::apply_source_to_state() 37 0.02861 0.02888 0.02978 0.20% MLCellLinOp::defineAuxData() 19 0.02023 0.02176 0.02842 0.19% MLMG::prepareForSolve() 18 0.0211 0.0221 0.02659 0.18% MLMG::mgVcycle_up::5 128 0.0154 0.0161 0.02566 0.17% MLMG::computeMLResidual() 18 0.01112 0.01824 0.02533 0.17% VisMF::Write(FabArray) 4 0.01416 0.01531 0.02502 0.17% Castro::save_data_for_retry() 3 0.01942 0.01979 0.02158 0.14% Castro::sum_integrated_quantities() 2 0.02112 0.02128 0.0213 0.14% FabArray::setDomainBndry() 187 0.0006619 0.00915 0.01861 0.12% FabArrayBase::getFB() 6380 0.01387 0.01459 0.0186 0.12% Castro::problem_post_timestep() 2 0.009387 0.01405 0.01733 0.11% Castro::wd_update() 2 0.009111 0.01402 0.01732 0.11% MLMG::ResNormInf() 146 0.01277 0.01479 0.0167 0.11% FabArrayBase::getCPC() 2045 0.01081 0.01209 0.01615 0.11% Gravity::test_composite_phi() 1 0.01496 0.01504 0.01507 0.10% FabArray::norminf() 244 0.01104 0.01142 0.01485 0.10% amrex::Add() 128 0.01457 0.01461 0.01471 0.10% MLLinOp::define() 19 0.004994 0.0102 0.01393 0.09% Other 2549 0.08475 0.1034 0.1497 0.99% -------------------------------------------------------------------------------------------- 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 408 KiB 411 KiB 417 KiB 1176 KiB 1176 KiB 1176 KiB Castro::do_advance_ctu() 960 2204 B 2244 B 2304 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 4 B 5 B 6 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 1588 KiB 1589 KiB 1589 KiB 2329 KiB 2329 KiB 2329 KiB Castro::construct_ctu_hydro_source() 6400 408 KiB 411 KiB 417 KiB 1176 KiB 1176 KiB 1176 KiB StateData::define() 192 1049 KiB 1049 KiB 1050 KiB 1164 KiB 1164 KiB 1164 KiB Castro::initMFs() 256 1028 KiB 1029 KiB 1029 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 674 KiB 675 KiB 675 KiB 803 KiB 803 KiB 803 KiB Amr::writePlotFile() 64 4485 B 4547 B 5046 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 51 KiB 51 KiB 53 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 13 KiB 15 KiB 16 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 435 KiB 435 KiB 435 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 2394 B 4870 B 6189 B 225 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1452 B 1461 B 1482 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 18 KiB 18 KiB 19 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 131 KiB 131 KiB 131 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 128 KiB 128 KiB 128 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 9357 B 9898 B 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 1280 498 B 499 B 509 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 419 B 420 B 426 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 2204 B 2244 B 2304 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 18 B 20 B 36 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 51 B 51 B 52 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 3 B 4 B 5 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 109 B 111 B 113 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 127 B 167 B 196 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 225 B 233 B 246 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 1093 B 1101 B 1118 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 118 B 118 B 118 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 1628 B 1693 B 1725 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 423 B 472 B 655 B 4608 B 5151 B 7152 B BndryData::define() 4864 439 B 441 B 442 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 21 B 30 B 40 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 0 B 18 B 29 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 15 B 24 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 5 B 7 B 648 B 648 B 648 B MLCGSolver::bicgstab 640 0 B 0 B 2 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 0 B 0 B 0 B 8 B 8 B 8 B ----------------------------------------------------------------------------------------------------------------------- AMReX (26.01-7-gc2a1abe7911b) finalized