[petsc-dev] getting ready for PETSc release
Barry Smith
bsmith at mcs.anl.gov
Sun Jun 8 17:12:02 CDT 2014
Peter and Lisandro,
Sorry I don’t remember this.
In looking at the code this would happen if snes->ksp doesn’t yet exist. We could
1) have SNESComputeJacobian() call SNESGetKSP() to force the existence of KSP. or
2) only call the KSPSetReusePreconditioner() in the various places if snes->ksp exists.
The issue is if the Jacobian is ever computed but KSP is not used? For example, for some weird non-Newton nonlinear solver. If the Jacobian computation always goes hand in hand with having/using a KSP then I would do 1) otherwise 2.
Barry
On Jun 8, 2014, at 3:32 PM, Lisandro Dalcin <dalcinl at gmail.com> wrote:
> On 5 June 2014 16:37, Barry Smith <bsmith at mcs.anl.gov> wrote:
>>
>> We’ll be making a PETSc release in a few days so please get anything you need into next for testing and movement over to master.
>>
>> Barry
>>
>> Tentative plan for release Wed June 11
>>
>
> Barry, I'm still getting this error in petsc4py, we discussed the
> issue some time ago.
>
> Traceback (most recent call last):
> File "test/test_snes.py", line 155, in testCompJac
> self.snes.computeJacobian(x, J)
> File "SNES.pyx", line 236, in petsc4py.PETSc.SNES.computeJacobian
> (src/petsc4py.PETSc.c:141969)
> Error: error code 85
> [0] SNESComputeJacobian() line 2209 in
> /Users/dalcinl/Devel/petsc-dev/src/snes/interface/snes.c
> [0] KSPSetReusePreconditioner() line 192 in
> /Users/dalcinl/Devel/petsc-dev/src/ksp/ksp/interface/itfunc.c
> [0] Null argument, when expecting valid pointer
> [0] Null Object: Parameter # 1
>
>
>
>
> --
> Lisandro Dalcin
> ---------------
> CIMEC (UNL/CONICET)
> Predio CONICET-Santa Fe
> Colectora RN 168 Km 472, Paraje El Pozo
> 3000 Santa Fe, Argentina
> Tel: +54-342-4511594 (ext 1016)
> Tel/Fax: +54-342-4511169
More information about the petsc-dev
mailing list