Initializing AMReX (25.06-22-g1d307c2c661a-dirty)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (25.06-22-g1d307c2c661a-dirty) initialized Starting run at 06:53:50 UTC on 2025-08-12. Successfully read inputs file ... Castro git describe: 25.08-12-gf8a308cc1 AMReX git describe: 25.06-22-g1d307c2c6-dirty Microphysics git describe: 25.08-3-gd0f79a2d 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.177049207 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001207974 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010618234 Gravity::multilevel_solve_for_new_phi() time = 0.109270149 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010664374 TIME= 0 MASS = 2.549996568e+33 TIME= 0 XMOM = 0 TIME= 0 YMOM = 3.48170636e+26 TIME= 0 ZMOM = 0 TIME= 0 ANG MOM X = 0 TIME= 0 ANG MOM Y = 0 TIME= 0 ANG MOM Z = 2.076918743e+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 = -3.641005016e-07 TIME= 0 CENTER OF MASS Y-VEL = 1.365376881e-07 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.048774475 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004123121 seconds [Level 0 step 1] ADVANCE at time 0 with dt = 0.25  Beginning subcycle 1 starting at time 0 with dt = 0.25 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -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.239071148 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.769409 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477044 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035502 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 = -0.01295059363 at index 21, 83, 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.7670713549 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.7670713549 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 = -0.01295059363 at index 21, 172, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035502 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477044 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.769409 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.239071148 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.00444886 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.503960463e+28 zmom added: 0 eint added: 0 ener added: 6.931398274e+47 Castro::do_old_sources() time = 0.065507438 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -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 Castro::construct_old_gravity() time = 0.005194945 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.251980232e+28 zmom added: 0 eint added: 0 ener added: 3.465699137e+47 Castro::do_old_sources() time = 0.065072873 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.125 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.0625   Beginning subcycle 1 starting at time 0 with dt = 0.0625 Estimated number of subcycles remaining: 4 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.1060744788 at index 20, 78, 0 Invalid density = -0.1060744788 at index 20, 177, 0 Castro::construct_old_gravity() time = 0.0052505 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -2.93534348e+40 ymom added: 1.625990116e+28 zmom added: 0 eint added: 0 ener added: 1.732849568e+47 Castro::do_old_sources() time = 0.061729626 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.0625 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.03125   Beginning subcycle 1 starting at time 0 with dt = 0.03125 Estimated number of subcycles remaining: 8 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Castro::construct_old_gravity() time = 0.005814149 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.129950579e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.059243378 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.337342614 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.015308882 Gravity::solve_for_phi() time = 0.092445587 on level 0 Castro::construct_new_gravity() time = 0.094630761 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.673145426e+44 ener added: 3.81935025e+42 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.752993387e+35 ymom added: -3.588950123e+27 zmom added: 0 eint added: 0 ener added: 6.222804276e+43 Castro::do_new_sources() time = 0.060173306 on level 0 Castro::construct_old_gravity() time = 0.002920391 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.129950579e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.058588843 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.337375085 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.015097625 Gravity::solve_for_phi() time = 0.092605351 on level 0 Castro::construct_new_gravity() time = 0.09472744 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.837104563e+44 ener added: 3.792000243e+42 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.742618664e+35 ymom added: -3.590636008e+27 zmom added: 0 eint added: 0 ener added: 4.583540443e+43 Castro::do_new_sources() time = 0.059594188 on level 0 Castro::construct_old_gravity() time = 0.003375554 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.129950579e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.058996191 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.343030263 on level 0 Ending SDC iteration 3 of 3.  Subcycle completed  Beginning subcycle 2 starting at time 0.03125 with dt = 0.03125 Estimated number of subcycles remaining: 7 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.015233009 Gravity::solve_for_phi() time = 0.092914317 on level 0 Castro::construct_new_gravity() time = 0.095572593 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.837201172e+44 ener added: 3.792181427e+42 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.742608901e+35 ymom added: -3.595677134e+27 zmom added: 0 eint added: 0 ener added: 4.58255499e+43 Castro::do_new_sources() time = 0.059614645 on level 0 Castro::construct_old_gravity() time = 0.003522435 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.985021797e+46 ener added: 5.427586462e+40 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467706592e+40 ymom added: 9.380508782e+26 zmom added: 0 eint added: 0 ener added: 8.721813599e+46 Castro::do_old_sources() time = 0.058700415 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.985569545e+46 ener added: 2.746016564e+42 Castro::construct_ctu_hydro_source() time = 0.346013812 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.015915424 Gravity::solve_for_phi() time = 0.093882427 on level 0 Castro::construct_new_gravity() time = 0.096004661 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.269601218e+44 ener added: 4.892522312e+43 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.763301424e+35 ymom added: 4.919525284e+24 zmom added: 0 eint added: 0 ener added: 6.012622799e+43 Castro::do_new_sources() time = 0.060053668 on level 0 Castro::construct_old_gravity() time = 0.003452191 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.985021797e+46 ener added: 5.427586462e+40 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467706592e+40 ymom added: 9.380508782e+26 zmom added: 0 eint added: 0 ener added: 8.721813599e+46 Castro::do_old_sources() time = 0.059167112 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.986269155e+46 ener added: 2.7526843e+42 Castro::construct_ctu_hydro_source() time = 0.353406863 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 3 starting at time 0.0625 with dt = 0.03125 Estimated number of subcycles remaining: 6 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.015470501 Gravity::solve_for_phi() time = 0.092997967 on level 0 Castro::construct_new_gravity() time = 0.095086588 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.343102187e+44 ener added: 4.888810292e+43 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.757196357e+35 ymom added: 6.54874172e+24 zmom added: 0 eint added: 0 ener added: 5.281428219e+43 Castro::do_new_sources() time = 0.059563902 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.910635222e+29 ener added: 1.910635222e+29 Castro::construct_old_gravity() time = 0.002975557 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.984353231e+46 ener added: 6.990279267e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467741736e+40 ymom added: 9.545224924e+26 zmom added: 0 eint added: 0 ener added: 8.780735363e+46 Castro::do_old_sources() time = 0.058987714 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.987321945e+46 ener added: 2.18047792e+43 Castro::construct_ctu_hydro_source() time = 0.347885452 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.015760343 Gravity::solve_for_phi() time = 0.094133339 on level 0 Castro::construct_new_gravity() time = 0.09625909 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.031202547e+43 ener added: 1.184458535e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.781734557e+35 ymom added: 5.277244545e+23 zmom added: 0 eint added: 0 ener added: 5.994705329e+43 Castro::do_new_sources() time = 0.061353026 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.910635222e+29 ener added: 1.910635222e+29 Castro::construct_old_gravity() time = 0.0028441 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.984353231e+46 ener added: 6.990279267e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467741736e+40 ymom added: 9.545224924e+26 zmom added: 0 eint added: 0 ener added: 8.780735363e+46 Castro::do_old_sources() time = 0.058764783 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.988028457e+46 ener added: 2.185423907e+43 Castro::construct_ctu_hydro_source() time = 0.354399231 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 4 starting at time 0.09375 with dt = 0.03125 Estimated number of subcycles remaining: 5 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.015649573 Gravity::solve_for_phi() time = 0.093242414 on level 0 Castro::construct_new_gravity() time = 0.095344122 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.773882032e+43 ener added: 1.183651189e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.77560788e+35 ymom added: 2.905435979e+24 zmom added: 0 eint added: 0 ener added: 5.260150257e+43 Castro::do_new_sources() time = 0.059677642 on level 0 Castro::construct_old_gravity() time = 0.002884777 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.983917963e+46 ener added: 1.154261437e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467777248e+40 ymom added: 9.587537328e+26 zmom added: 0 eint added: 0 ener added: 8.84101861e+46 Castro::do_old_sources() time = 0.058512358 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.984973203e+46 ener added: 2.291971287e+43 Castro::construct_ctu_hydro_source() time = 0.346285212 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.015288487 Gravity::solve_for_phi() time = 0.09258315 on level 0 Castro::construct_new_gravity() time = 0.094688132 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.612861023e+43 ener added: 1.488715317e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.801491314e+35 ymom added: 3.612610359e+23 zmom added: 0 eint added: 0 ener added: 4.637706889e+43 Castro::do_new_sources() time = 0.05968941 on level 0 Castro::construct_old_gravity() time = 0.002955949 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.983917963e+46 ener added: 1.154261437e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467777248e+40 ymom added: 9.587537328e+26 zmom added: 0 eint added: 0 ener added: 8.84101861e+46 Castro::do_old_sources() time = 0.058774297 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.985683204e+46 ener added: 2.297906432e+43 Castro::construct_ctu_hydro_source() time = 0.352137259 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 5 starting at time 0.125 with dt = 0.03125 Estimated number of subcycles remaining: 4 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.015234554 Gravity::solve_for_phi() time = 0.092382909 on level 0 Castro::construct_new_gravity() time = 0.094526801 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.357824616e+43 ener added: 1.487804348e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.795347886e+35 ymom added: -1.042462401e+25 zmom added: 0 eint added: 0 ener added: 3.901837814e+43 Castro::do_new_sources() time = 0.059821707 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.371823034e+29 ener added: 2.371823034e+29 Castro::construct_old_gravity() time = 0.00356744 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.983564871e+46 ener added: 1.521590942e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467813155e+40 ymom added: 9.391842461e+26 zmom added: 0 eint added: 0 ener added: 8.902573541e+46 Castro::do_old_sources() time = 0.058598548 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.992394606e+46 ener added: 2.321200641e+43 Castro::construct_ctu_hydro_source() time = 0.347013282 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.015631942 Gravity::solve_for_phi() time = 0.092941153 on level 0 Castro::construct_new_gravity() time = 0.095063075 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.122723364e+43 ener added: 1.619563371e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.826309838e+35 ymom added: 1.288025458e+24 zmom added: 0 eint added: 0 ener added: 2.081560288e+43 Castro::do_new_sources() time = 0.059686625 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.371823034e+29 ener added: 2.371823034e+29 Castro::construct_old_gravity() time = 0.002881174 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.983564871e+46 ener added: 1.521590942e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467813155e+40 ymom added: 9.391842461e+26 zmom added: 0 eint added: 0 ener added: 8.902573541e+46 Castro::do_old_sources() time = 0.058605405 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.993114585e+46 ener added: 2.326635407e+43 Castro::construct_ctu_hydro_source() time = 0.351964128 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 6 starting at time 0.15625 with dt = 0.03125 Estimated number of subcycles remaining: 3 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.015774649 Gravity::solve_for_phi() time = 0.093131472 on level 0 Castro::construct_new_gravity() time = 0.095245098 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.874849314e+43 ener added: 1.618719464e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.820136286e+35 ymom added: 2.815711015e+24 zmom added: 0 eint added: 0 ener added: 1.337821578e+43 Castro::do_new_sources() time = 0.05951615 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.138949907e+29 ener added: 5.138949907e+29 Castro::construct_old_gravity() time = 0.002855247 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.983077706e+46 ener added: 2.015109106e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467849558e+40 ymom added: 9.409976348e+26 zmom added: 0 eint added: 0 ener added: 8.964911027e+46 Castro::do_old_sources() time = 0.05857239 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.000666253e+46 ener added: 2.405830666e+43 Castro::construct_ctu_hydro_source() time = 0.347450402 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.016072777 Gravity::solve_for_phi() time = 0.09329223 on level 0 Castro::construct_new_gravity() time = 0.095400843 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.911535206e+43 ener added: 1.713555602e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.844610083e+35 ymom added: -2.677581796e+24 zmom added: 0 eint added: 0 ener added: -1.419860559e+43 Castro::do_new_sources() time = 0.059767391 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.138949907e+29 ener added: 5.138949907e+29 Castro::construct_old_gravity() time = 0.002873986 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.983077706e+46 ener added: 2.015109106e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467849558e+40 ymom added: 9.409976348e+26 zmom added: 0 eint added: 0 ener added: 8.964911027e+46 Castro::do_old_sources() time = 0.05859163 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.001390031e+46 ener added: 2.412276394e+43 Castro::construct_ctu_hydro_source() time = 0.352686883 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 7 starting at time 0.1875 with dt = 0.03125 Estimated number of subcycles remaining: 2 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.015684757 Gravity::solve_for_phi() time = 0.093933765 on level 0 Castro::construct_new_gravity() time = 0.096149046 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.667367417e+43 ener added: 1.712648092e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.83840701e+35 ymom added: 5.636144397e+24 zmom added: 0 eint added: 0 ener added: -2.166675325e+43 Castro::do_new_sources() time = 0.060149955 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: 3.038898262e+29 ener added: 3.038898262e+29 Castro::construct_old_gravity() time = 0.003147013 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.982914015e+46 ener added: 2.174744554e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467886326e+40 ymom added: 9.585270592e+26 zmom added: 0 eint added: 0 ener added: 9.027797292e+46 Castro::do_old_sources() time = 0.059505892 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.9938659e+46 ener added: 2.675763102e+43 Castro::construct_ctu_hydro_source() time = 0.391850174 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.01651908 Gravity::solve_for_phi() time = 0.102542862 on level 0 Castro::construct_new_gravity() time = 0.104692957 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.153618601e+44 ener added: 1.743191878e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.87433057e+35 ymom added: -3.282044706e+23 zmom added: 0 eint added: 0 ener added: -4.056999371e+43 Castro::do_new_sources() time = 0.060226587 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: 3.038898262e+29 ener added: 3.038898262e+29 Castro::construct_old_gravity() time = 0.002818178 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.982914015e+46 ener added: 2.174744554e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467886326e+40 ymom added: 9.585270592e+26 zmom added: 0 eint added: 0 ener added: 9.027797292e+46 Castro::do_old_sources() time = 0.062219122 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.99459825e+46 ener added: 2.682276214e+43 Castro::construct_ctu_hydro_source() time = 0.355797852 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed  Beginning subcycle 8 starting at time 0.21875 with dt = 0.03125 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016108761 Gravity::solve_for_phi() time = 0.093148688 on level 0 Castro::construct_new_gravity() time = 0.095292954 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.230005451e+44 ener added: 1.742291039e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.86809364e+35 ymom added: -8.38456169e+24 zmom added: 0 eint added: 0 ener added: -4.811917635e+43 Castro::do_new_sources() time = 0.059647335 on level 0 Castro::construct_old_gravity() time = 0.002857829 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.982547556e+46 ener added: 2.529310204e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467923688e+40 ymom added: 9.382019939e+26 zmom added: 0 eint added: 0 ener added: 9.091359292e+46 Castro::do_old_sources() time = 0.058604016 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.993424638e+46 ener added: 2.936823048e+43 Castro::construct_ctu_hydro_source() time = 0.369872542 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.016073821 Gravity::solve_for_phi() time = 0.093322139 on level 0 Castro::construct_new_gravity() time = 0.095438808 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.584986492e+43 ener added: 1.766412994e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.939839844e+35 ymom added: 3.243085182e+24 zmom added: 0 eint added: 0 ener added: -1.081736549e+43 Castro::do_new_sources() time = 0.059797944 on level 0 Castro::construct_old_gravity() time = 0.002947949 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.982547556e+46 ener added: 2.529310204e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467923688e+40 ymom added: 9.382019939e+26 zmom added: 0 eint added: 0 ener added: 9.091359292e+46 Castro::do_old_sources() time = 0.058668779 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.994157717e+46 ener added: 2.94327526e+43 Castro::construct_ctu_hydro_source() time = 0.353111899 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002651547848 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.01592728 Gravity::solve_for_phi() time = 0.093477889 on level 0 Castro::construct_new_gravity() time = 0.095654652 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.350350804e+43 ener added: 1.765529791e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.933604547e+35 ymom added: 1.827555829e+24 zmom added: 0 eint added: 0 ener added: -1.838318421e+43 Castro::do_new_sources() time = 0.059930281 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.581215356e+29 ener added: 1.581215356e+29 TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -4.341809003e+38 TIME= 0.25 YMOM = 1.164437349e+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.469212319e+38 TIME= 0.25 RHO*e = 1.325894512e+50 TIME= 0.25 RHO*K = 3.170941987e+49 TIME= 0.25 RHO*E = 1.642988711e+50 TIME= 0.25 RHO*PHI = -5.449011378e+50 TIME= 0.25 TOTAL ENERGY= -1.081516979e+50 TIME= 0.25 CENTER OF MASS X-LOC = 316424388.8 TIME= 0.25 CENTER OF MASS X-VEL = -170267.2489 TIME= 0.25 CENTER OF MASS Y-LOC = 2.427336677e-06 TIME= 0.25 CENTER OF MASS Y-VEL = 4.566427124e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 57600050.72 TIME= 0.25 MAXIMUM DENSITY = 4458760.788 TIME= 0.25 MAXIMUM T_S / T_E = 0 STEP = 1 TIME = 0.25 DT = 0.25 [STEP 1] Coarse TimeStep time: 12.41800488 [STEP 1] FAB kilobyte spread across MPI nodes: [6998 ... 6998] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.076054206 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004134466 seconds Ending run at 06:54:05 UTC on 2025-08-12. Run time = 14.90452112 Run time without initialization = 12.49905768 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.8956 ... 14.8986 ... 14.9046 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.7085 5.8207 5.9362 39.83% Amr::InitAmr() 1 1.487 1.49 1.496 10.04% Castro::computeTemp() 80 0.8258 0.9482 1.091 7.32% FabArray::ParallelCopy_finish() 1429 0.2158 0.5554 0.8028 5.39% amrex::Copy() 724 0.7598 0.7703 0.7822 5.25% MultiFab::AddProduct() 3439 0.5731 0.59 0.7241 4.86% Castro::normalize_species() 79 0.5443 0.5521 0.56 3.76% Castro::reset_internal_energy(Fab) 160 0.5151 0.524 0.5333 3.58% MLPoisson::Fsmooth() 6144 0.3528 0.386 0.4708 3.16% Castro::initMFs() 1 0.3573 0.3627 0.3873 2.60% FabArray::setVal() 5013 0.3536 0.3607 0.3697 2.48% MultiFab::Subtract() 159 0.3419 0.3469 0.3523 2.36% check_for_negative_density() 20 0.09026 0.2059 0.3333 2.24% Gravity::fill_multipole_BCs() 19 0.1972 0.2547 0.2846 1.91% FabArray::ParallelCopy_nowait() 1429 0.2074 0.2295 0.2518 1.69% FabArray::sum() 3439 0.1784 0.1815 0.1848 1.24% Castro::initData() 1 0.003476 0.02877 0.1694 1.14% FillBoundary_finish() 6380 0.1051 0.1325 0.1591 1.07% FillBoundary_nowait() 6380 0.08947 0.1222 0.1565 1.05% TagBoxArray::collate() 2 0.0002667 0.03037 0.1448 0.97% MLCellLinOp::applyBC() 7148 0.06392 0.08068 0.1291 0.87% MLPoisson::Fapply() 986 0.09513 0.1003 0.117 0.79% Castro::construct_new_gravity_source() 17 0.07309 0.07452 0.07586 0.51% FabArray::Saxpy() 91 0.07231 0.07305 0.07391 0.50% FabArray::Xpay() 915 0.05025 0.05293 0.05929 0.40% Castro::construct_old_gravity_source() 20 0.04865 0.04965 0.05189 0.35% MultiFab::contains_nan() 37 0.0473 0.04808 0.04922 0.33% Amr::writePlotFile() 2 0.03493 0.04356 0.04502 0.30% MLMG::addInterpCorrection() 768 0.03357 0.03644 0.0447 0.30% AmrLevel::derive() 18 0.04072 0.04125 0.0422 0.28% amrex::average_down 768 0.02521 0.02764 0.03487 0.23% StateData::FillBoundary(geom) 143 0 0.0175 0.03356 0.23% VisMF::Write(FabArray) 4 0.01389 0.01495 0.02473 0.17% Castro::reset_internal_energy(MultiFab) 80 0.02338 0.0239 0.02458 0.16% MLCellLinOp::defineAuxData() 19 0.01584 0.01676 0.02053 0.14% Castro::enforce_min_density() 79 0.01854 0.01881 0.01924 0.13% FabArray::FillBoundary() 6380 0.01359 0.01424 0.01749 0.12% Castro::construct_new_source() 102 0.01623 0.01645 0.01673 0.11% MLCGSolver::bicgstab 8 0 0.001035 0.01656 0.11% FabArray::norminf() 244 0.01151 0.01199 0.01507 0.10% Castro::wd_update() 2 0.004998 0.01017 0.015 0.10% amrex::Add() 128 0.01399 0.01424 0.01448 0.10% Castro::evaluate_source_change() 196 0.01356 0.01404 0.01443 0.10% FabArray::setDomainBndry() 147 0.0004857 0.007557 0.01432 0.10% FabArrayBase::getFB() 6380 0.01067 0.01115 0.01376 0.09% FillPatchIterator::Initialize 147 0.01278 0.0131 0.01338 0.09% MLMG::prepareForSolve() 18 0.0102 0.01062 0.01186 0.08% MLCellLinOp::smooth() 1544 0.00976 0.01025 0.0114 0.08% Castro::construct_old_source() 120 0.01033 0.01048 0.0107 0.07% MLMG::ResNormInf() 146 0.003568 0.006024 0.009406 0.06% Gravity::actual_solve_with_mlmg() 19 0.007338 0.007806 0.009275 0.06% MLMG::solve() 18 0.0007858 0.004622 0.008929 0.06% FabArrayBase::getCPC() 2005 0.007665 0.007953 0.008449 0.06% Gravity::get_old_grav_vector() 20 0.006612 0.006774 0.007249 0.05% Gravity::get_new_grav_vector() 18 0.006039 0.00626 0.006713 0.05% Castro::initialize_do_advance() 20 0.004588 0.004692 0.0049 0.03% MLCellLinOp::apply() 986 0.002026 0.002251 0.004874 0.03% FabArray::ParallelCopy() 1429 0.003589 0.00376 0.004111 0.03% AmrLevel::FillPatch() 147 0.002416 0.002494 0.002565 0.02% Castro::sum_integrated_quantities() 2 0.001563 0.001971 0.002276 0.02% MLCellLinOp::correctionResidual() 768 0.001922 0.001996 0.0021 0.01% MLMG::mgVcycle() 128 0.001552 0.00165 0.001876 0.01% MLMG:computeResOfCorrection() 768 0.001483 0.001548 0.001658 0.01% Castro::save_data_for_retry() 3 0.001267 0.001309 0.00141 0.01% MLMG::actualBottomSolve() 128 0.0001323 0.0002055 0.001239 0.01% StateDataPhysBCFunct::() 147 0.0007144 0.0008679 0.001027 0.01% Castro::subcycle_advance_ctu() 1 0.0002226 0.0002971 0.0008929 0.01% main() 1 0.0006635 0.000719 0.0007879 0.01% FillPatchIterator::FillFromLevel0() 147 0.0007009 0.0007204 0.0007401 0.00% Gravity::solve_for_phi() 17 0.0006564 0.0006838 0.0007193 0.00% MLMG::mgVcycle_down::2 128 0.0005244 0.0005719 0.0006153 0.00% MLMG::mgVcycle_down::3 128 0.0004947 0.0005379 0.0006135 0.00% MLMG::mgVcycle_down::4 128 0.0004821 0.0005028 0.0006113 0.00% AmrMesh::MakeNewGrids() 2 0.0005633 0.0005822 0.0006089 0.00% Castro::Castro() 1 0.0003981 0.0004537 0.0006045 0.00% MLMG::mgVcycle_down::0 128 0.0005467 0.00057 0.0005962 0.00% MLMG::mgVcycle_down::1 128 0.0005487 0.0005685 0.0005944 0.00% MLMG::mgVcycle_down::5 128 0.000498 0.0005181 0.0005819 0.00% FillPatchSingleLevel 147 0.0004551 0.0004752 0.0005022 0.00% Castro::clean_state() 79 0.0004348 0.0004576 0.0004901 0.00% MLCellLinOp::solutionResidual() 147 0.0004594 0.0004705 0.0004885 0.00% Castro::do_advance_ctu() 20 0.000387 0.0004143 0.0004463 0.00% MLMG::oneIter() 128 0.0003904 0.0004038 0.0004264 0.00% MLMG::mgVcycle_up::1 128 0.000351 0.0003801 0.0004043 0.00% MLMG::mgVcycle_up::0 128 0.0003626 0.0003752 0.0003911 0.00% MLMG::mgVcycle_up::2 128 0.0003079 0.0003403 0.0003827 0.00% MLMG::mgVcycle_up::4 128 0.0003093 0.0003218 0.0003597 0.00% MLMG::mgVcycle_up::3 128 0.0003039 0.0003182 0.0003571 0.00% MLMG::mgVcycle_up::5 128 0.0003073 0.0003193 0.0003419 0.00% Castro::do_old_sources() 20 0.0002782 0.000295 0.0003125 0.00% Castro::problem_post_timestep() 2 1.162e-05 3.112e-05 0.0003048 0.00% Castro::construct_new_gravity() 17 0.0001689 0.0001868 0.0002895 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002203 0.0002381 0.0002636 0.00% MLMG::computeResidual() 128 0.0002344 0.0002429 0.0002488 0.00% MLMG::mgVcycle_bottom 128 0.0002212 0.0002294 0.0002475 0.00% Castro::do_new_sources() 17 0.0002136 0.0002244 0.0002461 0.00% Amr::defBaseLevel() 1 0.0001149 0.0001359 0.0001574 0.00% Amr::coarseTimeStep() 1 0.0001211 0.0001462 0.000151 0.00% TagBoxArray::mapPRD 2 0.0001229 0.0001305 0.0001403 0.00% MLPoisson::define() 19 0.000119 0.0001277 0.0001362 0.00% Castro::check_for_nan() 37 0.0001066 0.0001161 0.0001312 0.00% Castro::initialize_advance() 1 0.0001041 0.0001107 0.0001303 0.00% Castro::retry_advance_ctu() 11 6.287e-05 7.343e-05 0.00013 0.00% Gravity::test_composite_phi() 1 0.0001052 0.0001114 0.0001283 0.00% Castro::construct_old_gravity() 20 0.0001094 0.0001158 0.0001256 0.00% Castro::apply_source_to_state() 37 0.0001031 0.0001067 0.0001125 0.00% MLLinOp::define() 19 7.414e-05 7.935e-05 8.375e-05 0.00% Amr::grid_places() 2 5.311e-05 6.197e-05 7.837e-05 0.00% Castro::create_source_corrector() 17 5.699e-05 6.179e-05 6.745e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.118e-05 6.232e-05 6.325e-05 0.00% Castro::expand_state() 20 4.9e-05 5.244e-05 6.268e-05 0.00% Amr::FinalizeInit() 1 8.51e-06 1.57e-05 5.143e-05 0.00% MLMG::computeMLResidual() 18 4.202e-05 4.409e-05 4.851e-05 0.00% Castro::derive() 18 4.107e-05 4.28e-05 4.524e-05 0.00% Amr::regrid() 1 1.662e-05 2.037e-05 4.324e-05 0.00% Gravity::actual_multilevel_solve() 1 3.234e-05 3.3e-05 3.358e-05 0.00% Amr::init() 1 6.097e-06 8.394e-06 2.871e-05 0.00% Castro::post_regrid() 2 1.062e-05 1.275e-05 2.83e-05 0.00% Castro::post_init() 1 1.686e-05 1.901e-05 2.343e-05 0.00% Amr::timeStep() 1 1.439e-05 1.701e-05 2.131e-05 0.00% Amr::InitializeInit() 1 8.297e-06 1.139e-05 1.739e-05 0.00% Castro::advance() 1 6.93e-06 7.914e-06 1.074e-05 0.00% Amr::bldFineLevels() 1 6.872e-06 7.981e-06 1.051e-05 0.00% Castro::post_timestep() 1 5.799e-06 6.247e-06 7.393e-06 0.00% Amr::initialInit() 1 2.949e-06 3.205e-06 3.435e-06 0.00% Other 2521 0.04948 0.06432 0.09879 0.66% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.9 14.9 14.9 100.00% Amr::coarseTimeStep() 1 12.5 12.5 12.5 83.86% Amr::timeStep() 1 12.42 12.42 12.42 83.32% Castro::advance() 1 12.35 12.36 12.36 82.91% Castro::subcycle_advance_ctu() 1 12.3 12.31 12.31 82.58% Castro::do_advance_ctu() 20 12.26 12.27 12.27 82.33% Castro::construct_ctu_hydro_source() 20 6.785 6.809 6.853 45.98% Castro::clean_state() 79 3.63 3.796 3.93 26.37% Castro::computeTemp() 80 2.315 2.461 2.583 17.33% Gravity::solve_phi_with_mlmg() 19 1.62 1.679 1.713 11.49% Castro::construct_new_gravity() 17 1.538 1.596 1.63 10.94% Castro::reset_internal_energy(MultiFab) 80 1.486 1.513 1.593 10.69% Gravity::solve_for_phi() 17 1.5 1.559 1.593 10.69% Amr::InitAmr() 1 1.487 1.49 1.496 10.04% Castro::initialize_do_advance() 20 1.397 1.431 1.486 9.97% Gravity::actual_solve_with_mlmg() 19 1.417 1.421 1.425 9.56% MLMG::solve() 18 1.361 1.367 1.369 9.19% MLMG::oneIter() 128 1.224 1.227 1.231 8.26% MLMG::mgVcycle() 128 1.209 1.213 1.217 8.16% Castro::do_old_sources() 20 1.157 1.176 1.192 8.00% Castro::evaluate_source_change() 196 0.968 0.9908 1.132 7.60% FabArray::ParallelCopy() 1429 0.471 0.8086 1.05 7.05% Castro::do_new_sources() 17 0.918 0.9775 1.017 6.82% Amr::init() 1 0.9093 0.9093 0.9093 6.10% Castro::enforce_min_density() 79 0.7948 0.8095 0.88 5.90% Amr::initialInit() 1 0.8559 0.8559 0.856 5.74% MLCellLinOp::smooth() 1544 0.6304 0.6875 0.8318 5.58% FabArray::ParallelCopy_finish() 1429 0.2158 0.5554 0.8028 5.39% amrex::Copy() 724 0.7598 0.7703 0.7822 5.25% AmrLevel::FillPatch() 147 0.5173 0.6037 0.7305 4.90% MultiFab::AddProduct() 3439 0.5731 0.59 0.7241 4.86% Amr::InitializeInit() 1 0.443 0.4746 0.6185 4.15% Amr::defBaseLevel() 1 0.443 0.4745 0.6185 4.15% Castro::normalize_species() 79 0.5443 0.5521 0.56 3.76% FillPatchIterator::Initialize 147 0.3384 0.4277 0.5568 3.74% Castro::reset_internal_energy(Fab) 160 0.5151 0.524 0.5333 3.58% FillPatchIterator::FillFromLevel0() 147 0.3196 0.4071 0.53 3.56% FillPatchSingleLevel 147 0.3188 0.4064 0.5293 3.55% MLPoisson::Fsmooth() 6144 0.3528 0.386 0.4708 3.16% MLMG::mgVcycle_up::1 128 0.1073 0.2879 0.4684 3.14% MLCellLinOp::applyBC() 7148 0.3193 0.3616 0.4405 2.96% Amr::FinalizeInit() 1 0.2375 0.3814 0.413 2.77% MLMG::addInterpCorrection() 768 0.06824 0.3292 0.4064 2.73% Castro::Castro() 1 0.3627 0.3684 0.3929 2.64% Castro::initMFs() 1 0.3574 0.3629 0.3875 2.60% FabArray::setVal() 5013 0.3536 0.3607 0.3697 2.48% MultiFab::Subtract() 159 0.3419 0.3469 0.3523 2.36% check_for_negative_density() 20 0.09026 0.2059 0.3333 2.24% FabArray::FillBoundary() 6380 0.2407 0.2809 0.3322 2.23% MLMG::mgVcycle_up::2 128 0.01356 0.1289 0.3051 2.05% Gravity::fill_multipole_BCs() 19 0.1993 0.2567 0.2866 1.92% MLMG::mgVcycle_down::0 128 0.2771 0.281 0.2832 1.90% Castro::expand_state() 20 0.1971 0.2269 0.2807 1.88% FabArray::ParallelCopy_nowait() 1429 0.2265 0.2495 0.2722 1.83% Castro::initData() 1 0.07842 0.1044 0.2453 1.65% Amr::bldFineLevels() 1 0.04664 0.1904 0.2215 1.49% MLMG::mgVcycle_up::0 128 0.2171 0.219 0.2206 1.48% MLCellLinOp::apply() 986 0.1599 0.1709 0.211 1.42% MLMG::mgVcycle_up::3 128 0.01353 0.02777 0.1969 1.32% FabArray::sum() 3439 0.1784 0.1815 0.1848 1.24% Amr::grid_places() 2 0.002738 0.1466 0.1782 1.20% AmrMesh::MakeNewGrids() 2 0.002664 0.1466 0.1782 1.20% TagBoxArray::mapPRD 2 0.0007986 0.1148 0.1764 1.18% FillBoundary_nowait() 6380 0.1005 0.1341 0.1689 1.13% FillBoundary_finish() 6380 0.1051 0.1325 0.1591 1.07% MLMG:computeResOfCorrection() 768 0.1175 0.1292 0.1555 1.04% MLCellLinOp::correctionResidual() 768 0.116 0.1277 0.1538 1.03% Castro::post_init() 1 0.143 0.1476 0.1488 1.00% TagBoxArray::collate() 2 0.0002667 0.03037 0.1448 0.97% MLMG::mgVcycle_down::1 128 0.126 0.1297 0.1347 0.90% Castro::construct_new_source() 102 0.1299 0.1318 0.1337 0.90% Amr::writePlotFile() 2 0.1252 0.1252 0.1253 0.84% MLPoisson::Fapply() 986 0.09513 0.1003 0.117 0.79% Gravity::multilevel_solve_for_new_phi() 1 0.1092 0.1092 0.1093 0.73% Gravity::actual_multilevel_solve() 1 0.1091 0.1091 0.1092 0.73% MLCellLinOp::solutionResidual() 147 0.09121 0.09685 0.104 0.70% Castro::create_source_corrector() 17 0.05211 0.07265 0.09473 0.64% Castro::post_regrid() 2 0.08584 0.08706 0.09169 0.62% Castro::construct_old_source() 120 0.08285 0.08421 0.0867 0.58% MLMG::mgVcycle_down::2 128 0.01586 0.04723 0.0817 0.55% MLMG::computeResidual() 128 0.0773 0.07954 0.08025 0.54% Castro::construct_new_gravity_source() 17 0.07309 0.07452 0.07586 0.51% FabArray::Saxpy() 91 0.07231 0.07305 0.07391 0.50% MLMG::mgVcycle_bottom 128 0.0003538 0.004897 0.07289 0.49% MLMG::actualBottomSolve() 128 0.0001323 0.004668 0.07264 0.49% Castro::derive() 18 0.06882 0.06973 0.07111 0.48% AmrLevel::derive() 18 0.06878 0.06969 0.07107 0.48% MLCGSolver::bicgstab 8 0 0.004422 0.07075 0.47% MLMG::mgVcycle_down::3 128 0.0152 0.02207 0.06876 0.46% Castro::post_timestep() 1 0.05988 0.06228 0.06601 0.44% Castro::construct_old_gravity() 20 0.04194 0.0573 0.06426 0.43% Gravity::get_old_grav_vector() 20 0.04181 0.05719 0.06415 0.43% amrex::average_down 768 0.03433 0.04149 0.05935 0.40% FabArray::Xpay() 915 0.05025 0.05293 0.05929 0.40% Castro::construct_old_gravity_source() 20 0.04865 0.04965 0.05189 0.35% Castro::initialize_advance() 1 0.04893 0.04972 0.05157 0.35% Castro::check_for_nan() 37 0.04741 0.04819 0.04935 0.33% MultiFab::contains_nan() 37 0.0473 0.04808 0.04922 0.33% Amr::regrid() 1 0.04439 0.04491 0.04673 0.31% MLPoisson::define() 19 0.03061 0.03597 0.04164 0.28% MLMG::mgVcycle_down::4 128 0.01506 0.01692 0.04066 0.27% Gravity::get_new_grav_vector() 18 0.03615 0.03711 0.03829 0.26% Castro::retry_advance_ctu() 11 0.03525 0.03631 0.03827 0.26% StateDataPhysBCFunct::() 147 0.0007144 0.01837 0.03448 0.23% StateData::FillBoundary(geom) 143 0 0.0175 0.03356 0.23% MLMG::mgVcycle_up::4 128 0.01358 0.01496 0.03239 0.22% Castro::apply_source_to_state() 37 0.02949 0.02978 0.03014 0.20% MLMG::mgVcycle_down::5 128 0.015 0.0161 0.02845 0.19% MLCellLinOp::defineAuxData() 19 0.02 0.0218 0.02786 0.19% MLMG::prepareForSolve() 18 0.02051 0.02165 0.02484 0.17% VisMF::Write(FabArray) 4 0.01389 0.01495 0.02473 0.17% MLMG::mgVcycle_up::5 128 0.01369 0.01455 0.02437 0.16% MLMG::computeMLResidual() 18 0.01104 0.01679 0.02301 0.15% Castro::sum_integrated_quantities() 2 0.02117 0.02136 0.02138 0.14% MLMG::ResNormInf() 146 0.01399 0.01629 0.01952 0.13% Castro::save_data_for_retry() 3 0.01809 0.01853 0.01912 0.13% Castro::problem_post_timestep() 2 0.01045 0.01394 0.01836 0.12% Castro::wd_update() 2 0.01044 0.01391 0.01835 0.12% FabArrayBase::getCPC() 2005 0.01006 0.01156 0.01571 0.11% FabArrayBase::getFB() 6380 0.01102 0.0119 0.01542 0.10% FabArray::norminf() 244 0.01151 0.01199 0.01507 0.10% Gravity::test_composite_phi() 1 0.0144 0.01445 0.01448 0.10% amrex::Add() 128 0.01399 0.01424 0.01448 0.10% FabArray::setDomainBndry() 147 0.0004857 0.007557 0.01432 0.10% MLLinOp::define() 19 0.003951 0.007477 0.01116 0.07% Other 2521 0.08295 0.1056 0.1488 1.00% -------------------------------------------------------------------------------------------- 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 12 B 16 B 27 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 1559 KiB 1560 KiB 1560 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 402 KiB 410 KiB 419 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1033 KiB 1033 KiB 1033 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 983 KiB 983 KiB 983 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 660 KiB 660 KiB 660 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 4120 B 4198 B 4681 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 48 KiB 49 KiB 53 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 11 KiB 12 KiB 15 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 412 KiB 412 KiB 412 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 2147 B 3844 B 5361 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1287 B 1302 B 1330 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 16 KiB 17 KiB 18 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 130 KiB 130 KiB 130 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 192 127 KiB 127 KiB 127 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 6102 9604 B 9 KiB 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 255 B 258 B 261 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 429 B 434 B 439 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 17 B 20 B 36 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 50 B 50 B 51 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 3 B 8 B 9 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 109 B 112 B 116 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 128 B 175 B 197 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 225 B 232 B 237 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1059 B 1084 B 1240 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 119 B 119 B 120 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1649 B 1713 B 1750 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 434 B 484 B 673 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 451 B 452 B 453 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 24 B 33 B 44 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 19 B 30 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 16 B 25 B 1024 B 1024 B 1024 B amrex::average_down 5632 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 64 0 B 5 B 7 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 (25.06-22-g1d307c2c661a-dirty) finalized