Cheers Matthew...that worked now. Matthew Knepley wrote: > I think you need to tell use that the IS is a permutation: > > http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/IS/ISSetPermutation.html > > Matt > > On 10/25/07, Tim Stitt <timothy.stitt at ichec.ie> wrote: > >> 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) >> >> >> > > >
-- 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)
