Initializing AMReX (24.10-16-gc0602940a47c)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.10-16-gc0602940a47c) initialized Starting run at 08:49:56 UTC on 2024-10-23. Successfully read inputs file ... Castro git describe: 24.10-6-gb6503500c AMReX git describe: 24.10-16-gc0602940a 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.407295028 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001944259 new finest: 0 Gravity::fill_multipole_BCs() time = 0.02103065 Gravity::multilevel_solve_for_new_phi() time = 0.190524655 ... norm of composite residual at level 0 1.193955541e-11 Gravity::fill_multipole_BCs() time = 0.021301891 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.405535019 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.006543397 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 = -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 1Invalid density = -618.7284155 at index , 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 1910, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = , 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.99403820 Invalid density = -1.089487306 at index 6, 191, 0 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071149 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477042 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035501 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.767071355 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.767071355 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Castro::construct_old_gravity() time = 0.00826777 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.113010805 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.072667312 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.11627327 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, 177, 0 Invalid density = -0.1060744788 at index 20, 78, 0 Castro::construct_old_gravity() time = 0.061818873 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.112354297 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.0614826 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.109296885 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.664837371 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.06108751 Gravity::solve_for_phi() time = 0.219338718 on level 0 Castro::construct_new_gravity() time = 0.222981637 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.112672788 on level 0 Castro::construct_old_gravity() time = 0.012175265 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.109867057 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.67173532 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.062056205 Gravity::solve_for_phi() time = 0.212132651 on level 0 Castro::construct_new_gravity() time = 0.215989471 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.111825299 on level 0 Castro::construct_old_gravity() time = 0.006839418 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.1096245 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.66846914 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.061351166 Gravity::solve_for_phi() time = 0.200931354 on level 0 Castro::construct_new_gravity() time = 0.204433013 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.113389733 on level 0 Castro::construct_old_gravity() time = 0.007946154 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.11009863 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.674724994 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.063038897 Gravity::solve_for_phi() time = 0.214641966 on level 0 Castro::construct_new_gravity() time = 0.220460258 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.117604006 on level 0 Castro::construct_old_gravity() time = 0.020103576 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.116377477 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.738315837 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.071913124 Gravity::solve_for_phi() time = 0.223958804 on level 0 Castro::construct_new_gravity() time = 0.227765292 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.117913049 on level 0 Castro::construct_old_gravity() time = 0.020683634 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.116519808 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.720851537 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.07173094 Gravity::solve_for_phi() time = 0.222879749 on level 0 Castro::construct_new_gravity() time = 0.226639759 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.117667462 on level 0 Castro::construct_old_gravity() time = 0.016028151 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.109559588 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.687263045 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.071268136 Gravity::solve_for_phi() time = 0.221568905 on level 0 Castro::construct_new_gravity() time = 0.225348778 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.116817723 on level 0 Castro::construct_old_gravity() time = 0.012866615 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.110342323 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.691256329 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.071539595 Gravity::solve_for_phi() time = 0.222948818 on level 0 Castro::construct_new_gravity() time = 0.226509127 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.111692231 on level 0 Castro::construct_old_gravity() time = 0.00682477 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.108735213 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.686096016 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.071854876 Gravity::solve_for_phi() time = 0.219645607 on level 0 Castro::construct_new_gravity() time = 0.223382097 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.118046629 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.021095821 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.116203301 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.731240285 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.065149524 Gravity::solve_for_phi() time = 0.220871222 on level 0 Castro::construct_new_gravity() time = 0.224882236 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.121876095 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.018760484 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.116909134 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.732873552 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.0719048 Gravity::solve_for_phi() time = 0.222740503 on level 0 Castro::construct_new_gravity() time = 0.226555508 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.118529415 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.048556631 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.314423645 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.683630225 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.06584952 Gravity::solve_for_phi() time = 0.210308452 on level 0 Castro::construct_new_gravity() time = 0.214071458 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.111576814 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.006797177 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.108519181 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.729773862 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.07314563 Gravity::solve_for_phi() time = 0.236434418 on level 0 Castro::construct_new_gravity() time = 0.240004513 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.118766226 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.017920985 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.117276825 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.721954924 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.073859104 Gravity::solve_for_phi() time = 0.22220777 on level 0 Castro::construct_new_gravity() time = 0.226067818 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.118950459 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.010537176 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.110291339 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.689836062 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.074329645 Gravity::solve_for_phi() time = 0.22452443 on level 0 Castro::construct_new_gravity() time = 0.228365988 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.118424798 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.018443048 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.117296926 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.690822981 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.074690391 Gravity::solve_for_phi() time = 0.230353964 on level 0 Castro::construct_new_gravity() time = 0.236163993 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.119254099 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.006722992 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.109003754 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.692763078 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.001338865299 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.074697244 Gravity::solve_for_phi() time = 0.222129312 on level 0 Castro::construct_new_gravity() time = 0.225852022 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.118612085 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: 24.60300992 [STEP 1] FAB kilobyte spread across MPI nodes: [6973 ... 6973] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.095725778 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.006165216 seconds Ending run at 08:50:25 UTC on 2024-10-23. Run time = 28.81709255 Run time without initialization = 24.70630425 Average number of zones advanced per microsecond: 0.001 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 28.8167 ... 28.8169 ... 28.8172 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 6.0800 8.3988 11.5793 40.18% check_for_negative_density() 20 0.4393 3.673 6.041 20.96% FabArray::ParallelCopy_finish() 1429 1.037 3.096 4.857 16.85% Amr::InitAmr() 1 2.134 2.134 2.134 7.41% Castro::computeTemp() 80 0.8921 1.372 1.941 6.74% amrex::Copy() 724 0.8153 1.105 1.458 5.06% Castro::normalize_species() 79 0.5789 0.8735 1.256 4.36% Gravity::fill_multipole_BCs() 19 0.4873 0.8829 1.209 4.20% Castro::reset_internal_energy(Fab) 160 0.5456 0.7859 1.086 3.77% FillBoundary_finish() 6380 0.3643 0.7022 1.026 3.56% MultiFab::AddProduct() 3439 0.5785 0.7655 0.9914 3.44% MLPoisson::Fsmooth() 6144 0.3624 0.554 0.8551 2.97% Castro::initMFs() 1 0.3651 0.5375 0.6977 2.42% MultiFab::Subtract() 159 0.3618 0.4571 0.5741 1.99% FabArray::setVal() 5013 0.3733 0.4659 0.5713 1.98% FabArray::ParallelCopy_nowait() 1429 0.2105 0.3173 0.4338 1.51% Amr::writePlotFile() 2 0.327 0.3844 0.418 1.45% FabArray::sum() 3439 0.1986 0.2908 0.4082 1.42% Castro::sum_integrated_quantities() 2 0.005158 0.2804 0.3059 1.06% Castro::problem_post_timestep() 2 1.4e-05 0.01779 0.2844 0.99% FillBoundary_nowait() 6380 0.09226 0.163 0.2324 0.81% MLPoisson::Fapply() 986 0.1001 0.1555 0.2322 0.81% MLCellLinOp::applyBC() 7148 0.07227 0.1094 0.2144 0.74% Castro::initData() 1 0.003544 0.03108 0.1718 0.60% Castro::construct_new_gravity_source() 17 0.07205 0.1113 0.1568 0.54% FabArray::Saxpy() 264 0.0716 0.0953 0.1319 0.46% Castro::wd_update() 2 0.008323 0.0809 0.1249 0.43% MLMG::ResNormInf() 146 0.03369 0.07775 0.1146 0.40% FabArray::Xpay() 945 0.05344 0.07728 0.1118 0.39% Castro::construct_old_gravity_source() 20 0.04779 0.07305 0.1032 0.36% AmrLevel::derive() 18 0.04242 0.05951 0.08537 0.30% MultiFab::contains_nan() 37 0.04768 0.06332 0.08176 0.28% MultiFab::max() 2 0.0002903 0.04006 0.08027 0.28% MLMG::addInterpCorrection() 768 0.03514 0.05134 0.07672 0.27% StateData::FillBoundary(geom) 143 0 0.02448 0.06362 0.22% amrex::average_down 768 0.02609 0.03752 0.05737 0.20% VisMF::Write(FabArray) 4 0.01411 0.02065 0.053 0.18% Castro::construct_new_source() 102 0.01647 0.02514 0.03543 0.12% Castro::reset_internal_energy(MultiFab) 80 0.02572 0.02937 0.03335 0.12% MLCellLinOp::defineAuxData() 19 0.01635 0.02255 0.03297 0.11% FabArray::norminf() 282 0.01231 0.0173 0.0294 0.10% Castro::enforce_min_density() 79 0.02 0.02425 0.02912 0.10% MLCGSolver::bicgstab 8 0 0.001726 0.02761 0.10% FabArray::FillBoundary() 6380 0.01417 0.0186 0.02709 0.09% MLMG::solve() 18 0.002829 0.01269 0.02641 0.09% amrex::Add() 128 0.01461 0.02003 0.02605 0.09% Castro::evaluate_source_change() 196 0.0152 0.02005 0.02536 0.09% FabArray::setDomainBndry() 147 0.0006515 0.009734 0.023 0.08% Castro::construct_old_source() 120 0.01071 0.01598 0.02215 0.08% MLCellLinOp::smooth() 3072 0.0118 0.01623 0.02161 0.07% FillPatchIterator::Initialize 147 0.01401 0.01595 0.01811 0.06% MLMG::prepareForSolve() 18 0.01065 0.01366 0.01806 0.06% Gravity::actual_solve_with_mlmg() 19 0.007984 0.01007 0.0138 0.05% Gravity::get_old_grav_vector() 20 0.006917 0.009885 0.01324 0.05% Gravity::get_new_grav_vector() 18 0.00633 0.009162 0.012 0.04% Castro::save_data_for_retry() 3 0.001466 0.002096 0.007619 0.03% MLCellLinOp::apply() 986 0.002086 0.002967 0.007414 0.03% Castro::initialize_do_advance() 20 0.004741 0.005412 0.006314 0.02% FabArray::ParallelCopy() 1429 0.003718 0.004856 0.006184 0.02% AmrLevel::FillPatch() 147 0.002798 0.003399 0.004135 0.01% MLCellLinOp::correctionResidual() 768 0.001968 0.002533 0.003262 0.01% MLMG::mgVcycle() 128 0.001699 0.00214 0.002602 0.01% MLMG:computeResOfCorrection() 768 0.001487 0.001936 0.002579 0.01% StateDataPhysBCFunct::() 147 0.001112 0.001281 0.001671 0.01% MLMG::actualBottomSolve() 128 0.0001423 0.0002651 0.001669 0.01% main() 1 0.0007259 0.0009486 0.001324 0.00% Castro::subcycle_advance_ctu() 1 0.0003192 0.0004396 0.001307 0.00% MLMG::mgVcycle_down::2 128 0.0006542 0.000909 0.001306 0.00% MLMG::mgVcycle_down::1 128 0.0007182 0.0009388 0.001221 0.00% MLMG::mgVcycle_down::3 128 0.0006138 0.0008382 0.001203 0.00% MLMG::mgVcycle_down::0 128 0.0007053 0.0009278 0.001171 0.00% FillPatchIterator::FillFromLevel0() 147 0.0007825 0.0009333 0.001054 0.00% Gravity::solve_for_phi() 17 0.0007546 0.0008954 0.001053 0.00% MLMG::mgVcycle_down::4 128 0.0005902 0.000779 0.001035 0.00% MLMG::mgVcycle_down::5 128 0.0006202 0.0007954 0.001024 0.00% AmrMesh::MakeNewGrids() 2 0.0005846 0.0007938 0.001006 0.00% MLMG::mgVcycle_up::0 128 0.0005062 0.0006822 0.0008915 0.00% Castro::clean_state() 79 0.0005555 0.0007146 0.0008892 0.00% MLMG::mgVcycle_up::1 128 0.0005044 0.0006605 0.0008571 0.00% FillPatchSingleLevel 147 0.0005443 0.0006733 0.0007796 0.00% Castro::do_advance_ctu() 20 0.0004747 0.0006003 0.0007713 0.00% MLCellLinOp::solutionResidual() 147 0.0004785 0.0006004 0.0007485 0.00% MLMG::mgVcycle_up::2 128 0.000414 0.0005781 0.0007466 0.00% MLMG::mgVcycle_up::3 128 0.0004103 0.0005394 0.0007169 0.00% MLMG::mgVcycle_up::4 128 0.0004162 0.0005329 0.0006866 0.00% MLMG::mgVcycle_up::5 128 0.0004183 0.0005228 0.0006842 0.00% MLMG::oneIter() 128 0.0004299 0.0005468 0.0006694 0.00% Castro::Castro() 1 0.0003972 0.000487 0.0005627 0.00% Castro::do_old_sources() 20 0.000338 0.00042 0.0005468 0.00% MLMG::computeResidual() 128 0.0002581 0.0003308 0.0004245 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002687 0.0003273 0.0004113 0.00% Gravity::update_max_rhs() 2 0.0002374 0.0003028 0.0004109 0.00% Castro::construct_new_gravity() 17 0.0001808 0.0002288 0.0004002 0.00% Castro::do_new_sources() 17 0.0002393 0.0003128 0.0003921 0.00% MLMG::mgVcycle_bottom 128 0.0002292 0.0002888 0.0003593 0.00% Castro::check_for_nan() 37 0.0001276 0.0001806 0.0002453 0.00% Castro::retry_advance_ctu() 11 7.702e-05 0.0001127 0.0002402 0.00% Castro::initialize_advance() 1 0.0001172 0.00015 0.0002218 0.00% Castro::construct_old_gravity() 20 0.0001285 0.0001543 0.0002184 0.00% TagBoxArray::mapPRD 2 0.0001298 0.0001714 0.0002182 0.00% Amr::coarseTimeStep() 1 0.0001564 0.0001991 0.0002068 0.00% Castro::apply_source_to_state() 37 0.0001109 0.0001528 0.0001985 0.00% Castro::create_source_corrector() 17 7.215e-05 0.0001022 0.0001913 0.00% MLPoisson::define() 19 0.0001265 0.0001574 0.0001878 0.00% Gravity::test_composite_phi() 1 0.0001178 0.0001486 0.0001841 0.00% Gravity::multilevel_solve_for_new_phi() 1 7.377e-05 9.31e-05 0.0001423 0.00% Amr::defBaseLevel() 1 9.581e-05 0.0001216 0.0001418 0.00% Amr::grid_places() 2 6.01e-05 7.465e-05 9.124e-05 0.00% Castro::expand_state() 20 5.813e-05 7.332e-05 8.891e-05 0.00% Amr::FinalizeInit() 1 1.126e-05 2.093e-05 8.866e-05 0.00% Castro::derive() 18 4.448e-05 6.222e-05 8.111e-05 0.00% MLMG::computeMLResidual() 18 4.225e-05 5.371e-05 6.479e-05 0.00% Gravity::actual_multilevel_solve() 1 3.633e-05 4.734e-05 6.325e-05 0.00% Amr::regrid() 1 1.802e-05 2.241e-05 5.417e-05 0.00% Castro::post_init() 1 1.734e-05 2.648e-05 3.828e-05 0.00% Amr::timeStep() 1 1.552e-05 2.129e-05 3.368e-05 0.00% Castro::advance() 1 8.098e-06 1.196e-05 2.885e-05 0.00% Amr::init() 1 7.909e-06 1.227e-05 1.984e-05 0.00% Castro::post_regrid() 2 1.159e-05 1.542e-05 1.813e-05 0.00% Amr::InitializeInit() 1 8.312e-06 1.063e-05 1.629e-05 0.00% Amr::bldFineLevels() 1 7.964e-06 9.519e-06 1.165e-05 0.00% Castro::post_timestep() 1 6.338e-06 8.207e-06 1.054e-05 0.00% Amr::initialInit() 1 3.902e-06 5.897e-06 9.284e-06 0.00% Other 10861 0.07299 0.1114 0.1764 0.61% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 28.82 28.82 28.82 100.00% Amr::coarseTimeStep() 1 24.71 24.71 24.71 85.73% Amr::timeStep() 1 24.6 24.6 24.6 85.38% Castro::advance() 1 24.42 24.44 24.48 84.96% Castro::subcycle_advance_ctu() 1 24.32 24.37 24.42 84.73% Castro::do_advance_ctu() 20 24.24 24.32 24.38 84.59% Castro::construct_ctu_hydro_source() 20 12.91 13.19 13.57 47.09% Castro::clean_state() 79 3.885 5.452 7.237 25.11% check_for_negative_density() 20 0.4393 3.673 6.041 20.96% FabArray::ParallelCopy() 1429 1.492 3.445 5.112 17.74% FabArray::ParallelCopy_finish() 1429 1.037 3.096 4.857 16.85% Castro::computeTemp() 80 2.493 3.515 4.649 16.13% AmrLevel::FillPatch() 147 1.47 3.007 4.294 14.90% FillPatchIterator::Initialize 147 1.15 2.757 4.108 14.25% FillPatchIterator::FillFromLevel0() 147 1.109 2.731 4.078 14.15% FillPatchSingleLevel 147 1.108 2.73 4.077 14.15% Gravity::solve_phi_with_mlmg() 19 3.231 3.628 3.945 13.69% Castro::construct_new_gravity() 17 3.09 3.483 3.801 13.19% Gravity::solve_for_phi() 17 3.02 3.416 3.732 12.95% Castro::initialize_do_advance() 20 2.7 3.241 3.61 12.53% Castro::reset_internal_energy(MultiFab) 80 1.579 2.143 2.843 9.86% Gravity::actual_solve_with_mlmg() 19 2.733 2.741 2.75 9.54% MLMG::solve() 18 2.645 2.662 2.673 9.27% Castro::do_old_sources() 20 2.198 2.318 2.414 8.38% MLMG::oneIter() 128 2.303 2.33 2.36 8.19% MLMG::mgVcycle() 128 2.287 2.309 2.333 8.10% Amr::InitAmr() 1 2.134 2.134 2.134 7.41% Amr::init() 1 1.976 1.976 1.976 6.86% Castro::do_new_sources() 17 0.9708 1.39 1.896 6.58% Castro::evaluate_source_change() 196 1.004 1.341 1.748 6.07% MLCellLinOp::smooth() 3072 1.265 1.374 1.625 5.64% Amr::initialInit() 1 1.564 1.564 1.564 5.43% amrex::Copy() 724 0.8153 1.105 1.458 5.06% Castro::enforce_min_density() 79 0.836 1.102 1.431 4.97% Castro::expand_state() 20 0.362 0.9056 1.319 4.58% MLCellLinOp::applyBC() 7148 0.72 1.008 1.259 4.37% Castro::normalize_species() 79 0.5789 0.8735 1.256 4.36% Gravity::fill_multipole_BCs() 19 0.4909 0.8859 1.211 4.20% FabArray::FillBoundary() 6380 0.596 0.8987 1.163 4.04% Castro::create_source_corrector() 17 0.1806 0.7007 1.128 3.91% Amr::FinalizeInit() 1 0.7062 0.8765 1.111 3.85% Castro::reset_internal_energy(Fab) 160 0.5456 0.7859 1.086 3.77% FillBoundary_finish() 6380 0.3643 0.7022 1.026 3.56% MultiFab::AddProduct() 3439 0.5785 0.7655 0.9914 3.44% MLMG::mgVcycle_up::1 128 0.1962 0.533 0.9111 3.16% Amr::InitializeInit() 1 0.4529 0.6871 0.8574 2.98% Amr::defBaseLevel() 1 0.4529 0.6871 0.8574 2.98% MLPoisson::Fsmooth() 6144 0.3624 0.554 0.8551 2.97% MLMG::addInterpCorrection() 768 0.1161 0.5934 0.7691 2.67% Castro::Castro() 1 0.3708 0.5445 0.7062 2.45% Castro::initMFs() 1 0.3653 0.5377 0.6979 2.42% Castro::post_init() 1 0.5343 0.5734 0.6131 2.13% MLMG::mgVcycle_down::0 128 0.5075 0.5481 0.5922 2.06% MLMG::mgVcycle_up::2 128 0.0141 0.236 0.5839 2.03% MultiFab::Subtract() 159 0.3618 0.4571 0.5741 1.99% FabArray::setVal() 5013 0.3733 0.4659 0.5713 1.98% Amr::writePlotFile() 2 0.5018 0.502 0.5021 1.74% FabArray::ParallelCopy_nowait() 1429 0.2305 0.3447 0.4676 1.62% Amr::bldFineLevels() 1 0.08906 0.2393 0.4535 1.57% MLMG::mgVcycle_up::0 128 0.4208 0.4338 0.4478 1.55% MLCellLinOp::apply() 986 0.3381 0.3594 0.4147 1.44% Amr::grid_places() 2 0.004797 0.1751 0.4092 1.42% AmrMesh::MakeNewGrids() 2 0.004715 0.175 0.4092 1.42% FabArray::sum() 3439 0.1986 0.2908 0.4082 1.42% TagBoxArray::mapPRD 2 0.002017 0.172 0.4068 1.41% Castro::construct_old_gravity() 20 0.142 0.2727 0.3764 1.31% Gravity::get_old_grav_vector() 20 0.1419 0.2726 0.3763 1.31% MLMG::mgVcycle_down::1 128 0.2732 0.3134 0.3607 1.25% MLMG::mgVcycle_up::3 128 0.01401 0.0428 0.3438 1.19% Castro::sum_integrated_quantities() 2 0.04127 0.3078 0.3256 1.13% Castro::problem_post_timestep() 2 0.01488 0.1037 0.3189 1.11% MLMG:computeResOfCorrection() 768 0.1992 0.24 0.3016 1.05% MLCellLinOp::correctionResidual() 768 0.1977 0.238 0.2991 1.04% Castro::construct_new_source() 102 0.1313 0.1893 0.2587 0.90% Castro::initData() 1 0.08038 0.1405 0.2501 0.87% FillBoundary_nowait() 6380 0.1035 0.1779 0.25 0.87% MLPoisson::Fapply() 986 0.1001 0.1555 0.2322 0.81% MLCellLinOp::solutionResidual() 147 0.1858 0.199 0.2109 0.73% Gravity::multilevel_solve_for_new_phi() 1 0.1899 0.1902 0.1905 0.66% Gravity::actual_multilevel_solve() 1 0.1897 0.1901 0.1904 0.66% Castro::post_timestep() 1 0.1199 0.1604 0.186 0.65% MLMG::computeResidual() 128 0.1465 0.1571 0.169 0.59% Castro::post_regrid() 2 0.08829 0.128 0.167 0.58% Castro::construct_old_source() 120 0.08302 0.1201 0.1645 0.57% Castro::construct_new_gravity_source() 17 0.07205 0.1113 0.1568 0.54% MLMG::mgVcycle_down::2 128 0.01647 0.07668 0.153 0.53% Castro::derive() 18 0.07214 0.1002 0.141 0.49% AmrLevel::derive() 18 0.07209 0.1002 0.1409 0.49% FabArray::Saxpy() 264 0.0716 0.0953 0.1319 0.46% Castro::wd_update() 2 0.01486 0.08595 0.1284 0.45% MLMG::ResNormInf() 146 0.05363 0.09225 0.1247 0.43% MLMG::mgVcycle_bottom 128 0.0003745 0.00817 0.1239 0.43% MLMG::mgVcycle_down::3 128 0.01571 0.03217 0.1236 0.43% MLMG::actualBottomSolve() 128 0.0001423 0.007881 0.1235 0.43% MLCGSolver::bicgstab 8 0 0.007569 0.1211 0.42% FabArray::Xpay() 945 0.05344 0.07728 0.1118 0.39% Castro::construct_old_gravity_source() 20 0.04779 0.07305 0.1032 0.36% amrex::average_down 768 0.03533 0.06545 0.1025 0.36% Castro::initialize_advance() 1 0.0515 0.06971 0.09721 0.34% Amr::regrid() 1 0.04615 0.06612 0.08619 0.30% Castro::check_for_nan() 37 0.04782 0.0635 0.08196 0.28% MultiFab::contains_nan() 37 0.04768 0.06332 0.08176 0.28% Gravity::update_max_rhs() 2 0.002957 0.04205 0.08172 0.28% MultiFab::max() 2 0.0002903 0.04006 0.08027 0.28% Castro::retry_advance_ctu() 11 0.03719 0.04837 0.076 0.26% Gravity::get_new_grav_vector() 18 0.06302 0.06733 0.07086 0.25% MLMG::mgVcycle_down::4 128 0.01563 0.02303 0.06854 0.24% MLPoisson::define() 19 0.0402 0.05225 0.06585 0.23% StateDataPhysBCFunct::() 147 0.001112 0.02576 0.06529 0.23% StateData::FillBoundary(geom) 143 0 0.02448 0.06362 0.22% MLMG::computeMLResidual() 18 0.02156 0.04062 0.05696 0.20% MLMG::mgVcycle_up::4 128 0.01415 0.01991 0.05395 0.19% VisMF::Write(FabArray) 4 0.01411 0.02065 0.053 0.18% Castro::apply_source_to_state() 37 0.02926 0.03865 0.04952 0.17% Castro::save_data_for_retry() 3 0.01944 0.02615 0.04584 0.16% MLMG::mgVcycle_down::5 128 0.01563 0.0213 0.04447 0.15% MLCellLinOp::defineAuxData() 19 0.02048 0.02902 0.04395 0.15% MLMG::prepareForSolve() 18 0.02131 0.02798 0.03867 0.13% MLMG::mgVcycle_up::5 128 0.01427 0.01894 0.0377 0.13% FabArray::norminf() 282 0.01231 0.0173 0.0294 0.10% Gravity::test_composite_phi() 1 0.02664 0.02696 0.02727 0.09% amrex::Add() 128 0.01461 0.02003 0.02605 0.09% FabArray::setDomainBndry() 147 0.0006515 0.009734 0.023 0.08% Other 10861 0.1221 0.1857 0.2823 0.98% -------------------------------------------------------------------------------------------- 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 6 B 7 B 8 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 1594 KiB 1595 KiB 1595 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 221 KiB 307 KiB 426 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1063 KiB 1063 KiB 1063 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 1011 KiB 1012 KiB 1012 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 674 KiB 675 KiB 675 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 2230 B 3116 B 4656 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 26 KiB 35 KiB 47 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 13 KiB 28 KiB 40 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 421 KiB 422 KiB 423 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 2765 B 12 KiB 19 KiB 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 695 B 962 B 1372 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 9455 B 12 KiB 15 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 133 KiB 133 KiB 134 KiB 157 KiB 157 KiB 157 KiB MLMG::prepareForSolve() 6102 6102 9647 B 10 KiB 11 KiB 102 KiB 108 KiB 126 KiB Castro::buildMetrics() 160 160 108 KiB 108 KiB 108 KiB 117 KiB 117 KiB 117 KiB Castro::construct_old_source() 640 640 135 B 182 B 238 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 223 B 305 B 410 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 9 B 12 B 26 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 33 KiB 34 KiB 34 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 48 B 49 B 50 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 3 B 24 B 47 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 99 B 106 B 112 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 226 B 435 B 602 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 64 33 KiB 33 KiB 33 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 2432 175 B 208 B 244 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 567 B 758 B 989 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 107 B 107 B 108 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1715 B 1941 B 2121 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 433 B 484 B 673 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 452 B 453 B 455 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 34 B 72 B 102 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 15 B 36 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 26 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 3 B 9 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-16-gc0602940a47c) finalized