Initializing AMReX (25.11-46-gce6502b2c0cb)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (25.11-46-gce6502b2c0cb) initialized Starting run at 07:53:19 UTC on 2025-11-21. Successfully read inputs file ... Castro git describe: 25.11-7-gc9f724d8e AMReX git describe: 25.11-46-gce6502b2c Microphysics git describe: 25.11-7-gb1a752a3 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.174626314 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001224786 new finest: 0 Gravity::fill_multipole_BCs() time = 0.011036065 Gravity::multilevel_solve_for_new_phi() time = 0.10810184 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010924486 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.051728656 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004317251 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 = -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 = -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 = -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.002810201 on level 0 Castro::do_old_sources() time = 0.007690317 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.065561805 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.004450588 on level 0 Castro::do_old_sources() time = 0.008335444 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.064555345 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.005008045 on level 0 Castro::do_old_sources() time = 0.008735422 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.061170657 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.003409749 on level 0 Castro::do_old_sources() time = 0.007947392 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.059932342 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.333874956 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.014084786 Gravity::solve_for_phi() time = 0.09033908 on level 0 Castro::construct_new_gravity() time = 0.093036589 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.060827926 on level 0 Castro::construct_old_gravity() time = 0.003353035 on level 0 Castro::do_old_sources() time = 0.00778203 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.060065446 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.332919617 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.013856333 Gravity::solve_for_phi() time = 0.090257929 on level 0 Castro::construct_new_gravity() time = 0.092685626 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.060901321 on level 0 Castro::construct_old_gravity() time = 0.003043398 on level 0 Castro::do_old_sources() time = 0.00777297 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.060183172 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.335465093 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.013681775 Gravity::solve_for_phi() time = 0.09070098 on level 0 Castro::construct_new_gravity() time = 0.09286671 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.061328864 on level 0 Castro::construct_old_gravity() time = 0.00395517 on level 0 Castro::do_old_sources() time = 0.00789001 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.059798598 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.341163905 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.015530274 Gravity::solve_for_phi() time = 0.092697604 on level 0 Castro::construct_new_gravity() time = 0.094866263 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.060655289 on level 0 Castro::construct_old_gravity() time = 0.002814681 on level 0 Castro::do_old_sources() time = 0.007806622 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.060315461 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.345230355 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.014189257 Gravity::solve_for_phi() time = 0.092528604 on level 0 Castro::construct_new_gravity() time = 0.094662432 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.060691269 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.002794767 on level 0 Castro::do_old_sources() time = 0.007700393 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.059755548 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.339820421 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.014180535 Gravity::solve_for_phi() time = 0.090625719 on level 0 Castro::construct_new_gravity() time = 0.092805216 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.061094873 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.002849979 on level 0 Castro::do_old_sources() time = 0.007755095 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.060232314 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.345598154 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.014212481 Gravity::solve_for_phi() time = 0.091383213 on level 0 Castro::construct_new_gravity() time = 0.093524768 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.061142543 on level 0 Castro::construct_old_gravity() time = 0.003605274 on level 0 Castro::do_old_sources() time = 0.00782281 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.059842615 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.340355436 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.01463093 Gravity::solve_for_phi() time = 0.091902389 on level 0 Castro::construct_new_gravity() time = 0.094043846 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.061527045 on level 0 Castro::construct_old_gravity() time = 0.002796827 on level 0 Castro::do_old_sources() time = 0.007674208 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.06018983 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.346307104 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.016877802 Gravity::solve_for_phi() time = 0.095969545 on level 0 Castro::construct_new_gravity() time = 0.098199182 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.060370489 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.003084165 on level 0 Castro::do_old_sources() time = 0.008030869 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.060373608 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.341299962 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.014433681 Gravity::solve_for_phi() time = 0.090868867 on level 0 Castro::construct_new_gravity() time = 0.093045884 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.060570548 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.002715474 on level 0 Castro::do_old_sources() time = 0.007710801 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.059916264 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.345466729 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.01486961 Gravity::solve_for_phi() time = 0.091407664 on level 0 Castro::construct_new_gravity() time = 0.09362072 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.060187036 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.002801018 on level 0 Castro::do_old_sources() time = 0.007701519 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.060187392 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.341309168 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.014407035 Gravity::solve_for_phi() time = 0.091532322 on level 0 Castro::construct_new_gravity() time = 0.09368774 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.060618886 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.002756707 on level 0 Castro::do_old_sources() time = 0.007741396 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.059850976 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.347336107 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.014618424 Gravity::solve_for_phi() time = 0.091242542 on level 0 Castro::construct_new_gravity() time = 0.094104313 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.060315478 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.002777985 on level 0 Castro::do_old_sources() time = 0.007702003 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.060213014 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.364792016 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.018951838 Gravity::solve_for_phi() time = 0.095939258 on level 0 Castro::construct_new_gravity() time = 0.098089596 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.06073149 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.002716223 on level 0 Castro::do_old_sources() time = 0.008041614 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.059994654 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.348659065 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.014617199 Gravity::solve_for_phi() time = 0.091459975 on level 0 Castro::construct_new_gravity() time = 0.093670207 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.060811984 on level 0 Castro::construct_old_gravity() time = 0.002806585 on level 0 Castro::do_old_sources() time = 0.007699682 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.059696142 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.341393845 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.01482129 Gravity::solve_for_phi() time = 0.092092055 on level 0 Castro::construct_new_gravity() time = 0.094219163 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.060585861 on level 0 Castro::construct_old_gravity() time = 0.002823839 on level 0 Castro::do_old_sources() time = 0.007712494 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.060192573 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.354831602 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.0025912558 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.01866274 Gravity::solve_for_phi() time = 0.096405684 on level 0 Castro::construct_new_gravity() time = 0.09855409 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.067635383 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.137357411e+28 TIME= 0.25 ZMOM = 0 TIME= 0.25 ANG MOM X = 0 TIME= 0.25 ANG MOM Y = 0 TIME= 0.25 ANG MOM Z = -1.469316165e+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.460231145e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 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.70958976 [STEP 1] FAB kilobyte spread across MPI nodes: [7155 ... 7155] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.049589993 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004225396 seconds Ending run at 07:53:34 UTC on 2025-11-21. Run time = 15.05145897 Run time without initialization = 12.76435483 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]: 15.0439 ... 15.0457 ... 15.0516 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.6873 5.7192 5.7961 38.51% Amr::InitAmr() 1 1.368 1.37 1.376 9.14% Castro::computeTemp() 80 0.8268 0.9342 1.089 7.24% amrex::Copy() 766 0.8713 0.8757 0.898 5.97% FabArray::ParallelCopy_finish() 1469 0.1941 0.5322 0.7306 4.85% MultiFab::AddProduct() 3598 0.5938 0.596 0.6045 4.02% Castro::normalize_species() 79 0.5499 0.5521 0.563 3.74% Castro::reset_internal_energy(Fab) 160 0.5128 0.5163 0.5252 3.49% MLPoisson::Fsmooth() 6144 0.3499 0.3793 0.4551 3.02% FabArray::setVal() 5232 0.3735 0.3764 0.3852 2.56% Castro::initMFs() 1 0.358 0.3623 0.3835 2.55% MultiFab::Subtract() 159 0.3496 0.352 0.3574 2.37% FabArray::ParallelCopy_nowait() 1469 0.2335 0.254 0.2762 1.84% Castro::do_advance_ctu() 20 0.2701 0.2727 0.2762 1.83% Gravity::fill_multipole_BCs() 19 0.2033 0.2522 0.2737 1.82% FabArray::sum() 3598 0.2052 0.2061 0.2095 1.39% Castro::initData() 1 0.003343 0.02858 0.1693 1.13% TagBoxArray::collate() 2 0.0003063 0.03216 0.1649 1.10% check_for_negative_density() 20 0.05703 0.1181 0.1649 1.10% FillBoundary_finish() 6380 0.09338 0.1285 0.1575 1.05% FillBoundary_nowait() 6380 0.08664 0.1164 0.1502 1.00% MLCellLinOp::applyBC() 7148 0.06736 0.08281 0.1298 0.86% MLPoisson::Fapply() 986 0.09589 0.1 0.1164 0.77% Castro::construct_old_gravity_source() 40 0.09891 0.09919 0.1002 0.67% FabArray::Saxpy() 111 0.08521 0.08564 0.0871 0.58% Castro::construct_new_gravity_source() 17 0.07321 0.074 0.07569 0.50% FabArray::Xpay() 915 0.05051 0.05245 0.05793 0.38% MultiFab::contains_nan() 37 0.04757 0.0478 0.04831 0.32% MLMG::addInterpCorrection() 768 0.03367 0.03613 0.04375 0.29% AmrLevel::derive() 18 0.04025 0.04059 0.04127 0.27% StateData::FillBoundary(geom) 193 0 0.01981 0.03809 0.25% amrex::average_down 768 0.0247 0.02686 0.03356 0.22% VisMF::Write(FabArray) 4 0.01404 0.01519 0.02702 0.18% Castro::reset_internal_energy(MultiFab) 80 0.02408 0.0245 0.02514 0.17% Castro::construct_old_source() 240 0.02029 0.02048 0.02129 0.14% Castro::wd_update() 2 0.005326 0.01792 0.02083 0.14% FabArray::FillBoundary() 6380 0.01665 0.0171 0.02027 0.13% MLCellLinOp::defineAuxData() 19 0.0154 0.01638 0.01949 0.13% Castro::enforce_min_density() 79 0.01874 0.01891 0.0191 0.13% Amr::writePlotFile() 2 0.006755 0.01762 0.01889 0.13% MLCGSolver::bicgstab 8 0 0.001053 0.01685 0.11% FabArray::setDomainBndry() 187 0.0006436 0.008874 0.01675 0.11% Castro::construct_new_source() 102 0.01598 0.01614 0.01666 0.11% FabArrayBase::getFB() 6380 0.01309 0.01359 0.01647 0.11% Castro::evaluate_source_change() 196 0.0142 0.0146 0.01571 0.10% FillPatchIterator::Initialize 187 0.01499 0.01522 0.01552 0.10% FabArray::norminf() 244 0.01095 0.01138 0.01477 0.10% amrex::Add() 128 0.0141 0.01421 0.01454 0.10% MLCellLinOp::smooth() 1544 0.01076 0.01114 0.01259 0.08% MLMG::prepareForSolve() 18 0.01004 0.01043 0.01157 0.08% Gravity::actual_solve_with_mlmg() 19 0.007321 0.007794 0.009445 0.06% FabArrayBase::getCPC() 2045 0.008385 0.008563 0.009032 0.06% MLMG::ResNormInf() 146 0.003579 0.005875 0.008903 0.06% Gravity::get_new_grav_vector() 18 0.006034 0.006239 0.007226 0.05% Gravity::get_old_grav_vector() 20 0.006517 0.006592 0.006766 0.04% MLCellLinOp::apply() 986 0.002425 0.002693 0.005806 0.04% MLMG::solve() 18 0.0007563 0.002872 0.005188 0.03% Castro::initialize_do_advance() 20 0.004726 0.004894 0.005144 0.03% FabArray::ParallelCopy() 1469 0.004223 0.004373 0.004749 0.03% AmrLevel::FillPatch() 187 0.002974 0.003097 0.003457 0.02% Castro::sum_integrated_quantities() 2 0.001553 0.002567 0.002888 0.02% MLCellLinOp::correctionResidual() 768 0.002276 0.002319 0.002432 0.02% MLMG::mgVcycle() 128 0.002143 0.002206 0.002308 0.02% MLMG:computeResOfCorrection() 768 0.001799 0.00184 0.001959 0.01% Castro::save_data_for_retry() 3 0.001276 0.0013 0.001322 0.01% StateDataPhysBCFunct::() 187 0.0009615 0.001141 0.001296 0.01% MLMG::actualBottomSolve() 128 0.000179 0.0002516 0.001249 0.01% Castro::subcycle_advance_ctu() 1 0.0002561 0.0003305 0.001137 0.01% FillPatchIterator::FillFromLevel0() 187 0.0009015 0.0009323 0.001021 0.01% main() 1 0.0006756 0.00076 0.0008488 0.01% Gravity::solve_for_phi() 17 0.0006434 0.0006723 0.0007632 0.01% Castro::Castro() 1 0.0004415 0.000519 0.0007276 0.00% FillPatchSingleLevel 187 0.0006529 0.0006728 0.0007067 0.00% MLMG::mgVcycle_down::2 128 0.0005919 0.0006315 0.0007007 0.00% MLMG::mgVcycle_down::0 128 0.0006335 0.0006476 0.0006834 0.00% MLMG::mgVcycle_down::3 128 0.0005629 0.0005985 0.0006716 0.00% MLMG::mgVcycle_down::1 128 0.0006167 0.0006311 0.0006703 0.00% MLMG::mgVcycle_down::5 128 0.0005847 0.0006018 0.0006554 0.00% MLMG::mgVcycle_down::4 128 0.0005602 0.0005745 0.0006281 0.00% AmrMesh::MakeNewGrids() 2 0.0005463 0.0005715 0.0006089 0.00% Castro::clean_state() 79 0.000465 0.0004898 0.0005572 0.00% MLCellLinOp::solutionResidual() 147 0.0005054 0.0005173 0.0005556 0.00% MLMG::oneIter() 128 0.0004452 0.0004586 0.0005019 0.00% Castro::do_old_sources() 40 0.000447 0.0004635 0.0004993 0.00% MLMG::mgVcycle_up::1 128 0.0004091 0.0004224 0.0004577 0.00% MLMG::mgVcycle_up::2 128 0.0003909 0.0004081 0.0004531 0.00% MLMG::mgVcycle_up::0 128 0.0004139 0.0004232 0.0004513 0.00% MLMG::mgVcycle_up::4 128 0.0003853 0.0003952 0.0004316 0.00% MLMG::mgVcycle_up::3 128 0.0003809 0.0003917 0.0004224 0.00% MLMG::mgVcycle_up::5 128 0.0003748 0.0003818 0.0003991 0.00% Castro::construct_new_gravity() 17 0.0001689 0.0001869 0.0003489 0.00% MLMG::computeResidual() 128 0.0002923 0.0002975 0.0003041 0.00% MLMG::mgVcycle_bottom 128 0.000274 0.0002812 0.0002975 0.00% Castro::problem_post_timestep() 2 1.158e-05 2.784e-05 0.0002588 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002064 0.0002219 0.0002519 0.00% Castro::do_new_sources() 17 0.0002133 0.0002243 0.0002447 0.00% Amr::coarseTimeStep() 1 0.0001423 0.0001709 0.0001797 0.00% Amr::defBaseLevel() 1 0.0001298 0.0001448 0.0001717 0.00% MLPoisson::define() 19 0.0001275 0.0001347 0.0001516 0.00% Castro::check_for_nan() 37 0.0001166 0.0001255 0.0001424 0.00% Gravity::test_composite_phi() 1 0.0001042 0.0001107 0.0001404 0.00% Castro::construct_old_gravity() 20 0.0001115 0.0001255 0.0001359 0.00% Castro::initialize_advance() 1 0.0001093 0.0001211 0.0001346 0.00% Castro::retry_advance_ctu() 11 7.079e-05 8.001e-05 0.0001297 0.00% TagBoxArray::mapPRD 2 0.0001216 0.0001248 0.0001295 0.00% Castro::apply_source_to_state() 37 0.0001156 0.0001209 0.0001287 0.00% MLLinOp::define() 19 7.845e-05 8.336e-05 9.203e-05 0.00% Amr::FinalizeInit() 1 1.192e-05 2.066e-05 9.199e-05 0.00% Amr::grid_places() 2 5.26e-05 5.849e-05 8.788e-05 0.00% Castro::create_source_corrector() 17 5.1e-05 5.633e-05 7.172e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.77e-05 5.93e-05 6.502e-05 0.00% Castro::expand_state() 20 5.51e-05 5.848e-05 6.405e-05 0.00% MLMG::computeMLResidual() 18 5.093e-05 5.438e-05 5.766e-05 0.00% Castro::derive() 18 4.67e-05 4.84e-05 5.223e-05 0.00% Amr::regrid() 1 1.619e-05 2.056e-05 4.655e-05 0.00% Gravity::actual_multilevel_solve() 1 3.031e-05 3.155e-05 3.54e-05 0.00% Amr::InitializeInit() 1 9.953e-06 1.217e-05 2.823e-05 0.00% Amr::timeStep() 1 1.566e-05 1.829e-05 2.511e-05 0.00% Castro::advance() 1 8.177e-06 9.994e-06 1.41e-05 0.00% Amr::bldFineLevels() 1 7.45e-06 8.449e-06 1.352e-05 0.00% Castro::post_init() 1 1.139e-05 1.227e-05 1.318e-05 0.00% Castro::post_regrid() 2 1.084e-05 1.157e-05 1.284e-05 0.00% Amr::init() 1 6.386e-06 7.106e-06 8.048e-06 0.00% Castro::post_timestep() 1 6.029e-06 6.681e-06 7.671e-06 0.00% Amr::initialInit() 1 3.197e-06 3.44e-06 3.843e-06 0.00% Other 2549 0.04785 0.06103 0.09807 0.65% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 15.04 15.05 15.05 100.00% Amr::coarseTimeStep() 1 12.76 12.76 12.76 84.80% Amr::timeStep() 1 12.71 12.71 12.71 84.44% Castro::advance() 1 12.64 12.64 12.65 84.02% Castro::subcycle_advance_ctu() 1 12.58 12.59 12.59 83.67% Castro::do_advance_ctu() 20 12.54 12.55 12.55 83.41% Castro::construct_ctu_hydro_source() 20 6.6 6.637 6.686 44.42% Castro::clean_state() 79 3.757 3.867 4.017 26.69% Castro::computeTemp() 80 2.39 2.497 2.649 17.60% Gravity::solve_phi_with_mlmg() 19 1.609 1.661 1.684 11.19% Castro::construct_new_gravity() 17 1.528 1.58 1.603 10.65% Castro::reset_internal_energy(MultiFab) 80 1.554 1.562 1.59 10.57% Gravity::solve_for_phi() 17 1.49 1.542 1.565 10.40% Castro::initialize_do_advance() 20 1.405 1.438 1.472 9.78% Gravity::actual_solve_with_mlmg() 19 1.401 1.405 1.409 9.36% Amr::InitAmr() 1 1.368 1.37 1.376 9.14% Castro::do_old_sources() 40 1.308 1.337 1.366 9.08% MLMG::solve() 18 1.347 1.352 1.354 9.00% MLMG::oneIter() 128 1.213 1.216 1.22 8.11% MLMG::mgVcycle() 128 1.198 1.201 1.206 8.01% Castro::evaluate_source_change() 196 1.017 1.021 1.038 6.89% Castro::do_new_sources() 17 0.9428 0.99 1.032 6.86% FabArray::ParallelCopy() 1469 0.4715 0.8123 1.008 6.70% Amr::init() 1 0.9111 0.9111 0.9111 6.05% amrex::Copy() 766 0.8713 0.8757 0.898 5.97% Castro::enforce_min_density() 79 0.8426 0.846 0.8637 5.74% Amr::initialInit() 1 0.8545 0.8545 0.8546 5.68% MLCellLinOp::smooth() 1544 0.6234 0.6795 0.8126 5.40% FabArray::ParallelCopy_finish() 1469 0.1941 0.5322 0.7306 4.85% AmrLevel::FillPatch() 187 0.5893 0.662 0.727 4.83% Amr::InitializeInit() 1 0.4455 0.4762 0.6187 4.11% Amr::defBaseLevel() 1 0.4455 0.4762 0.6186 4.11% MultiFab::AddProduct() 3598 0.5938 0.596 0.6045 4.02% Castro::normalize_species() 79 0.5499 0.5521 0.563 3.74% Castro::reset_internal_energy(Fab) 160 0.5128 0.5163 0.5252 3.49% FillPatchIterator::Initialize 187 0.3734 0.4461 0.5118 3.40% FillPatchIterator::FillFromLevel0() 187 0.3415 0.422 0.4801 3.19% FillPatchSingleLevel 187 0.3406 0.4211 0.4791 3.18% MLMG::mgVcycle_up::1 128 0.105 0.2806 0.4559 3.03% MLPoisson::Fsmooth() 6144 0.3499 0.3793 0.4551 3.02% MLCellLinOp::applyBC() 7148 0.3113 0.3591 0.4389 2.92% Amr::FinalizeInit() 1 0.2358 0.3782 0.4089 2.72% MLMG::addInterpCorrection() 768 0.06733 0.3236 0.3975 2.64% Castro::Castro() 1 0.364 0.3683 0.3909 2.60% FabArray::setVal() 5232 0.3735 0.3764 0.3852 2.56% Castro::initMFs() 1 0.3582 0.3624 0.3837 2.55% MultiFab::Subtract() 159 0.3496 0.352 0.3574 2.37% FabArray::FillBoundary() 6380 0.2385 0.2763 0.3192 2.12% MLMG::mgVcycle_up::2 128 0.01517 0.1276 0.3012 2.00% FabArray::ParallelCopy_nowait() 1469 0.2545 0.2757 0.2987 1.98% MLMG::mgVcycle_down::0 128 0.276 0.2791 0.281 1.87% Gravity::fill_multipole_BCs() 19 0.2055 0.2544 0.2758 1.83% Castro::expand_state() 20 0.2024 0.2277 0.2549 1.69% Castro::initData() 1 0.07968 0.1062 0.248 1.65% Amr::bldFineLevels() 1 0.04781 0.1897 0.22 1.46% MLCellLinOp::apply() 986 0.1571 0.1707 0.2137 1.42% MLMG::mgVcycle_up::0 128 0.2079 0.2101 0.2132 1.42% FabArray::sum() 3598 0.2052 0.2061 0.2095 1.39% MLMG::mgVcycle_up::3 128 0.01514 0.02917 0.1979 1.31% Amr::grid_places() 2 0.002733 0.1451 0.1759 1.17% AmrMesh::MakeNewGrids() 2 0.002679 0.1451 0.1758 1.17% TagBoxArray::mapPRD 2 0.0008545 0.1116 0.174 1.16% Castro::construct_old_source() 240 0.1652 0.1658 0.1686 1.12% FillBoundary_nowait() 6380 0.1002 0.1307 0.165 1.10% TagBoxArray::collate() 2 0.0003063 0.03216 0.1649 1.10% check_for_negative_density() 20 0.05703 0.1181 0.1649 1.10% FillBoundary_finish() 6380 0.09338 0.1285 0.1575 1.05% MLMG:computeResOfCorrection() 768 0.1189 0.1304 0.1554 1.03% MLCellLinOp::correctionResidual() 768 0.1171 0.1285 0.1535 1.02% Castro::post_init() 1 0.1415 0.1441 0.145 0.96% Castro::construct_new_source() 102 0.1283 0.1294 0.1317 0.88% MLMG::mgVcycle_down::1 128 0.1228 0.1264 0.1304 0.87% MLPoisson::Fapply() 986 0.09589 0.1 0.1164 0.77% Gravity::multilevel_solve_for_new_phi() 1 0.108 0.108 0.1081 0.72% Gravity::actual_multilevel_solve() 1 0.1079 0.108 0.108 0.72% Amr::writePlotFile() 2 0.1018 0.1019 0.1019 0.68% MLCellLinOp::solutionResidual() 147 0.08989 0.09554 0.1004 0.67% Castro::construct_old_gravity_source() 40 0.09891 0.09919 0.1002 0.67% Castro::post_regrid() 2 0.08806 0.08891 0.09155 0.61% FabArray::Saxpy() 111 0.08521 0.08564 0.0871 0.58% MLMG::mgVcycle_bottom 128 0.000454 0.005472 0.08057 0.54% MLMG::actualBottomSolve() 128 0.000179 0.005191 0.08027 0.53% MLMG::mgVcycle_down::2 128 0.01779 0.04746 0.07939 0.53% MLMG::computeResidual() 128 0.0763 0.07834 0.07896 0.52% MLCGSolver::bicgstab 8 0 0.004893 0.07829 0.52% Castro::construct_new_gravity_source() 17 0.07321 0.074 0.07569 0.50% Castro::derive() 18 0.07123 0.07188 0.07517 0.50% AmrLevel::derive() 18 0.07118 0.07184 0.07512 0.50% Castro::post_timestep() 1 0.06393 0.07225 0.07444 0.49% Castro::create_source_corrector() 17 0.05171 0.06011 0.07293 0.48% MLMG::mgVcycle_down::3 128 0.01713 0.02362 0.06745 0.45% Castro::construct_old_gravity() 20 0.04272 0.05112 0.05876 0.39% Gravity::get_old_grav_vector() 20 0.0426 0.05099 0.05862 0.39% amrex::average_down 768 0.03447 0.04107 0.05862 0.39% FabArray::Xpay() 915 0.05051 0.05245 0.05793 0.38% Castro::initialize_advance() 1 0.05023 0.05071 0.05201 0.35% Castro::check_for_nan() 37 0.04768 0.04793 0.04844 0.32% MultiFab::contains_nan() 37 0.04757 0.0478 0.04831 0.32% Amr::regrid() 1 0.04529 0.04578 0.04744 0.32% MLMG::mgVcycle_down::4 128 0.01707 0.01886 0.04116 0.27% MLPoisson::define() 19 0.03022 0.03536 0.04102 0.27% StateDataPhysBCFunct::() 187 0.0009615 0.02095 0.03934 0.26% Castro::retry_advance_ctu() 11 0.037 0.03752 0.03922 0.26% Gravity::get_new_grav_vector() 18 0.03675 0.03762 0.0392 0.26% StateData::FillBoundary(geom) 193 0 0.01981 0.03809 0.25% MLMG::mgVcycle_up::4 128 0.01516 0.01656 0.03295 0.22% MLMG::mgVcycle_down::5 128 0.01703 0.01811 0.02961 0.20% Castro::apply_source_to_state() 37 0.0285 0.02867 0.02922 0.19% VisMF::Write(FabArray) 4 0.01404 0.01519 0.02702 0.18% MLCellLinOp::defineAuxData() 19 0.0199 0.02162 0.02663 0.18% MLMG::mgVcycle_up::5 128 0.01527 0.01614 0.02518 0.17% MLMG::prepareForSolve() 18 0.02046 0.02143 0.02413 0.16% Castro::problem_post_timestep() 2 0.009506 0.02139 0.02402 0.16% Castro::wd_update() 2 0.009248 0.02137 0.024 0.16% Castro::sum_integrated_quantities() 2 0.02182 0.02196 0.02197 0.15% MLMG::computeMLResidual() 18 0.0109 0.01676 0.02132 0.14% Castro::save_data_for_retry() 3 0.01961 0.01988 0.02121 0.14% MLMG::ResNormInf() 146 0.01332 0.01561 0.01862 0.12% FabArrayBase::getFB() 6380 0.01344 0.01428 0.01794 0.12% FabArray::setDomainBndry() 187 0.0006436 0.008874 0.01675 0.11% FabArrayBase::getCPC() 2045 0.01073 0.01199 0.01547 0.10% FabArray::norminf() 244 0.01095 0.01138 0.01477 0.10% Gravity::test_composite_phi() 1 0.01459 0.01465 0.01469 0.10% amrex::Add() 128 0.0141 0.01421 0.01454 0.10% MLLinOp::define() 19 0.003556 0.007291 0.01089 0.07% Other 2549 0.08081 0.09988 0.1462 0.97% -------------------------------------------------------------------------------------------- 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 9 B 11 B 13 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 1601 KiB 1603 KiB 1603 KiB 2329 KiB 2329 KiB 2329 KiB Castro::construct_ctu_hydro_source() 6400 6400 403 KiB 405 KiB 411 KiB 1176 KiB 1176 KiB 1176 KiB StateData::define() 192 192 1057 KiB 1058 KiB 1058 KiB 1164 KiB 1164 KiB 1164 KiB Castro::initMFs() 256 256 1036 KiB 1037 KiB 1037 KiB 1142 KiB 1142 KiB 1142 KiB Castro::swap_state_time_levels() 160 160 680 KiB 681 KiB 681 KiB 803 KiB 803 KiB 803 KiB Amr::writePlotFile() 64 64 4311 B 4378 B 4967 B 752 KiB 752 KiB 752 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 53 KiB 53 KiB 54 KiB 672 KiB 672 KiB 672 KiB FillPatchIterator::Initialize 5984 5984 12 KiB 13 KiB 14 KiB 525 KiB 525 KiB 525 KiB Castro::initialize_do_advance() 640 640 438 KiB 438 KiB 439 KiB 525 KiB 525 KiB 525 KiB FabArray::ParallelCopy_nowait() 13984 13984 2155 B 3871 B 4976 B 225 KiB 299 KiB 378 KiB AmrLevel::derive() 1152 1152 1365 B 1379 B 1447 B 352 KiB 352 KiB 352 KiB Castro::enforce_min_density() 2528 2528 18 KiB 18 KiB 19 KiB 336 KiB 336 KiB 336 KiB Castro::initialize_advance() 32 32 132 KiB 132 KiB 132 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 192 129 KiB 129 KiB 129 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 6102 9404 B 9949 B 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 1280 1280 491 B 493 B 506 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 413 B 415 B 420 B 112 KiB 112 KiB 112 KiB Castro::do_advance_ctu() 960 960 2045 B 2055 B 2075 B 89 KiB 89 KiB 89 KiB Amr::writeSmallPlotFile() 64 64 18 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 6 B 7 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 110 B 112 B 117 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 129 B 154 B 178 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 64 31 KiB 32 KiB 32 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 2432 220 B 228 B 243 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1102 B 1106 B 1124 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 117 B 117 B 117 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1621 B 1677 B 1703 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 425 B 474 B 658 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 442 B 443 B 444 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 23 B 31 B 42 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 15 B 24 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.11-46-gce6502b2c0cb) finalized