MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (22.03-34-g8b691cbf8ddb) initialized Starting run at 07:57:03 UTC on 2022-03-30. Successfully read inputs file ... Castro git describe: 22.02-19-ga92ac477d AMReX git describe: 22.03-34-g8b691cbf8 Microphysics git describe: 22.03-14-g4f1304ed 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 integrator.retry_burn = 1 Initializing the data at level 0 Done initializing the level 0 data grid_places() time: 0.436502685 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001507068 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.020739684 Gravity::fill_multipole_BCs() time = 0.021475474 ... 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.084608371 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.009429189 seconds [Level 0 step 1] ADVANCE with dt = 0.25 [1m[32m Beginning subcycle 1 starting at time 0 with dt = 0.25[0m 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.107173229 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.498593537 Advance was unsuccessful with reason: negative density (density = -9.081732e+02 at index (14,70); started at 4.477055e+03); proceeding to a retry. [1m[32m Subcycle completed[0m [1m[31m Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125 [0m [1m[32m Beginning subcycle 1 starting at time 0 with dt = 0.125[0m 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.106430231 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.502549773 Advance was unsuccessful with reason: negative density (density = -1.592544e+00 at index (20,79); started at 2.790161e+00); proceeding to a retry. [1m[32m Subcycle completed[0m [1m[31m Timestep 0.125 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.0625 [0m [1m[32m Beginning subcycle 1 starting at time 0 with dt = 0.0625[0m 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.105696545 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.498735078 Advance was unsuccessful with reason: negative density (density = -1.060745e-01 at index (20,78); started at 2.068339e-01); proceeding to a retry. [1m[32m Subcycle completed[0m [1m[31m Timestep 0.0625 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.03125 [0m [1m[32m Beginning subcycle 1 starting at time 0 with dt = 0.03125[0m 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.105202069 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.494822285 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.06348257 Gravity::solve_for_phi() time = 0.190565425 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.112958454 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.104934638 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.496355937 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.063878679 Gravity::solve_for_phi() time = 0.191044595 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.11288396 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.104983749 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.503653856 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.063324227 Gravity::solve_for_phi() time = 0.190395792 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.11232102 Ending SDC iteration 3 of 3. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 2 starting at time 0.03125 with dt = 0.03125[0m 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.104804839 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.507483124 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.064947698 Gravity::solve_for_phi() time = 0.192635822 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.112164297 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.103642642 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.516878402 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.065028169 Gravity::solve_for_phi() time = 0.19324268 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.112500856 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 3 starting at time 0.0625 with dt = 0.03125[0m 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.105390155 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.506813739 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.066611946 Gravity::solve_for_phi() time = 0.194542228 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.112233688 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.104283926 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.516091431 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.066795186 Gravity::solve_for_phi() time = 0.194128101 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.112209248 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 4 starting at time 0.09375 with dt = 0.03125[0m 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.105394021 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.515849906 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.07144611 Gravity::solve_for_phi() time = 0.198786417 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.11212783 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.103358508 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.516866485 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.066742609 Gravity::solve_for_phi() time = 0.194052152 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.112845565 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 5 starting at time 0.125 with dt = 0.03125[0m 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.104686297 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.505709929 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.067608858 Gravity::solve_for_phi() time = 0.195035299 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.113670043 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.104663843 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.515295187 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.06795053 Gravity::solve_for_phi() time = 0.19517304 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.113283004 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 6 starting at time 0.15625 with dt = 0.03125[0m 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.10539266 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.507161806 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.067339658 Gravity::solve_for_phi() time = 0.19450291 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.113058464 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.104142491 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.51680992 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.068181455 Gravity::solve_for_phi() time = 0.195671974 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.113831038 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 7 starting at time 0.1875 with dt = 0.03125[0m 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.104260512 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.512085802 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.066872517 Gravity::solve_for_phi() time = 0.191336588 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.111225245 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.1038743 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.51697165 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.067814895 Gravity::solve_for_phi() time = 0.194999204 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.113160398 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 8 starting at time 0.21875 with dt = 0.03125[0m 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.105084156 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.505859502 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.068445914 Gravity::solve_for_phi() time = 0.195761625 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.112902711 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.103618043 ... Entering construct_ctu_hydro_source() ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.515947192 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.069282645 Gravity::solve_for_phi() time = 0.196480907 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.113089219 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m Subcycling complete Zones advanced per microsecond at this level: 0.001519042417 [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: 21.75064452 [STEP 1] FAB kilobyte spread across MPI nodes: [7068 ... 7068] STEP = 1 TIME = 0.25 DT = 0.25 PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.084056393 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.009463385 seconds Ending run at 07:57:29 UTC on 2022-03-30. Run time = 25.83791226 Run time without initialization = 21.84656272 Average number of zones advanced per microsecond: 0.001 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 25.84 ... 25.84 ... 25.84 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 4.6994 7.3955 10.1874 39.43% MultiFab::min() 23 0.003435 2.86 5.586 21.62% FabArray::ParallelCopy_finish() 1583 0.18 2.489 4.304 16.66% Castro::reset_internal_energy() 80 0.9069 1.492 3.051 11.81% Amr::InitAmr() 1 2.525 2.526 2.526 9.78% Castro::computeTemp() 80 0.9115 1.335 2.033 7.87% Gravity::fill_multipole_BCs() 19 0.3969 0.8329 1.18 4.57% Castro::enforce_min_density() 79 0.1905 0.3723 1.158 4.48% Castro::normalize_species() 79 0.4854 0.7119 1.084 4.20% MultiFab::AddProduct() 3439 0.5568 0.7284 1.043 4.04% FillBoundary_finish() 6388 0.1407 0.593 0.8608 3.33% MLPoisson::Fsmooth() 6144 0.2821 0.4344 0.7697 2.98% Amr::defBaseLevel() 1 0.3675 0.4908 0.726 2.81% FabArray::setVal() 5021 0.3626 0.4492 0.6363 2.46% MultiFab::Subtract() 159 0.3358 0.4414 0.6315 2.44% FabArray::ParallelCopy_nowait() 1583 0.2298 0.3308 0.4984 1.93% MultiFab::sum() 3439 0.1874 0.256 0.4031 1.56% Castro::do_old_sources() 20 0.0002971 0.0347 0.2663 1.03% MLPoisson::Fapply() 994 0.09813 0.1465 0.2399 0.93% FillBoundary_nowait() 6388 0.08447 0.1394 0.2349 0.91% Castro::do_new_sources() 17 0.0002475 0.03064 0.2244 0.87% Castro::do_advance_ctu() 20 0.1158 0.1452 0.1932 0.75% Castro::derive() 156 0.09803 0.1303 0.1888 0.73% MLCellLinOp::applyBC() 7156 0.05402 0.08204 0.1704 0.66% Castro::construct_new_gravity_source() 17 0.07048 0.1049 0.1625 0.63% Castro::initData() 1 0.003543 0.02632 0.1356 0.52% Castro::wd_update() 2 0.007803 0.07447 0.128 0.50% Castro::expand_state() 20 0.06827 0.08662 0.1217 0.47% MultiFab::Saxpy() 91 0.07078 0.08435 0.1116 0.43% MLMG::ResNormInf() 146 0.02152 0.07464 0.1061 0.41% MultiFab::Xpay() 923 0.05309 0.07018 0.1047 0.41% Castro::construct_old_gravity_source() 20 0.04525 0.06632 0.1024 0.40% Amr::writePlotFile() 2 0.02807 0.07838 0.101 0.39% MultiFab::contains_nan() 37 0.0479 0.06671 0.09961 0.39% MultiFab::max() 8 0.0009811 0.04416 0.08217 0.32% MLMG::addInterpCorrection() 768 0.0402 0.05018 0.07117 0.28% StateData::FillBoundary(geom) 143 0 0.02338 0.06355 0.25% amrex::average_down 768 0.02625 0.03491 0.05428 0.21% TagBoxArray::collate() 2 8.788e-05 0.004058 0.05368 0.21% VisMF::Write(FabArray) 4 0.01288 0.01726 0.04201 0.16% Castro::problem_post_timestep() 2 0.0005884 0.02447 0.03743 0.14% Gravity::get_old_grav_vector() 20 0.01933 0.02532 0.03714 0.14% Castro::construct_new_source() 102 0.016 0.0229 0.03472 0.13% Castro::retry_advance_ctu() 11 0.01984 0.0249 0.03412 0.13% Gravity::get_new_grav_vector() 18 0.0175 0.02289 0.03281 0.13% Castro::create_source_corrector() 17 0.01869 0.02349 0.03175 0.12% MLCellLinOp::defineAuxData() 19 0.01525 0.01974 0.02954 0.11% MLMG::oneIter() 128 0.0149 0.01992 0.02905 0.11% FabArray::setDomainBndry() 301 0.001234 0.009175 0.02501 0.10% MLCGSolver::bicgstab 8 0 0.001553 0.02484 0.10% FabArray::FillBoundary() 6388 0.01377 0.0173 0.0241 0.09% MLCellLinOp::smooth() 3072 0.0116 0.0149 0.02361 0.09% Castro::evaluate_source_change() 196 0.01309 0.01641 0.02304 0.09% MLMG::solve() 18 0.00103 0.01133 0.02154 0.08% FillPatchIterator::Initialize 301 0.01696 0.01804 0.02126 0.08% Castro::construct_old_source() 120 0.01019 0.01424 0.02099 0.08% Castro::volWgtSum() 66 0.01017 0.01311 0.01935 0.07% Amr::writeSmallPlotFile() 2 0.01 0.01612 0.01715 0.07% FabArrayBase::getFB() 6388 0.01005 0.01236 0.0163 0.06% MLLinOp::makeSubCommunicator() 19 0.0003824 0.01033 0.01488 0.06% Castro::post_init() 1 0.0004599 0.009434 0.01404 0.05% Castro::post_timestep() 1 0.0003038 0.009341 0.01387 0.05% FabArray::mult() 59 0.007175 0.008906 0.013 0.05% Gravity::actual_solve_with_mlmg() 19 0.007147 0.008552 0.01175 0.05% MLCellLinOp::compGrad() 18 0.00592 0.007682 0.01171 0.05% FabArrayBase::getCPC() 2005 0.007012 0.008523 0.01109 0.04% MLMG::prepareForSolve() 18 0.006891 0.008254 0.01095 0.04% MultiFab::LinComb() 203 0 0.0006382 0.01021 0.04% Gravity::solve_for_phi() 17 0.002538 0.004459 0.01005 0.04% FabArrayBase::CPC::define() 576 0.002226 0.004095 0.009037 0.03% Castro::volInBoundary() 14 0.003653 0.00488 0.007256 0.03% FabArray::ParallelCopy() 1583 0.003908 0.004994 0.00724 0.03% MLCGSolver::ParallelAllReduce 189 0 0.0004406 0.007049 0.03% MultiFab::Dot() 142 0 0.0003951 0.006322 0.02% CGSolver::sxay() 203 0 0.0003834 0.006134 0.02% BndryData::define() 19 0.003186 0.004057 0.005628 0.02% MLCellLinOp::apply() 994 0.002158 0.002891 0.00546 0.02% Castro::initialize_do_advance() 20 0.004351 0.004525 0.004934 0.02% Castro::construct_new_gravity() 17 0.002274 0.002997 0.004487 0.02% MultiFab::Copy() 18 0.002051 0.002579 0.003775 0.01% Castro::locWgtSum() 12 0.001876 0.002376 0.003463 0.01% MLCellLinOp::prepareForSolve() 19 0.001782 0.002336 0.003455 0.01% MLCellLinOp::correctionResidual() 776 0.001915 0.002399 0.003201 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.001297 0.002001 0.003188 0.01% MLPoisson::normalize() 79 0 0.0001949 0.003118 0.01% MLLinOp::defineGrids() 19 0.001821 0.002223 0.002975 0.01% MLCellLinOp::defineBC() 19 0.001687 0.002011 0.002617 0.01% MLMG:computeResOfCorrection() 768 0.001494 0.001871 0.002611 0.01% MLMG::mgVcycle() 128 0.001647 0.001976 0.002601 0.01% MLMG::MLRhsNormInf() 18 0.001151 0.001514 0.002248 0.01% StateDataPhysBCFunct::() 301 0.00126 0.001538 0.002126 0.01% FabArrayBase::FB::FB() 75 0.0002864 0.0008238 0.001879 0.01% StateData::define() 6 0.0012 0.001436 0.001732 0.01% Castro::volProductSum() 4 0.0007868 0.001044 0.001561 0.01% Gravity::test_composite_phi() 1 0.0006954 0.001055 0.001451 0.01% FillPatchSingleLevel 301 0.000878 0.001048 0.0014 0.01% main() 1 0.001146 0.001236 0.001338 0.01% MLMG::mgVcycle_down::1 128 0.0007282 0.0009231 0.001335 0.01% MLMG::actualBottomSolve() 128 0.0001516 0.0002524 0.001331 0.01% MLMG::mgVcycle_down::2 128 0.0006805 0.0009116 0.001329 0.01% MLMG::mgVcycle_down::3 128 0.0006264 0.0008232 0.001296 0.01% MLMG::mgVcycle_down::4 128 0.0006144 0.0007869 0.001287 0.00% MLCellLinOp::setLevelBC() 19 0.000664 0.0008406 0.001265 0.00% MLMG::mgVcycle_down::0 128 0.0007408 0.0009277 0.00125 0.00% Gravity::applyMetricTerms() 2 0.0004822 0.0006983 0.001163 0.00% MLMG::mgVcycle_down::5 128 0.0006027 0.000755 0.001152 0.00% Castro::subcycle_advance_ctu() 1 0.0003108 0.0004248 0.00106 0.00% MLMG::mgVcycle_up::0 128 0.0005257 0.0006768 0.0009805 0.00% MLMG::mgVcycle_up::1 128 0.0005297 0.0006781 0.0009798 0.00% AmrMesh::MakeNewGrids() 2 0.0004944 0.0006403 0.0009589 0.00% MLMG::mgVcycle_up::2 128 0.0004365 0.0005896 0.0009221 0.00% Castro::enforce_consistent_e() 1 0.0003784 0.0005599 0.0009072 0.00% MLMG::mgVcycle_up::3 128 0.0004335 0.0005525 0.0008855 0.00% MLCellLinOp::solutionResidual() 147 0.0004453 0.0005532 0.0008426 0.00% Castro::swap_state_time_levels() 8 0.0005113 0.0006834 0.0008071 0.00% Castro::clean_state() 79 0.00045 0.0005752 0.000783 0.00% MLMG::mgVcycle_up::5 128 0.000456 0.0005451 0.0007358 0.00% MLMG::mgVcycle_up::4 128 0.0004315 0.0005227 0.0007144 0.00% Gravity::update_max_rhs() 2 0.0003994 0.0004935 0.0006584 0.00% makeSFC 95 0.0004129 0.0004856 0.0006471 0.00% MLLinOp::makeAgglomeratedDMap 19 0.0003791 0.0004401 0.0005727 0.00% MLMG::computeResidual() 128 0.0002722 0.0003548 0.0005644 0.00% MLMG::compResidual() 1 0.0002548 0.0003185 0.0004551 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002672 0.0003141 0.000399 0.00% DistributionMapping::Distribute() 95 0.0002265 0.0002866 0.0003877 0.00% MLMG::mgVcycle_bottom 128 0.0002347 0.0002829 0.0003751 0.00% DistributionMapping::LeastUsedCPUs() 1 2.939e-05 0.0002582 0.0003661 0.00% Castro::apply_tagging_restrictions() 2 0.0001428 0.0001998 0.0003108 0.00% OwnerMask() 2 0.000151 0.0001964 0.0002826 0.00% Castro::finalize_do_advance() 17 0.0001365 0.0001829 0.0002711 0.00% Castro::apply_problem_tags() 2 0.0001657 0.000203 0.00027 0.00% Gravity::actual_multilevel_solve() 1 0.0001389 0.0001759 0.0002522 0.00% MLLinOp::define() 19 0.0001614 0.000189 0.0002261 0.00% Castro::check_for_nan() 37 0.0001115 0.0001483 0.0002094 0.00% TagBoxArray::mapPRD 2 0.0001138 0.0001376 0.000187 0.00% MLPoisson::define() 19 0.0001235 0.0001497 0.0001838 0.00% Castro::enforce_speed_limit() 79 0.0001045 0.0001345 0.0001805 0.00% Castro::construct_old_gravity() 20 9.957e-05 0.000121 0.0001702 0.00% Amr::coarseTimeStep() 1 0.0001188 0.0001612 0.0001665 0.00% Castro::apply_source_to_state() 37 8.968e-05 0.0001175 0.0001629 0.00% Castro::initialize_advance() 1 9.89e-05 0.0001122 0.0001318 0.00% Amr::grid_places() 2 7.965e-05 9.796e-05 0.0001234 0.00% MLPoisson::prepareForSolve() 19 7.678e-05 8.953e-05 0.0001168 0.00% Amr::initSubcycle() 1 5.707e-05 7.578e-05 0.0001088 0.00% MultiFab::Swap() 5 7.335e-05 8.407e-05 0.000108 0.00% knapsack() 1 6.104e-05 7.411e-05 0.0001026 0.00% Gravity::multilevel_solve_for_new_phi() 1 4.972e-05 5.885e-05 7.371e-05 0.00% MLMG::getGradSolution() 18 4.513e-05 5.409e-05 6.928e-05 0.00% MLMG::MLResNormInf() 18 3.606e-05 4.617e-05 6.922e-05 0.00% MLMG::buildFineMask() 18 4.331e-05 5.307e-05 6.717e-05 0.00% MLMG::computeMLResidual() 18 4.466e-05 5.419e-05 6.482e-05 0.00% Amr::FinalizeInit() 1 9.19e-06 1.461e-05 5.989e-05 0.00% Gravity::swapTimeLevels() 8 3.515e-05 4.202e-05 5.527e-05 0.00% Castro::estTimeStep() 19 2.893e-05 3.736e-05 5.494e-05 0.00% Amr::regrid() 1 1.189e-05 1.663e-05 5.12e-05 0.00% Castro::construct_old_rotation_source() 20 3.075e-05 3.584e-05 4.99e-05 0.00% AmrLevel::AmrLevel(dm) 1 2.959e-05 3.621e-05 4.711e-05 0.00% Castro::construct_new_rotation_source() 17 2.288e-05 2.87e-05 4.304e-05 0.00% DistributionMapping::KnapSackDoIt() 1 1.744e-05 2.111e-05 3.472e-05 0.00% Castro::finalize_advance() 1 2.041e-05 2.486e-05 3.077e-05 0.00% Amr::timeStep() 1 1.356e-05 1.682e-05 2.84e-05 0.00% Amr::init() 1 8.317e-06 1.187e-05 2.235e-05 0.00% Amr::InitializeInit() 1 7.297e-06 1.17e-05 1.855e-05 0.00% Castro::computeInitialDt() 2 7.443e-06 8.971e-06 1.833e-05 0.00% Castro::post_regrid() 2 1.135e-05 1.379e-05 1.83e-05 0.00% Amr::bldFineLevels() 1 7.525e-06 9.588e-06 1.698e-05 0.00% BoxList::parallelComplementIn() 4 7.607e-06 9.448e-06 1.265e-05 0.00% BoxList::complementIn 4 6.346e-06 7.956e-06 1.082e-05 0.00% Castro::errorEst() 2 5.967e-06 7.559e-06 1.057e-05 0.00% Castro::advance() 1 6.582e-06 7.822e-06 9.553e-06 0.00% Castro::gwstrain() 2 4.494e-06 6.002e-06 8.468e-06 0.00% Gravity::set_mass_offset() 2 4.898e-06 5.971e-06 7.818e-06 0.00% Amr::initialInit() 1 3.066e-06 3.747e-06 4.871e-06 0.00% Castro::FluxRegFineAdd() 1 1.654e-06 2.116e-06 2.948e-06 0.00% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 25.84 25.84 25.84 100.00% Amr::coarseTimeStep() 1 21.85 21.85 21.85 84.55% Amr::timeStep() 1 21.75 21.75 21.75 84.18% Castro::advance() 1 21.51 21.54 21.57 83.49% Castro::subcycle_advance_ctu() 1 21.46 21.47 21.48 83.15% Castro::do_advance_ctu() 20 21.41 21.42 21.44 83.00% Castro::construct_ctu_hydro_source() 20 4.699 7.395 10.19 39.43% Castro::clean_state() 79 3.774 5.426 7.4 28.64% MultiFab::min() 23 0.003435 2.86 5.586 21.62% Castro::computeTemp() 80 2.514 3.632 4.878 18.88% FabArray::ParallelCopy() 1583 0.6642 2.832 4.567 17.68% FabArray::ParallelCopy_finish() 1583 0.18 2.489 4.304 16.66% Castro::reset_internal_energy() 80 1.525 2.297 3.67 14.20% FillPatchIterator::Initialize 301 0.6121 2.189 3.651 14.13% FillPatchSingleLevel 301 0.568 2.162 3.63 14.05% Gravity::solve_phi_with_mlmg() 19 2.72 3.158 3.503 13.56% Castro::construct_new_gravity() 17 2.582 3.015 3.357 12.99% Gravity::solve_for_phi() 17 2.518 2.956 3.298 12.77% Amr::InitAmr() 1 2.525 2.526 2.526 9.78% Castro::initialize_do_advance() 20 2.185 2.302 2.46 9.52% Gravity::actual_solve_with_mlmg() 19 2.312 2.324 2.334 9.03% MLMG::solve() 18 2.239 2.254 2.262 8.75% Castro::do_old_sources() 20 1.075 1.553 2.096 8.11% MLMG::oneIter() 128 1.936 1.961 2.001 7.74% MLMG::mgVcycle() 128 1.921 1.941 1.972 7.63% Castro::do_new_sources() 17 0.9786 1.421 1.917 7.42% Castro::evaluate_source_change() 196 0.9636 1.257 1.841 7.12% Castro::enforce_min_density() 79 0.7995 1.171 1.783 6.90% Amr::init() 1 1.465 1.465 1.465 5.67% Amr::initialInit() 1 1.368 1.368 1.368 5.30% MLCellLinOp::smooth() 3072 1.018 1.116 1.35 5.23% Gravity::fill_multipole_BCs() 19 0.3969 0.8329 1.18 4.57% Castro::expand_state() 20 0.2846 0.7214 1.146 4.44% Castro::normalize_species() 79 0.4854 0.7119 1.084 4.20% MLCellLinOp::applyBC() 7156 0.482 0.8449 1.073 4.15% Castro::create_source_corrector() 17 0.08969 0.617 1.071 4.14% MultiFab::AddProduct() 3439 0.5568 0.7284 1.043 4.04% FabArray::FillBoundary() 6388 0.3823 0.7629 1.012 3.92% Amr::FinalizeInit() 1 0.4816 0.6858 0.9155 3.54% Amr::InitializeInit() 1 0.4526 0.6824 0.8863 3.43% Amr::defBaseLevel() 1 0.4526 0.6824 0.8863 3.43% FillBoundary_finish() 6388 0.1407 0.593 0.8608 3.33% MLPoisson::Fsmooth() 6144 0.2821 0.4344 0.7697 2.98% MLMG::mgVcycle_up::1 128 0.1664 0.4295 0.7633 2.95% MLMG::addInterpCorrection() 768 0.1916 0.5013 0.6624 2.56% FabArray::setVal() 5021 0.3626 0.4492 0.6363 2.46% MultiFab::Subtract() 159 0.3358 0.4414 0.6315 2.44% MLMG::mgVcycle_down::0 128 0.433 0.467 0.5279 2.04% FabArray::ParallelCopy_nowait() 1583 0.2355 0.3382 0.5087 1.97% MLMG::mgVcycle_up::2 128 0.01001 0.1992 0.4923 1.91% Amr::bldFineLevels() 1 0.08689 0.2695 0.481 1.86% Castro::initData() 1 0.07884 0.1843 0.4468 1.73% Amr::grid_places() 2 0.004144 0.2082 0.438 1.70% AmrMesh::MakeNewGrids() 2 0.004021 0.2082 0.438 1.69% TagBoxArray::mapPRD 2 0.002112 0.2028 0.4364 1.69% MultiFab::sum() 3439 0.1874 0.256 0.4031 1.56% Castro::post_init() 1 0.3113 0.3539 0.3916 1.52% MLMG::mgVcycle_up::0 128 0.359 0.3715 0.3904 1.51% MLCellLinOp::apply() 994 0.3045 0.3243 0.3539 1.37% Castro::construct_old_gravity() 20 0.07792 0.238 0.3468 1.34% Gravity::get_old_grav_vector() 20 0.07778 0.2379 0.3467 1.34% amrex::average_down 768 0.03477 0.07741 0.3381 1.31% MLMG::mgVcycle_down::5 128 0.01201 0.03476 0.3148 1.22% MLMG::mgVcycle_down::1 128 0.2161 0.2641 0.3068 1.19% MLMG:computeResOfCorrection() 768 0.1812 0.2127 0.2859 1.11% MLCellLinOp::correctionResidual() 776 0.1797 0.2112 0.2833 1.10% MLMG::mgVcycle_up::3 128 0.009937 0.03238 0.2749 1.06% Castro::derive() 156 0.1431 0.1871 0.2672 1.03% Castro::construct_new_source() 102 0.1269 0.1764 0.2604 1.01% FillBoundary_nowait() 6388 0.09481 0.1525 0.2523 0.98% MLPoisson::Fapply() 994 0.09813 0.1465 0.2399 0.93% Castro::post_timestep() 1 0.1791 0.2096 0.237 0.92% Castro::problem_post_timestep() 2 0.1344 0.1866 0.2329 0.90% MLCellLinOp::solutionResidual() 147 0.1751 0.1847 0.1919 0.74% Gravity::multilevel_solve_for_new_phi() 1 0.1811 0.1813 0.1817 0.70% Gravity::actual_multilevel_solve() 1 0.1811 0.1812 0.1816 0.70% Amr::writePlotFile() 2 0.1712 0.1713 0.1713 0.66% Castro::post_regrid() 2 0.08587 0.1236 0.1661 0.64% Castro::construct_new_gravity_source() 17 0.07048 0.1049 0.1625 0.63% Castro::construct_old_source() 120 0.07922 0.1092 0.1623 0.63% MLMG::computeResidual() 128 0.1387 0.1459 0.1604 0.62% Castro::wd_update() 2 0.0548 0.107 0.1533 0.59% MLMG::mgVcycle_down::2 128 0.01298 0.05924 0.1193 0.46% MLMG::mgVcycle_up::5 128 0.01001 0.01862 0.1181 0.46% MultiFab::Saxpy() 91 0.07078 0.08435 0.1116 0.43% MLMG::ResNormInf() 146 0.02152 0.07464 0.1061 0.41% MLMG::mgVcycle_down::3 128 0.01222 0.02567 0.105 0.41% MultiFab::Xpay() 923 0.05309 0.07018 0.1047 0.41% Castro::construct_old_gravity_source() 20 0.04525 0.06632 0.1024 0.40% Castro::volWgtSum() 66 0.05381 0.07041 0.1024 0.40% Castro::check_for_nan() 37 0.04802 0.06686 0.0998 0.39% MultiFab::contains_nan() 37 0.0479 0.06671 0.09961 0.39% Castro::initialize_advance() 1 0.05016 0.06999 0.09421 0.36% MLMG::mgVcycle_bottom 128 0.0003923 0.006087 0.09046 0.35% MLMG::actualBottomSolve() 128 0.0001516 0.005805 0.09017 0.35% MLCGSolver::bicgstab 8 0 0.005519 0.0883 0.34% Amr::regrid() 1 0.04451 0.06272 0.08436 0.33% Gravity::update_max_rhs() 2 0.002782 0.04533 0.08305 0.32% MultiFab::max() 8 0.0009811 0.04416 0.08217 0.32% Gravity::get_new_grav_vector() 18 0.05181 0.05955 0.06742 0.26% MLPoisson::define() 19 0.03431 0.04717 0.06625 0.26% StateDataPhysBCFunct::() 301 0.001773 0.02492 0.06554 0.25% StateData::FillBoundary(geom) 143 0 0.02338 0.06355 0.25% Castro::retry_advance_ctu() 11 0.03754 0.04657 0.06243 0.24% MLMG::mgVcycle_down::4 128 0.01211 0.01747 0.06132 0.24% TagBoxArray::collate() 2 8.788e-05 0.004058 0.05368 0.21% MLMG::computeMLResidual() 18 0.02032 0.03758 0.04706 0.18% Castro::apply_source_to_state() 37 0.02887 0.03444 0.04556 0.18% MLMG::mgVcycle_up::4 128 0.009973 0.0139 0.04438 0.17% VisMF::Write(FabArray) 4 0.01288 0.01726 0.04201 0.16% MLCellLinOp::defineAuxData() 19 0.01923 0.02565 0.03963 0.15% Castro::volInBoundary() 14 0.02063 0.0272 0.03938 0.15% MLMG::prepareForSolve() 18 0.01763 0.02206 0.03101 0.12% Gravity::test_composite_phi() 1 0.02675 0.02714 0.02733 0.11% FabArray::setDomainBndry() 301 0.001234 0.009175 0.02501 0.10% Amr::writeSmallPlotFile() 2 0.0208 0.0208 0.02081 0.08% FabArrayBase::getCPC() 2005 0.009238 0.01262 0.02013 0.08% MLLinOp::define() 19 0.004533 0.01395 0.01865 0.07% MLLinOp::defineGrids() 19 0.004335 0.01376 0.01843 0.07% FabArrayBase::getFB() 6388 0.01034 0.01318 0.01782 0.07% CGSolver::sxay() 203 0 0.001022 0.01635 0.06% MLMG::getGradSolution() 18 0.008544 0.01058 0.01553 0.06% MLCellLinOp::compGrad() 18 0.008498 0.01053 0.01547 0.06% MLLinOp::makeSubCommunicator() 19 0.0003824 0.01033 0.01488 0.06% FabArray::mult() 59 0.007175 0.008906 0.013 0.05% MLCGSolver::ParallelAllReduce 189 0 0.0006547 0.01047 0.04% MLCellLinOp::defineBC() 19 0.005855 0.007421 0.01024 0.04% MultiFab::LinComb() 203 0 0.0006382 0.01021 0.04% Castro::locWgtSum() 12 0.005178 0.006639 0.009392 0.04% FabArrayBase::CPC::define() 576 0.002226 0.004095 0.009037 0.03% BndryData::define() 19 0.00411 0.00541 0.007685 0.03% MultiFab::Dot() 142 0 0.0003951 0.006322 0.02% Castro::volProductSum() 4 0.002989 0.003861 0.005536 0.02% MultiFab::Copy() 18 0.002051 0.002579 0.003775 0.01% MLPoisson::prepareForSolve() 19 0.001865 0.002425 0.003553 0.01% MLCellLinOp::prepareForSolve() 19 0.001782 0.002336 0.003455 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.001297 0.002001 0.003188 0.01% MLPoisson::normalize() 79 0 0.0001949 0.003118 0.01% MLMG::MLResNormInf() 18 0.001214 0.001583 0.002348 0.01% MLMG::compResidual() 1 0.001665 0.002122 0.002341 0.01% MLMG::MLRhsNormInf() 18 0.001151 0.001514 0.002248 0.01% Gravity::swapTimeLevels() 8 0.001214 0.001508 0.002159 0.01% FabArrayBase::FB::FB() 75 0.0002864 0.0008238 0.001879 0.01% AmrLevel::AmrLevel(dm) 1 0.001231 0.001472 0.001777 0.01% StateData::define() 6 0.0012 0.001436 0.001732 0.01% MLLinOp::makeAgglomeratedDMap 19 0.001018 0.001212 0.001607 0.01% MLCellLinOp::setLevelBC() 19 0.000664 0.0008406 0.001265 0.00% Gravity::applyMetricTerms() 2 0.0004822 0.0006983 0.001163 0.00% makeSFC 95 0.0006394 0.0007722 0.001035 0.00% Castro::enforce_consistent_e() 1 0.0003784 0.0005599 0.0009072 0.00% Castro::swap_state_time_levels() 8 0.0005113 0.0006834 0.0008071 0.00% Castro::errorEst() 2 0.0003255 0.0004104 0.0005902 0.00% DistributionMapping::KnapSackDoIt() 1 0.0001667 0.0003534 0.0004469 0.00% DistributionMapping::Distribute() 95 0.0002265 0.0002866 0.0003877 0.00% DistributionMapping::LeastUsedCPUs() 1 2.939e-05 0.0002582 0.0003661 0.00% Castro::apply_tagging_restrictions() 2 0.0001428 0.0001998 0.0003108 0.00% OwnerMask() 2 0.000151 0.0001964 0.0002826 0.00% Castro::finalize_do_advance() 17 0.0001365 0.0001829 0.0002711 0.00% Castro::apply_problem_tags() 2 0.0001657 0.000203 0.00027 0.00% Castro::enforce_speed_limit() 79 0.0001045 0.0001345 0.0001805 0.00% Amr::initSubcycle() 1 5.707e-05 7.578e-05 0.0001088 0.00% MultiFab::Swap() 5 7.335e-05 8.407e-05 0.000108 0.00% knapsack() 1 6.104e-05 7.411e-05 0.0001026 0.00% MLMG::buildFineMask() 18 4.331e-05 5.307e-05 6.717e-05 0.00% Castro::estTimeStep() 19 2.893e-05 3.736e-05 5.494e-05 0.00% Castro::construct_old_rotation_source() 20 3.075e-05 3.584e-05 4.99e-05 0.00% Castro::construct_new_rotation_source() 17 2.288e-05 2.87e-05 4.304e-05 0.00% Castro::finalize_advance() 1 2.223e-05 2.698e-05 3.372e-05 0.00% BoxList::parallelComplementIn() 4 1.395e-05 1.74e-05 2.326e-05 0.00% Castro::computeInitialDt() 2 1.09e-05 1.303e-05 2.273e-05 0.00% BoxList::complementIn 4 6.346e-06 7.956e-06 1.082e-05 0.00% Castro::gwstrain() 2 4.494e-06 6.002e-06 8.468e-06 0.00% Gravity::set_mass_offset() 2 4.898e-06 5.971e-06 7.818e-06 0.00% Castro::FluxRegFineAdd() 1 1.654e-06 2.116e-06 2.948e-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]::integrator.retry_burn(nvals = 1) :: [1] [TOP]::mg.v(nvals = 1) :: [0] AMReX (22.03-34-g8b691cbf8ddb) finalized