Initializing AMReX (25.06-14-g65bf80d763f6)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (25.06-14-g65bf80d763f6) initialized Starting run at 12:52:53 UTC on 2025-06-13. Successfully read inputs file ... Castro git describe: 25.06-3-gb5799437e AMReX git describe: 25.06-14-g65bf80d76 Microphysics git describe: 25.06-11-gf3ca3777 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.18738446 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001250827 new finest: 0 Gravity::fill_multipole_BCs() time = 0.011779978 Gravity::multilevel_solve_for_new_phi() time = 0.10876304 ... norm of composite residual at level 0 1.197718084e-11 Gravity::fill_multipole_BCs() time = 0.010738444 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.049062611 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004086952 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 = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071148 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.769409 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477044 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035502 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -0.01295059363 at index 21, 83, 0 Invalid density = -0.01295059363 at index 21, 172, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035502 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477044 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.769409 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071148 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Castro::construct_old_gravity() time = 0.005651903 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.064756002 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 = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Castro::construct_old_gravity() time = 0.005600483 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.063311959 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, 177, 0 Invalid density = -0.1060744788 at index 20, 78, 0 Castro::construct_old_gravity() time = 0.007085178 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.060743695 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.00453072 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.059188873 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.367771604 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.028955134 Gravity::solve_for_phi() time = 0.117097906 on level 0 Castro::construct_new_gravity() time = 0.119258543 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.066397381 on level 0 Castro::construct_old_gravity() time = 0.005310285 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.058834635 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.332306122 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.014306262 Gravity::solve_for_phi() time = 0.092111473 on level 0 Castro::construct_new_gravity() time = 0.09422016 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.061021568 on level 0 Castro::construct_old_gravity() time = 0.007723787 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.060009285 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.356753281 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.022923601 Gravity::solve_for_phi() time = 0.106646639 on level 0 Castro::construct_new_gravity() time = 0.108759491 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.060422311 on level 0 Castro::construct_old_gravity() time = 0.003096845 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.059672528 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.348196917 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.015440026 Gravity::solve_for_phi() time = 0.093261534 on level 0 Castro::construct_new_gravity() time = 0.095444103 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.060250487 on level 0 Castro::construct_old_gravity() time = 0.004511098 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.058940779 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.344606773 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.014767342 Gravity::solve_for_phi() time = 0.090552551 on level 0 Castro::construct_new_gravity() time = 0.092666524 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.058966499 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.003381203 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.058758245 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.339450615 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.014762956 Gravity::solve_for_phi() time = 0.091396735 on level 0 Castro::construct_new_gravity() time = 0.093921909 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.059419313 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.074198974 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.073570536 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.40467892 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.014604817 Gravity::solve_for_phi() time = 0.090951645 on level 0 Castro::construct_new_gravity() time = 0.093048644 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.060991596 on level 0 Castro::construct_old_gravity() time = 0.004440145 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.059574678 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.348320931 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.019528906 Gravity::solve_for_phi() time = 0.09549099 on level 0 Castro::construct_new_gravity() time = 0.097599173 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.060778414 on level 0 Castro::construct_old_gravity() time = 0.002730907 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.059098227 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.345094371 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.014799132 Gravity::solve_for_phi() time = 0.112571627 on level 0 Castro::construct_new_gravity() time = 0.122349405 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.068179385 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.006585872 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.061640092 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.34800058 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.015838379 Gravity::solve_for_phi() time = 0.09258065 on level 0 Castro::construct_new_gravity() time = 0.094689311 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.058909849 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.00295573 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.059911967 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.349235846 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.014637585 Gravity::solve_for_phi() time = 0.090331927 on level 0 Castro::construct_new_gravity() time = 0.092508391 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.058984433 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.003523418 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.059401147 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.341901937 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.015457379 Gravity::solve_for_phi() time = 0.095738235 on level 0 Castro::construct_new_gravity() time = 0.097872939 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.060809655 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.002722455 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.058754495 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.348548416 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.015601483 Gravity::solve_for_phi() time = 0.091082333 on level 0 Castro::construct_new_gravity() time = 0.093252678 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.060039529 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.003052415 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.058994916 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.347155692 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.014902498 Gravity::solve_for_phi() time = 0.092896647 on level 0 Castro::construct_new_gravity() time = 0.094998836 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.059852025 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.003665952 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.059831369 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.354709442 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.016502804 Gravity::solve_for_phi() time = 0.09491128 on level 0 Castro::construct_new_gravity() time = 0.09702944 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.060317743 on level 0 Castro::construct_old_gravity() time = 0.002858319 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.059161357 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.341195124 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.014936341 Gravity::solve_for_phi() time = 0.090786702 on level 0 Castro::construct_new_gravity() time = 0.092900496 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.05952957 on level 0 Castro::construct_old_gravity() time = 0.003936132 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.060176885 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.347647636 on level 0 Ending SDC iteration 2 of 2. [1m[32m Subcycle completed[0m Subcycling complete Zones advanced per microsecond at this level: 0.002638335846 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.014944441 Gravity::solve_for_phi() time = 0.091723041 on level 0 Castro::construct_new_gravity() time = 0.093870408 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.061230093 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.51708166 [STEP 1] FAB kilobyte spread across MPI nodes: [6998 ... 6998] PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.051709211 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004187933 seconds Ending run at 12:53:08 UTC on 2025-06-13. Run time = 14.89949902 Run time without initialization = 12.5738916 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.8915 ... 14.8934 ... 14.8995 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 20 5.6551 5.7146 5.8658 39.37% Amr::InitAmr() 1 1.407 1.409 1.415 9.50% Castro::computeTemp() 80 0.8277 0.9385 1.1 7.38% FabArray::ParallelCopy_finish() 1429 0.2824 0.6712 0.8922 5.99% amrex::Copy() 724 0.7755 0.7844 0.8007 5.37% MultiFab::AddProduct() 3439 0.5753 0.5809 0.595 3.99% Castro::normalize_species() 79 0.5479 0.554 0.5915 3.97% Castro::reset_internal_energy(Fab) 160 0.5095 0.5164 0.5341 3.58% MLPoisson::Fsmooth() 6144 0.3527 0.3828 0.4683 3.14% Castro::initMFs() 1 0.3492 0.359 0.3739 2.51% FabArray::setVal() 5013 0.3466 0.3523 0.3644 2.45% MultiFab::Subtract() 159 0.3388 0.3414 0.3499 2.35% check_for_negative_density() 20 0.07524 0.2499 0.3219 2.16% Gravity::fill_multipole_BCs() 19 0.2098 0.2779 0.3006 2.02% FabArray::ParallelCopy_nowait() 1429 0.2023 0.2186 0.2375 1.59% FabArray::sum() 3439 0.1917 0.1939 0.199 1.34% FillBoundary_finish() 6380 0.1328 0.1528 0.1786 1.20% TagBoxArray::collate() 2 0.0002627 0.03724 0.1666 1.12% Castro::initData() 1 0.003418 0.02803 0.1631 1.09% FillBoundary_nowait() 6380 0.09217 0.1211 0.1543 1.04% MLCellLinOp::applyBC() 7148 0.06178 0.07899 0.1278 0.86% MLPoisson::Fapply() 986 0.09285 0.09807 0.1148 0.77% Castro::construct_new_gravity_source() 17 0.07134 0.07308 0.07975 0.54% FabArray::Saxpy() 91 0.06706 0.06743 0.06844 0.46% FabArray::Xpay() 915 0.05088 0.05343 0.05997 0.40% Castro::construct_old_gravity_source() 20 0.04824 0.04865 0.0513 0.34% Castro::wd_update() 2 0.006164 0.04476 0.05006 0.34% MultiFab::contains_nan() 37 0.04639 0.04689 0.04819 0.32% MLMG::addInterpCorrection() 768 0.0336 0.0363 0.04473 0.30% AmrLevel::derive() 18 0.04005 0.04088 0.04449 0.30% amrex::average_down 768 0.02516 0.02756 0.03486 0.23% StateData::FillBoundary(geom) 143 0 0.01765 0.03464 0.23% VisMF::Write(FabArray) 4 0.01415 0.01546 0.02876 0.19% Castro::reset_internal_energy(MultiFab) 80 0.02324 0.02388 0.02519 0.17% Amr::writePlotFile() 2 0.006986 0.01893 0.02097 0.14% MLCellLinOp::defineAuxData() 19 0.01581 0.01678 0.02041 0.14% Castro::enforce_min_density() 79 0.01854 0.01886 0.01916 0.13% FabArray::FillBoundary() 6380 0.01356 0.01418 0.01705 0.11% MLCGSolver::bicgstab 8 0 0.001051 0.01682 0.11% Castro::construct_new_source() 102 0.0158 0.01593 0.01624 0.11% amrex::Add() 128 0.01426 0.01452 0.0153 0.10% Castro::evaluate_source_change() 196 0.01353 0.01407 0.01504 0.10% MLMG::ResNormInf() 146 0.006626 0.01175 0.01488 0.10% FabArray::setDomainBndry() 147 0.0004747 0.007613 0.01483 0.10% FabArray::norminf() 244 0.01082 0.01139 0.01409 0.09% FillPatchIterator::Initialize 147 0.01284 0.01312 0.01407 0.09% FabArrayBase::getFB() 6380 0.01007 0.01059 0.01308 0.09% MLMG::prepareForSolve() 18 0.01023 0.01069 0.01177 0.08% MLCellLinOp::smooth() 1544 0.009521 0.01015 0.01168 0.08% MLMG::solve() 18 0.001248 0.004767 0.01072 0.07% Castro::construct_old_source() 120 0.01009 0.01021 0.01041 0.07% Gravity::actual_solve_with_mlmg() 19 0.007297 0.007815 0.009338 0.06% FabArrayBase::getCPC() 2005 0.007341 0.007578 0.008103 0.05% Gravity::get_new_grav_vector() 18 0.006014 0.006252 0.007743 0.05% Gravity::get_old_grav_vector() 20 0.006517 0.006609 0.006695 0.04% Castro::sum_integrated_quantities() 2 0.003056 0.004727 0.005226 0.04% Castro::initialize_do_advance() 20 0.004595 0.004836 0.005052 0.03% MLCellLinOp::apply() 986 0.001914 0.002161 0.004588 0.03% FabArray::ParallelCopy() 1429 0.003478 0.003629 0.003987 0.03% AmrLevel::FillPatch() 147 0.002371 0.002514 0.002836 0.02% MLLinOp::defineGrids() 19 0.001988 0.002123 0.002415 0.02% MLCellLinOp::correctionResidual() 768 0.001867 0.001931 0.002046 0.01% MLMG::mgVcycle() 128 0.001679 0.001767 0.001935 0.01% MLMG:computeResOfCorrection() 768 0.001477 0.001544 0.00166 0.01% Castro::save_data_for_retry() 3 0.001315 0.001357 0.001614 0.01% MLMG::actualBottomSolve() 128 0.0001342 0.000207 0.001237 0.01% Castro::subcycle_advance_ctu() 1 0.0002261 0.0003104 0.001036 0.01% StateDataPhysBCFunct::() 147 0.0007278 0.0008941 0.0009975 0.01% main() 1 0.0006827 0.0007257 0.0008 0.01% FillPatchIterator::FillFromLevel0() 147 0.0006733 0.0007061 0.0007605 0.01% Gravity::solve_for_phi() 17 0.0006547 0.0006936 0.0007467 0.01% AmrMesh::MakeNewGrids() 2 0.0005475 0.0005962 0.0006361 0.00% Castro::Castro() 1 0.0003974 0.0004804 0.0006312 0.00% MLMG::mgVcycle_down::0 128 0.0005517 0.0005736 0.0006235 0.00% MLMG::mgVcycle_down::2 128 0.0005496 0.0005875 0.0006228 0.00% MLMG::mgVcycle_down::1 128 0.0005488 0.0005708 0.0006201 0.00% MLMG::mgVcycle_down::3 128 0.0005129 0.0005516 0.0006194 0.00% MLMG::mgVcycle_down::4 128 0.0004911 0.0005166 0.0006009 0.00% MLMG::mgVcycle_down::5 128 0.0005076 0.0005289 0.00059 0.00% Castro::clean_state() 79 0.0004198 0.000459 0.000501 0.00% FillPatchSingleLevel 147 0.0004693 0.0004822 0.0004968 0.00% MLCellLinOp::solutionResidual() 147 0.0004371 0.0004509 0.0004734 0.00% MLMG::oneIter() 128 0.0003815 0.0004081 0.0004564 0.00% Castro::do_advance_ctu() 20 0.000381 0.0004148 0.000443 0.00% MLMG::mgVcycle_up::0 128 0.0003593 0.0003755 0.0004095 0.00% MLMG::mgVcycle_up::1 128 0.0003575 0.0003719 0.0004052 0.00% MLMG::mgVcycle_up::2 128 0.0003155 0.0003458 0.0003878 0.00% MLMG::mgVcycle_up::4 128 0.0003158 0.0003306 0.000367 0.00% MLMG::mgVcycle_up::3 128 0.0003155 0.0003322 0.000367 0.00% MLMG::computeResidual() 128 0.0002356 0.000251 0.0003649 0.00% MLMG::mgVcycle_up::5 128 0.0003155 0.0003234 0.0003421 0.00% Castro::do_old_sources() 20 0.0002745 0.0002895 0.0003063 0.00% Castro::construct_new_gravity() 17 0.0001652 0.0001819 0.0003027 0.00% Castro::problem_post_timestep() 2 1.138e-05 3.021e-05 0.0002731 0.00% MLMG::mgVcycle_bottom 128 0.0002299 0.0002368 0.0002576 0.00% Gravity::solve_phi_with_mlmg() 19 0.0002174 0.000237 0.0002566 0.00% Castro::do_new_sources() 17 0.0002078 0.0002211 0.0002385 0.00% Amr::coarseTimeStep() 1 0.0001317 0.0001661 0.0001812 0.00% Amr::defBaseLevel() 1 0.0001033 0.0001295 0.0001484 0.00% MLPoisson::define() 19 0.0001221 0.0001332 0.0001437 0.00% Castro::check_for_nan() 37 0.0001046 0.000118 0.0001435 0.00% TagBoxArray::mapPRD 2 0.0001263 0.0001337 0.0001425 0.00% Castro::retry_advance_ctu() 11 6.475e-05 7.579e-05 0.0001361 0.00% Castro::initialize_advance() 1 0.00011 0.0001219 0.0001337 0.00% Castro::construct_old_gravity() 20 0.0001074 0.0001148 0.0001316 0.00% Gravity::test_composite_phi() 1 0.0001052 0.0001097 0.0001247 0.00% Castro::apply_source_to_state() 37 9.822e-05 0.0001029 0.0001113 0.00% MLLinOp::define() 19 8.384e-05 8.85e-05 9.842e-05 0.00% Amr::grid_places() 2 5.386e-05 6.366e-05 9.789e-05 0.00% Castro::create_source_corrector() 17 5.85e-05 6.631e-05 7.746e-05 0.00% Amr::FinalizeInit() 1 8.32e-06 1.906e-05 7.471e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.883e-05 6.034e-05 6.441e-05 0.00% Castro::expand_state() 20 5.076e-05 5.343e-05 5.6e-05 0.00% MLMG::computeMLResidual() 18 4.127e-05 4.347e-05 4.891e-05 0.00% Castro::derive() 18 4.064e-05 4.27e-05 4.535e-05 0.00% Amr::regrid() 1 1.57e-05 1.974e-05 4.37e-05 0.00% Amr::timeStep() 1 1.415e-05 1.779e-05 3.808e-05 0.00% Gravity::actual_multilevel_solve() 1 3.091e-05 3.178e-05 3.508e-05 0.00% Amr::init() 1 6.336e-06 8.202e-06 2.388e-05 0.00% Castro::post_regrid() 2 1.11e-05 1.361e-05 2.342e-05 0.00% Castro::advance() 1 6.96e-06 1.104e-05 2.213e-05 0.00% Amr::InitializeInit() 1 7.181e-06 9.27e-06 1.774e-05 0.00% Castro::post_init() 1 1.187e-05 1.301e-05 1.643e-05 0.00% Castro::post_timestep() 1 5.89e-06 7.436e-06 1.268e-05 0.00% Amr::bldFineLevels() 1 6.713e-06 7.447e-06 8.314e-06 0.00% Amr::initialInit() 1 2.942e-06 3.173e-06 3.757e-06 0.00% Other 2502 0.04774 0.06193 0.09786 0.66% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 14.89 14.89 14.9 100.00% Amr::coarseTimeStep() 1 12.57 12.57 12.57 84.39% Amr::timeStep() 1 12.52 12.52 12.52 84.01% Castro::advance() 1 12.41 12.42 12.42 83.36% Castro::subcycle_advance_ctu() 1 12.36 12.37 12.37 83.03% Castro::do_advance_ctu() 20 12.32 12.33 12.33 82.77% Castro::construct_ctu_hydro_source() 20 6.714 6.743 6.799 45.64% Castro::clean_state() 79 3.663 3.783 3.95 26.51% Castro::computeTemp() 80 2.331 2.448 2.612 17.53% Gravity::solve_phi_with_mlmg() 19 1.656 1.723 1.75 11.74% Castro::construct_new_gravity() 17 1.583 1.645 1.671 11.21% Gravity::solve_for_phi() 17 1.536 1.603 1.63 10.94% Castro::reset_internal_energy(MultiFab) 80 1.492 1.509 1.544 10.36% Castro::initialize_do_advance() 20 1.4 1.493 1.542 10.35% Gravity::actual_solve_with_mlmg() 19 1.437 1.442 1.446 9.70% Amr::InitAmr() 1 1.407 1.409 1.415 9.50% MLMG::solve() 18 1.382 1.388 1.389 9.32% MLMG::oneIter() 128 1.238 1.241 1.243 8.34% MLMG::mgVcycle() 128 1.223 1.226 1.228 8.24% Castro::do_old_sources() 20 1.134 1.176 1.199 8.05% FabArray::ParallelCopy() 1429 0.5319 0.9133 1.13 7.59% Castro::do_new_sources() 17 0.9208 0.9676 1.018 6.83% Castro::evaluate_source_change() 196 0.9792 0.9888 1.013 6.80% Amr::init() 1 0.9106 0.9107 0.9107 6.11% FabArray::ParallelCopy_finish() 1429 0.2824 0.6712 0.8922 5.99% Amr::initialInit() 1 0.8569 0.857 0.8572 5.75% MLCellLinOp::smooth() 1544 0.6368 0.6964 0.8481 5.69% Castro::enforce_min_density() 79 0.8009 0.8076 0.8255 5.54% amrex::Copy() 724 0.7755 0.7844 0.8007 5.37% AmrLevel::FillPatch() 147 0.6118 0.7076 0.7751 5.20% Amr::InitializeInit() 1 0.435 0.47 0.6209 4.17% Amr::defBaseLevel() 1 0.4349 0.4699 0.6209 4.17% FillPatchIterator::Initialize 147 0.4329 0.5289 0.598 4.01% MultiFab::AddProduct() 3439 0.5753 0.5809 0.595 3.99% Castro::normalize_species() 79 0.5479 0.554 0.5915 3.97% FillPatchIterator::FillFromLevel0() 147 0.4052 0.5081 0.5797 3.89% FillPatchSingleLevel 147 0.4045 0.5074 0.579 3.89% Castro::reset_internal_energy(Fab) 160 0.5095 0.5164 0.5341 3.58% MLMG::mgVcycle_up::1 128 0.1091 0.2922 0.4777 3.21% MLPoisson::Fsmooth() 6144 0.3527 0.3828 0.4683 3.14% MLCellLinOp::applyBC() 7148 0.3332 0.3784 0.4599 3.09% Amr::FinalizeInit() 1 0.2362 0.3871 0.4221 2.83% MLMG::addInterpCorrection() 768 0.06769 0.3322 0.4137 2.78% Castro::Castro() 1 0.3546 0.3649 0.3813 2.56% Castro::initMFs() 1 0.3493 0.3592 0.3741 2.51% FabArray::setVal() 5013 0.3466 0.3523 0.3644 2.45% MultiFab::Subtract() 159 0.3388 0.3414 0.3499 2.35% FabArray::FillBoundary() 6380 0.2578 0.2994 0.3458 2.32% check_for_negative_density() 20 0.07524 0.2499 0.3219 2.16% MLMG::mgVcycle_up::2 128 0.0133 0.1298 0.3065 2.06% Gravity::fill_multipole_BCs() 19 0.2119 0.28 0.3027 2.03% MLMG::mgVcycle_down::0 128 0.2792 0.2808 0.284 1.91% Castro::expand_state() 20 0.2004 0.238 0.2692 1.81% FabArray::ParallelCopy_nowait() 1429 0.2215 0.2384 0.2584 1.73% Castro::initData() 1 0.07814 0.1032 0.2377 1.60% Amr::bldFineLevels() 1 0.04576 0.1967 0.2317 1.55% MLMG::mgVcycle_up::0 128 0.2172 0.2179 0.2194 1.47% MLCellLinOp::apply() 986 0.1591 0.173 0.2091 1.40% FabArray::sum() 3439 0.1917 0.1939 0.199 1.34% MLMG::mgVcycle_up::3 128 0.01326 0.0275 0.195 1.31% Amr::grid_places() 2 0.002689 0.1536 0.1886 1.27% AmrMesh::MakeNewGrids() 2 0.002624 0.1535 0.1885 1.27% Castro::create_source_corrector() 17 0.06445 0.1256 0.1804 1.21% FillBoundary_finish() 6380 0.1328 0.1528 0.1786 1.20% TagBoxArray::mapPRD 2 0.000946 0.1149 0.1765 1.18% TagBoxArray::collate() 2 0.0002627 0.03724 0.1666 1.12% FillBoundary_nowait() 6380 0.1026 0.1324 0.1662 1.12% MLMG:computeResOfCorrection() 768 0.119 0.1304 0.1555 1.04% MLCellLinOp::correctionResidual() 768 0.1175 0.1289 0.1538 1.03% Castro::post_init() 1 0.1439 0.1474 0.1481 0.99% Castro::construct_old_gravity() 20 0.04731 0.08311 0.1471 0.99% Gravity::get_old_grav_vector() 20 0.0472 0.08299 0.147 0.99% MLMG::mgVcycle_down::1 128 0.1343 0.1388 0.1442 0.97% Castro::construct_new_source() 102 0.1253 0.1273 0.134 0.90% MLPoisson::Fapply() 986 0.09285 0.09807 0.1148 0.77% Gravity::multilevel_solve_for_new_phi() 1 0.1087 0.1087 0.1088 0.73% Gravity::actual_multilevel_solve() 1 0.1086 0.1086 0.1087 0.73% MLCellLinOp::solutionResidual() 147 0.09051 0.09824 0.1058 0.71% Castro::post_timestep() 1 0.09706 0.1008 0.1032 0.69% Amr::writePlotFile() 2 0.1011 0.1012 0.1012 0.68% Castro::post_regrid() 2 0.08534 0.08608 0.08955 0.60% MLMG::mgVcycle_down::2 128 0.01558 0.04879 0.086 0.58% Castro::construct_old_source() 120 0.08091 0.08145 0.08466 0.57% MLMG::computeResidual() 128 0.07836 0.07993 0.08239 0.55% Castro::construct_new_gravity_source() 17 0.07134 0.07308 0.07975 0.54% Castro::derive() 18 0.06883 0.06987 0.07329 0.49% AmrLevel::derive() 18 0.06879 0.06983 0.07324 0.49% MLMG::mgVcycle_bottom 128 0.0003672 0.004861 0.07218 0.48% MLMG::actualBottomSolve() 128 0.0001342 0.004625 0.07192 0.48% MLCGSolver::bicgstab 8 0 0.004378 0.07005 0.47% FabArray::Saxpy() 91 0.06706 0.06743 0.06844 0.46% MLMG::mgVcycle_down::3 128 0.01484 0.02187 0.0681 0.46% amrex::average_down 768 0.03393 0.04206 0.06054 0.41% FabArray::Xpay() 915 0.05088 0.05343 0.05997 0.40% Castro::problem_post_timestep() 2 0.009922 0.04857 0.0535 0.36% Castro::wd_update() 2 0.009903 0.04854 0.05349 0.36% Castro::initialize_advance() 1 0.04918 0.05007 0.05324 0.36% Castro::construct_old_gravity_source() 20 0.04824 0.04865 0.0513 0.34% Castro::check_for_nan() 37 0.0465 0.04701 0.04832 0.32% MultiFab::contains_nan() 37 0.04639 0.04689 0.04819 0.32% Gravity::get_new_grav_vector() 18 0.03564 0.04091 0.04622 0.31% Amr::regrid() 1 0.04404 0.04434 0.04553 0.31% MLPoisson::define() 19 0.031 0.03632 0.04202 0.28% MLMG::mgVcycle_down::4 128 0.01475 0.0167 0.04033 0.27% Castro::retry_advance_ctu() 11 0.03595 0.03625 0.03756 0.25% StateDataPhysBCFunct::() 147 0.0007278 0.01854 0.03564 0.24% StateData::FillBoundary(geom) 143 0 0.01765 0.03464 0.23% MLMG::mgVcycle_up::4 128 0.01337 0.01478 0.03167 0.21% VisMF::Write(FabArray) 4 0.01415 0.01546 0.02876 0.19% Castro::apply_source_to_state() 37 0.02734 0.02757 0.02797 0.19% MLMG::mgVcycle_down::5 128 0.0147 0.01584 0.02778 0.19% MLCellLinOp::defineAuxData() 19 0.01994 0.02172 0.02757 0.19% MLMG::prepareForSolve() 18 0.02053 0.02165 0.02458 0.16% Castro::sum_integrated_quantities() 2 0.02435 0.02452 0.02454 0.16% MLMG::ResNormInf() 146 0.01665 0.02133 0.02448 0.16% MLMG::mgVcycle_up::5 128 0.01346 0.01438 0.024 0.16% MLMG::computeMLResidual() 18 0.01118 0.01771 0.02349 0.16% Castro::save_data_for_retry() 3 0.01861 0.01884 0.01941 0.13% amrex::Add() 128 0.01426 0.01452 0.0153 0.10% FabArrayBase::getCPC() 2005 0.009622 0.0111 0.01522 0.10% FabArray::setDomainBndry() 147 0.0004747 0.007613 0.01483 0.10% Gravity::test_composite_phi() 1 0.01471 0.01477 0.01479 0.10% FabArrayBase::getFB() 6380 0.01038 0.01133 0.01468 0.10% FabArray::norminf() 244 0.01082 0.01139 0.01409 0.09% MLLinOp::define() 19 0.004368 0.008008 0.01174 0.08% MLLinOp::defineGrids() 19 0.004284 0.007919 0.01165 0.08% Other 2502 0.07964 0.09737 0.1417 0.95% -------------------------------------------------------------------------------------------- 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 10 B 11 B 13 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- Cpu Memory Usage: ------------------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max ------------------------------------------------------------------------------------------------------------------------------- Castro::save_data_for_retry() 416 416 1570 KiB 1570 KiB 1571 KiB 2297 KiB 2297 KiB 2297 KiB Castro::construct_ctu_hydro_source() 6400 6400 398 KiB 402 KiB 413 KiB 1156 KiB 1156 KiB 1156 KiB StateData::define() 192 192 1039 KiB 1039 KiB 1040 KiB 1148 KiB 1148 KiB 1148 KiB Castro::initMFs() 256 256 989 KiB 989 KiB 989 KiB 1093 KiB 1093 KiB 1093 KiB Castro::swap_state_time_levels() 160 160 664 KiB 664 KiB 664 KiB 787 KiB 787 KiB 787 KiB Amr::writePlotFile() 64 64 4125 B 4230 B 4884 B 736 KiB 736 KiB 736 KiB Castro::reset_internal_energy(MultiFab) 5120 5120 49 KiB 49 KiB 50 KiB 640 KiB 640 KiB 640 KiB FillPatchIterator::Initialize 4704 4704 12 KiB 14 KiB 15 KiB 500 KiB 500 KiB 500 KiB Castro::initialize_do_advance() 640 640 414 KiB 414 KiB 415 KiB 500 KiB 500 KiB 500 KiB FabArray::ParallelCopy_nowait() 12704 12704 2222 B 4413 B 5700 B 215 KiB 285 KiB 360 KiB AmrLevel::derive() 1152 1152 1280 B 1298 B 1359 B 336 KiB 336 KiB 336 KiB Castro::enforce_min_density() 2528 2528 17 KiB 17 KiB 17 KiB 320 KiB 320 KiB 320 KiB Castro::initialize_advance() 32 32 131 KiB 131 KiB 131 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 9736 B 10 KiB 11 KiB 102 KiB 108 KiB 126 KiB Castro::construct_old_source() 640 640 244 B 245 B 249 B 112 KiB 112 KiB 112 KiB Castro::construct_new_source() 544 544 409 B 410 B 417 B 112 KiB 112 KiB 112 KiB Amr::writeSmallPlotFile() 64 64 19 B 20 B 36 B 80 KiB 80 KiB 80 KiB Castro::Castro() 128 128 33 KiB 33 KiB 33 KiB 74 KiB 74 KiB 74 KiB Gravity::test_composite_phi() 96 96 52 B 52 B 52 B 52 KiB 52 KiB 52 KiB Gravity::update_max_rhs() 192 192 3 B 7 B 8 B 49 KiB 49 KiB 49 KiB Gravity::get_new_grav_vector() 576 576 108 B 123 B 140 B 45 KiB 45 KiB 45 KiB Gravity::get_old_grav_vector() 640 640 144 B 255 B 454 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 231 B 236 B 244 B 17 KiB 17 KiB 17 KiB Castro::evaluate_source_change() 6272 6272 1071 B 1082 B 1108 B 16 KiB 16 KiB 16 KiB Gravity::actual_multilevel_solve() 32 32 119 B 119 B 119 B 16 KiB 16 KiB 16 KiB Gravity::solve_for_phi() 544 544 1688 B 1762 B 1791 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 440 B 492 B 683 B 4608 B 5151 B 7152 B BndryData::define() 4864 4864 458 B 459 B 460 B 4864 B 4864 B 4864 B FillBoundary_nowait() 91296 91296 26 B 36 B 49 B 2592 B 3392 B 4192 B AmrMesh::MakeNewGrids() 128 128 0 B 19 B 30 B 2916 B 2916 B 2916 B OwnerMask() 64 64 0 B 0 B 0 B 2592 B 2592 B 2592 B MLMG::addInterpCorrection() 5632 5632 1 B 16 B 26 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 6 B 8 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-14-g65bf80d763f6) finalized