MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (23.06-14-g7c7a8fb1fff5) initialized Starting run at 07:42:45 UTC on 2023-06-09. Successfully read inputs file ... Castro git describe: 23.06-16-g0603698ef AMReX git describe: 23.06-14-g7c7a8fb1f Microphysics git describe: 23.06-8-gc1297423 reading extern runtime parameters ... 13 Species: He4 C12 O16 Ne20 Mg24 Si28 S32 Ar36 Ca40 Ti44 Cr48 Fe52 Ni56 Successfully read inputs file ... 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 Warning: the following castro.* parameters are ignored castro.dtnuc_mode = 1 Initializing the data at level 0 Done initializing the level 0 data grid_places() time: 0.137466332 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.000998148 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010475536 Gravity::fill_multipole_BCs() time = 0.012466025 ... norm of composite residual at level 0 1.193955541e-11 TIME= 0 MASS = 2.549996568e+33 TIME= 0 XMOM = 0 TIME= 0 YMOM = 0 TIME= 0 ZMOM = 0 TIME= 0 ANG MOM X = 0 TIME= 0 ANG MOM Y = 0 TIME= 0 ANG MOM Z = -3.115378115e+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 = -4.854673355e-07 TIME= 0 CENTER OF MASS Y-VEL = 0 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.050522103 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004429539 seconds [Level 0 step 1] ADVANCE with dt = 0.25  Beginning subcycle 1 starting at time 0 with dt = 0.25 Estimated number of subcycles remaining: 1 Beginning SDC iteration 1 of 2. 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.498822529e+28 zmom added: 0 eint added: 0 ener added: 6.931398274e+47 Castro::do_old_sources() time = 0.065531655 ... Entering construct_ctu_hydro_source() >>> RESETTING NEG. DENSITY AT 0, 62, 0 >>> FROM -9.947567811e-05 TO 1e-05 >>> IN GRID ((0,32) (31,63) (0,0)) >>> RESETTING NEG. DENSITY AT 0, 63, 0 >>> FROM -1.942940124 TO 1e-05 >>> IN GRID ((0,32) (31,63) (0,0)) >>> RESETTING NEG. DENSITY AT 1, 63, 0 >>> FROM -1.164962232 TO 1e-05 >>> IN GRID ((0,32) (31,63) (0,0)) >>> RESETTING NEG. DENSITY AT 2, 63, 0 >>> FROM -0.767071355 TO 1e-05 >>> IN GRID ((0,32) (31,63) (0,0)) >>> RESETTING NEG. DENSITY AT 3, 63, 0 >>> FROM -0.08878802337 TO 1e-05 >>> IN GRID ((0,32) (31,63) (0,0)) >>> RESETTING NEG. DENSITY AT 4, 63, 0 >>> FROM -0.01174181879 TO 1e-05 >>> IN GRID ((0,32) (31,63) (0,0)) >>> RESETTING NEG. DENSITY AT 5, 63, 0 >>> FROM -0.0008017363671 TO 1e-05 >>> IN GRID ((0,32) (31,63) (0,0)) >>> RESETTING NEG. DENSITY AT 0, 192, 0 >>> FROM -1.942940124 TO 1e-05 >>> IN GRID ((0,192) (31,223) (0,0)) >>> RESETTING NEG. DENSITY AT 1, 192, 0 >>> FROM -1.164962232 TO 1e-05 >>> IN GRID ((0,192) (31,223) (0,0)) >>> RESETTING NEG. DENSITY AT 2, 192, 0 >>> FROM -0.767071355 TO 1e-05 >>> IN GRID ((0,192) (31,223) (0,0)) >>> RESETTING NEG. DENSITY AT 3, 192, 0 >>> FROM -0.08878802337 TO 1e-05 >>> IN GRID ((0,192) (31,223) (0,0)) >>> RESETTING NEG. DENSITY AT 4, 192, 0 >>> FROM -0.01174181879 TO 1e-05 >>> IN GRID ((0,192) (31,223) (0,0)) >>> RESETTING NEG. DENSITY AT 5, 192, 0 >>> FROM -0.0008017363671 TO 1e-05 >>> IN GRID ((0,192) (31,223) (0,0)) >>> RESETTING NEG. DENSITY AT 0, 193, 0 >>> FROM -9.947567811e-05 TO 1e-05 >>> IN GRID ((0,192) (31,223) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 174, 0 >>> FROM -8.076810125 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 175, 0 >>> FROM -15.15828486 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 177, 0 >>> FROM -108.2035501 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 178, 0 >>> FROM -61.58489855 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 18, 179, 0 >>> FROM -80.16477042 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 18, 180, 0 >>> FROM -160.2179808 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 17, 181, 0 >>> FROM -552.7694089 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 16, 183, 0 >>> FROM -415.2911996 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 15, 184, 0 >>> FROM -764.1508347 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 14, 185, 0 >>> FROM -908.1731819 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 13, 186, 0 >>> FROM -859.0333271 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 12, 187, 0 >>> FROM -606.7538937 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 10, 188, 0 >>> FROM -907.0826581 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 11, 188, 0 >>> FROM -241.113047 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 8, 189, 0 >>> FROM -820.8014889 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 9, 189, 0 >>> FROM -456.4708539 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 10, 189, 0 >>> FROM -6.239071149 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 5, 190, 0 >>> FROM -625.0368072 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 6, 190, 0 >>> FROM -691.1269281 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 7, 190, 0 >>> FROM -318.5467376 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 8, 190, 0 >>> FROM -11.08517908 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 0, 191, 0 >>> FROM -618.7284155 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 1, 191, 0 >>> FROM -540.6137805 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 2, 191, 0 >>> FROM -443.905673 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 3, 191, 0 >>> FROM -302.5298758 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 4, 191, 0 >>> FROM -167.9940382 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 5, 191, 0 >>> FROM -55.18215035 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 6, 191, 0 >>> FROM -1.089487306 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 0, 64, 0 >>> FROM -618.7284155 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 1, 64, 0 >>> FROM -540.6137805 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 2, 64, 0 >>> FROM -443.905673 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 3, 64, 0 >>> FROM -302.5298758 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 4, 64, 0 >>> FROM -167.9940382 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 5, 64, 0 >>> FROM -55.18215035 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 6, 64, 0 >>> FROM -1.089487306 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 5, 65, 0 >>> FROM -625.0368072 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 6, 65, 0 >>> FROM -691.1269281 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 7, 65, 0 >>> FROM -318.5467376 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 8, 65, 0 >>> FROM -11.08517908 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 8, 66, 0 >>> FROM -820.8014889 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 9, 66, 0 >>> FROM -456.4708539 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 10, 66, 0 >>> FROM -6.239071149 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 10, 67, 0 >>> FROM -907.0826581 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 11, 67, 0 >>> FROM -241.113047 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 12, 68, 0 >>> FROM -606.7538937 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 13, 69, 0 >>> FROM -859.0333271 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 14, 70, 0 >>> FROM -908.1731819 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 15, 71, 0 >>> FROM -764.1508347 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 16, 72, 0 >>> FROM -415.2911996 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 17, 74, 0 >>> FROM -552.7694089 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 18, 75, 0 >>> FROM -160.2179808 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 18, 76, 0 >>> FROM -80.16477042 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 77, 0 >>> FROM -61.58489855 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 78, 0 >>> FROM -108.2035501 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 80, 0 >>> FROM -15.15828486 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 81, 0 >>> FROM -8.076810125 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) Contributions to the state from negative density resets: mass added: -8.63794634e+28 xmom added: -6.442755578e+36 ymom added: -5.625519073e+23 zmom added: 0 eint added: -4.036942199e+44 ener added: -6.630395009e+43 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 4.399017972e+48 ener added: 4.595479623e+45 ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.353049772 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.010474424 Gravity::solve_for_phi() time = 0.0795973 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 4.016301068e+46 ener added: 5.32095639e+52 Contributions to the state from new-time sources: mass added: 0 xmom added: 9.681452345e+35 ymom added: -2.452736705e+28 zmom added: 0 eint added: 0 ener added: 1.024838906e+46 Castro::do_new_sources() time = 0.065764572 Ending SDC iteration 1 of 2. Beginning SDC iteration 2 of 2. 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.498822529e+28 zmom added: 0 eint added: 0 ener added: 6.931398274e+47 Castro::do_old_sources() time = 0.066841654 ... Entering construct_ctu_hydro_source() >>> RESETTING NEG. DENSITY AT 0, 62, 0 >>> FROM -7.355747832e-05 TO 1e-05 >>> IN GRID ((0,32) (31,63) (0,0)) >>> RESETTING NEG. DENSITY AT 0, 63, 0 >>> FROM -0.5625596674 TO 1e-05 >>> IN GRID ((0,32) (31,63) (0,0)) >>> RESETTING NEG. DENSITY AT 1, 63, 0 >>> FROM -0.3360676128 TO 1e-05 >>> IN GRID ((0,32) (31,63) (0,0)) >>> RESETTING NEG. DENSITY AT 2, 63, 0 >>> FROM -0.4666974183 TO 1e-05 >>> IN GRID ((0,32) (31,63) (0,0)) >>> RESETTING NEG. DENSITY AT 3, 63, 0 >>> FROM -0.02178000789 TO 1e-05 >>> RESETTING NEG. DENSITY AT 0, 148, 0 >>> FROM >>> IN GRID ((0,32) (31,63) (0,0)) >>> RESETTING NEG. DENSITY AT 4, 63, 0 >>> FROM -0.00173451199 TO 1e-05 -0.0002386989058 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 1, >>> IN GRID ((0,32) (31,63) (0,0)) >>> RESETTING NEG. DENSITY AT 0, 192, 0 >>> FROM -0.5625596674 TO 1e-05 148, 0 >>> FROM -0.0002442862726 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT >>> IN GRID ((0,192) (31,223) (0,0)) >>> RESETTING NEG. DENSITY AT 1, 192, 0 >>> FROM -0.3360676128 TO 1e-05 >>> IN GRID 2, 148, 0 >>> FROM -0.0002433077831 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT ((0,192) (31,223) (0,0)) >>> RESETTING NEG. DENSITY AT 2, 192, 0 >>> FROM -0.4666974183 TO 1e-05 >>> IN GRID 3, 148, 0 >>> FROM -0.0002173903377 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT ((0,192) (31,223) (0,0)) >>> RESETTING NEG. DENSITY AT 3, 192, 0 >>> FROM -0.02178000789 TO 1e-05 >>> IN GRID ((0,192) 4, 148, 0 >>> FROM -0.0001938467994 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 1, 149, 0 >>> FROM -0.1835531121 TO 1e-05 >>> IN GRID ((0(31,223) (0,0)) >>> RESETTING NEG. DENSITY AT 4, 192, 0 >>> FROM -0.00173451199 TO 1e-05 >>> IN GRID ((0,192) (31,223,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 2, 149, 0 >>> FROM -0.01182645204 TO 1e-05 >>> IN GRID ((0,) (0,0)) >>> RESETTING NEG. DENSITY AT 0, 193, 0 >>> FROM -7.355747832e-05 TO 1e-05 >>> IN GRID ((0,192) (31,223128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 5, 149, 0 >>> FROM -0.000127100675 TO 1e-05 >>> IN GRID ((0) (0,0)) ,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 6, 149, 0 >>> FROM -0.0001998888189 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 7, 149, 0 >>> FROM -0.0002963218592 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 9, 149, 0 >>> FROM -1.218500506e-05 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 0, 150, 0 >>> FROM -2891.682129 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 2, 150, 0 >>> FROM -1947.592872 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 4, 150, 0 >>> FROM -15.63589649 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 6, 150, 0 >>> FROM -0.01430210716 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 7, 150, 0 >>> FROM -0.005826589474 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 8, 150, 0 >>> FROM -0.0003869262728 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 9, 150, 0 >>> FROM -0.0002265300772 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 11, 150, 0 >>> FROM -4.094001771e-05 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 5, 151, 0 >>> FROM -11291.48442 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 8, 151, 0 >>> FROM -0.02280560839 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 9, 151, 0 >>> FROM -0.02663416158 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 10, 151, 0 >>> FROM -0.0003820609327 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 11, 151, 0 >>> FROM -0.0001918215002 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 13, 151, 0 >>> FROM -2.414090548e-05 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 8, 152, 0 >>> FROM -4718.917246 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 11, 152, 0 >>> FROM -0.03043351756 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 12, 152, 0 >>> FROM -0.0004128595986 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 13, 152, 0 >>> FROM -5.585694285e-05 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 10, 153, 0 >>> FROM -5397.328835 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 12, 153, 0 >>> FROM -0.1039031694 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 13, 153, 0 >>> FROM -0.001435731739 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 14, 153, 0 >>> FROM -0.0004780182213 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 14, 154, 0 >>> FROM -0.009816209719 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 15, 154, 0 >>> FROM -0.0005032745505 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 13, 155, 0 >>> FROM -1809.897411 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 15, 155, 0 >>> FROM -0.01850143321 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 16, 155, 0 >>> FROM -0.0005249404306 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 16, 156, 0 >>> FROM -0.01364865427 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 17, 156, 0 >>> FROM -0.0004610255684 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 17, 157, 0 >>> FROM -0.003298985043 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 18, 157, 0 >>> FROM -0.0001129243674 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 157, 0 >>> FROM -5.302322498e-05 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 17, 158, 0 >>> FROM -0.1301738217 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 18, 158, 0 >>> FROM -0.0004334535201 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 17, 159, 0 >>> FROM -0.2777479179 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 18, 159, 0 >>> FROM -0.01012067404 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 159, 0 >>> FROM -0.0001027182743 TO 1e-05 >>> IN GRID ((0,128) (31,159) (0,0)) >>> RESETTING NEG. DENSITY AT 18, 160, 0 >>> FROM -0.1400920109 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 160, 0 >>> FROM -0.001828334082 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 160, 0 >>> FROM -0.0001790829998 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING SMALL DENSITY AT 21, 160, 0 >>> FROM 7.190266466e-06 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 161, 0 >>> FROM -7.375343036e-05 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 162, 0 >>> FROM -0.0004830068255 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 21, 162, 0 >>> FROM -0.0001938110725 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 21, 163, 0 >>> FROM -0.0001577330399 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING SMALL DENSITY AT 21, 164, 0 >>> FROM 4.485106849e-06 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 22, 164, 0 >>> FROM -7.191625508e-05 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 165, 0 >>> FROM -3.221894241 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 21, 165, 0 >>> FROM -0.0001055752926 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 22, 165, 0 >>> FROM -1.893312962e-05 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 21, 166, 0 >>> FROM -0.001425551306 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 22, 167, 0 >>> FROM -0.0002667408166 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 22, 168, 0 >>> FROM -0.0001725195728 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 21, 169, 0 >>> FROM -0.06741775584 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING SMALL DENSITY AT 22, 172, 0 >>> FROM 3.357381342e-06 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 174, 0 >>> FROM -2.041828578 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 175, 0 >>> FROM -21.90201134 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 177, 0 >>> FROM -11.96023318 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 178, 0 >>> FROM -74.37700663 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 18, 180, 0 >>> FROM -167.30368 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 17, 181, 0 >>> FROM -252.4011455 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 182, 0 >>> FROM -2.009505659e-05 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 16, 183, 0 >>> FROM -380.0962691 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 15, 184, 0 >>> FROM -660.7610919 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 14, 185, 0 >>> FROM -737.8680307 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 13, 186, 0 >>> FROM -741.903731 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 12, 187, 0 >>> FROM -516.4803492 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 10, 188, 0 >>> FROM -480.1358683 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 11, 188, 0 >>> FROM -205.4944334 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 8, 189, 0 >>> FROM -272.0201172 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 9, 189, 0 >>> FROM -369.0327843 TO 1e-05 >>> IN GRID ((0,160) (31, >>> RESETTING NEG. DENSITY AT 17, 96, 0 >>> FROM -0.2777479187 TO 1e-05 >>> IN GRID ((0,96) 191) (0,0)) >>> RESETTING NEG. DENSITY AT 6, 190, 0 >>> FROM -327.5782885 TO 1e-05 >>> IN GRID ((0,160) (31,191)(31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 18, 96, 0 >>> FROM -0.01012067404 TO 1e-05 >>> IN GRID ((0 (0,0)) >>> RESETTING NEG. DENSITY AT 7, 190, 0 >>> FROM -240.7258411 TO 1e-05 >>> IN GRID ((0,160) (31,191) ,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 96, 0 >>> FROM -0.0001027182743 TO 1e-05 >>> IN GRID ((0,(0,0)) >>> RESETTING NEG. DENSITY AT 0, 191, 0 >>> FROM -180.9630155 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 17, 97, 0 >>> FROM -0.1301738217 TO 1e-05 >>> IN GRID ((0,96) (31,) >>> RESETTING NEG. DENSITY AT 1, 191, 0 >>> FROM -152.908186 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) 127) (0,0)) >>> RESETTING NEG. DENSITY AT 18, 97, 0 >>> FROM -0.0004334535201 TO 1e-05 >>> IN GRID ((0,96) (31,127) >>> RESETTING NEG. DENSITY AT 2, 191, 0 >>> FROM -154.1842501 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) (0,0)) >>> RESETTING NEG. DENSITY AT 17, 98, 0 >>> FROM -0.003298985043 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0, >>> RESETTING NEG. DENSITY AT 3, 191, 0 >>> FROM -125.9846333 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) 0)) >>> RESETTING NEG. DENSITY AT 18, 98, 0 >>> FROM -0.0001129243674 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 4, 191, 0 >>> FROM -89.75187256 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 98, 0 >>> FROM -5.302322498e-05 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 5, 191, 0 >>> FROM -40.09380837 TO 1e-05 >>> IN GRID ((0,160) (31,191) (0,0)) >>> RESETTING NEG. DENSITY AT 16, 99, 0 >>> FROM -0.01364865427 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 17, 99, 0 >>> FROM -0.0004610255684 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 13, 100, 0 >>> FROM -1809.897411 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 15, 100, 0 >>> FROM -0.01850143321 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 16, 100, 0 >>> FROM -0.0005249404306 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 14, 101, 0 >>> FROM -0.009816209719 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 15, 101, 0 >>> FROM -0.0005032745505 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 10, 102, 0 >>> FROM -5397.328835 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 12, 102, 0 >>> FROM -0.1039031694 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 13, 102, 0 >>> FROM -0.001435731739 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 14, 102, 0 >>> FROM -0.0004780182213 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 8, 103, 0 >>> FROM -4718.917246 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 11, 103, 0 >>> FROM -0.03043351756 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 12, 103, 0 >>> FROM -0.0004128595986 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 13, 103, 0 >>> FROM -5.585694285e-05 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 5, 104, 0 >>> FROM -11291.48442 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 8, 104, 0 >>> FROM -0.02280560839 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 9, 104, 0 >>> FROM -0.02663416158 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 10, 104, 0 >>> FROM -0.0003820609327 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 11, 104, 0 >>> FROM -0.0001918215002 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 13, 104, 0 >>> FROM -2.414090548e-05 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 0, 105, 0 >>> FROM -2891.682129 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 2, 105, 0 >>> FROM -1947.592872 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 4, 105, 0 >>> FROM -15.63589649 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 6, 105, 0 >>> FROM -0.01430210716 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 7, 105, 0 >>> FROM -0.005826589474 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 8, 105, 0 >>> FROM -0.0003869262728 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 9, 105, 0 >>> FROM -0.0002265300772 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 11, 105, 0 >>> FROM -4.094001771e-05 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 1, 106, 0 >>> FROM -0.1835531122 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 2, 106, 0 >>> FROM -0.01182645204 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 5, 106, 0 >>> FROM -0.000127100675 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 6, 106, 0 >>> FROM -0.0001998888189 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 7, 106, 0 >>> FROM -0.0002963218592 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 9, 106, 0 >>> FROM -1.218500506e-05 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 0, 107, 0 >>> FROM -0.0002386989058 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 1, 107, 0 >>> FROM -0.0002442862726 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 2, 107, 0 >>> FROM -0.0002433077831 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 3, 107, 0 >>> FROM -0.0002173903378 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 4, 107, 0 >>> FROM -0.0001938467994 TO 1e-05 >>> IN GRID ((0,96) (31,127) (0,0)) >>> RESETTING NEG. DENSITY AT 0, 64, 0 >>> FROM -180.9630155 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 1, 64, 0 >>> FROM -152.908186 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 2, 64, 0 >>> FROM -154.1842501 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 3, 64, 0 >>> FROM -125.9846333 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 4, 64, 0 >>> FROM -89.75187256 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 5, 64, 0 >>> FROM -40.09380837 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 6, 65, 0 >>> FROM -327.5782885 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 7, 65, 0 >>> FROM -240.7258411 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 8, 66, 0 >>> FROM -272.0201172 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 9, 66, 0 >>> FROM -369.0327843 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 10, 67, 0 >>> FROM -480.1358683 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 11, 67, 0 >>> FROM -205.4944334 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 12, 68, 0 >>> FROM -516.4803492 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 13, 69, 0 >>> FROM -741.903731 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 14, 70, 0 >>> FROM -737.8680307 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 15, 71, 0 >>> FROM -660.7610919 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 16, 72, 0 >>> FROM -380.0962691 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 73, 0 >>> FROM -2.009505659e-05 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 17, 74, 0 >>> FROM -252.4011455 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 18, 75, 0 >>> FROM -167.30368 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 77, 0 >>> FROM -74.37700663 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 78, 0 >>> FROM -11.96023318 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 80, 0 >>> FROM -21.90201134 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 81, 0 >>> FROM -2.041828578 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING SMALL DENSITY AT 22, 83, 0 >>> FROM 3.357381341e-06 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 21, 86, 0 >>> FROM -0.06741775584 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 22, 87, 0 >>> FROM -0.0001725195728 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 22, 88, 0 >>> FROM -0.0002667408166 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 21, 89, 0 >>> FROM -0.001425551306 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 90, 0 >>> FROM -3.221894241 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 21, 90, 0 >>> FROM -0.0001055752926 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 22, 90, 0 >>> FROM -1.893312962e-05 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING SMALL DENSITY AT 21, 91, 0 >>> FROM 4.485106854e-06 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 22, 91, 0 >>> FROM -7.191625508e-05 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 21, 92, 0 >>> FROM -0.0001577330399 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 93, 0 >>> FROM -0.0004830068255 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 21, 93, 0 >>> FROM -0.0001938110725 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 94, 0 >>> FROM -7.375343036e-05 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 18, 95, 0 >>> FROM -0.1400920109 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 19, 95, 0 >>> FROM -0.001828334082 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING NEG. DENSITY AT 20, 95, 0 >>> FROM -0.0001790829998 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) >>> RESETTING SMALL DENSITY AT 21, 95, 0 >>> FROM 7.190266465e-06 TO 1e-05 >>> IN GRID ((0,64) (31,95) (0,0)) Contributions to the state from negative density resets: mass added: -2.107823144e+29 xmom added: 2.556641392e+38 ymom added: -3.179499255e+26 zmom added: 0 eint added: -1.86234351e+49 ener added: -4.899585579e+49 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 3.658755439e+49 ener added: 1.443380131e+49 ... Leaving construct_ctu_hydro_source() Castro::construct_ctu_hydro_source() time = 0.338890218 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.010537428 Gravity::solve_for_phi() time = 0.079900007 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -3.392134076e+47 ener added: 3.848343937e+52 Contributions to the state from new-time sources: mass added: 0 xmom added: 1.609831707e+37 ymom added: -2.494195305e+28 zmom added: 0 eint added: 0 ener added: -1.925807569e+46 Castro::do_new_sources() time = 0.079268872 Ending SDC iteration 2 of 2.  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.024153703 [Level 0 step 1] Advanced 32768 cells TIME= 0.25 MASS = 2.55020735e+33 TIME= 0.25 XMOM = -3.731657684e+38 TIME= 0.25 YMOM = 4.061990754e+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 = -6.697647564e+38 TIME= 0.25 RHO*e = 1.909328145e+50 TIME= 0.25 RHO*K = 3.852415681e+52 TIME= 0.25 RHO*E = 3.871508969e+52 TIME= 0.25 RHO*PHI = -5.448837462e+50 TIME= 0.25 TOTAL ENERGY= 3.844264781e+52 TIME= 0.25 CENTER OF MASS X-LOC = 316437217.8 TIME= 0.25 CENTER OF MASS X-VEL = -146327.6186 TIME= 0.25 CENTER OF MASS Y-LOC = 3.033920063e-06 TIME= 0.25 CENTER OF MASS Y-VEL = 1.592808033e-05 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 5412040692 TIME= 0.25 MAXIMUM DENSITY = 4532127.98 TIME= 0.25 MAXIMUM T_S / T_E = 0 [STEP 1] Coarse TimeStep time: 1.427646413 [STEP 1] FAB kilobyte spread across MPI nodes: [5036 ... 5036] STEP = 1 TIME = 0.25 DT = 0.25 PLOTFILE: file = wdmerger_retry-SDC_plt00001 Write plotfile time = 0.05278937 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004413893 seconds Ending run at 07:42:49 UTC on 2023-06-09. Run time = 3.895719875 Run time without initialization = 1.48628468 Average number of zones advanced per microsecond: 0.022 Average number of zones advanced per microsecond per rank: 0.001 TinyProfiler total time across processes [min...avg...max]: 3.896 ... 3.896 ... 3.896 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Amr::InitAmr() 1 1.5303 1.5304 1.5310 39.30% Castro::construct_ctu_hydro_source() 2 0.5344 0.5476 0.5783 14.84% Castro::initMFs() 1 0.3623 0.3836 0.4334 11.12% FabArray::ParallelCopy_finish() 277 0.01151 0.1438 0.2222 5.70% Castro::computeTemp() 14 0.1595 0.1692 0.2067 5.31% amrex::Copy() 259 0.1605 0.1625 0.1667 4.28% Castro::initData() 1 0.003303 0.02316 0.1313 3.37% Castro::reset_internal_energy(MultiFab) 14 0.004507 0.02541 0.1212 3.11% MultiFab::AddProduct() 568 0.09359 0.09463 0.09893 2.54% Castro::normalize_species() 13 0.08472 0.08617 0.09796 2.51% Castro::reset_internal_energy(Fab) 28 0.0885 0.08945 0.09515 2.44% MLPoisson::Fsmooth() 1104 0.05642 0.06085 0.07298 1.87% check_for_negative_density() 2 0.002393 0.04003 0.06805 1.75% MultiFab::Subtract() 27 0.05728 0.05859 0.06669 1.71% Castro::enforce_min_density() 13 0.00273 0.01365 0.0647 1.66% FabArray::setVal() 852 0.05973 0.0611 0.06363 1.63% Castro::derive() 18 0.04199 0.04276 0.04812 1.24% Gravity::fill_multipole_BCs() 4 0.04081 0.04266 0.04349 1.12% Castro::wd_update() 2 0.004966 0.03102 0.03648 0.94% FabArray::sum() 568 0.03077 0.03121 0.03252 0.83% VisMF::Write(FabArray) 4 0.01239 0.01367 0.02872 0.74% FabArray::ParallelCopy_nowait() 277 0.02135 0.02302 0.02585 0.66% FillBoundary_finish() 1147 0.01615 0.02042 0.02447 0.63% FillBoundary_nowait() 1147 0.01334 0.01779 0.02269 0.58% Amr::writePlotFile() 2 0.007021 0.01931 0.02183 0.56% MLPoisson::Fapply() 178 0.01732 0.0182 0.02081 0.53% FabArray::Saxpy() 27 0.007543 0.008384 0.01972 0.51% MLCellLinOp::applyBC() 1285 0.008619 0.01136 0.01886 0.48% Castro::do_old_sources() 2 8.489e-05 0.002013 0.01169 0.30% FabArray::Xpay() 171 0.009017 0.009563 0.01089 0.28% Castro::volWgtSum() 46 0.007349 0.007532 0.009095 0.23% Castro::construct_new_gravity_source() 2 0.008027 0.008163 0.00833 0.21% MLMG::addInterpCorrection() 138 0.006114 0.006541 0.007722 0.20% Castro::construct_new_source() 12 0.001821 0.002119 0.00606 0.16% Amr::writeSmallPlotFile() 2 0.003995 0.005683 0.005965 0.15% MLMG::ResNormInf() 26 0.0009546 0.004253 0.005921 0.15% amrex::average_down 138 0.004507 0.004875 0.005894 0.15% MultiFab::contains_nan() 4 0.004798 0.004908 0.005487 0.14% Castro::volInBoundary() 14 0.003451 0.003739 0.0054 0.14% Castro::construct_old_gravity_source() 2 0.004613 0.004678 0.005363 0.14% Castro::do_new_sources() 2 4.695e-05 0.0009185 0.005161 0.13% Castro::sum_integrated_quantities() 2 0.002456 0.004047 0.004365 0.11% MultiFab::max() 2 0.0001529 0.002632 0.004098 0.11% MLCGSolver::bicgstab 1 0 0.0002449 0.003918 0.10% MLCellLinOp::defineAuxData() 4 0.003071 0.003228 0.003654 0.09% StateData::FillBoundary(geom) 16 0 0.00189 0.003629 0.09% amrex::Add() 24 0.002572 0.002636 0.003168 0.08% FabArray::FillBoundary() 1147 0.002501 0.002597 0.003125 0.08% FillPatchIterator::Initialize 45 0.002587 0.002685 0.002813 0.07% MLCellLinOp::compGrad() 3 0.0009261 0.001052 0.002588 0.07% FabArray::norminf() 43 0.001994 0.002124 0.002525 0.06% MLCellLinOp::smooth() 552 0.002115 0.002201 0.002372 0.06% FabArrayBase::getFB() 1147 0.001844 0.001923 0.002368 0.06% Castro::evaluate_source_change() 31 0.002171 0.002205 0.002343 0.06% Castro::buildMetrics() 1 0.000569 0.0007584 0.001871 0.05% MLMG::prepareForSolve() 3 0.001716 0.001746 0.00183 0.05% MLLinOp::makeSubCommunicator() 4 0.0001683 0.0009174 0.001789 0.05% StateData::define() 6 0.001345 0.001431 0.001722 0.04% Gravity::actual_solve_with_mlmg() 4 0.001356 0.001435 0.001711 0.04% FabArrayBase::CPC::define() 126 0.0008177 0.001147 0.001629 0.04% FabArrayBase::getCPC() 373 0.00135 0.001391 0.00148 0.04% main() 1 0.000887 0.001166 0.001415 0.04% FabArray::setDomainBndry() 45 0.00013 0.0007531 0.001329 0.03% FabArray::mult() 11 0.001167 0.001185 0.00122 0.03% Castro::locWgtSum() 6 0.001159 0.001174 0.001189 0.03% amrex::Dot() 24 0 6.977e-05 0.001116 0.03% Gravity::get_new_grav_vector() 3 0.001002 0.001046 0.001104 0.03% MLCGSolver::ParallelAllReduce 32 0 6.872e-05 0.0011 0.03% TagBoxArray::collate() 2 0.0001247 0.000509 0.001047 0.03% Castro::construct_old_source() 12 0.0009661 0.0009822 0.001005 0.03% MLCellLinOp::apply() 178 0.0003835 0.0004296 0.0009407 0.02% DistributionMapping::LeastUsedCPUs() 1 2.779e-05 0.0008138 0.0009264 0.02% MLMG::solve() 3 0.0001171 0.0004134 0.0009029 0.02% Castro::Castro() 1 0.0004731 0.0006243 0.0007802 0.02% FabArray::ParallelCopy() 277 0.0006944 0.0007244 0.0007723 0.02% Castro::enforce_consistent_e() 1 0.000358 0.0003934 0.0007699 0.02% BndryData::define() 4 0.0006503 0.0006617 0.0007057 0.02% Gravity::get_old_grav_vector() 2 0.0006649 0.0006731 0.000704 0.02% FabArray::LinComb() 12 0 4.284e-05 0.0006855 0.02% Amr::coarseTimeStep() 1 0.0005777 0.0006021 0.0006065 0.02% MLPoisson::normalize() 13 0 3.583e-05 0.0005732 0.01% AmrMesh::MakeNewGrids() 2 0.0005078 0.0005244 0.0005428 0.01% Castro::apply_tagging_restrictions() 2 0.0004058 0.0004655 0.0005216 0.01% Gravity::applyMetricTerms() 2 0.0004716 0.0004752 0.0004849 0.01% Castro::initialize_do_advance() 2 0.0004159 0.0004336 0.0004715 0.01% MLLinOp::defineGrids() 4 0.0003826 0.0004183 0.0004569 0.01% MLCellLinOp::prepareForSolve() 4 0.0003545 0.0003848 0.0004503 0.01% MLCellLinOp::correctionResidual() 139 0.000346 0.0003626 0.0004426 0.01% Castro::swap_state_time_levels() 1 0.0003753 0.0003886 0.0004248 0.01% Castro::volProductSum() 2 0.0003963 0.0004079 0.0004189 0.01% FabArrayBase::FB::FB() 14 0.0001612 0.0002682 0.0003913 0.01% Castro::problem_post_timestep() 2 1.149e-05 3.733e-05 0.0003821 0.01% MLCellLinOp::defineBC() 4 0.0003275 0.0003432 0.0003753 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 128 0.000274 0.0003242 0.0003727 0.01% MLMG::mgVcycle() 23 0.0002872 0.0003049 0.00034 0.01% MLMG:computeResOfCorrection() 138 0.0002574 0.0002651 0.0002762 0.01% MLMG::actualBottomSolve() 23 2.5e-05 3.716e-05 0.0002106 0.01% Amr::initSubcycle() 1 6.434e-05 8.194e-05 0.0002005 0.01% StateDataPhysBCFunct::() 45 0.000168 0.0001831 0.0001996 0.01% Gravity::solve_for_phi() 2 9.434e-05 0.0001259 0.0001979 0.01% MLCellLinOp::setLevelBC() 4 0.0001175 0.0001448 0.000182 0.00% Gravity::update_max_rhs() 2 0.0001639 0.0001684 0.0001776 0.00% OwnerMask() 2 0.000136 0.0001551 0.0001713 0.00% MLMG::compResidual() 1 0.0001533 0.0001563 0.0001608 0.00% FillPatchSingleLevel 45 0.0001372 0.0001435 0.0001581 0.00% MLMG::mgVcycle_down::1 23 0.0001265 0.0001308 0.0001528 0.00% TagBoxArray::mapPRD 2 0.0001094 0.0001215 0.0001389 0.00% MLMG::mgVcycle_down::2 23 0.0001195 0.0001296 0.0001388 0.00% MLMG::mgVcycle_down::0 23 0.0001276 0.0001297 0.0001387 0.00% knapsack() 1 5.801e-05 7.292e-05 0.0001369 0.00% Amr::grid_places() 2 7.804e-05 9.569e-05 0.0001366 0.00% MLMG::mgVcycle_down::4 23 0.000107 0.0001118 0.0001352 0.00% Amr::defBaseLevel() 1 9.74e-05 0.0001075 0.0001346 0.00% MLMG::mgVcycle_down::3 23 0.00011 0.00012 0.0001343 0.00% MLMG::mgVcycle_down::5 23 0.0001101 0.0001131 0.0001283 0.00% Castro::initialize_advance() 1 0.0001011 0.0001099 0.0001201 0.00% Gravity::test_composite_phi() 1 9.254e-05 9.667e-05 0.000109 0.00% MLMG::mgVcycle_up::0 23 9.048e-05 9.384e-05 0.000108 0.00% Castro::subcycle_advance_ctu() 1 2.632e-05 3.372e-05 0.0001007 0.00% Castro::clean_state() 13 7.587e-05 8.168e-05 0.0001005 0.00% AmrLevel::AmrLevel(dm) 1 4.623e-05 7.937e-05 9.861e-05 0.00% Castro::apply_problem_tags() 2 9.465e-05 9.709e-05 9.853e-05 0.00% makeSFC 20 8.351e-05 9.004e-05 9.815e-05 0.00% MLMG::mgVcycle_up::2 23 7.409e-05 8.282e-05 9.795e-05 0.00% MLMG::mgVcycle_up::1 23 8.746e-05 9.151e-05 9.529e-05 0.00% MLMG::mgVcycle_up::3 23 7.384e-05 7.693e-05 9.146e-05 0.00% MLCellLinOp::solutionResidual() 27 8.269e-05 8.541e-05 8.971e-05 0.00% MLMG::mgVcycle_up::5 23 7.73e-05 7.912e-05 8.931e-05 0.00% Castro::do_advance_ctu() 2 6.981e-05 7.511e-05 8.645e-05 0.00% MLMG::mgVcycle_up::4 23 7.437e-05 7.647e-05 8.491e-05 0.00% MLMG::oneIter() 23 7.19e-05 7.675e-05 8.315e-05 0.00% MLLinOp::makeAgglomeratedDMap 4 7.166e-05 7.601e-05 8.037e-05 0.00% Amr::FinalizeInit() 1 8.949e-06 1.442e-05 6.917e-05 0.00% Gravity::solve_phi_with_mlmg() 4 4.845e-05 5.033e-05 5.699e-05 0.00% DistributionMapping::Distribute() 20 4.553e-05 4.844e-05 5.685e-05 0.00% MLMG::computeResidual() 23 4.464e-05 4.755e-05 5.581e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 4.932e-05 5.049e-05 5.156e-05 0.00% MLMG::mgVcycle_bottom 23 4.133e-05 4.268e-05 4.93e-05 0.00% Castro::expand_state() 2 4.181e-05 4.448e-05 4.884e-05 0.00% Castro::create_source_corrector() 2 3.345e-05 3.603e-05 4.105e-05 0.00% MLPoisson::define() 4 3.153e-05 3.382e-05 3.862e-05 0.00% Castro::construct_new_gravity() 2 1.854e-05 2.104e-05 3.755e-05 0.00% DistributionMapping::KnapSackDoIt() 1 1.692e-05 2.035e-05 3.59e-05 0.00% Amr::regrid() 1 1.153e-05 1.428e-05 3.531e-05 0.00% Castro::finalize_advance() 1 2.437e-05 2.656e-05 3.296e-05 0.00% Castro::finalize_do_advance() 2 2.433e-05 2.645e-05 3.243e-05 0.00% Gravity::actual_multilevel_solve() 1 2.744e-05 2.85e-05 3.066e-05 0.00% Castro::post_init() 1 1.228e-05 1.39e-05 3.064e-05 0.00% MLPoisson::prepareForSolve() 4 2.203e-05 2.343e-05 2.725e-05 0.00% MLLinOp::define() 4 1.517e-05 1.67e-05 2.665e-05 0.00% Castro::computeInitialDt() 2 7.233e-06 8.905e-06 2.644e-05 0.00% Castro::enforce_speed_limit() 13 1.802e-05 1.923e-05 2.221e-05 0.00% Amr::InitializeInit() 1 6.748e-06 9.66e-06 2.178e-05 0.00% Castro::post_timestep() 1 6.286e-06 8.934e-06 2.152e-05 0.00% Amr::timeStep() 1 1.135e-05 1.312e-05 1.903e-05 0.00% Castro::apply_source_to_state() 4 1.071e-05 1.196e-05 1.712e-05 0.00% Castro::post_regrid() 2 1.102e-05 1.301e-05 1.557e-05 0.00% Gravity::set_mass_offset() 2 4.602e-06 5.599e-06 1.34e-05 0.00% Castro::check_for_nan() 4 1.114e-05 1.185e-05 1.305e-05 0.00% MLMG::getGradSolution() 3 8.008e-06 8.785e-06 1.209e-05 0.00% Castro::construct_old_gravity() 2 9.182e-06 1.031e-05 1.123e-05 0.00% Amr::bldFineLevels() 1 7.201e-06 8.729e-06 1.071e-05 0.00% BoxList::parallelComplementIn() 4 7.055e-06 7.574e-06 1.017e-05 0.00% Amr::init() 1 6.682e-06 7.721e-06 9.751e-06 0.00% Castro::estTimeStep() 3 3.934e-06 4.733e-06 9.46e-06 0.00% Castro::gwstrain() 2 4.33e-06 5.535e-06 8.665e-06 0.00% BoxList::complementIn 4 6.127e-06 6.451e-06 8.347e-06 0.00% Castro::advance() 1 6.926e-06 7.522e-06 8.053e-06 0.00% MLMG::computeMLResidual() 3 7.268e-06 7.568e-06 8.046e-06 0.00% Castro::errorEst() 2 6.34e-06 6.613e-06 7.348e-06 0.00% MLMG::MLResNormInf() 3 5.921e-06 6.223e-06 6.658e-06 0.00% Gravity::swapTimeLevels() 1 5.546e-06 6.064e-06 6.567e-06 0.00% MLMG::MLRhsNormInf() 3 5.824e-06 6.017e-06 6.383e-06 0.00% Amr::initialInit() 1 3.01e-06 3.516e-06 5.161e-06 0.00% Castro::retry_advance_ctu() 1 3.699e-06 4.173e-06 4.944e-06 0.00% Castro::construct_old_rotation_source() 2 3.751e-06 3.967e-06 4.292e-06 0.00% Castro::construct_new_rotation_source() 2 3.222e-06 3.442e-06 4.04e-06 0.00% Castro::FluxRegCrseInit 1 1.768e-06 2.23e-06 3.183e-06 0.00% Castro::FluxRegFineAdd() 1 1.36e-06 1.944e-06 2.686e-06 0.00% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 3.895 3.896 3.896 100.00% Amr::InitAmr() 1 1.53 1.531 1.531 39.30% Amr::coarseTimeStep() 1 1.486 1.486 1.486 38.13% Amr::timeStep() 1 1.428 1.428 1.428 36.65% Castro::advance() 1 1.335 1.339 1.357 34.82% Castro::subcycle_advance_ctu() 1 1.285 1.289 1.306 33.53% Castro::do_advance_ctu() 2 1.285 1.289 1.306 33.53% Amr::init() 1 0.8775 0.8782 0.8783 22.54% Amr::initialInit() 1 0.8216 0.8223 0.8225 21.11% Castro::clean_state() 13 0.6149 0.641 0.7164 18.39% Castro::construct_ctu_hydro_source() 2 0.6301 0.6564 0.692 17.76% Amr::InitializeInit() 1 0.4519 0.5123 0.5874 15.08% Amr::defBaseLevel() 1 0.4519 0.5123 0.5874 15.08% Castro::computeTemp() 14 0.4249 0.4575 0.5527 14.19% Castro::Castro() 1 0.3678 0.3894 0.4404 11.31% Castro::initMFs() 1 0.3624 0.3838 0.4336 11.13% Castro::reset_internal_energy(MultiFab) 14 0.2646 0.2883 0.3815 9.79% Amr::FinalizeInit() 1 0.235 0.31 0.3705 9.51% Gravity::solve_phi_with_mlmg() 4 0.2801 0.2802 0.2803 7.20% FabArray::ParallelCopy() 277 0.05085 0.1833 0.2616 6.72% Gravity::actual_solve_with_mlmg() 4 0.2358 0.2368 0.2386 6.12% MLMG::solve() 3 0.2247 0.2254 0.2257 5.79% FabArray::ParallelCopy_finish() 277 0.01151 0.1438 0.2222 5.70% Castro::initData() 1 0.07938 0.1204 0.2089 5.36% Castro::enforce_min_density() 13 0.1329 0.1456 0.2041 5.24% MLMG::oneIter() 23 0.1966 0.1977 0.1997 5.13% MLMG::mgVcycle() 23 0.1939 0.195 0.1965 5.04% Amr::bldFineLevels() 1 0.04695 0.1227 0.1831 4.70% Castro::evaluate_source_change() 31 0.1587 0.1606 0.168 4.31% amrex::Copy() 259 0.1605 0.1625 0.1667 4.28% Castro::construct_new_gravity() 2 0.1637 0.1638 0.164 4.21% Gravity::solve_for_phi() 2 0.1594 0.1595 0.1595 4.09% Castro::initialize_do_advance() 2 0.1421 0.1471 0.1492 3.83% Castro::do_new_sources() 2 0.1163 0.1217 0.1451 3.72% Castro::post_init() 1 0.1407 0.1431 0.1446 3.71% Amr::grid_places() 2 0.002946 0.07794 0.1385 3.55% AmrMesh::MakeNewGrids() 2 0.002809 0.07785 0.1384 3.55% TagBoxArray::mapPRD 2 0.0007396 0.07607 0.137 3.52% Castro::do_old_sources() 2 0.1183 0.129 0.1322 3.39% MLCellLinOp::smooth() 552 0.09807 0.106 0.1274 3.27% Gravity::multilevel_solve_for_new_phi() 1 0.1058 0.106 0.108 2.77% Gravity::actual_multilevel_solve() 1 0.1058 0.106 0.108 2.77% Amr::writePlotFile() 2 0.1038 0.104 0.104 2.67% MultiFab::AddProduct() 568 0.09359 0.09463 0.09893 2.54% Castro::normalize_species() 13 0.08472 0.08617 0.09796 2.51% Castro::reset_internal_energy(Fab) 28 0.0885 0.08945 0.09515 2.44% Castro::post_timestep() 1 0.07096 0.08853 0.09218 2.37% Castro::post_regrid() 2 0.08742 0.0889 0.09132 2.34% Castro::derive() 18 0.07225 0.07347 0.07938 2.04% MLMG::mgVcycle_up::1 23 0.01618 0.0457 0.07547 1.94% MLPoisson::Fsmooth() 1104 0.05642 0.06085 0.07298 1.87% FillPatchIterator::Initialize 45 0.05038 0.06203 0.06999 1.80% check_for_negative_density() 2 0.002393 0.04003 0.06805 1.75% MLMG::addInterpCorrection() 138 0.01125 0.05552 0.06746 1.73% FillPatchSingleLevel 45 0.04641 0.0586 0.06679 1.71% MultiFab::Subtract() 27 0.05728 0.05859 0.06669 1.71% MLCellLinOp::applyBC() 1285 0.04873 0.05435 0.0666 1.71% FabArray::setVal() 852 0.05973 0.0611 0.06363 1.63% MLMG::mgVcycle_up::2 23 0.001989 0.02132 0.05094 1.31% Castro::initialize_advance() 1 0.04925 0.04971 0.05026 1.29% FabArray::FillBoundary() 1147 0.037 0.043 0.04942 1.27% Amr::regrid() 1 0.04476 0.04578 0.04834 1.24% MLMG::mgVcycle_down::0 23 0.04564 0.04602 0.04649 1.19% Gravity::fill_multipole_BCs() 4 0.04128 0.04314 0.04396 1.13% FabArray::ParallelCopy_nowait() 277 0.03704 0.03876 0.04233 1.09% Castro::problem_post_timestep() 2 0.008906 0.0348 0.04 1.03% Castro::wd_update() 2 0.008787 0.03476 0.03998 1.03% MLMG::mgVcycle_up::0 23 0.03542 0.03637 0.0379 0.97% MLCellLinOp::apply() 178 0.02758 0.02964 0.03598 0.92% MLMG::mgVcycle_up::3 23 0.001996 0.004436 0.03435 0.88% FabArray::sum() 568 0.03077 0.03121 0.03252 0.83% Castro::construct_new_source() 12 0.01445 0.01556 0.0292 0.75% VisMF::Write(FabArray) 4 0.01239 0.01367 0.02872 0.74% MLCellLinOp::correctionResidual() 139 0.01957 0.02123 0.02596 0.67% MLMG:computeResOfCorrection() 138 0.01983 0.02145 0.02546 0.65% FillBoundary_nowait() 1147 0.01536 0.01998 0.02501 0.64% FillBoundary_finish() 1147 0.01615 0.02042 0.02447 0.63% Castro::sum_integrated_quantities() 2 0.02385 0.0241 0.02413 0.62% Castro::expand_state() 2 0.01863 0.02241 0.02391 0.61% MLPoisson::Fapply() 178 0.01732 0.0182 0.02081 0.53% MLMG::mgVcycle_down::1 23 0.01971 0.02022 0.02074 0.53% FabArray::Saxpy() 27 0.007543 0.008384 0.01972 0.51% MLCellLinOp::solutionResidual() 27 0.01675 0.01814 0.01914 0.49% Gravity::test_composite_phi() 1 0.01386 0.01586 0.01607 0.41% MLMG::computeResidual() 23 0.01446 0.01499 0.01596 0.41% MLMG::mgVcycle_bottom 23 6.699e-05 0.001034 0.01552 0.40% MLMG::actualBottomSolve() 23 2.5e-05 0.000991 0.01547 0.40% MLCGSolver::bicgstab 1 0 0.0009484 0.01517 0.39% MLMG::mgVcycle_down::2 23 0.002337 0.007109 0.01225 0.31% Castro::create_source_corrector() 2 0.005823 0.00977 0.01174 0.30% FabArray::Xpay() 171 0.009017 0.009563 0.01089 0.28% MLMG::mgVcycle_down::3 23 0.002217 0.003301 0.01064 0.27% amrex::average_down 138 0.005991 0.007074 0.009738 0.25% Amr::writeSmallPlotFile() 2 0.009321 0.009323 0.009324 0.24% Castro::volWgtSum() 46 0.007349 0.007532 0.009095 0.23% Castro::construct_old_source() 12 0.007844 0.007933 0.008622 0.22% MLPoisson::define() 4 0.006307 0.007325 0.00853 0.22% Castro::construct_new_gravity_source() 2 0.008027 0.008163 0.00833 0.21% MLMG::ResNormInf() 26 0.002961 0.00597 0.007613 0.20% Castro::construct_old_gravity() 2 0.003988 0.00527 0.007479 0.19% Gravity::get_old_grav_vector() 2 0.003978 0.005259 0.007468 0.19% Castro::apply_source_to_state() 4 0.003007 0.003273 0.007132 0.18% MLMG::mgVcycle_down::4 23 0.002193 0.002483 0.006304 0.16% Gravity::get_new_grav_vector() 3 0.005927 0.006068 0.006243 0.16% Castro::check_for_nan() 4 0.004809 0.00492 0.005499 0.14% MultiFab::contains_nan() 4 0.004798 0.004908 0.005487 0.14% Castro::volInBoundary() 14 0.003451 0.003739 0.0054 0.14% Castro::construct_old_gravity_source() 2 0.004613 0.004678 0.005363 0.14% Gravity::update_max_rhs() 2 0.001442 0.00389 0.00535 0.14% MLCellLinOp::defineAuxData() 4 0.00405 0.004387 0.005184 0.13% MLMG::mgVcycle_up::4 23 0.001998 0.00221 0.004907 0.13% MLMG::mgVcycle_down::5 23 0.002178 0.002337 0.00419 0.11% MultiFab::max() 2 0.0001529 0.002632 0.004098 0.11% MLMG::prepareForSolve() 3 0.003463 0.003589 0.003954 0.10% StateDataPhysBCFunct::() 45 0.000168 0.002073 0.003811 0.10% StateData::FillBoundary(geom) 16 0 0.00189 0.003629 0.09% MLMG::mgVcycle_up::5 23 0.002005 0.002151 0.003579 0.09% amrex::Add() 24 0.002572 0.002636 0.003168 0.08% MLMG::getGradSolution() 3 0.001205 0.001427 0.003139 0.08% MLCellLinOp::compGrad() 3 0.001196 0.001418 0.003127 0.08% MLMG::computeMLResidual() 3 0.00178 0.002507 0.003117 0.08% Castro::finalize_do_advance() 2 0.002408 0.002506 0.003083 0.08% FabArrayBase::getCPC() 373 0.002189 0.002538 0.003058 0.08% FabArrayBase::getFB() 1147 0.002012 0.002191 0.002733 0.07% FabArray::norminf() 43 0.001994 0.002124 0.002525 0.06% MLLinOp::define() 4 0.0008528 0.001567 0.002452 0.06% MLLinOp::defineGrids() 4 0.0008351 0.00155 0.002436 0.06% Castro::buildMetrics() 1 0.0006282 0.0008379 0.002193 0.06% AmrLevel::AmrLevel(dm) 1 0.001423 0.00151 0.001809 0.05% MLLinOp::makeSubCommunicator() 4 0.0001683 0.0009174 0.001789 0.05% StateData::define() 6 0.001345 0.001431 0.001722 0.04% MLCGSolver::ParallelAllReduce 32 0 0.0001066 0.001706 0.04% FabArrayBase::CPC::define() 126 0.0008177 0.001147 0.001629 0.04% MLCellLinOp::defineBC() 4 0.001275 0.001337 0.001413 0.04% FabArray::setDomainBndry() 45 0.00013 0.0007531 0.001329 0.03% FabArray::mult() 11 0.001167 0.001185 0.00122 0.03% MLMG::compResidual() 1 0.0009391 0.001082 0.001212 0.03% Castro::locWgtSum() 6 0.001159 0.001174 0.001189 0.03% amrex::Dot() 24 0 6.977e-05 0.001116 0.03% BndryData::define() 4 0.0009174 0.000994 0.001056 0.03% TagBoxArray::collate() 2 0.0001247 0.000509 0.001047 0.03% DistributionMapping::KnapSackDoIt() 1 0.0002006 0.000907 0.001006 0.03% DistributionMapping::LeastUsedCPUs() 1 2.779e-05 0.0008138 0.0009264 0.02% Castro::enforce_consistent_e() 1 0.000358 0.0003934 0.0007699 0.02% FabArray::LinComb() 12 0 4.284e-05 0.0006855 0.02% Castro::errorEst() 2 0.0005094 0.0005692 0.0006255 0.02% MLPoisson::normalize() 13 0 3.583e-05 0.0005732 0.01% Castro::apply_tagging_restrictions() 2 0.0004058 0.0004655 0.0005216 0.01% Gravity::applyMetricTerms() 2 0.0004716 0.0004752 0.0004849 0.01% MLPoisson::prepareForSolve() 4 0.0003798 0.0004082 0.0004737 0.01% MLCellLinOp::prepareForSolve() 4 0.0003545 0.0003848 0.0004503 0.01% Castro::swap_state_time_levels() 1 0.0003753 0.0003886 0.0004248 0.01% Castro::volProductSum() 2 0.0003963 0.0004079 0.0004189 0.01% FabArrayBase::FB::FB() 14 0.0001612 0.0002682 0.0003913 0.01% FabArray::setVal(val, thecmd, scomp, ncomp) 128 0.000274 0.0003242 0.0003727 0.01% Castro::finalize_advance() 1 0.000279 0.0002921 0.0003387 0.01% MLLinOp::makeAgglomeratedDMap 4 0.0002026 0.0002145 0.0002257 0.01% MLMG::MLResNormInf() 3 0.0002043 0.0002074 0.0002127 0.01% MLMG::MLRhsNormInf() 3 0.0001986 0.0002024 0.0002092 0.01% Amr::initSubcycle() 1 6.434e-05 8.194e-05 0.0002005 0.01% MLCellLinOp::setLevelBC() 4 0.0001175 0.0001448 0.000182 0.00% OwnerMask() 2 0.000136 0.0001551 0.0001713 0.00% Gravity::swapTimeLevels() 1 0.0001377 0.0001414 0.000147 0.00% makeSFC 20 0.0001297 0.0001385 0.0001458 0.00% knapsack() 1 5.801e-05 7.292e-05 0.0001369 0.00% Castro::apply_problem_tags() 2 9.465e-05 9.709e-05 9.853e-05 0.00% DistributionMapping::Distribute() 20 4.553e-05 4.844e-05 5.685e-05 0.00% Castro::computeInitialDt() 2 1.014e-05 1.205e-05 2.967e-05 0.00% Castro::enforce_speed_limit() 13 1.802e-05 1.923e-05 2.221e-05 0.00% BoxList::parallelComplementIn() 4 1.32e-05 1.403e-05 1.851e-05 0.00% Gravity::set_mass_offset() 2 4.602e-06 5.599e-06 1.34e-05 0.00% Castro::estTimeStep() 3 3.934e-06 4.733e-06 9.46e-06 0.00% Castro::gwstrain() 2 4.33e-06 5.535e-06 8.665e-06 0.00% BoxList::complementIn 4 6.127e-06 6.451e-06 8.347e-06 0.00% Castro::retry_advance_ctu() 1 3.699e-06 4.173e-06 4.944e-06 0.00% Castro::construct_old_rotation_source() 2 3.751e-06 3.967e-06 4.292e-06 0.00% Castro::construct_new_rotation_source() 2 3.222e-06 3.442e-06 4.04e-06 0.00% Castro::FluxRegCrseInit 1 1.768e-06 2.23e-06 3.183e-06 0.00% Castro::FluxRegFineAdd() 1 1.36e-06 1.944e-06 2.686e-06 0.00% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::castro.dtnuc_mode(nvals = 1) :: [1] [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 47 B 54 B 71 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (23.06-14-g7c7a8fb1fff5) finalized