<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">David:</div><div class="gmail_quote">I'll patch petsc-maint (v3.7), then merge it to petsc-master.</div><div class="gmail_quote">It might take 1-2 days for regression tests.</div><div class="gmail_quote"><br></div><div class="gmail_quote">Hong</div><div class="gmail_quote"><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class="">On Tue, Sep 20, 2016 at 2:01 PM, Hong <span dir="ltr"><<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@mcs.anl.gov</a>></span> wrote:<br></span><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">David :</div><span class=""><div class="gmail_quote">This is a bug in PETSc. A change </div><div class="gmail_quote"><div class="gmail_quote">$ git diff ../../../pc/impls/factor/chole<wbr>sky/cholesky.c</div><div class="gmail_quote">diff --git a/src/ksp/pc/impls/factor/chol<wbr>esky/cholesky.c b/src/ksp/pc/impls/factor/chol<wbr>esky/cholesky.c</div><div class="gmail_quote">index 953d551..cc28369 100644</div><div class="gmail_quote">--- a/src/ksp/pc/impls/factor/chol<wbr>esky/cholesky.c</div><div class="gmail_quote">+++ b/src/ksp/pc/impls/factor/chol<wbr>esky/cholesky.c</div><div class="gmail_quote">@@ -141,9 +141,7 @@ static PetscErrorCode PCSetUp_Cholesky(PC pc)</div><div class="gmail_quote"><br></div><div class="gmail_quote">     ierr = MatCholeskyFactorNumeric(((PC_<wbr>Factor*)dir)->fact,pc->pmat,&(<wbr>(PC_Factor*)dir)->info);CHKERR<wbr>Q(ierr);</div><div class="gmail_quote">     ierr = MatFactorGetError(((PC_Factor*<wbr>)dir)->fact,&err);CHKERRQ(ierr<wbr>);</div><div class="gmail_quote">-    if (err) { /* FactorNumeric() fails */</div><div class="gmail_quote">-      pc->failedreason = (PCFailedReason)err;</div><div class="gmail_quote">-    }</div><div class="gmail_quote">+    pc->failedreason = (PCFailedReason)err;</div><div class="gmail_quote">   }</div><div class="gmail_quote"><br></div><div class="gmail_quote">fixed the problem. I'll fix this problem in petsc-release, including other routines. </div><div class="gmail_quote">Thanks for reporting the bug and sending matrix.dat. Let us know whenever you encounter problem using PETSc.</div></div></span></div></div></blockquote><div><br></div><div><br></div><div>OK, great, thanks for the fix.</div><div><br></div><div>Will this fix be included in the next patch release of 3.7?</div><div><br></div><div>Thanks,</div><div>David</div><div><br></div><div><br></div><div> </div></div></div></div>
</blockquote></div><br></div></div>