Initializing AMReX (24.10-6-ge1222803739e)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.10-6-ge1222803739e) initialized Starting run at 07:39:01 UTC on 2024-10-08. Successfully read inputs file ... Castro git describe: 24.10-1-g68c1a1537 AMReX git describe: 24.10-6-ge12228037 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.20203989 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001238626 new finest: 0 Gravity::fill_multipole_BCs() time = 0.011496194 Gravity::multilevel_solve_for_new_phi() time = 0.116313542 ... norm of composite residual at level 0 1.193955541e-11 Gravity::fill_multipole_BCs() time = 0.011695666 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.053291471 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.00733996 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 = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477042 at index 18, 179, 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.7694089 at index 17, 181, 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, 63Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347, 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 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 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 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0Invalid density = -9.947567811e-05 at index 0, 193, 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.005580707 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.067656795 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.010071716 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.067084928 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.010376557 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.065409474 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.007365599 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.062343566 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.34969966 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.018115368 Gravity::solve_for_phi() time = 0.097429174 on level 0 Castro::construct_new_gravity() time = 0.099683816 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.060919157 on level 0 Castro::construct_old_gravity() time = 0.004432562 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.060052443 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.349887343 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.018009631 Gravity::solve_for_phi() time = 0.099580071 on level 0 Castro::construct_new_gravity() time = 0.101818229 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.061229515 on level 0 Castro::construct_old_gravity() time = 0.005348745 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.061496741 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.354589737 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.018150489 Gravity::solve_for_phi() time = 0.097556988 on level 0 Castro::construct_new_gravity() time = 0.099702579 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.060633772 on level 0 Castro::construct_old_gravity() time = 0.005030604 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.062957054 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.358789085 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.018467004 Gravity::solve_for_phi() time = 0.098199279 on level 0 Castro::construct_new_gravity() time = 0.100338997 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.060851112 on level 0 Castro::construct_old_gravity() time = 0.004861594 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.060578254 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.362973118 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.01599687 Gravity::solve_for_phi() time = 0.093556793 on level 0 Castro::construct_new_gravity() time = 0.09567362 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.060825136 on level 0 Castro::construct_old_gravity() time = 0.005423779 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.061847787 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.360801226 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.015795065 Gravity::solve_for_phi() time = 0.093583843 on level 0 Castro::construct_new_gravity() time = 0.095759445 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.061047447 on level 0 Castro::construct_old_gravity() time = 0.00373085 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.060879633 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.364316064 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.017107808 Gravity::solve_for_phi() time = 0.097824713 on level 0 Castro::construct_new_gravity() time = 0.099945236 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.060690801 on level 0 Castro::construct_old_gravity() time = 0.002783455 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.05994631 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.358681669 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.018213571 Gravity::solve_for_phi() time = 0.097339314 on level 0 Castro::construct_new_gravity() time = 0.09951601 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.061342901 on level 0 Castro::construct_old_gravity() time = 0.00608908 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.062390651 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.356355596 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.015846768 Gravity::solve_for_phi() time = 0.092824134 on level 0 Castro::construct_new_gravity() time = 0.094965105 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.060241624 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.002765693 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.06006081 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.352229582 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.015925656 Gravity::solve_for_phi() time = 0.093270056 on level 0 Castro::construct_new_gravity() time = 0.095426356 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.060509785 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.002771391 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.060087781 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.357033828 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.016144532 Gravity::solve_for_phi() time = 0.093498516 on level 0 Castro::construct_new_gravity() time = 0.095699526 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.060356845 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.002871536 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.059909416 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.353066259 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.016104867 Gravity::solve_for_phi() time = 0.093470989 on level 0 Castro::construct_new_gravity() time = 0.095597112 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.060976806 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.003013255 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.059826441 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.356998684 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.016199199 Gravity::solve_for_phi() time = 0.093279947 on level 0 Castro::construct_new_gravity() time = 0.095394007 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.060453104 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.002809016 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.059912093 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.352323055 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.016005079 Gravity::solve_for_phi() time = 0.095516106 on level 0 Castro::construct_new_gravity() time = 0.097728132 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.060481217 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.002851869 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.059937404 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.4375695 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.016147674 Gravity::solve_for_phi() time = 0.101402262 on level 0 Castro::construct_new_gravity() time = 0.103622294 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.061034262 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.007013453 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.062904111 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.357799732 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.017008918 Gravity::solve_for_phi() time = 0.094141521 on level 0 Castro::construct_new_gravity() time = 0.096323189 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.060480224 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.002762748 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.059988985 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.358442489 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.00258440801 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.016647218 Gravity::solve_for_phi() time = 0.094053634 on level 0 Castro::construct_new_gravity() time = 0.096188045 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.060391183 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: 12.74154469 [STEP 1] FAB kilobyte spread across MPI nodes: [6973 ... 6973] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.050407 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.006711293 seconds Ending run at 07:39:16 UTC on 2024-10-08. Run time = 15.41052423 Run time without initialization = 12.79963347 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]: 15.4105 ... 15.4105 ... 15.4106 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.7622 5.8786 6.0295 39.13% Amr::InitAmr() 1 1.663 1.663 1.663 10.79% Castro::computeTemp() 80 0.8339 0.9602 1.14 7.40% FabArray::ParallelCopy_finish() 1429 0.3366 0.6353 0.8764 5.69% amrex::Copy() 724 0.7768 0.7918 0.8098 5.25% MultiFab::AddProduct() 3439 0.5495 0.5595 0.5743 3.73% Castro::normalize_species() 79 0.5474 0.558 0.5699 3.70% Castro::reset_internal_energy(Fab) 160 0.516 0.527 0.5424 3.52% check_for_negative_density() 20 0.1957 0.3732 0.5002 3.25% MLPoisson::Fsmooth() 6144 0.3491 0.381 0.4502 2.92% Castro::initMFs() 1 0.3518 0.3643 0.3794 2.46% FabArray::setVal() 5013 0.3543 0.3639 0.3735 2.42% MultiFab::Subtract() 159 0.3394 0.346 0.3541 2.30% Gravity::fill_multipole_BCs() 19 0.2145 0.2754 0.3061 1.99% FabArray::ParallelCopy_nowait() 1429 0.2064 0.2199 0.2439 1.58% FillBoundary_finish() 6380 0.1248 0.1495 0.2037 1.32% FabArray::sum() 3439 0.1832 0.1866 0.1918 1.24% Castro::initData() 1 0.003344 0.02929 0.1752 1.14% TagBoxArray::collate() 2 0.0001971 0.03355 0.1748 1.13% FillBoundary_nowait() 6380 0.09361 0.1223 0.1598 1.04% MLCellLinOp::applyBC() 7148 0.06434 0.07997 0.1237 0.80% MLPoisson::Fapply() 986 0.09618 0.1017 0.1152 0.75% FabArray::Saxpy() 264 0.06889 0.07038 0.07642 0.50% Castro::construct_new_gravity_source() 17 0.06823 0.07023 0.07285 0.47% FabArray::Xpay() 945 0.05219 0.05493 0.0612 0.40% Castro::construct_old_gravity_source() 20 0.04531 0.04634 0.0485 0.31% MultiFab::contains_nan() 37 0.04587 0.04685 0.04766 0.31% MLMG::addInterpCorrection() 768 0.03408 0.0371 0.04432 0.29% AmrLevel::derive() 18 0.04043 0.04155 0.04375 0.28% amrex::average_down 768 0.02517 0.02768 0.03435 0.22% StateData::FillBoundary(geom) 143 0 0.01789 0.03424 0.22% VisMF::Write(FabArray) 4 0.01591 0.01845 0.03406 0.22% Castro::reset_internal_energy(MultiFab) 80 0.02329 0.0241 0.02508 0.16% Amr::writePlotFile() 2 0.007966 0.01999 0.02414 0.16% MLCellLinOp::defineAuxData() 19 0.01616 0.0171 0.0202 0.13% Castro::enforce_min_density() 79 0.01879 0.01928 0.01989 0.13% MLCGSolver::bicgstab 8 0 0.001098 0.01756 0.11% FabArray::FillBoundary() 6380 0.01363 0.01422 0.01672 0.11% Castro::construct_new_source() 102 0.01592 0.01616 0.01663 0.11% FabArray::norminf() 282 0.01131 0.0118 0.01584 0.10% amrex::Add() 128 0.01413 0.01437 0.01456 0.09% Castro::evaluate_source_change() 196 0.01371 0.01408 0.01449 0.09% Castro::wd_update() 2 0.004799 0.01015 0.01429 0.09% FabArray::setDomainBndry() 147 0.0004773 0.007575 0.01421 0.09% FillPatchIterator::Initialize 147 0.01289 0.01348 0.01391 0.09% MLCellLinOp::smooth() 3072 0.01187 0.0123 0.01378 0.09% FabArrayBase::getFB() 6380 0.01049 0.01101 0.01331 0.09% MLMG::prepareForSolve() 18 0.01029 0.01082 0.01155 0.07% Castro::construct_old_source() 120 0.01015 0.01039 0.01108 0.07% MLMG::ResNormInf() 146 0.005548 0.008356 0.0106 0.07% Gravity::actual_solve_with_mlmg() 19 0.007544 0.007957 0.009134 0.06% Amr::writeSmallPlotFile() 2 0.006678 0.008376 0.008855 0.06% FabArrayBase::getCPC() 2005 0.007566 0.007793 0.008103 0.05% Gravity::get_old_grav_vector() 20 0.006507 0.006709 0.007415 0.05% Gravity::get_new_grav_vector() 18 0.006044 0.006192 0.00636 0.04% Castro::initialize_do_advance() 20 0.004384 0.004917 0.005573 0.04% MLCellLinOp::apply() 986 0.001996 0.00222 0.004562 0.03% MLMG::solve() 18 0.0008605 0.00274 0.004299 0.03% Castro::sum_integrated_quantities() 2 0.00235 0.00377 0.004269 0.03% FabArray::ParallelCopy() 1429 0.003457 0.003608 0.003999 0.03% AmrLevel::FillPatch() 147 0.00242 0.002547 0.002738 0.02% MLCellLinOp::correctionResidual() 768 0.001895 0.001953 0.002038 0.01% MLMG::mgVcycle() 128 0.001606 0.001669 0.001785 0.01% MLMG:computeResOfCorrection() 768 0.001473 0.001521 0.001608 0.01% Castro::save_data_for_retry() 3 0.001393 0.001453 0.001545 0.01% MLMG::actualBottomSolve() 128 0.0001327 0.0001978 0.001094 0.01% Castro::subcycle_advance_ctu() 1 0.0002553 0.0003447 0.001022 0.01% StateDataPhysBCFunct::() 147 0.0007586 0.0008966 0.0009937 0.01% MLMG::mgVcycle_down::2 128 0.0006375 0.0007456 0.000954 0.01% main() 1 0.0007592 0.000815 0.0008681 0.01% MLMG::mgVcycle_down::5 128 0.0005881 0.0006467 0.0008534 0.01% FillPatchIterator::FillFromLevel0() 147 0.000696 0.0007239 0.0007702 0.00% Gravity::solve_for_phi() 17 0.0006692 0.0007022 0.0007626 0.00% MLMG::mgVcycle_down::3 128 0.0005883 0.0006375 0.0007523 0.00% MLMG::mgVcycle_down::1 128 0.0006757 0.000691 0.0007225 0.00% MLMG::mgVcycle_down::0 128 0.0006765 0.0006946 0.0007184 0.00% MLMG::mgVcycle_down::4 128 0.0005707 0.0005956 0.0006911 0.00% AmrMesh::MakeNewGrids() 2 0.0005616 0.0005861 0.0006027 0.00% MLMG::mgVcycle_up::1 128 0.0004797 0.0004974 0.0005463 0.00% MLMG::mgVcycle_up::0 128 0.0004893 0.0005055 0.0005352 0.00% FillPatchSingleLevel 147 0.0004623 0.0004881 0.0005293 0.00% MLMG::mgVcycle_up::2 128 0.0004137 0.0004548 0.0005257 0.00% Castro::clean_state() 79 0.0004477 0.0004772 0.000517 0.00% MLMG::mgVcycle_up::3 128 0.0004074 0.0004273 0.0004846 0.00% MLCellLinOp::solutionResidual() 147 0.0004401 0.0004555 0.0004749 0.00% MLMG::mgVcycle_up::4 128 0.0004092 0.0004283 0.0004718 0.00% Castro::do_advance_ctu() 20 0.000392 0.0004276 0.0004649 0.00% Castro::Castro() 1 0.0003731 0.0004125 0.0004539 0.00% MLMG::mgVcycle_up::5 128 0.0004093 0.0004261 0.0004523 0.00% MLMG::oneIter() 128 0.0003903 0.0004042 0.0004339 0.00% Castro::do_old_sources() 20 0.0002786 0.0002918 0.0003128 0.00% Castro::construct_new_gravity() 17 0.000158 0.0001826 0.000302 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002282 0.0002415 0.00028 0.00% Castro::problem_post_timestep() 2 1.128e-05 2.823e-05 0.0002534 0.00% MLMG::computeResidual() 128 0.0002385 0.0002453 0.0002517 0.00% MLMG::mgVcycle_bottom 128 0.0002245 0.0002328 0.0002471 0.00% Castro::do_new_sources() 17 0.0002087 0.0002206 0.0002341 0.00% Gravity::update_max_rhs() 2 0.0002005 0.0002142 0.0002253 0.00% Amr::coarseTimeStep() 1 0.0001371 0.0001601 0.0001647 0.00% Castro::initialize_advance() 1 0.0001285 0.0001498 0.0001634 0.00% Castro::retry_advance_ctu() 11 6.413e-05 8.192e-05 0.0001616 0.00% TagBoxArray::mapPRD 2 0.000125 0.0001325 0.0001484 0.00% Amr::defBaseLevel() 1 0.0001063 0.0001243 0.0001411 0.00% MLPoisson::define() 19 0.0001223 0.0001296 0.0001363 0.00% Castro::construct_old_gravity() 20 0.0001124 0.0001189 0.0001328 0.00% Castro::check_for_nan() 37 0.0001072 0.0001169 0.0001276 0.00% Gravity::test_composite_phi() 1 0.0001075 0.000112 0.0001232 0.00% Castro::apply_source_to_state() 37 9.759e-05 0.0001053 0.0001127 0.00% Castro::create_source_corrector() 17 5.647e-05 6.491e-05 7.677e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.28e-05 6.6e-05 6.84e-05 0.00% Amr::grid_places() 2 5.364e-05 5.596e-05 6.045e-05 0.00% Castro::expand_state() 20 5.01e-05 5.318e-05 5.831e-05 0.00% MLMG::computeMLResidual() 18 3.95e-05 4.539e-05 5.397e-05 0.00% Castro::derive() 18 4.071e-05 4.348e-05 5e-05 0.00% Amr::FinalizeInit() 1 9.12e-06 1.257e-05 4.978e-05 0.00% Amr::regrid() 1 1.621e-05 2.063e-05 4.694e-05 0.00% Gravity::actual_multilevel_solve() 1 3.311e-05 3.517e-05 3.726e-05 0.00% Castro::advance() 1 6.964e-06 9.528e-06 2.789e-05 0.00% Amr::timeStep() 1 1.488e-05 1.712e-05 2.029e-05 0.00% Castro::post_init() 1 1.266e-05 1.396e-05 1.725e-05 0.00% Amr::InitializeInit() 1 7.96e-06 9.472e-06 1.405e-05 0.00% Castro::post_regrid() 2 1.113e-05 1.185e-05 1.281e-05 0.00% Amr::bldFineLevels() 1 7.419e-06 8.105e-06 9.573e-06 0.00% Amr::init() 1 6.933e-06 7.848e-06 8.812e-06 0.00% Castro::post_timestep() 1 5.649e-06 6.403e-06 6.999e-06 0.00% Amr::initialInit() 1 2.967e-06 3.416e-06 4.081e-06 0.00% Other 2474 0.0463 0.06248 0.0913 0.59% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 15.41 15.41 15.41 100.00% Amr::coarseTimeStep() 1 12.8 12.8 12.8 83.05% Amr::timeStep() 1 12.74 12.74 12.74 82.68% Castro::advance() 1 12.67 12.68 12.68 82.29% Castro::subcycle_advance_ctu() 1 12.63 12.63 12.63 81.96% Castro::do_advance_ctu() 20 12.58 12.59 12.59 81.71% Castro::construct_ctu_hydro_source() 20 6.992 7.034 7.08 45.95% Castro::clean_state() 79 3.643 3.806 4.028 26.14% Castro::computeTemp() 80 2.326 2.473 2.678 17.38% Gravity::solve_phi_with_mlmg() 19 1.658 1.723 1.754 11.38% Amr::InitAmr() 1 1.663 1.663 1.663 10.79% Castro::construct_new_gravity() 17 1.568 1.632 1.663 10.79% Gravity::solve_for_phi() 17 1.53 1.594 1.626 10.55% Castro::reset_internal_energy(MultiFab) 80 1.486 1.513 1.549 10.05% Castro::initialize_do_advance() 20 1.421 1.46 1.502 9.74% Gravity::actual_solve_with_mlmg() 19 1.44 1.444 1.449 9.40% MLMG::solve() 18 1.384 1.389 1.391 9.02% MLMG::oneIter() 128 1.243 1.246 1.248 8.10% MLMG::mgVcycle() 128 1.228 1.231 1.233 8.00% Castro::do_old_sources() 20 1.143 1.195 1.217 7.89% FabArray::ParallelCopy() 1429 0.6059 0.8792 1.111 7.21% Castro::do_new_sources() 17 0.9172 0.972 1.027 6.67% Castro::evaluate_source_change() 196 0.9481 0.9655 0.991 6.43% Amr::init() 1 0.9473 0.9473 0.9473 6.15% Amr::initialInit() 1 0.8862 0.8862 0.8863 5.75% FabArray::ParallelCopy_finish() 1429 0.3366 0.6353 0.8764 5.69% MLCellLinOp::smooth() 3072 0.6413 0.6982 0.8453 5.49% Castro::enforce_min_density() 79 0.7877 0.8015 0.8209 5.33% amrex::Copy() 724 0.7768 0.7918 0.8098 5.25% AmrLevel::FillPatch() 147 0.5388 0.6648 0.7723 5.01% Amr::InitializeInit() 1 0.4363 0.4767 0.6366 4.13% Amr::defBaseLevel() 1 0.4363 0.4767 0.6366 4.13% FillPatchIterator::Initialize 147 0.3558 0.4837 0.5942 3.86% MultiFab::AddProduct() 3439 0.5495 0.5595 0.5743 3.73% Castro::normalize_species() 79 0.5474 0.558 0.5699 3.70% FillPatchIterator::FillFromLevel0() 147 0.3277 0.4626 0.5673 3.68% FillPatchSingleLevel 147 0.327 0.4619 0.5665 3.68% Castro::reset_internal_energy(Fab) 160 0.516 0.527 0.5424 3.52% check_for_negative_density() 20 0.1957 0.3732 0.5002 3.25% MLMG::mgVcycle_up::1 128 0.1132 0.2961 0.4798 3.11% MLCellLinOp::applyBC() 7148 0.33 0.3777 0.4764 3.09% MLPoisson::Fsmooth() 6144 0.3491 0.381 0.4502 2.92% Amr::FinalizeInit() 1 0.2496 0.4095 0.4499 2.92% MLMG::addInterpCorrection() 768 0.06754 0.332 0.4109 2.67% Castro::Castro() 1 0.3572 0.3698 0.385 2.50% Castro::initMFs() 1 0.352 0.3644 0.3796 2.46% FabArray::setVal() 5013 0.3543 0.3639 0.3735 2.42% MultiFab::Subtract() 159 0.3394 0.346 0.3541 2.30% FabArray::FillBoundary() 6380 0.2537 0.2977 0.3527 2.29% MLMG::mgVcycle_up::2 128 0.01367 0.1302 0.3087 2.00% Gravity::fill_multipole_BCs() 19 0.2167 0.2774 0.3082 2.00% MLMG::mgVcycle_down::0 128 0.2819 0.2846 0.2872 1.86% Castro::expand_state() 20 0.2003 0.2389 0.2741 1.78% FabArray::ParallelCopy_nowait() 1429 0.2263 0.2403 0.2654 1.72% Castro::initData() 1 0.07738 0.1051 0.2522 1.64% Amr::bldFineLevels() 1 0.04709 0.2063 0.2458 1.60% MLMG::mgVcycle_up::0 128 0.2171 0.2182 0.2189 1.42% MLCellLinOp::apply() 986 0.1637 0.1742 0.212 1.38% FillBoundary_finish() 6380 0.1248 0.1495 0.2037 1.32% Amr::grid_places() 2 0.002956 0.1628 0.2033 1.32% AmrMesh::MakeNewGrids() 2 0.002902 0.1628 0.2032 1.32% TagBoxArray::mapPRD 2 0.001197 0.1278 0.2008 1.30% MLMG::mgVcycle_up::3 128 0.01358 0.02774 0.1936 1.26% FabArray::sum() 3439 0.1832 0.1866 0.1918 1.24% TagBoxArray::collate() 2 0.0001971 0.03355 0.1748 1.13% FillBoundary_nowait() 6380 0.1046 0.134 0.1724 1.12% Castro::post_init() 1 0.1515 0.1594 0.1617 1.05% MLMG:computeResOfCorrection() 768 0.1224 0.1324 0.1576 1.02% MLCellLinOp::correctionResidual() 768 0.121 0.1309 0.1561 1.01% MLMG::mgVcycle_down::1 128 0.129 0.1336 0.1387 0.90% Castro::construct_new_source() 102 0.1235 0.1262 0.1295 0.84% Castro::create_source_corrector() 17 0.05641 0.08442 0.1203 0.78% Gravity::multilevel_solve_for_new_phi() 1 0.1162 0.1163 0.1163 0.75% Gravity::actual_multilevel_solve() 1 0.1162 0.1162 0.1163 0.75% MLPoisson::Fapply() 986 0.09618 0.1017 0.1152 0.75% Amr::writePlotFile() 2 0.1041 0.1041 0.1042 0.68% MLCellLinOp::solutionResidual() 147 0.09162 0.09894 0.1032 0.67% Castro::post_regrid() 2 0.08494 0.08721 0.09507 0.62% MLMG::mgVcycle_down::2 128 0.01596 0.04824 0.08468 0.55% Castro::construct_old_source() 120 0.07865 0.08026 0.08368 0.54% MLMG::computeResidual() 128 0.07932 0.08078 0.08233 0.53% Castro::construct_old_gravity() 20 0.05641 0.0689 0.07914 0.51% Gravity::get_old_grav_vector() 20 0.05629 0.06878 0.07903 0.51% FabArray::Saxpy() 264 0.06889 0.07038 0.07642 0.50% MLMG::mgVcycle_bottom 128 0.0003588 0.005101 0.07607 0.49% MLMG::actualBottomSolve() 128 0.0001327 0.004868 0.07582 0.49% Castro::derive() 18 0.06888 0.07099 0.07515 0.49% AmrLevel::derive() 18 0.06884 0.07094 0.0751 0.49% MLCGSolver::bicgstab 8 0 0.00464 0.07425 0.48% Castro::construct_new_gravity_source() 17 0.06823 0.07023 0.07285 0.47% MLMG::mgVcycle_down::3 128 0.01526 0.02239 0.06999 0.45% Castro::post_timestep() 1 0.05985 0.06358 0.06677 0.43% FabArray::Xpay() 945 0.05219 0.05493 0.0612 0.40% amrex::average_down 768 0.03405 0.04188 0.0592 0.38% Castro::initialize_advance() 1 0.04848 0.05004 0.05183 0.34% Castro::construct_old_gravity_source() 20 0.04531 0.04634 0.0485 0.31% Castro::check_for_nan() 37 0.04598 0.04697 0.04778 0.31% MultiFab::contains_nan() 37 0.04587 0.04685 0.04766 0.31% Amr::regrid() 1 0.04375 0.0447 0.04658 0.30% Castro::retry_advance_ctu() 11 0.03608 0.03818 0.04437 0.29% MLPoisson::define() 19 0.03121 0.03646 0.04281 0.28% MLMG::mgVcycle_down::4 128 0.01512 0.01707 0.0396 0.26% Gravity::get_new_grav_vector() 18 0.03627 0.03688 0.03808 0.25% StateDataPhysBCFunct::() 147 0.0007586 0.01879 0.03519 0.23% StateData::FillBoundary(geom) 143 0 0.01789 0.03424 0.22% VisMF::Write(FabArray) 4 0.01591 0.01845 0.03406 0.22% MLMG::mgVcycle_up::4 128 0.0137 0.01517 0.03132 0.20% Castro::apply_source_to_state() 37 0.02804 0.02847 0.02891 0.19% MLMG::mgVcycle_down::5 128 0.01508 0.01637 0.02775 0.18% MLCellLinOp::defineAuxData() 19 0.0204 0.02212 0.02719 0.18% MLMG::prepareForSolve() 18 0.0207 0.022 0.02424 0.16% MLMG::mgVcycle_up::5 128 0.01379 0.01478 0.02373 0.15% Castro::sum_integrated_quantities() 2 0.02326 0.02343 0.02344 0.15% Castro::save_data_for_retry() 3 0.01882 0.01953 0.02285 0.15% MLMG::computeMLResidual() 18 0.01175 0.01763 0.02203 0.14% MLMG::ResNormInf() 146 0.01573 0.01847 0.02055 0.13% Castro::problem_post_timestep() 2 0.00932 0.01379 0.01763 0.11% Castro::wd_update() 2 0.009306 0.01377 0.01762 0.11% FabArray::norminf() 282 0.01131 0.0118 0.01584 0.10% Gravity::test_composite_phi() 1 0.01535 0.0154 0.01543 0.10% FabArrayBase::getFB() 6380 0.01083 0.01176 0.01495 0.10% FabArrayBase::getCPC() 2005 0.01001 0.01143 0.01494 0.10% amrex::Add() 128 0.01413 0.01437 0.01456 0.09% Amr::writeSmallPlotFile() 2 0.01428 0.01429 0.01429 0.09% FabArray::setDomainBndry() 147 0.0004773 0.007575 0.01421 0.09% Gravity::update_max_rhs() 2 0.001591 0.009461 0.01171 0.08% Other 2474 0.07679 0.1029 0.144 0.93% -------------------------------------------------------------------------------------------- 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 11 B 14 B 21 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 1540 KiB 1541 KiB 1541 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 392 KiB 400 KiB 410 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1023 KiB 1024 KiB 1024 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 974 KiB 974 KiB 975 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 653 KiB 653 KiB 654 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 3973 B 4185 B 4810 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 47 KiB 47 KiB 49 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 11 KiB 13 KiB 15 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 409 KiB 409 KiB 409 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 2051 B 4072 B 5413 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1241 B 1277 B 1331 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 16 KiB 16 KiB 16 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 9430 B 9982 B 11 KiB 102 KiB 108 KiB 126 KiB Castro::buildMetrics() 160 160 104 KiB 104 KiB 104 KiB 117 KiB 117 KiB 117 KiB Castro::construct_old_source() 640 640 240 B 244 B 254 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 404 B 410 B 419 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 29 B 32 B 47 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 33 KiB 33 KiB 33 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 52 B 52 B 52 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 3 B 11 B 14 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 106 B 107 B 111 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 167 B 204 B 235 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 224 B 229 B 232 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1002 B 1020 B 1048 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 123 B 123 B 123 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1625 B 1693 B 1727 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 426 B 476 B 660 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 443 B 444 B 445 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 24 B 34 B 44 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 21 B 33 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 15 B 25 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 6 B 8 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-6-ge1222803739e) finalized