<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META content="MSHTML 6.00.2900.2180" name=GENERATOR></HEAD>
<BODY>
<DIV dir=ltr align=left><SPAN class=940290913-17052006><FONT face=Arial
size=2>Doing this, I realize that I wrote <FONT face="Times New Roman">MAT
SYMMETRIC</FONT> instead of <FONT
face="Times New Roman">MAT_SYMMETRIC</FONT>. Everything is
ok.</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=940290913-17052006><FONT face=Arial
size=2>Sorry for all those mails...</FONT></SPAN></DIV>
<DIV> </DIV>
<DIV align=left><FONT face=Arial size=2>Faïza</FONT></DIV>
<DIV> </DIV><BR>
<DIV class=OutlookMessageHeader lang=fr dir=ltr align=left>
<HR tabIndex=-1>
<FONT face=Tahoma size=2><B>De :</B> Matthew Knepley
[mailto:knepley@gmail.com] <BR><B>Envoyé :</B> mercredi 17 mai 2006
13:18<BR><B>À :</B> petsc-users@mcs.anl.gov<BR><B>Cc :</B> Boulahya
Faiza<BR><B>Objet :</B> Re: Petsc + BlockSolve95<BR></FONT><BR></DIV>
<DIV></DIV>My guess is somehow you are passing the wrong integer for the
option<BR>because I am looking at the code and it is fine. However, it is easy
to<BR>clear up. Go to<BR><BR>
$PETSC_DIR/src/mat/impls/rowbs/mpi/mpirowbs.c:1411<BR><BR>and replace it
with<BR><BR> SETERRQ1(PETSC_ERR_SUP,"unknown option %d",op);<BR><BR>and
then in that directory<BR><BR> make<BR> make shared<BR><BR>This will
tell us what integer was passed in for the option. I have already made<BR>this
fix in the dev version.<BR><BR> Matt<BR><BR>
<DIV><SPAN class=gmail_quote>On 5/17/06, <B class=gmail_sendername>Boulahya
Faïza</B> <<A href="mailto:F.Boulahya@brgm.fr">F.Boulahya@brgm.fr</A>>
wrote:</SPAN>
<BLOCKQUOTE class=gmail_quote
style="PADDING-LEFT: 1ex; MARGIN: 0pt 0pt 0pt 0.8ex; BORDER-LEFT: rgb(204,204,204) 1px solid">
<DIV>
<DIV>
<P><FONT size=2>Ok but the error does not come from parallell ICC. It's from
the MatSetOption().</FONT> </P>
<P><FONT size=2>Moreover, when I tried with MATMPIROWBS , it's the same
problem :</FONT> <BR><FONT size=2>[0]PETSC ERROR: MatSetOption_MPIRowbs() line
1411 in src/mat/impls/rowbs/mpi/mpirowbs.c</FONT> </P></DIV>
<DIV><SPAN class=q><BR><FONT size=2>[0]PETSC ERROR: No support for this
operation for this object type!</FONT> <BR><FONT size=2>[0]PETSC ERROR:
unknown option!</FONT> <BR><FONT size=2>[0]PETSC ERROR: MatSetOption() line
4137 in src/mat/interface/matrix.c</FONT> </SPAN></DIV>
<DIV>
<P></P><BR>
<P><FONT size=2>Faïza</FONT> </P>
<P><FONT size=2>-----Message d'origine-----</FONT> <BR><FONT size=2>De :
Satish Balay [<A onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:balay@mcs.anl.gov" target=_blank>mailto:balay@mcs.anl.gov</A>]
</FONT><BR><FONT size=2>Envoyé : mardi 16 mai 2006 15:39</FONT> </P></DIV>
<DIV><SPAN class=q><BR><FONT size=2>À : '<A
onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:petsc-users@mcs.anl.gov"
target=_blank>petsc-users@mcs.anl.gov</A>'</FONT> <BR></SPAN></DIV>
<DIV><FONT size=2>Objet : RE: Petsc + BlockSolve95</FONT>
<P></P></DIV>
<DIV><SPAN class=e id=q_10b41ff7d5e86adf_5>
<P><FONT size=2>> [0]PETSC ERROR: MatSetOption_MPIAIJ() line 1251 in
</FONT><BR><FONT size=2>> src/mat/impls/aij/mpi/mpiaij.c [0]PETSC ERROR: No
support for this operation for this object type!</FONT> </P>
<P><FONT size=2>AIJ matrix does not support parallell ICC. BlockSolve provides
parallel ICC - and to use it one must use MatType MATMPIROWBS</FONT></P>
<P><FONT size=2>Satish</FONT> </P>
<P><FONT size=2>On Tue, 16 May 2006, Boulahya Faïza wrote:</FONT> </P>
<P><FONT size=2>> Is there something new about it?</FONT> </P>
<P><FONT size=2>> De : Vaz, Guilherme [<A
onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:G.Vaz@marin.nl" target=_blank>mailto:G.Vaz@marin.nl</A>]</FONT>
</P>
<P><FONT size=2>> I have exactly the same problem as Faiza...
</FONT><BR><FONT size=2>> In sequential it runs ok but in parallel
not.</FONT> </P>
<P><FONT size=2>> -----Original Message-----</FONT> <BR><FONT size=2>>
From: <A onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:owner-petsc-users@mcs.anl.gov"
target=_blank>owner-petsc-users@mcs.anl.gov</A> </FONT><BR><FONT size=2>>
[<A onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:owner-petsc-users@mcs.anl.gov"
target=_blank>mailto:owner-petsc-users@mcs.anl.gov</A>]</FONT> <BR><FONT
size=2>> On Behalf Of Boulahya Faïza</FONT> </P>
<P><FONT size=2>> I tried something else :</FONT> <BR><FONT
size=2>> </FONT><BR><FONT size=2>> - creation of
the matrix with MatCreateMPIAIJ</FONT> <BR><FONT size=2>> </FONT><BR><FONT
size=2>> - initialization</FONT> <BR><FONT
size=2>> </FONT><BR><FONT size=2>> - adding of
the options MAT_SYMMETRIC and MAT_SYMMETRY_ETERNAL</FONT> <BR><FONT
size=2>> </FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>>
</FONT><BR><FONT size=2>> When solving in sequential CG + ICC, everything
is ok. When I tried in </FONT><BR><FONT size=2>> parallel the same code the
options lead to the same error :</FONT> <BR><FONT size=2>> </FONT><BR><FONT
size=2>> [0]PETSC ERROR: MatSetOption_MPIAIJ() line 1251 in
</FONT><BR><FONT size=2>> src/mat/impls/aij/mpi/mpiaij.c [0]PETSC ERROR: No
support for this </FONT><BR><FONT size=2>> operation for this object
type!</FONT> <BR><FONT size=2>> [0]PETSC ERROR: unknown option!</FONT>
<BR><FONT size=2>> [0]PETSC ERROR: MatSetOption() line 4137 in
src/mat/interface/matrix.c</FONT> <BR><FONT size=2>> </FONT><BR><FONT
size=2>> Can MatSetOption be only used in sequential?</FONT> <BR><FONT
size=2>> </FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>>
</FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>>
</FONT><BR><FONT size=2>> Faïza</FONT> <BR><FONT size=2>>
</FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>>
</FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>>
</FONT><BR><FONT size=2>> _____</FONT> <BR><FONT size=2>>
</FONT><BR><FONT size=2>> De : Matthew Knepley [<A
onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:knepley@gmail.com" target=_blank>mailto:knepley@gmail.com</A>]
Envoyé : mardi 9 mai </FONT><BR><FONT size=2>> 2006 18:19 À : <A
onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:petsc-users@mcs.anl.gov"
target=_blank>petsc-users@mcs.anl.gov</A> Objet : Re: Petsc + </FONT><BR><FONT
size=2>> BlockSolve95</FONT> <BR><FONT size=2>> </FONT><BR><FONT
size=2>> I believe there is a problem with the option that you specified.
All </FONT><BR><FONT size=2>> these are integers, and it is complaining
that the integer does not </FONT><BR><FONT size=2>> match MAT_SYMMETRIC. I
will fix the error message to print the </FONT><BR><FONT size=2>> offending
option, but please check the code.</FONT> <BR><FONT size=2>>
</FONT><BR><FONT size=2>> Thanks,</FONT> <BR><FONT size=2>>
</FONT><BR><FONT size=2>> Matt</FONT>
<BR><FONT size=2>> </FONT><BR><FONT size=2>> On 5/9/06, Boulahya Faïza
<<A onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:F.Boulahya@brgm.fr" target=_blank>F.Boulahya@brgm.fr</A>
</FONT><BR><FONT size=2>> <<A
onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:F.Boulahya@brgm.fr"
target=_blank>mailto:F.Boulahya@brgm.fr</A>> ></FONT> <BR><FONT
size=2>> wrote: </FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>>
Hi All,</FONT> <BR><FONT size=2>> </FONT><BR><FONT size=2>>
</FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>> Has anyone used
Conjugate Gradient Solver + Icomplete Cholesky </FONT><BR><FONT size=2>>
Preconditionner in parallel case? I tried as said in the manual : I
</FONT><BR><FONT size=2>> use MATMPIROWBS for the storage of the matrice.
However I get this message :</FONT> <BR><FONT size=2>> </FONT><BR><FONT
size=2>> </FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>>
PETSC ERROR: To use incomplete Cholesky</FONT> <BR><FONT
size=2>>
preconditioning with a MATMPIROWBS matrix you </FONT><BR><FONT size=2>>
must declare it to be</FONT> <BR><FONT
size=2>>
symmetric using the option </FONT><BR><FONT size=2>>
MatSetOption(A,MAT_SYMMETRIC)!</FONT> <BR><FONT size=2>> </FONT><BR><FONT
size=2>> </FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>> So
I tried adding this option (even if in the namual it is written
</FONT><BR><FONT size=2>> that it is not required). Then I obtained this
message</FONT> <BR><FONT size=2>> </FONT><BR><FONT size=2>> PETSC ERROR:
MatSetOption_MPIRowbs() line 1411 in </FONT><BR><FONT size=2>>
src/mat/impls/rowbs/mpi/mpirowbs.c</FONT> <BR><FONT size=2>> PETSC ERROR:
No support for this operation for this object type!</FONT> <BR><FONT
size=2>> PETSC ERROR: unknown option!</FONT> <BR><FONT size=2>> PETSC
ERROR: MatSetOption() line 4137 in src/mat/interface/matrix.c</FONT> <BR><FONT
size=2>> </FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>>
</FONT><BR><FONT size=2>> In advance thanks,</FONT> <BR><FONT size=2>>
</FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>>
</FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>>
</FONT><BR><FONT size=2>> Faïza</FONT> <BR><FONT size=2>>
</FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>>
</FONT><BR><FONT size=2>> </FONT><BR><FONT size=2>> Pensez à
visiter le site BRGM sur.... <A
onclick="return top.js.OpenExtLink(window,event,this)"
href="http://www.brgm.fr/" target=_blank>http://www.brgm.fr/</A></FONT>
<BR><FONT size=2>> </FONT><BR><FONT size=2>> ***</FONT> <BR><FONT
size=2>> Le contenu de cet e-mail et de ses pièces jointes est destiné à
</FONT><BR><FONT size=2>> l'usage exclusif du</FONT> <BR><FONT size=2>>
(des) destinataire(s) expressément désigné(s) comme tel(s). En cas de
</FONT><BR><FONT size=2>> réception de cet e-mail par erreur, le
signaler à son expéditeur et </FONT><BR><FONT size=2>> ne pas en divulguer
le contenu.</FONT> <BR><FONT size=2>> L'absence de virus a été vérifié
à l'émission du message. Il convient </FONT><BR><FONT size=2>>
néanmoins de vérifier l'absence de corruption à sa réception.</FONT> <BR><FONT
size=2>> </FONT><BR><FONT size=2>> The contents of this email and
any attachments are confidential. They </FONT><BR><FONT size=2>> are
intended for the named recipient(s) only. If you have received
</FONT><BR><FONT size=2>> this email in error please notify the system
manager or the sender </FONT><BR><FONT size=2>> immediately and do
not disclose the contents to anyone or make copies. </FONT><BR><FONT
size=2>> eSafe scanned this email for viruses, vandals and malicious
content.</FONT> <BR><FONT size=2>> ***</FONT> <BR><FONT size=2>> Pensez
à visiter le site BRGM sur.... <A
onclick="return top.js.OpenExtLink(window,event,this)"
href="http://www.brgm.fr/" target=_blank>http://www.brgm.fr/</A></FONT>
<BR><FONT size=2>> </FONT><BR><FONT size=2>> ***</FONT> <BR><FONT
size=2>> Le contenu de cet e-mail et de ses pièces jointes est destiné à
</FONT><BR><FONT size=2>> l'usage exclusif du</FONT> <BR><FONT size=2>>
(des) destinataire(s) expressément désigné(s) comme tel(s). En cas de
</FONT><BR><FONT size=2>> réception de cet e-mail par erreur, le
signaler à son expéditeur et ne pas en divulguer le contenu.</FONT> <BR><FONT
size=2>> L'absence de virus a été vérifié à l'émission du message. Il
convient </FONT><BR><FONT size=2>> néanmoins de vérifier l'absence de
corruption à sa réception.</FONT> <BR><FONT size=2>> </FONT><BR><FONT
size=2>> The contents of this email and any attachments are confidential.
They </FONT><BR><FONT size=2>> are intended for the named recipient(s)
only. If you have received </FONT><BR><FONT size=2>> this email in error
please notify the system manager or the sender </FONT><BR><FONT
size=2>> immediately and do not disclose the contents to anyone or make
copies. </FONT><BR><FONT size=2>> eSafe scanned this email for viruses,
vandals and malicious content.</FONT> <BR><FONT size=2>> ***</FONT>
</P><FONT color=blue size=3></FONT><PRE><FONT color=blue size=3>Pensez à visiter le site BRGM sur.... <A onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.brgm.fr/" target=_blank>http://www.brgm.fr/
</A><BR><BR>***<BR>Le contenu de cet e-mail et de ses pièces jointes est destiné à l'usage exclusif du <BR>(des) destinataire(s) expressément désigné(s) comme tel(s). En cas de réception de cet <BR> e-mail par erreur, le signaler à son expéditeur et ne pas en divulguer le contenu.
<BR>L'absence de virus a été vérifié à l'émission du message. Il convient néanmoins de <BR>vérifier l'absence de corruption à sa réception.<BR><BR>The contents of this email and any attachments are confidential. They are intended for
<BR>the named recipient(s) only. If you have received this email in error please notify the <BR>system manager or the sender immediately and do not disclose the contents to <BR>anyone or make copies. eSafe scanned this email for viruses, vandals and malicious
<BR>content.<BR>***<BR></FONT></PRE></SPAN></DIV>
<DIV></DIV></DIV></BLOCKQUOTE></DIV><BR><BR clear=all><BR>-- <BR>"Failure has a
thousand explanations. Success doesn't need one" -- Sir Alec
Guiness
<FONT SIZE=3 COLOR=BLUE><PRE>
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.
***
</PRE></FONT>
</BODY></HTML>