<div dir="ltr"><div><div><div><div><div><div>Greetings,<br><br></div>We have a test that has started failing upon switching from 3.5.4 to 3.6.0 (actually went straight to 3.6.3 but checked this is repeatable with 3.6.0). I've attached the matrix generated with -mat_view binary and a small PETSc program that runs in serial that reproduces the behavior by loading the matrix and solving a linear system (RHS doesn't matter here). For context, this matrix is the Jacobian of a Taylor-Hood approximation of the axisymmetric incompressible Navier-Stokes equations for flow between concentric cylinders (for which there is an exact solution). The matrix is for a two element case, hopefully small enough for debugging.<br><br></div>Using the following command line options with the test program works with PETSc 3.5.4 and gives a NAN residual with PETSc 3.6.0:<br><br>PETSC_OPTIONS="-pc_type asm -pc_asm_overlap 12 -sub_pc_type ilu -sub_pc_factor_mat_ordering_type 1wd -sub_pc_factor_levels 4"<br><br></div>If I remove the mat ordering option, all is well again in PETSc 3.6.x:<br><br>PETSC_OPTIONS="-pc_type asm -pc_asm_overlap 12 -sub_pc_type ilu -sub_pc_factor_levels 4"<br><br></div>Those options are nothing special. They were arrived at through trial/error to get decent behavior for the solver on up to 4 processors to keep the time to something reasonable for the test suite without getting really fancy. Specifically, we'd noticed this mat ordering on some problems in the test suite behaved noticeably better.</div><div><br></div><div>As always, thanks for your time.<br></div><div><br></div>Best,<br><br></div>Paul<br><br><br></div>