Invalid argument! - not a permutation!

[Tim Stitt]

Hi all,

Following on from some help I got a few days ago with respect to 
MatLUFactorSymbolic() usage in Fortran...I think I now have all the 
parts in place but keep getting the following runtime error upon execution:

[0]PETSC ERROR: --------------------- Error Message 
------------------------------------
[0]PETSC ERROR: Invalid argument!
[0]PETSC ERROR: not a permutation!
...
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: ISInvertPermutation() line 183 in 
src/vec/is/interface/index.c
[0]PETSC ERROR: MatLUFactorSymbolic_SeqAIJ() line 276 in 
src/mat/impls/aij/seq/aijfact.c
[0]PETSC ERROR: MatLUFactorSymbolic() line 2178 in 
src/mat/interface/matrix.c
[localhost.localdomain:24730] MPI_ABORT invoked on rank 0 in 
communicator MPI_COMM_WORLD with errorcode 62

Here is the portion of  the offending code (I think):

  call MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY,error);CHKERRQ(error)
  call MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY,error);CHKERRQ(error)
  call ISCreateGeneral(PETSC_COMM_WORLD,order,columnIndices, 
indexSet,error);CHKERRQ(error)
  call MatFactorInfoInitialize(info,error);CHKERRQ(error)
  call 
MatLUFactorSymbolic(A,indexSet,indexSet,info,factorMat,error);CHKERRQ(error)

order - number of rows/columns in square matrix A
columnIndices - 1D array containing values 0-->order-1

All types are set correctly as far as I can tell.

Can anyone suggest why the execution falls over during the 
MatLUFactorSymbolic() call?

Thanks again,

Tim.

-- 
Dr. Timothy Stitt <timothy_dot_stitt_at_ichec.ie>
HPC Application Consultant - ICHEC (www.ichec.ie)

Dublin Institute for Advanced Studies
5 Merrion Square - Dublin 2 - Ireland

+353-1-6621333 (tel) / +353-1-6621477 (fax)