MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (22.01-39-g874ccc6f3f8d) initialized Starting run at 07:56:28 UTC on 2022-01-26. Successfully read inputs file ... Castro git describe: 21.12-27-g958e0cc7a AMReX git describe: 22.01-39-g874ccc6f3 Microphysics git describe: 22.01-20-g1a290c72 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.use_stopping_criterion = 1 castro.dtnuc_mode = 1 Warning: the following integrator.* parameters are ignored integrator.do_constant_volume_burn = 1 Initializing the data at level 0 Done initializing the level 0 data grid_places() time: 0.149175879 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.000957466 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.01036963 Gravity::fill_multipole_BCs() time = 0.010170662 ... 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.046628241 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.006448705 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.069154059 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.27072506 Advance was unsuccessful with reason: negative density (density = -9.081732e+02 at index (14,70); started at 4.477055e+03); 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.068403324 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.262386557 Advance was unsuccessful with reason: negative density (density = -1.592544e+00 at index (20,79); started at 2.790161e+00); 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.064226374 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.31528064 Advance was unsuccessful with reason: negative density (density = -1.060745e-01 at index (20,78); started at 2.068339e-01); 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.063477368 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.275262648 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015915659 Gravity::solve_for_phi() time = 0.082445845 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.070261472 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.064000998 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.26800437 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015025343 Gravity::solve_for_phi() time = 0.081657034 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.070548886 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.064142657 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.273379681 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014828271 Gravity::solve_for_phi() time = 0.081301669 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.07068938 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.063796772 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.275590233 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.015597435 Gravity::solve_for_phi() time = 0.081752125 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.069642494 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.064325368 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.279581345 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.014938603 Gravity::solve_for_phi() time = 0.081117698 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.070691103 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.063669736 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.284075084 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.016757592 Gravity::solve_for_phi() time = 0.085060697 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.069919969 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.063353813 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.275016429 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.01443902 Gravity::solve_for_phi() time = 0.079997156 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.06914635 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.0636465 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.267665874 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.015158006 Gravity::solve_for_phi() time = 0.080274807 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.068840597 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.063397506 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.273811543 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.014119732 Gravity::solve_for_phi() time = 0.079419175 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.06908163 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.063504246 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.268536592 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.014827271 Gravity::solve_for_phi() time = 0.080155354 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.070442275 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.063277182 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.274843263 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.014178744 Gravity::solve_for_phi() time = 0.080012059 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.069443961 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.063212585 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.26885636 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.014823587 Gravity::solve_for_phi() time = 0.080343525 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.068426186 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.06342678 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.273572885 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.014403176 Gravity::solve_for_phi() time = 0.079835634 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.069010675 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.064112807 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.273513785 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.015719492 Gravity::solve_for_phi() time = 0.08120519 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.068892941 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.063577063 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.28069623 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.01561569 Gravity::solve_for_phi() time = 0.08189045 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.0712661 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.063811694 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.276166581 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.015972967 Gravity::solve_for_phi() time = 0.082274044 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.069601031 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.064596017 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.280665869 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.01591849 Gravity::solve_for_phi() time = 0.082259704 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.070976623 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002724064375 [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: 12.12507793 [STEP 1] FAB kilobyte spread across MPI nodes: [8620 ... 8620] STEP = 1 TIME = 0.25 DT = 0.25 PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.047192763 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.006471518 seconds Ending run at 07:56:43 UTC on 2022-01-26. Run time = 14.45643759 Run time without initialization = 12.18100958 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]: 14.46 ... 14.46 ... 14.46 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.2856 5.4027 5.5077 38.10% Amr::InitAmr() 1 1.385 1.385 1.386 9.58% Castro::computeTemp() 80 0.8514 0.9646 1.114 7.70% Castro::reset_internal_energy() 80 0.8733 0.9202 1.047 7.24% MultiFab::AddProduct() 3920 0.6069 0.6174 0.6749 4.67% FabArray::ParallelCopy_finish() 1583 0.1495 0.4577 0.6368 4.40% FabArray::setVal() 5502 0.4534 0.4598 0.4657 3.22% Castro::normalize_species() 79 0.4445 0.449 0.4529 3.13% Amr::defBaseLevel() 1 0.3539 0.3606 0.3735 2.58% MLPoisson::Fsmooth() 6144 0.2653 0.2902 0.35 2.42% MultiFab::Subtract() 159 0.3351 0.3383 0.3415 2.36% FabArray::ParallelCopy_nowait() 1583 0.2969 0.3189 0.3406 2.36% Gravity::fill_multipole_BCs() 19 0.1967 0.2473 0.278 1.92% MultiFab::min() 23 0.002079 0.1347 0.2718 1.88% Castro::enforce_min_density() 79 0.1862 0.1968 0.2357 1.63% MultiFab::sum() 3920 0.2113 0.2137 0.2162 1.50% MultiFab::Saxpy() 91 0.201 0.203 0.2059 1.42% Castro::do_advance_ctu() 20 0.1518 0.1532 0.1546 1.07% FillBoundary_nowait() 6388 0.07899 0.1066 0.1357 0.94% Castro::initData() 1 0.003388 0.02293 0.1308 0.90% FillBoundary_finish() 6388 0.07308 0.09185 0.123 0.85% MLPoisson::Fapply() 994 0.09243 0.09729 0.1059 0.73% Castro::derive() 156 0.09509 0.09654 0.09955 0.69% MLCellLinOp::applyBC() 7156 0.04536 0.06069 0.09242 0.64% Castro::construct_new_gravity_source() 17 0.07934 0.08055 0.08209 0.57% Castro::expand_state() 20 0.06652 0.06741 0.06838 0.47% Castro::construct_old_gravity_source() 20 0.06337 0.06412 0.06502 0.45% Castro::do_old_sources() 20 0.000305 0.01446 0.06263 0.43% Castro::do_new_sources() 17 0.0002509 0.01198 0.05951 0.41% MultiFab::Xpay() 923 0.04967 0.05206 0.0569 0.39% Castro::create_source_corrector() 17 0.05136 0.05189 0.05338 0.37% MultiFab::contains_nan() 37 0.04454 0.04513 0.04558 0.32% StateData::FillBoundary(geom) 143 0 0.02209 0.04264 0.29% MLMG::addInterpCorrection() 768 0.02904 0.03157 0.038 0.26% amrex::average_down 768 0.02496 0.02716 0.03272 0.23% Amr::writePlotFile() 2 0.01742 0.02598 0.02754 0.19% VisMF::Write(FabArray) 4 0.01294 0.01439 0.02678 0.19% MLCGSolver::bicgstab 8 0 0.001534 0.02454 0.17% Castro::retry_advance_ctu() 11 0.02387 0.02395 0.02407 0.17% FillPatchIterator::Initialize 301 0.0203 0.02059 0.02097 0.15% TagBoxArray::collate() 2 6.667e-05 0.003151 0.02087 0.14% Gravity::get_old_grav_vector() 20 0.01872 0.01891 0.01914 0.13% Castro::construct_new_source() 102 0.01691 0.0172 0.0184 0.13% MLCellLinOp::defineAuxData() 19 0.01581 0.01648 0.018 0.12% Gravity::get_new_grav_vector() 18 0.01731 0.01749 0.01771 0.12% FabArray::setDomainBndry() 301 0.0008409 0.008875 0.0165 0.11% FabArray::FillBoundary() 6388 0.01325 0.0138 0.01606 0.11% Castro::evaluate_source_change() 196 0.01359 0.01414 0.01506 0.10% MLMG::ResNormInf() 146 0.01068 0.01275 0.01443 0.10% MLMG::oneIter() 128 0.01408 0.01425 0.01443 0.10% Castro::wd_update() 2 0.003892 0.009192 0.01289 0.09% FabArrayBase::getFB() 6388 0.009576 0.01015 0.01226 0.08% MLCellLinOp::smooth() 3072 0.01071 0.01131 0.01223 0.08% Castro::construct_old_source() 120 0.01137 0.01153 0.01204 0.08% Amr::writeSmallPlotFile() 2 0.006879 0.01048 0.01105 0.08% Castro::volWgtSum() 66 0.009872 0.009974 0.01009 0.07% MultiFab::LinComb() 203 0 0.000615 0.00984 0.07% MLLinOp::makeSubCommunicator() 19 0.0003164 0.00392 0.008092 0.06% Gravity::actual_solve_with_mlmg() 19 0.006811 0.007144 0.007713 0.05% FabArray::mult() 59 0.007228 0.00735 0.00757 0.05% FabArrayBase::getCPC() 2005 0.007025 0.007227 0.007538 0.05% MLMG::prepareForSolve() 18 0.006743 0.007031 0.007334 0.05% MLCGSolver::ParallelAllReduce 189 0 0.0004356 0.00697 0.05% MultiFab::Dot() 142 0 0.000379 0.006065 0.04% FabArrayBase::CPC::define() 576 0.002193 0.003376 0.006054 0.04% CGSolver::sxay() 203 0 0.0003694 0.00591 0.04% MLCellLinOp::compGrad() 18 0.005553 0.005687 0.005798 0.04% MLCellLinOp::apply() 994 0.001998 0.00229 0.005076 0.04% Castro::initialize_do_advance() 20 0.004344 0.004452 0.004623 0.03% FabArray::ParallelCopy() 1583 0.003731 0.003937 0.004216 0.03% Castro::volInBoundary() 14 0.003487 0.003569 0.003983 0.03% MLMG::solve() 18 0.0007631 0.001763 0.003791 0.03% BndryData::define() 19 0.00321 0.003309 0.003422 0.02% Gravity::solve_for_phi() 17 0.002431 0.002692 0.003403 0.02% MLPoisson::normalize() 79 0 0.0001919 0.00307 0.02% Castro::construct_new_gravity() 17 0.00222 0.002296 0.002425 0.02% MLCellLinOp::correctionResidual() 776 0.001829 0.001933 0.002279 0.02% MLCellLinOp::prepareForSolve() 19 0.001744 0.001903 0.002094 0.01% MultiFab::Copy() 18 0.001988 0.002022 0.002055 0.01% MLLinOp::defineGrids() 19 0.00181 0.001904 0.002023 0.01% Castro::locWgtSum() 12 0.001781 0.001797 0.001821 0.01% MLMG::mgVcycle() 128 0.001572 0.001675 0.001807 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.001301 0.001552 0.001795 0.01% MLCellLinOp::defineBC() 19 0.001662 0.001727 0.00179 0.01% MultiFab::max() 8 0.0005198 0.001156 0.00165 0.01% StateData::define() 6 0.001403 0.001435 0.001616 0.01% MLMG:computeResOfCorrection() 768 0.001445 0.001507 0.001602 0.01% Castro::problem_post_timestep() 2 0.0002769 0.0007737 0.001551 0.01% StateDataPhysBCFunct::() 301 0.001145 0.001288 0.001401 0.01% FabArrayBase::FB::FB() 75 0.000273 0.000681 0.001286 0.01% MLMG::actualBottomSolve() 128 0.0001463 0.0002207 0.001281 0.01% main() 1 0.0009401 0.00114 0.001253 0.01% MLMG::MLRhsNormInf() 18 0.001135 0.001155 0.001178 0.01% Castro::swap_state_time_levels() 8 0.0005148 0.0009672 0.001012 0.01% Castro::subcycle_advance_ctu() 1 0.0002829 0.0003838 0.0008909 0.01% FillPatchSingleLevel 301 0.0008403 0.0008652 0.0008863 0.01% MLMG::mgVcycle_down::2 128 0.0006541 0.0007255 0.0008235 0.01% MLCellLinOp::setLevelBC() 19 0.0005339 0.000689 0.0008172 0.01% MLMG::mgVcycle_down::0 128 0.0007136 0.0007405 0.0007942 0.01% Castro::volProductSum() 4 0.0007576 0.0007739 0.000788 0.01% MLMG::mgVcycle_down::1 128 0.0006989 0.0007258 0.0007758 0.01% MLMG::mgVcycle_down::3 128 0.0005972 0.0006546 0.0007526 0.01% Castro::post_timestep() 1 0.0002309 0.000466 0.0007499 0.01% MLMG::mgVcycle_down::4 128 0.0005961 0.0006299 0.0007462 0.01% MLMG::mgVcycle_down::5 128 0.0005885 0.0006152 0.0007085 0.00% Castro::post_init() 1 0.0004312 0.0006463 0.0006997 0.00% DistributionMapping::LeastUsedCPUs() 1 7.06e-05 0.0006043 0.0006493 0.00% Gravity::test_composite_phi() 1 0.0003841 0.0004383 0.0005677 0.00% MLMG::mgVcycle_up::1 128 0.0005027 0.0005288 0.0005639 0.00% MLMG::mgVcycle_up::0 128 0.000493 0.0005145 0.0005442 0.00% AmrMesh::MakeNewGrids() 2 0.0004872 0.0005248 0.0005432 0.00% MLMG::mgVcycle_up::2 128 0.000416 0.0004688 0.0005343 0.00% MLMG::mgVcycle_up::3 128 0.0004213 0.000444 0.0005158 0.00% Castro::clean_state() 79 0.000434 0.0004627 0.0005156 0.00% MLMG::mgVcycle_up::5 128 0.0004372 0.0004551 0.0004987 0.00% MLMG::mgVcycle_up::4 128 0.0004174 0.0004336 0.0004866 0.00% Gravity::applyMetricTerms() 2 0.0004626 0.0004684 0.0004742 0.00% makeSFC 95 0.0003967 0.0004271 0.0004589 0.00% MLCellLinOp::solutionResidual() 147 0.0004183 0.0004344 0.0004527 0.00% Gravity::update_max_rhs() 2 0.0003824 0.000392 0.0004061 0.00% MLLinOp::makeAgglomeratedDMap 19 0.0003734 0.0003857 0.0003967 0.00% Castro::enforce_consistent_e() 1 0.0003611 0.0003705 0.0003959 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002434 0.0002582 0.0002999 0.00% Castro::initialize_advance() 1 0.0002254 0.000272 0.0002856 0.00% MLMG::computeResidual() 128 0.0002625 0.0002722 0.0002798 0.00% MLMG::mgVcycle_bottom 128 0.0002177 0.0002321 0.0002768 0.00% MLMG::compResidual() 1 0.000251 0.0002564 0.0002623 0.00% DistributionMapping::Distribute() 95 0.0002229 0.0002301 0.0002381 0.00% OwnerMask() 2 0.0001375 0.000163 0.0002295 0.00% MLLinOp::define() 19 0.0001564 0.0001859 0.000216 0.00% Castro::apply_tagging_restrictions() 2 0.0001393 0.0001563 0.000182 0.00% Castro::finalize_do_advance() 17 0.0001426 0.0001573 0.0001802 0.00% Castro::apply_problem_tags() 2 0.0001599 0.0001634 0.0001743 0.00% Castro::check_for_nan() 37 0.0001215 0.000133 0.0001742 0.00% Amr::grid_places() 2 8.322e-05 0.0001199 0.0001724 0.00% knapsack() 1 5.904e-05 7.164e-05 0.0001642 0.00% Amr::initSubcycle() 1 5.748e-05 6.562e-05 0.0001639 0.00% TagBoxArray::mapPRD 2 0.0001133 0.0001216 0.0001455 0.00% Amr::coarseTimeStep() 1 0.0001142 0.0001389 0.0001419 0.00% MLPoisson::define() 19 0.000114 0.0001239 0.000138 0.00% Gravity::actual_multilevel_solve() 1 0.0001332 0.0001352 0.000138 0.00% Castro::construct_old_gravity() 20 9.546e-05 0.0001063 0.0001279 0.00% Castro::apply_source_to_state() 37 0.000106 0.0001149 0.0001264 0.00% Castro::enforce_speed_limit() 79 0.0001031 0.0001085 0.0001156 0.00% Amr::FinalizeInit() 1 8.891e-06 2.568e-05 9.074e-05 0.00% MultiFab::Swap() 5 7.106e-05 7.472e-05 9.061e-05 0.00% MLPoisson::prepareForSolve() 19 6.852e-05 7.404e-05 8.007e-05 0.00% AmrLevel::AmrLevel(dm) 1 3.037e-05 3.595e-05 5.686e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 4.896e-05 5.041e-05 5.251e-05 0.00% MLMG::computeMLResidual() 18 4.198e-05 4.635e-05 5.241e-05 0.00% MLMG::buildFineMask() 18 4.277e-05 4.654e-05 5.116e-05 0.00% DistributionMapping::KnapSackDoIt() 1 1.712e-05 2.118e-05 4.96e-05 0.00% Castro::finalize_advance() 1 1.988e-05 2.766e-05 4.734e-05 0.00% MLMG::getGradSolution() 18 4.098e-05 4.297e-05 4.556e-05 0.00% Gravity::swapTimeLevels() 8 3.475e-05 3.727e-05 4.052e-05 0.00% Castro::estTimeStep() 19 2.744e-05 3.043e-05 3.751e-05 0.00% MLMG::MLResNormInf() 18 3.347e-05 3.548e-05 3.726e-05 0.00% Castro::construct_old_rotation_source() 20 3.043e-05 3.207e-05 3.487e-05 0.00% Amr::regrid() 1 1.166e-05 1.544e-05 3.348e-05 0.00% Castro::construct_new_rotation_source() 17 2.204e-05 2.409e-05 2.535e-05 0.00% Amr::timeStep() 1 1.259e-05 1.379e-05 1.923e-05 0.00% Amr::bldFineLevels() 1 7.659e-06 9.42e-06 1.757e-05 0.00% Gravity::set_mass_offset() 2 4.016e-06 8.587e-06 1.636e-05 0.00% Castro::post_regrid() 2 1.094e-05 1.197e-05 1.603e-05 0.00% Amr::InitializeInit() 1 6.144e-06 8.027e-06 1.571e-05 0.00% Amr::init() 1 7.529e-06 8.843e-06 1.382e-05 0.00% Castro::computeInitialDt() 2 7.044e-06 9.911e-06 1.316e-05 0.00% Amr::initialInit() 1 2.854e-06 5.427e-06 1.147e-05 0.00% BoxList::parallelComplementIn() 4 7.232e-06 7.542e-06 7.929e-06 0.00% Castro::advance() 1 6.056e-06 6.588e-06 7.042e-06 0.00% Castro::errorEst() 2 5.669e-06 6.047e-06 6.848e-06 0.00% BoxList::complementIn 4 6.219e-06 6.408e-06 6.626e-06 0.00% Castro::gwstrain() 2 4.24e-06 4.737e-06 6.256e-06 0.00% Castro::FluxRegFineAdd() 1 1.589e-06 1.882e-06 2.788e-06 0.00% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.46 14.46 14.46 100.00% Amr::coarseTimeStep() 1 12.18 12.18 12.18 84.26% Amr::timeStep() 1 12.13 12.13 12.13 83.87% Castro::advance() 1 12.02 12.03 12.03 83.21% Castro::subcycle_advance_ctu() 1 11.97 11.98 11.98 82.87% Castro::do_advance_ctu() 20 11.93 11.93 11.93 82.55% Castro::construct_ctu_hydro_source() 20 5.286 5.403 5.508 38.10% Castro::clean_state() 79 3.61 3.681 3.801 26.29% Castro::computeTemp() 80 2.4 2.491 2.673 18.49% Castro::reset_internal_energy() 80 1.471 1.526 1.646 11.38% Gravity::solve_phi_with_mlmg() 19 1.398 1.449 1.484 10.26% Castro::construct_new_gravity() 17 1.331 1.382 1.417 9.80% Amr::InitAmr() 1 1.385 1.385 1.386 9.59% Gravity::solve_for_phi() 17 1.295 1.346 1.381 9.55% Castro::do_old_sources() 20 1.201 1.245 1.281 8.86% Castro::initialize_do_advance() 20 1.233 1.248 1.26 8.71% Gravity::actual_solve_with_mlmg() 19 1.196 1.2 1.204 8.33% Castro::do_new_sources() 17 1.128 1.155 1.187 8.21% MLMG::solve() 18 1.146 1.149 1.15 7.96% Castro::evaluate_source_change() 196 1.052 1.068 1.121 7.75% MLMG::oneIter() 128 1.027 1.028 1.029 7.12% MLMG::mgVcycle() 128 1.013 1.013 1.015 7.02% FabArray::ParallelCopy() 1583 0.4721 0.7867 0.9763 6.75% Amr::init() 1 0.889 0.8894 0.8894 6.15% Amr::initialInit() 1 0.8337 0.8341 0.8342 5.77% Castro::enforce_min_density() 79 0.7771 0.7957 0.828 5.73% MultiFab::AddProduct() 3920 0.6069 0.6174 0.6749 4.67% MLCellLinOp::smooth() 3072 0.486 0.5309 0.6466 4.47% FabArray::ParallelCopy_finish() 1583 0.1495 0.4577 0.6368 4.40% Amr::InitializeInit() 1 0.4388 0.4936 0.5862 4.05% Amr::defBaseLevel() 1 0.4388 0.4936 0.5862 4.05% FillPatchIterator::Initialize 301 0.4235 0.4719 0.524 3.62% FillPatchSingleLevel 301 0.3865 0.4425 0.4969 3.44% FabArray::setVal() 5502 0.4534 0.4598 0.4657 3.22% Castro::normalize_species() 79 0.4445 0.449 0.4529 3.13% MLMG::mgVcycle_up::1 128 0.08635 0.2335 0.4094 2.83% Amr::FinalizeInit() 1 0.2479 0.3406 0.3953 2.73% MLMG::addInterpCorrection() 768 0.1401 0.289 0.3624 2.51% MLPoisson::Fsmooth() 6144 0.2653 0.2902 0.35 2.42% FabArray::ParallelCopy_nowait() 1583 0.3024 0.3251 0.3476 2.40% MultiFab::Subtract() 159 0.3351 0.3383 0.3415 2.36% MLCellLinOp::applyBC() 7156 0.2535 0.2837 0.3411 2.36% Gravity::fill_multipole_BCs() 19 0.1967 0.2473 0.278 1.92% MLMG::mgVcycle_up::2 128 0.009636 0.1137 0.2748 1.90% MultiFab::min() 23 0.002079 0.1347 0.2718 1.88% FabArray::FillBoundary() 6388 0.1964 0.223 0.263 1.82% MLMG::mgVcycle_down::0 128 0.2362 0.2376 0.2396 1.66% Castro::initData() 1 0.07704 0.1255 0.224 1.55% MultiFab::sum() 3920 0.2113 0.2137 0.2162 1.50% Castro::construct_new_source() 102 0.2096 0.212 0.2151 1.49% MultiFab::Saxpy() 91 0.201 0.203 0.2059 1.42% Castro::expand_state() 20 0.1778 0.1854 0.1967 1.36% Amr::bldFineLevels() 1 0.04438 0.137 0.1917 1.33% MLMG::mgVcycle_up::3 128 0.009593 0.02288 0.186 1.29% MLMG::mgVcycle_up::0 128 0.1722 0.1727 0.1737 1.20% MLCellLinOp::apply() 994 0.1443 0.1522 0.1711 1.18% amrex::average_down 768 0.03302 0.0475 0.166 1.15% Castro::post_init() 1 0.1613 0.1619 0.1624 1.12% Castro::create_source_corrector() 17 0.1299 0.139 0.1554 1.08% Amr::grid_places() 2 0.002715 0.09535 0.1502 1.04% AmrMesh::MakeNewGrids() 2 0.002586 0.09523 0.15 1.04% TagBoxArray::mapPRD 2 0.0007011 0.09106 0.1489 1.03% FillBoundary_nowait() 6388 0.08884 0.1174 0.1469 1.02% Castro::derive() 156 0.1387 0.1407 0.1442 1.00% MLMG::mgVcycle_down::5 128 0.01165 0.02089 0.1441 1.00% Castro::construct_old_source() 120 0.1413 0.1425 0.1439 1.00% MLMG:computeResOfCorrection() 768 0.1093 0.1183 0.1403 0.97% MLCellLinOp::correctionResidual() 776 0.1078 0.1171 0.1387 0.96% FillBoundary_finish() 6388 0.07308 0.09185 0.123 0.85% MLMG::mgVcycle_down::1 128 0.105 0.1082 0.111 0.77% MLMG::mgVcycle_up::5 128 0.009702 0.01607 0.1063 0.74% MLPoisson::Fapply() 994 0.09243 0.09729 0.1059 0.73% Castro::post_timestep() 1 0.0959 0.0982 0.1005 0.70% Amr::writePlotFile() 2 0.09607 0.09616 0.09616 0.67% Gravity::multilevel_solve_for_new_phi() 1 0.09227 0.0923 0.09232 0.64% Gravity::actual_multilevel_solve() 1 0.09222 0.09225 0.09227 0.64% MLCellLinOp::solutionResidual() 147 0.08383 0.08808 0.09107 0.63% MLMG::mgVcycle_bottom 128 0.0003656 0.00585 0.08791 0.61% MLMG::actualBottomSolve() 128 0.0001463 0.005618 0.08763 0.61% MLCGSolver::bicgstab 8 0 0.005343 0.08549 0.59% Castro::post_regrid() 2 0.08275 0.08328 0.08389 0.58% Castro::apply_source_to_state() 37 0.08192 0.08275 0.08361 0.58% Castro::construct_old_gravity() 20 0.04379 0.06575 0.08292 0.57% Gravity::get_old_grav_vector() 20 0.04369 0.06564 0.0828 0.57% Castro::construct_new_gravity_source() 17 0.07934 0.08055 0.08209 0.57% Castro::problem_post_timestep() 2 0.07304 0.07599 0.07906 0.55% MLMG::computeResidual() 128 0.07249 0.0737 0.07419 0.51% MLMG::mgVcycle_down::2 128 0.0124 0.03874 0.06746 0.47% Castro::construct_old_gravity_source() 20 0.06337 0.06412 0.06502 0.45% MultiFab::Xpay() 923 0.04967 0.05206 0.0569 0.39% MLMG::mgVcycle_down::3 128 0.01168 0.01756 0.05585 0.39% Castro::volWgtSum() 66 0.05243 0.053 0.05359 0.37% Castro::initialize_advance() 1 0.0491 0.04992 0.05053 0.35% Castro::check_for_nan() 37 0.04466 0.04527 0.04574 0.32% MultiFab::contains_nan() 37 0.04454 0.04513 0.04558 0.32% Castro::retry_advance_ctu() 11 0.04373 0.04405 0.04433 0.31% StateDataPhysBCFunct::() 301 0.001145 0.02338 0.04393 0.30% Amr::regrid() 1 0.04243 0.04263 0.04289 0.30% StateData::FillBoundary(geom) 143 0 0.02209 0.04264 0.29% MLPoisson::define() 19 0.03005 0.03463 0.04045 0.28% Castro::wd_update() 2 0.03091 0.03384 0.03693 0.26% Gravity::get_new_grav_vector() 18 0.03536 0.03561 0.03613 0.25% MLMG::mgVcycle_down::4 128 0.01169 0.01335 0.03306 0.23% VisMF::Write(FabArray) 4 0.01294 0.01439 0.02678 0.19% MLMG::mgVcycle_up::4 128 0.009658 0.01085 0.02459 0.17% MLCellLinOp::defineAuxData() 19 0.01999 0.02135 0.02445 0.17% Castro::volInBoundary() 14 0.01998 0.02051 0.0226 0.16% TagBoxArray::collate() 2 6.667e-05 0.003151 0.02087 0.14% MLMG::prepareForSolve() 18 0.01689 0.01794 0.01966 0.14% MLMG::computeMLResidual() 18 0.01029 0.01399 0.01696 0.12% FabArray::setDomainBndry() 301 0.0008409 0.008875 0.0165 0.11% CGSolver::sxay() 203 0 0.0009844 0.01575 0.11% MLMG::ResNormInf() 146 0.01068 0.01275 0.01443 0.10% Amr::writeSmallPlotFile() 2 0.01442 0.01442 0.01442 0.10% Gravity::test_composite_phi() 1 0.0136 0.01363 0.01364 0.09% FabArrayBase::getCPC() 2005 0.009225 0.0106 0.0135 0.09% FabArrayBase::getFB() 6388 0.009849 0.01083 0.01278 0.09% MLLinOp::define() 19 0.003349 0.007053 0.01133 0.08% MLLinOp::defineGrids() 19 0.003151 0.006867 0.01114 0.08% MLCGSolver::ParallelAllReduce 189 0 0.0006445 0.01031 0.07% MultiFab::LinComb() 203 0 0.000615 0.00984 0.07% MLLinOp::makeSubCommunicator() 19 0.0003164 0.00392 0.008092 0.06% MLMG::getGradSolution() 18 0.007216 0.007457 0.007693 0.05% MLCellLinOp::compGrad() 18 0.007174 0.007414 0.007649 0.05% FabArray::mult() 59 0.007228 0.00735 0.00757 0.05% MLCellLinOp::defineBC() 19 0.005789 0.006101 0.006347 0.04% MultiFab::Dot() 142 0 0.000379 0.006065 0.04% FabArrayBase::CPC::define() 576 0.002193 0.003376 0.006054 0.04% Castro::locWgtSum() 12 0.005016 0.005061 0.005159 0.04% BndryData::define() 19 0.004121 0.004374 0.004567 0.03% MLPoisson::normalize() 79 0 0.0001919 0.00307 0.02% Castro::volProductSum() 4 0.002902 0.002956 0.003003 0.02% Gravity::update_max_rhs() 2 0.001386 0.002021 0.002536 0.02% MLPoisson::prepareForSolve() 19 0.001823 0.001977 0.002163 0.01% MLCellLinOp::prepareForSolve() 19 0.001744 0.001903 0.002094 0.01% MultiFab::Copy() 18 0.001988 0.002022 0.002055 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.001301 0.001552 0.001795 0.01% AmrLevel::AmrLevel(dm) 1 0.001434 0.001471 0.001673 0.01% MultiFab::max() 8 0.0005198 0.001156 0.00165 0.01% StateData::define() 6 0.001403 0.001435 0.001616 0.01% FabArrayBase::FB::FB() 75 0.000273 0.000681 0.001286 0.01% MLMG::MLResNormInf() 18 0.001171 0.001203 0.001235 0.01% Gravity::swapTimeLevels() 8 0.001127 0.00117 0.001211 0.01% MLMG::MLRhsNormInf() 18 0.001135 0.001155 0.001178 0.01% MLMG::compResidual() 1 0.000918 0.001063 0.001177 0.01% MLLinOp::makeAgglomeratedDMap 19 0.0009982 0.001043 0.001077 0.01% Castro::swap_state_time_levels() 8 0.0005148 0.0009672 0.001012 0.01% MLCellLinOp::setLevelBC() 19 0.0005339 0.000689 0.0008172 0.01% DistributionMapping::KnapSackDoIt() 1 0.0002844 0.0006971 0.0007459 0.01% makeSFC 95 0.0006215 0.0006572 0.0006882 0.00% DistributionMapping::LeastUsedCPUs() 1 7.06e-05 0.0006043 0.0006493 0.00% Gravity::applyMetricTerms() 2 0.0004626 0.0004684 0.0004742 0.00% Castro::enforce_consistent_e() 1 0.0003611 0.0003705 0.0003959 0.00% Castro::errorEst() 2 0.0003084 0.0003257 0.0003511 0.00% DistributionMapping::Distribute() 95 0.0002229 0.0002301 0.0002381 0.00% OwnerMask() 2 0.0001375 0.000163 0.0002295 0.00% Castro::apply_tagging_restrictions() 2 0.0001393 0.0001563 0.000182 0.00% Castro::finalize_do_advance() 17 0.0001426 0.0001573 0.0001802 0.00% Castro::apply_problem_tags() 2 0.0001599 0.0001634 0.0001743 0.00% knapsack() 1 5.904e-05 7.164e-05 0.0001642 0.00% Amr::initSubcycle() 1 5.748e-05 6.562e-05 0.0001639 0.00% Castro::enforce_speed_limit() 79 0.0001031 0.0001085 0.0001156 0.00% MultiFab::Swap() 5 7.106e-05 7.472e-05 9.061e-05 0.00% MLMG::buildFineMask() 18 4.277e-05 4.654e-05 5.116e-05 0.00% Castro::finalize_advance() 1 2.162e-05 2.954e-05 4.898e-05 0.00% Castro::estTimeStep() 19 2.744e-05 3.043e-05 3.751e-05 0.00% Castro::construct_old_rotation_source() 20 3.043e-05 3.207e-05 3.487e-05 0.00% Castro::construct_new_rotation_source() 17 2.204e-05 2.409e-05 2.535e-05 0.00% Gravity::set_mass_offset() 2 4.016e-06 8.587e-06 1.636e-05 0.00% Castro::computeInitialDt() 2 1.072e-05 1.346e-05 1.625e-05 0.00% BoxList::parallelComplementIn() 4 1.345e-05 1.395e-05 1.455e-05 0.00% BoxList::complementIn 4 6.219e-06 6.408e-06 6.626e-06 0.00% Castro::gwstrain() 2 4.24e-06 4.737e-06 6.256e-06 0.00% Castro::FluxRegFineAdd() 1 1.589e-06 1.882e-06 2.788e-06 0.00% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::castro.center(nvals = 2) :: [0.0e0, 0.0e0] [TOP]::castro.dtnuc_mode(nvals = 1) :: [1] [TOP]::integrator.do_constant_volume_burn(nvals = 1) :: [1] [TOP]::mg.v(nvals = 1) :: [0] AMReX (22.01-39-g874ccc6f3f8d) finalized