[petsc-users] MatMatSolve in sequential call

[Alexander Grayver]

Hello,

Trying to call this routine on 1 proc:

call MatMatSolve(F,B,solutions,ierr);CHKERRQ(ierr)

I get following crash:

[0]PETSC ERROR: --------------------- Error Message 
------------------------------------
[0]PETSC ERROR: Null argument, when expecting valid pointer!
[0]PETSC ERROR: Null Object: Parameter # 1!
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 2, Fri Jul 13 
15:42:00 CDT 2012
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: /app on a openmpi-i named glic1 by ag Mon Aug 13 
15:06:36 2012
[0]PETSC ERROR: Libraries linked from 
/home/lib/petsc-3.3-p2/openmpi-intel-complex-debug-f-ds/lib
[0]PETSC ERROR: Configure run at Wed Jul 18 12:14:30 2012
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: ISGetIndices() line 370 in 
/home/lib/petsc-3.3-p2/src/vec/is/interface/index.c
[0]PETSC ERROR: MatMatSolve_SeqAIJ() line 1109 in 
/home/lib/petsc-3.3-p2/src/mat/impls/aij/seq/aijfact.c
[0]PETSC ERROR: MatMatSolve() line 3122 in 
/home/lib/petsc-3.3-p2/src/mat/interface/matrix.c

The same code works fine on nproc > 1.

Any ideas?

-- 
Regards,
Alexander