Initializing AMReX (24.05)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.05) initialized Starting run at 21:38:43 UTC on 2024-05-01. Successfully read inputs file ... Castro git describe: 24.05-1-g7d1a1c16d AMReX git describe: 24.05 Microphysics git describe: 24.05 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 grid_places() time: 0.175520093 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001090808 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010527789 Gravity::multilevel_solve_for_new_phi() time = 0.099897724 Gravity::fill_multipole_BCs() time = 0.010614254 ... norm of composite residual at level 0 1.193955541e-11 TIME= 0 MASS = 2.549996568e+33 TIME= 0 XMOM = 0 TIME= 0 YMOM = 0 TIME= 0 ZMOM = 0 TIME= 0 ANG MOM X = 0 TIME= 0 ANG MOM Y = 0 TIME= 0 ANG MOM Z = -3.115378115e+34 TIME= 0 RHO*e = 1.325915365e+50 TIME= 0 RHO*K = 3.09960185e+49 TIME= 0 RHO*E = 1.63587555e+50 TIME= 0 RHO*PHI = -5.434813569e+50 TIME= 0 TOTAL ENERGY= -1.081531234e+50 TIME= 0 CENTER OF MASS X-LOC = 316439254.4 TIME= 0 CENTER OF MASS X-VEL = 0 TIME= 0 CENTER OF MASS Y-LOC = -4.854673355e-07 TIME= 0 CENTER OF MASS Y-VEL = 0 TIME= 0 CENTER OF MASS Z-LOC = 0 TIME= 0 CENTER OF MASS Z-VEL = 0 TIME= 0 MAXIMUM TEMPERATURE = 10000000 TIME= 0 MAXIMUM DENSITY = 4476474 TIME= 0 MAXIMUM T_S / T_E = 0 INITIAL GRIDS Level 0 32 grids 32768 cells 100 % of domain smallest grid: 32 x 32 biggest grid: 32 x 32 PLOTFILE: file = wdmerger_retry-SDC_plt00000 Write plotfile time = 0.048837826 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.018856301 seconds [Level 0 step 1] ADVANCE at time 0 with dt = 0.25  Beginning subcycle 1 starting at time 0 with dt = 0.25 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. Castro::construct_old_gravity() time = 0.002717459 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.065008611 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.767071355 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.767071355 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071149 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477042 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035501 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035501 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477042 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.7694089 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071149 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. Castro::construct_old_gravity() time = 0.004775788 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.064049969 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.125 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.0625   Beginning subcycle 1 starting at time 0 with dt = 0.0625 Estimated number of subcycles remaining: 4 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. Castro::construct_old_gravity() time = 0.00365421 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.060473243 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.1060744788 at index 20, 78, 0 Invalid density = -0.1060744788 at index 20, 177, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.0625 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.03125   Beginning subcycle 1 starting at time 0 with dt = 0.03125 Estimated number of subcycles remaining: 8 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. Castro::construct_old_gravity() time = 0.005768549 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.058448457 on level 0 ... Entering construct_ctu_hydro_source() on level 0 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.317411453 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014122761 Gravity::solve_for_phi() time = 0.083444492 on level 0 Castro::construct_new_gravity() time = 0.085577386 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.059092266 on level 0 Ending SDC iteration 1 of 3. Beginning SDC iteration 2 of 3. Castro::construct_old_gravity() time = 0.002518629 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.058131469 on level 0 ... Entering construct_ctu_hydro_source() on level 0 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.316911639 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014069712 Gravity::solve_for_phi() time = 0.083451266 on level 0 Castro::construct_new_gravity() time = 0.085562585 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.058957486 on level 0 Ending SDC iteration 2 of 3. Beginning SDC iteration 3 of 3. Castro::construct_old_gravity() time = 0.002832631 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.058362405 on level 0 ... Entering construct_ctu_hydro_source() on level 0 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.32157459 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01377285 Gravity::solve_for_phi() time = 0.083038143 on level 0 Castro::construct_new_gravity() time = 0.085199206 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.05895581 on level 0 Ending SDC iteration 3 of 3.  Subcycle completed  Beginning subcycle 2 starting at time 0.03125 with dt = 0.03125 Estimated number of subcycles remaining: 7 Beginning SDC iteration 1 of 2. Castro::construct_old_gravity() time = 0.002520339 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.058131617 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.985598081e+46 ener added: 2.745852676e+42 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.325465234 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014878854 Gravity::solve_for_phi() time = 0.084524543 on level 0 Castro::construct_new_gravity() time = 0.086558848 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.274072206e+44 ener added: 4.891864863e+43 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.763054894e+35 ymom added: -3.010508633e+24 zmom added: 0 eint added: 0 ener added: 5.968608016e+43 Castro::do_new_sources() time = 0.058898666 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.002929949 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.058100978 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.986297782e+46 ener added: 2.752525031e+42 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.330696426 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014505044 Gravity::solve_for_phi() time = 0.083793301 on level 0 Castro::construct_new_gravity() time = 0.085959117 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.347573489e+44 ener added: 4.888153391e+43 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.756949739e+35 ymom added: -4.309159416e+24 zmom added: 0 eint added: 0 ener added: 5.2374099e+43 Castro::do_new_sources() time = 0.059022481 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 3 starting at time 0.0625 with dt = 0.03125 Estimated number of subcycles remaining: 6 Beginning SDC iteration 1 of 2. Castro::construct_old_gravity() time = 0.002590349 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.984353173e+46 ener added: 6.990361455e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46774173e+40 ymom added: 9.528602194e+26 zmom added: 0 eint added: 0 ener added: 8.780735352e+46 Castro::do_old_sources() time = 0.058063225 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.98733239e+46 ener added: 2.180535114e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.396494265 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.017180516 Gravity::solve_for_phi() time = 0.086335936 on level 0 Castro::construct_new_gravity() time = 0.088458012 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.058042977e+43 ener added: 1.184398423e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.782398841e+35 ymom added: 1.223092919e+24 zmom added: 0 eint added: 0 ener added: 5.96849108e+43 Castro::do_new_sources() time = 0.058833578 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.002687167 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.984353173e+46 ener added: 6.990361455e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46774173e+40 ymom added: 9.528602194e+26 zmom added: 0 eint added: 0 ener added: 8.780735352e+46 Castro::do_old_sources() time = 0.058375394 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.988039066e+46 ener added: 2.185480971e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.330814791 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014534566 Gravity::solve_for_phi() time = 0.083774321 on level 0 Castro::construct_new_gravity() time = 0.085811904 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.800848844e+43 ener added: 1.183591058e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.776272322e+35 ymom added: 1.448585919e+24 zmom added: 0 eint added: 0 ener added: 5.233809214e+43 Castro::do_new_sources() time = 0.058886003 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 4 starting at time 0.09375 with dt = 0.03125 Estimated number of subcycles remaining: 5 Beginning SDC iteration 1 of 2. Castro::construct_old_gravity() time = 0.002692996 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.983917922e+46 ener added: 1.154265495e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467777256e+40 ymom added: 9.528602194e+26 zmom added: 0 eint added: 0 ener added: 8.841018465e+46 Castro::do_old_sources() time = 0.058569048 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.98498831e+46 ener added: 2.292089986e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.325476353 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01465826 Gravity::solve_for_phi() time = 0.083666665 on level 0 Castro::construct_new_gravity() time = 0.085695815 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -4.614759295e+43 ener added: 1.488589688e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.803591267e+35 ymom added: -6.871043233e+23 zmom added: 0 eint added: 0 ener added: 4.637081169e+43 Castro::do_new_sources() time = 0.059422447 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.00281359 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.983917922e+46 ener added: 1.154265495e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467777256e+40 ymom added: 9.528602194e+26 zmom added: 0 eint added: 0 ener added: 8.841018465e+46 Castro::do_old_sources() time = 0.057845073 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.985698801e+46 ener added: 2.298024982e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.329892628 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014567466 Gravity::solve_for_phi() time = 0.083804555 on level 0 Castro::construct_new_gravity() time = 0.085907684 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -5.359798244e+43 ener added: 1.487684259e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.797448342e+35 ymom added: -5.254813304e+24 zmom added: 0 eint added: 0 ener added: 3.901081726e+43 Castro::do_new_sources() time = 0.059189148 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 5 starting at time 0.125 with dt = 0.03125 Estimated number of subcycles remaining: 4 Beginning SDC iteration 1 of 2. Castro::construct_old_gravity() time = 0.002615002 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.98356532e+46 ener added: 1.521222674e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467813205e+40 ymom added: 9.429621393e+26 zmom added: 0 eint added: 0 ener added: 8.902573838e+46 Castro::do_old_sources() time = 0.058011287 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.992363162e+46 ener added: 2.321339475e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.325994454 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014942345 Gravity::solve_for_phi() time = 0.084366592 on level 0 Castro::construct_new_gravity() time = 0.08640478 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.068921344e+43 ener added: 1.620777103e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.829776437e+35 ymom added: 5.133212367e+24 zmom added: 0 eint added: 0 ener added: 2.123241032e+43 Castro::do_new_sources() time = 0.059828788 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.002686941 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.98356532e+46 ener added: 1.521222674e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467813205e+40 ymom added: 9.429621393e+26 zmom added: 0 eint added: 0 ener added: 8.902573838e+46 Castro::do_old_sources() time = 0.058326862 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.993083003e+46 ener added: 2.326773931e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.331014342 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014757266 Gravity::solve_for_phi() time = 0.083713456 on level 0 Castro::construct_new_gravity() time = 0.085829265 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.821028634e+43 ener added: 1.619931887e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.823603759e+35 ymom added: 8.029203612e+24 zmom added: 0 eint added: 0 ener added: 1.379533843e+43 Castro::do_new_sources() time = 0.058806685 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 6 starting at time 0.15625 with dt = 0.03125 Estimated number of subcycles remaining: 3 Beginning SDC iteration 1 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 4.920709246e+29 ener added: 4.920709246e+29 Castro::construct_old_gravity() time = 0.002703291 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.983073759e+46 ener added: 2.019352467e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467849677e+40 ymom added: 9.585270592e+26 zmom added: 0 eint added: 0 ener added: 8.96490842e+46 Castro::do_old_sources() time = 0.058306692 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.000609588e+46 ener added: 2.406064512e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.328453831 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014733463 Gravity::solve_for_phi() time = 0.083918069 on level 0 Castro::construct_new_gravity() time = 0.086045256 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -8.854195607e+43 ener added: 1.71395957e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.848617495e+35 ymom added: -5.170991299e+24 zmom added: 0 eint added: 0 ener added: -1.366537222e+43 Castro::do_new_sources() time = 0.058971496 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 4.920709246e+29 ener added: 4.920709246e+29 Castro::construct_old_gravity() time = 0.002712261 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.983073759e+46 ener added: 2.019352467e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467849677e+40 ymom added: 9.585270592e+26 zmom added: 0 eint added: 0 ener added: 8.96490842e+46 Castro::do_old_sources() time = 0.058668403 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.001333674e+46 ener added: 2.41251015e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.346270702 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01678217 Gravity::solve_for_phi() time = 0.090634529 on level 0 Castro::construct_new_gravity() time = 0.092741579 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -9.610031177e+43 ener added: 1.713054017e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.842414694e+35 ymom added: -7.880449068e+24 zmom added: 0 eint added: 0 ener added: -2.113374984e+43 Castro::do_new_sources() time = 0.061569509 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 7 starting at time 0.1875 with dt = 0.03125 Estimated number of subcycles remaining: 2 Beginning SDC iteration 1 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.075345155e+29 ener added: 2.075345155e+29 Castro::construct_old_gravity() time = 0.002619578 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.982914566e+46 ener added: 2.174746794e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467886525e+40 ymom added: 9.439443915e+26 zmom added: 0 eint added: 0 ener added: 9.027789929e+46 Castro::do_old_sources() time = 0.058285144 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.993854135e+46 ener added: 2.676087223e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.326071985 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01508361 Gravity::solve_for_phi() time = 0.084372021 on level 0 Castro::construct_new_gravity() time = 0.086907641 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.152849071e+44 ener added: 1.743211419e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.877263956e+35 ymom added: 7.059937892e+24 zmom added: 0 eint added: 0 ener added: -4.049549507e+43 Castro::do_new_sources() time = 0.059149191 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.075345155e+29 ener added: 2.075345155e+29 Castro::construct_old_gravity() time = 0.003214893 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.982914566e+46 ener added: 2.174746794e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467886525e+40 ymom added: 9.439443915e+26 zmom added: 0 eint added: 0 ener added: 9.027789929e+46 Castro::do_old_sources() time = 0.058140733 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.994586452e+46 ener added: 2.682600444e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.332649031 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015162194 Gravity::solve_for_phi() time = 0.084711796 on level 0 Castro::construct_new_gravity() time = 0.086843237 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.229235664e+44 ener added: 1.74231031e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.871026436e+35 ymom added: 8.771795742e+23 zmom added: 0 eint added: 0 ener added: -4.804462625e+43 Castro::do_new_sources() time = 0.058745672 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 8 starting at time 0.21875 with dt = 0.03125 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.383563436e+30 ener added: 1.383563436e+30 Castro::construct_old_gravity() time = 0.002595405 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.982548267e+46 ener added: 2.529322566e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467923946e+40 ymom added: 9.471933797e+26 zmom added: 0 eint added: 0 ener added: 9.091353867e+46 Castro::do_old_sources() time = 0.058185035 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.993555556e+46 ener added: 2.937230463e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.328576194 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015234591 Gravity::solve_for_phi() time = 0.084460572 on level 0 Castro::construct_new_gravity() time = 0.086569859 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -8.731399969e+43 ener added: 1.766195089e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.939762636e+35 ymom added: 1.676440101e+23 zmom added: 0 eint added: 0 ener added: -1.225922477e+43 Castro::do_new_sources() time = 0.058942529 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.383563436e+30 ener added: 1.383563436e+30 Castro::construct_old_gravity() time = 0.002691848 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.982548267e+46 ener added: 2.529322566e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467923946e+40 ymom added: 9.471933797e+26 zmom added: 0 eint added: 0 ener added: 9.091353867e+46 Castro::do_old_sources() time = 0.058214997 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.994288712e+46 ener added: 2.943682879e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.33200881 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015125747 Gravity::solve_for_phi() time = 0.086263402 on level 0 Castro::construct_new_gravity() time = 0.08836929 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -9.496809054e+43 ener added: 1.76531212e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.933525523e+35 ymom added: 5.748300601e+24 zmom added: 0 eint added: 0 ener added: -1.982551467e+43 Castro::do_new_sources() time = 0.058888415 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002793100299 [Level 0 step 1] Advanced 32768 cells Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 8.037844726e+29 ener added: 8.037844726e+29 TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -4.340632835e+38 TIME= 0.25 YMOM = 1.156700224e+28 TIME= 0.25 ZMOM = 0 TIME= 0.25 ANG MOM X = 0 TIME= 0.25 ANG MOM Y = 0 TIME= 0.25 ANG MOM Z = -1.469523857e+38 TIME= 0.25 RHO*e = 1.325894499e+50 TIME= 0.25 RHO*K = 3.170941995e+49 TIME= 0.25 RHO*E = 1.642988698e+50 TIME= 0.25 RHO*PHI = -5.449011351e+50 TIME= 0.25 TOTAL ENERGY= -1.081516977e+50 TIME= 0.25 CENTER OF MASS X-LOC = 316424375.9 TIME= 0.25 CENTER OF MASS X-VEL = -170221.1246 TIME= 0.25 CENTER OF MASS Y-LOC = 1.213668339e-07 TIME= 0.25 CENTER OF MASS Y-VEL = 4.536085416e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 58251090.19 TIME= 0.25 MAXIMUM DENSITY = 4458791.937 TIME= 0.25 MAXIMUM T_S / T_E = 0 [STEP 1] Coarse TimeStep time: 11.7906367 [STEP 1] FAB kilobyte spread across MPI nodes: [6973 ... 6973] STEP = 1 TIME = 0.25 DT = 0.25 PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.086210177 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.018644649 seconds Ending run at 21:38:57 UTC on 2024-05-01. Run time = 14.18795917 Run time without initialization = 11.89638021 Average number of zones advanced per microsecond: 0.003 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 14.19 ... 14.19 ... 14.19 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.4154 5.4654 5.5975 39.45% Amr::InitAmr() 1 1.375 1.375 1.375 9.69% Castro::computeTemp() 80 0.8455 0.9589 1.115 7.86% amrex::Copy() 724 0.7686 0.7733 0.7861 5.54% FabArray::ParallelCopy_finish() 1429 0.05851 0.4769 0.6875 4.85% MultiFab::AddProduct() 3439 0.5586 0.5603 0.5672 4.00% Castro::normalize_species() 79 0.541 0.5441 0.5515 3.89% Castro::reset_internal_energy(Fab) 160 0.5232 0.5267 0.5324 3.75% MLPoisson::Fsmooth() 6144 0.3155 0.34 0.4157 2.93% Castro::initMFs() 1 0.3646 0.37 0.391 2.76% FabArray::setVal() 5013 0.3432 0.3468 0.355 2.50% MultiFab::Subtract() 159 0.3323 0.3333 0.3376 2.38% Gravity::fill_multipole_BCs() 19 0.1948 0.2353 0.2719 1.92% check_for_negative_density() 20 0.04267 0.1943 0.2481 1.75% FabArray::ParallelCopy_nowait() 1429 0.1825 0.199 0.2182 1.54% Castro::reset_internal_energy(MultiFab) 80 0.02275 0.06178 0.2029 1.43% FabArray::sum() 3439 0.19 0.1907 0.1921 1.35% Castro::initData() 1 0.003278 0.02862 0.1701 1.20% FillBoundary_nowait() 6380 0.07689 0.1028 0.1301 0.92% FillBoundary_finish() 6380 0.09268 0.1082 0.1258 0.89% MLPoisson::Fapply() 986 0.0951 0.09939 0.1151 0.81% MLCellLinOp::applyBC() 7148 0.04607 0.06125 0.1032 0.73% FabArray::Saxpy() 264 0.06824 0.06888 0.07634 0.54% Castro::do_new_sources() 17 0.0003664 0.01277 0.07214 0.51% Castro::construct_new_gravity_source() 17 0.06947 0.07022 0.0717 0.51% VisMF::Write(FabArray) 4 0.01199 0.01628 0.06085 0.43% FabArray::Xpay() 945 0.05026 0.05245 0.05994 0.42% AmrLevel::derive() 18 0.04075 0.04224 0.05685 0.40% Amr::writePlotFile() 2 0.02309 0.05254 0.05674 0.40% MultiFab::contains_nan() 37 0.04534 0.04599 0.04704 0.33% Castro::enforce_min_density() 79 0.01718 0.02423 0.04596 0.32% Castro::construct_old_gravity_source() 20 0.0447 0.045 0.04541 0.32% MLMG::addInterpCorrection() 768 0.03318 0.03553 0.04274 0.30% Castro::do_old_sources() 20 0.0005266 0.008714 0.04183 0.29% Amr::writeSmallPlotFile() 2 0.01838 0.0327 0.03491 0.25% amrex::average_down 768 0.02429 0.02624 0.03236 0.23% StateData::FillBoundary(geom) 143 0 0.01663 0.0319 0.22% MLCellLinOp::defineAuxData() 19 0.01447 0.01534 0.01744 0.12% FabArray::FillBoundary() 6380 0.01371 0.01428 0.01728 0.12% MLCGSolver::bicgstab 8 0 0.001029 0.01647 0.12% Castro::construct_new_source() 102 0.0156 0.01568 0.01586 0.11% FabArray::norminf() 282 0.01096 0.0113 0.01562 0.11% amrex::Add() 128 0.01396 0.01407 0.01426 0.10% Castro::evaluate_source_change() 196 0.01256 0.01291 0.01403 0.10% MLCellLinOp::smooth() 3072 0.01124 0.01171 0.01338 0.09% FillPatchIterator::Initialize 147 0.0124 0.01262 0.01302 0.09% FabArrayBase::getFB() 6380 0.009676 0.01011 0.01301 0.09% Castro::construct_old_gravity() 20 0.0003789 0.003083 0.01281 0.09% Castro::wd_update() 2 0.005128 0.008771 0.01188 0.08% MLMG::prepareForSolve() 18 0.01007 0.01028 0.01107 0.08% Castro::construct_old_source() 120 0.009923 0.01002 0.01016 0.07% Gravity::actual_solve_with_mlmg() 19 0.006732 0.007016 0.007947 0.06% FabArrayBase::getCPC() 2005 0.007103 0.007359 0.007814 0.06% MLMG::ResNormInf() 146 0.002853 0.005644 0.007022 0.05% MLMG::solve() 18 0.0009559 0.002653 0.006783 0.05% Gravity::get_old_grav_vector() 20 0.00635 0.006468 0.006774 0.05% Gravity::get_new_grav_vector() 18 0.005981 0.006167 0.006611 0.05% Castro::initialize_do_advance() 20 0.004424 0.00464 0.00483 0.03% MLCellLinOp::apply() 986 0.002021 0.002236 0.004747 0.03% FabArray::ParallelCopy() 1429 0.003545 0.003701 0.004022 0.03% AmrLevel::FillPatch() 147 0.002143 0.00223 0.002516 0.02% Castro::sum_integrated_quantities() 2 0.00159 0.002113 0.002357 0.02% MLLinOp::defineGrids() 19 0.001776 0.001873 0.00214 0.02% MLCellLinOp::correctionResidual() 768 0.001817 0.00187 0.001979 0.01% MLMG::mgVcycle() 128 0.001525 0.001596 0.001771 0.01% Castro::construct_new_gravity() 17 0.0003013 0.0005584 0.00168 0.01% MLMG:computeResOfCorrection() 768 0.001368 0.00142 0.001512 0.01% Castro::save_data_for_retry() 3 0.001262 0.001353 0.001502 0.01% Gravity::solve_for_phi() 17 0.0008183 0.0009805 0.001476 0.01% main() 1 0.0009811 0.001177 0.001271 0.01% MLMG::actualBottomSolve() 128 0.0001284 0.0001955 0.001083 0.01% StateDataPhysBCFunct::() 147 0.0006632 0.0008156 0.001011 0.01% Castro::subcycle_advance_ctu() 1 0.0001877 0.0002567 0.0007956 0.01% FillPatchIterator::FillFromLevel0() 147 0.0006852 0.000712 0.0007816 0.01% MLMG::mgVcycle_down::2 128 0.0006179 0.0006656 0.0007621 0.01% MLMG::mgVcycle_down::0 128 0.00065 0.0006641 0.0007249 0.01% MLMG::mgVcycle_down::3 128 0.0005745 0.0006223 0.0007189 0.01% MLMG::mgVcycle_down::1 128 0.0006442 0.0006573 0.000713 0.01% MLMG::mgVcycle_down::4 128 0.0005582 0.0005823 0.0006956 0.00% MLMG::mgVcycle_down::5 128 0.0005817 0.000599 0.0006666 0.00% Castro::Castro() 1 0.0005403 0.0005618 0.000637 0.00% Castro::do_advance_ctu() 20 0.0004094 0.0004556 0.0005824 0.00% AmrMesh::MakeNewGrids() 2 0.0005125 0.0005271 0.0005483 0.00% MLMG::mgVcycle_up::1 128 0.0004521 0.0004735 0.0005272 0.00% MLMG::mgVcycle_up::0 128 0.0004647 0.0004762 0.0005229 0.00% Castro::clean_state() 79 0.0004375 0.0004606 0.0005065 0.00% FillPatchSingleLevel 147 0.0004675 0.0004807 0.000506 0.00% MLMG::mgVcycle_up::2 128 0.0003877 0.0004266 0.0004878 0.00% MLCellLinOp::solutionResidual() 147 0.0004445 0.0004574 0.0004827 0.00% MLMG::mgVcycle_up::3 128 0.0003862 0.0004031 0.0004576 0.00% MLMG::mgVcycle_up::5 128 0.0003982 0.0004104 0.0004534 0.00% MLMG::oneIter() 128 0.0003839 0.0003986 0.0004453 0.00% MLMG::mgVcycle_up::4 128 0.0003895 0.0004016 0.0004436 0.00% Castro::problem_post_timestep() 2 1.09e-05 3.448e-05 0.00037 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002129 0.0002322 0.0002543 0.00% MLMG::mgVcycle_bottom 128 0.0002281 0.0002371 0.0002527 0.00% MLMG::computeResidual() 128 0.0002326 0.0002398 0.000248 0.00% MLPoisson::define() 19 0.0001319 0.000156 0.0001929 0.00% Castro::initialize_advance() 1 0.0001158 0.0001222 0.0001556 0.00% Amr::coarseTimeStep() 1 0.0001249 0.0001468 0.0001539 0.00% Castro::retry_advance_ctu() 11 6.687e-05 9.002e-05 0.000143 0.00% Amr::defBaseLevel() 1 0.0001055 0.0001186 0.0001371 0.00% Castro::check_for_nan() 37 0.0001034 0.0001129 0.0001276 0.00% TagBoxArray::mapPRD 2 0.0001117 0.000118 0.0001255 0.00% Amr::grid_places() 2 8.658e-05 9.943e-05 0.0001242 0.00% Gravity::test_composite_phi() 1 9.28e-05 9.962e-05 0.0001216 0.00% Castro::apply_source_to_state() 37 9.133e-05 0.0001003 0.0001094 0.00% MLLinOp::define() 19 8.007e-05 8.861e-05 0.0001005 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.058e-05 6.892e-05 9.262e-05 0.00% Castro::create_source_corrector() 17 5.816e-05 6.274e-05 7.344e-05 0.00% Castro::expand_state() 20 5.441e-05 5.767e-05 6.746e-05 0.00% Amr::FinalizeInit() 1 1.099e-05 1.845e-05 6.18e-05 0.00% Castro::derive() 18 4.155e-05 4.457e-05 5.524e-05 0.00% MLMG::computeMLResidual() 18 3.979e-05 4.19e-05 4.52e-05 0.00% Amr::regrid() 1 1.051e-05 1.267e-05 3.134e-05 0.00% Gravity::actual_multilevel_solve() 1 2.649e-05 2.786e-05 2.96e-05 0.00% Amr::timeStep() 1 1.179e-05 1.631e-05 2.365e-05 0.00% Castro::post_init() 1 1.202e-05 1.314e-05 1.74e-05 0.00% Castro::post_regrid() 2 1.128e-05 1.214e-05 1.474e-05 0.00% Amr::InitializeInit() 1 6.714e-06 7.914e-06 1.224e-05 0.00% Castro::advance() 1 6.979e-06 8.315e-06 1.191e-05 0.00% Amr::init() 1 6.143e-06 7.562e-06 1.001e-05 0.00% Amr::bldFineLevels() 1 6.951e-06 7.68e-06 8.78e-06 0.00% Castro::post_timestep() 1 5.644e-06 6.128e-06 6.831e-06 0.00% Amr::initialInit() 1 2.941e-06 3.481e-06 6.663e-06 0.00% Other 2587 0.04344 0.05915 0.0989 0.70% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.19 14.19 14.19 100.00% Amr::coarseTimeStep() 1 11.9 11.9 11.9 83.84% Amr::timeStep() 1 11.79 11.79 11.79 83.10% Castro::advance() 1 11.73 11.73 11.73 82.69% Castro::subcycle_advance_ctu() 1 11.68 11.68 11.68 82.34% Castro::do_advance_ctu() 20 11.64 11.64 11.64 82.07% Castro::construct_ctu_hydro_source() 20 6.405 6.434 6.471 45.61% Castro::clean_state() 79 3.662 3.766 3.909 27.55% Castro::computeTemp() 80 2.374 2.491 2.699 19.02% Castro::reset_internal_energy(MultiFab) 80 1.488 1.532 1.674 11.80% Gravity::solve_phi_with_mlmg() 19 1.465 1.508 1.547 10.90% Castro::construct_new_gravity() 17 1.393 1.434 1.473 10.38% Gravity::solve_for_phi() 17 1.355 1.398 1.437 10.13% Castro::initialize_do_advance() 20 1.351 1.384 1.413 9.96% Amr::InitAmr() 1 1.375 1.375 1.375 9.69% Gravity::actual_solve_with_mlmg() 19 1.266 1.269 1.274 8.98% MLMG::solve() 18 1.215 1.219 1.221 8.61% Castro::do_old_sources() 20 1.132 1.155 1.175 8.28% MLMG::oneIter() 128 1.087 1.088 1.09 7.69% MLMG::mgVcycle() 128 1.073 1.074 1.076 7.58% Castro::do_new_sources() 17 0.9324 0.9719 1.005 7.08% Castro::evaluate_source_change() 196 0.9558 0.9587 0.9697 6.83% Amr::init() 1 0.9159 0.9159 0.9159 6.46% FabArray::ParallelCopy() 1429 0.2839 0.6992 0.9066 6.39% Amr::initialInit() 1 0.8474 0.8475 0.8475 5.97% Castro::enforce_min_density() 79 0.784 0.7942 0.816 5.75% amrex::Copy() 724 0.7686 0.7733 0.7861 5.54% MLCellLinOp::smooth() 3072 0.5429 0.5912 0.7154 5.04% FabArray::ParallelCopy_finish() 1429 0.05851 0.4769 0.6875 4.85% Amr::InitializeInit() 1 0.4503 0.5179 0.6246 4.40% Amr::defBaseLevel() 1 0.4503 0.5179 0.6246 4.40% AmrLevel::FillPatch() 147 0.4769 0.5321 0.608 4.29% MultiFab::AddProduct() 3439 0.5586 0.5603 0.5672 4.00% Castro::normalize_species() 79 0.541 0.5441 0.5515 3.89% Castro::reset_internal_energy(Fab) 160 0.5232 0.5267 0.5324 3.75% FillPatchIterator::Initialize 147 0.2995 0.3562 0.4305 3.03% MLPoisson::Fsmooth() 6144 0.3155 0.34 0.4157 2.93% MLMG::mgVcycle_up::1 128 0.09457 0.2551 0.4151 2.93% FillPatchIterator::FillFromLevel0() 147 0.2758 0.3377 0.4138 2.92% FillPatchSingleLevel 147 0.2751 0.337 0.4131 2.91% Castro::Castro() 1 0.3705 0.376 0.3979 2.80% Amr::FinalizeInit() 1 0.2229 0.3295 0.3972 2.80% Castro::initMFs() 1 0.3648 0.3702 0.3911 2.76% MLMG::addInterpCorrection() 768 0.06252 0.2977 0.3657 2.58% MLCellLinOp::applyBC() 7148 0.2672 0.2973 0.3608 2.54% FabArray::setVal() 5013 0.3432 0.3468 0.355 2.50% MultiFab::Subtract() 159 0.3323 0.3333 0.3376 2.38% Gravity::fill_multipole_BCs() 19 0.1969 0.2374 0.274 1.93% FabArray::FillBoundary() 6380 0.2094 0.236 0.2712 1.91% MLMG::mgVcycle_up::2 128 0.0107 0.1133 0.2711 1.91% MLMG::mgVcycle_down::0 128 0.2584 0.2594 0.2619 1.85% Castro::initData() 1 0.07723 0.1402 0.2523 1.78% check_for_negative_density() 20 0.04267 0.1943 0.2481 1.75% Castro::expand_state() 20 0.178 0.2046 0.2435 1.72% FabArray::ParallelCopy_nowait() 1429 0.2011 0.2186 0.2386 1.68% Amr::bldFineLevels() 1 0.04654 0.1526 0.2205 1.55% MLMG::mgVcycle_up::0 128 0.1931 0.1935 0.1939 1.37% FabArray::sum() 3439 0.19 0.1907 0.1921 1.35% MLCellLinOp::apply() 986 0.1471 0.1577 0.1921 1.35% MLMG::mgVcycle_up::3 128 0.01065 0.02339 0.1766 1.24% Amr::grid_places() 2 0.002354 0.109 0.1766 1.24% AmrMesh::MakeNewGrids() 2 0.002247 0.1089 0.1765 1.24% TagBoxArray::mapPRD 2 0.0007583 0.106 0.1751 1.23% MLMG:computeResOfCorrection() 768 0.1119 0.121 0.1451 1.02% MLCellLinOp::correctionResidual() 768 0.1105 0.1196 0.1436 1.01% FillBoundary_nowait() 6380 0.08693 0.1136 0.1414 1.00% Amr::writePlotFile() 2 0.1356 0.1357 0.1357 0.96% Castro::post_init() 1 0.1326 0.1336 0.1349 0.95% Castro::construct_new_source() 102 0.1237 0.1245 0.1258 0.89% FillBoundary_finish() 6380 0.09268 0.1082 0.1258 0.89% MLMG::mgVcycle_down::1 128 0.1121 0.1157 0.1189 0.84% MLPoisson::Fapply() 986 0.0951 0.09939 0.1151 0.81% Gravity::multilevel_solve_for_new_phi() 1 0.09986 0.09989 0.09992 0.70% Gravity::actual_multilevel_solve() 1 0.09979 0.09982 0.09985 0.70% MLCellLinOp::solutionResidual() 147 0.08564 0.09121 0.0953 0.67% Castro::derive() 18 0.06893 0.07102 0.09148 0.64% AmrLevel::derive() 18 0.06889 0.07098 0.09143 0.64% Castro::post_regrid() 2 0.08563 0.08689 0.08793 0.62% Castro::construct_old_source() 120 0.07733 0.07777 0.07864 0.55% MLMG::computeResidual() 128 0.07422 0.07481 0.07693 0.54% FabArray::Saxpy() 264 0.06824 0.06888 0.07634 0.54% Castro::construct_new_gravity_source() 17 0.06947 0.07022 0.0717 0.51% MLMG::mgVcycle_bottom 128 0.0003616 0.004822 0.07157 0.50% MLMG::actualBottomSolve() 128 0.0001284 0.004585 0.07132 0.50% MLCGSolver::bicgstab 8 0 0.004361 0.06978 0.49% MLMG::mgVcycle_down::2 128 0.01247 0.03976 0.06928 0.49% Castro::post_timestep() 1 0.05882 0.06071 0.06305 0.44% VisMF::Write(FabArray) 4 0.01199 0.01628 0.06085 0.43% Castro::create_source_corrector() 17 0.04746 0.05179 0.06073 0.43% FabArray::Xpay() 945 0.05026 0.05245 0.05994 0.42% Castro::construct_old_gravity() 20 0.04028 0.0519 0.05933 0.42% MLMG::mgVcycle_down::3 128 0.01186 0.01791 0.05905 0.42% Gravity::get_old_grav_vector() 20 0.03893 0.04882 0.05894 0.42% amrex::average_down 768 0.03249 0.03852 0.0542 0.38% Castro::initialize_advance() 1 0.04879 0.04969 0.05122 0.36% Castro::check_for_nan() 37 0.04545 0.0461 0.04716 0.33% MultiFab::contains_nan() 37 0.04534 0.04599 0.04704 0.33% Castro::construct_old_gravity_source() 20 0.0447 0.045 0.04541 0.32% Amr::regrid() 1 0.04403 0.04467 0.04527 0.32% MLPoisson::define() 19 0.02876 0.03349 0.03997 0.28% Castro::retry_advance_ctu() 11 0.03535 0.03645 0.03883 0.27% Amr::writeSmallPlotFile() 2 0.03789 0.03789 0.03789 0.27% Gravity::get_new_grav_vector() 18 0.03483 0.03582 0.03716 0.26% MLMG::mgVcycle_down::4 128 0.01179 0.01338 0.03458 0.24% StateDataPhysBCFunct::() 147 0.0006632 0.01745 0.03274 0.23% StateData::FillBoundary(geom) 143 0 0.01663 0.0319 0.22% Castro::apply_source_to_state() 37 0.02784 0.02795 0.02818 0.20% MLMG::mgVcycle_up::4 128 0.01076 0.01195 0.02729 0.19% MLCellLinOp::defineAuxData() 19 0.01857 0.02015 0.02409 0.17% MLMG::prepareForSolve() 18 0.02013 0.0209 0.02323 0.16% MLMG::mgVcycle_down::5 128 0.0118 0.01264 0.023 0.16% Castro::sum_integrated_quantities() 2 0.02092 0.02119 0.02121 0.15% Castro::save_data_for_retry() 3 0.01859 0.01899 0.02034 0.14% MLMG::computeMLResidual() 18 0.0106 0.01589 0.01986 0.14% MLMG::mgVcycle_up::5 128 0.01085 0.01153 0.01937 0.14% MLMG::ResNormInf() 146 0.01253 0.01532 0.01671 0.12% FabArray::norminf() 282 0.01096 0.0113 0.01562 0.11% Castro::problem_post_timestep() 2 0.009381 0.0125 0.01533 0.11% Castro::wd_update() 2 0.009012 0.01246 0.01532 0.11% FabArrayBase::getFB() 6380 0.01003 0.01079 0.01444 0.10% amrex::Add() 128 0.01396 0.01407 0.01426 0.10% Gravity::test_composite_phi() 1 0.01412 0.01414 0.01417 0.10% FabArrayBase::getCPC() 2005 0.009346 0.01071 0.01414 0.10% MLLinOp::define() 19 0.00363 0.007476 0.01171 0.08% MLLinOp::defineGrids() 19 0.00355 0.007387 0.01162 0.08% Other 2587 0.06787 0.08596 0.131 0.92% -------------------------------------------------------------------------------------------- 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 13 B 17 B 22 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.05) finalized