MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (22.01-18-g405b490356e5) initialized Starting run at 08:07:17 UTC on 2022-01-12. Successfully read inputs file ... Castro git describe: 21.12-23-g727e39bb6 AMReX git describe: 22.01-18-g405b49035 Microphysics git describe: 22.01-14-geeef682d 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.149192798 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001009505 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.01143598557 Gravity::fill_multipole_BCs() time = 0.01142692566 ... 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.064582846 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.008227002 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.07521605492 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.2947971821 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.07487225533 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.294907093 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.07257604599 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.2952935696 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.07161045074 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.3161077499 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01479411125 Gravity::solve_for_phi() time = 0.1101796627 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.07688951492 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.07090544701 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.2948725224 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01427268982 Gravity::solve_for_phi() time = 0.1098246574 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.07658600807 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.07081246376 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.299333334 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01438212395 Gravity::solve_for_phi() time = 0.1092765331 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.07635974884 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.07054138184 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.3005955219 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.01524162292 Gravity::solve_for_phi() time = 0.110480547 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.07578516006 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.07023692131 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.305524826 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.01477265358 Gravity::solve_for_phi() time = 0.1101891994 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.0769302845 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.0704369545 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.300665617 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.01574206352 Gravity::solve_for_phi() time = 0.1313610077 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.07782053947 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.07022333145 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.3053565025 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.0149371624 Gravity::solve_for_phi() time = 0.1105959415 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.07640242577 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.07076406479 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.3010056019 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.01600790024 Gravity::solve_for_phi() time = 0.1119582653 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.07565879822 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.07030391693 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.306112051 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.01506018639 Gravity::solve_for_phi() time = 0.11074543 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.07636237144 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.07045340538 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.3112981319 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.01604747772 Gravity::solve_for_phi() time = 0.1139543056 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.07623457909 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.07044005394 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.3066015244 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.01525449753 Gravity::solve_for_phi() time = 0.1110956669 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.07643723488 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.06983375549 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.3007442951 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.01576900482 Gravity::solve_for_phi() time = 0.1114058495 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.07595849037 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.07024407387 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.3062412739 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.01524615288 Gravity::solve_for_phi() time = 0.1111328602 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.07645893097 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.07009673119 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.3016545773 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.01607847214 Gravity::solve_for_phi() time = 0.1116688251 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.0759601593 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.07010316849 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.3060510159 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.01547574997 Gravity::solve_for_phi() time = 0.1111228466 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.07643675804 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.07025694847 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.3020327091 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.01629209518 Gravity::solve_for_phi() time = 0.111702919 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.0757598877 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.07076454163 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.3061521053 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.01575994492 Gravity::solve_for_phi() time = 0.1115438938 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.07683181763 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002404174757 [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: 13.7857958 [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.05043177 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.006841412 seconds Ending run at 08:07:33 UTC on 2022-01-12. Run time = 16.36302519 Run time without initialization = 13.84577036 Average number of zones advanced per microsecond: 0.002 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 16.36 ... 16.36 ... 16.36 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.9026 5.9611 6.0494 36.97% Amr::InitAmr() 1 1.524 1.524 1.524 9.31% Castro::computeTemp() 80 0.8879 0.9969 1.173 7.17% Castro::reset_internal_energy() 80 0.934 0.9711 1.081 6.61% FabArray::ParallelCopy_finish() 1583 0.2843 0.619 0.8417 5.14% MultiFab::AddProduct() 3920 0.695 0.7087 0.7161 4.38% Castro::normalize_species() 79 0.5117 0.5146 0.5198 3.18% FabArray::setVal() 5502 0.5054 0.5112 0.5162 3.15% MLPoisson::Fsmooth() 6144 0.3134 0.3413 0.4161 2.54% FabArray::ParallelCopy_nowait() 1583 0.3422 0.362 0.3858 2.36% Amr::defBaseLevel() 1 0.3677 0.377 0.3842 2.35% MultiFab::Subtract() 159 0.3622 0.3643 0.3667 2.24% Gravity::fill_multipole_BCs() 19 0.2077 0.2616 0.2822 1.72% MultiFab::sum() 3920 0.2381 0.2405 0.2421 1.48% Castro::enforce_min_density() 79 0.2008 0.2104 0.2379 1.45% MultiFab::Saxpy() 91 0.213 0.2151 0.2177 1.33% MultiFab::min() 23 0.00212 0.1328 0.2086 1.27% FillBoundary_nowait() 6388 0.1175 0.1469 0.1806 1.10% Castro::do_advance_ctu() 20 0.1649 0.1671 0.1698 1.04% FillBoundary_finish() 6388 0.1085 0.1239 0.148 0.90% MLCellLinOp::applyBC() 7156 0.08084 0.09828 0.1362 0.83% Castro::initData() 1 0.003908 0.02354 0.1336 0.82% MLPoisson::Fapply() 994 0.1113 0.1165 0.1291 0.79% Castro::derive() 156 0.1025 0.1051 0.1147 0.70% Castro::construct_new_gravity_source() 17 0.09357 0.09445 0.09534 0.58% Castro::construct_old_gravity_source() 20 0.07035 0.07209 0.08353 0.51% FabArray::FillBoundary() 6388 0.06764 0.06933 0.08174 0.50% Castro::do_old_sources() 20 0.001004 0.0175 0.07986 0.49% Castro::expand_state() 20 0.07228 0.07332 0.07446 0.46% Castro::do_new_sources() 17 0.0008469 0.01443 0.06578 0.40% MultiFab::Xpay() 923 0.05692 0.05928 0.06511 0.40% Castro::create_source_corrector() 17 0.05528 0.05613 0.05724 0.35% MLCellLinOp::smooth() 3072 0.05392 0.05508 0.05597 0.34% MLCGSolver::bicgstab 8 0 0.003371 0.05393 0.33% MultiFab::contains_nan() 37 0.05101 0.05127 0.05179 0.32% Castro::post_timestep() 1 0.05023 0.05047 0.05066 0.31% Castro::evaluate_source_change() 196 0.04779 0.04857 0.04983 0.30% StateData::FillBoundary(geom) 143 0 0.02464 0.04831 0.30% MLMG::addInterpCorrection() 768 0.03539 0.03803 0.04538 0.28% VisMF::Write(FabArray) 4 0.01453 0.01713 0.04362 0.27% Amr::writePlotFile() 2 0.01914 0.03984 0.0435 0.27% amrex::average_down 768 0.03229 0.03459 0.04043 0.25% FabArrayBase::getFB() 6388 0.02428 0.02538 0.03118 0.19% Castro::retry_advance_ctu() 11 0.02571 0.02599 0.02638 0.16% CGSolver::sxay() 203 0 0.001532 0.02451 0.15% FillPatchIterator::Initialize 301 0.02291 0.02321 0.0239 0.15% MLCellLinOp::apply() 994 0.009633 0.01075 0.02375 0.15% TagBoxArray::collate() 2 8.858e-05 0.003288 0.02325 0.14% Gravity::get_old_grav_vector() 20 0.02077 0.02104 0.02133 0.13% MLCellLinOp::defineAuxData() 19 0.01867 0.01927 0.02093 0.13% Gravity::get_new_grav_vector() 18 0.01898 0.01938 0.02005 0.12% Castro::construct_new_source() 102 0.01923 0.01946 0.01981 0.12% MultiFab::LinComb() 203 0 0.001218 0.01949 0.12% FabArray::setDomainBndry() 301 0.001626 0.01049 0.01928 0.12% MLCGSolver::ParallelAllReduce 189 0 0.001162 0.01859 0.11% MLMG::oneIter() 128 0.01608 0.01637 0.01761 0.11% FabArray::ParallelCopy() 1583 0.01638 0.01687 0.01744 0.11% MLMG::ResNormInf() 146 0.01439 0.01578 0.01736 0.11% FabArrayBase::getCPC() 2005 0.01391 0.01432 0.01524 0.09% Amr::writeSmallPlotFile() 2 0.008067 0.0127 0.01333 0.08% Castro::construct_old_source() 120 0.013 0.0131 0.0132 0.08% MultiFab::Dot() 142 0 0.000807 0.01291 0.08% Castro::wd_update() 2 0.004354 0.008821 0.01141 0.07% Castro::volWgtSum() 66 0.01081 0.01089 0.01095 0.07% MLCellLinOp::correctionResidual() 776 0.008155 0.008372 0.009669 0.06% MLMG::prepareForSolve() 18 0.008363 0.008634 0.009293 0.06% FabArray::mult() 59 0.007958 0.00817 0.008539 0.05% Gravity::actual_solve_with_mlmg() 19 0.007081 0.007345 0.007925 0.05% MLLinOp::makeSubCommunicator() 19 0.000327 0.003966 0.007877 0.05% FabArrayBase::CPC::define() 576 0.003423 0.004752 0.007782 0.05% MLCellLinOp::compGrad() 18 0.006595 0.006653 0.006746 0.04% MLMG::mgVcycle() 128 0.006358 0.006515 0.006674 0.04% MLPoisson::normalize() 79 0 0.0004125 0.0066 0.04% MLMG:computeResOfCorrection() 768 0.005487 0.005555 0.005699 0.03% MultiFab::max() 8 0.001238 0.003092 0.005578 0.03% MLMG::solve() 18 0.00327 0.004067 0.005283 0.03% Castro::initialize_do_advance() 20 0.004436 0.004634 0.005084 0.03% Castro::volInBoundary() 14 0.003974 0.00405 0.004368 0.03% BndryData::define() 19 0.003698 0.003791 0.004066 0.02% Gravity::solve_for_phi() 17 0.002759 0.002952 0.003814 0.02% FillPatchSingleLevel 301 0.003301 0.003384 0.003461 0.02% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.002774 0.003071 0.003451 0.02% MLMG::mgVcycle_down::0 128 0.002759 0.002877 0.002988 0.02% MLMG::mgVcycle_down::2 128 0.002672 0.002802 0.002983 0.02% Castro::construct_new_gravity() 17 0.002499 0.002586 0.002883 0.02% StateDataPhysBCFunct::() 301 0.001978 0.002508 0.002873 0.02% MLMG::mgVcycle_down::1 128 0.002632 0.002765 0.002859 0.02% MLMG::mgVcycle_down::3 128 0.002674 0.002736 0.002845 0.02% MLMG::mgVcycle_down::4 128 0.002629 0.002697 0.002832 0.02% MLMG::mgVcycle_down::5 128 0.002689 0.002738 0.002826 0.02% MLMG::actualBottomSolve() 128 0.0004355 0.0006068 0.002812 0.02% MLLinOp::defineGrids() 19 0.001985 0.002086 0.002384 0.01% MultiFab::Copy() 18 0.002178 0.002243 0.002322 0.01% MLCellLinOp::prepareForSolve() 19 0.001844 0.002028 0.002296 0.01% MLMG::mgVcycle_up::2 128 0.001779 0.001897 0.002043 0.01% MLMG::mgVcycle_up::5 128 0.001782 0.001835 0.002033 0.01% MLCellLinOp::defineBC() 19 0.001787 0.001845 0.002015 0.01% MLMG::mgVcycle_up::0 128 0.001756 0.001858 0.001996 0.01% Castro::locWgtSum() 12 0.001915 0.001943 0.001994 0.01% MLMG::mgVcycle_up::3 128 0.001791 0.00186 0.00199 0.01% MLMG::mgVcycle_up::1 128 0.0018 0.001868 0.001953 0.01% MLMG::mgVcycle_up::4 128 0.001784 0.001837 0.001907 0.01% StateData::define() 6 0.001473 0.001535 0.001844 0.01% MLCellLinOp::solutionResidual() 147 0.001534 0.001609 0.001683 0.01% Castro::clean_state() 79 0.001525 0.001573 0.001674 0.01% FabArrayBase::FB::FB() 75 0.0004389 0.0008825 0.001529 0.01% Castro::subcycle_advance_ctu() 1 0.0004298 0.0005338 0.001494 0.01% Castro::problem_post_timestep() 2 0.0004399 0.001021 0.001437 0.01% MLMG::MLRhsNormInf() 18 0.001327 0.001356 0.001395 0.01% Castro::post_init() 1 0.0004608 0.0009527 0.001375 0.01% Castro::swap_state_time_levels() 8 0.0005594 0.001052 0.001233 0.01% main() 1 0.0008615 0.001007 0.001147 0.01% makeSFC 95 0.0009006 0.0009672 0.001043 0.01% MLMG::computeResidual() 128 0.000906 0.0009514 0.001009 0.01% MLMG::mgVcycle_bottom 128 0.0008812 0.0009195 0.0009818 0.01% MLCellLinOp::setLevelBC() 19 0.0006249 0.0007542 0.0009265 0.01% Castro::volProductSum() 4 0.0008389 0.0008784 0.0009098 0.01% Gravity::applyMetricTerms() 2 0.0005143 0.0005676 0.0008167 0.00% MLLinOp::makeAgglomeratedDMap 19 0.0006696 0.0007176 0.000804 0.00% AmrMesh::MakeNewGrids() 2 0.0005708 0.0005839 0.0006062 0.00% Gravity::update_max_rhs() 2 0.0004518 0.0004837 0.0005966 0.00% Gravity::test_composite_phi() 1 0.0004525 0.0005062 0.000584 0.00% DistributionMapping::Distribute() 95 0.0004135 0.0004536 0.0005126 0.00% Gravity::solve_phi_with_mlmg() 19 0.0004307 0.0004572 0.0005024 0.00% Castro::enforce_consistent_e() 1 0.0004061 0.0004245 0.0004468 0.00% Castro::initialize_advance() 1 0.0002761 0.0003201 0.0003854 0.00% MLPoisson::define() 19 0.0003097 0.0003272 0.0003523 0.00% Castro::check_for_nan() 37 0.000313 0.0003318 0.0003516 0.00% Castro::enforce_speed_limit() 79 0.0003052 0.000327 0.000342 0.00% Castro::apply_source_to_state() 37 0.000297 0.0003129 0.0003276 0.00% MLMG::compResidual() 1 0.0002723 0.0002819 0.0002898 0.00% Castro::construct_old_gravity() 20 0.0002043 0.0002195 0.0002389 0.00% MLLinOp::define() 19 0.0002003 0.0002162 0.0002306 0.00% Castro::finalize_do_advance() 17 0.0001742 0.0001981 0.0002261 0.00% Castro::apply_tagging_restrictions() 2 0.00016 0.0001827 0.0002124 0.00% DistributionMapping::LeastUsedCPUs() 1 2.2e-05 0.0001809 0.0002036 0.00% MLPoisson::prepareForSolve() 19 0.0001469 0.0001689 0.0002005 0.00% Castro::apply_problem_tags() 2 0.0001674 0.0001746 0.0001883 0.00% OwnerMask() 2 0.0001517 0.0001651 0.0001763 0.00% Amr::coarseTimeStep() 1 0.0001356 0.0001611 0.000169 0.00% MLMG::computeMLResidual() 18 0.000132 0.0001445 0.0001619 0.00% Gravity::actual_multilevel_solve() 1 0.0001468 0.0001506 0.0001552 0.00% MLMG::getGradSolution() 18 0.0001448 0.0001476 0.0001518 0.00% Castro::construct_old_rotation_source() 20 0.0001342 0.0001397 0.0001515 0.00% MLMG::MLResNormInf() 18 0.0001215 0.0001308 0.0001465 0.00% TagBoxArray::mapPRD 2 0.0001371 0.0001414 0.000146 0.00% Castro::construct_new_rotation_source() 17 0.0001119 0.0001179 0.0001313 0.00% Gravity::swapTimeLevels() 8 9.925e-05 0.0001073 0.0001312 0.00% Amr::grid_places() 2 9.954e-05 0.0001109 0.0001213 0.00% Amr::initSubcycle() 1 4.408e-05 5.391e-05 0.0001144 0.00% MLMG::buildFineMask() 18 8.225e-05 9.404e-05 0.0001117 0.00% MultiFab::Swap() 5 8.213e-05 9.115e-05 0.0001065 0.00% Castro::estTimeStep() 19 7.459e-05 8.198e-05 0.000101 0.00% knapsack() 1 6.221e-05 6.669e-05 8.383e-05 0.00% AmrLevel::AmrLevel(dm) 1 4.4e-05 5.322e-05 7.733e-05 0.00% Amr::FinalizeInit() 1 2.183e-05 2.633e-05 6.013e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.379e-05 5.63e-05 5.867e-05 0.00% Castro::finalize_advance() 1 2.808e-05 3.32e-05 5.775e-05 0.00% Amr::regrid() 1 1.825e-05 2.232e-05 4.304e-05 0.00% DistributionMapping::KnapSackDoIt() 1 2.662e-05 3.392e-05 4.258e-05 0.00% Amr::timeStep() 1 2.096e-05 2.42e-05 3.504e-05 0.00% BoxList::parallelComplementIn() 4 2.737e-05 2.892e-05 3.125e-05 0.00% Castro::post_regrid() 2 2.004e-05 2.242e-05 2.704e-05 0.00% Castro::advance() 1 1.696e-05 1.931e-05 2.504e-05 0.00% Amr::InitializeInit() 1 1.225e-05 1.517e-05 2.5e-05 0.00% Castro::errorEst() 2 2.117e-05 2.246e-05 2.4e-05 0.00% Amr::init() 1 1.942e-05 2.051e-05 2.254e-05 0.00% Castro::FluxRegFineAdd() 1 2.458e-06 5.156e-06 2.179e-05 0.00% Castro::computeInitialDt() 2 1.692e-05 1.833e-05 1.992e-05 0.00% Amr::bldFineLevels() 1 1.292e-05 1.646e-05 1.754e-05 0.00% BoxList::complementIn 4 1.492e-05 1.607e-05 1.679e-05 0.00% Castro::gwstrain() 2 6.292e-06 1.006e-05 1.463e-05 0.00% Amr::initialInit() 1 1.208e-05 1.257e-05 1.304e-05 0.00% Gravity::set_mass_offset() 2 7e-06 9.956e-06 1.233e-05 0.00% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 16.36 16.36 16.36 100.00% Amr::coarseTimeStep() 1 13.85 13.85 13.85 84.61% Amr::timeStep() 1 13.79 13.79 13.79 84.25% Castro::advance() 1 13.63 13.63 13.63 83.29% Castro::subcycle_advance_ctu() 1 13.56 13.57 13.57 82.96% Castro::do_advance_ctu() 20 13.51 13.52 13.52 82.65% Castro::construct_ctu_hydro_source() 20 5.903 5.961 6.049 36.97% Castro::clean_state() 79 3.937 4.019 4.189 25.60% Castro::computeTemp() 80 2.567 2.662 2.824 17.26% Gravity::solve_phi_with_mlmg() 19 1.97 2.028 2.051 12.53% Castro::construct_new_gravity() 17 1.867 1.923 1.946 11.89% Gravity::solve_for_phi() 17 1.826 1.883 1.906 11.65% Castro::reset_internal_energy() 80 1.626 1.665 1.769 10.81% Gravity::actual_solve_with_mlmg() 19 1.761 1.764 1.769 10.81% MLMG::solve() 18 1.698 1.701 1.703 10.41% MLMG::oneIter() 128 1.548 1.549 1.551 9.48% MLMG::mgVcycle() 128 1.531 1.533 1.535 9.38% Amr::InitAmr() 1 1.524 1.524 1.524 9.31% Castro::do_old_sources() 20 1.319 1.37 1.419 8.67% Castro::initialize_do_advance() 20 1.359 1.371 1.384 8.46% Castro::do_new_sources() 17 1.226 1.266 1.298 7.94% Castro::evaluate_source_change() 196 1.237 1.253 1.262 7.71% FabArray::ParallelCopy() 1583 0.6724 1.009 1.23 7.52% Amr::init() 1 0.9927 0.9929 0.9929 6.07% Castro::enforce_min_density() 79 0.8814 0.897 0.9247 5.65% Amr::initialInit() 1 0.9173 0.9175 0.9176 5.61% MLCellLinOp::smooth() 3072 0.725 0.7755 0.9066 5.54% FabArray::ParallelCopy_finish() 1583 0.2843 0.619 0.8417 5.14% MultiFab::AddProduct() 3920 0.695 0.7087 0.7161 4.38% FillPatchIterator::Initialize 301 0.4825 0.5375 0.6243 3.82% Amr::InitializeInit() 1 0.4594 0.5166 0.6069 3.71% Amr::defBaseLevel() 1 0.4594 0.5165 0.6069 3.71% MLMG::mgVcycle_up::1 128 0.1151 0.3335 0.5995 3.66% FillPatchSingleLevel 301 0.4399 0.5038 0.5936 3.63% MLMG::addInterpCorrection() 768 0.1674 0.4515 0.5423 3.31% MLCellLinOp::applyBC() 7156 0.4209 0.4647 0.52 3.18% Castro::normalize_species() 79 0.5117 0.5146 0.5198 3.18% FabArray::setVal() 5502 0.5054 0.5112 0.5162 3.15% Amr::FinalizeInit() 1 0.3106 0.4009 0.4581 2.80% MLMG::mgVcycle_up::2 128 0.03143 0.1937 0.4503 2.75% MLPoisson::Fsmooth() 6144 0.3134 0.3413 0.4161 2.54% FabArray::ParallelCopy_nowait() 1583 0.3528 0.3732 0.3979 2.43% FabArray::FillBoundary() 6388 0.3279 0.3664 0.3955 2.42% MultiFab::Subtract() 159 0.3622 0.3643 0.3667 2.24% MLMG::mgVcycle_up::3 128 0.0314 0.05377 0.3475 2.12% MLMG::mgVcycle_down::0 128 0.2911 0.2925 0.2946 1.80% Gravity::fill_multipole_BCs() 19 0.2077 0.2616 0.2822 1.72% MLCellLinOp::apply() 994 0.1973 0.2105 0.2741 1.68% MLMG::mgVcycle_up::5 128 0.03083 0.04577 0.2588 1.58% MultiFab::sum() 3920 0.2381 0.2405 0.2421 1.48% Castro::construct_new_source() 102 0.233 0.2352 0.238 1.45% Castro::expand_state() 20 0.1988 0.21 0.2333 1.43% MLMG::mgVcycle_bottom 128 0.001325 0.01503 0.2199 1.34% MLMG::actualBottomSolve() 128 0.0004355 0.01411 0.2189 1.34% Castro::initData() 1 0.08406 0.1318 0.2181 1.33% MultiFab::Saxpy() 91 0.213 0.2151 0.2177 1.33% Castro::post_init() 1 0.2128 0.2147 0.2171 1.33% MLCGSolver::bicgstab 8 0 0.01345 0.2152 1.31% MLMG::mgVcycle_up::0 128 0.2115 0.2118 0.2125 1.30% MultiFab::min() 23 0.00212 0.1328 0.2086 1.27% FillBoundary_nowait() 6388 0.1427 0.1732 0.2077 1.27% amrex::average_down 768 0.05267 0.06874 0.2035 1.24% MLMG:computeResOfCorrection() 768 0.1609 0.172 0.1979 1.21% Amr::bldFineLevels() 1 0.05029 0.1397 0.1957 1.20% MLCellLinOp::correctionResidual() 776 0.1554 0.1671 0.1924 1.18% MLMG::mgVcycle_down::5 128 0.03792 0.04859 0.1888 1.15% Castro::derive() 156 0.156 0.159 0.1693 1.03% Castro::construct_old_source() 120 0.1546 0.1567 0.1678 1.03% Castro::create_source_corrector() 17 0.1441 0.1506 0.161 0.98% Castro::post_timestep() 1 0.1561 0.1582 0.1604 0.98% MLMG::mgVcycle_down::1 128 0.1441 0.1472 0.151 0.92% Amr::grid_places() 2 0.002762 0.09308 0.1502 0.92% AmrMesh::MakeNewGrids() 2 0.002647 0.09297 0.1501 0.92% TagBoxArray::mapPRD 2 0.0007689 0.08847 0.1481 0.90% FillBoundary_finish() 6388 0.1085 0.1239 0.148 0.90% Gravity::multilevel_solve_for_new_phi() 1 0.1326 0.1326 0.1327 0.81% Gravity::actual_multilevel_solve() 1 0.1326 0.1326 0.1326 0.81% MLPoisson::Fapply() 994 0.1113 0.1165 0.1291 0.79% Amr::writePlotFile() 2 0.1176 0.1177 0.1177 0.72% MLCellLinOp::solutionResidual() 147 0.1026 0.1085 0.1115 0.68% MLMG::mgVcycle_down::2 128 0.03863 0.06778 0.09985 0.61% Castro::construct_new_gravity_source() 17 0.09357 0.09445 0.09534 0.58% Castro::post_regrid() 2 0.09066 0.09306 0.09498 0.58% Castro::construct_old_gravity() 20 0.04519 0.07129 0.09286 0.57% Gravity::get_old_grav_vector() 20 0.04499 0.07107 0.09264 0.57% MLMG::computeResidual() 128 0.08957 0.09072 0.09189 0.56% Castro::apply_source_to_state() 37 0.0866 0.08754 0.08838 0.54% MLMG::mgVcycle_down::3 128 0.03825 0.0446 0.08774 0.54% Castro::problem_post_timestep() 2 0.08201 0.08543 0.08759 0.54% Castro::construct_old_gravity_source() 20 0.07035 0.07209 0.08353 0.51% Castro::initialize_advance() 1 0.05465 0.0575 0.06718 0.41% MultiFab::Xpay() 923 0.05692 0.05928 0.06511 0.40% MLMG::mgVcycle_down::4 128 0.03715 0.03938 0.06216 0.38% Castro::volWgtSum() 66 0.06029 0.06086 0.06158 0.38% Castro::check_for_nan() 37 0.05133 0.0516 0.05212 0.32% MultiFab::contains_nan() 37 0.05101 0.05127 0.05179 0.32% StateDataPhysBCFunct::() 301 0.001978 0.02715 0.05107 0.31% MLPoisson::define() 19 0.0397 0.04466 0.05064 0.31% Amr::regrid() 1 0.04646 0.04763 0.04859 0.30% StateData::FillBoundary(geom) 143 0 0.02464 0.04831 0.30% Castro::retry_advance_ctu() 11 0.04706 0.04758 0.04813 0.29% MLMG::mgVcycle_up::4 128 0.03123 0.03277 0.04796 0.29% CGSolver::sxay() 203 0 0.00275 0.044 0.27% VisMF::Write(FabArray) 4 0.01453 0.01713 0.04362 0.27% Gravity::get_new_grav_vector() 18 0.03896 0.04004 0.04115 0.25% Castro::wd_update() 2 0.03418 0.03763 0.03981 0.24% FabArrayBase::getFB() 6388 0.02474 0.02627 0.0319 0.19% MLCellLinOp::defineAuxData() 19 0.02732 0.0288 0.03184 0.19% Castro::volInBoundary() 14 0.02333 0.02392 0.02559 0.16% MLCGSolver::ParallelAllReduce 189 0 0.001598 0.02558 0.16% MLMG::prepareForSolve() 18 0.02075 0.02182 0.02395 0.15% TagBoxArray::collate() 2 8.858e-05 0.003288 0.02325 0.14% FabArrayBase::getCPC() 2005 0.01763 0.01907 0.0217 0.13% MLMG::computeMLResidual() 18 0.01274 0.01803 0.02088 0.13% MultiFab::LinComb() 203 0 0.001218 0.01949 0.12% FabArray::setDomainBndry() 301 0.001626 0.01049 0.01928 0.12% MLMG::ResNormInf() 146 0.01439 0.01578 0.01736 0.11% Amr::writeSmallPlotFile() 2 0.01706 0.01706 0.01706 0.10% Gravity::test_composite_phi() 1 0.01566 0.01567 0.01569 0.10% MultiFab::Dot() 142 0 0.000807 0.01291 0.08% MLLinOp::define() 19 0.0046 0.008406 0.01269 0.08% MLLinOp::defineGrids() 19 0.004376 0.00819 0.01247 0.08% MLMG::getGradSolution() 18 0.008886 0.00909 0.009307 0.06% MLCellLinOp::compGrad() 18 0.008739 0.008942 0.009157 0.06% FabArray::mult() 59 0.007958 0.00817 0.008539 0.05% MLLinOp::makeSubCommunicator() 19 0.000327 0.003966 0.007877 0.05% FabArrayBase::CPC::define() 576 0.003423 0.004752 0.007782 0.05% MLCellLinOp::defineBC() 19 0.006888 0.00713 0.007724 0.05% MLPoisson::normalize() 79 0 0.0004125 0.0066 0.04% Gravity::update_max_rhs() 2 0.002264 0.00418 0.006585 0.04% Castro::locWgtSum() 12 0.005854 0.00592 0.005994 0.04% BndryData::define() 19 0.005083 0.005285 0.005708 0.03% MultiFab::max() 8 0.001238 0.003092 0.005578 0.03% Castro::volProductSum() 4 0.003469 0.003542 0.003601 0.02% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.002774 0.003071 0.003451 0.02% MLPoisson::prepareForSolve() 19 0.00201 0.002197 0.002472 0.02% MultiFab::Copy() 18 0.002178 0.002243 0.002322 0.01% MLCellLinOp::prepareForSolve() 19 0.001844 0.002028 0.002296 0.01% MLLinOp::makeAgglomeratedDMap 19 0.002003 0.002138 0.002282 0.01% AmrLevel::AmrLevel(dm) 1 0.00152 0.001589 0.001919 0.01% StateData::define() 6 0.001473 0.001535 0.001844 0.01% makeSFC 95 0.001323 0.001421 0.001555 0.01% MLMG::MLResNormInf() 18 0.001421 0.001465 0.001539 0.01% Gravity::swapTimeLevels() 8 0.001334 0.00141 0.001537 0.01% FabArrayBase::FB::FB() 75 0.0004389 0.0008825 0.001529 0.01% MLMG::MLRhsNormInf() 18 0.001327 0.001356 0.001395 0.01% MLMG::compResidual() 1 0.001096 0.001247 0.001384 0.01% Castro::swap_state_time_levels() 8 0.0005594 0.001052 0.001233 0.01% MLCellLinOp::setLevelBC() 19 0.0006249 0.0007542 0.0009265 0.01% Gravity::applyMetricTerms() 2 0.0005143 0.0005676 0.0008167 0.00% DistributionMapping::Distribute() 95 0.0004135 0.0004536 0.0005126 0.00% Castro::enforce_consistent_e() 1 0.0004061 0.0004245 0.0004468 0.00% Castro::errorEst() 2 0.0003498 0.0003798 0.0004166 0.00% Castro::enforce_speed_limit() 79 0.0003052 0.000327 0.000342 0.00% DistributionMapping::KnapSackDoIt() 1 0.0001389 0.0002815 0.0003104 0.00% Castro::finalize_do_advance() 17 0.0001742 0.0001981 0.0002261 0.00% Castro::apply_tagging_restrictions() 2 0.00016 0.0001827 0.0002124 0.00% DistributionMapping::LeastUsedCPUs() 1 2.2e-05 0.0001809 0.0002036 0.00% Castro::apply_problem_tags() 2 0.0001674 0.0001746 0.0001883 0.00% OwnerMask() 2 0.0001517 0.0001651 0.0001763 0.00% Castro::construct_old_rotation_source() 20 0.0001342 0.0001397 0.0001515 0.00% Castro::construct_new_rotation_source() 17 0.0001119 0.0001179 0.0001313 0.00% Amr::initSubcycle() 1 4.408e-05 5.391e-05 0.0001144 0.00% MLMG::buildFineMask() 18 8.225e-05 9.404e-05 0.0001117 0.00% MultiFab::Swap() 5 8.213e-05 9.115e-05 0.0001065 0.00% Castro::estTimeStep() 19 7.459e-05 8.198e-05 0.000101 0.00% knapsack() 1 6.221e-05 6.669e-05 8.383e-05 0.00% Castro::finalize_advance() 1 3.267e-05 3.836e-05 6.2e-05 0.00% BoxList::parallelComplementIn() 4 4.317e-05 4.499e-05 4.779e-05 0.00% Castro::computeInitialDt() 2 2.629e-05 2.81e-05 3.129e-05 0.00% Castro::FluxRegFineAdd() 1 2.458e-06 5.156e-06 2.179e-05 0.00% BoxList::complementIn 4 1.492e-05 1.607e-05 1.679e-05 0.00% Castro::gwstrain() 2 6.292e-06 1.006e-05 1.463e-05 0.00% Gravity::set_mass_offset() 2 7e-06 9.956e-06 1.233e-05 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-18-g405b490356e5) finalized