Initializing AMReX (24.06-21-g259db7cfb99e)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.06-21-g259db7cfb99e) initialized Starting run at 07:24:14 UTC on 2024-06-24. Successfully read inputs file ... Castro git describe: 24.06-16-g830d937a6 AMReX git describe: 24.06-21-g259db7cfb Microphysics git describe: 24.06-31-g48149cdd 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.173082315 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001083775 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010561122 Gravity::multilevel_solve_for_new_phi() time = 0.10053377 Gravity::fill_multipole_BCs() time = 0.010511938 ... 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.04846042 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.003931351 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.002694673 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.065357239 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.007544675 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.064561784 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.004857759 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.060963824 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.006925994 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.059393515 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.322625249 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015799305 Gravity::solve_for_phi() time = 0.086466443 on level 0 Castro::construct_new_gravity() time = 0.088545087 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.061290042 on level 0 Ending SDC iteration 1 of 3. Beginning SDC iteration 2 of 3. Castro::construct_old_gravity() time = 0.012577367 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.068134307 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.335806869 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01548512 Gravity::solve_for_phi() time = 0.086405311 on level 0 Castro::construct_new_gravity() time = 0.08844726 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.060908389 on level 0 Ending SDC iteration 2 of 3. Beginning SDC iteration 3 of 3. Castro::construct_old_gravity() time = 0.003050559 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.05899914 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.327058837 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015426259 Gravity::solve_for_phi() time = 0.086540745 on level 0 Castro::construct_new_gravity() time = 0.088592256 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.059851715 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.003365241 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.058805473 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.329826032 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015934522 Gravity::solve_for_phi() time = 0.087102288 on level 0 Castro::construct_new_gravity() time = 0.089144311 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.059854814 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.002705871 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.060196634 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.340261603 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016059674 Gravity::solve_for_phi() time = 0.086555935 on level 0 Castro::construct_new_gravity() time = 0.088622726 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.061005927 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.002661955 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.058676083 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.329303468 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015981859 Gravity::solve_for_phi() time = 0.087557357 on level 0 Castro::construct_new_gravity() time = 0.089601112 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.06037956 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.003001476 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.058826133 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.331486864 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014876855 Gravity::solve_for_phi() time = 0.084421637 on level 0 Castro::construct_new_gravity() time = 0.086462168 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.059816983 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.002588929 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.058953105 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.323194671 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015345641 Gravity::solve_for_phi() time = 0.084570275 on level 0 Castro::construct_new_gravity() time = 0.086625232 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.058567851 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.003863843 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.059171819 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.328301092 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015021987 Gravity::solve_for_phi() time = 0.084574875 on level 0 Castro::construct_new_gravity() time = 0.086617291 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.05887558 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.002540674 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.058440935 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.323190621 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015786443 Gravity::solve_for_phi() time = 0.084833975 on level 0 Castro::construct_new_gravity() time = 0.0868825 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.059374734 on level 0 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. Castro::construct_old_gravity() time = 0.004649976 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.058680992 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.327498646 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015642453 Gravity::solve_for_phi() time = 0.08512708 on level 0 Castro::construct_new_gravity() time = 0.087172761 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.05958761 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.004069661 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.059924587 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.32408094 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015294253 Gravity::solve_for_phi() time = 0.084281414 on level 0 Castro::construct_new_gravity() time = 0.086315142 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.063091088 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.004995666 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.059825609 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.329186155 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01528698 Gravity::solve_for_phi() time = 0.084682569 on level 0 Castro::construct_new_gravity() time = 0.086715616 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.060491453 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.003002808 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.059229321 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.323301385 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015469778 Gravity::solve_for_phi() time = 0.086384261 on level 0 Castro::construct_new_gravity() time = 0.088455505 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.061855179 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.002822367 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.058201373 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.33200564 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015932277 Gravity::solve_for_phi() time = 0.085776885 on level 0 Castro::construct_new_gravity() time = 0.087876813 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.060211836 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.004386104 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.057808598 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.322698478 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01513678 Gravity::solve_for_phi() time = 0.084563491 on level 0 Castro::construct_new_gravity() time = 0.086647384 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.058849177 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.002534007 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.059415736 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.340313115 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015987281 Gravity::solve_for_phi() time = 0.08520092 on level 0 Castro::construct_new_gravity() time = 0.08723979 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.059173041 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002782261805 [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.83749703 [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.050816752 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.003913486 seconds Ending run at 07:24:28 UTC on 2024-06-24. Run time = 14.24610726 Run time without initialization = 11.89309188 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.25 ... 14.25 ... 14.25 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.3339 5.4360 5.5218 38.76% Amr::InitAmr() 1 1.451 1.451 1.451 10.19% Castro::computeTemp() 80 0.8505 0.9602 1.112 7.81% amrex::Copy() 724 0.7558 0.7741 0.8897 6.25% FabArray::ParallelCopy_finish() 1429 0.06184 0.543 0.8244 5.79% Castro::normalize_species() 79 0.541 0.5541 0.5883 4.13% MultiFab::AddProduct() 3439 0.5446 0.551 0.56 3.93% Castro::reset_internal_energy(Fab) 160 0.5183 0.5255 0.5323 3.74% MLPoisson::Fsmooth() 6144 0.3119 0.3392 0.4128 2.90% Castro::initMFs() 1 0.359 0.3668 0.3892 2.73% MultiFab::Subtract() 159 0.346 0.3508 0.3554 2.49% FabArray::setVal() 5013 0.3355 0.3426 0.3528 2.48% Gravity::fill_multipole_BCs() 19 0.1949 0.2357 0.2831 1.99% Castro::reset_internal_energy(MultiFab) 80 0.02209 0.08245 0.2553 1.79% check_for_negative_density() 20 0.05147 0.1572 0.2528 1.77% FabArray::ParallelCopy_nowait() 1429 0.1803 0.1984 0.2182 1.53% FabArray::sum() 3439 0.1825 0.1848 0.1877 1.32% Castro::initData() 1 0.003312 0.02811 0.1665 1.17% FillBoundary_finish() 6380 0.09108 0.1108 0.1352 0.95% FillBoundary_nowait() 6380 0.07638 0.1029 0.1321 0.93% MLPoisson::Fapply() 986 0.09569 0.1005 0.1162 0.82% MLCellLinOp::applyBC() 7148 0.04831 0.06363 0.107 0.75% FabArray::Saxpy() 264 0.07176 0.07304 0.08054 0.57% Castro::enforce_min_density() 79 0.01689 0.03372 0.07959 0.56% Castro::do_new_sources() 17 0.000374 0.01339 0.07948 0.56% Castro::construct_new_gravity_source() 17 0.06711 0.06861 0.0708 0.50% Castro::do_old_sources() 20 0.0005142 0.01189 0.06739 0.47% FabArray::Xpay() 945 0.04964 0.05211 0.05945 0.42% Castro::construct_old_gravity_source() 20 0.04541 0.04598 0.04677 0.33% MultiFab::contains_nan() 37 0.04392 0.04447 0.04529 0.32% MLMG::addInterpCorrection() 768 0.03305 0.03567 0.04281 0.30% AmrLevel::derive() 18 0.04086 0.04118 0.04159 0.29% amrex::average_down 768 0.02491 0.02707 0.03301 0.23% StateData::FillBoundary(geom) 143 0 0.01688 0.03268 0.23% VisMF::Write(FabArray) 4 0.01204 0.01343 0.02718 0.19% Amr::writePlotFile() 2 0.007214 0.01928 0.02094 0.15% MLCellLinOp::defineAuxData() 19 0.01471 0.01552 0.01762 0.12% FabArray::FillBoundary() 6380 0.01334 0.014 0.01681 0.12% Castro::construct_new_source() 102 0.01556 0.01583 0.01661 0.12% MLCGSolver::bicgstab 8 0 0.001015 0.01624 0.11% FabArray::norminf() 282 0.01082 0.01124 0.01545 0.11% amrex::Add() 128 0.01418 0.01432 0.01446 0.10% Castro::construct_old_gravity() 20 0.0003677 0.002989 0.01382 0.10% Castro::evaluate_source_change() 196 0.01244 0.01281 0.01339 0.09% FabArrayBase::getFB() 6380 0.0101 0.01059 0.01326 0.09% Castro::wd_update() 2 0.005165 0.008463 0.01269 0.09% FillPatchIterator::Initialize 147 0.01206 0.01243 0.01268 0.09% MLCellLinOp::smooth() 3072 0.01079 0.01128 0.01221 0.09% FabArray::setDomainBndry() 147 0.0004751 0.005901 0.01149 0.08% MLMG::prepareForSolve() 18 0.01 0.01023 0.01081 0.08% Castro::construct_old_source() 120 0.009838 0.01003 0.01043 0.07% FabArrayBase::getCPC() 2005 0.007392 0.007605 0.007913 0.06% Gravity::actual_solve_with_mlmg() 19 0.006605 0.006912 0.007849 0.06% MLMG::ResNormInf() 146 0.001965 0.004736 0.007099 0.05% Gravity::get_old_grav_vector() 20 0.006393 0.006509 0.006656 0.05% Gravity::get_new_grav_vector() 18 0.005915 0.006118 0.006502 0.05% MLMG::solve() 18 0.000789 0.00295 0.006436 0.05% Castro::initialize_do_advance() 20 0.004278 0.004504 0.004976 0.03% MLCellLinOp::apply() 986 0.00194 0.002158 0.004552 0.03% Castro::sum_integrated_quantities() 2 0.001654 0.003711 0.004239 0.03% FabArray::ParallelCopy() 1429 0.003456 0.00359 0.003925 0.03% Castro::construct_new_gravity() 17 0.0002908 0.0008073 0.003126 0.02% AmrLevel::FillPatch() 147 0.00214 0.002242 0.002453 0.02% MLCellLinOp::correctionResidual() 768 0.00181 0.001871 0.001961 0.01% MLMG::mgVcycle() 128 0.001511 0.001589 0.001753 0.01% MLMG:computeResOfCorrection() 768 0.001352 0.001406 0.001499 0.01% Gravity::solve_for_phi() 17 0.0008079 0.000971 0.001369 0.01% Castro::save_data_for_retry() 3 0.001233 0.001268 0.0013 0.01% main() 1 0.0009752 0.001138 0.001261 0.01% MLMG::actualBottomSolve() 128 0.0001305 0.0001918 0.001056 0.01% StateDataPhysBCFunct::() 147 0.0006924 0.0008305 0.0009153 0.01% Castro::subcycle_advance_ctu() 1 0.0001855 0.000245 0.0008014 0.01% FillPatchIterator::FillFromLevel0() 147 0.0006877 0.0007144 0.00077 0.01% MLMG::mgVcycle_down::3 128 0.0005583 0.000606 0.0006975 0.00% MLMG::mgVcycle_down::1 128 0.0006413 0.0006561 0.0006916 0.00% MLMG::mgVcycle_down::2 128 0.0005971 0.0006466 0.000688 0.00% MLMG::mgVcycle_down::0 128 0.0006316 0.0006535 0.0006861 0.00% MLMG::mgVcycle_down::4 128 0.0005435 0.0005672 0.0006637 0.00% MLMG::mgVcycle_down::5 128 0.0005629 0.0005842 0.0006569 0.00% Castro::Castro() 1 0.0004839 0.0005397 0.0006419 0.00% AmrMesh::MakeNewGrids() 2 0.0005321 0.0005503 0.0005684 0.00% FillPatchSingleLevel 147 0.0004608 0.0004822 0.0005042 0.00% Castro::clean_state() 79 0.0004182 0.0004446 0.000493 0.00% MLMG::mgVcycle_up::0 128 0.0004396 0.0004561 0.0004893 0.00% MLMG::mgVcycle_up::1 128 0.0004403 0.0004604 0.0004838 0.00% MLMG::mgVcycle_up::2 128 0.0003743 0.0004157 0.0004695 0.00% Castro::do_advance_ctu() 20 0.000402 0.0004285 0.0004606 0.00% MLCellLinOp::solutionResidual() 147 0.0004294 0.0004416 0.0004582 0.00% MLMG::mgVcycle_up::3 128 0.0003687 0.0003881 0.0004508 0.00% MLMG::mgVcycle_up::5 128 0.0003881 0.0003996 0.0004443 0.00% MLMG::mgVcycle_up::4 128 0.0003758 0.0003908 0.0004275 0.00% MLMG::oneIter() 128 0.0003708 0.000382 0.0004098 0.00% Castro::problem_post_timestep() 2 1.018e-05 2.653e-05 0.0002555 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002183 0.0002293 0.0002535 0.00% MLMG::computeResidual() 128 0.0002312 0.0002381 0.000253 0.00% MLMG::mgVcycle_bottom 128 0.0002226 0.0002293 0.0002455 0.00% Castro::retry_advance_ctu() 11 6.142e-05 7.686e-05 0.0001608 0.00% MLPoisson::define() 19 0.0001211 0.0001308 0.0001521 0.00% Amr::coarseTimeStep() 1 0.0001146 0.000139 0.0001445 0.00% Castro::check_for_nan() 37 0.0001053 0.0001139 0.0001366 0.00% TagBoxArray::mapPRD 2 0.0001112 0.0001221 0.0001335 0.00% Amr::defBaseLevel() 1 9.572e-05 0.0001086 0.00013 0.00% Amr::grid_places() 2 8.173e-05 9.827e-05 0.0001252 0.00% Castro::initialize_advance() 1 0.0001012 0.00011 0.0001203 0.00% Gravity::test_composite_phi() 1 9.186e-05 9.663e-05 0.000111 0.00% Castro::apply_source_to_state() 37 9.376e-05 9.969e-05 0.0001089 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.16e-05 7.153e-05 0.0001087 0.00% MLLinOp::define() 19 7.499e-05 8.183e-05 8.938e-05 0.00% Amr::FinalizeInit() 1 8.286e-06 1.741e-05 6.989e-05 0.00% Castro::create_source_corrector() 17 5.362e-05 6.034e-05 6.623e-05 0.00% Castro::expand_state() 20 4.965e-05 5.374e-05 5.918e-05 0.00% Castro::derive() 18 4.03e-05 4.225e-05 5.223e-05 0.00% MLMG::computeMLResidual() 18 4.172e-05 4.372e-05 4.835e-05 0.00% Amr::regrid() 1 1.066e-05 1.253e-05 2.997e-05 0.00% Gravity::actual_multilevel_solve() 1 2.643e-05 2.769e-05 2.929e-05 0.00% Amr::timeStep() 1 1.203e-05 1.473e-05 1.885e-05 0.00% Castro::advance() 1 6.465e-06 9.528e-06 1.498e-05 0.00% Amr::InitializeInit() 1 6.452e-06 8.126e-06 1.381e-05 0.00% Castro::post_init() 1 1.151e-05 1.249e-05 1.322e-05 0.00% Castro::post_regrid() 2 1.062e-05 1.138e-05 1.262e-05 0.00% Amr::bldFineLevels() 1 6.845e-06 7.603e-06 9.476e-06 0.00% Amr::init() 1 5.685e-06 6.576e-06 7.692e-06 0.00% Castro::post_timestep() 1 5.644e-06 6.012e-06 6.423e-06 0.00% Amr::initialInit() 1 2.918e-06 3.267e-06 3.854e-06 0.00% Other 2461 0.04774 0.06005 0.09159 0.64% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.25 14.25 14.25 100.00% Amr::coarseTimeStep() 1 11.89 11.89 11.89 83.48% Amr::timeStep() 1 11.84 11.84 11.84 83.09% Castro::advance() 1 11.77 11.78 11.78 82.67% Castro::subcycle_advance_ctu() 1 11.72 11.73 11.73 82.33% Castro::do_advance_ctu() 20 11.68 11.69 11.69 82.08% Castro::construct_ctu_hydro_source() 20 6.34 6.373 6.409 44.99% Castro::clean_state() 79 3.647 3.81 4.005 28.11% Castro::computeTemp() 80 2.355 2.516 2.773 19.47% Castro::reset_internal_energy(MultiFab) 80 1.474 1.556 1.736 12.19% Gravity::solve_phi_with_mlmg() 19 1.474 1.518 1.566 10.99% Castro::initialize_do_advance() 20 1.397 1.441 1.492 10.47% Castro::construct_new_gravity() 17 1.403 1.443 1.49 10.46% Gravity::solve_for_phi() 17 1.364 1.407 1.455 10.21% Amr::InitAmr() 1 1.451 1.451 1.451 10.19% Gravity::actual_solve_with_mlmg() 19 1.275 1.278 1.283 9.01% MLMG::solve() 18 1.224 1.228 1.23 8.63% Castro::do_old_sources() 20 1.159 1.185 1.201 8.43% MLMG::oneIter() 128 1.096 1.098 1.101 7.73% MLMG::mgVcycle() 128 1.082 1.083 1.086 7.62% FabArray::ParallelCopy() 1429 0.2914 0.7647 1.041 7.31% Castro::do_new_sources() 17 0.9241 0.9807 1.023 7.18% Castro::evaluate_source_change() 196 0.9331 0.9442 0.9597 6.74% Amr::init() 1 0.9011 0.9012 0.9012 6.33% amrex::Copy() 724 0.7558 0.7741 0.8897 6.25% Amr::initialInit() 1 0.8479 0.848 0.8481 5.95% Castro::enforce_min_density() 79 0.7772 0.8045 0.8452 5.93% FabArray::ParallelCopy_finish() 1429 0.06184 0.543 0.8244 5.79% AmrLevel::FillPatch() 147 0.483 0.5923 0.7285 5.11% MLCellLinOp::smooth() 3072 0.545 0.5947 0.7227 5.07% Amr::InitializeInit() 1 0.4507 0.5151 0.6219 4.37% Amr::defBaseLevel() 1 0.4507 0.5151 0.6219 4.37% Castro::normalize_species() 79 0.541 0.5541 0.5883 4.13% MultiFab::AddProduct() 3439 0.5446 0.551 0.56 3.93% FillPatchIterator::Initialize 147 0.3073 0.4161 0.5553 3.90% FillPatchIterator::FillFromLevel0() 147 0.2832 0.3978 0.539 3.78% FillPatchSingleLevel 147 0.2825 0.3971 0.5383 3.78% Castro::reset_internal_energy(Fab) 160 0.5183 0.5255 0.5323 3.74% MLMG::mgVcycle_up::1 128 0.09387 0.2572 0.4207 2.95% MLPoisson::Fsmooth() 6144 0.3119 0.3392 0.4128 2.90% Amr::FinalizeInit() 1 0.2261 0.3329 0.3973 2.79% Castro::Castro() 1 0.3642 0.3723 0.3949 2.77% Castro::initMFs() 1 0.3591 0.367 0.3894 2.73% MLCellLinOp::applyBC() 7148 0.2685 0.3026 0.3741 2.63% MLMG::addInterpCorrection() 768 0.06285 0.3015 0.3712 2.61% MultiFab::Subtract() 159 0.346 0.3508 0.3554 2.49% FabArray::setVal() 5013 0.3355 0.3426 0.3528 2.48% Gravity::fill_multipole_BCs() 19 0.197 0.2378 0.2851 2.00% MLMG::mgVcycle_up::2 128 0.01076 0.1156 0.2755 1.93% FabArray::FillBoundary() 6380 0.2075 0.239 0.269 1.89% Castro::expand_state() 20 0.1782 0.2155 0.2652 1.86% MLMG::mgVcycle_down::0 128 0.2598 0.2617 0.2644 1.86% Castro::initData() 1 0.07741 0.1406 0.2551 1.79% check_for_negative_density() 20 0.05147 0.1572 0.2528 1.77% FabArray::ParallelCopy_nowait() 1429 0.199 0.2181 0.2389 1.68% Amr::bldFineLevels() 1 0.04691 0.1535 0.2175 1.53% MLCellLinOp::apply() 986 0.1491 0.1592 0.1952 1.37% MLMG::mgVcycle_up::0 128 0.1932 0.1937 0.1942 1.36% FabArray::sum() 3439 0.1825 0.1848 0.1877 1.32% MLMG::mgVcycle_up::3 128 0.01073 0.02356 0.1773 1.24% Amr::grid_places() 2 0.002958 0.1098 0.1742 1.22% AmrMesh::MakeNewGrids() 2 0.002833 0.1097 0.1741 1.22% TagBoxArray::mapPRD 2 0.0008229 0.1077 0.1726 1.21% MLMG:computeResOfCorrection() 768 0.1124 0.1221 0.1461 1.03% MLCellLinOp::correctionResidual() 768 0.111 0.1207 0.1446 1.01% FillBoundary_nowait() 6380 0.08693 0.1142 0.144 1.01% Castro::post_init() 1 0.1337 0.1356 0.1376 0.97% FillBoundary_finish() 6380 0.09108 0.1108 0.1352 0.95% Castro::construct_new_source() 102 0.1223 0.1244 0.1274 0.89% Castro::create_source_corrector() 17 0.04749 0.08274 0.1255 0.88% MLMG::mgVcycle_down::1 128 0.1135 0.1163 0.12 0.84% MLPoisson::Fapply() 986 0.09569 0.1005 0.1162 0.82% Gravity::multilevel_solve_for_new_phi() 1 0.1005 0.1005 0.1006 0.71% Gravity::actual_multilevel_solve() 1 0.1004 0.1005 0.1005 0.71% Amr::writePlotFile() 2 0.09985 0.09993 0.09994 0.70% MLCellLinOp::solutionResidual() 147 0.08672 0.09138 0.09563 0.67% Castro::post_regrid() 2 0.08551 0.08747 0.08943 0.63% Castro::construct_old_source() 120 0.07855 0.07957 0.0809 0.57% FabArray::Saxpy() 264 0.07176 0.07304 0.08054 0.57% Castro::construct_old_gravity() 20 0.04256 0.05751 0.07711 0.54% MLMG::computeResidual() 128 0.07506 0.07563 0.07629 0.54% Castro::derive() 18 0.06882 0.06973 0.07448 0.52% AmrLevel::derive() 18 0.06878 0.06969 0.07443 0.52% MLMG::mgVcycle_bottom 128 0.0003536 0.004854 0.07223 0.51% MLMG::actualBottomSolve() 128 0.0001305 0.004625 0.07198 0.51% MLMG::mgVcycle_down::2 128 0.01254 0.04065 0.07094 0.50% Castro::construct_new_gravity_source() 17 0.06711 0.06861 0.0708 0.50% MLCGSolver::bicgstab 8 0 0.004403 0.07045 0.49% Castro::post_timestep() 1 0.05999 0.06207 0.06512 0.46% Gravity::get_old_grav_vector() 20 0.04002 0.05452 0.06329 0.44% MLMG::mgVcycle_down::3 128 0.012 0.01823 0.06031 0.42% FabArray::Xpay() 945 0.04964 0.05211 0.05945 0.42% amrex::average_down 768 0.0332 0.03962 0.0549 0.39% Castro::initialize_advance() 1 0.04822 0.04946 0.05056 0.35% Castro::construct_old_gravity_source() 20 0.04541 0.04598 0.04677 0.33% Amr::regrid() 1 0.04392 0.04481 0.04562 0.32% Castro::check_for_nan() 37 0.04404 0.04459 0.04541 0.32% MultiFab::contains_nan() 37 0.04392 0.04447 0.04529 0.32% MLPoisson::define() 19 0.02904 0.03375 0.03998 0.28% Castro::retry_advance_ctu() 11 0.034 0.03539 0.03946 0.28% Gravity::get_new_grav_vector() 18 0.03414 0.03511 0.03701 0.26% MLMG::mgVcycle_down::4 128 0.01191 0.01357 0.035 0.25% StateDataPhysBCFunct::() 147 0.0006924 0.01771 0.03354 0.24% StateData::FillBoundary(geom) 143 0 0.01688 0.03268 0.23% Castro::apply_source_to_state() 37 0.02928 0.02963 0.03045 0.21% MLMG::mgVcycle_up::4 128 0.01083 0.01203 0.02731 0.19% VisMF::Write(FabArray) 4 0.01204 0.01343 0.02718 0.19% MLCellLinOp::defineAuxData() 19 0.01874 0.02042 0.0247 0.17% MLMG::mgVcycle_down::5 128 0.01186 0.0128 0.0232 0.16% MLMG::prepareForSolve() 18 0.02004 0.02089 0.02316 0.16% Castro::save_data_for_retry() 3 0.01809 0.01878 0.02294 0.16% Castro::sum_integrated_quantities() 2 0.0227 0.02286 0.02288 0.16% MLMG::mgVcycle_up::5 128 0.0109 0.01164 0.01977 0.14% MLMG::computeMLResidual() 18 0.01067 0.01526 0.01946 0.14% MLMG::ResNormInf() 146 0.01155 0.01429 0.01663 0.12% Castro::problem_post_timestep() 2 0.009181 0.01216 0.0161 0.11% Castro::wd_update() 2 0.009081 0.01214 0.01609 0.11% FabArray::norminf() 282 0.01082 0.01124 0.01545 0.11% FabArrayBase::getCPC() 2005 0.009672 0.01114 0.0149 0.10% FabArrayBase::getFB() 6380 0.01045 0.01131 0.01485 0.10% amrex::Add() 128 0.01418 0.01432 0.01446 0.10% Gravity::test_composite_phi() 1 0.014 0.01403 0.01406 0.10% MLLinOp::define() 19 0.004002 0.007362 0.01149 0.08% FabArray::setDomainBndry() 147 0.0004751 0.005901 0.01149 0.08% Other 2461 0.07911 0.09688 0.1414 0.99% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Pinned Memory Usage: --------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max --------------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 16 11 B 14 B 17 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.06-21-g259db7cfb99e) finalized