superlu_dist doesn't work in peysc-3.0.0-p1

Ravi Kannan rxk at
Tue Mar 31 10:41:43 CDT 2009

Hi all,

Thanks for the reply. As you suggested I used the following


when setting PC type.

It works only in serial mode. For the parallel run, the results are wrong.

Have you guys seen the same thing or is there something else I overlooked?

The version of my superlu_dist is 2.3.



  -----Original Message-----
  From: Ravi Kannan [mailto:rxk at]
  Sent: Wednesday, March 25, 2009 10:27 AM
  To: PETSc users list
  Subject: superlu_dist doesn't work in peysc-3.0.0-p1


  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)
  _ierr = MatSetType(_A,MAT_SOLVER_SUPERLU_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