MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (23.10-4-ge470d3350ed3) initialized Starting run at 07:20:14 UTC on 2023-10-03. Successfully read inputs file ... Castro git describe: 23.10-4-g104a57d85 AMReX git describe: 23.10-4-ge470d3350 Microphysics git describe: 23.10-1-g4803fc8b 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 Initializing the data at level 0 Done initializing the level 0 data grid_places() time: 0.136658915 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001004583 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.011193779 Gravity::multilevel_solve_for_new_phi() time = 0.098095698 Gravity::fill_multipole_BCs() time = 0.011533658 ... 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 = 316439254.4 TIME= 0 CENTER OF MASS X-VEL = 0 TIME= 0 CENTER OF MASS Y-LOC = -4.854673355e-07 TIME= 0 CENTER OF MASS Y-VEL = 0 TIME= 0 CENTER OF MASS Z-LOC = 0 TIME= 0 CENTER OF MASS Z-VEL = 0 TIME= 0 MAXIMUM TEMPERATURE = 10000000 TIME= 0 MAXIMUM DENSITY = 4476474 TIME= 0 MAXIMUM T_S / T_E = 0 INITIAL GRIDS Level 0 32 grids 32768 cells 100 % of domain smallest grid: 32 x 32 biggest grid: 32 x 32 PLOTFILE: file = wdmerger_retry-SDC_plt00000 Write plotfile time = 0.050035302 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.007525995 seconds [Level 0 step 1] ADVANCE at time 0 with dt = 0.25  Beginning subcycle 1 starting at time 0 with dt = 0.25 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. Castro::construct_old_gravity() time = 0.004367444 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -5.111789233e+47 ener added: 3.959219145e+48 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.174137392e+41 ymom added: 6.498822529e+28 zmom added: 0 eint added: 0 ener added: 6.931398274e+47 Castro::do_old_sources() time = 0.064937061 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.767071355 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.767071355 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071149 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477042 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035501 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035501 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477042 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.7694089 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071149 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. Castro::construct_old_gravity() time = 0.016685145 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -5.024842118e+47 ener added: 6.151153052e+47 Contributions to the state from old-time sources: mass added: 0 xmom added: -5.870686959e+40 ymom added: 3.249411264e+28 zmom added: 0 eint added: 0 ener added: 3.465699137e+47 Castro::do_old_sources() time = 0.070168013 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.125 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.0625   Beginning subcycle 1 starting at time 0 with dt = 0.0625 Estimated number of subcycles remaining: 4 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. Castro::construct_old_gravity() time = 0.002616341 on level 0 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.060234491 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.1060744788 at index 20, 78, 0 Invalid density = -0.1060744788 at index 20, 177, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.0625 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.03125   Beginning subcycle 1 starting at time 0 with dt = 0.03125 Estimated number of subcycles remaining: 8 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. Castro::construct_old_gravity() time = 0.003094403 on level 0 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.059464106 on level 0 ... Entering construct_ctu_hydro_source() on level 0 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.319111915 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014824957 Gravity::solve_for_phi() time = 0.084860308 on level 0 Castro::construct_new_gravity() time = 0.086980074 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.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.060004791 on level 0 Ending SDC iteration 1 of 3. Beginning SDC iteration 2 of 3. Castro::construct_old_gravity() time = 0.002811527 on level 0 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.059091959 on level 0 ... Entering construct_ctu_hydro_source() on level 0 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.315084868 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014299608 Gravity::solve_for_phi() time = 0.083563079 on level 0 Castro::construct_new_gravity() time = 0.085672907 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.838148647e+44 ener added: 3.819761762e+42 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.742583287e+35 ymom added: -3.588335034e+27 zmom added: 0 eint added: 0 ener added: 4.570264853e+43 Castro::do_new_sources() time = 0.059535934 on level 0 Ending SDC iteration 2 of 3. Beginning SDC iteration 3 of 3. Castro::construct_old_gravity() time = 0.002712228 on level 0 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.05923001 on level 0 ... Entering construct_ctu_hydro_source() on level 0 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.322430546 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014385739 Gravity::solve_for_phi() time = 0.084019649 on level 0 Castro::construct_new_gravity() time = 0.086116752 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.838245244e+44 ener added: 3.819942812e+42 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.742573524e+35 ymom added: -3.587065898e+27 zmom added: 0 eint added: 0 ener added: 4.569279535e+43 Castro::do_new_sources() time = 0.058913949 on level 0 Ending SDC iteration 3 of 3.  Subcycle completed  Beginning subcycle 2 starting at time 0.03125 with dt = 0.03125 Estimated number of subcycles remaining: 7 Beginning SDC iteration 1 of 2. Castro::construct_old_gravity() time = 0.002559512 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.985021785e+46 ener added: 5.435221686e+40 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467706591e+40 ymom added: 9.489312105e+26 zmom added: 0 eint added: 0 ener added: 8.721813579e+46 Castro::do_old_sources() time = 0.05913356 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.985598081e+46 ener added: 2.745852676e+42 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.323290334 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014855673 Gravity::solve_for_phi() time = 0.084938441 on level 0 Castro::construct_new_gravity() time = 0.087043477 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.274072206e+44 ener added: 4.891864863e+43 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.763054894e+35 ymom added: -9.059860097e+24 zmom added: 0 eint added: 0 ener added: 5.968608016e+43 Castro::do_new_sources() time = 0.059486435 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.002718299 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.985021785e+46 ener added: 5.435221686e+40 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467706591e+40 ymom added: 9.489312105e+26 zmom added: 0 eint added: 0 ener added: 8.721813579e+46 Castro::do_old_sources() time = 0.059130314 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.986297782e+46 ener added: 2.752525031e+42 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.327738661 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014911842 Gravity::solve_for_phi() time = 0.084074517 on level 0 Castro::construct_new_gravity() time = 0.086174211 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.347573489e+44 ener added: 4.888153391e+43 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.756949739e+35 ymom added: -3.751920171e+24 zmom added: 0 eint added: 0 ener added: 5.2374099e+43 Castro::do_new_sources() time = 0.059663684 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 3 starting at time 0.0625 with dt = 0.03125 Estimated number of subcycles remaining: 6 Beginning SDC iteration 1 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.881810984e+29 ener added: 1.881810984e+29 Castro::construct_old_gravity() time = 0.003201653 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.984353173e+46 ener added: 6.990361455e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46774173e+40 ymom added: 9.41450982e+26 zmom added: 0 eint added: 0 ener added: 8.780735352e+46 Castro::do_old_sources() time = 0.059133176 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.98733239e+46 ener added: 2.180535114e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.324222145 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014902858 Gravity::solve_for_phi() time = 0.0855276 on level 0 Castro::construct_new_gravity() time = 0.087621594 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.058042977e+43 ener added: 1.184398423e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.782398841e+35 ymom added: -3.291489439e+24 zmom added: 0 eint added: 0 ener added: 5.968491079e+43 Castro::do_new_sources() time = 0.059867832 on level 0 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.881810984e+29 ener added: 1.881810984e+29 Castro::construct_old_gravity() time = 0.003048463 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.984353173e+46 ener added: 6.990361455e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46774173e+40 ymom added: 9.41450982e+26 zmom added: 0 eint added: 0 ener added: 8.780735352e+46 Castro::do_old_sources() time = 0.059245086 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.988039066e+46 ener added: 2.185480971e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.329257318 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014935943 Gravity::solve_for_phi() time = 0.084962544 on level 0 Castro::construct_new_gravity() time = 0.087054185 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.800848844e+43 ener added: 1.183591058e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.776272322e+35 ymom added: 1.291567233e+24 zmom added: 0 eint added: 0 ener added: 5.233809214e+43 Castro::do_new_sources() time = 0.060293645 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 4 starting at time 0.09375 with dt = 0.03125 Estimated number of subcycles remaining: 5 Beginning SDC iteration 1 of 2. Castro::construct_old_gravity() time = 0.00403044 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.983917922e+46 ener added: 1.154265495e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467777256e+40 ymom added: 9.44246623e+26 zmom added: 0 eint added: 0 ener added: 8.841018465e+46 Castro::do_old_sources() time = 0.059706847 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.98498831e+46 ener added: 2.292089986e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.32672384 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015652073 Gravity::solve_for_phi() time = 0.085104578 on level 0 Castro::construct_new_gravity() time = 0.087198411 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -4.614759295e+43 ener added: 1.488589688e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.803591267e+35 ymom added: 6.776595903e+23 zmom added: 0 eint added: 0 ener added: 4.637081169e+43 Castro::do_new_sources() time = 0.059841719 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.003085222 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.983917922e+46 ener added: 1.154265495e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467777256e+40 ymom added: 9.44246623e+26 zmom added: 0 eint added: 0 ener added: 8.841018465e+46 Castro::do_old_sources() time = 0.059050874 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.985698801e+46 ener added: 2.298024982e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.329721704 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015241916 Gravity::solve_for_phi() time = 0.084213639 on level 0 Castro::construct_new_gravity() time = 0.086313899 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -5.359798245e+43 ener added: 1.487684259e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.797448342e+35 ymom added: 5.252452121e+24 zmom added: 0 eint added: 0 ener added: 3.901081726e+43 Castro::do_new_sources() time = 0.059556759 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 5 starting at time 0.125 with dt = 0.03125 Estimated number of subcycles remaining: 4 Beginning SDC iteration 1 of 2. Castro::construct_old_gravity() time = 0.002732679 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.98356532e+46 ener added: 1.521222674e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467813205e+40 ymom added: 9.540691453e+26 zmom added: 0 eint added: 0 ener added: 8.902573838e+46 Castro::do_old_sources() time = 0.059381557 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.992363162e+46 ener added: 2.321339475e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.324974993 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015383341 Gravity::solve_for_phi() time = 0.084681629 on level 0 Castro::construct_new_gravity() time = 0.086776867 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.068921344e+43 ener added: 1.620777103e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.829776437e+35 ymom added: -4.179294337e+24 zmom added: 0 eint added: 0 ener added: 2.123241033e+43 Castro::do_new_sources() time = 0.059659406 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.003143976 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.98356532e+46 ener added: 1.521222674e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467813205e+40 ymom added: 9.540691453e+26 zmom added: 0 eint added: 0 ener added: 8.902573838e+46 Castro::do_old_sources() time = 0.059301079 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.993083003e+46 ener added: 2.326773931e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.330458966 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015516624 Gravity::solve_for_phi() time = 0.085110383 on level 0 Castro::construct_new_gravity() time = 0.087242324 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.821028634e+43 ener added: 1.619931887e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.823603759e+35 ymom added: 1.205384045e+24 zmom added: 0 eint added: 0 ener added: 1.379533843e+43 Castro::do_new_sources() time = 0.059780293 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 6 starting at time 0.15625 with dt = 0.03125 Estimated number of subcycles remaining: 3 Beginning SDC iteration 1 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 4.908356001e+29 ener added: 4.908356001e+29 Castro::construct_old_gravity() time = 0.00294969 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.983073759e+46 ener added: 2.019352467e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467849677e+40 ymom added: 9.586026171e+26 zmom added: 0 eint added: 0 ener added: 8.96490842e+46 Castro::do_old_sources() time = 0.059169238 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.000609588e+46 ener added: 2.406064512e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.325939031 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015243877 Gravity::solve_for_phi() time = 0.084371173 on level 0 Castro::construct_new_gravity() time = 0.086499808 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -8.854195608e+43 ener added: 1.71395957e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.848617495e+35 ymom added: -3.296211805e+24 zmom added: 0 eint added: 0 ener added: -1.366537222e+43 Castro::do_new_sources() time = 0.059544322 on level 0 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.908356001e+29 ener added: 4.908356001e+29 Castro::construct_old_gravity() time = 0.002840191 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.983073759e+46 ener added: 2.019352467e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467849677e+40 ymom added: 9.586026171e+26 zmom added: 0 eint added: 0 ener added: 8.96490842e+46 Castro::do_old_sources() time = 0.058933468 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.001333674e+46 ener added: 2.41251015e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.329926966 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015386611 Gravity::solve_for_phi() time = 0.086008424 on level 0 Castro::construct_new_gravity() time = 0.088125929 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -9.610031177e+43 ener added: 1.713054017e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.842414694e+35 ymom added: -1.986935698e+24 zmom added: 0 eint added: 0 ener added: -2.113374984e+43 Castro::do_new_sources() time = 0.059599912 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 7 starting at time 0.1875 with dt = 0.03125 Estimated number of subcycles remaining: 2 Beginning SDC iteration 1 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.009461182e+29 ener added: 2.009461182e+29 Castro::construct_old_gravity() time = 0.00328278 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.982914566e+46 ener added: 2.174746794e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467886525e+40 ymom added: 9.502912521e+26 zmom added: 0 eint added: 0 ener added: 9.027789929e+46 Castro::do_old_sources() time = 0.058931152 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.993854135e+46 ener added: 2.676087223e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.324421918 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015498936 Gravity::solve_for_phi() time = 0.08468034 on level 0 Castro::construct_new_gravity() time = 0.086797624 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.152849071e+44 ener added: 1.743211419e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.877263956e+35 ymom added: 2.991619167e+24 zmom added: 0 eint added: 0 ener added: -4.049549507e+43 Castro::do_new_sources() time = 0.059813462 on level 0 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: 2.009461182e+29 ener added: 2.009461182e+29 Castro::construct_old_gravity() time = 0.00339714 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.982914566e+46 ener added: 2.174746794e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467886525e+40 ymom added: 9.502912521e+26 zmom added: 0 eint added: 0 ener added: 9.027789929e+46 Castro::do_old_sources() time = 0.059197139 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.994586452e+46 ener added: 2.682600444e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.329787144 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015402598 Gravity::solve_for_phi() time = 0.084496661 on level 0 Castro::construct_new_gravity() time = 0.086618967 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.229235664e+44 ener added: 1.74231031e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.871026436e+35 ymom added: -4.264296934e+24 zmom added: 0 eint added: 0 ener added: -4.804462625e+43 Castro::do_new_sources() time = 0.059760447 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 8 starting at time 0.21875 with dt = 0.03125 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 9.20728525e+29 ener added: 9.20728525e+29 Castro::construct_old_gravity() time = 0.003212142 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.982548267e+46 ener added: 2.529322566e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467923946e+40 ymom added: 9.435666022e+26 zmom added: 0 eint added: 0 ener added: 9.091353867e+46 Castro::do_old_sources() time = 0.059244309 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.993555556e+46 ener added: 2.937230463e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.325746735 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01571946 Gravity::solve_for_phi() time = 0.084961736 on level 0 Castro::construct_new_gravity() time = 0.087064691 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -8.731399969e+43 ener added: 1.766195089e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.939762635e+35 ymom added: 4.946678891e+24 zmom added: 0 eint added: 0 ener added: -1.225922477e+43 Castro::do_new_sources() time = 0.059701741 on level 0 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: 9.20728525e+29 ener added: 9.20728525e+29 Castro::construct_old_gravity() time = 0.002745011 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.982548267e+46 ener added: 2.529322566e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467923946e+40 ymom added: 9.435666022e+26 zmom added: 0 eint added: 0 ener added: 9.091353867e+46 Castro::do_old_sources() time = 0.059177488 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.994288712e+46 ener added: 2.943682879e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.32882233 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015327234 Gravity::solve_for_phi() time = 0.084440197 on level 0 Castro::construct_new_gravity() time = 0.086568776 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -9.496809054e+43 ener added: 1.76531212e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.933525523e+35 ymom added: -4.247768651e+24 zmom added: 0 eint added: 0 ener added: -1.982551467e+43 Castro::do_new_sources() time = 0.059893411 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002808440132 [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: 4.117748323e+27 ener added: 4.117748323e+27 TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -4.340632835e+38 TIME= 0.25 YMOM = 1.160568787e+28 TIME= 0.25 ZMOM = 0 TIME= 0.25 ANG MOM X = 0 TIME= 0.25 ANG MOM Y = 0 TIME= 0.25 ANG MOM Z = -1.471393084e+38 TIME= 0.25 RHO*e = 1.325894499e+50 TIME= 0.25 RHO*K = 3.170941995e+49 TIME= 0.25 RHO*E = 1.642988698e+50 TIME= 0.25 RHO*PHI = -5.449011351e+50 TIME= 0.25 TOTAL ENERGY= -1.081516977e+50 TIME= 0.25 CENTER OF MASS X-LOC = 316424375.9 TIME= 0.25 CENTER OF MASS X-VEL = -170221.1246 TIME= 0.25 CENTER OF MASS Y-LOC = 7.282010032e-07 TIME= 0.25 CENTER OF MASS Y-VEL = 4.55125627e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 58251090.19 TIME= 0.25 MAXIMUM DENSITY = 4458791.937 TIME= 0.25 MAXIMUM T_S / T_E = 0 [STEP 1] Coarse TimeStep time: 11.72617626 [STEP 1] FAB kilobyte spread across MPI nodes: [6973 ... 6973] STEP = 1 TIME = 0.25 DT = 0.25 PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.049873356 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.00729417 seconds Ending run at 07:20:28 UTC on 2023-10-03. Run time = 14.09515149 Run time without initialization = 11.78424451 Average number of zones advanced per microsecond: 0.003 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 14.1 ... 14.1 ... 14.1 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.3309 5.3974 5.4806 38.88% Amr::InitAmr() 1 1.451 1.451 1.451 10.30% Castro::computeTemp() 80 0.8531 0.9643 1.136 8.06% amrex::Copy() 795 0.7889 0.7945 0.7994 5.67% FabArray::ParallelCopy_finish() 1429 0.04999 0.4788 0.6836 4.85% MultiFab::AddProduct() 3439 0.5497 0.5537 0.5575 3.95% Castro::normalize_species() 79 0.5342 0.5395 0.5465 3.88% Castro::reset_internal_energy(Fab) 160 0.521 0.5276 0.5364 3.81% MLPoisson::Fsmooth() 6144 0.313 0.3398 0.4069 2.89% Castro::initMFs() 1 0.3589 0.368 0.3973 2.82% FabArray::setVal() 5021 0.3446 0.35 0.3549 2.52% MultiFab::Subtract() 159 0.3441 0.3466 0.3489 2.48% Gravity::fill_multipole_BCs() 19 0.1957 0.2367 0.2752 1.95% FabArray::ParallelCopy_nowait() 1429 0.1825 0.1967 0.2132 1.51% check_for_negative_density() 20 0.04961 0.1393 0.1966 1.39% FabArray::sum() 3439 0.1904 0.1918 0.194 1.38% Castro::reset_internal_energy(MultiFab) 80 0.02243 0.055 0.1633 1.16% Castro::initData() 1 0.003301 0.02261 0.1289 0.91% FillBoundary_nowait() 6388 0.0738 0.09867 0.1263 0.90% FillBoundary_finish() 6388 0.08021 0.09348 0.1216 0.86% MLPoisson::Fapply() 994 0.09438 0.09897 0.1134 0.80% MLCellLinOp::applyBC() 7156 0.04631 0.06198 0.1035 0.73% Castro::do_new_sources() 17 0.0003695 0.01456 0.07621 0.54% FabArray::Saxpy() 192 0.06857 0.06946 0.07324 0.52% Castro::construct_new_gravity_source() 17 0.06883 0.06953 0.07073 0.50% FabArray::Xpay() 953 0.04961 0.05203 0.05926 0.42% Castro::do_old_sources() 20 0.0005145 0.01288 0.0549 0.39% MultiFab::contains_nan() 37 0.04494 0.04553 0.04594 0.33% Castro::construct_old_gravity_source() 20 0.04431 0.0446 0.04502 0.32% MLMG::addInterpCorrection() 768 0.03353 0.03601 0.04286 0.30% Castro::enforce_min_density() 79 0.01714 0.02379 0.04266 0.30% AmrLevel::derive() 18 0.04075 0.0414 0.0424 0.30% StateData::FillBoundary(geom) 143 0 0.0177 0.03399 0.24% amrex::average_down 768 0.02491 0.02702 0.03274 0.23% VisMF::Write(FabArray) 4 0.01523 0.01713 0.02868 0.20% MLCGSolver::bicgstab 8 0 0.001369 0.0219 0.16% Amr::writePlotFile() 2 0.008827 0.01827 0.02063 0.15% FabArray::FillBoundary() 6388 0.01417 0.01485 0.018 0.13% MLCellLinOp::defineAuxData() 19 0.01525 0.01591 0.01789 0.13% Castro::construct_old_gravity() 20 0.0003494 0.00349 0.0174 0.12% Castro::construct_new_source() 102 0.01534 0.01556 0.01625 0.12% Castro::evaluate_source_change() 196 0.01296 0.01343 0.01506 0.11% amrex::Add() 136 0.01439 0.01449 0.01481 0.11% MLCellLinOp::smooth() 3072 0.01178 0.01233 0.01458 0.10% FabArray::norminf() 244 0.01084 0.0112 0.01394 0.10% FillPatchIterator::Initialize 147 0.01261 0.013 0.01381 0.10% FabArrayBase::getFB() 6388 0.01019 0.01065 0.01326 0.09% MLMG::prepareForSolve() 18 0.01051 0.01084 0.0117 0.08% FabArray::setDomainBndry() 147 0.0004564 0.005913 0.01105 0.08% Castro::wd_update() 2 0.004719 0.007853 0.01075 0.08% Castro::construct_old_source() 120 0.009835 0.009904 0.01027 0.07% Amr::writeSmallPlotFile() 2 0.007346 0.009033 0.009522 0.07% FabArrayBase::getCPC() 2005 0.007311 0.007577 0.008308 0.06% MLLinOp::makeSubCommunicator() 19 0.0005529 0.004136 0.008195 0.06% Gravity::actual_solve_with_mlmg() 19 0.006848 0.007128 0.008018 0.06% Castro::volWgtSum() 46 0.007273 0.007323 0.007432 0.05% FabArray::mult() 59 0.006882 0.006998 0.007284 0.05% MLMG::solve() 18 0.001544 0.003201 0.006996 0.05% Gravity::get_old_grav_vector() 20 0.006396 0.006529 0.006741 0.05% Gravity::get_new_grav_vector() 18 0.00604 0.006215 0.00665 0.05% FabArrayBase::CPC::define() 576 0.002282 0.003456 0.006628 0.05% amrex::Dot() 142 0 0.0003823 0.006117 0.04% MLCellLinOp::compGrad() 18 0.005578 0.005716 0.005833 0.04% Castro::initialize_do_advance() 20 0.004862 0.005206 0.005646 0.04% Castro::volInBoundary() 14 0.003369 0.003631 0.005417 0.04% MLMG::ResNormInf() 146 0.003183 0.004247 0.005294 0.04% MLCellLinOp::apply() 994 0.002093 0.00236 0.005221 0.04% MLCGSolver::ParallelAllReduce 189 0 0.0002881 0.00461 0.03% FabArray::ParallelCopy() 1429 0.003789 0.003958 0.004275 0.03% FabArray::LinComb() 71 0 0.0002301 0.003681 0.03% BndryData::define() 19 0.003128 0.003204 0.003441 0.02% MLPoisson::normalize() 79 0 0.0001959 0.003135 0.02% AmrLevel::FillPatch() 147 0.002192 0.002312 0.002612 0.02% MLCellLinOp::correctionResidual() 776 0.001903 0.001997 0.002431 0.02% Castro::sum_integrated_quantities() 2 0.001609 0.002046 0.002302 0.02% MLLinOp::defineGrids() 19 0.001819 0.001978 0.002204 0.02% MLCellLinOp::prepareForSolve() 19 0.001732 0.001849 0.002156 0.02% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.001325 0.001608 0.001902 0.01% MLMG::mgVcycle() 128 0.001665 0.00174 0.001894 0.01% MLCellLinOp::defineBC() 19 0.001584 0.001641 0.001759 0.01% MLMG:computeResOfCorrection() 768 0.001469 0.001513 0.001615 0.01% StateData::define() 6 0.001258 0.001428 0.001549 0.01% FabArrayBase::FB::FB() 75 0.0002896 0.0007061 0.0015 0.01% Castro::save_data_for_retry() 3 0.001394 0.001463 0.001499 0.01% main() 1 0.001044 0.001263 0.001435 0.01% TagBoxArray::collate() 2 0.0001021 0.0006865 0.001425 0.01% Castro::buildMetrics() 1 0.0005599 0.0006302 0.00134 0.01% Gravity::solve_for_phi() 17 0.0008111 0.0009288 0.001331 0.01% Castro::construct_new_gravity() 17 0.0002948 0.0005227 0.001265 0.01% Castro::locWgtSum() 6 0.001147 0.001175 0.001213 0.01% MLMG::actualBottomSolve() 128 0.0001356 0.000204 0.001164 0.01% MultiFab::max() 2 0.000169 0.000812 0.001074 0.01% StateDataPhysBCFunct::() 147 0.0006679 0.0008453 0.001033 0.01% Castro::subcycle_advance_ctu() 1 0.0001972 0.0002758 0.000984 0.01% MLCellLinOp::setLevelBC() 19 0.0005443 0.0006954 0.0008586 0.01% FillPatchIterator::FillFromLevel0() 147 0.0007573 0.0007902 0.0008379 0.01% MLMG::mgVcycle_down::2 128 0.0006716 0.0007231 0.0008346 0.01% MLMG::mgVcycle_down::0 128 0.0006906 0.0007116 0.0007899 0.01% MLMG::mgVcycle_down::1 128 0.0006894 0.0007059 0.0007774 0.01% MLMG::mgVcycle_down::3 128 0.0006218 0.0006706 0.0007527 0.01% MLMG::mgVcycle_down::4 128 0.0006029 0.0006296 0.0007176 0.01% MLMG::mgVcycle_down::5 128 0.0006251 0.0006468 0.0007152 0.01% Castro::swap_state_time_levels() 8 0.0006278 0.0006744 0.0007107 0.01% Gravity::applyMetricTerms() 2 0.000476 0.0004926 0.0005817 0.00% Castro::clean_state() 79 0.0004525 0.0004924 0.0005805 0.00% Castro::do_advance_ctu() 20 0.0004784 0.0005136 0.0005762 0.00% MLMG::mgVcycle_up::1 128 0.0004841 0.0005068 0.0005709 0.00% MLMG::mgVcycle_up::0 128 0.0004888 0.0005048 0.0005528 0.00% AmrMesh::MakeNewGrids() 2 0.0004882 0.0005165 0.0005515 0.00% Castro::Castro() 1 0.0004605 0.0004927 0.0005469 0.00% Castro::apply_tagging_restrictions() 2 0.0004219 0.0004797 0.000541 0.00% MLMG::mgVcycle_up::2 128 0.0004192 0.0004599 0.000535 0.00% FillPatchSingleLevel 147 0.000493 0.0005055 0.0005196 0.00% MLMG::mgVcycle_up::5 128 0.0004248 0.0004359 0.0004915 0.00% MLCellLinOp::solutionResidual() 147 0.0004402 0.0004538 0.0004904 0.00% MLMG::mgVcycle_up::4 128 0.0004149 0.0004301 0.0004861 0.00% MLMG::mgVcycle_up::3 128 0.0004108 0.000429 0.0004796 0.00% makeSFC 95 0.000399 0.0004236 0.0004619 0.00% MLMG::oneIter() 128 0.0003914 0.0004042 0.0004546 0.00% Castro::enforce_consistent_e() 1 0.0003608 0.0003867 0.0004293 0.00% Castro::volProductSum() 2 0.0003953 0.0003985 0.0004016 0.00% MLLinOp::makeAgglomeratedDMap 19 0.0003408 0.0003586 0.0003917 0.00% Castro::problem_post_timestep() 2 1.085e-05 3.471e-05 0.0003721 0.00% MLMG::mgVcycle_bottom 128 0.0002364 0.0002445 0.0002769 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002335 0.0002501 0.0002751 0.00% MLMG::computeResidual() 128 0.000252 0.0002586 0.0002669 0.00% DistributionMapping::Distribute() 95 0.0002196 0.0002319 0.0002479 0.00% Castro::finalize_do_advance() 17 0.000172 0.0001973 0.0002177 0.00% Gravity::update_max_rhs() 2 0.0001673 0.000174 0.0001876 0.00% Castro::retry_advance_ctu() 11 5.581e-05 6.918e-05 0.0001819 0.00% MLMG::compResidual() 1 0.0001593 0.0001608 0.0001691 0.00% OwnerMask() 2 0.000138 0.0001533 0.0001684 0.00% Castro::initialize_advance() 1 0.0001279 0.0001442 0.000158 0.00% MLPoisson::define() 19 0.0001311 0.0001376 0.0001465 0.00% Amr::coarseTimeStep() 1 0.0001164 0.0001388 0.0001437 0.00% Castro::check_for_nan() 37 0.0001097 0.0001214 0.0001422 0.00% Amr::defBaseLevel() 1 0.0001111 0.0001223 0.0001347 0.00% Amr::grid_places() 2 8.783e-05 0.0001037 0.0001295 0.00% Gravity::test_composite_phi() 1 9.444e-05 0.000101 0.000129 0.00% Castro::enforce_speed_limit() 79 0.0001088 0.0001172 0.0001276 0.00% TagBoxArray::mapPRD 2 0.0001116 0.0001193 0.0001273 0.00% Castro::apply_source_to_state() 37 9.968e-05 0.0001071 0.0001218 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.579e-05 8.689e-05 0.0001077 0.00% Amr::FinalizeInit() 1 9.478e-06 1.949e-05 0.0001019 0.00% Castro::apply_problem_tags() 2 9.143e-05 9.38e-05 9.821e-05 0.00% MLPoisson::prepareForSolve() 19 7.862e-05 8.413e-05 9.673e-05 0.00% MLLinOp::define() 19 7.981e-05 8.506e-05 9.142e-05 0.00% Castro::create_source_corrector() 17 6.062e-05 6.97e-05 9.02e-05 0.00% MultiFab::Swap() 5 7.136e-05 7.716e-05 8.947e-05 0.00% Castro::expand_state() 20 5.876e-05 6.465e-05 7.083e-05 0.00% Amr::initSubcycle() 1 5.706e-05 6.196e-05 6.329e-05 0.00% knapsack() 1 5.829e-05 6.081e-05 6.281e-05 0.00% Gravity::swapTimeLevels() 8 3.371e-05 4.136e-05 5.404e-05 0.00% MLMG::computeMLResidual() 18 4.104e-05 4.62e-05 5.249e-05 0.00% Castro::derive() 18 4.334e-05 4.636e-05 5.074e-05 0.00% MLMG::getGradSolution() 18 4.062e-05 4.25e-05 4.63e-05 0.00% MLMG::MLResNormInf() 18 3.535e-05 3.774e-05 4.428e-05 0.00% Gravity::actual_multilevel_solve() 1 2.773e-05 3.31e-05 3.925e-05 0.00% AmrLevel::AmrLevel(dm) 1 3.431e-05 3.561e-05 3.887e-05 0.00% Amr::timeStep() 1 1.255e-05 1.69e-05 3.725e-05 0.00% Castro::construct_old_rotation_source() 20 3e-05 3.217e-05 3.679e-05 0.00% MLMG::MLRhsNormInf() 18 3.202e-05 3.294e-05 3.527e-05 0.00% DistributionMapping::LeastUsedCPUs() 1 1.584e-05 2.551e-05 3.155e-05 0.00% Amr::regrid() 1 1.148e-05 1.401e-05 3.142e-05 0.00% Castro::estTimeStep() 10 1.306e-05 1.507e-05 2.446e-05 0.00% Gravity::set_mass_offset() 2 4.406e-06 9.355e-06 2.434e-05 0.00% Castro::construct_new_rotation_source() 17 2.192e-05 2.287e-05 2.426e-05 0.00% Castro::finalize_advance() 1 1.751e-05 1.894e-05 2.418e-05 0.00% Castro::post_regrid() 2 1.15e-05 1.428e-05 2.296e-05 0.00% DistributionMapping::KnapSackDoIt() 1 1.677e-05 1.791e-05 1.985e-05 0.00% Castro::computeInitialDt() 2 7.557e-06 1.21e-05 1.941e-05 0.00% Castro::post_init() 1 1.254e-05 1.312e-05 1.432e-05 0.00% Castro::advance() 1 7.272e-06 8.636e-06 1.304e-05 0.00% Amr::InitializeInit() 1 6.579e-06 8.376e-06 1.291e-05 0.00% Amr::bldFineLevels() 1 7.556e-06 9.376e-06 1.134e-05 0.00% Amr::init() 1 7.508e-06 8.157e-06 1.014e-05 0.00% Castro::gwstrain() 2 4.256e-06 5.476e-06 9.934e-06 0.00% Castro::errorEst() 2 6.454e-06 6.885e-06 7.927e-06 0.00% BoxList::parallelComplementIn() 4 7.229e-06 7.543e-06 7.807e-06 0.00% Castro::post_timestep() 1 5.753e-06 6.427e-06 7.613e-06 0.00% BoxList::complementIn 4 6.219e-06 6.496e-06 6.814e-06 0.00% Amr::initialInit() 1 3.08e-06 3.459e-06 4.828e-06 0.00% Castro::FluxRegCrseInit 1 1.695e-06 2.047e-06 2.631e-06 0.00% Castro::FluxRegFineAdd() 1 1.476e-06 1.67e-06 2.067e-06 0.00% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.1 14.1 14.1 100.00% Amr::coarseTimeStep() 1 11.78 11.78 11.78 83.60% Amr::timeStep() 1 11.73 11.73 11.73 83.19% Castro::advance() 1 11.66 11.67 11.67 82.78% Castro::subcycle_advance_ctu() 1 11.61 11.62 11.62 82.42% Castro::do_advance_ctu() 20 11.57 11.58 11.58 82.15% Castro::construct_ctu_hydro_source() 20 6.28 6.317 6.368 45.18% Castro::clean_state() 79 3.689 3.8 3.965 28.13% Castro::computeTemp() 80 2.397 2.506 2.656 18.85% Castro::reset_internal_energy(MultiFab) 80 1.499 1.541 1.645 11.67% Gravity::solve_phi_with_mlmg() 19 1.463 1.506 1.548 10.98% Castro::construct_new_gravity() 17 1.39 1.433 1.474 10.46% Amr::InitAmr() 1 1.451 1.451 1.451 10.30% Gravity::solve_for_phi() 17 1.354 1.397 1.438 10.21% Castro::initialize_do_advance() 20 1.364 1.387 1.425 10.11% Gravity::actual_solve_with_mlmg() 19 1.263 1.266 1.271 9.01% MLMG::solve() 18 1.212 1.215 1.217 8.63% Castro::do_old_sources() 20 1.15 1.183 1.201 8.52% MLMG::oneIter() 128 1.086 1.087 1.088 7.72% MLMG::mgVcycle() 128 1.071 1.072 1.073 7.61% Castro::do_new_sources() 17 0.938 0.9814 1.015 7.20% Castro::evaluate_source_change() 196 0.95 0.9569 0.9666 6.86% FabArray::ParallelCopy() 1429 0.2721 0.6992 0.9011 6.39% Amr::init() 1 0.8589 0.8589 0.859 6.09% Castro::enforce_min_density() 79 0.7935 0.8036 0.8172 5.80% Amr::initialInit() 1 0.8005 0.8005 0.8006 5.68% amrex::Copy() 795 0.7889 0.7945 0.7994 5.67% MLCellLinOp::smooth() 3072 0.5311 0.5783 0.7026 4.98% FabArray::ParallelCopy_finish() 1429 0.04999 0.4788 0.6836 4.85% AmrLevel::FillPatch() 147 0.4674 0.5494 0.6282 4.46% Amr::InitializeInit() 1 0.4447 0.4953 0.5793 4.11% Amr::defBaseLevel() 1 0.4447 0.4953 0.5793 4.11% MultiFab::AddProduct() 3439 0.5497 0.5537 0.5575 3.95% Castro::normalize_species() 79 0.5342 0.5395 0.5465 3.88% Castro::reset_internal_energy(Fab) 160 0.521 0.5276 0.5364 3.81% FillPatchIterator::Initialize 147 0.2881 0.3686 0.4484 3.18% FillPatchIterator::FillFromLevel0() 147 0.2644 0.3497 0.4315 3.06% FillPatchSingleLevel 147 0.2636 0.3489 0.4306 3.06% MLMG::mgVcycle_up::1 128 0.08919 0.2554 0.4214 2.99% MLPoisson::Fsmooth() 6144 0.313 0.3398 0.4069 2.89% Castro::Castro() 1 0.3642 0.3733 0.4028 2.86% Castro::initMFs() 1 0.3591 0.3681 0.3975 2.82% MLMG::addInterpCorrection() 768 0.0628 0.3083 0.3756 2.66% MLCellLinOp::applyBC() 7156 0.2488 0.2803 0.3578 2.54% Amr::FinalizeInit() 1 0.2212 0.3052 0.3558 2.52% FabArray::setVal() 5021 0.3446 0.35 0.3549 2.52% MultiFab::Subtract() 159 0.3441 0.3466 0.3489 2.48% MLMG::mgVcycle_up::2 128 0.0111 0.1177 0.2823 2.00% Gravity::fill_multipole_BCs() 19 0.1979 0.2388 0.2772 1.97% MLMG::mgVcycle_down::0 128 0.2549 0.2563 0.2574 1.83% FabArray::FillBoundary() 6388 0.1898 0.2184 0.2543 1.80% FabArray::ParallelCopy_nowait() 1429 0.2019 0.2165 0.2339 1.66% Castro::expand_state() 20 0.176 0.1859 0.2118 1.50% Castro::initData() 1 0.07853 0.1203 0.2056 1.46% check_for_negative_density() 20 0.04961 0.1393 0.1966 1.39% FabArray::sum() 3439 0.1904 0.1918 0.194 1.38% MLCellLinOp::apply() 994 0.1439 0.1539 0.1923 1.36% MLMG::mgVcycle_up::0 128 0.1881 0.1888 0.19 1.35% MLMG::mgVcycle_up::3 128 0.01103 0.02449 0.1878 1.33% Amr::bldFineLevels() 1 0.04642 0.1303 0.1808 1.28% MLCellLinOp::correctionResidual() 776 0.1085 0.1178 0.1447 1.03% MLMG:computeResOfCorrection() 768 0.11 0.119 0.142 1.01% FillBoundary_nowait() 6388 0.08439 0.11 0.1381 0.98% Amr::grid_places() 2 0.003053 0.08706 0.1377 0.98% AmrMesh::MakeNewGrids() 2 0.002964 0.08695 0.1376 0.98% TagBoxArray::mapPRD 2 0.0007485 0.085 0.135 0.96% Castro::post_init() 1 0.1315 0.1321 0.1323 0.94% Castro::construct_new_source() 102 0.1228 0.124 0.1256 0.89% FillBoundary_finish() 6388 0.08021 0.09348 0.1216 0.86% MLMG::mgVcycle_down::1 128 0.1086 0.1114 0.1141 0.81% MLPoisson::Fapply() 994 0.09438 0.09897 0.1134 0.80% Amr::writePlotFile() 2 0.1005 0.1006 0.1006 0.71% Gravity::multilevel_solve_for_new_phi() 1 0.09805 0.09808 0.09812 0.70% Gravity::actual_multilevel_solve() 1 0.09797 0.09799 0.09803 0.70% MLCellLinOp::solutionResidual() 147 0.08454 0.08896 0.0926 0.66% Castro::create_source_corrector() 17 0.04701 0.06652 0.09006 0.64% Castro::post_regrid() 2 0.08573 0.08598 0.08657 0.61% MLMG::mgVcycle_bottom 128 0.0003733 0.005623 0.08423 0.60% MLMG::actualBottomSolve() 128 0.0001356 0.005378 0.08395 0.60% MLCGSolver::bicgstab 8 0 0.005145 0.08232 0.58% Castro::construct_old_source() 120 0.07701 0.07749 0.07846 0.56% MLMG::computeResidual() 128 0.07233 0.07366 0.07431 0.53% Castro::construct_old_gravity() 20 0.04587 0.06241 0.07356 0.52% FabArray::Saxpy() 192 0.06857 0.06946 0.07324 0.52% Castro::derive() 18 0.06951 0.07051 0.07224 0.51% AmrLevel::derive() 18 0.06947 0.07047 0.07219 0.51% MLMG::mgVcycle_down::2 128 0.01301 0.04074 0.07184 0.51% Castro::construct_new_gravity_source() 17 0.06883 0.06953 0.07073 0.50% Gravity::get_old_grav_vector() 20 0.03936 0.05892 0.06504 0.46% Castro::post_timestep() 1 0.05845 0.06057 0.06292 0.45% FabArray::Xpay() 953 0.04961 0.05203 0.05926 0.42% MLMG::mgVcycle_down::3 128 0.01237 0.01842 0.05814 0.41% Castro::initialize_advance() 1 0.04876 0.0502 0.05627 0.40% amrex::average_down 768 0.03331 0.03937 0.05462 0.39% Castro::check_for_nan() 37 0.04505 0.04565 0.04608 0.33% MultiFab::contains_nan() 37 0.04494 0.04553 0.04594 0.33% Castro::construct_old_gravity_source() 20 0.04431 0.0446 0.04502 0.32% Amr::regrid() 1 0.044 0.04423 0.04451 0.32% MLPoisson::define() 19 0.02976 0.03423 0.04003 0.28% Castro::retry_advance_ctu() 11 0.03623 0.03659 0.03834 0.27% Gravity::get_new_grav_vector() 18 0.03489 0.03536 0.03624 0.26% StateDataPhysBCFunct::() 147 0.0006679 0.01854 0.03491 0.25% MLMG::mgVcycle_down::4 128 0.01231 0.014 0.03482 0.25% StateData::FillBoundary(geom) 143 0 0.0177 0.03399 0.24% Castro::apply_source_to_state() 37 0.02798 0.02825 0.02869 0.20% VisMF::Write(FabArray) 4 0.01523 0.01713 0.02868 0.20% MLMG::mgVcycle_up::4 128 0.0111 0.01234 0.02722 0.19% MLCellLinOp::defineAuxData() 19 0.01949 0.02096 0.02494 0.18% MLMG::prepareForSolve() 18 0.0208 0.0217 0.02377 0.17% MLMG::mgVcycle_down::5 128 0.01228 0.01322 0.02328 0.17% Castro::sum_integrated_quantities() 2 0.02109 0.02135 0.02137 0.15% Castro::save_data_for_retry() 3 0.01933 0.01961 0.02063 0.15% MLMG::mgVcycle_up::5 128 0.01117 0.01195 0.01992 0.14% MLMG::computeMLResidual() 18 0.01046 0.01487 0.0183 0.13% Amr::writeSmallPlotFile() 2 0.0152 0.0152 0.0152 0.11% Gravity::test_composite_phi() 1 0.01508 0.01511 0.01512 0.11% MLMG::ResNormInf() 146 0.01285 0.01387 0.01503 0.11% amrex::Add() 136 0.01439 0.01449 0.01481 0.11% FabArrayBase::getFB() 6388 0.01052 0.01136 0.01476 0.10% FabArrayBase::getCPC() 2005 0.009759 0.01103 0.01443 0.10% Castro::problem_post_timestep() 2 0.009173 0.01152 0.01416 0.10% Castro::wd_update() 2 0.009161 0.01148 0.01414 0.10% FabArray::norminf() 244 0.01084 0.0112 0.01394 0.10% MLLinOp::define() 19 0.003832 0.007213 0.0114 0.08% MLLinOp::defineGrids() 19 0.003742 0.007128 0.01132 0.08% FabArray::setDomainBndry() 147 0.0004564 0.005913 0.01105 0.08% MLLinOp::makeSubCommunicator() 19 0.0005529 0.004136 0.008195 0.06% MLCGSolver::ParallelAllReduce 189 0 0.0004954 0.007926 0.06% MLMG::getGradSolution() 18 0.007239 0.007416 0.007595 0.05% MLCellLinOp::compGrad() 18 0.007197 0.007374 0.007551 0.05% Castro::volWgtSum() 46 0.007273 0.007323 0.007432 0.05% FabArray::mult() 59 0.006882 0.006998 0.007284 0.05% FabArrayBase::CPC::define() 576 0.002282 0.003456 0.006628 0.05% MLCellLinOp::defineBC() 19 0.005673 0.00592 0.006264 0.04% amrex::Dot() 142 0 0.0003823 0.006117 0.04% Castro::volInBoundary() 14 0.003369 0.003631 0.005417 0.04% BndryData::define() 19 0.004059 0.004279 0.004517 0.03% FabArray::LinComb() 71 0 0.0002301 0.003681 0.03% MLPoisson::normalize() 79 0 0.0001959 0.003135 0.02% Gravity::update_max_rhs() 2 0.001414 0.002057 0.002318 0.02% MLPoisson::prepareForSolve() 19 0.00182 0.001933 0.002241 0.02% MLCellLinOp::prepareForSolve() 19 0.001732 0.001849 0.002156 0.02% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.001325 0.001608 0.001902 0.01% AmrLevel::AmrLevel(dm) 1 0.001293 0.001464 0.001585 0.01% StateData::define() 6 0.001258 0.001428 0.001549 0.01% FabArrayBase::FB::FB() 75 0.0002896 0.0007061 0.0015 0.01% TagBoxArray::collate() 2 0.0001021 0.0006865 0.001425 0.01% Castro::buildMetrics() 1 0.0006197 0.0006911 0.001405 0.01% Gravity::swapTimeLevels() 8 0.00115 0.001233 0.001308 0.01% MLMG::MLResNormInf() 18 0.001254 0.00127 0.001288 0.01% MLMG::compResidual() 1 0.0009466 0.001107 0.001274 0.01% MLMG::MLRhsNormInf() 18 0.001198 0.001211 0.001229 0.01% Castro::locWgtSum() 6 0.001147 0.001175 0.001213 0.01% MultiFab::max() 2 0.000169 0.000812 0.001074 0.01% MLLinOp::makeAgglomeratedDMap 19 0.0009652 0.001014 0.001056 0.01% MLCellLinOp::setLevelBC() 19 0.0005443 0.0006954 0.0008586 0.01% Castro::swap_state_time_levels() 8 0.0006278 0.0006744 0.0007107 0.01% makeSFC 95 0.0006244 0.0006555 0.0006933 0.00% Castro::errorEst() 2 0.0005206 0.0005804 0.0006418 0.00% Gravity::applyMetricTerms() 2 0.000476 0.0004926 0.0005817 0.00% Castro::apply_tagging_restrictions() 2 0.0004219 0.0004797 0.000541 0.00% Castro::enforce_consistent_e() 1 0.0003608 0.0003867 0.0004293 0.00% Castro::volProductSum() 2 0.0003953 0.0003985 0.0004016 0.00% Castro::finalize_advance() 1 0.0002611 0.0002707 0.0002827 0.00% DistributionMapping::Distribute() 95 0.0002196 0.0002319 0.0002479 0.00% Castro::finalize_do_advance() 17 0.0001827 0.0002088 0.0002299 0.00% OwnerMask() 2 0.000138 0.0001533 0.0001684 0.00% Castro::enforce_speed_limit() 79 0.0001088 0.0001172 0.0001276 0.00% DistributionMapping::KnapSackDoIt() 1 9.431e-05 0.0001042 0.0001129 0.00% Castro::apply_problem_tags() 2 9.143e-05 9.38e-05 9.821e-05 0.00% MultiFab::Swap() 5 7.136e-05 7.716e-05 8.947e-05 0.00% Amr::initSubcycle() 1 5.706e-05 6.196e-05 6.329e-05 0.00% knapsack() 1 5.829e-05 6.081e-05 6.281e-05 0.00% Castro::construct_old_rotation_source() 20 3e-05 3.217e-05 3.679e-05 0.00% DistributionMapping::LeastUsedCPUs() 1 1.584e-05 2.551e-05 3.155e-05 0.00% Castro::estTimeStep() 10 1.306e-05 1.507e-05 2.446e-05 0.00% Gravity::set_mass_offset() 2 4.406e-06 9.355e-06 2.434e-05 0.00% Castro::construct_new_rotation_source() 17 2.192e-05 2.287e-05 2.426e-05 0.00% Castro::computeInitialDt() 2 1.072e-05 1.569e-05 2.255e-05 0.00% BoxList::parallelComplementIn() 4 1.355e-05 1.404e-05 1.45e-05 0.00% Castro::gwstrain() 2 4.256e-06 5.476e-06 9.934e-06 0.00% BoxList::complementIn 4 6.219e-06 6.496e-06 6.814e-06 0.00% Castro::FluxRegCrseInit 1 1.695e-06 2.047e-06 2.631e-06 0.00% Castro::FluxRegFineAdd() 1 1.476e-06 1.67e-06 2.067e-06 0.00% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Pinned Memory Usage: --------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max --------------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 16 12 B 13 B 18 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (23.10-4-ge470d3350ed3) finalized