On 10/25/07, Tim Stitt <timothy.stitt at ichec.ie> wrote: > OK...that is useful thanks. > > If you don't mind just one more question...I notice that the > MatLUFactorSymbolic() routine now executes with one processor but > complains when I run > 2 processors. The actual error message is: > > [1]PETSC ERROR: No support for this operation for this object type! > [1]PETSC ERROR: Matrix type mpiaij symbolic LU! > > Is this because the default LU factorisation routine isn't parallel...I > just assumed parallel superlu would be picked up automatically (it is > compiled into my PETSc build) ?
It is only valid for the superlu matrix type, -mat_type superlu_dist. Matt > Tim. > > Barry Smith wrote: > > I have added this information to the error message that gets printed. > > > > Barry > > > > > > On Thu, 25 Oct 2007, 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) > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener
