[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Invalid argument! - not a permutation!
- To: petsc-users@xxxxxxxxxxx
- Subject: Re: Invalid argument! - not a permutation!
- From: "Matthew Knepley" <knepley@xxxxxxxxx>
- Date: Thu, 25 Oct 2007 14:00:09 -0500
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=beta; h=domainkey-signature:received:received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; bh=8jBz27dyN2BLxWZvszLKKIyE6/AGvDFamEcHvLacRSo=; b=UREdUnHtfAAPabuAkRhDhhKrFUTM/v3Drmv3aENHdNDWz2EkfZdiZb65FqcRR2ig3HExpkTbKjpg/ahhWqMqf9QTCix2X9bVpQHecGfO0FgSbrJ4s/cvBWyxQte558ctDq17kWqlNMQ97o2+tyg3XEKM71ObNU/D5vv7BC8P0ps=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=beta; h=received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=j9uMfNeCLv7px1VSO2DEqdLTUx/kED3JfQL1V/hjzkiSUIKJ9MijK2Fcn9DlNHOPSAzaz6VyPAysZ4SxWNQu4xuA6m4UedBGoOfFOer3AY/fEkb2F88GG3Sg++m9IHR5PP3ZRC8xDXV826GSNtL3Q8JRj5z4YOU/UbdJoUvWm1Y=
- In-reply-to: <4720E65C.9070504@ichec.ie>
- References: <4720E65C.9070504@ichec.ie>
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
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@xxxxxxxx> 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)
>
>
--
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