[petsc-users] Issue updating MUMPS ictnl after failed solve

David Knezevic david.knezevic at akselos.com
Wed Sep 21 16:48:48 CDT 2016


On Tue, Sep 20, 2016 at 2:24 PM, Hong <hzhang at mcs.anl.gov> wrote:

> David:
> I'll patch petsc-maint (v3.7), then merge it to petsc-master.
> It might take 1-2 days for regression tests.
>
>
I pulled the patched version and it works well for me now, thanks!

David




> 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/20160921/88118564/attachment.html>


More information about the petsc-users mailing list