Initializing AMReX (25.03-15-g546cc5e16431-dirty)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (25.03-15-g546cc5e16431-dirty) initialized Starting run at 06:56:49 UTC on 2025-06-09. Successfully read inputs file ... Castro git describe: 25.06 AMReX git describe: 25.03-15-g546cc5e16-dirty Microphysics git describe: 25.06-7-gd774002a 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.160851357 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001217183 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010858249 Gravity::multilevel_solve_for_new_phi() time = 0.112570714 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.011443447 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.049800979 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004104394 seconds [Level 0 step 1] ADVANCE at time 0 with dt = 0.25 [1m[32m Beginning subcycle 1 starting at time 0 with dt = 0.25[0m Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.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 = -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 = -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.239071148Invalid density = -108.2035502 at index 19, 177, 0 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 = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477044 at index 18, 179, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.769409 at index Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 17, 181, 0 Invalid density = -415.2911996Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 -415.2911996 at index 16, 72, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, Invalid density = -552.769409 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -80.16477044 at index 18, 76, 0 Invalid density = 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-61.58489855 at index 19, 77, 0 Invalid density = -108.2035502 at index 19, 78, 0 , 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 = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 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 Invalid density = -0.01295059363 at index 21, 83, 0 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.004989057 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.063221144 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry. [1m[32m Subcycle completed[0m [1m[31m Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125 [0m [1m[32m Beginning subcycle 1 starting at time 0 with dt = 0.125[0m Estimated number of subcycles remaining: 2 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Castro::construct_old_gravity() time = 0.005427119 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.063579998 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry. [1m[32m Subcycle completed[0m [1m[31m Timestep 0.125 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.0625 [0m [1m[32m Beginning subcycle 1 starting at time 0 with dt = 0.0625[0m Estimated number of subcycles remaining: 4 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.1060744788 at index 20, 78, 0 Invalid density = -0.1060744788 at index 20, 177, 0 Castro::construct_old_gravity() time = 0.005463029 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.060038877 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry. [1m[32m Subcycle completed[0m [1m[31m Timestep 0.0625 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.03125 [0m [1m[32m Beginning subcycle 1 starting at time 0 with dt = 0.03125[0m Estimated number of subcycles remaining: 8 Adding an SDC iteration due to the retry. Beginning SDC iteration 1 of 3. ... Entering construct_ctu_hydro_source() on level 0 Castro::construct_old_gravity() time = 0.005219567 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.059007036 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.332170361 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.013918508 Gravity::solve_for_phi() time = 0.090260533 on level 0 Castro::construct_new_gravity() time = 0.092490054 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.05885448 on level 0 Castro::construct_old_gravity() time = 0.003680338 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.058862792 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.334712282 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.014129363 Gravity::solve_for_phi() time = 0.089896105 on level 0 Castro::construct_new_gravity() time = 0.092061163 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.058512809 on level 0 Castro::construct_old_gravity() time = 0.003360799 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.058398164 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.341089106 on level 0 Ending SDC iteration 3 of 3. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 2 starting at time 0.03125 with dt = 0.03125[0m Estimated number of subcycles remaining: 7 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.013777353 Gravity::solve_for_phi() time = 0.094467306 on level 0 Castro::construct_new_gravity() time = 0.096635794 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.058878892 on level 0 Castro::construct_old_gravity() time = 0.004515643 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.058388507 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.342240247 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.014181882 Gravity::solve_for_phi() time = 0.089551891 on level 0 Castro::construct_new_gravity() time = 0.091820135 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.059074148 on level 0 Castro::construct_old_gravity() time = 0.00292411 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.059382684 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.350417012 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 3 starting at time 0.0625 with dt = 0.03125[0m Estimated number of subcycles remaining: 6 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.014857846 Gravity::solve_for_phi() time = 0.093677124 on level 0 Castro::construct_new_gravity() time = 0.095838809 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.059347408 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.00348755 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.058837727 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.344559483 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.014357612 Gravity::solve_for_phi() time = 0.090084932 on level 0 Castro::construct_new_gravity() time = 0.092229049 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.058898648 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.003826608 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.058299909 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.35002589 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 4 starting at time 0.09375 with dt = 0.03125[0m Estimated number of subcycles remaining: 5 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.014424883 Gravity::solve_for_phi() time = 0.092543316 on level 0 Castro::construct_new_gravity() time = 0.094730245 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.058760615 on level 0 Castro::construct_old_gravity() time = 0.004058879 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.058493642 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.342512765 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.014705199 Gravity::solve_for_phi() time = 0.09154106 on level 0 Castro::construct_new_gravity() time = 0.093679467 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.058651264 on level 0 Castro::construct_old_gravity() time = 0.00333529 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.058320127 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.351610862 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 5 starting at time 0.125 with dt = 0.03125[0m Estimated number of subcycles remaining: 4 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.014315053 Gravity::solve_for_phi() time = 0.093418437 on level 0 Castro::construct_new_gravity() time = 0.095583681 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.058483885 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.002811523 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.058237672 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.339475607 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.014892746 Gravity::solve_for_phi() time = 0.090707842 on level 0 Castro::construct_new_gravity() time = 0.09287438 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.058356726 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.002872887 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.058379913 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.343814957 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 6 starting at time 0.15625 with dt = 0.03125[0m Estimated number of subcycles remaining: 3 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.014681654 Gravity::solve_for_phi() time = 0.090458745 on level 0 Castro::construct_new_gravity() time = 0.092614425 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.06093994 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.002766991 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.058577186 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.339645922 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.014830327 Gravity::solve_for_phi() time = 0.090403936 on level 0 Castro::construct_new_gravity() time = 0.09254677 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.058422112 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.002791042 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.058162009 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.345414578 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 7 starting at time 0.1875 with dt = 0.03125[0m Estimated number of subcycles remaining: 2 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.014679478 Gravity::solve_for_phi() time = 0.090890514 on level 0 Castro::construct_new_gravity() time = 0.093082979 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.058759531 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.003194045 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.058496891 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.342311601 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.016455947 Gravity::solve_for_phi() time = 0.092519688 on level 0 Castro::construct_new_gravity() time = 0.09467393 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.058789589 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.002824812 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.058403515 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.34589702 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m [1m[32m Beginning subcycle 8 starting at time 0.21875 with dt = 0.03125[0m Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.014906553 Gravity::solve_for_phi() time = 0.09103874 on level 0 Castro::construct_new_gravity() time = 0.093183516 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.058563082 on level 0 Castro::construct_old_gravity() time = 0.0029206 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.05861011 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.340789165 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.014982668 Gravity::solve_for_phi() time = 0.112902833 on level 0 Castro::construct_new_gravity() time = 0.1150483 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.058883022 on level 0 Castro::construct_old_gravity() time = 0.002808595 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.05894766 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.347795672 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m Subcycling complete Zones advanced per microsecond at this level: 0.002704463489 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.015753265 Gravity::solve_for_phi() time = 0.092883508 on level 0 Castro::construct_new_gravity() time = 0.095085842 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.058646378 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.17644655 [STEP 1] FAB kilobyte spread across MPI nodes: [6998 ... 6998] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.04781568 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004110307 seconds Ending run at 06:57:04 UTC on 2025-06-09. Run time = 14.47325635 Run time without initialization = 12.22946041 Average number of zones advanced per microsecond: 0.003 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 14.4642 ... 14.4669 ... 14.4735 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.6503 5.6988 5.7911 40.01% Amr::InitAmr() 1 1.345 1.348 1.354 9.36% Castro::computeTemp() 80 0.8169 0.9236 1.079 7.45% amrex::Copy() 724 0.8029 0.8084 0.8218 5.68% FabArray::ParallelCopy_finish() 1429 0.2212 0.5288 0.7203 4.98% MultiFab::AddProduct() 3439 0.5462 0.5501 0.5673 3.92% Castro::normalize_species() 79 0.5429 0.5472 0.5602 3.87% Castro::reset_internal_energy(Fab) 160 0.5037 0.5085 0.5192 3.59% MLPoisson::Fsmooth() 6144 0.3472 0.3754 0.4506 3.11% Castro::initMFs() 1 0.35 0.3565 0.3714 2.57% FabArray::setVal() 5013 0.344 0.3486 0.3578 2.47% MultiFab::Subtract() 159 0.3345 0.3372 0.3514 2.43% Gravity::fill_multipole_BCs() 19 0.1964 0.2457 0.2668 1.84% check_for_negative_density() 20 0.07739 0.1764 0.2476 1.71% FabArray::ParallelCopy_nowait() 1429 0.2082 0.2259 0.246 1.70% FabArray::sum() 3439 0.1925 0.1948 0.2003 1.38% FillBoundary_finish() 6380 0.1254 0.1561 0.19 1.31% Castro::initData() 1 0.003325 0.02721 0.158 1.09% FillBoundary_nowait() 6380 0.09203 0.1206 0.1563 1.08% TagBoxArray::collate() 2 0.0002248 0.02662 0.1404 0.97% MLCellLinOp::applyBC() 7148 0.06221 0.07797 0.1245 0.86% MLPoisson::Fapply() 986 0.0966 0.1013 0.1158 0.80% FabArray::Saxpy() 264 0.07246 0.07317 0.07992 0.55% Castro::construct_new_gravity_source() 17 0.06832 0.0693 0.07318 0.51% FabArray::Xpay() 945 0.0518 0.05431 0.06092 0.42% Castro::construct_old_gravity_source() 20 0.04578 0.04608 0.04722 0.33% MultiFab::contains_nan() 37 0.04394 0.04441 0.04572 0.32% MLMG::addInterpCorrection() 768 0.03313 0.03584 0.04335 0.30% AmrLevel::derive() 18 0.04005 0.04038 0.04229 0.29% StateData::FillBoundary(geom) 143 0 0.01744 0.03378 0.23% amrex::average_down 768 0.02493 0.02715 0.03368 0.23% VisMF::Write(FabArray) 4 0.0136 0.01467 0.02604 0.18% Castro::reset_internal_energy(MultiFab) 80 0.02322 0.02359 0.02412 0.17% MLCellLinOp::defineAuxData() 19 0.01557 0.01647 0.02011 0.14% Castro::enforce_min_density() 79 0.01859 0.01875 0.01915 0.13% FabArray::FillBoundary() 6380 0.01371 0.01431 0.01735 0.12% Amr::writePlotFile() 2 0.006568 0.01592 0.0172 0.12% MLCGSolver::bicgstab 8 0 0.001067 0.01708 0.12% Castro::construct_new_source() 102 0.01596 0.01606 0.01654 0.11% amrex::Add() 128 0.01429 0.01444 0.01503 0.10% Castro::evaluate_source_change() 196 0.01386 0.01423 0.01498 0.10% FabArray::setDomainBndry() 147 0.0004527 0.007847 0.01486 0.10% FabArray::norminf() 244 0.01097 0.01134 0.01403 0.10% FillPatchIterator::Initialize 147 0.0128 0.01298 0.01346 0.09% FabArrayBase::getFB() 6380 0.01033 0.01067 0.01314 0.09% MLMG::ResNormInf() 146 0.005922 0.009173 0.0129 0.09% Castro::wd_update() 2 0.005213 0.008873 0.01185 0.08% MLMG::prepareForSolve() 18 0.01021 0.01052 0.01176 0.08% MLCellLinOp::smooth() 1544 0.01002 0.01057 0.01169 0.08% Castro::construct_old_source() 120 0.01015 0.01032 0.01108 0.08% MLMG::solve() 18 0.001462 0.004589 0.009018 0.06% Gravity::actual_solve_with_mlmg() 19 0.007245 0.007609 0.00896 0.06% FabArrayBase::getCPC() 2005 0.007387 0.007569 0.007892 0.05% Gravity::get_old_grav_vector() 20 0.006609 0.006793 0.007804 0.05% Gravity::get_new_grav_vector() 18 0.006179 0.006241 0.006401 0.04% Castro::initialize_do_advance() 20 0.00437 0.004608 0.005058 0.03% MLCellLinOp::apply() 986 0.002002 0.002225 0.004691 0.03% FabArray::ParallelCopy() 1429 0.0035 0.00366 0.004002 0.03% Castro::sum_integrated_quantities() 2 0.002222 0.003445 0.003741 0.03% AmrLevel::FillPatch() 147 0.002336 0.002439 0.002541 0.02% MLCellLinOp::correctionResidual() 768 0.001902 0.00194 0.001997 0.01% MLMG::mgVcycle() 128 0.00161 0.001691 0.001833 0.01% MLMG:computeResOfCorrection() 768 0.001442 0.001487 0.001561 0.01% Castro::save_data_for_retry() 3 0.001311 0.001343 0.001489 0.01% MLMG::actualBottomSolve() 128 0.0001301 0.0002011 0.0012 0.01% StateDataPhysBCFunct::() 147 0.0006931 0.0008538 0.0009728 0.01% main() 1 0.0007019 0.0007805 0.0009501 0.01% Castro::subcycle_advance_ctu() 1 0.0002126 0.0002895 0.0009325 0.01% Gravity::solve_for_phi() 17 0.000648 0.0006785 0.0007424 0.01% FillPatchIterator::FillFromLevel0() 147 0.0006738 0.0006985 0.0007197 0.00% MLMG::mgVcycle_down::2 128 0.0005316 0.0005713 0.0006078 0.00% MLMG::mgVcycle_down::3 128 0.0004996 0.0005363 0.0006022 0.00% AmrMesh::MakeNewGrids() 2 0.0005445 0.0005653 0.0005976 0.00% MLMG::mgVcycle_down::0 128 0.0005451 0.0005611 0.0005913 0.00% MLMG::mgVcycle_down::1 128 0.0005453 0.0005601 0.0005904 0.00% Castro::Castro() 1 0.0004738 0.0005213 0.0005883 0.00% MLMG::mgVcycle_down::5 128 0.0005065 0.0005227 0.0005814 0.00% MLMG::mgVcycle_down::4 128 0.0004839 0.0005053 0.0005808 0.00% FillPatchSingleLevel 147 0.0004624 0.0004779 0.0004932 0.00% Castro::clean_state() 79 0.0004308 0.0004533 0.0004789 0.00% MLCellLinOp::solutionResidual() 147 0.0004442 0.0004547 0.000469 0.00% Castro::do_advance_ctu() 20 0.0003905 0.000411 0.0004425 0.00% MLMG::oneIter() 128 0.0003821 0.0003924 0.0004082 0.00% MLMG::mgVcycle_up::1 128 0.0003625 0.000382 0.0003979 0.00% MLMG::mgVcycle_up::0 128 0.000369 0.00038 0.000393 0.00% Castro::do_old_sources() 20 0.0002782 0.0002988 0.0003875 0.00% MLMG::mgVcycle_up::2 128 0.0003152 0.0003423 0.0003776 0.00% MLMG::mgVcycle_up::3 128 0.0003109 0.000321 0.0003643 0.00% MLMG::mgVcycle_up::4 128 0.0003155 0.0003244 0.0003536 0.00% MLMG::mgVcycle_up::5 128 0.000308 0.0003162 0.0003397 0.00% Castro::construct_new_gravity() 17 0.0001552 0.0001716 0.0002788 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002169 0.0002311 0.0002649 0.00% Castro::problem_post_timestep() 2 1.066e-05 2.729e-05 0.0002466 0.00% MLMG::computeResidual() 128 0.0002277 0.0002368 0.0002456 0.00% MLMG::mgVcycle_bottom 128 0.0002189 0.0002275 0.000243 0.00% Castro::do_new_sources() 17 0.0002082 0.0002209 0.0002329 0.00% Amr::coarseTimeStep() 1 0.0001403 0.0001655 0.0001731 0.00% Amr::defBaseLevel() 1 0.0001258 0.0001475 0.0001721 0.00% TagBoxArray::mapPRD 2 0.0001281 0.0001335 0.0001399 0.00% MLPoisson::define() 19 0.000121 0.0001251 0.0001333 0.00% Castro::check_for_nan() 37 0.0001034 0.0001119 0.0001316 0.00% Castro::construct_old_gravity() 20 0.0001068 0.0001138 0.0001316 0.00% Castro::retry_advance_ctu() 11 6.013e-05 7.06e-05 0.0001292 0.00% Gravity::test_composite_phi() 1 0.000106 0.00011 0.0001249 0.00% Castro::apply_source_to_state() 37 0.0001 0.0001046 0.0001197 0.00% Castro::initialize_advance() 1 0.0001087 0.0001126 0.0001191 0.00% Amr::grid_places() 2 5.703e-05 6.538e-05 8.614e-05 0.00% MLLinOp::define() 19 7.075e-05 7.609e-05 8.133e-05 0.00% Castro::create_source_corrector() 17 5.897e-05 6.445e-05 7.244e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.841e-05 6.019e-05 6.311e-05 0.00% Castro::expand_state() 20 5.142e-05 5.397e-05 5.783e-05 0.00% Amr::FinalizeInit() 1 1.702e-05 3.106e-05 5.735e-05 0.00% MLMG::computeMLResidual() 18 4.122e-05 4.398e-05 4.79e-05 0.00% Castro::derive() 18 4.022e-05 4.162e-05 4.485e-05 0.00% Amr::regrid() 1 1.641e-05 2.018e-05 4.442e-05 0.00% Amr::timeStep() 1 1.44e-05 2.373e-05 3.687e-05 0.00% Gravity::actual_multilevel_solve() 1 3.133e-05 3.218e-05 3.293e-05 0.00% Amr::InitializeInit() 1 7.575e-06 1.17e-05 2.793e-05 0.00% Castro::post_regrid() 2 1.086e-05 1.187e-05 1.685e-05 0.00% Castro::advance() 1 6.789e-06 8.608e-06 1.493e-05 0.00% Amr::bldFineLevels() 1 6.967e-06 8.411e-06 1.46e-05 0.00% Castro::post_init() 1 1.175e-05 1.263e-05 1.416e-05 0.00% Amr::initialInit() 1 3.379e-06 7.114e-06 1.095e-05 0.00% Amr::init() 1 6.045e-06 7.168e-06 8.773e-06 0.00% Castro::post_timestep() 1 5.45e-06 6.253e-06 7.524e-06 0.00% Other 2421 0.0487 0.06252 0.09125 0.63% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.46 14.47 14.47 100.00% Amr::coarseTimeStep() 1 12.23 12.23 12.23 84.49% Amr::timeStep() 1 12.18 12.18 12.18 84.13% Castro::advance() 1 12.11 12.11 12.12 83.72% Castro::subcycle_advance_ctu() 1 12.06 12.07 12.07 83.38% Castro::do_advance_ctu() 20 12.02 12.03 12.03 83.12% Castro::construct_ctu_hydro_source() 20 6.618 6.644 6.694 46.25% Castro::clean_state() 79 3.617 3.733 3.888 26.86% Castro::computeTemp() 80 2.305 2.415 2.57 17.76% Gravity::solve_phi_with_mlmg() 19 1.622 1.675 1.701 11.75% Castro::construct_new_gravity() 17 1.536 1.588 1.613 11.15% Gravity::solve_for_phi() 17 1.498 1.551 1.576 10.89% Castro::reset_internal_energy(MultiFab) 80 1.481 1.492 1.528 10.56% Castro::initialize_do_advance() 20 1.367 1.406 1.441 9.96% Gravity::actual_solve_with_mlmg() 19 1.422 1.426 1.43 9.88% MLMG::solve() 18 1.368 1.374 1.375 9.50% Amr::InitAmr() 1 1.345 1.348 1.354 9.36% MLMG::oneIter() 128 1.225 1.229 1.231 8.51% MLMG::mgVcycle() 128 1.21 1.214 1.216 8.40% Castro::do_old_sources() 20 1.129 1.151 1.18 8.15% Castro::do_new_sources() 17 0.9101 0.9561 0.9974 6.89% Castro::evaluate_source_change() 196 0.9492 0.9559 0.9852 6.81% FabArray::ParallelCopy() 1429 0.471 0.7786 0.9678 6.69% Amr::init() 1 0.8895 0.8896 0.8896 6.15% MLCellLinOp::smooth() 1544 0.6401 0.6949 0.8356 5.77% Amr::initialInit() 1 0.8352 0.8352 0.8353 5.77% Castro::enforce_min_density() 79 0.7911 0.7971 0.8227 5.68% amrex::Copy() 724 0.8029 0.8084 0.8218 5.68% FabArray::ParallelCopy_finish() 1429 0.2212 0.5288 0.7203 4.98% AmrLevel::FillPatch() 147 0.5256 0.5994 0.6725 4.65% Amr::InitializeInit() 1 0.4357 0.4664 0.5951 4.11% Amr::defBaseLevel() 1 0.4357 0.4664 0.5951 4.11% MultiFab::AddProduct() 3439 0.5462 0.5501 0.5673 3.92% Castro::normalize_species() 79 0.5429 0.5472 0.5602 3.87% Castro::reset_internal_energy(Fab) 160 0.5037 0.5085 0.5192 3.59% FillPatchIterator::Initialize 147 0.3412 0.4148 0.4894 3.38% MLMG::mgVcycle_up::1 128 0.1149 0.2914 0.4679 3.23% FillPatchIterator::FillFromLevel0() 147 0.3136 0.394 0.4663 3.22% FillPatchSingleLevel 147 0.3129 0.3933 0.4656 3.22% MLCellLinOp::applyBC() 7148 0.3347 0.3804 0.4651 3.21% MLPoisson::Fsmooth() 6144 0.3472 0.3754 0.4506 3.11% MLMG::addInterpCorrection() 768 0.06538 0.3232 0.4009 2.77% Amr::FinalizeInit() 1 0.2401 0.3688 0.3995 2.76% Castro::Castro() 1 0.3561 0.3629 0.3782 2.61% Castro::initMFs() 1 0.3501 0.3568 0.3726 2.57% FabArray::setVal() 5013 0.344 0.3486 0.3578 2.47% MultiFab::Subtract() 159 0.3345 0.3372 0.3514 2.43% FabArray::FillBoundary() 6380 0.2614 0.3024 0.3476 2.40% MLMG::mgVcycle_up::2 128 0.01329 0.1264 0.2982 2.06% MLMG::mgVcycle_down::0 128 0.2756 0.2784 0.281 1.94% Gravity::fill_multipole_BCs() 19 0.1985 0.2478 0.2689 1.86% FabArray::ParallelCopy_nowait() 1429 0.2278 0.2461 0.2672 1.85% Castro::expand_state() 20 0.1969 0.2274 0.2601 1.80% check_for_negative_density() 20 0.07739 0.1764 0.2476 1.71% Castro::initData() 1 0.077 0.1016 0.232 1.60% MLMG::mgVcycle_up::0 128 0.2167 0.219 0.2206 1.52% MLCellLinOp::apply() 986 0.1615 0.1729 0.2105 1.45% Amr::bldFineLevels() 1 0.04533 0.1744 0.205 1.42% FabArray::sum() 3439 0.1925 0.1948 0.2003 1.38% MLMG::mgVcycle_up::3 128 0.01323 0.02695 0.1913 1.32% FillBoundary_finish() 6380 0.1254 0.1561 0.19 1.31% FillBoundary_nowait() 6380 0.1027 0.132 0.1684 1.16% Amr::grid_places() 2 0.002619 0.1314 0.162 1.12% AmrMesh::MakeNewGrids() 2 0.002533 0.1313 0.162 1.12% TagBoxArray::mapPRD 2 0.0008965 0.1033 0.1602 1.11% MLMG:computeResOfCorrection() 768 0.1198 0.1301 0.1547 1.07% MLCellLinOp::correctionResidual() 768 0.1183 0.1287 0.1531 1.06% Castro::post_init() 1 0.1474 0.1517 0.1527 1.05% TagBoxArray::collate() 2 0.0002248 0.02662 0.1404 0.97% MLMG::mgVcycle_down::1 128 0.1303 0.1362 0.1403 0.97% Castro::construct_new_source() 102 0.1246 0.1256 0.1294 0.89% MLPoisson::Fapply() 986 0.0966 0.1013 0.1158 0.80% Gravity::multilevel_solve_for_new_phi() 1 0.1124 0.1124 0.1126 0.78% Gravity::actual_multilevel_solve() 1 0.1124 0.1124 0.1125 0.78% MLCellLinOp::solutionResidual() 147 0.09298 0.09925 0.105 0.73% Amr::writePlotFile() 2 0.09795 0.09802 0.09804 0.68% Castro::post_regrid() 2 0.08472 0.08573 0.09006 0.62% Castro::create_source_corrector() 17 0.05558 0.06761 0.08782 0.61% MLMG::computeResidual() 128 0.08043 0.08164 0.08319 0.57% Castro::construct_old_source() 120 0.07958 0.08013 0.08238 0.57% MLMG::mgVcycle_down::2 128 0.01557 0.04661 0.0808 0.56% FabArray::Saxpy() 264 0.07246 0.07317 0.07992 0.55% Castro::derive() 18 0.06959 0.07012 0.07346 0.51% AmrLevel::derive() 18 0.06954 0.07008 0.07342 0.51% Castro::construct_new_gravity_source() 17 0.06832 0.0693 0.07318 0.51% MLMG::mgVcycle_bottom 128 0.000351 0.004907 0.0731 0.51% MLMG::actualBottomSolve() 128 0.0001301 0.00468 0.07286 0.50% MLCGSolver::bicgstab 8 0 0.004439 0.07103 0.49% MLMG::mgVcycle_down::3 128 0.01488 0.02156 0.06725 0.46% Castro::construct_old_gravity() 20 0.04949 0.05869 0.06484 0.45% Gravity::get_old_grav_vector() 20 0.04937 0.05858 0.06473 0.45% Castro::post_timestep() 1 0.05918 0.06167 0.06396 0.44% FabArray::Xpay() 945 0.0518 0.05431 0.06092 0.42% amrex::average_down 768 0.03388 0.04093 0.05842 0.40% Castro::initialize_advance() 1 0.04858 0.04903 0.05148 0.36% Castro::construct_old_gravity_source() 20 0.04578 0.04608 0.04722 0.33% Amr::regrid() 1 0.04378 0.04426 0.04657 0.32% Castro::check_for_nan() 37 0.04404 0.04452 0.04586 0.32% MultiFab::contains_nan() 37 0.04394 0.04441 0.04572 0.32% MLPoisson::define() 19 0.03015 0.03526 0.04103 0.28% MLMG::mgVcycle_down::4 128 0.01477 0.01646 0.03928 0.27% Castro::retry_advance_ctu() 11 0.03599 0.03671 0.039 0.27% Gravity::get_new_grav_vector() 18 0.03628 0.03683 0.03744 0.26% StateDataPhysBCFunct::() 147 0.0006931 0.01829 0.03469 0.24% StateData::FillBoundary(geom) 143 0 0.01744 0.03378 0.23% MLMG::mgVcycle_up::4 128 0.01335 0.01454 0.03117 0.22% Castro::apply_source_to_state() 37 0.02956 0.02969 0.02995 0.21% MLMG::mgVcycle_down::5 128 0.0148 0.01571 0.02747 0.19% MLCellLinOp::defineAuxData() 19 0.01976 0.02139 0.02715 0.19% VisMF::Write(FabArray) 4 0.0136 0.01467 0.02604 0.18% MLMG::prepareForSolve() 18 0.02047 0.02137 0.02433 0.17% MLMG::mgVcycle_up::5 128 0.01346 0.01418 0.02344 0.16% Castro::sum_integrated_quantities() 2 0.02285 0.023 0.02301 0.16% MLMG::computeMLResidual() 18 0.01095 0.01709 0.02259 0.16% MLMG::ResNormInf() 146 0.01593 0.0189 0.02258 0.16% Castro::save_data_for_retry() 3 0.01933 0.01946 0.02009 0.14% Castro::problem_post_timestep() 2 0.009682 0.01253 0.01523 0.11% Castro::wd_update() 2 0.009668 0.0125 0.01522 0.11% Gravity::test_composite_phi() 1 0.01503 0.01517 0.0152 0.11% amrex::Add() 128 0.01429 0.01444 0.01503 0.10% FabArrayBase::getCPC() 2005 0.009692 0.0111 0.01489 0.10% FabArray::setDomainBndry() 147 0.0004527 0.007847 0.01486 0.10% FabArrayBase::getFB() 6380 0.01063 0.01141 0.01477 0.10% FabArray::norminf() 244 0.01097 0.01134 0.01403 0.10% MLLinOp::define() 19 0.003841 0.007408 0.01114 0.08% Other 2421 0.08132 0.1028 0.1408 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 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 1569 KiB 1571 KiB 1571 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 409 KiB 413 KiB 420 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1040 KiB 1041 KiB 1041 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 990 KiB 991 KiB 991 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 664 KiB 665 KiB 665 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 4272 B 4343 B 4875 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 50 KiB 50 KiB 51 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 11 KiB 13 KiB 14 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 416 KiB 416 KiB 416 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 2059 B 3857 B 5038 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1329 B 1339 B 1398 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 17 KiB 17 KiB 18 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 132 KiB 132 KiB 132 KiB 157 KiB 157 KiB 157 KiB Castro::buildMetrics() 192 192 128 KiB 128 KiB 128 KiB 142 KiB 142 KiB 142 KiB MLMG::prepareForSolve() 6102 6102 9927 B 10 KiB 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 260 B 262 B 271 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 438 B 439 B 444 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 19 B 21 B 37 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 33 KiB 33 KiB 33 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 54 B 54 B 55 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 3 B 8 B 9 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 113 B 114 B 117 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 155 B 184 B 204 B 45 KiB 45 KiB 45 KiB Gravity::multilevel_solve_for_new_phi() 64 64 31 KiB 31 KiB 31 KiB 37 KiB 37 KiB 37 KiB Gravity::fill_multipole_BCs() 2432 2432 229 B 236 B 245 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1068 B 1076 B 1109 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 127 B 127 B 127 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1695 B 1755 B 1784 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 448 B 501 B 695 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 466 B 468 B 469 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 27 B 38 B 50 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 18 B 28 B 2916 B 2916 B 2916 B OwnerMask() 64 64 0 B 0 B 0 B 2592 B 2592 B 2592 B MLMG::addInterpCorrection() 5632 5632 1 B 16 B 25 B 1024 B 1024 B 1024 B amrex::average_down 5632 5632 0 B 0 B 1 B 1024 B 1024 B 1024 B TagBoxArray::mapPRD 64 64 0 B 5 B 7 B 648 B 648 B 648 B MLCGSolver::bicgstab 640 640 0 B 0 B 2 B 0 B 36 B 576 B MLCellLinOp::applyBC() 114375 114375 0 B 0 B 0 B 8 B 8 B 8 B ------------------------------------------------------------------------------------------------------------------------------- AMReX (25.03-15-g546cc5e16431-dirty) finalized