<!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.2800.1611" name=GENERATOR></HEAD>
<BODY>
<DIV><FONT size=2>
<P>Hi all,</P>
<P>Thanks for the reply. As you suggested I used the following</P>
<P>PCFactorSetMatSolverPackage(pc,MAT_SOLVER_SUPERLU_DIST)</P>
<P>when setting PC type.</P>
<P>It works only in serial mode. For the parallel run, the results are
wrong.</P>
<P>Have you guys seen the same thing or is there something else I
overlooked?</P>
<P>The version of my superlu_dist is 2.3.</P>
<P>Thanks,</P>
<P>XG, RAVI</P></FONT></DIV>
<BLOCKQUOTE dir=ltr style="MARGIN-RIGHT: 0px">
<DIV class=OutlookMessageHeader dir=ltr align=left><FONT face=Tahoma
size=2>-----Original Message-----<BR><B>From:</B> Ravi Kannan
[mailto:rxk@cfdrc.com]<BR><B>Sent:</B> Wednesday, March 25, 2009 10:27
AM<BR><B>To:</B> PETSc users list<BR><B>Subject:</B> superlu_dist doesn't work
in peysc-3.0.0-p1<BR><BR></FONT></DIV>
<DIV><FONT face=Arial color=#0000ff size=2><SPAN
class=062342016-25032009>Hi,</SPAN></FONT></DIV>
<DIV><FONT face=Arial color=#0000ff size=2><SPAN
class=062342016-25032009></SPAN></FONT> </DIV>
<DIV><FONT face=Arial color=#0000ff size=2><SPAN
class=062342016-25032009>After I upgrade the petsc from 2.3.3 to 3.0.0, I have
made the change for the superlu from</SPAN></FONT></DIV>
<DIV><FONT face=Arial color=#0000ff size=2><SPAN
class=062342016-25032009>_ierr =
MatSetType(_A,MATSUPERLU_DIST) </SPAN></FONT></DIV>
<DIV><FONT face=Arial color=#0000ff size=2><SPAN
class=062342016-25032009>to</SPAN></FONT></DIV>
<DIV><FONT face=Arial color=#0000ff size=2>_ierr =
MatSetType(_A,MAT_SOLVER_SUPERLU_DIST)</FONT></DIV>
<DIV><FONT face=Arial color=#0000ff size=2></FONT> </DIV>
<DIV><SPAN class=062342016-25032009><FONT face=Arial color=#0000ff size=2>Is
this the only change I need to do?</FONT></SPAN></DIV>
<DIV><SPAN class=062342016-25032009><FONT face=Arial color=#0000ff
size=2></FONT></SPAN> </DIV>
<DIV><SPAN class=062342016-25032009><FONT face=Arial color=#0000ff
size=2>Ravi, X.G</FONT></SPAN></DIV>
<BLOCKQUOTE>
<DIV class=OutlookMessageHeader dir=ltr align=left><FONT face=Tahoma
size=2>-----Original Message-----<BR><B>From:</B>
petsc-users-bounces@mcs.anl.gov
[mailto:petsc-users-bounces@mcs.anl.gov]<B>On Behalf Of </B>Matthew
Knepley<BR><B>Sent:</B> Wednesday, March 25, 2009 7:08 AM<BR><B>To:</B>
PETSc users list<BR><B>Subject:</B> Re: Petsc parallel vectors with two
communicators<BR><BR></FONT></DIV>On Wed, Mar 25, 2009 at 1:05 AM, Khan,
Irfan <SPAN dir=ltr><<A
href="mailto:irfan.khan@gatech.edu">irfan.khan@gatech.edu</A>></SPAN>
wrote:<BR>
<DIV class=gmail_quote>
<BLOCKQUOTE class=gmail_quote
style="PADDING-LEFT: 1ex; MARGIN: 0pt 0pt 0pt 0.8ex; BORDER-LEFT: rgb(204,204,204) 1px solid">Hi<BR>Can
the petsc parallel vectors be used with two different communicators? For
instance, I have created two different communicators called FEA_Comm and
FSI_Comm. The total number of processes are x+y. FSI_Comm works on x+y but
FEA_Comm works only on x.<BR><BR>Now I am trying to create parallel
vectors a1 and a2 such that a1 has entries from x+y processes but a2 has
entries from only y processes.<BR><BR>After splitting the communicators I
assign PETSC_COMM_WORLD to FEA_Comm which works on only x processes.
Subsequently petsc is initialized (PetscInitialize()). But when the
parallel vectors are created, the processes hang.</BLOCKQUOTE>
<DIV><BR>PETSC_COMM_WORLD should encompass all processes you wish to use in
PETSc, so that means x+y. You can create Vec<BR>objects on subcommunicators,
like x.<BR><BR> Matt<BR> </DIV>
<BLOCKQUOTE class=gmail_quote
style="PADDING-LEFT: 1ex; MARGIN: 0pt 0pt 0pt 0.8ex; BORDER-LEFT: rgb(204,204,204) 1px solid"><BR>Any
suggestions will be helpful<BR><BR>Thankyou<BR>Irfan<BR>Graduate Research
Assistant<BR>Woodruff school of Mechanical Engineering<BR>Atlanta, GA
(30307)<BR></BLOCKQUOTE></DIV>-- <BR>What most experimenters take for
granted before they begin their experiments is infinitely more interesting
than any results to which their experiments lead.<BR>-- Norbert
Wiener<BR></BLOCKQUOTE></BLOCKQUOTE></BODY></HTML>