Initializing AMReX (24.09-38-g1bfdb84c124a)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.09-38-g1bfdb84c124a) initialized Starting run at 07:37:21 UTC on 2024-09-23. Successfully read inputs file ... Castro git describe: 24.09-12-gfb34ee127 AMReX git describe: 24.09-38-g1bfdb84c1 Microphysics git describe: 24.09-5-g397d3d20 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.220938435 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001166962 new finest: 0 Gravity::fill_multipole_BCs() time = 0.012285411 Gravity::multilevel_solve_for_new_phi() time = 0.118570015 ... norm of composite residual at level 0 1.193955541e-11 Gravity::fill_multipole_BCs() time = 0.01216114 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.049171287 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004246094 seconds [Level 0 step 1] ADVANCE at time 0 with dt = 0.25 [1m[32m Beginning subcycle 1 starting at time 0 with dt = 0.25[0m Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. ... 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 = Invalid density = -8.076810125 at index 20, 174, 0 -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0Invalid density = -108.2035501 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, Invalid density = -415.2911996 at index 16, 72, 0 178, 0 Invalid density = -80.16477042 at index 18, 179, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.7694089 at index 17, 181, Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477042 at index 18, 76, 0 Invalid density = -415.2911996 at index 16, 0 Invalid density = -61.58489855 at index 19, 77, 0 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -108.2035501 at index 19, 78, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 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 Castro::construct_old_gravity() time = 0.004677489 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.065958436 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry. [1m[32m Subcycle completed[0m [1m[31m Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125 [0m [1m[32m Beginning subcycle 1 starting at time 0 with dt = 0.125[0m 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.007517227 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.06499144 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry. [1m[32m Subcycle completed[0m [1m[31m Timestep 0.125 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.0625 [0m [1m[32m Beginning subcycle 1 starting at time 0 with dt = 0.0625[0m 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.007444307 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.06112014 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry. [1m[32m Subcycle completed[0m [1m[31m Timestep 0.0625 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.03125 [0m [1m[32m Beginning subcycle 1 starting at time 0 with dt = 0.03125[0m 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.005984193 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.060518346 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.32295802 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.015733093 Gravity::solve_for_phi() time = 0.100595618 on level 0 Castro::construct_new_gravity() time = 0.102818824 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.061380148 on level 0 Castro::construct_old_gravity() time = 0.006183582 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.06062893 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.328822168 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.015302897 Gravity::solve_for_phi() time = 0.099605739 on level 0 Castro::construct_new_gravity() time = 0.101836886 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.059900023 on level 0 Castro::construct_old_gravity() time = 0.004282405 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.05945259 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.332826896 on level 0 Ending SDC iteration 3 of 3. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 2 starting at time 0.03125 with dt = 0.03125[0m 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.01460243 Gravity::solve_for_phi() time = 0.098953764 on level 0 Castro::construct_new_gravity() time = 0.101499062 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.061025748 on level 0 Castro::construct_old_gravity() time = 0.006389097 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.060214172 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.985598081e+46 ener added: 2.745852676e+42 Castro::construct_ctu_hydro_source() time = 0.335438679 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.014649033 Gravity::solve_for_phi() time = 0.096675227 on level 0 Castro::construct_new_gravity() time = 0.098899109 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.061400474 on level 0 Castro::construct_old_gravity() time = 0.005536675 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.060723444 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.986297782e+46 ener added: 2.752525031e+42 Castro::construct_ctu_hydro_source() time = 0.337840999 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 3 starting at time 0.0625 with dt = 0.03125[0m 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.014603133 Gravity::solve_for_phi() time = 0.091497299 on level 0 Castro::construct_new_gravity() time = 0.093674558 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.058712432 on level 0 Castro::construct_old_gravity() time = 0.008708285 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.060601368 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.98733239e+46 ener added: 2.180535114e+43 Castro::construct_ctu_hydro_source() time = 0.335217943 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.014396886 Gravity::solve_for_phi() time = 0.096598649 on level 0 Castro::construct_new_gravity() time = 0.098784188 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.059696012 on level 0 Castro::construct_old_gravity() time = 0.002955251 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.060060055 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.988039066e+46 ener added: 2.185480971e+43 Castro::construct_ctu_hydro_source() time = 0.331136786 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 4 starting at time 0.09375 with dt = 0.03125[0m 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.015646457 Gravity::solve_for_phi() time = 0.095423609 on level 0 Castro::construct_new_gravity() time = 0.097633828 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.060631884 on level 0 Castro::construct_old_gravity() time = 0.002855362 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.060106875 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.98498831e+46 ener added: 2.292089986e+43 Castro::construct_ctu_hydro_source() time = 0.326990051 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.016205901 Gravity::solve_for_phi() time = 0.09573349 on level 0 Castro::construct_new_gravity() time = 0.098072988 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.060571859 on level 0 Castro::construct_old_gravity() time = 0.003912002 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.060790812 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.985698801e+46 ener added: 2.298024982e+43 Castro::construct_ctu_hydro_source() time = 0.36071689 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 5 starting at time 0.125 with dt = 0.03125[0m 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.020753727 Gravity::solve_for_phi() time = 0.109174188 on level 0 Castro::construct_new_gravity() time = 0.111384704 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.061016097 on level 0 Castro::construct_old_gravity() time = 0.002856339 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.060255841 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.992363162e+46 ener added: 2.321339475e+43 Castro::construct_ctu_hydro_source() time = 0.336539971 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.014975975 Gravity::solve_for_phi() time = 0.095214661 on level 0 Castro::construct_new_gravity() time = 0.097346138 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.062295295 on level 0 Castro::construct_old_gravity() time = 0.009614739 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.065782303 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.993083003e+46 ener added: 2.326773931e+43 Castro::construct_ctu_hydro_source() time = 0.348221102 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 6 starting at time 0.15625 with dt = 0.03125[0m 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.016540832 Gravity::solve_for_phi() time = 0.097642701 on level 0 Castro::construct_new_gravity() time = 0.099830143 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.060397596 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.920709246e+29 ener added: 4.920709246e+29 Castro::construct_old_gravity() time = 0.002842356 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.060062932 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.000609588e+46 ener added: 2.406064512e+43 Castro::construct_ctu_hydro_source() time = 0.326221813 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.016387567 Gravity::solve_for_phi() time = 0.0960451 on level 0 Castro::construct_new_gravity() time = 0.098264145 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.060466698 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.920709246e+29 ener added: 4.920709246e+29 Castro::construct_old_gravity() time = 0.002867139 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.060105854 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.001333674e+46 ener added: 2.41251015e+43 Castro::construct_ctu_hydro_source() time = 0.331516951 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 7 starting at time 0.1875 with dt = 0.03125[0m 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.016613603 Gravity::solve_for_phi() time = 0.096877141 on level 0 Castro::construct_new_gravity() time = 0.09907558 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.060514806 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.075345155e+29 ener added: 2.075345155e+29 Castro::construct_old_gravity() time = 0.002990378 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.065750536 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.993854135e+46 ener added: 2.676087223e+43 Castro::construct_ctu_hydro_source() time = 0.348548708 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.01643023 Gravity::solve_for_phi() time = 0.100608075 on level 0 Castro::construct_new_gravity() time = 0.102765912 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.064695857 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.075345155e+29 ener added: 2.075345155e+29 Castro::construct_old_gravity() time = 0.003944382 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.060033878 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.994586452e+46 ener added: 2.682600444e+43 Castro::construct_ctu_hydro_source() time = 0.330942602 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 8 starting at time 0.21875 with dt = 0.03125[0m 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.016912575 Gravity::solve_for_phi() time = 0.096318393 on level 0 Castro::construct_new_gravity() time = 0.098451749 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.060498647 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.383563436e+30 ener added: 1.383563436e+30 Castro::construct_old_gravity() time = 0.002781636 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.060009364 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.993555556e+46 ener added: 2.937230463e+43 Castro::construct_ctu_hydro_source() time = 0.327696421 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.017103781 Gravity::solve_for_phi() time = 0.096699038 on level 0 Castro::construct_new_gravity() time = 0.098816619 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.060570094 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.383563436e+30 ener added: 1.383563436e+30 Castro::construct_old_gravity() time = 0.002854895 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.060088963 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.994288712e+46 ener added: 2.943682879e+43 Castro::construct_ctu_hydro_source() time = 0.331831848 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m Subcycling complete Zones advanced per microsecond at this level: 0.002696380757 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.017037761 Gravity::solve_for_phi() time = 0.096757941 on level 0 Castro::construct_new_gravity() time = 0.098891893 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.06048691 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.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 TIME = 0.25 DT = 0.25 [STEP 1] Coarse TimeStep time: 12.21311508 [STEP 1] FAB kilobyte spread across MPI nodes: [6973 ... 6973] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.049336596 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004510681 seconds Ending run at 07:37:35 UTC on 2024-09-23. Run time = 14.58166085 Run time without initialization = 12.26796718 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.5816 ... 14.5816 ... 14.5817 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.2730 5.3678 5.4348 37.27% Amr::InitAmr() 1 1.347 1.347 1.347 9.24% Castro::computeTemp() 80 0.8468 0.9597 1.132 7.76% FabArray::ParallelCopy_finish() 1429 0.2607 0.6577 0.9319 6.39% amrex::Copy() 724 0.7906 0.8034 0.8127 5.57% Castro::normalize_species() 79 0.5469 0.5597 0.5888 4.04% MultiFab::AddProduct() 3439 0.5435 0.5513 0.5587 3.83% Castro::reset_internal_energy(Fab) 160 0.5203 0.5295 0.5415 3.71% MLPoisson::Fsmooth() 6144 0.3507 0.3841 0.4657 3.19% check_for_negative_density() 20 0.246 0.3147 0.4177 2.86% Castro::initMFs() 1 0.3614 0.3702 0.4172 2.86% FabArray::setVal() 5013 0.35 0.359 0.3705 2.54% MultiFab::Subtract() 159 0.3506 0.3563 0.3607 2.47% Gravity::fill_multipole_BCs() 19 0.207 0.2648 0.2948 2.02% FabArray::ParallelCopy_nowait() 1429 0.2066 0.2232 0.2514 1.72% FillBoundary_finish() 6380 0.1431 0.1699 0.2062 1.41% TagBoxArray::collate() 2 0.0002582 0.05129 0.1908 1.31% FabArray::sum() 3439 0.1811 0.1838 0.1866 1.28% Castro::initData() 1 0.0033 0.02806 0.1634 1.12% FillBoundary_nowait() 6380 0.09282 0.1232 0.1588 1.09% MLCellLinOp::applyBC() 7148 0.06408 0.08358 0.1315 0.90% MLPoisson::Fapply() 986 0.09377 0.09899 0.1145 0.79% FabArray::Saxpy() 264 0.06861 0.07002 0.07783 0.53% Castro::construct_new_gravity_source() 17 0.06846 0.07035 0.07265 0.50% FabArray::Xpay() 945 0.05214 0.05505 0.06246 0.43% MultiFab::contains_nan() 37 0.04563 0.04639 0.04774 0.33% Castro::construct_old_gravity_source() 20 0.04458 0.04547 0.04738 0.32% MLMG::addInterpCorrection() 768 0.03398 0.03695 0.0453 0.31% AmrLevel::derive() 18 0.0408 0.04141 0.04215 0.29% amrex::average_down 768 0.02529 0.02806 0.03456 0.24% StateData::FillBoundary(geom) 143 0 0.01779 0.03439 0.24% Castro::reset_internal_energy(MultiFab) 80 0.02356 0.02416 0.02489 0.17% VisMF::Write(FabArray) 4 0.01383 0.01475 0.02468 0.17% MLCellLinOp::defineAuxData() 19 0.01612 0.01708 0.02052 0.14% Castro::enforce_min_density() 79 0.0185 0.0188 0.01927 0.13% Castro::construct_new_source() 102 0.01591 0.01637 0.01823 0.12% MLCGSolver::bicgstab 8 0 0.001126 0.01801 0.12% Amr::writePlotFile() 2 0.007043 0.01587 0.01748 0.12% FabArray::FillBoundary() 6380 0.01384 0.01448 0.01747 0.12% FabArray::norminf() 282 0.01151 0.01203 0.01624 0.11% amrex::Add() 128 0.01469 0.0149 0.0151 0.10% Castro::evaluate_source_change() 196 0.01383 0.01415 0.01448 0.10% FabArray::setDomainBndry() 147 0.0004755 0.007582 0.01441 0.10% Castro::wd_update() 2 0.005218 0.009957 0.01439 0.10% FillPatchIterator::Initialize 147 0.01288 0.01329 0.01423 0.10% FabArrayBase::getFB() 6380 0.01081 0.01129 0.01381 0.09% MLCellLinOp::smooth() 3072 0.01184 0.01243 0.0134 0.09% MLMG::ResNormInf() 146 0.005604 0.008795 0.01256 0.09% MLMG::prepareForSolve() 18 0.01037 0.01085 0.01186 0.08% Castro::construct_old_source() 120 0.01019 0.01036 0.01046 0.07% MLMG::solve() 18 0.001577 0.005038 0.009724 0.07% Gravity::actual_solve_with_mlmg() 19 0.007562 0.007912 0.009345 0.06% FabArrayBase::getCPC() 2005 0.007601 0.00776 0.008132 0.06% Gravity::get_old_grav_vector() 20 0.006475 0.006693 0.007296 0.05% Gravity::get_new_grav_vector() 18 0.006055 0.006165 0.006385 0.04% Castro::initialize_do_advance() 20 0.004332 0.004818 0.005704 0.04% MLCellLinOp::apply() 986 0.002051 0.002285 0.004804 0.03% FabArray::ParallelCopy() 1429 0.003538 0.003695 0.00402 0.03% Castro::sum_integrated_quantities() 2 0.001951 0.00233 0.002625 0.02% AmrLevel::FillPatch() 147 0.002389 0.002503 0.002596 0.02% MLCellLinOp::correctionResidual() 768 0.00188 0.001948 0.002064 0.01% MLMG::mgVcycle() 128 0.001594 0.001674 0.001834 0.01% MLMG:computeResOfCorrection() 768 0.00147 0.00152 0.00161 0.01% Castro::save_data_for_retry() 3 0.001337 0.001412 0.001574 0.01% MLMG::actualBottomSolve() 128 0.0001339 0.0002003 0.001114 0.01% Castro::subcycle_advance_ctu() 1 0.0002554 0.0003436 0.00102 0.01% StateDataPhysBCFunct::() 147 0.0007151 0.0008681 0.000992 0.01% main() 1 0.000715 0.0007709 0.0008997 0.01% Gravity::solve_for_phi() 17 0.0006682 0.0007037 0.0008227 0.01% MLMG::mgVcycle_down::3 128 0.0005975 0.000649 0.0007606 0.01% MLMG::mgVcycle_down::2 128 0.000643 0.000696 0.0007435 0.01% FillPatchIterator::FillFromLevel0() 147 0.0006844 0.0007037 0.000727 0.00% MLMG::mgVcycle_down::0 128 0.0006809 0.0006985 0.0007208 0.00% MLMG::mgVcycle_down::4 128 0.0005777 0.0006068 0.0007105 0.00% MLMG::mgVcycle_down::1 128 0.0006756 0.0006911 0.0007057 0.00% MLMG::mgVcycle_down::5 128 0.0005981 0.0006249 0.0007003 0.00% AmrMesh::MakeNewGrids() 2 0.000548 0.0005731 0.0006337 0.00% Castro::Castro() 1 0.000468 0.0005181 0.000576 0.00% MLMG::mgVcycle_up::1 128 0.000483 0.0005046 0.0005382 0.00% MLMG::mgVcycle_up::0 128 0.0004881 0.000502 0.00053 0.00% Castro::clean_state() 79 0.0004485 0.0004763 0.0005101 0.00% MLMG::mgVcycle_up::2 128 0.0004169 0.0004554 0.0005045 0.00% FillPatchSingleLevel 147 0.0004697 0.0004807 0.0004984 0.00% Castro::do_advance_ctu() 20 0.0003973 0.0004258 0.0004978 0.00% MLMG::mgVcycle_up::3 128 0.0004133 0.0004299 0.0004935 0.00% MLMG::mgVcycle_up::4 128 0.0004194 0.0004299 0.0004704 0.00% MLMG::mgVcycle_up::5 128 0.0004175 0.0004281 0.0004695 0.00% MLCellLinOp::solutionResidual() 147 0.000446 0.0004555 0.0004655 0.00% MLMG::oneIter() 128 0.0003971 0.0004064 0.000414 0.00% Castro::problem_post_timestep() 2 1.079e-05 3.404e-05 0.0003493 0.00% Castro::do_old_sources() 20 0.0002895 0.0002987 0.0003207 0.00% Castro::construct_new_gravity() 17 0.000154 0.0001745 0.0003123 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002196 0.0002328 0.000263 0.00% Castro::do_new_sources() 17 0.0002184 0.0002295 0.0002589 0.00% MLMG::computeResidual() 128 0.0002347 0.0002423 0.0002534 0.00% MLMG::mgVcycle_bottom 128 0.0002235 0.0002328 0.0002533 0.00% Gravity::test_composite_phi() 1 0.0001065 0.0001164 0.0001792 0.00% Amr::coarseTimeStep() 1 0.0001492 0.0001737 0.0001789 0.00% Castro::retry_advance_ctu() 11 6.13e-05 8.085e-05 0.00017 0.00% TagBoxArray::mapPRD 2 0.0001294 0.0001367 0.0001462 0.00% MLPoisson::define() 19 0.0001203 0.0001273 0.000137 0.00% Castro::check_for_nan() 37 0.0001077 0.0001152 0.0001327 0.00% Castro::construct_old_gravity() 20 0.0001109 0.0001186 0.0001302 0.00% Castro::initialize_advance() 1 0.0001046 0.0001114 0.0001157 0.00% Castro::apply_source_to_state() 37 9.826e-05 0.0001041 0.0001148 0.00% Amr::defBaseLevel() 1 9.772e-05 0.0001065 0.0001142 0.00% Castro::create_source_corrector() 17 5.576e-05 6.354e-05 0.0001124 0.00% MLLinOp::define() 19 7.269e-05 7.767e-05 8.322e-05 0.00% Amr::grid_places() 2 5.474e-05 6.003e-05 6.894e-05 0.00% Amr::FinalizeInit() 1 1.402e-05 2.696e-05 6.597e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.855e-05 6.144e-05 6.493e-05 0.00% Castro::expand_state() 20 5.062e-05 5.365e-05 6.489e-05 0.00% MLMG::computeMLResidual() 18 4.086e-05 4.411e-05 5.067e-05 0.00% Amr::regrid() 1 1.609e-05 1.985e-05 4.538e-05 0.00% Castro::derive() 18 3.989e-05 4.164e-05 4.407e-05 0.00% Gravity::actual_multilevel_solve() 1 3.135e-05 3.325e-05 3.744e-05 0.00% Castro::advance() 1 6.751e-06 8.841e-06 2.344e-05 0.00% Amr::timeStep() 1 1.391e-05 1.692e-05 2.335e-05 0.00% Castro::post_regrid() 2 1.201e-05 1.423e-05 2.113e-05 0.00% Amr::InitializeInit() 1 8.582e-06 9.877e-06 1.479e-05 0.00% Castro::post_init() 1 1.194e-05 1.269e-05 1.407e-05 0.00% Amr::bldFineLevels() 1 7.427e-06 9.048e-06 1.053e-05 0.00% Amr::init() 1 6.112e-06 7.289e-06 7.913e-06 0.00% Castro::post_timestep() 1 5.816e-06 6.253e-06 7.625e-06 0.00% Amr::initialInit() 1 2.985e-06 3.183e-06 3.768e-06 0.00% Other 2459 0.05022 0.06221 0.09076 0.62% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.58 14.58 14.58 100.00% Amr::coarseTimeStep() 1 12.27 12.27 12.27 84.13% Amr::timeStep() 1 12.21 12.21 12.21 83.76% Castro::advance() 1 12.15 12.15 12.15 83.35% Castro::subcycle_advance_ctu() 1 12.1 12.1 12.1 83.00% Castro::do_advance_ctu() 20 12.06 12.06 12.06 82.72% Castro::construct_ctu_hydro_source() 20 6.438 6.468 6.533 44.80% Castro::clean_state() 79 3.671 3.815 4.006 27.47% Castro::computeTemp() 80 2.355 2.479 2.666 18.28% Gravity::solve_phi_with_mlmg() 19 1.696 1.758 1.786 12.25% Castro::construct_new_gravity() 17 1.603 1.665 1.693 11.61% Gravity::solve_for_phi() 17 1.565 1.627 1.656 11.35% Castro::reset_internal_energy(MultiFab) 80 1.496 1.52 1.542 10.57% Castro::initialize_do_advance() 20 1.418 1.468 1.525 10.46% Gravity::actual_solve_with_mlmg() 19 1.486 1.49 1.494 10.24% MLMG::solve() 18 1.429 1.435 1.437 9.85% Amr::InitAmr() 1 1.347 1.347 1.347 9.24% MLMG::oneIter() 128 1.285 1.288 1.292 8.86% MLMG::mgVcycle() 128 1.269 1.273 1.276 8.75% Castro::do_old_sources() 20 1.175 1.196 1.216 8.34% FabArray::ParallelCopy() 1429 0.5108 0.9051 1.175 8.06% Castro::do_new_sources() 17 0.9184 0.9723 1.023 7.02% Amr::init() 1 0.9667 0.9667 0.9667 6.63% Castro::evaluate_source_change() 196 0.9375 0.9515 0.9645 6.61% FabArray::ParallelCopy_finish() 1429 0.2607 0.6577 0.9319 6.39% Amr::initialInit() 1 0.9129 0.9129 0.9129 6.26% MLCellLinOp::smooth() 3072 0.6634 0.7224 0.8739 5.99% Castro::enforce_min_density() 79 0.7914 0.8031 0.8153 5.59% amrex::Copy() 724 0.7906 0.8034 0.8127 5.57% AmrLevel::FillPatch() 147 0.5588 0.6697 0.7746 5.31% Amr::InitializeInit() 1 0.4472 0.4827 0.6663 4.57% Amr::defBaseLevel() 1 0.4471 0.4827 0.6663 4.57% FillPatchIterator::Initialize 147 0.3755 0.4874 0.5927 4.06% Castro::normalize_species() 79 0.5469 0.5597 0.5888 4.04% FillPatchIterator::FillFromLevel0() 147 0.3477 0.4666 0.573 3.93% FillPatchSingleLevel 147 0.347 0.4659 0.5722 3.92% MultiFab::AddProduct() 3439 0.5435 0.5513 0.5587 3.83% Castro::reset_internal_energy(Fab) 160 0.5203 0.5295 0.5415 3.71% MLMG::mgVcycle_up::1 128 0.1155 0.3062 0.4995 3.43% MLCellLinOp::applyBC() 7148 0.3616 0.4032 0.4924 3.38% Amr::FinalizeInit() 1 0.2466 0.4301 0.4657 3.19% MLPoisson::Fsmooth() 6144 0.3507 0.3841 0.4657 3.19% MLMG::addInterpCorrection() 768 0.0687 0.3471 0.4298 2.95% Castro::Castro() 1 0.3673 0.3764 0.4246 2.91% check_for_negative_density() 20 0.246 0.3147 0.4177 2.86% Castro::initMFs() 1 0.3616 0.3704 0.4173 2.86% FabArray::setVal() 5013 0.35 0.359 0.3705 2.54% FabArray::FillBoundary() 6380 0.2842 0.3196 0.3609 2.47% MultiFab::Subtract() 159 0.3506 0.3563 0.3607 2.47% MLMG::mgVcycle_up::2 128 0.01395 0.1366 0.3245 2.23% Gravity::fill_multipole_BCs() 19 0.2092 0.267 0.297 2.04% MLMG::mgVcycle_down::0 128 0.2858 0.2891 0.2935 2.01% Castro::expand_state() 20 0.2044 0.2416 0.2904 1.99% FabArray::ParallelCopy_nowait() 1429 0.2266 0.2437 0.2728 1.87% Amr::bldFineLevels() 1 0.04725 0.2304 0.265 1.82% Castro::initData() 1 0.07807 0.1046 0.2401 1.65% MLMG::mgVcycle_up::0 128 0.221 0.2232 0.2258 1.55% Amr::grid_places() 2 0.002968 0.1865 0.2221 1.52% AmrMesh::MakeNewGrids() 2 0.002912 0.1864 0.222 1.52% TagBoxArray::mapPRD 2 0.001314 0.1338 0.2203 1.51% MLCellLinOp::apply() 986 0.1649 0.1764 0.2149 1.47% FillBoundary_finish() 6380 0.1431 0.1699 0.2062 1.41% MLMG::mgVcycle_up::3 128 0.01387 0.02872 0.2044 1.40% TagBoxArray::collate() 2 0.0002582 0.05129 0.1908 1.31% FabArray::sum() 3439 0.1811 0.1838 0.1866 1.28% FillBoundary_nowait() 6380 0.1041 0.1352 0.1717 1.18% MLMG:computeResOfCorrection() 768 0.1225 0.1332 0.1587 1.09% Castro::post_init() 1 0.1539 0.1563 0.1581 1.08% MLCellLinOp::correctionResidual() 768 0.1211 0.1317 0.1571 1.08% MLMG::mgVcycle_down::1 128 0.1402 0.1447 0.1517 1.04% Castro::construct_new_source() 102 0.1233 0.1261 0.1304 0.89% Gravity::multilevel_solve_for_new_phi() 1 0.1185 0.1185 0.1186 0.81% Gravity::actual_multilevel_solve() 1 0.1184 0.1184 0.1185 0.81% MLPoisson::Fapply() 986 0.09377 0.09899 0.1145 0.79% Castro::create_source_corrector() 17 0.06138 0.08789 0.112 0.77% MLCellLinOp::solutionResidual() 147 0.09359 0.1003 0.1064 0.73% Amr::writePlotFile() 2 0.09886 0.09894 0.09895 0.68% Castro::post_regrid() 2 0.08544 0.08729 0.08966 0.61% MLMG::mgVcycle_down::2 128 0.01629 0.04962 0.08913 0.61% MLMG::computeResidual() 128 0.07815 0.08096 0.08367 0.57% Castro::construct_old_source() 120 0.07797 0.07935 0.08318 0.57% MLMG::mgVcycle_bottom 128 0.0003584 0.005258 0.07857 0.54% MLMG::actualBottomSolve() 128 0.0001339 0.005026 0.07832 0.54% FabArray::Saxpy() 264 0.06861 0.07002 0.07783 0.53% Castro::construct_old_gravity() 20 0.05349 0.06736 0.07677 0.53% MLCGSolver::bicgstab 8 0 0.004795 0.07672 0.53% Gravity::get_old_grav_vector() 20 0.05337 0.06725 0.07664 0.53% MLMG::mgVcycle_down::3 128 0.01552 0.02302 0.07396 0.51% Castro::construct_new_gravity_source() 17 0.06846 0.07035 0.07265 0.50% Castro::derive() 18 0.0701 0.07106 0.07213 0.49% AmrLevel::derive() 18 0.07006 0.07101 0.07209 0.49% Castro::post_timestep() 1 0.05955 0.06306 0.06656 0.46% amrex::average_down 768 0.0343 0.04283 0.06391 0.44% FabArray::Xpay() 945 0.05214 0.05505 0.06246 0.43% Castro::initialize_advance() 1 0.04915 0.05026 0.05227 0.36% Castro::check_for_nan() 37 0.04574 0.04651 0.04786 0.33% MultiFab::contains_nan() 37 0.04563 0.04639 0.04774 0.33% Castro::construct_old_gravity_source() 20 0.04458 0.04547 0.04738 0.32% Amr::regrid() 1 0.044 0.04502 0.04689 0.32% MLPoisson::define() 19 0.03138 0.0366 0.04285 0.29% MLMG::mgVcycle_down::4 128 0.01548 0.01746 0.04158 0.29% Castro::retry_advance_ctu() 11 0.03588 0.03758 0.03961 0.27% Gravity::get_new_grav_vector() 18 0.03688 0.0375 0.03831 0.26% StateDataPhysBCFunct::() 147 0.0007151 0.01866 0.03531 0.24% StateData::FillBoundary(geom) 143 0 0.01779 0.03439 0.24% MLMG::mgVcycle_up::4 128 0.01398 0.01541 0.03298 0.23% MLMG::mgVcycle_down::5 128 0.01543 0.01661 0.02901 0.20% Castro::apply_source_to_state() 37 0.02799 0.02838 0.02894 0.20% MLCellLinOp::defineAuxData() 19 0.02039 0.02213 0.02775 0.19% MLMG::mgVcycle_up::5 128 0.0141 0.01503 0.02503 0.17% MLMG::prepareForSolve() 18 0.02083 0.02209 0.02483 0.17% VisMF::Write(FabArray) 4 0.01383 0.01475 0.02468 0.17% MLMG::computeMLResidual() 18 0.01155 0.01884 0.02397 0.16% MLMG::ResNormInf() 146 0.01585 0.01902 0.02281 0.16% Castro::sum_integrated_quantities() 2 0.02189 0.02215 0.02217 0.15% Castro::save_data_for_retry() 3 0.01893 0.01968 0.02156 0.15% Castro::problem_post_timestep() 2 0.009746 0.01367 0.01782 0.12% Castro::wd_update() 2 0.009732 0.01364 0.0178 0.12% FabArray::norminf() 282 0.01151 0.01203 0.01624 0.11% Gravity::test_composite_phi() 1 0.01597 0.01602 0.01606 0.11% FabArrayBase::getFB() 6380 0.01119 0.01205 0.0155 0.11% FabArrayBase::getCPC() 2005 0.01006 0.01139 0.01536 0.11% amrex::Add() 128 0.01469 0.0149 0.0151 0.10% FabArray::setDomainBndry() 147 0.0004755 0.007582 0.01441 0.10% MLLinOp::define() 19 0.004185 0.007812 0.01145 0.08% Other 2459 0.08425 0.1018 0.1401 0.96% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Pinned Memory Usage: --------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max --------------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 16 11 B 13 B 16 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 1560 KiB 1561 KiB 1562 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 381 KiB 388 KiB 393 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1041 KiB 1042 KiB 1042 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 991 KiB 991 KiB 992 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 661 KiB 662 KiB 662 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 4272 B 4362 B 4868 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 50 KiB 50 KiB 51 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 12 KiB 14 KiB 16 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 414 KiB 414 KiB 414 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 2263 B 4391 B 6165 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1332 B 1349 B 1372 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 17 KiB 17 KiB 17 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 131 KiB 131 KiB 131 KiB 157 KiB 157 KiB 157 KiB MLMG::prepareForSolve() 6102 6102 10 KiB 10 KiB 12 KiB 102 KiB 108 KiB 126 KiB Castro::buildMetrics() 160 160 106 KiB 106 KiB 106 KiB 117 KiB 117 KiB 117 KiB Castro::construct_old_source() 640 640 254 B 258 B 274 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 425 B 431 B 447 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 18 B 20 B 37 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 33 KiB 33 KiB 33 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 57 B 57 B 58 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 3 B 5 B 8 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 114 B 115 B 118 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 167 B 210 B 241 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 64 31 KiB 31 KiB 31 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 2432 235 B 243 B 277 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1048 B 1063 B 1078 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 132 B 132 B 133 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1757 B 1826 B 1859 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 465 B 519 B 720 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 483 B 485 B 486 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 29 B 40 B 54 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 23 B 39 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 17 B 27 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 7 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.09-38-g1bfdb84c124a) finalized