[petsc-users] strange error

[Juha Jäykkä]

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.

Any help is appreciated...

-Juha