Initializing AMReX (24.03-13-g3fe7aadc23c9)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.03-13-g3fe7aadc23c9) initialized Starting run at 08:11:54 UTC on 2024-03-15. Successfully read inputs file ... Castro git describe: 24.03-7-ge19c0bdaa AMReX git describe: 24.03-13-g3fe7aadc2 Microphysics git describe: 24.03-10-ga5faa5f9 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.160672916 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.00105532 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010345031 Gravity::multilevel_solve_for_new_phi() time = 0.100980485 Gravity::fill_multipole_BCs() time = 0.010337239 ... 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.050686877 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004117003 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.005541646 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.064205965 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 = -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, 0Invalid 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, 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 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 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.007311087 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.062886135 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.002545154 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.059580767 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.1060744788 at index 20, 177, 0 Invalid density = -0.1060744788 at index 20, 78, 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.002680285 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.058577197 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.318363355 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014215988 Gravity::solve_for_phi() time = 0.083476878 on level 0 Castro::construct_new_gravity() time = 0.085630788 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.05985014 on level 0 Ending SDC iteration 1 of 3. Beginning SDC iteration 2 of 3. Castro::construct_old_gravity() time = 0.002239063 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.058394169 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.31779203 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.013892922 Gravity::solve_for_phi() time = 0.083632232 on level 0 Castro::construct_new_gravity() time = 0.085764482 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.059204906 on level 0 Ending SDC iteration 2 of 3. Beginning SDC iteration 3 of 3. Castro::construct_old_gravity() time = 0.002549789 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.058344621 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.323241641 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014279546 Gravity::solve_for_phi() time = 0.083449804 on level 0 Castro::construct_new_gravity() time = 0.085572579 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.058965122 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.002537621 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.059769435 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.325795284 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01450229 Gravity::solve_for_phi() time = 0.084219789 on level 0 Castro::construct_new_gravity() time = 0.086336079 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.05909069 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.002689789 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.058361884 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.33191648 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015730567 Gravity::solve_for_phi() time = 0.084888849 on level 0 Castro::construct_new_gravity() time = 0.08705809 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.060113381 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.002540199 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.058401102 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.326233711 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014476043 Gravity::solve_for_phi() time = 0.083836824 on level 0 Castro::construct_new_gravity() time = 0.085948943 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.059125344 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.002450197 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.05821451 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.331755487 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014349754 Gravity::solve_for_phi() time = 0.084342097 on level 0 Castro::construct_new_gravity() time = 0.086459196 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.059360359 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.003394289 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.058426454 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.326229118 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014968681 Gravity::solve_for_phi() time = 0.084594934 on level 0 Castro::construct_new_gravity() time = 0.086744698 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.060409755 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.002674755 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.059402165 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.332795724 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014775822 Gravity::solve_for_phi() time = 0.084299754 on level 0 Castro::construct_new_gravity() time = 0.086415457 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.059360924 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.00266375 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.05894255 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.327065383 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015047097 Gravity::solve_for_phi() time = 0.084151587 on level 0 Castro::construct_new_gravity() time = 0.087259941 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.059198826 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.002605923 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.058494681 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.333302273 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015277459 Gravity::solve_for_phi() time = 0.087608959 on level 0 Castro::construct_new_gravity() time = 0.089749901 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.060762343 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.002790811 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.064282313 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.339707331 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015899937 Gravity::solve_for_phi() time = 0.094557038 on level 0 Castro::construct_new_gravity() time = 0.096921732 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.059348427 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.002508042 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.058359405 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.332348123 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014825842 Gravity::solve_for_phi() time = 0.085168708 on level 0 Castro::construct_new_gravity() time = 0.087326689 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.065072029 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.002946171 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.060460359 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.329804039 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015084111 Gravity::solve_for_phi() time = 0.083936343 on level 0 Castro::construct_new_gravity() time = 0.086127846 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.059478699 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.00230948 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.058179135 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.332295153 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015095129 Gravity::solve_for_phi() time = 0.084156735 on level 0 Castro::construct_new_gravity() time = 0.08630886 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.059035334 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.002511469 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.058313108 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.336114254 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.025114861 Gravity::solve_for_phi() time = 0.095661867 on level 0 Castro::construct_new_gravity() time = 0.0977788 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.059511264 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.002674896 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.058776529 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.330981927 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015527129 Gravity::solve_for_phi() time = 0.085278483 on level 0 Castro::construct_new_gravity() time = 0.087475825 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.059665636 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.00279691738 [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.77439902 [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.047416991 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004176937 seconds Ending run at 08:12:08 UTC on 2024-03-15. Run time = 14.02114964 Run time without initialization = 11.82679629 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.02 ... 14.02 ... 14.02 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.3849 5.4252 5.5408 39.52% Amr::InitAmr() 1 1.303 1.304 1.304 9.30% Castro::computeTemp() 80 0.8587 0.9597 1.123 8.01% amrex::Copy() 724 0.7832 0.7875 0.7985 5.69% FabArray::ParallelCopy_finish() 1429 0.04691 0.4797 0.6751 4.82% MultiFab::AddProduct() 3439 0.5569 0.5614 0.5719 4.08% Castro::normalize_species() 79 0.546 0.5496 0.5571 3.97% Castro::reset_internal_energy(Fab) 160 0.5174 0.5226 0.532 3.79% MLPoisson::Fsmooth() 6144 0.3142 0.3388 0.413 2.95% Castro::initMFs() 1 0.3604 0.3656 0.3862 2.75% FabArray::setVal() 5013 0.3499 0.3531 0.3661 2.61% MultiFab::Subtract() 159 0.3407 0.3434 0.3546 2.53% Gravity::fill_multipole_BCs() 19 0.1954 0.2398 0.2806 2.00% FabArray::ParallelCopy_nowait() 1429 0.187 0.2007 0.2197 1.57% check_for_negative_density() 20 0.03436 0.1629 0.2163 1.54% Castro::reset_internal_energy(MultiFab) 80 0.02274 0.05975 0.2118 1.51% FabArray::sum() 3439 0.1786 0.1802 0.183 1.31% Castro::initData() 1 0.003328 0.02562 0.1484 1.06% FillBoundary_finish() 6380 0.09901 0.115 0.1383 0.99% FillBoundary_nowait() 6380 0.07726 0.1027 0.1314 0.94% MLPoisson::Fapply() 986 0.0946 0.0992 0.1146 0.82% MLCellLinOp::applyBC() 7148 0.04695 0.06205 0.1031 0.74% FabArray::Saxpy() 264 0.07065 0.07162 0.07947 0.57% Castro::construct_new_gravity_source() 17 0.06943 0.0703 0.07141 0.51% Castro::do_new_sources() 17 0.00039 0.01254 0.06816 0.49% Castro::enforce_min_density() 79 0.01748 0.02777 0.06321 0.45% FabArray::Xpay() 945 0.05097 0.05328 0.06079 0.43% MultiFab::contains_nan() 37 0.04546 0.04587 0.04652 0.33% Castro::construct_old_gravity_source() 20 0.04477 0.0451 0.04645 0.33% Castro::do_old_sources() 20 0.0005471 0.009429 0.04331 0.31% MLMG::addInterpCorrection() 768 0.0327 0.03507 0.04183 0.30% AmrLevel::derive() 18 0.04056 0.04076 0.04104 0.29% StateData::FillBoundary(geom) 143 0 0.01692 0.03281 0.23% amrex::average_down 768 0.02431 0.02621 0.03219 0.23% VisMF::Write(FabArray) 4 0.01216 0.01344 0.02568 0.18% Amr::writePlotFile() 2 0.007377 0.01802 0.01921 0.14% MLCellLinOp::defineAuxData() 19 0.01471 0.01545 0.01776 0.13% FabArray::FillBoundary() 6380 0.01362 0.01419 0.01688 0.12% MLCGSolver::bicgstab 8 0 0.001039 0.01663 0.12% Castro::construct_new_source() 102 0.01538 0.01554 0.01591 0.11% FabArray::norminf() 282 0.01076 0.01121 0.01523 0.11% TagBoxArray::collate() 2 8.579e-05 0.002082 0.01478 0.11% amrex::Add() 128 0.0141 0.0142 0.01446 0.10% Castro::evaluate_source_change() 196 0.01276 0.01301 0.01381 0.10% MLCellLinOp::smooth() 3072 0.01119 0.01171 0.01287 0.09% FillPatchIterator::Initialize 147 0.01236 0.01253 0.01282 0.09% Castro::construct_old_gravity() 20 0.0003944 0.002482 0.01275 0.09% FabArrayBase::getFB() 6380 0.009931 0.01036 0.0127 0.09% FabArray::setDomainBndry() 147 0.0004524 0.005949 0.01167 0.08% Castro::wd_update() 2 0.005215 0.008676 0.01081 0.08% MLMG::prepareForSolve() 18 0.009762 0.01004 0.01061 0.08% Castro::construct_old_source() 120 0.009834 0.00987 0.009935 0.07% MLLinOp::makeSubCommunicator() 19 0.001011 0.004777 0.008877 0.06% MLMG::ResNormInf() 146 0.004518 0.006598 0.007917 0.06% FabArrayBase::getCPC() 2005 0.00723 0.007396 0.007796 0.06% Gravity::actual_solve_with_mlmg() 19 0.006598 0.006859 0.007737 0.06% FabArray::mult() 59 0.006927 0.007066 0.0075 0.05% MultiFab::max() 2 0.0001974 0.005929 0.007473 0.05% Gravity::get_new_grav_vector() 18 0.006123 0.006308 0.00745 0.05% Castro::volWgtSum() 46 0.007203 0.007288 0.007419 0.05% MLMG::solve() 18 0.001226 0.003171 0.00708 0.05% Gravity::get_old_grav_vector() 20 0.006423 0.006488 0.00658 0.05% FabArrayBase::CPC::define() 576 0.002213 0.003349 0.006322 0.05% amrex::Dot() 142 0 0.000379 0.006064 0.04% MLCellLinOp::compGrad() 18 0.005768 0.005825 0.005885 0.04% Amr::writeSmallPlotFile() 2 0.003688 0.005321 0.005591 0.04% Castro::volInBoundary() 14 0.00339 0.003646 0.005385 0.04% MLCGSolver::ParallelAllReduce 227 0 0.0003118 0.004988 0.04% Castro::initialize_do_advance() 20 0.00435 0.004494 0.004876 0.03% MLCellLinOp::apply() 986 0.002002 0.002217 0.004696 0.03% FabArray::ParallelCopy() 1429 0.003526 0.003667 0.003996 0.03% BndryData::define() 19 0.00311 0.003182 0.003415 0.02% MLPoisson::normalize() 79 0 0.0001946 0.003113 0.02% AmrLevel::FillPatch() 147 0.002148 0.002229 0.002384 0.02% Castro::sum_integrated_quantities() 2 0.001634 0.002112 0.00222 0.02% MLLinOp::defineGrids() 19 0.001786 0.001899 0.002177 0.02% MLCellLinOp::prepareForSolve() 19 0.001657 0.001786 0.002077 0.01% MLCellLinOp::correctionResidual() 768 0.001814 0.001861 0.001975 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.001323 0.001544 0.001781 0.01% MLMG::mgVcycle() 128 0.001572 0.001631 0.00178 0.01% MLCellLinOp::defineBC() 19 0.001524 0.001561 0.001684 0.01% Gravity::solve_for_phi() 17 0.0008102 0.0009948 0.001519 0.01% MLMG:computeResOfCorrection() 768 0.001382 0.001432 0.001505 0.01% FabArrayBase::FB::FB() 75 0.000272 0.0006713 0.001426 0.01% StateData::define() 6 0.001339 0.001358 0.001386 0.01% Castro::save_data_for_retry() 3 0.001252 0.001293 0.001341 0.01% main() 1 0.001021 0.001145 0.001261 0.01% Castro::locWgtSum() 6 0.001106 0.00111 0.001127 0.01% MLMG::actualBottomSolve() 128 0.0001322 0.0001948 0.001081 0.01% StateDataPhysBCFunct::() 147 0.0007189 0.0008521 0.001024 0.01% Castro::construct_new_gravity() 17 0.0003146 0.0004341 0.0009878 0.01% MLCellLinOp::setLevelBC() 19 0.0005536 0.0006744 0.0008649 0.01% Castro::create_source_corrector() 17 5.619e-05 0.0001095 0.0008271 0.01% Castro::subcycle_advance_ctu() 1 0.0001798 0.0002351 0.0008262 0.01% DistributionMapping::LeastUsedCPUs() 1 2.198e-05 0.0007269 0.0007911 0.01% MLMG::mgVcycle_down::2 128 0.0006192 0.0006713 0.0007564 0.01% FillPatchIterator::FillFromLevel0() 147 0.0006826 0.0007074 0.000745 0.01% MLMG::mgVcycle_down::0 128 0.0006598 0.0006739 0.0007121 0.01% MLMG::mgVcycle_down::3 128 0.000582 0.0006243 0.0007105 0.01% MLMG::mgVcycle_down::1 128 0.0006558 0.0006705 0.0007087 0.01% MLMG::mgVcycle_down::4 128 0.0005689 0.0005914 0.000688 0.00% MLMG::mgVcycle_down::5 128 0.0005919 0.0006075 0.0006716 0.00% Castro::Castro() 1 0.0005172 0.0005763 0.000643 0.00% Castro::buildMetrics() 1 0.0005797 0.0005969 0.0006277 0.00% Castro::clean_state() 79 0.0004371 0.0004775 0.0005742 0.00% Castro::apply_tagging_restrictions() 2 0.0004145 0.0004821 0.00056 0.00% Castro::do_advance_ctu() 20 0.0004346 0.0004589 0.0005343 0.00% Gravity::applyMetricTerms() 2 0.0004871 0.0004994 0.0005334 0.00% AmrMesh::MakeNewGrids() 2 0.0004973 0.0005121 0.0005234 0.00% FillPatchSingleLevel 147 0.0004696 0.0004879 0.0005094 0.00% MLMG::mgVcycle_up::1 128 0.0004555 0.0004728 0.0005091 0.00% MLMG::mgVcycle_up::0 128 0.0004535 0.0004663 0.000493 0.00% MLMG::mgVcycle_up::2 128 0.0003942 0.0004281 0.0004829 0.00% MLCellLinOp::solutionResidual() 147 0.0004324 0.000446 0.0004679 0.00% MLMG::mgVcycle_up::5 128 0.000401 0.0004155 0.000461 0.00% makeSFC 95 0.0003676 0.0003969 0.0004521 0.00% Castro::swap_state_time_levels() 8 0.0004187 0.000434 0.0004478 0.00% MLMG::mgVcycle_up::3 128 0.0003905 0.0004028 0.0004386 0.00% MLMG::oneIter() 128 0.0003944 0.0004044 0.0004384 0.00% MLMG::mgVcycle_up::4 128 0.0003924 0.000401 0.0004313 0.00% Castro::problem_post_timestep() 2 1.112e-05 3.629e-05 0.0003988 0.00% Castro::volProductSum() 2 0.0003843 0.0003874 0.000391 0.00% Castro::enforce_consistent_e() 1 0.0003628 0.0003676 0.0003806 0.00% MLLinOp::makeAgglomeratedDMap 19 0.000328 0.0003447 0.0003753 0.00% DistributionMapping::Distribute() 95 0.0002138 0.0002331 0.0002863 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002216 0.0002368 0.0002637 0.00% MLMG::computeResidual() 128 0.0002407 0.0002453 0.000257 0.00% MLMG::mgVcycle_bottom 128 0.0002264 0.0002343 0.0002477 0.00% Castro::enforce_speed_limit() 79 0.0001082 0.00012 0.000212 0.00% Castro::finalize_do_advance() 17 0.0001572 0.00017 0.0001915 0.00% Castro::retry_advance_ctu() 11 5.271e-05 6.321e-05 0.0001847 0.00% Gravity::update_max_rhs() 2 0.0001567 0.0001623 0.0001677 0.00% OwnerMask() 2 0.0001338 0.0001491 0.0001644 0.00% Amr::coarseTimeStep() 1 0.0001295 0.0001582 0.0001642 0.00% MLMG::compResidual() 1 0.0001514 0.0001524 0.0001592 0.00% Amr::initSubcycle() 1 5.536e-05 6.326e-05 0.000153 0.00% MLPoisson::define() 19 0.0001244 0.0001333 0.0001444 0.00% Gravity::test_composite_phi() 1 9.112e-05 9.974e-05 0.0001318 0.00% Amr::defBaseLevel() 1 0.0001002 0.0001166 0.0001314 0.00% Castro::initialize_advance() 1 0.0001109 0.0001143 0.0001256 0.00% TagBoxArray::mapPRD 2 0.0001109 0.0001148 0.0001224 0.00% Castro::check_for_nan() 37 0.0001041 0.0001107 0.0001221 0.00% Amr::grid_places() 2 7.901e-05 9.422e-05 0.0001208 0.00% Castro::apply_source_to_state() 37 9.544e-05 0.0001025 0.000112 0.00% Castro::apply_problem_tags() 2 9.184e-05 9.353e-05 0.0001009 0.00% MLLinOp::define() 19 7.774e-05 8.177e-05 9.892e-05 0.00% MLPoisson::prepareForSolve() 19 7.744e-05 8.43e-05 9.741e-05 0.00% Amr::FinalizeInit() 1 8.761e-06 1.879e-05 9.346e-05 0.00% MultiFab::Swap() 5 6.448e-05 6.724e-05 8.354e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.965e-05 6.56e-05 8.23e-05 0.00% AmrLevel::AmrLevel(dm) 1 4.932e-05 5.645e-05 6.552e-05 0.00% Castro::expand_state() 20 5.098e-05 5.522e-05 6.426e-05 0.00% knapsack() 1 5.527e-05 5.786e-05 6.16e-05 0.00% Castro::derive() 18 4.23e-05 4.466e-05 4.928e-05 0.00% MLMG::computeMLResidual() 18 3.91e-05 4.19e-05 4.565e-05 0.00% MLMG::getGradSolution() 18 3.872e-05 4.131e-05 4.324e-05 0.00% MLMG::MLRhsNormInf() 18 3.21e-05 3.276e-05 3.826e-05 0.00% Gravity::swapTimeLevels() 8 3.181e-05 3.449e-05 3.724e-05 0.00% MLMG::MLResNormInf() 18 3.369e-05 3.494e-05 3.633e-05 0.00% Castro::construct_old_rotation_source() 20 2.997e-05 3.077e-05 3.372e-05 0.00% DistributionMapping::KnapSackDoIt() 1 1.617e-05 1.959e-05 3.275e-05 0.00% Amr::regrid() 1 1.074e-05 1.261e-05 3.054e-05 0.00% Castro::post_init() 1 1.192e-05 1.366e-05 2.886e-05 0.00% Gravity::actual_multilevel_solve() 1 2.642e-05 2.713e-05 2.808e-05 0.00% Castro::computeInitialDt() 2 6.704e-06 8.475e-06 2.432e-05 0.00% Castro::construct_new_rotation_source() 17 2.132e-05 2.19e-05 2.372e-05 0.00% Castro::finalize_advance() 1 1.69e-05 1.818e-05 2.227e-05 0.00% Amr::InitializeInit() 1 6.398e-06 9.837e-06 1.896e-05 0.00% Amr::timeStep() 1 1.232e-05 1.319e-05 1.691e-05 0.00% Castro::post_regrid() 2 1.075e-05 1.163e-05 1.507e-05 0.00% Castro::estTimeStep() 10 1.226e-05 1.288e-05 1.364e-05 0.00% Castro::advance() 1 6.907e-06 8.773e-06 1.358e-05 0.00% Amr::initialInit() 1 2.939e-06 4.526e-06 1.223e-05 0.00% Amr::init() 1 6.639e-06 7.92e-06 1.031e-05 0.00% Amr::bldFineLevels() 1 6.945e-06 7.594e-06 9.09e-06 0.00% Castro::gwstrain() 2 4.162e-06 4.959e-06 8.569e-06 0.00% Castro::errorEst() 2 6.181e-06 6.665e-06 7.397e-06 0.00% BoxList::parallelComplementIn() 4 6.976e-06 7.208e-06 7.377e-06 0.00% Castro::post_timestep() 1 5.715e-06 6.328e-06 7.28e-06 0.00% BoxList::complementIn 4 5.873e-06 6.085e-06 6.462e-06 0.00% Gravity::set_mass_offset() 2 4.35e-06 4.706e-06 5.012e-06 0.00% Castro::FluxRegCrseInit 1 1.542e-06 1.932e-06 2.152e-06 0.00% Castro::FluxRegFineAdd() 1 1.409e-06 1.528e-06 1.737e-06 0.00% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.02 14.02 14.02 100.00% Amr::coarseTimeStep() 1 11.83 11.83 11.83 84.35% Amr::timeStep() 1 11.77 11.77 11.77 83.98% Castro::advance() 1 11.71 11.71 11.72 83.56% Castro::subcycle_advance_ctu() 1 11.66 11.66 11.66 83.19% Castro::do_advance_ctu() 20 11.62 11.63 11.63 82.93% Castro::construct_ctu_hydro_source() 20 6.337 6.368 6.399 45.64% Castro::clean_state() 79 3.706 3.792 3.942 28.11% Castro::computeTemp() 80 2.39 2.495 2.709 19.32% Castro::reset_internal_energy(MultiFab) 80 1.486 1.535 1.677 11.96% Gravity::solve_phi_with_mlmg() 19 1.476 1.523 1.568 11.19% Castro::construct_new_gravity() 17 1.404 1.449 1.495 10.66% Gravity::solve_for_phi() 17 1.365 1.412 1.457 10.39% Castro::initialize_do_advance() 20 1.36 1.396 1.428 10.18% Amr::InitAmr() 1 1.304 1.304 1.305 9.30% Gravity::actual_solve_with_mlmg() 19 1.276 1.28 1.284 9.16% MLMG::solve() 18 1.225 1.23 1.231 8.78% Castro::do_old_sources() 20 1.141 1.168 1.188 8.47% MLMG::oneIter() 128 1.096 1.097 1.099 7.84% MLMG::mgVcycle() 128 1.082 1.083 1.085 7.74% Castro::do_new_sources() 17 0.9406 0.9794 1.017 7.25% Castro::evaluate_source_change() 196 0.9471 0.9552 0.9703 6.92% FabArray::ParallelCopy() 1429 0.2728 0.7037 0.8979 6.40% Amr::init() 1 0.8892 0.89 0.89 6.35% Castro::enforce_min_density() 79 0.7882 0.8043 0.8351 5.96% Amr::initialInit() 1 0.8336 0.8343 0.8344 5.95% amrex::Copy() 724 0.7832 0.7875 0.7985 5.69% MLCellLinOp::smooth() 3072 0.5473 0.5962 0.7232 5.16% FabArray::ParallelCopy_finish() 1429 0.04691 0.4797 0.6751 4.82% AmrLevel::FillPatch() 147 0.4653 0.5442 0.6117 4.36% Amr::InitializeInit() 1 0.4466 0.505 0.6061 4.32% Amr::defBaseLevel() 1 0.4466 0.505 0.6061 4.32% MultiFab::AddProduct() 3439 0.5569 0.5614 0.5719 4.08% Castro::normalize_species() 79 0.546 0.5496 0.5571 3.97% Castro::reset_internal_energy(Fab) 160 0.5174 0.5226 0.532 3.79% FillPatchIterator::Initialize 147 0.2881 0.365 0.4322 3.08% MLMG::mgVcycle_up::1 128 0.09375 0.2558 0.4179 2.98% FillPatchIterator::FillFromLevel0() 147 0.2639 0.3466 0.4156 2.96% FillPatchSingleLevel 147 0.2632 0.3459 0.4149 2.96% MLPoisson::Fsmooth() 6144 0.3142 0.3388 0.413 2.95% Castro::Castro() 1 0.3663 0.3716 0.3921 2.80% Amr::FinalizeInit() 1 0.2283 0.3293 0.3877 2.77% Castro::initMFs() 1 0.3605 0.3658 0.3864 2.76% MLCellLinOp::applyBC() 7148 0.2744 0.305 0.372 2.65% MLMG::addInterpCorrection() 768 0.06232 0.2988 0.3683 2.63% FabArray::setVal() 5013 0.3499 0.3531 0.3661 2.61% MultiFab::Subtract() 159 0.3407 0.3434 0.3546 2.53% Gravity::fill_multipole_BCs() 19 0.1975 0.2419 0.2827 2.02% FabArray::FillBoundary() 6380 0.2161 0.2429 0.2801 2.00% MLMG::mgVcycle_up::2 128 0.01067 0.1133 0.2708 1.93% MLMG::mgVcycle_down::0 128 0.2625 0.2637 0.2658 1.90% Castro::expand_state() 20 0.1765 0.2074 0.2407 1.72% FabArray::ParallelCopy_nowait() 1429 0.2059 0.2204 0.2402 1.71% Castro::initData() 1 0.07769 0.131 0.2373 1.69% check_for_negative_density() 20 0.03436 0.1629 0.2163 1.54% Amr::bldFineLevels() 1 0.04571 0.1469 0.2046 1.46% MLMG::mgVcycle_up::0 128 0.1945 0.1948 0.1953 1.39% MLCellLinOp::apply() 986 0.1482 0.1589 0.1921 1.37% FabArray::sum() 3439 0.1786 0.1802 0.183 1.31% MLMG::mgVcycle_up::3 128 0.0106 0.0233 0.176 1.26% Amr::grid_places() 2 0.002281 0.1033 0.1618 1.15% AmrMesh::MakeNewGrids() 2 0.00216 0.1032 0.1617 1.15% TagBoxArray::mapPRD 2 0.0007164 0.09988 0.1599 1.14% MLMG:computeResOfCorrection() 768 0.1137 0.1231 0.1454 1.04% MLCellLinOp::correctionResidual() 768 0.1123 0.1217 0.1439 1.03% FillBoundary_nowait() 6380 0.08747 0.1137 0.1431 1.02% Castro::post_init() 1 0.1332 0.139 0.1405 1.00% FillBoundary_finish() 6380 0.09901 0.115 0.1383 0.99% Castro::construct_new_source() 102 0.1246 0.1256 0.1275 0.91% MLMG::mgVcycle_down::1 128 0.1132 0.117 0.1226 0.87% MLPoisson::Fapply() 986 0.0946 0.0992 0.1146 0.82% Gravity::multilevel_solve_for_new_phi() 1 0.1009 0.101 0.101 0.72% Gravity::actual_multilevel_solve() 1 0.1009 0.1009 0.1009 0.72% Amr::writePlotFile() 2 0.0987 0.09879 0.0988 0.70% MLCellLinOp::solutionResidual() 147 0.08624 0.09128 0.09622 0.69% Castro::post_regrid() 2 0.08542 0.08695 0.09268 0.66% Castro::construct_old_source() 120 0.07798 0.07836 0.07973 0.57% FabArray::Saxpy() 264 0.07065 0.07162 0.07947 0.57% MLMG::computeResidual() 128 0.07473 0.07528 0.07663 0.55% MLMG::mgVcycle_down::2 128 0.01255 0.04098 0.07229 0.52% MLMG::mgVcycle_bottom 128 0.0003611 0.00482 0.07158 0.51% Castro::construct_new_gravity_source() 17 0.06943 0.0703 0.07141 0.51% MLMG::actualBottomSolve() 128 0.0001322 0.004585 0.07133 0.51% Castro::derive() 18 0.06919 0.06953 0.07091 0.51% AmrLevel::derive() 18 0.06914 0.06949 0.07086 0.51% MLCGSolver::bicgstab 8 0 0.004361 0.06978 0.50% Castro::create_source_corrector() 17 0.04735 0.05531 0.06585 0.47% Castro::post_timestep() 1 0.0586 0.06132 0.06354 0.45% FabArray::Xpay() 945 0.05097 0.05328 0.06079 0.43% MLMG::mgVcycle_down::3 128 0.01191 0.01796 0.05944 0.42% Castro::construct_old_gravity() 20 0.04177 0.05147 0.05685 0.41% Gravity::get_old_grav_vector() 20 0.039 0.04899 0.05607 0.40% amrex::average_down 768 0.03261 0.0391 0.05465 0.39% Castro::initialize_advance() 1 0.04877 0.04974 0.05301 0.38% Amr::regrid() 1 0.0439 0.04467 0.04708 0.34% Castro::check_for_nan() 37 0.04556 0.04598 0.04664 0.33% MultiFab::contains_nan() 37 0.04546 0.04587 0.04652 0.33% Castro::construct_old_gravity_source() 20 0.04477 0.0451 0.04645 0.33% MLPoisson::define() 19 0.02947 0.03396 0.04011 0.29% Gravity::get_new_grav_vector() 18 0.03515 0.03702 0.03804 0.27% Castro::retry_advance_ctu() 11 0.03545 0.0358 0.03699 0.26% MLMG::mgVcycle_down::4 128 0.01185 0.01347 0.03463 0.25% StateDataPhysBCFunct::() 147 0.0007189 0.01778 0.03372 0.24% StateData::FillBoundary(geom) 143 0 0.01692 0.03281 0.23% Castro::apply_source_to_state() 37 0.02872 0.029 0.02934 0.21% MLMG::mgVcycle_up::4 128 0.01064 0.01179 0.02681 0.19% VisMF::Write(FabArray) 4 0.01216 0.01344 0.02568 0.18% MLCellLinOp::defineAuxData() 19 0.01888 0.02031 0.02453 0.17% MLMG::prepareForSolve() 18 0.02002 0.02089 0.02328 0.17% MLMG::mgVcycle_down::5 128 0.01184 0.01269 0.02289 0.16% Castro::sum_integrated_quantities() 2 0.02075 0.02104 0.02106 0.15% Castro::save_data_for_retry() 3 0.01868 0.01893 0.02009 0.14% MLMG::mgVcycle_up::5 128 0.0108 0.01147 0.01947 0.14% MLMG::computeMLResidual() 18 0.01068 0.01554 0.01935 0.14% MLMG::ResNormInf() 146 0.01416 0.01615 0.01739 0.12% FabArray::norminf() 282 0.01076 0.01121 0.01523 0.11% TagBoxArray::collate() 2 8.579e-05 0.002082 0.01478 0.11% amrex::Add() 128 0.0141 0.0142 0.01446 0.10% Castro::problem_post_timestep() 2 0.009415 0.01236 0.01426 0.10% Castro::wd_update() 2 0.009016 0.01232 0.01425 0.10% FabArrayBase::getFB() 6380 0.0102 0.01103 0.01412 0.10% FabArrayBase::getCPC() 2005 0.00948 0.01074 0.01409 0.10% Gravity::test_composite_phi() 1 0.01381 0.01383 0.01385 0.10% MLLinOp::define() 19 0.00384 0.007732 0.01198 0.09% MLLinOp::defineGrids() 19 0.00376 0.007651 0.0119 0.08% FabArray::setDomainBndry() 147 0.0004524 0.005949 0.01167 0.08% MLLinOp::makeSubCommunicator() 19 0.001011 0.004777 0.008877 0.06% Gravity::update_max_rhs() 2 0.001439 0.00717 0.008703 0.06% Amr::writeSmallPlotFile() 2 0.008686 0.008689 0.008694 0.06% MLCGSolver::ParallelAllReduce 227 0 0.0005155 0.008248 0.06% MLMG::getGradSolution() 18 0.007358 0.007647 0.00788 0.06% MLCellLinOp::compGrad() 18 0.007319 0.007606 0.007838 0.06% FabArray::mult() 59 0.006927 0.007066 0.0075 0.05% MultiFab::max() 2 0.0001974 0.005929 0.007473 0.05% Castro::volWgtSum() 46 0.007203 0.007288 0.007419 0.05% FabArrayBase::CPC::define() 576 0.002213 0.003349 0.006322 0.05% amrex::Dot() 142 0 0.000379 0.006064 0.04% MLCellLinOp::defineBC() 19 0.005656 0.005779 0.006008 0.04% Castro::volInBoundary() 14 0.00339 0.003646 0.005385 0.04% BndryData::define() 19 0.004009 0.004217 0.004366 0.03% MLPoisson::normalize() 79 0 0.0001946 0.003113 0.02% MLPoisson::prepareForSolve() 19 0.001749 0.001871 0.002164 0.02% MLCellLinOp::prepareForSolve() 19 0.001657 0.001786 0.002077 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 608 0.001323 0.001544 0.001781 0.01% MLMG::MLRhsNormInf() 18 0.001189 0.001226 0.001573 0.01% AmrLevel::AmrLevel(dm) 1 0.001397 0.001415 0.001444 0.01% FabArrayBase::FB::FB() 75 0.000272 0.0006713 0.001426 0.01% StateData::define() 6 0.001339 0.001358 0.001386 0.01% MLMG::MLResNormInf() 18 0.001223 0.001249 0.001377 0.01% Gravity::swapTimeLevels() 8 0.001144 0.001197 0.001281 0.01% MLMG::compResidual() 1 0.0009422 0.001107 0.001268 0.01% Castro::locWgtSum() 6 0.001106 0.00111 0.001127 0.01% MLLinOp::makeAgglomeratedDMap 19 0.0009155 0.0009746 0.001085 0.01% DistributionMapping::KnapSackDoIt() 1 9.556e-05 0.0008043 0.0008723 0.01% MLCellLinOp::setLevelBC() 19 0.0005536 0.0006744 0.0008649 0.01% DistributionMapping::LeastUsedCPUs() 1 2.198e-05 0.0007269 0.0007911 0.01% makeSFC 95 0.0005874 0.00063 0.0007117 0.01% Castro::buildMetrics() 1 0.0006383 0.0006563 0.0006883 0.00% Castro::errorEst() 2 0.0005131 0.0005823 0.0006631 0.00% Castro::apply_tagging_restrictions() 2 0.0004145 0.0004821 0.00056 0.00% Gravity::applyMetricTerms() 2 0.0004871 0.0004994 0.0005334 0.00% Castro::swap_state_time_levels() 8 0.0004187 0.000434 0.0004478 0.00% Castro::volProductSum() 2 0.0003843 0.0003874 0.000391 0.00% Castro::enforce_consistent_e() 1 0.0003628 0.0003676 0.0003806 0.00% Castro::finalize_advance() 1 0.0002548 0.0002674 0.0003368 0.00% DistributionMapping::Distribute() 95 0.0002138 0.0002331 0.0002863 0.00% Castro::enforce_speed_limit() 79 0.0001082 0.00012 0.000212 0.00% Castro::finalize_do_advance() 17 0.0001666 0.0001799 0.0002015 0.00% OwnerMask() 2 0.0001338 0.0001491 0.0001644 0.00% Amr::initSubcycle() 1 5.536e-05 6.326e-05 0.000153 0.00% Castro::apply_problem_tags() 2 9.184e-05 9.353e-05 0.0001009 0.00% MultiFab::Swap() 5 6.448e-05 6.724e-05 8.354e-05 0.00% knapsack() 1 5.527e-05 5.786e-05 6.16e-05 0.00% Castro::construct_old_rotation_source() 20 2.997e-05 3.077e-05 3.372e-05 0.00% Castro::computeInitialDt() 2 9.638e-06 1.151e-05 2.752e-05 0.00% Castro::construct_new_rotation_source() 17 2.132e-05 2.19e-05 2.372e-05 0.00% BoxList::parallelComplementIn() 4 1.287e-05 1.329e-05 1.384e-05 0.00% Castro::estTimeStep() 10 1.226e-05 1.288e-05 1.364e-05 0.00% Castro::gwstrain() 2 4.162e-06 4.959e-06 8.569e-06 0.00% BoxList::complementIn 4 5.873e-06 6.085e-06 6.462e-06 0.00% Gravity::set_mass_offset() 2 4.35e-06 4.706e-06 5.012e-06 0.00% Castro::FluxRegCrseInit 1 1.542e-06 1.932e-06 2.152e-06 0.00% Castro::FluxRegFineAdd() 1 1.409e-06 1.528e-06 1.737e-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 14 B 17 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.03-13-g3fe7aadc23c9) finalized