superlu_dist doesn't work in peysc-3.0.0-p1

Ravi Kannan rxk at
Wed Mar 25 11:26:45 CDT 2009


After I upgrade the petsc from 2.3.3 to 3.0.0, I have made the change for
the superlu from
_ierr = MatSetType(_A,MATSUPERLU_DIST)

Is this the only change I need to do?

Ravi, X.G
  -----Original Message-----
  From: petsc-users-bounces at
[mailto:petsc-users-bounces at]On Behalf Of Matthew Knepley
  Sent: Wednesday, March 25, 2009 7:08 AM
  To: PETSc users list
  Subject: Re: Petsc parallel vectors with two communicators

  On Wed, Mar 25, 2009 at 1:05 AM, Khan, Irfan <irfan.khan at>

    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.

    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.

    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.

  PETSC_COMM_WORLD should encompass all processes you wish to use in PETSc,
so that means x+y. You can create Vec
  objects on subcommunicators, like x.


    Any suggestions will be helpful

    Graduate Research Assistant
    Woodruff school of Mechanical Engineering
    Atlanta, GA (30307)

  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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <>

More information about the petsc-users mailing list