[petsc-users] strange error

[Matthew Knepley]

On Fri, Apr 15, 2011 at 9:17 AM, Juha Jäykkä <J.Jaykka at leeds.ac.uk> wrote:

> Hi list!
>
> I keep getting strange errors when running a PETSc code:
>
> [6]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> [6]PETSC ERROR: Object is in wrong state!
> [6]PETSC ERROR: Matrix must be set first!
> [6]PETSC ERROR:
> ------------------------------------------------------------------------
>
> This happens on one machine only, but it has OpenMPI just as the others,
> where
> it works correctly.
>
> More specifically, the error comes from
>
> [6]PETSC ERROR: PCSetUp() line 775 in src/ksp/pc/interface/precon.c
> [6]PETSC ERROR: PCApply() line 353 in src/ksp/pc/interface/precon.c
> [6]PETSC ERROR: KSPSolve_PREONLY() line 29 in
> src/ksp/ksp/impls/preonly/preonly.c
> [6]PETSC ERROR: KSPSolve() line 396 in src/ksp/ksp/interface/itfunc.c
> [6]PETSC ERROR: PCApply_BJacobi_Singleblock() line 777 in
> src/ksp/pc/impls/bjacobi/bjacobi.c
> [6]PETSC ERROR: PCApply() line 357 in src/ksp/pc/interface/precon.c
> [6]PETSC ERROR: KSPInitialResidual() line 65 in
> src/ksp/ksp/interface/itres.c
> [6]PETSC ERROR: KSPSolve_GMRES() line 240 in
> src/ksp/ksp/impls/gmres/gmres.c
> [6]PETSC ERROR: KSPSolve() line 396 in src/ksp/ksp/interface/itfunc.c
> [6]PETSC ERROR: SNES_KSPSolve() line 2944 in src/snes/interface/snes.c
> [6]PETSC ERROR: SNESSolve_LS() line 191 in src/snes/impls/ls/ls.c
> [6]PETSC ERROR: SNESSolve() line 2255 in src/snes/interface/snes.c
> [6]PETSC ERROR: TSStep_BEuler_Nonlinear() line 176 in
> src/ts/impls/implicit/beuler/beuler.c
> [6]PETSC ERROR: TSStep() line 1693 in src/ts/interface/ts.c
> [6]PETSC ERROR: TSSolve() line 1731 in src/ts/interface/ts.c
>
> Now, I do realize that the Jacobian and preconditioner matrices must be
> properly created before calling TSSolve, but, to the best of my knowledge,
> they are. If they were not, the code should never work.
>
> And this always comes from just one MPI rank, as if somehow its memory is
> corrupt or something. Funny thing is, it only happens under the batch queue
> system: interactively, on the same machine, it works fine.
>

I would try valgrind on it first.

   Matt


> Any help is appreciated...
>
> -Juha
>



-- 
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110415/c4cf5004/attachment.htm>