MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (22.04-39-g50466d11e7b2) initialized Starting run at 23:50:48 UTC on 2022-04-30. Successfully read inputs file ... Castro git describe: 22.04-49-g7653ce36f AMReX git describe: 22.04-39-g50466d11e Microphysics git describe: 22.04-3-g3c498521 reading extern runtime parameters ... 13 Species: He4 C12 O16 Ne20 Mg24 Si28 S32 Ar36 Ca40 Ti44 Cr48 Fe52 Ni56 Successfully read inputs file ... 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 Warning: the following castro.* parameters are ignored castro.center = 0.0e0 0.0e0 castro.dtnuc_mode = 1 Initializing the data at level 0 Done initializing the level 0 data grid_places() time: 0.141806184 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.000936651 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010612078 Gravity::fill_multipole_BCs() time = 0.010626709 ... norm of composite residual at level 0 1.193955541e-11 TIME= 0 MASS = 2.549996568e+33 TIME= 0 XMOM = 0 TIME= 0 YMOM = 0 TIME= 0 ZMOM = 0 TIME= 0 ANG MOM X = 0 TIME= 0 ANG MOM Y = 0 TIME= 0 ANG MOM Z = -3.115378115e+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 = 9.947183943e-17 TIME= 0 CENTER OF MASS X-VEL = 0 TIME= 0 CENTER OF MASS Y-LOC = -2.5098073e-31 TIME= 0 CENTER OF MASS Y-VEL = 0 TIME= 0 CENTER OF MASS Z-LOC = 0 TIME= 0 CENTER OF MASS Z-VEL = 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.052522034 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.006592996 seconds [Level 0 step 1] ADVANCE 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. 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.498822529e+28 zmom added: 0 eint added: 0 ener added: 6.931398274e+47 Castro::do_old_sources() time = 0.061777693 ... Entering construct_ctu_hydro_source() 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.239071149 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.7694089 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477042 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035501 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 ... Leaving construct_ctu_hydro_source() Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035501 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477042 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.7694089 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.239071149 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_ctu_hydro_source() time = 0.257906447 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.767071355 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.767071355 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 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. 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.249411264e+28 zmom added: 0 eint added: 0 ener added: 3.465699137e+47 Castro::do_old_sources() time = 0.060927155 ... Entering construct_ctu_hydro_source() Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.251623181 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 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. Contributions to the state from old-time sources: mass added: 0 xmom added: -2.93534348e+40 ymom added: 1.624705632e+28 zmom added: 0 eint added: 0 ener added: 1.732849568e+47 Castro::do_old_sources() time = 0.057450373 ... Entering construct_ctu_hydro_source() Invalid density = -0.1060744788 at index 20, 78, 0 ... Leaving construct_ctu_hydro_source() Invalid density = -0.1060744788 at index 20, 177, 0 Castro::construct_ctu_hydro_source() time = 0.253131344 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. Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.123528161e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.055443031 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.254188184 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015540256 Gravity::solve_for_phi() time = 0.083452181 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.674192048e+44 ener added: 3.847280946e+42 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.752957884e+35 ymom added: -3.586520465e+27 zmom added: 0 eint added: 0 ener added: 6.209484931e+43 Castro::do_new_sources() time = 0.059985618 Ending SDC iteration 1 of 3. Beginning SDC iteration 2 of 3. Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.123528161e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.055387593 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.253675812 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014731253 Gravity::solve_for_phi() time = 0.081451954 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.948673093e+46 ener added: 3.718993455e+42 Contributions to the state from new-time sources: mass added: 0 xmom added: -4.993343564e+36 ymom added: -3.766130952e+27 zmom added: 0 eint added: 0 ener added: 6.97362419e+46 Castro::do_new_sources() time = 0.060183185 Ending SDC iteration 2 of 3. Beginning SDC iteration 3 of 3. Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.123528161e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.055468478 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.259185201 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.013608174 Gravity::solve_for_phi() time = 0.080813807 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.948673093e+46 ener added: 3.718993455e+42 Contributions to the state from new-time sources: mass added: 0 xmom added: -4.993343564e+36 ymom added: -3.766130952e+27 zmom added: 0 eint added: 0 ener added: 6.97362419e+46 Castro::do_new_sources() time = 0.058646775 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. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 4.503787228e+25 ener added: 4.503787228e+25 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.994270532e+46 ener added: 5.341284592e+40 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.468670409e+40 ymom added: 5.931292302e+26 zmom added: 0 eint added: 0 ener added: 8.719870913e+46 Castro::do_old_sources() time = 0.055089812 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.260791889 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.997068055e+46 ener added: 2.415044041e+42 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014302817 Gravity::solve_for_phi() time = 0.079683666 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.543114243e+44 ener added: 4.768138875e+43 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.765646947e+35 ymom added: 1.797592414e+26 zmom added: 0 eint added: 0 ener added: 5.380039006e+43 Castro::do_new_sources() time = 0.058459386 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 4.503787228e+25 ener added: 4.503787228e+25 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.994270532e+46 ener added: 5.341284592e+40 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.468670409e+40 ymom added: 5.931292302e+26 zmom added: 0 eint added: 0 ener added: 8.719870913e+46 Castro::do_old_sources() time = 0.072385438 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.283896493 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.317346246e+44 ener added: 3.120883219e+42 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015147786 Gravity::solve_for_phi() time = 0.083023061 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.957879842e+46 ener added: 4.564181618e+43 Contributions to the state from new-time sources: mass added: 0 xmom added: -5.001437306e+36 ymom added: 1.471489396e+25 zmom added: 0 eint added: 0 ener added: 6.980886509e+46 Castro::do_new_sources() time = 0.060240701 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. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -7.002879254e+46 ener added: 6.901118887e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.469670696e+40 ymom added: 6.16552168e+26 zmom added: 0 eint added: 0 ener added: 8.772022862e+46 Castro::do_old_sources() time = 0.054806809 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.262318906 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.012688492e+46 ener added: 2.110657671e+43 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014021353 Gravity::solve_for_phi() time = 0.079447138 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.326633224e+44 ener added: 1.174460404e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.711471844e+35 ymom added: 1.66063198e+26 zmom added: 0 eint added: 0 ener added: 2.826321793e+43 Castro::do_new_sources() time = 0.057952355 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -7.002879254e+46 ener added: 6.901118887e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.469670696e+40 ymom added: 6.16552168e+26 zmom added: 0 eint added: 0 ener added: 8.772022862e+46 Castro::do_old_sources() time = 0.062409832 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.262089027 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.090986835e+44 ener added: 2.655092076e+43 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01349832 Gravity::solve_for_phi() time = 0.079036753 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.968859957e+46 ener added: 1.10301368e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -5.004822123e+36 ymom added: -8.771795742e+24 zmom added: 0 eint added: 0 ener added: 6.987657644e+46 Castro::do_new_sources() time = 0.058416104 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. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -7.011724564e+46 ener added: 1.143975262e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.47067166e+40 ymom added: 6.065029721e+26 zmom added: 0 eint added: 0 ener added: 8.82072643e+46 Castro::do_old_sources() time = 0.0549174 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.257322427 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.025000527e+46 ener added: 2.183068315e+43 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015400717 Gravity::solve_for_phi() time = 0.080963932 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.651474965e+44 ener added: 1.477761377e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.60533052e+35 ymom added: 1.690843319e+26 zmom added: 0 eint added: 0 ener added: -1.190247427e+43 Castro::do_new_sources() time = 0.058781762 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -7.011724564e+46 ener added: 1.143975262e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.47067166e+40 ymom added: 6.065029721e+26 zmom added: 0 eint added: 0 ener added: 8.82072643e+46 Castro::do_old_sources() time = 0.055204869 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.261472185 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.213015184e+44 ener added: 2.894439872e+43 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01389422 Gravity::solve_for_phi() time = 0.079369883 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.97482432e+46 ener added: 1.388215779e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -5.00305216e+36 ymom added: 1.271497176e+24 zmom added: 0 eint added: 0 ener added: 6.99302976e+46 Castro::do_new_sources() time = 0.059013849 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. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.25179549e+29 ener added: 1.25179549e+29 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -7.020672591e+46 ener added: 1.490478693e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.471672271e+40 ymom added: 6.057473935e+26 zmom added: 0 eint added: 0 ener added: 8.865830705e+46 Castro::do_old_sources() time = 0.055671259 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.258810916 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.046400343e+46 ener added: 2.16073557e+43 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014782834 Gravity::solve_for_phi() time = 0.082152133 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -2.248919276e+44 ener added: 1.60656791e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.450686215e+35 ymom added: 1.697312961e+26 zmom added: 0 eint added: 0 ener added: -6.185767806e+43 Castro::do_new_sources() time = 0.0607876 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.25179549e+29 ener added: 1.25179549e+29 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -7.020672591e+46 ener added: 1.490478693e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.471672271e+40 ymom added: 6.057473935e+26 zmom added: 0 eint added: 0 ener added: 8.865830705e+46 Castro::do_old_sources() time = 0.057651201 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.263934876 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.824816589e+44 ener added: 2.945461224e+43 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014566628 Gravity::solve_for_phi() time = 0.0800974 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.977981871e+46 ener added: 1.507699921e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -4.996348783e+36 ymom added: 3.069538214e+24 zmom added: 0 eint added: 0 ener added: 6.997250022e+46 Castro::do_new_sources() time = 0.059026262 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. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.258977451e+28 ener added: 6.258977451e+28 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -7.029400513e+46 ener added: 2.040436118e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.472671541e+40 ymom added: 6.126231591e+26 zmom added: 0 eint added: 0 ener added: 8.906898514e+46 Castro::do_old_sources() time = 0.055748606 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.257219296 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.067953295e+46 ener added: 2.185941597e+43 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014832772 Gravity::solve_for_phi() time = 0.081141028 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -3.036038551e+44 ener added: 1.710131682e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.236317743e+35 ymom added: 1.665059198e+26 zmom added: 0 eint added: 0 ener added: -1.282743578e+44 Castro::do_new_sources() time = 0.05871225 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.258977451e+28 ener added: 6.258977451e+28 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -7.029400513e+46 ener added: 2.040436118e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.472671541e+40 ymom added: 6.126231591e+26 zmom added: 0 eint added: 0 ener added: 8.906898514e+46 Castro::do_old_sources() time = 0.054977985 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.2635153 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 3.23178587e+44 ener added: 2.927676839e+43 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014243945 Gravity::solve_for_phi() time = 0.079744856 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.979285299e+46 ener added: 1.614875382e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -4.983606539e+36 ymom added: -1.663453594e+24 zmom added: 0 eint added: 0 ener added: 6.999733101e+46 Castro::do_new_sources() time = 0.058582105 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. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 3.75538647e+29 ener added: 3.75538647e+29 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -7.038531048e+46 ener added: 2.158325616e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.473668262e+40 ymom added: 6.057473935e+26 zmom added: 0 eint added: 0 ener added: 8.943835283e+46 Castro::do_old_sources() time = 0.055363958 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.257489714 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.076022112e+46 ener added: 2.369765675e+43 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015257285 Gravity::solve_for_phi() time = 0.081764014 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -3.768256697e+44 ener added: 1.742935385e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -9.789276502e+34 ymom added: 1.761383669e+26 zmom added: 0 eint added: 0 ener added: -1.824543748e+44 Castro::do_new_sources() time = 0.059010326 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 3.75538647e+29 ener added: 3.75538647e+29 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -7.038531048e+46 ener added: 2.158325616e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.473668262e+40 ymom added: 6.057473935e+26 zmom added: 0 eint added: 0 ener added: 8.943835283e+46 Castro::do_old_sources() time = 0.054937663 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.261550834 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.399341724e+44 ener added: 3.055404896e+43 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014687257 Gravity::solve_for_phi() time = 0.080095728 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.981378524e+46 ener added: 1.65468052e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -4.966651828e+36 ymom added: 1.133367956e+23 zmom added: 0 eint added: 0 ener added: 7.003642946e+46 Castro::do_new_sources() time = 0.058837574 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. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -7.047344521e+46 ener added: 2.557244057e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.474661592e+40 ymom added: 6.086941502e+26 zmom added: 0 eint added: 0 ener added: 8.976696152e+46 Castro::do_old_sources() time = 0.055001398 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.257059164 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.085134477e+46 ener added: 2.616682453e+43 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.0153746 Gravity::solve_for_phi() time = 0.081021071 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -3.966387031e+44 ener added: 1.764409733e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -7.086349894e+34 ymom added: 1.71437251e+26 zmom added: 0 eint added: 0 ener added: -1.814844154e+44 Castro::do_new_sources() time = 0.058171652 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -7.047344521e+46 ener added: 2.557244057e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.474661592e+40 ymom added: 6.086941502e+26 zmom added: 0 eint added: 0 ener added: 8.976696152e+46 Castro::do_old_sources() time = 0.054950941 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.261610597 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.20798565e+44 ener added: 3.358777738e+43 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014882963 Gravity::solve_for_phi() time = 0.080561105 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.988653412e+46 ener added: 1.671056473e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -4.948364053e+36 ymom added: 7.142579305e+23 zmom added: 0 eint added: 0 ener added: 7.013022432e+46 Castro::do_new_sources() time = 0.059483104 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002909546984 [Level 0 step 1] Advanced 32768 cells Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.548275694e+29 ener added: 1.548275694e+29 TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -2.890848341e+38 TIME= 0.25 YMOM = 8.975065285e+27 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.505143013e+38 TIME= 0.25 RHO*e = 1.331449006e+50 TIME= 0.25 RHO*K = 3.171006988e+49 TIME= 0.25 RHO*E = 1.648549705e+50 TIME= 0.25 RHO*PHI = -5.460195609e+50 TIME= 0.25 TOTAL ENERGY= -1.0815481e+50 TIME= 0.25 CENTER OF MASS X-LOC = 9.947183943e-17 TIME= 0.25 CENTER OF MASS X-VEL = -113366.7542 TIME= 0.25 CENTER OF MASS Y-LOC = -4.196397805e-29 TIME= 0.25 CENTER OF MASS Y-VEL = 3.519638182e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 [STEP 1] Coarse TimeStep time: 11.3570507 [STEP 1] FAB kilobyte spread across MPI nodes: [7064 ... 7064] STEP = 1 TIME = 0.25 DT = 0.25 PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.046443935 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.006530595 seconds Ending run at 23:51:02 UTC on 2022-04-30. Run time = 13.80810385 Run time without initialization = 11.41206775 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.81 ... 13.81 ... 13.81 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 4.9858 5.0692 5.1962 37.63% Amr::InitAmr() 1 1.5 1.5 1.501 10.87% Castro::computeTemp() 80 0.8496 0.9481 1.094 7.92% FabArray::ParallelCopy_finish() 1583 0.1346 0.4769 0.6695 4.85% MultiFab::AddProduct() 3439 0.5507 0.5549 0.5607 4.06% Castro::normalize_species() 79 0.5306 0.5346 0.5427 3.93% Castro::reset_internal_energy(MultiFab) 80 0.3623 0.3993 0.5296 3.84% Castro::reset_internal_energy(Fab) 160 0.505 0.5109 0.5167 3.74% MultiFab::Subtract() 159 0.3398 0.3437 0.3671 2.66% Castro::initMFs() 1 0.3527 0.3568 0.3654 2.65% MLPoisson::Fsmooth() 6144 0.2672 0.2896 0.3605 2.61% Castro::do_advance_ctu() 20 0.1369 0.2725 0.3604 2.61% FabArray::setVal() 5021 0.3501 0.3547 0.3603 2.61% Gravity::fill_multipole_BCs() 19 0.1999 0.2439 0.2648 1.92% FabArray::ParallelCopy_nowait() 1583 0.226 0.2432 0.2633 1.91% Castro::enforce_min_density() 79 0.1855 0.1941 0.2232 1.62% MultiFab::sum() 3439 0.1811 0.1824 0.184 1.33% Castro::initData() 1 0.003319 0.02277 0.1311 0.95% FillBoundary_nowait() 6388 0.07362 0.09911 0.1281 0.93% FillBoundary_finish() 6388 0.08341 0.09839 0.1128 0.82% MLPoisson::Fapply() 994 0.0935 0.09806 0.1078 0.78% Castro::derive() 156 0.09383 0.09531 0.09783 0.71% MLCellLinOp::applyBC() 7156 0.04529 0.06051 0.09324 0.68% Castro::do_old_sources() 20 0.0002787 0.01595 0.07253 0.53% Castro::construct_new_gravity_source() 17 0.06879 0.0695 0.07084 0.51% Castro::expand_state() 20 0.06661 0.06737 0.06829 0.49% MultiFab::Saxpy() 91 0.06597 0.06634 0.0669 0.48% MultiFab::Xpay() 923 0.05026 0.05298 0.06276 0.45% Castro::do_new_sources() 17 0.0002325 0.01132 0.05506 0.40% MultiFab::contains_nan() 37 0.04648 0.04674 0.04717 0.34% Castro::construct_old_gravity_source() 20 0.04442 0.04476 0.04531 0.33% MLMG::addInterpCorrection() 768 0.02915 0.03147 0.03834 0.28% Amr::writePlotFile() 2 0.02335 0.03184 0.03323 0.24% amrex::average_down 768 0.02518 0.02718 0.03308 0.24% StateData::FillBoundary(geom) 143 0 0.01707 0.03265 0.24% VisMF::Write(FabArray) 4 0.0129 0.01422 0.02578 0.19% MLCGSolver::bicgstab 8 0 0.001535 0.02457 0.18% Gravity::get_old_grav_vector() 20 0.01875 0.01904 0.02021 0.15% Castro::retry_advance_ctu() 11 0.0193 0.01974 0.02003 0.15% Castro::create_source_corrector() 17 0.01825 0.01842 0.01863 0.13% MLMG::ResNormInf() 146 0.01129 0.01512 0.01819 0.13% Gravity::get_new_grav_vector() 18 0.01696 0.0173 0.01782 0.13% MLCellLinOp::defineAuxData() 19 0.01488 0.01581 0.0175 0.13% FillPatchIterator::Initialize 301 0.0163 0.01659 0.01696 0.12% FabArray::FillBoundary() 6388 0.01302 0.01355 0.01572 0.11% Castro::construct_new_source() 102 0.01543 0.01557 0.0157 0.11% MLMG::oneIter() 128 0.01414 0.01422 0.01431 0.10% Castro::evaluate_source_change() 196 0.01254 0.01291 0.01356 0.10% MLCellLinOp::smooth() 3072 0.01103 0.01154 0.01296 0.09% FabArrayBase::getFB() 6388 0.009732 0.01026 0.01268 0.09% FabArray::setDomainBndry() 301 0.0008313 0.006366 0.01163 0.08% Amr::writeSmallPlotFile() 2 0.007066 0.01061 0.01113 0.08% Castro::volWgtSum() 66 0.01055 0.01068 0.0109 0.08% Castro::construct_old_source() 120 0.009813 0.009944 0.01076 0.08% Castro::wd_update() 2 0.004218 0.007811 0.01033 0.07% MultiFab::LinComb() 203 0 0.0006108 0.009773 0.07% MLLinOp::makeSubCommunicator() 19 0.0003385 0.004169 0.008032 0.06% Gravity::actual_solve_with_mlmg() 19 0.006741 0.007061 0.007769 0.06% FabArrayBase::getCPC() 2005 0.006972 0.007203 0.007585 0.05% FabArray::mult() 59 0.007141 0.007271 0.007547 0.05% MLMG::prepareForSolve() 18 0.006636 0.006884 0.00748 0.05% MLCGSolver::ParallelAllReduce 189 0 0.0004495 0.007192 0.05% MultiFab::Dot() 142 0 0.0003839 0.006142 0.04% FabArrayBase::CPC::define() 576 0.002171 0.003381 0.00614 0.04% CGSolver::sxay() 203 0 0.0003695 0.005912 0.04% MLCellLinOp::compGrad() 18 0.005564 0.00562 0.0057 0.04% MLMG::solve() 18 0.0009571 0.001819 0.005437 0.04% MLCellLinOp::apply() 994 0.002037 0.002306 0.005241 0.04% Castro::initialize_do_advance() 20 0.004248 0.004449 0.004663 0.03% MultiFab::max() 8 0.0005079 0.002542 0.004514 0.03% FabArray::ParallelCopy() 1583 0.003754 0.00393 0.004234 0.03% Castro::volInBoundary() 14 0.003574 0.003663 0.004013 0.03% Gravity::solve_for_phi() 17 0.002415 0.002654 0.003558 0.03% BndryData::define() 19 0.00313 0.003185 0.003311 0.02% MLPoisson::normalize() 79 0 0.0001969 0.003151 0.02% Castro::construct_new_gravity() 17 0.002221 0.002257 0.002388 0.02% MLCellLinOp::correctionResidual() 776 0.001846 0.001923 0.002249 0.02% MLLinOp::defineGrids() 19 0.001835 0.001989 0.002245 0.02% MLCellLinOp::prepareForSolve() 19 0.001707 0.001846 0.002088 0.02% MultiFab::Copy() 18 0.001989 0.002006 0.002024 0.01% Castro::locWgtSum() 12 0.001927 0.001952 0.001981 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.001276 0.001534 0.00185 0.01% MLMG::mgVcycle() 128 0.001562 0.001642 0.001784 0.01% MLCellLinOp::defineBC() 19 0.00163 0.001682 0.001773 0.01% MLMG:computeResOfCorrection() 768 0.001425 0.001507 0.001618 0.01% Castro::problem_post_timestep() 2 0.0003445 0.001145 0.001602 0.01% FabArrayBase::FB::FB() 75 0.000277 0.0006947 0.001321 0.01% StateDataPhysBCFunct::() 301 0.001101 0.001222 0.001316 0.01% MLMG::actualBottomSolve() 128 0.0001466 0.0002212 0.001274 0.01% Castro::post_init() 1 0.0003957 0.0009321 0.001271 0.01% StateData::define() 6 0.00115 0.001205 0.001258 0.01% main() 1 0.0009789 0.001136 0.001222 0.01% MLMG::MLRhsNormInf() 18 0.001089 0.001097 0.001106 0.01% FillPatchSingleLevel 301 0.0008409 0.0008554 0.0008766 0.01% Castro::volProductSum() 4 0.0008223 0.0008345 0.0008485 0.01% Castro::subcycle_advance_ctu() 1 0.0003231 0.0004001 0.0008457 0.01% MLCellLinOp::setLevelBC() 19 0.0005435 0.0006764 0.0008374 0.01% MLMG::mgVcycle_down::2 128 0.000659 0.0007175 0.0007979 0.01% MLMG::mgVcycle_down::3 128 0.0006007 0.0006477 0.000797 0.01% MLMG::mgVcycle_down::0 128 0.0007165 0.0007366 0.0007683 0.01% MLMG::mgVcycle_down::1 128 0.0007033 0.0007219 0.0007581 0.01% TagBoxArray::collate() 2 9.193e-05 0.0003748 0.0007318 0.01% MLMG::mgVcycle_down::4 128 0.0005873 0.0006181 0.0007128 0.01% MLMG::mgVcycle_down::5 128 0.0005843 0.0006082 0.0006741 0.00% Castro::swap_state_time_levels() 8 0.0004202 0.000626 0.0006711 0.00% Castro::Castro() 1 0.0004604 0.0005316 0.0006091 0.00% Castro::clean_state() 79 0.000517 0.000538 0.0005783 0.00% MLMG::mgVcycle_up::1 128 0.000505 0.0005332 0.0005744 0.00% Gravity::test_composite_phi() 1 0.0004005 0.0004684 0.0005671 0.00% Castro::buildMetrics() 1 0.0005307 0.0005512 0.0005666 0.00% MLMG::mgVcycle_up::2 128 0.0004195 0.0004675 0.0005485 0.00% MLMG::mgVcycle_up::0 128 0.0005176 0.0005287 0.0005482 0.00% MLMG::mgVcycle_up::3 128 0.0004228 0.000444 0.000538 0.00% MLMG::mgVcycle_up::5 128 0.0004433 0.0004596 0.0005312 0.00% makeSFC 95 0.0003974 0.0004322 0.0005213 0.00% AmrMesh::MakeNewGrids() 2 0.0004856 0.0004976 0.0005096 0.00% Gravity::applyMetricTerms() 2 0.0004749 0.0004865 0.0004951 0.00% MLMG::mgVcycle_up::4 128 0.0004267 0.0004377 0.0004922 0.00% Castro::apply_tagging_restrictions() 2 0.0003608 0.0004265 0.0004886 0.00% Castro::post_timestep() 1 0.0002157 0.000411 0.0004669 0.00% MLCellLinOp::solutionResidual() 147 0.0004242 0.0004314 0.0004421 0.00% MLLinOp::makeAgglomeratedDMap 19 0.0003651 0.0003829 0.0004136 0.00% Gravity::update_max_rhs() 2 0.0003813 0.0003915 0.0003989 0.00% Castro::enforce_consistent_e() 1 0.0003651 0.0003773 0.0003914 0.00% MLMG::mgVcycle_bottom 128 0.0002224 0.000232 0.0002931 0.00% MLMG::computeResidual() 128 0.0002676 0.0002732 0.0002825 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002403 0.0002543 0.0002811 0.00% MLMG::compResidual() 1 0.000248 0.0002538 0.0002577 0.00% DistributionMapping::Distribute() 95 0.0002218 0.0002334 0.0002488 0.00% Castro::finalize_do_advance() 17 0.0001386 0.0001491 0.0001726 0.00% MLLinOp::define() 19 0.0001499 0.0001583 0.0001721 0.00% OwnerMask() 2 0.0001331 0.0001518 0.0001687 0.00% Amr::defBaseLevel() 1 0.0001104 0.0001406 0.0001591 0.00% MLPoisson::define() 19 0.0001124 0.0001228 0.0001473 0.00% Gravity::actual_multilevel_solve() 1 0.0001325 0.000137 0.0001406 0.00% Castro::check_for_nan() 37 0.0001127 0.0001193 0.0001392 0.00% Amr::coarseTimeStep() 1 0.000107 0.0001323 0.0001356 0.00% Castro::enforce_speed_limit() 79 0.0001082 0.0001152 0.0001343 0.00% TagBoxArray::mapPRD 2 0.0001152 0.0001216 0.0001333 0.00% Castro::construct_old_gravity() 20 9.227e-05 0.0001012 0.0001152 0.00% Castro::initialize_advance() 1 9.493e-05 0.0001033 0.0001128 0.00% Castro::apply_source_to_state() 37 9.563e-05 0.0001024 0.0001079 0.00% Amr::grid_places() 2 8.168e-05 8.967e-05 9.983e-05 0.00% Castro::apply_problem_tags() 2 9.029e-05 9.226e-05 9.466e-05 0.00% MultiFab::Swap() 5 7.055e-05 7.778e-05 9.162e-05 0.00% MLPoisson::prepareForSolve() 19 7.1e-05 7.61e-05 8.272e-05 0.00% Castro::finalize_advance() 1 1.92e-05 2.946e-05 6.822e-05 0.00% knapsack() 1 5.914e-05 6.223e-05 6.457e-05 0.00% Amr::initSubcycle() 1 5.704e-05 5.868e-05 5.957e-05 0.00% DistributionMapping::LeastUsedCPUs() 1 1.469e-05 3.052e-05 5.459e-05 0.00% MLMG::buildFineMask() 18 4.367e-05 4.665e-05 5.297e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 4.982e-05 5.088e-05 5.244e-05 0.00% MLMG::computeMLResidual() 18 4.219e-05 4.548e-05 5.217e-05 0.00% Amr::FinalizeInit() 1 8.644e-06 1.248e-05 5.147e-05 0.00% MLMG::getGradSolution() 18 4.118e-05 4.248e-05 4.388e-05 0.00% MLMG::MLResNormInf() 18 3.491e-05 3.692e-05 4.133e-05 0.00% Gravity::swapTimeLevels() 8 3.192e-05 3.54e-05 3.778e-05 0.00% Castro::estTimeStep() 19 2.683e-05 2.928e-05 3.535e-05 0.00% Castro::construct_old_rotation_source() 20 3.021e-05 3.118e-05 3.429e-05 0.00% Amr::regrid() 1 1.189e-05 1.698e-05 3.415e-05 0.00% AmrLevel::AmrLevel(dm) 1 3.059e-05 3.145e-05 3.261e-05 0.00% Castro::construct_new_rotation_source() 17 2.312e-05 2.418e-05 2.834e-05 0.00% Amr::timeStep() 1 1.216e-05 1.459e-05 2.128e-05 0.00% DistributionMapping::KnapSackDoIt() 1 1.727e-05 1.815e-05 2.027e-05 0.00% Amr::init() 1 7.396e-06 8.724e-06 2.019e-05 0.00% Amr::InitializeInit() 1 6.672e-06 8.302e-06 1.581e-05 0.00% Amr::initialInit() 1 2.834e-06 3.612e-06 1.259e-05 0.00% Castro::post_regrid() 2 1.094e-05 1.148e-05 1.237e-05 0.00% Amr::bldFineLevels() 1 7.169e-06 7.898e-06 8.979e-06 0.00% BoxList::parallelComplementIn() 4 7.214e-06 7.42e-06 7.814e-06 0.00% Castro::advance() 1 6.279e-06 6.969e-06 7.599e-06 0.00% Castro::computeInitialDt() 2 6.546e-06 7.053e-06 7.365e-06 0.00% Castro::errorEst() 2 6.05e-06 6.372e-06 6.703e-06 0.00% BoxList::complementIn 4 6.187e-06 6.306e-06 6.474e-06 0.00% Castro::gwstrain() 2 4.439e-06 4.848e-06 6.368e-06 0.00% Gravity::set_mass_offset() 2 4.319e-06 4.672e-06 5.159e-06 0.00% Castro::FluxRegCrseInit 1 1.78e-06 2.089e-06 2.641e-06 0.00% Castro::FluxRegFineAdd() 1 1.41e-06 1.588e-06 1.955e-06 0.00% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 13.81 13.81 13.81 100.00% Amr::coarseTimeStep() 1 11.41 11.41 11.41 82.64% Amr::timeStep() 1 11.36 11.36 11.36 82.25% Castro::advance() 1 11.26 11.26 11.26 81.56% Castro::subcycle_advance_ctu() 1 11.21 11.21 11.21 81.21% Castro::do_advance_ctu() 20 11.17 11.17 11.18 80.94% Castro::construct_ctu_hydro_source() 20 4.986 5.069 5.196 37.63% Castro::clean_state() 79 3.673 3.757 3.874 28.06% Castro::computeTemp() 80 2.367 2.474 2.654 19.22% Castro::reset_internal_energy(MultiFab) 80 1.479 1.526 1.646 11.92% Amr::InitAmr() 1 1.501 1.501 1.501 10.87% Gravity::solve_phi_with_mlmg() 19 1.404 1.451 1.474 10.67% Castro::construct_new_gravity() 17 1.333 1.381 1.404 10.16% Gravity::solve_for_phi() 17 1.298 1.345 1.368 9.91% Castro::initialize_do_advance() 20 1.262 1.274 1.282 9.29% Gravity::actual_solve_with_mlmg() 19 1.202 1.205 1.21 8.77% MLMG::solve() 18 1.152 1.155 1.157 8.38% Castro::do_old_sources() 20 1.03 1.082 1.138 8.24% MLMG::oneIter() 128 1.027 1.029 1.032 7.47% MLMG::mgVcycle() 128 1.013 1.015 1.018 7.37% Castro::do_new_sources() 17 0.9408 0.9735 1.004 7.27% Castro::evaluate_source_change() 196 0.9451 0.9512 0.9603 6.95% FabArray::ParallelCopy() 1583 0.3892 0.7303 0.924 6.69% Amr::init() 1 0.8948 0.8948 0.8948 6.48% Amr::initialInit() 1 0.8335 0.8335 0.8336 6.04% Castro::enforce_min_density() 79 0.7886 0.8008 0.8266 5.99% FabArray::ParallelCopy_finish() 1583 0.1346 0.4769 0.6695 4.85% MLCellLinOp::smooth() 3072 0.4834 0.5279 0.6464 4.68% Amr::InitializeInit() 1 0.4367 0.4878 0.5768 4.18% Amr::defBaseLevel() 1 0.4367 0.4878 0.5768 4.18% MultiFab::AddProduct() 3439 0.5507 0.5549 0.5607 4.06% Castro::normalize_species() 79 0.5306 0.5346 0.5427 3.93% Castro::reset_internal_energy(Fab) 160 0.505 0.5109 0.5167 3.74% FillPatchIterator::Initialize 301 0.3209 0.4007 0.4831 3.50% FillPatchSingleLevel 301 0.2928 0.3777 0.4614 3.34% MLMG::mgVcycle_up::1 128 0.08447 0.2332 0.4098 2.97% Amr::FinalizeInit() 1 0.2567 0.3458 0.3969 2.87% Castro::Castro() 1 0.3581 0.3624 0.3722 2.70% MultiFab::Subtract() 159 0.3398 0.3437 0.3671 2.66% Castro::initMFs() 1 0.3529 0.3569 0.3656 2.65% MLMG::addInterpCorrection() 768 0.1402 0.293 0.3656 2.65% MLPoisson::Fsmooth() 6144 0.2672 0.2896 0.3605 2.61% FabArray::setVal() 5021 0.3501 0.3547 0.3603 2.61% MLCellLinOp::applyBC() 7156 0.2539 0.2825 0.328 2.38% MLMG::mgVcycle_up::2 128 0.009719 0.1158 0.279 2.02% FabArray::ParallelCopy_nowait() 1583 0.2316 0.2494 0.2707 1.96% Gravity::fill_multipole_BCs() 19 0.1999 0.2439 0.2648 1.92% FabArray::FillBoundary() 6388 0.1968 0.222 0.2443 1.77% MLMG::mgVcycle_down::0 128 0.2365 0.239 0.2427 1.76% Castro::initData() 1 0.07679 0.1239 0.212 1.54% Castro::expand_state() 20 0.1786 0.1855 0.2007 1.45% MLMG::mgVcycle_up::3 128 0.009694 0.02292 0.1864 1.35% Amr::bldFineLevels() 1 0.04633 0.1349 0.1852 1.34% MultiFab::sum() 3439 0.1811 0.1824 0.184 1.33% MLMG::mgVcycle_up::0 128 0.1708 0.172 0.1733 1.26% MLCellLinOp::apply() 994 0.1445 0.1545 0.173 1.25% Castro::post_init() 1 0.1657 0.1677 0.1697 1.23% amrex::average_down 768 0.03334 0.04762 0.1675 1.21% MLMG::mgVcycle_down::5 128 0.01168 0.02093 0.1454 1.05% Amr::grid_places() 2 0.002601 0.09167 0.1428 1.03% AmrMesh::MakeNewGrids() 2 0.002503 0.09158 0.1427 1.03% Castro::derive() 156 0.1375 0.1393 0.1422 1.03% TagBoxArray::mapPRD 2 0.0007107 0.09003 0.1413 1.02% MLMG:computeResOfCorrection() 768 0.1102 0.1192 0.1408 1.02% FillBoundary_nowait() 6388 0.08381 0.1101 0.1402 1.02% MLCellLinOp::correctionResidual() 776 0.1087 0.118 0.1392 1.01% Castro::construct_new_source() 102 0.1222 0.1232 0.1249 0.90% FillBoundary_finish() 6388 0.08341 0.09839 0.1128 0.82% MLMG::mgVcycle_down::1 128 0.1057 0.108 0.1125 0.81% MLPoisson::Fapply() 994 0.0935 0.09806 0.1078 0.78% MLMG::mgVcycle_up::5 128 0.009774 0.01612 0.107 0.78% Amr::writePlotFile() 2 0.101 0.101 0.1011 0.73% Castro::post_timestep() 1 0.09477 0.09661 0.09834 0.71% Gravity::multilevel_solve_for_new_phi() 1 0.0949 0.09493 0.09497 0.69% Gravity::actual_multilevel_solve() 1 0.09485 0.09488 0.09492 0.69% MLCellLinOp::solutionResidual() 147 0.08438 0.09039 0.09361 0.68% Castro::create_source_corrector() 17 0.04825 0.06773 0.08839 0.64% MLMG::mgVcycle_bottom 128 0.0003717 0.005875 0.08832 0.64% Castro::post_regrid() 2 0.08437 0.08633 0.08831 0.64% MLMG::actualBottomSolve() 128 0.0001466 0.005643 0.08802 0.64% MLCGSolver::bicgstab 8 0 0.005388 0.08621 0.62% Castro::construct_old_source() 120 0.07675 0.07733 0.07896 0.57% MLMG::computeResidual() 128 0.07348 0.07488 0.07854 0.57% Castro::problem_post_timestep() 2 0.07304 0.07482 0.07672 0.56% Castro::construct_new_gravity_source() 17 0.06879 0.0695 0.07084 0.51% MultiFab::Saxpy() 91 0.06597 0.06634 0.0669 0.48% MLMG::mgVcycle_down::2 128 0.01242 0.03783 0.06599 0.48% MultiFab::Xpay() 923 0.05026 0.05298 0.06276 0.45% Castro::construct_old_gravity() 20 0.04196 0.05296 0.06162 0.45% Gravity::get_old_grav_vector() 20 0.04186 0.05286 0.06152 0.45% MLMG::mgVcycle_down::3 128 0.01178 0.0177 0.05748 0.42% Castro::volWgtSum() 66 0.05278 0.0534 0.05426 0.39% Castro::initialize_advance() 1 0.04854 0.0496 0.05121 0.37% Castro::check_for_nan() 37 0.0466 0.04686 0.04729 0.34% MultiFab::contains_nan() 37 0.04648 0.04674 0.04717 0.34% Castro::construct_old_gravity_source() 20 0.04442 0.04476 0.04531 0.33% Amr::regrid() 1 0.04314 0.04419 0.04524 0.33% MLPoisson::define() 19 0.02946 0.03408 0.04098 0.30% Castro::retry_advance_ctu() 11 0.03633 0.03708 0.03767 0.27% Gravity::get_new_grav_vector() 18 0.03459 0.03532 0.03618 0.26% Castro::wd_update() 2 0.03043 0.03231 0.03423 0.25% StateDataPhysBCFunct::() 301 0.001101 0.01829 0.03395 0.25% MLMG::mgVcycle_down::4 128 0.01172 0.01331 0.03352 0.24% StateData::FillBoundary(geom) 143 0 0.01707 0.03265 0.24% Castro::apply_source_to_state() 37 0.02697 0.02713 0.02766 0.20% VisMF::Write(FabArray) 4 0.0129 0.01422 0.02578 0.19% MLMG::mgVcycle_up::4 128 0.009724 0.01085 0.02504 0.18% MLCellLinOp::defineAuxData() 19 0.01892 0.02068 0.02391 0.17% Castro::volInBoundary() 14 0.01989 0.02042 0.02229 0.16% MLMG::prepareForSolve() 18 0.01694 0.01784 0.01991 0.14% MLMG::ResNormInf() 146 0.01129 0.01512 0.01819 0.13% MLMG::computeMLResidual() 18 0.01035 0.01513 0.0179 0.13% CGSolver::sxay() 203 0 0.0009803 0.01568 0.11% Amr::writeSmallPlotFile() 2 0.01465 0.01465 0.01465 0.11% Gravity::test_composite_phi() 1 0.014 0.01404 0.01407 0.10% FabArrayBase::getCPC() 2005 0.009146 0.01058 0.01361 0.10% FabArrayBase::getFB() 6388 0.01002 0.01095 0.01321 0.10% FabArray::setDomainBndry() 301 0.0008313 0.006366 0.01163 0.08% MLLinOp::define() 19 0.003419 0.007365 0.01162 0.08% MLLinOp::defineGrids() 19 0.003268 0.007207 0.01145 0.08% MLCGSolver::ParallelAllReduce 189 0 0.0006601 0.01056 0.08% MultiFab::LinComb() 203 0 0.0006108 0.009773 0.07% MLLinOp::makeSubCommunicator() 19 0.0003385 0.004169 0.008032 0.06% FabArray::mult() 59 0.007141 0.007271 0.007547 0.05% MLMG::getGradSolution() 18 0.007137 0.007332 0.007521 0.05% MLCellLinOp::compGrad() 18 0.007096 0.007289 0.007477 0.05% MLCellLinOp::defineBC() 19 0.005756 0.005916 0.006319 0.05% MultiFab::Dot() 142 0 0.0003839 0.006142 0.04% FabArrayBase::CPC::define() 576 0.002171 0.003381 0.00614 0.04% Gravity::update_max_rhs() 2 0.00141 0.003427 0.005388 0.04% Castro::locWgtSum() 12 0.00514 0.0052 0.005276 0.04% BndryData::define() 19 0.004037 0.004234 0.004558 0.03% MultiFab::max() 8 0.0005079 0.002542 0.004514 0.03% MLPoisson::normalize() 79 0 0.0001969 0.003151 0.02% Castro::volProductSum() 4 0.002959 0.002995 0.003044 0.02% MLPoisson::prepareForSolve() 19 0.001788 0.001922 0.002165 0.02% MLCellLinOp::prepareForSolve() 19 0.001707 0.001846 0.002088 0.02% MultiFab::Copy() 18 0.001989 0.002006 0.002024 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.001276 0.001534 0.00185 0.01% FabArrayBase::FB::FB() 75 0.000277 0.0006947 0.001321 0.01% AmrLevel::AmrLevel(dm) 1 0.001181 0.001236 0.00129 0.01% StateData::define() 6 0.00115 0.001205 0.001258 0.01% Gravity::swapTimeLevels() 8 0.001166 0.001204 0.001253 0.01% MLMG::compResidual() 1 0.0009211 0.001051 0.001189 0.01% MLMG::MLResNormInf() 18 0.001141 0.001152 0.00117 0.01% MLLinOp::makeAgglomeratedDMap 19 0.0009878 0.001049 0.001169 0.01% MLMG::MLRhsNormInf() 18 0.001089 0.001097 0.001106 0.01% MLCellLinOp::setLevelBC() 19 0.0005435 0.0006764 0.0008374 0.01% makeSFC 95 0.0006228 0.0006656 0.0007554 0.01% TagBoxArray::collate() 2 9.193e-05 0.0003748 0.0007318 0.01% Castro::swap_state_time_levels() 8 0.0004202 0.000626 0.0006711 0.00% Castro::buildMetrics() 1 0.0005908 0.0006127 0.0006295 0.00% Castro::errorEst() 2 0.0004583 0.0005251 0.0005881 0.00% Gravity::applyMetricTerms() 2 0.0004749 0.0004865 0.0004951 0.00% Castro::apply_tagging_restrictions() 2 0.0003608 0.0004265 0.0004886 0.00% Castro::enforce_consistent_e() 1 0.0003651 0.0003773 0.0003914 0.00% DistributionMapping::Distribute() 95 0.0002218 0.0002334 0.0002488 0.00% Castro::finalize_do_advance() 17 0.0001386 0.0001491 0.0001726 0.00% OwnerMask() 2 0.0001331 0.0001518 0.0001687 0.00% DistributionMapping::KnapSackDoIt() 1 9.347e-05 0.0001109 0.0001393 0.00% Castro::enforce_speed_limit() 79 0.0001082 0.0001152 0.0001343 0.00% Castro::apply_problem_tags() 2 9.029e-05 9.226e-05 9.466e-05 0.00% MultiFab::Swap() 5 7.055e-05 7.778e-05 9.162e-05 0.00% Castro::finalize_advance() 1 2.258e-05 3.314e-05 7.204e-05 0.00% knapsack() 1 5.914e-05 6.223e-05 6.457e-05 0.00% Amr::initSubcycle() 1 5.704e-05 5.868e-05 5.957e-05 0.00% DistributionMapping::LeastUsedCPUs() 1 1.469e-05 3.052e-05 5.459e-05 0.00% MLMG::buildFineMask() 18 4.367e-05 4.665e-05 5.297e-05 0.00% Castro::estTimeStep() 19 2.683e-05 2.928e-05 3.535e-05 0.00% Castro::construct_old_rotation_source() 20 3.021e-05 3.118e-05 3.429e-05 0.00% Castro::construct_new_rotation_source() 17 2.312e-05 2.418e-05 2.834e-05 0.00% BoxList::parallelComplementIn() 4 1.344e-05 1.373e-05 1.427e-05 0.00% Castro::computeInitialDt() 2 9.486e-06 1.008e-05 1.052e-05 0.00% BoxList::complementIn 4 6.187e-06 6.306e-06 6.474e-06 0.00% Castro::gwstrain() 2 4.439e-06 4.848e-06 6.368e-06 0.00% Gravity::set_mass_offset() 2 4.319e-06 4.672e-06 5.159e-06 0.00% Castro::FluxRegCrseInit 1 1.78e-06 2.089e-06 2.641e-06 0.00% Castro::FluxRegFineAdd() 1 1.41e-06 1.588e-06 1.955e-06 0.00% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::castro.center(nvals = 2) :: [0.0e0, 0.0e0] [TOP]::castro.dtnuc_mode(nvals = 1) :: [1] [TOP]::mg.v(nvals = 1) :: [0] AMReX (22.04-39-g50466d11e7b2) finalized