Petsc + BlockSolve95

Barry Smith bsmith at mcs.anl.gov
Wed May 10 08:29:30 CDT 2006


      No, MatSetOption() should always work and accept these arguments.
Are you sure the values of MAT_SYMMETRIC and MAT_SYMMETRY_ETERNAL
are being set (by including the right Fortran include file) before you
use them?

       Barry
BTW: is this PETSc 2.3.1?



On Wed, 10 May 2006, Boulahya Faïza wrote:

> Thanks.
>
> I tried something else :
>    - creation of the matrix with MatCreateMPIAIJ
>    - initialization
>    - adding of the options MAT_SYMMETRIC and MAT_SYMMETRY_ETERNAL
>
> When solving in sequential CG + ICC, everything is ok. When I tried in
> parallel the same code the options lead to the same error :
> [0]PETSC ERROR: MatSetOption_MPIAIJ() line 1251 in
> src/mat/impls/aij/mpi/mpiaij.c
> [0]PETSC ERROR: No support for this operation for this object type!
> [0]PETSC ERROR: unknown option!
> [0]PETSC ERROR: MatSetOption() line 4137 in src/mat/interface/matrix.c
>
> Can MatSetOption be only used in sequential?
>
>
> Faïza
>
>
>  _____
>
> De : Matthew Knepley [mailto:knepley at gmail.com]
> Envoyé : mardi 9 mai 2006 18:19
> À : petsc-users at mcs.anl.gov
> Objet : Re: Petsc + BlockSolve95
>
>
> I believe there is a problem with the option that you specified. All these
> are integers, and it is complaining that the integer does not match
> MAT_SYMMETRIC. I will fix the error message to print the offending
> option, but please check the code.
>
>  Thanks,
>
>     Matt
>
>
> On 5/9/06, Boulahya Faïza <F.Boulahya at brgm.fr <mailto:F.Boulahya at brgm.fr> >
> wrote:
>
> Hi All,
>
> Has anyone used Conjugate Gradient Solver + Icomplete Cholesky
> Preconditionner in parallel case? I tried as said in the manual : I use
> MATMPIROWBS for the storage of the matrice. However I get this message :
>
> PETSC ERROR: To use incomplete Cholesky
>                        preconditioning with a MATMPIROWBS matrix you must
> declare it to be
>                        symmetric using the option
> MatSetOption(A,MAT_SYMMETRIC)!
>
>
> So I tried adding this option (even if in the namual it is written that it
> is not required). Then I obtained this message
> PETSC ERROR: MatSetOption_MPIRowbs() line 1411 in
> src/mat/impls/rowbs/mpi/mpirowbs.c
> PETSC ERROR: No support for this operation for this object type!
> PETSC ERROR: unknown option!
> PETSC ERROR: MatSetOption() line 4137 in src/mat/interface/matrix.c
>
>
> In advance thanks,
>
>
> Faïza
>
>
> Pensez à visiter le site BRGM sur.... http://www.brgm.fr/
>
> ***
> Le contenu de cet e-mail et de ses pièces jointes est destiné à l'usage exclusif du
> (des) destinataire(s) expressément désigné(s) comme tel(s). En cas de réception de cet
> e-mail par erreur, le signaler à son expéditeur et ne pas en divulguer le contenu.
> L'absence de virus a été vérifié à  l'émission du message. Il convient néanmoins de
> vérifier l'absence de corruption à sa réception.
>
> The contents of this email and any attachments are confidential. They are intended for
> the named recipient(s) only. If you have received this email in error please notify the
> system manager or  the sender immediately and do not disclose the contents to
> anyone or make copies. eSafe scanned this email for viruses, vandals and malicious
> content.
> ***


More information about the petsc-users mailing list