<div dir="ltr">Hi Folks,<div><br></div><div>We are dealing with a singular matrix that arises from Stokes' equation with Dirichlet velocity BCs. The way we are handling it is to pass nullspace info to KSP, i.e, vector of {u,p} = {0,1}, and letting the matrix to be singular. However, we are using ASM as the preconditioner and when the subpc covers the entire domain (which can be the case in the preconditioning step of our MG algorithm), LU decomposition fails. We are thinking of pinning one of the pressure DOFs to zero explicitly, so that we can avoid singular matrix for the LU decomposition. Is there another method of doing it, which avoids explicit pinning down of an arbitrary pressure DOF?</div><div><br></div><div>Thanks,</div><div>--Amneet</div><div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div>--Amneet <br><br></div><div><br></div><div><br></div></div></div>
</div></div>