[petsc-users] Issue updating MUMPS ictnl after failed solve
David Knezevic
david.knezevic at akselos.com
Tue Sep 20 13:08:38 CDT 2016
On Tue, Sep 20, 2016 at 2:01 PM, Hong <hzhang at mcs.anl.gov> wrote:
> David :
> This is a bug in PETSc. A change
> $ git diff ../../../pc/impls/factor/cholesky/cholesky.c
> diff --git a/src/ksp/pc/impls/factor/cholesky/cholesky.c
> b/src/ksp/pc/impls/factor/cholesky/cholesky.c
> index 953d551..cc28369 100644
> --- a/src/ksp/pc/impls/factor/cholesky/cholesky.c
> +++ b/src/ksp/pc/impls/factor/cholesky/cholesky.c
> @@ -141,9 +141,7 @@ static PetscErrorCode PCSetUp_Cholesky(PC pc)
>
> ierr = MatCholeskyFactorNumeric(((PC_Factor*)dir)->fact,pc->pmat,&(
> (PC_Factor*)dir)->info);CHKERRQ(ierr);
> ierr = MatFactorGetError(((PC_Factor*)dir)->fact,&err);CHKERRQ(ierr);
> - if (err) { /* FactorNumeric() fails */
> - pc->failedreason = (PCFailedReason)err;
> - }
> + pc->failedreason = (PCFailedReason)err;
> }
>
> fixed the problem. I'll fix this problem in petsc-release, including other
> routines.
> Thanks for reporting the bug and sending matrix.dat. Let us know whenever
> you encounter problem using PETSc.
>
OK, great, thanks for the fix.
Will this fix be included in the next patch release of 3.7?
Thanks,
David
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160920/9ee74f27/attachment.html>
More information about the petsc-users
mailing list