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 18:13:18 UTC on 2025-08-03. Successfully read inputs file ... Castro git describe: 25.07-83-g9c4cd7c4c AMReX git describe: 25.06-22-g1d307c2c6-dirty Microphysics git describe: joss_submit-1-ga350de6c 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.162266537 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001209893 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010826405 Gravity::multilevel_solve_for_new_phi() time = 0.109371548 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010482017 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.054726632 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004110888 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 = -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 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 Castro::construct_old_gravity() time = 0.016588666 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.074714839 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 = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Castro::construct_old_gravity() time = 0.00801564 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.064262959 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, 177, 0 Invalid density = -0.1060744788 at index 20, 78, 0 Castro::construct_old_gravity() time = 0.007429221 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.063229108 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.004237687 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.059928436 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.33878677 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.014702993 Gravity::solve_for_phi() time = 0.091186147 on level 0 Castro::construct_new_gravity() time = 0.093616719 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.060059849 on level 0 Castro::construct_old_gravity() time = 0.003486416 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.060413644 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.342108185 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.014384163 Gravity::solve_for_phi() time = 0.091125268 on level 0 Castro::construct_new_gravity() time = 0.093449534 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.060020354 on level 0 Castro::construct_old_gravity() time = 0.006892685 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.061962309 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.342105362 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.015282271 Gravity::solve_for_phi() time = 0.092823271 on level 0 Castro::construct_new_gravity() time = 0.095160885 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.059452035 on level 0 Castro::construct_old_gravity() time = 0.004329959 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.059808643 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.348062942 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.014665291 Gravity::solve_for_phi() time = 0.091535814 on level 0 Castro::construct_new_gravity() time = 0.093856204 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.059838108 on level 0 Castro::construct_old_gravity() time = 0.002876955 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.059607807 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.360666586 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.014981123 Gravity::solve_for_phi() time = 0.092332816 on level 0 Castro::construct_new_gravity() time = 0.0946963 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.059361036 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.006239573 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.063283208 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.385533066 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.018594419 Gravity::solve_for_phi() time = 0.100823314 on level 0 Castro::construct_new_gravity() time = 0.103197362 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.061794966 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.003723079 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.061487038 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.362725845 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.016165314 Gravity::solve_for_phi() time = 0.093163628 on level 0 Castro::construct_new_gravity() time = 0.09551978 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.061122119 on level 0 Castro::construct_old_gravity() time = 0.003388325 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.061313237 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.362489237 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.016303803 Gravity::solve_for_phi() time = 0.095111326 on level 0 Castro::construct_new_gravity() time = 0.097410886 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.061143073 on level 0 Castro::construct_old_gravity() time = 0.003497776 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.061917525 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.360173448 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.016384349 Gravity::solve_for_phi() time = 0.096013154 on level 0 Castro::construct_new_gravity() time = 0.098411577 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.061123536 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.0037717 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.061360279 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.360925399 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.016603693 Gravity::solve_for_phi() time = 0.093618301 on level 0 Castro::construct_new_gravity() time = 0.095936685 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.061205872 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.003105725 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.061314088 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.361230568 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.016214548 Gravity::solve_for_phi() time = 0.093091939 on level 0 Castro::construct_new_gravity() time = 0.095433494 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.061167175 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.005112656 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.061608359 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.373982126 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.017131771 Gravity::solve_for_phi() time = 0.095222194 on level 0 Castro::construct_new_gravity() time = 0.097570488 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.061248417 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.003012297 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.061275409 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.374367603 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.025019831 Gravity::solve_for_phi() time = 0.105208482 on level 0 Castro::construct_new_gravity() time = 0.107526395 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.060968238 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.003449693 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.062456512 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.356870885 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.015367593 Gravity::solve_for_phi() time = 0.09278577 on level 0 Castro::construct_new_gravity() time = 0.095102388 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.059605752 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.003037013 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.059909129 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.355240403 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.015282793 Gravity::solve_for_phi() time = 0.090950154 on level 0 Castro::construct_new_gravity() time = 0.093280487 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.05983899 on level 0 Castro::construct_old_gravity() time = 0.002983284 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.060233903 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.348175875 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.016408153 Gravity::solve_for_phi() time = 0.092500612 on level 0 Castro::construct_new_gravity() time = 0.09494081 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.059761798 on level 0 Castro::construct_old_gravity() time = 0.00296809 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.059685994 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.35459435 on level 0 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.002613597129 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.015607409 Gravity::solve_for_phi() time = 0.091688959 on level 0 Castro::construct_new_gravity() time = 0.094087567 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.060586098 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.59945629 [STEP 1] FAB kilobyte spread across MPI nodes: [7120 ... 7120] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.048538589 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004040134 seconds Ending run at 18:13:33 UTC on 2025-08-03. Run time = 15.28193033 Run time without initialization = 12.65295496 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.2713 ... 15.2763 ... 15.2820 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.7955 5.8795 6.0260 39.43% Amr::InitAmr() 1 1.696 1.701 1.706 11.16% Castro::computeTemp() 80 0.8668 0.9571 1.131 7.40% amrex::Copy() 724 0.7904 0.8115 0.9676 6.33% FabArray::ParallelCopy_finish() 1429 0.229 0.5843 0.8085 5.29% Castro::normalize_species() 79 0.5442 0.5651 0.5955 3.90% MultiFab::AddProduct() 3439 0.5545 0.5617 0.5854 3.83% Castro::reset_internal_energy(Fab) 160 0.5162 0.5269 0.5432 3.55% MLPoisson::Fsmooth() 6144 0.3528 0.3831 0.4668 3.05% Castro::initMFs() 1 0.3578 0.3635 0.3971 2.60% MultiFab::Subtract() 159 0.3527 0.3569 0.3657 2.39% FabArray::setVal() 5013 0.3401 0.3466 0.3623 2.37% check_for_negative_density() 20 0.1183 0.2534 0.3369 2.20% Gravity::fill_multipole_BCs() 19 0.2083 0.2662 0.2965 1.94% FabArray::ParallelCopy_nowait() 1429 0.2153 0.2321 0.2569 1.68% FabArray::sum() 3439 0.1899 0.1925 0.2 1.31% Castro::initData() 1 0.003393 0.02788 0.1639 1.07% FillBoundary_nowait() 6380 0.09275 0.1206 0.157 1.03% FillBoundary_finish() 6380 0.1132 0.1321 0.1552 1.02% MLCellLinOp::applyBC() 7148 0.0633 0.08015 0.1278 0.84% MLPoisson::Fapply() 986 0.09463 0.09972 0.1159 0.76% TagBoxArray::collate() 2 0.0002211 0.01468 0.07622 0.50% Castro::construct_new_gravity_source() 17 0.07111 0.07272 0.07446 0.49% FabArray::Saxpy() 91 0.06932 0.0703 0.07234 0.47% FabArray::Xpay() 915 0.05164 0.05429 0.06025 0.39% Castro::construct_old_gravity_source() 20 0.0487 0.04922 0.0501 0.33% AmrLevel::derive() 18 0.04056 0.04149 0.04776 0.31% MultiFab::contains_nan() 37 0.044 0.04485 0.0466 0.30% MLMG::addInterpCorrection() 768 0.03342 0.03621 0.04455 0.29% StateData::FillBoundary(geom) 143 0 0.01941 0.03763 0.25% amrex::average_down 768 0.02524 0.02753 0.03451 0.23% VisMF::Write(FabArray) 4 0.01409 0.01538 0.03059 0.20% Castro::reset_internal_energy(MultiFab) 80 0.02316 0.02367 0.0243 0.16% Amr::writePlotFile() 2 0.006925 0.0203 0.02271 0.15% MLCellLinOp::defineAuxData() 19 0.01567 0.01668 0.02065 0.14% Castro::enforce_min_density() 79 0.01832 0.01866 0.019 0.12% FabArray::FillBoundary() 6380 0.01359 0.01427 0.01745 0.11% FabArray::setDomainBndry() 147 0.0004784 0.008585 0.01744 0.11% Castro::construct_new_source() 102 0.0161 0.01634 0.01698 0.11% MLCGSolver::bicgstab 8 0 0.00106 0.01696 0.11% Castro::sum_integrated_quantities() 2 0.002097 0.0148 0.01593 0.10% amrex::Add() 128 0.01442 0.01465 0.01574 0.10% Castro::evaluate_source_change() 196 0.01408 0.01447 0.01517 0.10% FabArray::norminf() 244 0.01078 0.01127 0.01436 0.09% Castro::wd_update() 2 0.005416 0.01087 0.01406 0.09% FabArrayBase::getFB() 6380 0.01066 0.01114 0.01372 0.09% FillPatchIterator::Initialize 147 0.01293 0.01315 0.01352 0.09% MLMG::prepareForSolve() 18 0.01019 0.01063 0.01186 0.08% MLCellLinOp::smooth() 1544 0.01028 0.01079 0.01178 0.08% Castro::construct_old_source() 120 0.0103 0.01047 0.01072 0.07% Gravity::actual_solve_with_mlmg() 19 0.007332 0.007783 0.009341 0.06% FabArrayBase::getCPC() 2005 0.00749 0.007744 0.008176 0.05% MLMG::ResNormInf() 146 0.004174 0.005692 0.008078 0.05% MLMG::solve() 18 0.001075 0.004204 0.00773 0.05% Gravity::get_old_grav_vector() 20 0.00645 0.006675 0.007154 0.05% Gravity::get_new_grav_vector() 18 0.006023 0.006188 0.006444 0.04% Castro::initialize_do_advance() 20 0.004367 0.004653 0.004908 0.03% MLCellLinOp::apply() 986 0.002051 0.002271 0.004838 0.03% FabArray::ParallelCopy() 1429 0.003566 0.003739 0.004118 0.03% AmrLevel::FillPatch() 147 0.002364 0.002491 0.002638 0.02% MLLinOp::defineGrids() 19 0.00195 0.002087 0.002372 0.02% MLCellLinOp::correctionResidual() 768 0.001919 0.001988 0.002097 0.01% MLMG::mgVcycle() 128 0.001684 0.001771 0.001936 0.01% MLMG:computeResOfCorrection() 768 0.001503 0.001581 0.001707 0.01% Castro::save_data_for_retry() 3 0.001339 0.001369 0.001554 0.01% MLMG::actualBottomSolve() 128 0.0001309 0.0002033 0.001237 0.01% Castro::subcycle_advance_ctu() 1 0.0002165 0.0002959 0.001039 0.01% StateDataPhysBCFunct::() 147 0.0007272 0.0008802 0.001027 0.01% main() 1 0.000641 0.0006966 0.0007635 0.00% FillPatchIterator::FillFromLevel0() 147 0.0006976 0.0007207 0.0007591 0.00% Gravity::solve_for_phi() 17 0.0006544 0.0006867 0.0007454 0.00% AmrMesh::MakeNewGrids() 2 0.0005483 0.0005781 0.0006958 0.00% MLMG::mgVcycle_down::3 128 0.0005152 0.0005539 0.0006387 0.00% MLMG::mgVcycle_down::2 128 0.000543 0.0005856 0.0006304 0.00% MLMG::mgVcycle_down::4 128 0.000495 0.0005185 0.0006253 0.00% MLMG::mgVcycle_down::0 128 0.0005644 0.00058 0.0006046 0.00% MLMG::mgVcycle_down::1 128 0.0005717 0.0005836 0.0006023 0.00% MLMG::mgVcycle_down::5 128 0.0005157 0.0005369 0.0005976 0.00% Castro::Castro() 1 0.0003962 0.000432 0.0005339 0.00% FillPatchSingleLevel 147 0.0004614 0.0004767 0.0005076 0.00% Castro::clean_state() 79 0.0004373 0.0004641 0.0004972 0.00% MLCellLinOp::solutionResidual() 147 0.0004569 0.0004654 0.0004798 0.00% Castro::do_advance_ctu() 20 0.0003936 0.0004144 0.0004745 0.00% MLMG::oneIter() 128 0.000395 0.0004093 0.0004246 0.00% MLMG::mgVcycle_up::0 128 0.0003793 0.0003908 0.0004037 0.00% MLMG::mgVcycle_up::1 128 0.0003597 0.0003813 0.0004035 0.00% MLMG::mgVcycle_up::2 128 0.0003183 0.0003447 0.0003899 0.00% MLMG::mgVcycle_up::4 128 0.0003184 0.0003329 0.00038 0.00% MLMG::mgVcycle_up::3 128 0.0003148 0.0003266 0.0003743 0.00% MLMG::mgVcycle_up::5 128 0.000309 0.0003193 0.0003499 0.00% Castro::do_old_sources() 20 0.0002854 0.0002973 0.000309 0.00% Castro::construct_new_gravity() 17 0.0001559 0.0001723 0.0002891 0.00% Gravity::solve_phi_with_mlmg() 19 0.000221 0.000234 0.0002616 0.00% MLMG::mgVcycle_bottom 128 0.0002271 0.0002361 0.0002579 0.00% Castro::problem_post_timestep() 2 1.043e-05 2.807e-05 0.0002546 0.00% MLMG::computeResidual() 128 0.0002363 0.000246 0.0002531 0.00% Castro::do_new_sources() 17 0.0002134 0.0002258 0.0002384 0.00% Gravity::update_max_rhs() 2 0.0001886 0.0001962 0.0002047 0.00% Amr::defBaseLevel() 1 0.0001405 0.0001588 0.0001837 0.00% Amr::coarseTimeStep() 1 0.0001274 0.000154 0.0001605 0.00% Gravity::test_composite_phi() 1 0.0001118 0.0001164 0.0001437 0.00% MLPoisson::define() 19 0.0001205 0.0001275 0.0001401 0.00% Castro::retry_advance_ctu() 11 6.421e-05 7.699e-05 0.0001379 0.00% Castro::construct_old_gravity() 20 0.000114 0.0001224 0.0001376 0.00% TagBoxArray::mapPRD 2 0.0001195 0.0001275 0.0001321 0.00% Castro::check_for_nan() 37 0.0001076 0.0001184 0.0001317 0.00% Castro::initialize_advance() 1 0.0001111 0.0001146 0.0001202 0.00% Castro::apply_source_to_state() 37 0.0001011 0.000107 0.0001161 0.00% MLLinOp::define() 19 7.387e-05 7.868e-05 8.866e-05 0.00% Amr::grid_places() 2 5.227e-05 6.781e-05 8.681e-05 0.00% Amr::FinalizeInit() 1 8.546e-06 1.428e-05 8.474e-05 0.00% Castro::create_source_corrector() 17 5.832e-05 6.397e-05 7.181e-05 0.00% Amr::regrid() 1 1.641e-05 2.648e-05 6.419e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.862e-05 6.085e-05 6.337e-05 0.00% Castro::expand_state() 20 5.126e-05 5.397e-05 6.255e-05 0.00% Castro::derive() 18 4.082e-05 4.452e-05 5.595e-05 0.00% Amr::timeStep() 1 1.403e-05 1.906e-05 5.14e-05 0.00% MLMG::computeMLResidual() 18 4.252e-05 4.542e-05 5e-05 0.00% Gravity::actual_multilevel_solve() 1 3.119e-05 3.25e-05 3.367e-05 0.00% Castro::post_init() 1 1.205e-05 1.331e-05 1.904e-05 0.00% Amr::InitializeInit() 1 9.132e-06 1.279e-05 1.783e-05 0.00% Castro::advance() 1 7.088e-06 8.944e-06 1.61e-05 0.00% Castro::post_regrid() 2 1.091e-05 1.204e-05 1.412e-05 0.00% Amr::bldFineLevels() 1 6.995e-06 8.254e-06 1.231e-05 0.00% Amr::init() 1 6.301e-06 7.087e-06 7.903e-06 0.00% Castro::post_timestep() 1 5.914e-06 6.343e-06 6.812e-06 0.00% Amr::initialInit() 1 2.999e-06 3.286e-06 3.957e-06 0.00% Other 2499 0.04764 0.07104 0.1097 0.72% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 15.27 15.28 15.28 100.00% Amr::coarseTimeStep() 1 12.65 12.65 12.65 82.79% Amr::timeStep() 1 12.6 12.6 12.6 82.45% Castro::advance() 1 12.53 12.54 12.54 82.05% Castro::subcycle_advance_ctu() 1 12.48 12.49 12.49 81.72% Castro::do_advance_ctu() 20 12.44 12.45 12.45 81.47% Castro::construct_ctu_hydro_source() 20 6.88 6.92 6.971 45.62% Castro::clean_state() 79 3.699 3.827 4.013 26.26% Castro::computeTemp() 80 2.381 2.48 2.667 17.45% Gravity::solve_phi_with_mlmg() 19 1.625 1.687 1.718 11.24% Amr::InitAmr() 1 1.696 1.701 1.706 11.16% Castro::construct_new_gravity() 17 1.546 1.607 1.638 10.72% Castro::reset_internal_energy(MultiFab) 80 1.495 1.523 1.602 10.48% Gravity::solve_for_phi() 17 1.505 1.567 1.598 10.46% Castro::initialize_do_advance() 20 1.401 1.453 1.492 9.77% Gravity::actual_solve_with_mlmg() 19 1.411 1.417 1.422 9.30% MLMG::solve() 18 1.355 1.361 1.363 8.92% Castro::do_old_sources() 20 1.165 1.193 1.225 8.02% MLMG::oneIter() 128 1.221 1.223 1.225 8.02% MLMG::mgVcycle() 128 1.206 1.208 1.21 7.92% FabArray::ParallelCopy() 1429 0.4868 0.8407 1.069 6.99% Castro::do_new_sources() 17 0.9418 0.9774 1.028 6.73% Castro::evaluate_source_change() 196 0.9501 0.9628 1.003 6.56% amrex::Copy() 724 0.7904 0.8115 0.9676 6.33% Amr::init() 1 0.9226 0.9226 0.9226 6.04% Amr::initialInit() 1 0.8633 0.8633 0.8634 5.65% Castro::enforce_min_density() 79 0.7961 0.8089 0.8487 5.55% MLCellLinOp::smooth() 1544 0.6269 0.684 0.8283 5.42% FabArray::ParallelCopy_finish() 1429 0.229 0.5843 0.8085 5.29% AmrLevel::FillPatch() 147 0.5556 0.6505 0.7509 4.91% Amr::InitializeInit() 1 0.4447 0.477 0.6057 3.96% Amr::defBaseLevel() 1 0.4447 0.477 0.6057 3.96% Castro::normalize_species() 79 0.5442 0.5651 0.5955 3.90% MultiFab::AddProduct() 3439 0.5545 0.5617 0.5854 3.83% FillPatchIterator::Initialize 147 0.3676 0.4627 0.567 3.71% FillPatchIterator::FillFromLevel0() 147 0.338 0.4409 0.5479 3.59% FillPatchSingleLevel 147 0.3372 0.4402 0.5472 3.58% Castro::reset_internal_energy(Fab) 160 0.5162 0.5269 0.5432 3.55% MLMG::mgVcycle_up::1 128 0.1078 0.2881 0.4687 3.07% MLPoisson::Fsmooth() 6144 0.3528 0.3831 0.4668 3.05% MLCellLinOp::applyBC() 7148 0.3118 0.359 0.4365 2.86% Amr::FinalizeInit() 1 0.2576 0.3863 0.4186 2.74% MLMG::addInterpCorrection() 768 0.06768 0.3287 0.406 2.66% Castro::Castro() 1 0.3634 0.369 0.4034 2.64% Castro::initMFs() 1 0.3579 0.3636 0.3972 2.60% MultiFab::Subtract() 159 0.3527 0.3569 0.3657 2.39% FabArray::setVal() 5013 0.3401 0.3466 0.3623 2.37% check_for_negative_density() 20 0.1183 0.2534 0.3369 2.20% FabArray::FillBoundary() 6380 0.2345 0.2789 0.3157 2.07% MLMG::mgVcycle_up::2 128 0.01373 0.1284 0.3044 1.99% Gravity::fill_multipole_BCs() 19 0.2104 0.2684 0.2986 1.95% MLMG::mgVcycle_down::0 128 0.2773 0.2795 0.2834 1.85% FabArray::ParallelCopy_nowait() 1429 0.2348 0.2526 0.2782 1.82% Castro::expand_state() 20 0.1936 0.2289 0.2681 1.75% Castro::initData() 1 0.07796 0.1062 0.2398 1.57% MLMG::mgVcycle_up::0 128 0.209 0.2101 0.2112 1.38% Amr::bldFineLevels() 1 0.04547 0.1757 0.2077 1.36% MLCellLinOp::apply() 986 0.1566 0.1688 0.2057 1.35% FabArray::sum() 3439 0.1899 0.1925 0.2 1.31% MLMG::mgVcycle_up::3 128 0.01365 0.02765 0.1936 1.27% FillBoundary_nowait() 6380 0.104 0.1325 0.1696 1.11% Castro::post_init() 1 0.1563 0.1672 0.1693 1.11% Amr::grid_places() 2 0.00251 0.1312 0.1635 1.07% AmrMesh::MakeNewGrids() 2 0.002449 0.1311 0.1634 1.07% TagBoxArray::mapPRD 2 0.0008338 0.1151 0.1616 1.06% FillBoundary_finish() 6380 0.1132 0.1321 0.1552 1.02% MLMG:computeResOfCorrection() 768 0.1183 0.1288 0.1535 1.00% MLCellLinOp::correctionResidual() 768 0.1168 0.1272 0.1518 0.99% MLMG::mgVcycle_down::1 128 0.1304 0.1346 0.1397 0.91% Castro::construct_new_source() 102 0.1261 0.1281 0.1313 0.86% Castro::create_source_corrector() 17 0.05953 0.08425 0.1225 0.80% MLPoisson::Fapply() 986 0.09463 0.09972 0.1159 0.76% Gravity::multilevel_solve_for_new_phi() 1 0.1093 0.1093 0.1094 0.72% Gravity::actual_multilevel_solve() 1 0.1092 0.1092 0.1093 0.72% Amr::writePlotFile() 2 0.1036 0.1037 0.1037 0.68% MLCellLinOp::solutionResidual() 147 0.08904 0.09657 0.1037 0.68% Castro::post_regrid() 2 0.08581 0.08783 0.09871 0.65% Castro::derive() 18 0.06935 0.0714 0.08665 0.57% AmrLevel::derive() 18 0.06931 0.07136 0.08659 0.57% Castro::construct_old_source() 120 0.08184 0.08303 0.08519 0.56% Castro::construct_old_gravity() 20 0.05677 0.07293 0.0848 0.55% Gravity::get_old_grav_vector() 20 0.05664 0.0728 0.08468 0.55% MLMG::mgVcycle_down::2 128 0.01593 0.04741 0.08148 0.53% MLMG::computeResidual() 128 0.07726 0.07831 0.07922 0.52% TagBoxArray::collate() 2 0.0002211 0.01468 0.07622 0.50% Castro::construct_new_gravity_source() 17 0.07111 0.07272 0.07446 0.49% MLMG::mgVcycle_bottom 128 0.000358 0.004935 0.07343 0.48% MLMG::actualBottomSolve() 128 0.0001309 0.004699 0.07317 0.48% FabArray::Saxpy() 91 0.06932 0.0703 0.07234 0.47% MLCGSolver::bicgstab 8 0 0.004455 0.07128 0.47% MLMG::mgVcycle_down::3 128 0.01534 0.02226 0.06931 0.45% Castro::post_timestep() 1 0.06057 0.06399 0.06661 0.44% amrex::average_down 768 0.03421 0.04151 0.06058 0.40% FabArray::Xpay() 915 0.05164 0.05429 0.06025 0.39% Amr::regrid() 1 0.04408 0.04567 0.05558 0.36% Castro::initialize_advance() 1 0.04862 0.04965 0.05249 0.34% Castro::construct_old_gravity_source() 20 0.0487 0.04922 0.0501 0.33% Castro::check_for_nan() 37 0.04411 0.04497 0.04672 0.31% MultiFab::contains_nan() 37 0.044 0.04485 0.0466 0.30% MLPoisson::define() 19 0.03289 0.03813 0.04451 0.29% Gravity::get_new_grav_vector() 18 0.03938 0.04047 0.04294 0.28% Castro::retry_advance_ctu() 11 0.03592 0.03715 0.04258 0.28% MLMG::mgVcycle_down::4 128 0.01526 0.01716 0.04067 0.27% StateDataPhysBCFunct::() 147 0.0007272 0.02029 0.03862 0.25% StateData::FillBoundary(geom) 143 0 0.01941 0.03763 0.25% Castro::sum_integrated_quantities() 2 0.03502 0.03517 0.0352 0.23% MLMG::mgVcycle_up::4 128 0.01375 0.01513 0.03215 0.21% VisMF::Write(FabArray) 4 0.01409 0.01538 0.03059 0.20% Castro::apply_source_to_state() 37 0.02836 0.02868 0.02912 0.19% MLMG::mgVcycle_down::5 128 0.01522 0.01629 0.02834 0.19% MLCellLinOp::defineAuxData() 19 0.0198 0.02167 0.02787 0.18% Castro::save_data_for_retry() 3 0.01933 0.01988 0.02522 0.17% MLMG::mgVcycle_up::5 128 0.01378 0.01469 0.02465 0.16% MLMG::prepareForSolve() 18 0.02029 0.02147 0.02449 0.16% MLMG::computeMLResidual() 18 0.01102 0.01773 0.02385 0.16% Castro::problem_post_timestep() 2 0.01043 0.01473 0.01754 0.11% MLMG::ResNormInf() 146 0.01389 0.01529 0.01753 0.11% Castro::wd_update() 2 0.01042 0.0147 0.01753 0.11% FabArray::setDomainBndry() 147 0.0004784 0.008585 0.01744 0.11% amrex::Add() 128 0.01442 0.01465 0.01574 0.10% FabArrayBase::getFB() 6380 0.01097 0.01187 0.01539 0.10% FabArrayBase::getCPC() 2005 0.009784 0.01123 0.01526 0.10% Gravity::update_max_rhs() 2 0.001604 0.01247 0.01448 0.09% FabArray::norminf() 244 0.01078 0.01127 0.01436 0.09% Gravity::test_composite_phi() 1 0.01424 0.01429 0.01431 0.09% MLLinOp::define() 19 0.003923 0.009841 0.01359 0.09% MLLinOp::defineGrids() 19 0.003841 0.009762 0.01351 0.09% Other 2499 0.07748 0.1007 0.1443 0.94% -------------------------------------------------------------------------------------------- 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 12 B 15 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 1574 KiB 1574 KiB 1575 KiB 2346 KiB 2346 KiB 2346 KiB StateData::define() 192 192 1041 KiB 1042 KiB 1042 KiB 1173 KiB 1173 KiB 1173 KiB Castro::construct_ctu_hydro_source() 6400 6400 400 KiB 406 KiB 416 KiB 1164 KiB 1164 KiB 1164 KiB Castro::initMFs() 256 256 970 KiB 971 KiB 971 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 672 KiB 672 KiB 672 KiB 812 KiB 812 KiB 812 KiB Amr::writePlotFile() 64 64 4057 B 4180 B 4885 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 47 KiB 48 KiB 50 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 11 KiB 13 KiB 15 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 408 KiB 408 KiB 408 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 2056 B 4151 B 5683 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1260 B 1288 B 1461 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 16 KiB 16 KiB 17 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 143 KiB 143 KiB 143 KiB 175 KiB 175 KiB 175 KiB Castro::buildMetrics() 192 192 126 KiB 126 KiB 126 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 6102 9322 B 9867 B 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 241 B 246 B 265 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 405 B 409 B 421 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 17 B 19 B 35 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 32 KiB 32 KiB 32 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 49 B 49 B 49 B 52 KiB 52 KiB 52 KiB Gravity::get_new_grav_vector() 576 576 129 B 132 B 141 B 50 KiB 50 KiB 50 KiB Gravity::get_old_grav_vector() 640 640 187 B 241 B 281 B 50 KiB 50 KiB 50 KiB Gravity::update_max_rhs() 192 192 3 B 14 B 17 B 49 KiB 49 KiB 49 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 223 B 229 B 236 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1013 B 1027 B 1069 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 1613 B 1679 B 1713 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 421 B 470 B 653 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 438 B 439 B 440 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 21 B 32 B 43 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 19 B 27 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 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 4 B 6 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