[petsc-users] Error with SBAIJ during KSPSolve

[Tabrez Ali]

Hi

I am trying to solve a system with constraints (0 on some diagonals). It 
works fine with AIJ but gives the following error (see below) with SBAIJ 
Matrices during KSPSolve (sequential). With GMRES it just segfaults.

What did I miss?

Thanks in advance.
Tabrez

stali at x61:~/src$ ./a.out <mesh_dense.inp -pc_type lu 
-pc_factor_shift_type NONZERO -pc_factor_shift_amount 1

  Reading input file ...
  Calculating NZ structure ...
  Forming [K] ...
  Forming RHS ...
  Setting up solver ...
  Solving ...
[0]PETSC ERROR: --------------------- Error Message 
------------------------------------
[0]PETSC ERROR: Invalid argument!
[0]PETSC ERROR: Index set is not a permutation!
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 5, Mon Sep 27 
11:51:54 CDT 2010
[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: ./a.out on a linux-gnu named x61 by stali Tue Apr 26 
00:55:04 2011
[0]PETSC ERROR: Libraries linked from /opt/petsc-3.1-p5/lib
[0]PETSC ERROR: Configure run at Mon Feb  7 00:03:37 2011
[0]PETSC ERROR: Configure options --prefix=/opt/petsc-3.1-p5 
--with-mpi-dir=/opt/mpich2-1.0.8 --with-fortran-interfaces=1 
--with-parmetis=1 --download-parmetis=ifneeded --with-shared=1
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: ISSetPermutation() line 141 in src/vec/is/interface/index.c
[0]PETSC ERROR: MatGetOrdering() line 231 in src/mat/order/sorder.c
[0]PETSC ERROR: PCSetUp_LU() line 125 in src/ksp/pc/impls/factor/lu/lu.c
[0]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
[0]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c
  Done
  Cleaning up ...
  Finished