Initializing AMReX (24.02-34-g3525b4a3f27e)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.02-34-g3525b4a3f27e) initialized Starting run at 09:10:20 UTC on 2024-03-04. Successfully read inputs file ... Castro git describe: 24.02-27-g465a83f5a AMReX git describe: 24.02-34-g3525b4a3f Microphysics git describe: 24.03-1-gce331bc8 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.163532806 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.000953275 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.012656438 Gravity::multilevel_solve_for_new_phi() time = 0.099734184 Gravity::fill_multipole_BCs() time = 0.012800421 ... 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.047342258 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.003999703 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.002595293 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.491992098e+28 zmom added: 0 eint added: 0 ener added: 6.931398274e+47 Castro::do_old_sources() time = 0.06324004 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.003188026 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.245996049e+28 zmom added: 0 eint added: 0 ener added: 3.465699137e+47 Castro::do_old_sources() time = 0.062438546 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.002588891 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -2.93534348e+40 ymom added: 1.622998024e+28 zmom added: 0 eint added: 0 ener added: 1.732849568e+47 Castro::do_old_sources() time = 0.059089265 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.003792197 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.114990122e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.057868902 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.316762695 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014656151 Gravity::solve_for_phi() time = 0.082024068 on level 0 Castro::construct_new_gravity() time = 0.084153901 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.579871373e+27 zmom added: 0 eint added: 0 ener added: 6.209484931e+43 Castro::do_new_sources() time = 0.058574128 on level 0 Ending SDC iteration 1 of 3. Beginning SDC iteration 2 of 3. Castro::construct_old_gravity() time = 0.002579575 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.114990122e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.05785748 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.316583542 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014237027 Gravity::solve_for_phi() time = 0.08117976 on level 0 Castro::construct_new_gravity() time = 0.083267196 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.583821633e+27 zmom added: 0 eint added: 0 ener added: 4.570264853e+43 Castro::do_new_sources() time = 0.058573621 on level 0 Ending SDC iteration 2 of 3. Beginning SDC iteration 3 of 3. Castro::construct_old_gravity() time = 0.002580199 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.114990122e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.057749928 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.320338389 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014179542 Gravity::solve_for_phi() time = 0.081424829 on level 0 Castro::construct_new_gravity() time = 0.083567278 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.577296503e+27 zmom added: 0 eint added: 0 ener added: 4.569279535e+43 Castro::do_new_sources() time = 0.058567262 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.002519092 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.60340448e+26 zmom added: 0 eint added: 0 ener added: 8.721813579e+46 Castro::do_old_sources() time = 0.057932281 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.323396173 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014759021 Gravity::solve_for_phi() time = 0.083595391 on level 0 Castro::construct_new_gravity() time = 0.085736611 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: -3.010508633e+24 zmom added: 0 eint added: 0 ener added: 5.968608016e+43 Castro::do_new_sources() time = 0.058655228 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.002568914 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.60340448e+26 zmom added: 0 eint added: 0 ener added: 8.721813579e+46 Castro::do_old_sources() time = 0.057869794 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.373421935 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.028594318 Gravity::solve_for_phi() time = 0.100965553 on level 0 Castro::construct_new_gravity() time = 0.103065763 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: -4.309159416e+24 zmom added: 0 eint added: 0 ener added: 5.2374099e+43 Castro::do_new_sources() time = 0.061620038 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. Castro::construct_old_gravity() time = 0.008739962 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.528602194e+26 zmom added: 0 eint added: 0 ener added: 8.780735352e+46 Castro::do_old_sources() time = 0.05843425 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.341534039 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015111321 Gravity::solve_for_phi() time = 0.08269506 on level 0 Castro::construct_new_gravity() time = 0.084829403 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: 1.223092919e+24 zmom added: 0 eint added: 0 ener added: 5.96849108e+43 Castro::do_new_sources() time = 0.071462424 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.002981961 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.528602194e+26 zmom added: 0 eint added: 0 ener added: 8.780735352e+46 Castro::do_old_sources() time = 0.058018527 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.359018481 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.027171288 Gravity::solve_for_phi() time = 0.095039265 on level 0 Castro::construct_new_gravity() time = 0.097181202 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.448585919e+24 zmom added: 0 eint added: 0 ener added: 5.233809214e+43 Castro::do_new_sources() time = 0.068365468 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.007106799 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.528602194e+26 zmom added: 0 eint added: 0 ener added: 8.841018465e+46 Castro::do_old_sources() time = 0.057881074 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.354240641 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.041507527 Gravity::solve_for_phi() time = 0.109103898 on level 0 Castro::construct_new_gravity() time = 0.111244606 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.871043233e+23 zmom added: 0 eint added: 0 ener added: 4.637081169e+43 Castro::do_new_sources() time = 0.05838671 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.010116206 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.528602194e+26 zmom added: 0 eint added: 0 ener added: 8.841018465e+46 Castro::do_old_sources() time = 0.06736478 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.344824098 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01611219 Gravity::solve_for_phi() time = 0.087020627 on level 0 Castro::construct_new_gravity() time = 0.089128794 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.359798244e+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.254813304e+24 zmom added: 0 eint added: 0 ener added: 3.901081726e+43 Castro::do_new_sources() time = 0.059371321 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.002618924 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.429621393e+26 zmom added: 0 eint added: 0 ener added: 8.902573838e+46 Castro::do_old_sources() time = 0.058081271 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.328390475 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015563267 Gravity::solve_for_phi() time = 0.082622878 on level 0 Castro::construct_new_gravity() time = 0.084768554 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: 5.133212367e+24 zmom added: 0 eint added: 0 ener added: 2.123241032e+43 Castro::do_new_sources() time = 0.058493772 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.002565369 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.429621393e+26 zmom added: 0 eint added: 0 ener added: 8.902573838e+46 Castro::do_old_sources() time = 0.057880901 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.331229003 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.025695922 Gravity::solve_for_phi() time = 0.116608852 on level 0 Castro::construct_new_gravity() time = 0.118771638 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: 8.029203612e+24 zmom added: 0 eint added: 0 ener added: 1.379533843e+43 Castro::do_new_sources() time = 0.063016462 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.920709246e+29 ener added: 4.920709246e+29 Castro::construct_old_gravity() time = 0.003838169 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.585270592e+26 zmom added: 0 eint added: 0 ener added: 8.96490842e+46 Castro::do_old_sources() time = 0.058146295 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.336859276 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014706776 Gravity::solve_for_phi() time = 0.084101919 on level 0 Castro::construct_new_gravity() time = 0.08612856 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.854195607e+43 ener added: 1.71395957e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.848617495e+35 ymom added: -5.170991299e+24 zmom added: 0 eint added: 0 ener added: -1.366537222e+43 Castro::do_new_sources() time = 0.059026031 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.920709246e+29 ener added: 4.920709246e+29 Castro::construct_old_gravity() time = 0.00257916 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.585270592e+26 zmom added: 0 eint added: 0 ener added: 8.96490842e+46 Castro::do_old_sources() time = 0.058274846 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.341755723 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014828752 Gravity::solve_for_phi() time = 0.08497257 on level 0 Castro::construct_new_gravity() time = 0.087018333 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: -7.880449068e+24 zmom added: 0 eint added: 0 ener added: -2.113374984e+43 Castro::do_new_sources() time = 0.06058 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.075345155e+29 ener added: 2.075345155e+29 Castro::construct_old_gravity() time = 0.002829163 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.439443915e+26 zmom added: 0 eint added: 0 ener added: 9.027789929e+46 Castro::do_old_sources() time = 0.0590547 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.338340624 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015014766 Gravity::solve_for_phi() time = 0.084005188 on level 0 Castro::construct_new_gravity() time = 0.086030489 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: 7.059937892e+24 zmom added: 0 eint added: 0 ener added: -4.049549507e+43 Castro::do_new_sources() time = 0.060242301 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.075345155e+29 ener added: 2.075345155e+29 Castro::construct_old_gravity() time = 0.002590114 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.439443915e+26 zmom added: 0 eint added: 0 ener added: 9.027789929e+46 Castro::do_old_sources() time = 0.058791827 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.341004941 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015153261 Gravity::solve_for_phi() time = 0.084731581 on level 0 Castro::construct_new_gravity() time = 0.086745019 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: 8.771795742e+23 zmom added: 0 eint added: 0 ener added: -4.804462625e+43 Castro::do_new_sources() time = 0.060479133 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: 1.383563436e+30 ener added: 1.383563436e+30 Castro::construct_old_gravity() time = 0.002552975 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.471933797e+26 zmom added: 0 eint added: 0 ener added: 9.091353867e+46 Castro::do_old_sources() time = 0.057800548 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.338344938 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015244792 Gravity::solve_for_phi() time = 0.084993757 on level 0 Castro::construct_new_gravity() time = 0.087091784 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.939762636e+35 ymom added: 1.676440101e+23 zmom added: 0 eint added: 0 ener added: -1.225922477e+43 Castro::do_new_sources() time = 0.06072578 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.383563436e+30 ener added: 1.383563436e+30 Castro::construct_old_gravity() time = 0.002699759 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.471933797e+26 zmom added: 0 eint added: 0 ener added: 9.091353867e+46 Castro::do_old_sources() time = 0.059198988 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.345504677 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016681736 Gravity::solve_for_phi() time = 0.086495313 on level 0 Castro::construct_new_gravity() time = 0.089584959 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: 5.748300601e+24 zmom added: 0 eint added: 0 ener added: -1.982551467e+43 Castro::do_new_sources() time = 0.060633279 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002762189715 [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: 8.037844726e+29 ener added: 8.037844726e+29 TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -4.340632835e+38 TIME= 0.25 YMOM = 1.156700224e+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.469523857e+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 = 1.213668339e-07 TIME= 0.25 CENTER OF MASS Y-VEL = 4.536085416e-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.9212152 [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.046172284 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004056468 seconds Ending run at 09:10:34 UTC on 2024-03-04. Run time = 14.72413774 Run time without initialization = 11.97221849 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.72 ... 14.72 ... 14.72 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.3878 5.4649 5.5715 37.84% Amr::InitAmr() 1 1.863 1.863 1.863 12.65% Castro::computeTemp() 80 0.8433 0.956 1.112 7.55% amrex::Copy() 724 0.7586 0.7652 0.7857 5.34% FabArray::ParallelCopy_finish() 1429 0.05046 0.4883 0.7041 4.78% Castro::normalize_species() 79 0.5467 0.5578 0.5811 3.95% MultiFab::AddProduct() 3439 0.5427 0.5474 0.5627 3.82% Castro::reset_internal_energy(Fab) 160 0.5202 0.5258 0.538 3.65% MLPoisson::Fsmooth() 6144 0.3158 0.3447 0.4121 2.80% Castro::initMFs() 1 0.3633 0.3668 0.3716 2.52% FabArray::setVal() 5013 0.3434 0.3474 0.3536 2.40% MultiFab::Subtract() 159 0.3384 0.341 0.3476 2.36% Gravity::fill_multipole_BCs() 19 0.2081 0.2904 0.3366 2.29% check_for_negative_density() 20 0.1117 0.2281 0.2926 1.99% Castro::reset_internal_energy(MultiFab) 80 0.02218 0.07157 0.2726 1.85% FabArray::ParallelCopy_nowait() 1429 0.1864 0.2013 0.2174 1.48% FabArray::sum() 3439 0.1874 0.1892 0.1942 1.32% Castro::initData() 1 0.003305 0.02637 0.1546 1.05% FillBoundary_finish() 6380 0.09535 0.1218 0.1486 1.01% FillBoundary_nowait() 6380 0.07088 0.09424 0.1189 0.81% MLPoisson::Fapply() 986 0.09419 0.09911 0.1126 0.77% MLCellLinOp::applyBC() 7148 0.04626 0.06173 0.1022 0.69% FabArray::Saxpy() 264 0.06918 0.07062 0.07649 0.52% Castro::construct_new_gravity_source() 17 0.06875 0.06958 0.07073 0.48% Castro::do_new_sources() 17 0.0003502 0.01111 0.06928 0.47% Castro::enforce_min_density() 79 0.01718 0.02891 0.06904 0.47% FabArray::Xpay() 945 0.05023 0.05275 0.0593 0.40% MultiFab::contains_nan() 37 0.04536 0.04599 0.04708 0.32% Castro::construct_old_gravity_source() 20 0.0443 0.04473 0.04667 0.32% MLMG::addInterpCorrection() 768 0.03332 0.03587 0.04233 0.29% AmrLevel::derive() 18 0.04039 0.04073 0.0414 0.28% Castro::do_old_sources() 20 0.000461 0.008277 0.03681 0.25% amrex::average_down 768 0.0244 0.02656 0.03288 0.22% StateData::FillBoundary(geom) 143 0 0.01699 0.03245 0.22% VisMF::Write(FabArray) 4 0.01203 0.01309 0.02339 0.16% Castro::construct_new_source() 102 0.01535 0.01579 0.01743 0.12% MLCellLinOp::defineAuxData() 19 0.01421 0.01485 0.01694 0.12% FabArray::FillBoundary() 6380 0.0132 0.01369 0.01632 0.11% MLCGSolver::bicgstab 8 0 0.001006 0.0161 0.11% Amr::writePlotFile() 2 0.00638 0.01485 0.01605 0.11% FabArray::norminf() 282 0.01123 0.01171 0.01562 0.11% amrex::Add() 128 0.01432 0.0145 0.01519 0.10% TagBoxArray::collate() 2 0.0001007 0.002402 0.01504 0.10% Castro::construct_old_gravity() 20 0.0003477 0.003474 0.01414 0.10% Castro::evaluate_source_change() 196 0.01217 0.0124 0.01271 0.09% FillPatchIterator::Initialize 147 0.01213 0.01234 0.0126 0.09% FabArrayBase::getFB() 6380 0.009818 0.01025 0.01256 0.09% MLCellLinOp::smooth() 3072 0.01073 0.01122 0.01199 0.08% Castro::wd_update() 2 0.005219 0.008879 0.01182 0.08% Castro::construct_old_source() 120 0.009763 0.01001 0.01163 0.08% FabArray::setDomainBndry() 147 0.0004332 0.005946 0.0112 0.08% MLMG::ResNormInf() 146 0.005089 0.007655 0.01103 0.07% MLMG::prepareForSolve() 18 0.009753 0.009958 0.01038 0.07% MLMG::solve() 18 0.001259 0.003642 0.008248 0.06% FabArray::mult() 59 0.006953 0.007092 0.00789 0.05% MLLinOp::makeSubCommunicator() 19 0.0006941 0.004068 0.007791 0.05% Castro::volWgtSum() 46 0.007337 0.007455 0.007748 0.05% Gravity::get_old_grav_vector() 20 0.006316 0.006459 0.007568 0.05% Gravity::actual_solve_with_mlmg() 19 0.006488 0.006794 0.007539 0.05% FabArrayBase::getCPC() 2005 0.007087 0.007205 0.007462 0.05% Gravity::get_new_grav_vector() 18 0.00603 0.006153 0.006472 0.04% MLCellLinOp::compGrad() 18 0.005672 0.00578 0.006396 0.04% FabArrayBase::CPC::define() 576 0.002193 0.003316 0.006315 0.04% amrex::Dot() 142 0 0.0003758 0.006013 0.04% Castro::volInBoundary() 14 0.003465 0.00374 0.005488 0.04% Amr::writeSmallPlotFile() 2 0.003517 0.005078 0.005349 0.04% MLCGSolver::ParallelAllReduce 227 0 0.0003086 0.004938 0.03% Castro::initialize_do_advance() 20 0.004263 0.004384 0.004649 0.03% MLCellLinOp::apply() 986 0.001955 0.00217 0.004607 0.03% FabArray::ParallelCopy() 1429 0.003447 0.00357 0.003771 0.03% BndryData::define() 19 0.003044 0.003161 0.003413 0.02% Castro::sum_integrated_quantities() 2 0.002224 0.003104 0.00341 0.02% MLPoisson::normalize() 79 0 0.0001973 0.003157 0.02% MultiFab::max() 2 0.000174 0.001188 0.002903 0.02% AmrLevel::FillPatch() 147 0.002038 0.002136 0.002388 0.02% StateData::define() 6 0.001272 0.001489 0.002264 0.02% MLLinOp::defineGrids() 19 0.001724 0.001831 0.002199 0.01% Castro::construct_new_gravity() 17 0.0002564 0.0007294 0.002133 0.01% MLCellLinOp::prepareForSolve() 19 0.001643 0.001806 0.002055 0.01% MLCellLinOp::correctionResidual() 768 0.001786 0.001842 0.001908 0.01% MLMG::mgVcycle() 128 0.001585 0.001656 0.001842 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.001314 0.001537 0.001745 0.01% Gravity::solve_for_phi() 17 0.0007585 0.0009441 0.001668 0.01% MLCellLinOp::defineBC() 19 0.001478 0.001552 0.001645 0.01% MLMG:computeResOfCorrection() 768 0.001374 0.001426 0.001507 0.01% FabArrayBase::FB::FB() 75 0.0002726 0.0006648 0.001413 0.01% Castro::save_data_for_retry() 3 0.001257 0.001274 0.001296 0.01% main() 1 0.0009243 0.001078 0.001196 0.01% Castro::volProductSum() 2 0.0003931 0.0004476 0.00119 0.01% Castro::locWgtSum() 6 0.001116 0.001137 0.001165 0.01% MLMG::actualBottomSolve() 128 0.0001331 0.0001937 0.001062 0.01% StateDataPhysBCFunct::() 147 0.0006379 0.0007719 0.0008916 0.01% MLCellLinOp::setLevelBC() 19 0.0005318 0.0006694 0.000819 0.01% FillPatchIterator::FillFromLevel0() 147 0.0006754 0.0006935 0.0007229 0.00% Castro::subcycle_advance_ctu() 1 0.0001788 0.0002225 0.0007053 0.00% Castro::buildMetrics() 1 0.0005585 0.0006127 0.0007003 0.00% MLMG::mgVcycle_down::2 128 0.0006026 0.0006441 0.0006852 0.00% MLMG::mgVcycle_down::3 128 0.0005651 0.0006046 0.0006826 0.00% MLMG::mgVcycle_down::0 128 0.0006243 0.0006484 0.0006724 0.00% MLMG::mgVcycle_down::1 128 0.0006277 0.0006468 0.0006682 0.00% MLMG::mgVcycle_down::5 128 0.0005767 0.0005947 0.0006663 0.00% MLMG::mgVcycle_down::4 128 0.0005557 0.0005772 0.0006589 0.00% Castro::Castro() 1 0.0004568 0.0005182 0.0006418 0.00% Castro::apply_tagging_restrictions() 2 0.0004289 0.0004763 0.0005219 0.00% AmrMesh::MakeNewGrids() 2 0.0004807 0.0004977 0.000517 0.00% FillPatchSingleLevel 147 0.0004678 0.0004759 0.0004915 0.00% MLMG::mgVcycle_up::0 128 0.0004504 0.0004626 0.0004883 0.00% Gravity::applyMetricTerms() 2 0.0004768 0.0004804 0.0004868 0.00% MLMG::mgVcycle_up::1 128 0.0004297 0.000451 0.0004738 0.00% Castro::clean_state() 79 0.0004108 0.000437 0.0004711 0.00% makeSFC 95 0.000367 0.0003923 0.0004679 0.00% MLMG::mgVcycle_up::2 128 0.0003826 0.0004117 0.000463 0.00% Castro::swap_state_time_levels() 8 0.0004072 0.000415 0.0004489 0.00% Castro::do_advance_ctu() 20 0.0004014 0.0004244 0.0004419 0.00% MLMG::mgVcycle_up::5 128 0.0003921 0.0004045 0.0004391 0.00% MLCellLinOp::solutionResidual() 147 0.0004133 0.000424 0.0004336 0.00% MLMG::mgVcycle_up::4 128 0.0003817 0.0003919 0.0004298 0.00% MLMG::mgVcycle_up::3 128 0.0003779 0.0003936 0.000429 0.00% MLMG::oneIter() 128 0.0003665 0.0003779 0.0003957 0.00% Castro::enforce_consistent_e() 1 0.0003652 0.000377 0.0003903 0.00% Castro::problem_post_timestep() 2 1.084e-05 3.52e-05 0.0003683 0.00% MLLinOp::makeAgglomeratedDMap 19 0.0003284 0.0003406 0.0003612 0.00% DistributionMapping::Distribute() 95 0.0002076 0.0002182 0.0002792 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002111 0.0002213 0.0002599 0.00% MLMG::computeResidual() 128 0.0002372 0.0002428 0.0002503 0.00% MLMG::mgVcycle_bottom 128 0.0002189 0.0002274 0.0002427 0.00% MLMG::compResidual() 1 0.0001587 0.0001649 0.0001817 0.00% Castro::finalize_do_advance() 17 0.0001455 0.00016 0.0001722 0.00% Gravity::update_max_rhs() 2 0.0001576 0.000162 0.0001707 0.00% OwnerMask() 2 0.0001319 0.0001473 0.0001636 0.00% MLPoisson::define() 19 0.0001237 0.0001354 0.0001599 0.00% Amr::defBaseLevel() 1 9.177e-05 0.000122 0.0001487 0.00% Amr::coarseTimeStep() 1 0.0001075 0.0001281 0.0001328 0.00% Castro::retry_advance_ctu() 11 5.121e-05 5.831e-05 0.0001294 0.00% TagBoxArray::mapPRD 2 0.0001053 0.0001113 0.000121 0.00% Castro::check_for_nan() 37 9.981e-05 0.0001068 0.0001204 0.00% Castro::enforce_speed_limit() 79 9.922e-05 0.0001059 0.0001179 0.00% Castro::initialize_advance() 1 9.649e-05 0.0001035 0.0001139 0.00% Amr::grid_places() 2 7.279e-05 8.247e-05 0.0001139 0.00% Gravity::test_composite_phi() 1 9.266e-05 9.633e-05 0.0001101 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.657e-05 6.688e-05 0.00011 0.00% Castro::apply_source_to_state() 37 9.059e-05 9.569e-05 0.0001015 0.00% Castro::apply_problem_tags() 2 9.246e-05 9.501e-05 0.0001015 0.00% MLLinOp::define() 19 7.631e-05 8.155e-05 9.902e-05 0.00% DistributionMapping::LeastUsedCPUs() 1 3.153e-05 8.717e-05 9.901e-05 0.00% MLPoisson::prepareForSolve() 19 7.657e-05 7.983e-05 8.705e-05 0.00% MultiFab::Swap() 5 6.385e-05 6.789e-05 8.413e-05 0.00% Amr::initSubcycle() 1 5.638e-05 6.071e-05 7.206e-05 0.00% knapsack() 1 5.559e-05 5.973e-05 6.866e-05 0.00% AmrLevel::AmrLevel(dm) 1 3.141e-05 3.683e-05 6.633e-05 0.00% Castro::create_source_corrector() 17 5.406e-05 5.87e-05 6.554e-05 0.00% Castro::expand_state() 20 4.994e-05 5.287e-05 5.542e-05 0.00% Castro::derive() 18 3.978e-05 4.171e-05 5.026e-05 0.00% MLMG::computeMLResidual() 18 3.911e-05 4.237e-05 4.942e-05 0.00% Amr::FinalizeInit() 1 8.213e-06 1.682e-05 4.814e-05 0.00% MLMG::getGradSolution() 18 3.799e-05 3.936e-05 4.152e-05 0.00% MLMG::MLResNormInf() 18 3.293e-05 3.516e-05 4.097e-05 0.00% Gravity::swapTimeLevels() 8 3.228e-05 3.36e-05 3.493e-05 0.00% Amr::timeStep() 1 1.205e-05 1.848e-05 3.41e-05 0.00% MLMG::MLRhsNormInf() 18 3.099e-05 3.192e-05 3.284e-05 0.00% Castro::construct_old_rotation_source() 20 2.996e-05 3.053e-05 3.164e-05 0.00% Amr::regrid() 1 1.086e-05 1.273e-05 2.956e-05 0.00% Gravity::actual_multilevel_solve() 1 2.626e-05 2.698e-05 2.758e-05 0.00% Amr::InitializeInit() 1 6.368e-06 9.132e-06 2.646e-05 0.00% DistributionMapping::KnapSackDoIt() 1 1.603e-05 1.784e-05 2.33e-05 0.00% Castro::construct_new_rotation_source() 17 2.1e-05 2.172e-05 2.267e-05 0.00% Castro::finalize_advance() 1 1.676e-05 1.79e-05 2.243e-05 0.00% Castro::advance() 1 7.232e-06 8.651e-06 1.931e-05 0.00% Castro::post_regrid() 2 9.878e-06 1.101e-05 1.717e-05 0.00% Castro::computeInitialDt() 2 6.79e-06 8.139e-06 1.459e-05 0.00% Castro::estTimeStep() 10 1.205e-05 1.313e-05 1.429e-05 0.00% Castro::post_init() 1 1.179e-05 1.235e-05 1.38e-05 0.00% Amr::initialInit() 1 2.945e-06 4e-06 1.038e-05 0.00% Amr::bldFineLevels() 1 6.63e-06 7.622e-06 9.928e-06 0.00% Amr::init() 1 6.311e-06 7.158e-06 8.359e-06 0.00% BoxList::parallelComplementIn() 4 6.843e-06 7.202e-06 7.723e-06 0.00% Castro::errorEst() 2 6.021e-06 6.317e-06 6.689e-06 0.00% BoxList::complementIn 4 5.884e-06 6.1e-06 6.491e-06 0.00% Castro::gwstrain() 2 4.221e-06 4.752e-06 6.429e-06 0.00% Castro::post_timestep() 1 5.291e-06 5.857e-06 6.33e-06 0.00% Gravity::set_mass_offset() 2 3.974e-06 4.157e-06 4.4e-06 0.00% Castro::FluxRegCrseInit 1 1.599e-06 1.921e-06 2.312e-06 0.00% Castro::FluxRegFineAdd() 1 1.406e-06 1.49e-06 1.555e-06 0.00% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.72 14.72 14.72 100.00% Amr::coarseTimeStep() 1 11.97 11.97 11.97 81.31% Amr::timeStep() 1 11.92 11.92 11.92 80.96% Castro::advance() 1 11.86 11.86 11.86 80.57% Castro::subcycle_advance_ctu() 1 11.81 11.81 11.81 80.23% Castro::do_advance_ctu() 20 11.77 11.77 11.78 79.98% Castro::construct_ctu_hydro_source() 20 6.431 6.479 6.544 44.44% Castro::clean_state() 79 3.663 3.787 3.985 27.06% Castro::computeTemp() 80 2.364 2.493 2.755 18.71% Amr::InitAmr() 1 1.863 1.863 1.863 12.65% Castro::reset_internal_energy(MultiFab) 80 1.474 1.537 1.727 11.73% Gravity::solve_phi_with_mlmg() 19 1.485 1.57 1.62 11.00% Castro::construct_new_gravity() 17 1.411 1.494 1.543 10.48% Gravity::solve_for_phi() 17 1.373 1.457 1.507 10.24% Castro::initialize_do_advance() 20 1.354 1.398 1.433 9.73% Gravity::actual_solve_with_mlmg() 19 1.272 1.276 1.28 8.70% MLMG::solve() 18 1.224 1.227 1.229 8.34% Castro::do_old_sources() 20 1.139 1.16 1.173 7.96% MLMG::oneIter() 128 1.089 1.092 1.094 7.43% MLMG::mgVcycle() 128 1.074 1.077 1.079 7.33% Castro::do_new_sources() 17 0.933 0.979 1.037 7.04% Castro::evaluate_source_change() 196 0.9371 0.9454 0.9713 6.60% FabArray::ParallelCopy() 1429 0.2729 0.7123 0.9277 6.30% Amr::init() 1 0.8881 0.8882 0.8882 6.03% Amr::initialInit() 1 0.836 0.8361 0.8361 5.68% Castro::enforce_min_density() 79 0.7765 0.795 0.8297 5.64% amrex::Copy() 724 0.7586 0.7652 0.7857 5.34% MLCellLinOp::smooth() 3072 0.5549 0.6013 0.7226 4.91% FabArray::ParallelCopy_finish() 1429 0.05046 0.4883 0.7041 4.78% AmrLevel::FillPatch() 147 0.4621 0.5524 0.6458 4.39% Amr::InitializeInit() 1 0.4485 0.5074 0.6101 4.14% Amr::defBaseLevel() 1 0.4485 0.5074 0.61 4.14% Castro::normalize_species() 79 0.5467 0.5578 0.5811 3.95% MultiFab::AddProduct() 3439 0.5427 0.5474 0.5627 3.82% Castro::reset_internal_energy(Fab) 160 0.5202 0.5258 0.538 3.65% FillPatchIterator::Initialize 147 0.29 0.3791 0.4736 3.22% FillPatchIterator::FillFromLevel0() 147 0.2664 0.3608 0.457 3.10% FillPatchSingleLevel 147 0.2657 0.3601 0.4563 3.10% MLPoisson::Fsmooth() 6144 0.3158 0.3447 0.4121 2.80% MLMG::mgVcycle_up::1 128 0.09474 0.2514 0.4087 2.78% Amr::FinalizeInit() 1 0.2261 0.3286 0.3875 2.63% Castro::Castro() 1 0.3688 0.3723 0.3773 2.56% MLCellLinOp::applyBC() 7148 0.2654 0.3023 0.373 2.53% Castro::initMFs() 1 0.3635 0.367 0.3718 2.52% MLMG::addInterpCorrection() 768 0.06163 0.2927 0.3582 2.43% FabArray::setVal() 5013 0.3434 0.3474 0.3536 2.40% MultiFab::Subtract() 159 0.3384 0.341 0.3476 2.36% Gravity::fill_multipole_BCs() 19 0.2102 0.2925 0.3388 2.30% check_for_negative_density() 20 0.1117 0.2281 0.2926 1.99% FabArray::FillBoundary() 6380 0.2092 0.2406 0.275 1.87% MLMG::mgVcycle_up::2 128 0.01053 0.1111 0.2653 1.80% MLMG::mgVcycle_down::0 128 0.2585 0.2611 0.2636 1.79% FabArray::ParallelCopy_nowait() 1429 0.2051 0.2204 0.2374 1.61% Castro::initData() 1 0.07751 0.1333 0.2339 1.59% Amr::bldFineLevels() 1 0.04712 0.1495 0.2075 1.41% Castro::expand_state() 20 0.1731 0.193 0.2059 1.40% MLMG::mgVcycle_up::0 128 0.2022 0.2042 0.2058 1.40% FabArray::sum() 3439 0.1874 0.1892 0.1942 1.32% MLCellLinOp::apply() 986 0.1443 0.1567 0.1902 1.29% MLMG::mgVcycle_up::3 128 0.01046 0.02292 0.1728 1.17% Amr::grid_places() 2 0.003008 0.1056 0.1645 1.12% AmrMesh::MakeNewGrids() 2 0.002919 0.1055 0.1644 1.12% TagBoxArray::mapPRD 2 0.0006325 0.1019 0.1621 1.10% FillBoundary_finish() 6380 0.09535 0.1218 0.1486 1.01% MLMG:computeResOfCorrection() 768 0.1097 0.1193 0.1412 0.96% MLCellLinOp::correctionResidual() 768 0.1083 0.1178 0.1397 0.95% Castro::post_init() 1 0.1346 0.1356 0.1373 0.93% Castro::construct_new_source() 102 0.1235 0.1253 0.1334 0.91% FillBoundary_nowait() 6380 0.08099 0.1052 0.1303 0.89% Castro::create_source_corrector() 17 0.04671 0.07846 0.1201 0.82% MLMG::mgVcycle_down::1 128 0.1106 0.1145 0.119 0.81% MLPoisson::Fapply() 986 0.09419 0.09911 0.1126 0.77% Gravity::multilevel_solve_for_new_phi() 1 0.0997 0.09972 0.09975 0.68% Gravity::actual_multilevel_solve() 1 0.09963 0.09965 0.09969 0.68% MLCellLinOp::solutionResidual() 147 0.08714 0.09227 0.0984 0.67% Amr::writePlotFile() 2 0.09405 0.09412 0.09413 0.64% Castro::post_regrid() 2 0.0857 0.08739 0.08839 0.60% Castro::construct_old_source() 120 0.07707 0.07803 0.08427 0.57% MLMG::computeResidual() 128 0.07474 0.07587 0.07757 0.53% FabArray::Saxpy() 264 0.06918 0.07062 0.07649 0.52% Castro::construct_new_gravity_source() 17 0.06875 0.06958 0.07073 0.48% MLMG::mgVcycle_bottom 128 0.0003534 0.004758 0.07069 0.48% MLMG::actualBottomSolve() 128 0.0001331 0.004531 0.07046 0.48% Castro::derive() 18 0.06843 0.06874 0.06937 0.47% AmrLevel::derive() 18 0.06839 0.0687 0.06933 0.47% MLCGSolver::bicgstab 8 0 0.004308 0.06892 0.47% MLMG::mgVcycle_down::2 128 0.0123 0.03935 0.06804 0.46% Castro::post_timestep() 1 0.05811 0.06181 0.06454 0.44% Castro::construct_old_gravity() 20 0.04819 0.0552 0.05933 0.40% FabArray::Xpay() 945 0.05023 0.05275 0.0593 0.40% Gravity::get_old_grav_vector() 20 0.03827 0.05172 0.05888 0.40% MLMG::mgVcycle_down::3 128 0.01163 0.01767 0.05828 0.40% amrex::average_down 768 0.03244 0.03859 0.05409 0.37% Castro::initialize_advance() 1 0.04814 0.04875 0.04923 0.33% Castro::check_for_nan() 37 0.04547 0.04609 0.0472 0.32% MultiFab::contains_nan() 37 0.04536 0.04599 0.04708 0.32% Castro::construct_old_gravity_source() 20 0.0443 0.04473 0.04667 0.32% Amr::regrid() 1 0.04391 0.04481 0.04551 0.31% MLPoisson::define() 19 0.02815 0.03243 0.03776 0.26% Gravity::get_new_grav_vector() 18 0.03471 0.03585 0.03671 0.25% Castro::retry_advance_ctu() 11 0.03498 0.03524 0.03547 0.24% MLMG::mgVcycle_down::4 128 0.01165 0.01328 0.03431 0.23% StateDataPhysBCFunct::() 147 0.0006379 0.01776 0.0333 0.23% StateData::FillBoundary(geom) 143 0 0.01699 0.03245 0.22% Castro::apply_source_to_state() 37 0.02819 0.02838 0.02853 0.19% MLMG::mgVcycle_up::4 128 0.01054 0.01169 0.02662 0.18% MLCellLinOp::defineAuxData() 19 0.01832 0.01959 0.02356 0.16% VisMF::Write(FabArray) 4 0.01203 0.01309 0.02339 0.16% MLMG::prepareForSolve() 18 0.01982 0.02072 0.02277 0.15% MLMG::mgVcycle_down::5 128 0.01162 0.01251 0.02266 0.15% Castro::sum_integrated_quantities() 2 0.022 0.02227 0.02229 0.15% MLMG::ResNormInf() 146 0.01515 0.01768 0.02099 0.14% MLMG::computeMLResidual() 18 0.0111 0.01596 0.02024 0.14% MLMG::mgVcycle_up::5 128 0.01062 0.01134 0.01913 0.13% Castro::save_data_for_retry() 3 0.01832 0.01842 0.01857 0.13% Gravity::test_composite_phi() 1 0.01618 0.01621 0.01624 0.11% FabArray::norminf() 282 0.01123 0.01171 0.01562 0.11% Castro::problem_post_timestep() 2 0.009481 0.01265 0.01553 0.11% Castro::wd_update() 2 0.009112 0.01262 0.01552 0.11% amrex::Add() 128 0.01432 0.0145 0.01519 0.10% TagBoxArray::collate() 2 0.0001007 0.002402 0.01504 0.10% FabArrayBase::getFB() 6380 0.01012 0.01092 0.01397 0.09% FabArrayBase::getCPC() 2005 0.009322 0.01052 0.01378 0.09% FabArray::setDomainBndry() 147 0.0004332 0.005946 0.0112 0.08% MLLinOp::define() 19 0.003535 0.006932 0.01102 0.07% MLLinOp::defineGrids() 19 0.003454 0.006851 0.01093 0.07% Amr::writeSmallPlotFile() 2 0.008457 0.008458 0.00846 0.06% MLCGSolver::ParallelAllReduce 227 0 0.0005119 0.008191 0.06% MLMG::getGradSolution() 18 0.007228 0.007458 0.008118 0.06% MLCellLinOp::compGrad() 18 0.00719 0.007419 0.008078 0.05% FabArray::mult() 59 0.006953 0.007092 0.00789 0.05% MLLinOp::makeSubCommunicator() 19 0.0006941 0.004068 0.007791 0.05% Castro::volWgtSum() 46 0.007337 0.007455 0.007748 0.05% FabArrayBase::CPC::define() 576 0.002193 0.003316 0.006315 0.04% MLCellLinOp::defineBC() 19 0.00553 0.005768 0.006245 0.04% amrex::Dot() 142 0 0.0003758 0.006013 0.04% Castro::volInBoundary() 14 0.003465 0.00374 0.005488 0.04% BndryData::define() 19 0.003933 0.004215 0.00462 0.03% Gravity::update_max_rhs() 2 0.001421 0.002422 0.00413 0.03% MLPoisson::normalize() 79 0 0.0001973 0.003157 0.02% MultiFab::max() 2 0.000174 0.001188 0.002903 0.02% AmrLevel::AmrLevel(dm) 1 0.001307 0.001526 0.002306 0.02% StateData::define() 6 0.001272 0.001489 0.002264 0.02% MLPoisson::prepareForSolve() 19 0.001724 0.001886 0.002135 0.01% MLCellLinOp::prepareForSolve() 19 0.001643 0.001806 0.002055 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.001314 0.001537 0.001745 0.01% FabArrayBase::FB::FB() 75 0.0002726 0.0006648 0.001413 0.01% MLMG::MLResNormInf() 18 0.001289 0.001319 0.00136 0.01% MLMG::MLRhsNormInf() 18 0.00123 0.001256 0.001312 0.01% MLMG::compResidual() 1 0.0009555 0.001101 0.001241 0.01% Gravity::swapTimeLevels() 8 0.001142 0.001167 0.001239 0.01% Castro::volProductSum() 2 0.0003931 0.0004476 0.00119 0.01% Castro::locWgtSum() 6 0.001116 0.001137 0.001165 0.01% MLLinOp::makeAgglomeratedDMap 19 0.000903 0.0009511 0.001102 0.01% MLCellLinOp::setLevelBC() 19 0.0005318 0.0006694 0.000819 0.01% Castro::buildMetrics() 1 0.0006207 0.000675 0.0007625 0.01% makeSFC 95 0.0005746 0.0006105 0.0007471 0.01% Castro::errorEst() 2 0.0005318 0.0005777 0.0006217 0.00% Castro::apply_tagging_restrictions() 2 0.0004289 0.0004763 0.0005219 0.00% Gravity::applyMetricTerms() 2 0.0004768 0.0004804 0.0004868 0.00% Castro::swap_state_time_levels() 8 0.0004072 0.000415 0.0004489 0.00% Castro::enforce_consistent_e() 1 0.0003652 0.000377 0.0003903 0.00% DistributionMapping::Distribute() 95 0.0002076 0.0002182 0.0002792 0.00% Castro::finalize_advance() 1 0.0002506 0.0002573 0.0002637 0.00% Castro::finalize_do_advance() 17 0.0001548 0.00017 0.0001817 0.00% DistributionMapping::KnapSackDoIt() 1 0.0001168 0.0001647 0.0001797 0.00% OwnerMask() 2 0.0001319 0.0001473 0.0001636 0.00% Castro::enforce_speed_limit() 79 9.922e-05 0.0001059 0.0001179 0.00% Castro::apply_problem_tags() 2 9.246e-05 9.501e-05 0.0001015 0.00% DistributionMapping::LeastUsedCPUs() 1 3.153e-05 8.717e-05 9.901e-05 0.00% MultiFab::Swap() 5 6.385e-05 6.789e-05 8.413e-05 0.00% Amr::initSubcycle() 1 5.638e-05 6.071e-05 7.206e-05 0.00% knapsack() 1 5.559e-05 5.973e-05 6.866e-05 0.00% Castro::construct_old_rotation_source() 20 2.996e-05 3.053e-05 3.164e-05 0.00% Castro::construct_new_rotation_source() 17 2.1e-05 2.172e-05 2.267e-05 0.00% Castro::computeInitialDt() 2 9.682e-06 1.126e-05 1.761e-05 0.00% Castro::estTimeStep() 10 1.205e-05 1.313e-05 1.429e-05 0.00% BoxList::parallelComplementIn() 4 1.28e-05 1.33e-05 1.419e-05 0.00% BoxList::complementIn 4 5.884e-06 6.1e-06 6.491e-06 0.00% Castro::gwstrain() 2 4.221e-06 4.752e-06 6.429e-06 0.00% Gravity::set_mass_offset() 2 3.974e-06 4.157e-06 4.4e-06 0.00% Castro::FluxRegCrseInit 1 1.599e-06 1.921e-06 2.312e-06 0.00% Castro::FluxRegFineAdd() 1 1.406e-06 1.49e-06 1.555e-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 10 B 12 B 15 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.02-34-g3525b4a3f27e) finalized