<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, May 27, 2015 at 10:52 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
  The built in PETSc LU factorizations have special options to attempt to "shift" a singular matrix (by adding something to a zero pivot). These methods are ad hoc but sometimes perform satisfactory inside a preconditioner. See the manual page for PCLU in particular -pc_factor_shift_type <shifttype> - Sets shift type or PETSC_DECIDE for the default; use '-help' for a list of available types  See also <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#zeropivot" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html#zeropivot</a></blockquote><div><br></div><div>The other option is to use -pc_type svd for the coarse grid.</div><div><br></div><div>  Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="HOEnZb"><font color="#888888"><br>
  Barry<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
> On May 27, 2015, at 10:40 PM, Amneet Bhalla <<a href="mailto:mail2amneet@gmail.com">mail2amneet@gmail.com</a>> wrote:<br>
><br>
> Hi Folks,<br>
><br>
> 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?<br>
><br>
> Thanks,<br>
> --Amneet<br>
><br>
> --<br>
> --Amneet<br>
><br>
><br>
><br>
<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>