Initializing AMReX (24.10-18-g92679babfc2c)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.10-18-g92679babfc2c) initialized Starting run at 07:32:14 UTC on 2024-10-24. Successfully read inputs file ... Castro git describe: 24.10-7-g1686df30a AMReX git describe: 24.10-18-g92679babf Microphysics git describe: 24.10 reading extern runtime parameters ... 13 Species: He4 C12 O16 Ne20 Mg24 Si28 S32 Ar36 Ca40 Ti44 Cr48 Fe52 Ni56 Successfully read inputs file ... Creating primary with CO core with mass fractions C = 0.5and O = 0.5 and a He shell of solar mass =0. Creating secondary with CO core with mass fractions C = 0.5and O = 0.5 and a He shell of solar mass =0. 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 Now regridding at level lbase = 0 grid_places() time: 0.181037362 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001226192 new finest: 0 Gravity::fill_multipole_BCs() time = 0.013255983 Gravity::multilevel_solve_for_new_phi() time = 0.13126334 ... norm of composite residual at level 0 1.193955541e-11 Gravity::fill_multipole_BCs() time = 0.010988656 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.051864176 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.007225143 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. ... Entering construct_ctu_hydro_source() on level 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 = -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 Invalid density = -61.58489855 at index 19, 178, 0 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, Invalid density = -80.16477042 at index 18, 179, 0 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 = -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 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 Castro::construct_old_gravity() time = 0.0034267 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.067644221 on level 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. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Castro::construct_old_gravity() time = 0.040857703 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.12419156 on level 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. ... 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 Castro::construct_old_gravity() time = 0.076838003 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.108248769 on level 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. ... Entering construct_ctu_hydro_source() on level 0 Castro::construct_old_gravity() time = 0.005301524 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.08307648 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at level 0 Castro::construct_ctu_hydro_source() time = 0.360249197 on level 0 Ending SDC iteration 1 of 3. Beginning SDC iteration 2 of 3. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.019607913 Gravity::solve_for_phi() time = 0.101443308 on level 0 Castro::construct_new_gravity() time = 0.103667982 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.0623056 on level 0 Castro::construct_old_gravity() time = 0.005685132 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.085212796 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at level 0 Castro::construct_ctu_hydro_source() time = 0.352160554 on level 0 Ending SDC iteration 2 of 3. Beginning SDC iteration 3 of 3. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016916463 Gravity::solve_for_phi() time = 0.099618983 on level 0 Castro::construct_new_gravity() time = 0.101875374 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.062455951 on level 0 Castro::construct_old_gravity() time = 0.004737905 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.078012763 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at level 0 Castro::construct_ctu_hydro_source() time = 0.437920287 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. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.066023326 Gravity::solve_for_phi() time = 0.197814457 on level 0 Castro::construct_new_gravity() time = 0.201260513 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.110800008 on level 0 Castro::construct_old_gravity() time = 0.046648182 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.110660635 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.985562121e+46 ener added: 2.745854866e+42 Castro::construct_ctu_hydro_source() time = 0.658053516 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.066879206 Gravity::solve_for_phi() time = 0.176593975 on level 0 Castro::construct_new_gravity() time = 0.178907516 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.270502209e+44 ener added: 4.891857977e+43 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.763276272e+35 ymom added: -1.134784666e+25 zmom added: 0 eint added: 0 ener added: 6.004315505e+43 Castro::do_new_sources() time = 0.077821091 on level 0 Castro::construct_old_gravity() time = 0.017582796 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.066429159 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.986261807e+46 ener added: 2.752527215e+42 Castro::construct_ctu_hydro_source() time = 0.375182275 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. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.022185134 Gravity::solve_for_phi() time = 0.105359426 on level 0 Castro::construct_new_gravity() time = 0.107587679 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.344002121e+44 ener added: 4.888146542e+43 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.757171236e+35 ymom added: -7.092994457e+24 zmom added: 0 eint added: 0 ener added: 5.273131065e+43 Castro::do_new_sources() time = 0.06937159 on level 0 Castro::construct_old_gravity() time = 0.014780433 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.984353216e+46 ener added: 6.990361107e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467741735e+40 ymom added: 9.47571169e+26 zmom added: 0 eint added: 0 ener added: 8.780735301e+46 Castro::do_old_sources() time = 0.065397236 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.987305793e+46 ener added: 2.180540635e+43 Castro::construct_ctu_hydro_source() time = 0.39836977 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.057294144 Gravity::solve_for_phi() time = 0.15420261 on level 0 Castro::construct_new_gravity() time = 0.156458248 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.032152963e+43 ener added: 1.18439637e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.781732367e+35 ymom added: -1.228995877e+24 zmom added: 0 eint added: 0 ener added: 5.994406477e+43 Castro::do_new_sources() time = 0.06229 on level 0 Castro::construct_old_gravity() time = 0.003706424 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.984353216e+46 ener added: 6.990361107e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467741735e+40 ymom added: 9.47571169e+26 zmom added: 0 eint added: 0 ener added: 8.780735301e+46 Castro::do_old_sources() time = 0.061294426 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.988012461e+46 ener added: 2.185486484e+43 Castro::construct_ctu_hydro_source() time = 0.440891869 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. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.015756471 Gravity::solve_for_phi() time = 0.096951654 on level 0 Castro::construct_new_gravity() time = 0.099215785 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.774956503e+43 ener added: 1.18358902e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.775605302e+35 ymom added: 2.895991246e+24 zmom added: 0 eint added: 0 ener added: 5.259726767e+43 Castro::do_new_sources() time = 0.075445666 on level 0 Castro::construct_old_gravity() time = 0.004982816 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.983917951e+46 ener added: 1.154265389e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467777247e+40 ymom added: 9.523313144e+26 zmom added: 0 eint added: 0 ener added: 8.841018523e+46 Castro::do_old_sources() time = 0.092486836 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.984994198e+46 ener added: 2.292103419e+43 Castro::construct_ctu_hydro_source() time = 0.35678405 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.015529711 Gravity::solve_for_phi() time = 0.097324052 on level 0 Castro::construct_new_gravity() time = 0.09955781 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.621864779e+43 ener added: 1.488585404e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.801468831e+35 ymom added: -8.889854904e+23 zmom added: 0 eint added: 0 ener added: 4.63003051e+43 Castro::do_new_sources() time = 0.061632418 on level 0 Castro::construct_old_gravity() time = 0.002892722 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.983917951e+46 ener added: 1.154265389e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467777247e+40 ymom added: 9.523313144e+26 zmom added: 0 eint added: 0 ener added: 8.841018523e+46 Castro::do_old_sources() time = 0.061498848 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.985704698e+46 ener added: 2.298038423e+43 Castro::construct_ctu_hydro_source() time = 0.362036956 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. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.015164157 Gravity::solve_for_phi() time = 0.09651178 on level 0 Castro::construct_new_gravity() time = 0.099139988 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.366917123e+43 ener added: 1.48767998e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.795324993e+35 ymom added: 4.258393976e+24 zmom added: 0 eint added: 0 ener added: 3.89401763e+43 Castro::do_new_sources() time = 0.061845569 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.44944741e+29 ener added: 1.44944741e+29 Castro::construct_old_gravity() time = 0.006826912 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.983565234e+46 ener added: 1.521217471e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467813153e+40 ymom added: 9.596604272e+26 zmom added: 0 eint added: 0 ener added: 8.902573412e+46 Castro::do_old_sources() time = 0.061604652 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.992395521e+46 ener added: 2.321353363e+43 Castro::construct_ctu_hydro_source() time = 0.356256379 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.01655995 Gravity::solve_for_phi() time = 0.097478925 on level 0 Castro::construct_new_gravity() time = 0.099763615 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.10309639e+43 ener added: 1.620771775e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.826334308e+35 ymom added: -1.563103306e+24 zmom added: 0 eint added: 0 ener added: 2.089132415e+43 Castro::do_new_sources() time = 0.063403254 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.44944741e+29 ener added: 1.44944741e+29 Castro::construct_old_gravity() time = 0.003684073 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.983565234e+46 ener added: 1.521217471e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467813153e+40 ymom added: 9.596604272e+26 zmom added: 0 eint added: 0 ener added: 8.902573412e+46 Castro::do_old_sources() time = 0.062712208 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.993115374e+46 ener added: 2.326787816e+43 Castro::construct_ctu_hydro_source() time = 0.360846848 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. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.015858845 Gravity::solve_for_phi() time = 0.10541441 on level 0 Castro::construct_new_gravity() time = 0.107662289 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.855214541e+43 ener added: 1.619926565e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.820160807e+35 ymom added: -5.840386748e+24 zmom added: 0 eint added: 0 ener added: 1.345414315e+43 Castro::do_new_sources() time = 0.069700297 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.414226202e+29 ener added: 4.414226202e+29 Castro::construct_old_gravity() time = 0.00345473 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.983073451e+46 ener added: 2.019355568e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467849557e+40 ymom added: 9.476467269e+26 zmom added: 0 eint added: 0 ener added: 8.96491088e+46 Castro::do_old_sources() time = 0.089909578 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.000649661e+46 ener added: 2.406090033e+43 Castro::construct_ctu_hydro_source() time = 0.373818357 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016098623 Gravity::solve_for_phi() time = 0.097053559 on level 0 Castro::construct_new_gravity() time = 0.099281269 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.896789717e+43 ener added: 1.713953472e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.844643379e+35 ymom added: 5.511001686e+24 zmom added: 0 eint added: 0 ener added: -1.409065908e+43 Castro::do_new_sources() time = 0.060832919 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.414226202e+29 ener added: 4.414226202e+29 Castro::construct_old_gravity() time = 0.003050537 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.983073451e+46 ener added: 2.019355568e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467849557e+40 ymom added: 9.476467269e+26 zmom added: 0 eint added: 0 ener added: 8.96491088e+46 Castro::do_old_sources() time = 0.062981818 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.001373749e+46 ener added: 2.412535663e+43 Castro::construct_ctu_hydro_source() time = 0.36654504 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. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016199724 Gravity::solve_for_phi() time = 0.117223084 on level 0 Castro::construct_new_gravity() time = 0.120132698 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.652630993e+43 ener added: 1.713047936e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.838440231e+35 ymom added: 7.569953472e+24 zmom added: 0 eint added: 0 ener added: -2.155909539e+43 Castro::do_new_sources() time = 0.061642331 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.628569462e+29 ener added: 1.628569462e+29 Castro::construct_old_gravity() time = 0.012720758 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.982914011e+46 ener added: 2.174745257e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467886325e+40 ymom added: 9.629849732e+26 zmom added: 0 eint added: 0 ener added: 9.027797154e+46 Castro::do_old_sources() time = 0.104586646 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.993848698e+46 ener added: 2.676124277e+43 Castro::construct_ctu_hydro_source() time = 0.667062865 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.067189116 Gravity::solve_for_phi() time = 0.184697901 on level 0 Castro::construct_new_gravity() time = 0.188038508 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.152558515e+44 ener added: 1.743203049e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.874363167e+35 ymom added: -1.093700077e+25 zmom added: 0 eint added: 0 ener added: -4.046567664e+43 Castro::do_new_sources() time = 0.097701609 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.628569462e+29 ener added: 1.628569462e+29 Castro::construct_old_gravity() time = 0.003080377 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.982914011e+46 ener added: 2.174745257e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467886325e+40 ymom added: 9.629849732e+26 zmom added: 0 eint added: 0 ener added: 9.027797154e+46 Castro::do_old_sources() time = 0.097269254 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.994580994e+46 ener added: 2.682637513e+43 Castro::construct_ctu_hydro_source() time = 0.542926415 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. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.021344834 Gravity::solve_for_phi() time = 0.103742915 on level 0 Castro::construct_new_gravity() time = 0.106007361 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.228943912e+44 ener added: 1.742301947e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.868126277e+35 ymom added: -1.191216945e+25 zmom added: 0 eint added: 0 ener added: -4.801468886e+43 Castro::do_new_sources() time = 0.061413156 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.927278364e+29 ener added: 8.927278364e+29 Castro::construct_old_gravity() time = 0.002982292 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.982547542e+46 ener added: 2.529323823e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467923688e+40 ymom added: 9.382775517e+26 zmom added: 0 eint added: 0 ener added: 9.091359174e+46 Castro::do_old_sources() time = 0.061479894 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.99339064e+46 ener added: 2.937271748e+43 Castro::construct_ctu_hydro_source() time = 0.369332501 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.0182184 Gravity::solve_for_phi() time = 0.105962745 on level 0 Castro::construct_new_gravity() time = 0.108219387 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.568306138e+43 ener added: 1.766182377e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.939899642e+35 ymom added: 1.818111096e+23 zmom added: 0 eint added: 0 ener added: -1.062741416e+43 Castro::do_new_sources() time = 0.065577398 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.927278364e+29 ener added: 8.927278364e+29 Castro::construct_old_gravity() time = 0.005549414 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.982547542e+46 ener added: 2.529323823e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467923688e+40 ymom added: 9.382775517e+26 zmom added: 0 eint added: 0 ener added: 9.091359174e+46 Castro::do_old_sources() time = 0.065112309 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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.994123755e+46 ener added: 2.943724162e+43 Castro::construct_ctu_hydro_source() time = 0.375828361 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002189731144 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.018392019 Gravity::solve_for_phi() time = 0.105137893 on level 0 Castro::construct_new_gravity() time = 0.107383244 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.333668478e+43 ener added: 1.765299431e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.933664355e+35 ymom added: 1.216009369e+24 zmom added: 0 eint added: 0 ener added: -1.819323893e+43 Castro::do_new_sources() time = 0.063635979 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: 2.042403168e+29 ener added: 2.042403168e+29 TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -4.341806745e+38 TIME= 0.25 YMOM = 1.137357411e+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.468589243e+38 TIME= 0.25 RHO*e = 1.325894515e+50 TIME= 0.25 RHO*K = 3.170941989e+49 TIME= 0.25 RHO*E = 1.642988714e+50 TIME= 0.25 RHO*PHI = -5.449011382e+50 TIME= 0.25 TOTAL ENERGY= -1.081516977e+50 TIME= 0.25 CENTER OF MASS X-LOC = 316424388.6 TIME= 0.25 CENTER OF MASS X-VEL = -170267.1603 TIME= 0.25 CENTER OF MASS Y-LOC = -3.641005016e-07 TIME= 0.25 CENTER OF MASS Y-VEL = 4.460231145e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 58251071.63 TIME= 0.25 MAXIMUM DENSITY = 4458761.874 TIME= 0.25 MAXIMUM T_S / T_E = 0 STEP = 1 TIME = 0.25 DT = 0.25 [STEP 1] Coarse TimeStep time: 15.03339185 [STEP 1] FAB kilobyte spread across MPI nodes: [6973 ... 6973] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.0527232 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004909212 seconds Ending run at 07:32:32 UTC on 2024-10-24. Run time = 18.23084017 Run time without initialization = 15.0920454 Average number of zones advanced per microsecond: 0.002 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 18.2308 ... 18.2308 ... 18.2309 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.9521 6.1969 6.4579 35.42% Amr::InitAmr() 1 2.183 2.183 2.183 11.97% FabArray::ParallelCopy_finish() 1429 0.5848 1.133 1.461 8.01% check_for_negative_density() 20 0.8345 1.141 1.392 7.64% Castro::computeTemp() 80 0.8803 1.016 1.234 6.77% amrex::Copy() 724 0.837 0.8646 0.9292 5.10% MultiFab::AddProduct() 3439 0.5726 0.5965 0.6954 3.81% Castro::normalize_species() 79 0.5756 0.5979 0.6493 3.56% Castro::reset_internal_energy(Fab) 160 0.5398 0.5603 0.6031 3.31% Gravity::fill_multipole_BCs() 19 0.3333 0.4536 0.4999 2.74% MLPoisson::Fsmooth() 6144 0.3619 0.4004 0.4723 2.59% FabArray::setVal() 5013 0.373 0.3842 0.4041 2.22% Castro::initMFs() 1 0.3653 0.3718 0.3884 2.13% MultiFab::Subtract() 159 0.3493 0.3594 0.3828 2.10% FillBoundary_finish() 6380 0.2274 0.2775 0.342 1.88% FabArray::ParallelCopy_nowait() 1429 0.2248 0.2448 0.2684 1.47% FabArray::sum() 3439 0.1913 0.1997 0.217 1.19% Castro::initData() 1 0.003407 0.02924 0.1731 0.95% FillBoundary_nowait() 6380 0.09961 0.1289 0.173 0.95% TagBoxArray::collate() 2 0.0001859 0.02973 0.1549 0.85% MLCellLinOp::applyBC() 7148 0.06667 0.0839 0.1296 0.71% MLPoisson::Fapply() 986 0.09673 0.1047 0.1199 0.66% FabArray::Saxpy() 264 0.07129 0.07491 0.09363 0.51% Castro::construct_new_gravity_source() 17 0.07059 0.07294 0.07821 0.43% FabArray::Xpay() 945 0.05357 0.0572 0.06414 0.35% Castro::construct_old_gravity_source() 20 0.04764 0.05134 0.05879 0.32% MultiFab::contains_nan() 37 0.04749 0.04956 0.05282 0.29% MLMG::addInterpCorrection() 768 0.03452 0.03787 0.04477 0.25% AmrLevel::derive() 18 0.04155 0.04242 0.0442 0.24% StateData::FillBoundary(geom) 143 0 0.01892 0.03618 0.20% amrex::average_down 768 0.026 0.02871 0.03501 0.19% MLMG::ResNormInf() 146 0.02209 0.02766 0.03468 0.19% Castro::reset_internal_energy(MultiFab) 80 0.02486 0.02567 0.02709 0.15% VisMF::Write(FabArray) 4 0.01584 0.01693 0.02698 0.15% Castro::enforce_min_density() 79 0.0198 0.02031 0.02148 0.12% Amr::writePlotFile() 2 0.009701 0.01897 0.02107 0.12% MLCellLinOp::defineAuxData() 19 0.01638 0.01776 0.02086 0.11% Castro::wd_update() 2 0.00508 0.01379 0.01883 0.10% MLMG::solve() 18 0.005966 0.0105 0.0185 0.10% MLCGSolver::bicgstab 8 0 0.001132 0.01812 0.10% Castro::construct_new_source() 102 0.01632 0.0167 0.01773 0.10% FabArray::FillBoundary() 6380 0.01423 0.01493 0.01759 0.10% Castro::evaluate_source_change() 196 0.01459 0.01522 0.01702 0.09% FabArray::setDomainBndry() 147 0.0004967 0.008482 0.01656 0.09% amrex::Add() 128 0.01488 0.01523 0.01584 0.09% FabArray::norminf() 282 0.01089 0.01173 0.01559 0.09% FillPatchIterator::Initialize 147 0.01392 0.01453 0.01531 0.08% FabArrayBase::getFB() 6380 0.01104 0.01158 0.01385 0.08% MLCellLinOp::smooth() 3072 0.01189 0.0126 0.01366 0.07% Castro::construct_old_source() 120 0.01068 0.01136 0.01285 0.07% MLMG::prepareForSolve() 18 0.01064 0.01125 0.0128 0.07% Gravity::actual_solve_with_mlmg() 19 0.008032 0.008391 0.009743 0.05% FabArrayBase::getCPC() 2005 0.008001 0.008293 0.008816 0.05% Gravity::get_old_grav_vector() 20 0.006793 0.007222 0.008062 0.04% Gravity::get_new_grav_vector() 18 0.006216 0.006375 0.006819 0.04% Castro::initialize_do_advance() 20 0.004665 0.005166 0.00609 0.03% MLCellLinOp::apply() 986 0.002089 0.002335 0.004721 0.03% FabArray::ParallelCopy() 1429 0.003641 0.003812 0.004132 0.02% Castro::sum_integrated_quantities() 2 0.002183 0.003629 0.003919 0.02% AmrLevel::FillPatch() 147 0.002668 0.002834 0.003152 0.02% Castro::save_data_for_retry() 3 0.001819 0.001993 0.002304 0.01% MLCellLinOp::correctionResidual() 768 0.001933 0.002016 0.002142 0.01% MLMG::mgVcycle() 128 0.001695 0.001767 0.00192 0.01% MLMG:computeResOfCorrection() 768 0.001502 0.001563 0.001676 0.01% MLMG::actualBottomSolve() 128 0.0001411 0.0002112 0.001149 0.01% StateDataPhysBCFunct::() 147 0.0008607 0.001027 0.001148 0.01% Castro::subcycle_advance_ctu() 1 0.0002836 0.0003855 0.001049 0.01% Gravity::solve_for_phi() 17 0.0007495 0.000807 0.000895 0.00% main() 1 0.0007025 0.0007863 0.0008761 0.00% FillPatchIterator::FillFromLevel0() 147 0.0007382 0.0007864 0.0008407 0.00% MLMG::mgVcycle_down::3 128 0.0005957 0.0006563 0.0007976 0.00% MLMG::mgVcycle_down::2 128 0.0006578 0.0007146 0.0007949 0.00% MLMG::mgVcycle_down::0 128 0.0007006 0.0007256 0.0007738 0.00% MLMG::mgVcycle_down::1 128 0.0006862 0.0007138 0.0007623 0.00% Castro::Castro() 1 0.0003781 0.0004369 0.0007288 0.00% MLMG::mgVcycle_down::4 128 0.0005774 0.0006085 0.0007234 0.00% MLMG::mgVcycle_down::5 128 0.0006005 0.000628 0.000711 0.00% AmrMesh::MakeNewGrids() 2 0.0005825 0.0005966 0.0006216 0.00% Castro::clean_state() 79 0.0005204 0.0005665 0.0006173 0.00% Castro::do_advance_ctu() 20 0.0004396 0.0004927 0.0005853 0.00% FillPatchSingleLevel 147 0.0005047 0.0005374 0.0005769 0.00% MLMG::mgVcycle_up::0 128 0.000494 0.0005238 0.0005612 0.00% MLMG::mgVcycle_up::1 128 0.0004969 0.0005214 0.0005611 0.00% MLMG::mgVcycle_up::2 128 0.000416 0.0004668 0.0005327 0.00% MLCellLinOp::solutionResidual() 147 0.0004698 0.0004866 0.0005103 0.00% MLMG::mgVcycle_up::3 128 0.0004112 0.0004335 0.0005043 0.00% MLMG::mgVcycle_up::4 128 0.0004156 0.0004348 0.0004758 0.00% MLMG::mgVcycle_up::5 128 0.0004142 0.0004314 0.0004736 0.00% MLMG::oneIter() 128 0.0004198 0.0004367 0.0004725 0.00% Castro::do_old_sources() 20 0.000313 0.0003461 0.0004198 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002515 0.0002882 0.0003441 0.00% Castro::construct_new_gravity() 17 0.000178 0.0002013 0.0003154 0.00% Castro::problem_post_timestep() 2 1.304e-05 3.268e-05 0.0002777 0.00% Castro::do_new_sources() 17 0.0002292 0.0002509 0.0002766 0.00% MLMG::computeResidual() 128 0.0002465 0.0002555 0.0002688 0.00% MLMG::mgVcycle_bottom 128 0.0002324 0.0002453 0.0002575 0.00% Amr::coarseTimeStep() 1 0.0001746 0.0002131 0.0002203 0.00% Castro::check_for_nan() 37 0.0001297 0.0001443 0.0001702 0.00% Castro::initialize_advance() 1 0.0001325 0.000149 0.000161 0.00% Amr::defBaseLevel() 1 0.0001106 0.0001282 0.0001506 0.00% Castro::construct_old_gravity() 20 0.0001157 0.0001307 0.0001485 0.00% MLPoisson::define() 19 0.0001256 0.0001353 0.0001475 0.00% Castro::retry_advance_ctu() 11 7.426e-05 9.39e-05 0.0001415 0.00% Castro::apply_source_to_state() 37 0.0001111 0.0001237 0.0001415 0.00% TagBoxArray::mapPRD 2 0.0001254 0.0001329 0.0001405 0.00% Gravity::test_composite_phi() 1 0.0001108 0.0001162 0.000131 0.00% Castro::create_source_corrector() 17 6.564e-05 7.731e-05 9.27e-05 0.00% Amr::grid_places() 2 5.463e-05 6.027e-05 8.526e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.247e-05 6.726e-05 7.674e-05 0.00% Castro::expand_state() 20 5.777e-05 6.36e-05 7.105e-05 0.00% Amr::regrid() 1 1.646e-05 2.217e-05 6.028e-05 0.00% Castro::derive() 18 4.411e-05 4.898e-05 5.392e-05 0.00% Amr::FinalizeInit() 1 8.757e-06 1.332e-05 5.314e-05 0.00% MLMG::computeMLResidual() 18 4.225e-05 4.752e-05 5.294e-05 0.00% Gravity::actual_multilevel_solve() 1 3.279e-05 3.508e-05 3.824e-05 0.00% Amr::timeStep() 1 1.553e-05 1.964e-05 2.944e-05 0.00% Castro::advance() 1 7.968e-06 1.153e-05 2.732e-05 0.00% Castro::post_regrid() 2 1.083e-05 1.318e-05 1.925e-05 0.00% Castro::post_init() 1 1.215e-05 1.45e-05 1.758e-05 0.00% Amr::InitializeInit() 1 7.694e-06 8.913e-06 1.349e-05 0.00% Amr::init() 1 6.413e-06 7.993e-06 1.116e-05 0.00% Castro::post_timestep() 1 7.277e-06 8.81e-06 1.085e-05 0.00% Amr::bldFineLevels() 1 7.394e-06 8.445e-06 1.032e-05 0.00% Amr::initialInit() 1 3.088e-06 3.35e-06 3.863e-06 0.00% Other 2478 0.05451 0.07095 0.1008 0.55% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 18.23 18.23 18.23 100.00% Amr::coarseTimeStep() 1 15.09 15.09 15.09 82.78% Amr::timeStep() 1 15.03 15.03 15.03 82.46% Castro::advance() 1 14.96 14.96 14.97 82.10% Castro::subcycle_advance_ctu() 1 14.91 14.91 14.91 81.81% Castro::do_advance_ctu() 20 14.86 14.87 14.87 81.57% Castro::construct_ctu_hydro_source() 20 8.122 8.17 8.265 45.33% Castro::clean_state() 79 3.835 4.054 4.487 24.61% Castro::computeTemp() 80 2.45 2.632 2.957 16.22% Amr::InitAmr() 1 2.183 2.183 2.183 11.97% Gravity::solve_phi_with_mlmg() 19 2.008 2.129 2.175 11.93% Castro::construct_new_gravity() 17 1.908 2.029 2.074 11.38% Gravity::solve_for_phi() 17 1.866 1.987 2.033 11.15% Castro::initialize_do_advance() 20 1.608 1.694 1.78 9.76% FabArray::ParallelCopy() 1429 0.8792 1.403 1.724 9.45% Castro::reset_internal_energy(MultiFab) 80 1.56 1.616 1.723 9.45% Gravity::actual_solve_with_mlmg() 19 1.668 1.672 1.676 9.19% MLMG::solve() 18 1.609 1.615 1.617 8.87% Castro::do_old_sources() 20 1.427 1.51 1.564 8.58% FabArray::ParallelCopy_finish() 1429 0.5848 1.133 1.461 8.01% MLMG::oneIter() 128 1.416 1.421 1.427 7.83% MLMG::mgVcycle() 128 1.4 1.406 1.411 7.74% check_for_negative_density() 20 0.8345 1.141 1.392 7.64% AmrLevel::FillPatch() 147 0.8631 1.191 1.358 7.45% FillPatchIterator::Initialize 147 0.6483 0.993 1.168 6.40% FillPatchIterator::FillFromLevel0() 147 0.621 0.97 1.147 6.29% FillPatchSingleLevel 147 0.6201 0.9692 1.146 6.29% Castro::evaluate_source_change() 196 0.9894 1.029 1.131 6.20% Castro::do_new_sources() 17 0.9644 1.018 1.096 6.01% MLCellLinOp::smooth() 3072 0.7598 0.8253 0.9839 5.40% Amr::init() 1 0.9557 0.9557 0.9557 5.24% amrex::Copy() 724 0.837 0.8646 0.9292 5.10% Castro::enforce_min_density() 79 0.8238 0.8521 0.908 4.98% Amr::initialInit() 1 0.8961 0.8962 0.8962 4.92% MultiFab::AddProduct() 3439 0.5726 0.5965 0.6954 3.81% Castro::normalize_species() 79 0.5756 0.5979 0.6493 3.56% Amr::InitializeInit() 1 0.4532 0.4864 0.6326 3.47% Amr::defBaseLevel() 1 0.4532 0.4864 0.6326 3.47% MLCellLinOp::applyBC() 7148 0.4461 0.5176 0.6286 3.45% Castro::reset_internal_energy(Fab) 160 0.5398 0.5603 0.6031 3.31% MLMG::mgVcycle_up::1 128 0.1285 0.3283 0.5314 2.92% Gravity::fill_multipole_BCs() 19 0.3358 0.4559 0.5021 2.75% FabArray::FillBoundary() 6380 0.3627 0.4337 0.499 2.74% MLPoisson::Fsmooth() 6144 0.3619 0.4004 0.4723 2.59% MLMG::addInterpCorrection() 768 0.06833 0.3582 0.4508 2.47% Amr::FinalizeInit() 1 0.2635 0.4098 0.443 2.43% Castro::expand_state() 20 0.2432 0.3497 0.4266 2.34% FabArray::setVal() 5013 0.373 0.3842 0.4041 2.22% Castro::Castro() 1 0.3709 0.3774 0.394 2.16% Castro::initMFs() 1 0.3655 0.372 0.3886 2.13% MultiFab::Subtract() 159 0.3493 0.3594 0.3828 2.10% FillBoundary_finish() 6380 0.2274 0.2775 0.342 1.88% MLMG::mgVcycle_down::0 128 0.3238 0.3308 0.3369 1.85% MLMG::mgVcycle_up::2 128 0.01418 0.1404 0.3324 1.82% FabArray::ParallelCopy_nowait() 1429 0.2456 0.2662 0.2908 1.60% MLMG::mgVcycle_up::0 128 0.2536 0.2559 0.2581 1.42% Castro::initData() 1 0.08036 0.1071 0.2515 1.38% MLCellLinOp::apply() 986 0.1952 0.2102 0.2445 1.34% Amr::bldFineLevels() 1 0.04766 0.1938 0.2267 1.24% FabArray::sum() 3439 0.1913 0.1997 0.217 1.19% Castro::construct_old_gravity() 20 0.07608 0.1558 0.2042 1.12% Gravity::get_old_grav_vector() 20 0.07595 0.1557 0.2041 1.12% MLMG::mgVcycle_up::3 128 0.0141 0.02889 0.2007 1.10% Castro::create_source_corrector() 17 0.1032 0.1447 0.1914 1.05% FillBoundary_nowait() 6380 0.1115 0.1413 0.1865 1.02% Amr::grid_places() 2 0.002848 0.1491 0.1822 1.00% AmrMesh::MakeNewGrids() 2 0.002793 0.149 0.1822 1.00% MLMG::mgVcycle_down::1 128 0.1646 0.1705 0.1805 0.99% TagBoxArray::mapPRD 2 0.001153 0.1179 0.1805 0.99% MLMG:computeResOfCorrection() 768 0.1352 0.1485 0.174 0.95% MLCellLinOp::correctionResidual() 768 0.1337 0.147 0.1724 0.95% Castro::post_init() 1 0.1653 0.171 0.1722 0.94% TagBoxArray::collate() 2 0.0001859 0.02973 0.1549 0.85% Castro::construct_new_source() 102 0.1277 0.1313 0.1437 0.79% Gravity::multilevel_solve_for_new_phi() 1 0.1312 0.1312 0.1313 0.72% Gravity::actual_multilevel_solve() 1 0.1311 0.1311 0.1312 0.72% MLCellLinOp::solutionResidual() 147 0.1082 0.1211 0.1283 0.70% MLPoisson::Fapply() 986 0.09673 0.1047 0.1199 0.66% Amr::writePlotFile() 2 0.105 0.105 0.1051 0.58% Castro::construct_old_source() 120 0.08284 0.08811 0.1026 0.56% MLMG::computeResidual() 128 0.09121 0.094 0.09688 0.53% MLMG::mgVcycle_down::2 128 0.01652 0.05429 0.09671 0.53% Castro::post_regrid() 2 0.08847 0.08967 0.09536 0.52% FabArray::Saxpy() 264 0.07129 0.07491 0.09363 0.51% MLMG::mgVcycle_bottom 128 0.0003735 0.005287 0.0787 0.43% MLMG::actualBottomSolve() 128 0.0001411 0.005042 0.07844 0.43% Castro::construct_new_gravity_source() 17 0.07059 0.07294 0.07821 0.43% MLMG::mgVcycle_down::3 128 0.01578 0.0238 0.07698 0.42% MLCGSolver::bicgstab 8 0 0.004799 0.07679 0.42% Castro::derive() 18 0.07198 0.0732 0.07492 0.41% AmrLevel::derive() 18 0.07193 0.07315 0.07487 0.41% Castro::post_timestep() 1 0.06619 0.07065 0.07432 0.41% amrex::average_down 768 0.03534 0.04546 0.07045 0.39% FabArray::Xpay() 945 0.05357 0.0572 0.06414 0.35% Castro::construct_old_gravity_source() 20 0.04764 0.05134 0.05879 0.32% Castro::initialize_advance() 1 0.05109 0.05171 0.05314 0.29% Castro::check_for_nan() 37 0.04762 0.0497 0.05296 0.29% MultiFab::contains_nan() 37 0.04749 0.04956 0.05282 0.29% MLPoisson::define() 19 0.03387 0.03897 0.04824 0.26% Castro::retry_advance_ctu() 11 0.03959 0.04207 0.04753 0.26% Amr::regrid() 1 0.04553 0.04593 0.04715 0.26% MLMG::ResNormInf() 146 0.03217 0.03762 0.04427 0.24% Gravity::get_new_grav_vector() 18 0.04004 0.04096 0.04199 0.23% MLMG::mgVcycle_down::4 128 0.01562 0.01774 0.04107 0.23% Castro::apply_source_to_state() 37 0.029 0.03043 0.03866 0.21% StateDataPhysBCFunct::() 147 0.0008607 0.01995 0.03725 0.20% StateData::FillBoundary(geom) 143 0 0.01892 0.03618 0.20% MLMG::mgVcycle_up::4 128 0.01425 0.01576 0.03261 0.18% MLMG::computeMLResidual() 18 0.01539 0.02662 0.03158 0.17% MLMG::prepareForSolve() 18 0.02151 0.0232 0.02917 0.16% MLMG::mgVcycle_down::5 128 0.01556 0.01689 0.02876 0.16% MLCellLinOp::defineAuxData() 19 0.02072 0.02302 0.02817 0.15% VisMF::Write(FabArray) 4 0.01584 0.01693 0.02698 0.15% Castro::save_data_for_retry() 3 0.02091 0.02189 0.02486 0.14% MLMG::mgVcycle_up::5 128 0.01429 0.01533 0.02465 0.14% Castro::sum_integrated_quantities() 2 0.0237 0.02385 0.02387 0.13% Castro::problem_post_timestep() 2 0.009679 0.01757 0.02235 0.12% Castro::wd_update() 2 0.009665 0.01754 0.02233 0.12% FabArray::setDomainBndry() 147 0.0004967 0.008482 0.01656 0.09% FabArrayBase::getCPC() 2005 0.0104 0.01207 0.01593 0.09% amrex::Add() 128 0.01488 0.01523 0.01584 0.09% FabArray::norminf() 282 0.01089 0.01173 0.01559 0.09% FabArrayBase::getFB() 6380 0.01139 0.01237 0.01554 0.09% Gravity::test_composite_phi() 1 0.01479 0.01483 0.01485 0.08% Other 2478 0.0966 0.1261 0.1711 0.94% -------------------------------------------------------------------------------------------- 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 9 B 10 B 15 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- Cpu Memory Usage: ------------------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max ------------------------------------------------------------------------------------------------------------------------------- Castro::save_data_for_retry() 416 416 1544 KiB 1545 KiB 1546 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 343 KiB 356 KiB 372 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1010 KiB 1011 KiB 1011 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 961 KiB 962 KiB 962 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 651 KiB 651 KiB 651 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 3525 B 3617 B 4056 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 41 KiB 43 KiB 46 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 13 KiB 18 KiB 20 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 408 KiB 408 KiB 408 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 3179 B 6530 B 8671 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1093 B 1111 B 1138 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 14 KiB 14 KiB 15 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 129 KiB 129 KiB 129 KiB 157 KiB 157 KiB 157 KiB MLMG::prepareForSolve() 6102 6102 9253 B 9804 B 11 KiB 102 KiB 108 KiB 126 KiB Castro::buildMetrics() 160 160 103 KiB 103 KiB 103 KiB 117 KiB 117 KiB 117 KiB Castro::construct_old_source() 640 640 214 B 224 B 268 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 353 B 361 B 406 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 20 B 22 B 37 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 32 KiB 32 KiB 32 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 42 B 42 B 42 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 3 B 8 B 9 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 99 B 101 B 104 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 190 B 392 B 515 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 64 31 KiB 31 KiB 31 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 2432 211 B 218 B 227 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 884 B 919 B 1011 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 117 B 117 B 117 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1676 B 1784 B 1826 B 16 KiB 16 KiB 16 KiB MLMG::compResidual() 32 32 0 B 0 B 1 B 16 KiB 16 KiB 16 KiB MLCellLinOp::defineAuxData() 11704 11704 417 B 466 B 647 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 435 B 436 B 437 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 33 B 49 B 62 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 16 B 25 B 2916 B 2916 B 2916 B OwnerMask() 64 64 0 B 0 B 0 B 2592 B 2592 B 2592 B MLMG::addInterpCorrection() 5632 5632 1 B 14 B 23 B 1024 B 1024 B 1024 B amrex::average_down 5632 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 64 0 B 4 B 6 B 648 B 648 B 648 B MLCGSolver::bicgstab 640 640 0 B 0 B 2 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 114375 0 B 0 B 0 B 8 B 8 B 8 B ------------------------------------------------------------------------------------------------------------------------------- AMReX (24.10-18-g92679babfc2c) finalized