superlu_dist doesn't work in peysc-3.0.0-p1

[Ravi Kannan]

Hi all,

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

PCFactorSetMatSolverPackage(pc,MAT_SOLVER_SUPERLU_DIST)

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.

Thanks,

XG, RAVI

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


  Hi,

  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)
  to
  _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 mcs.anl.gov
[mailto:petsc-users-bounces at mcs.anl.gov]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 gatech.edu>
wrote:

      Hi
      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.

      Matt


      Any suggestions will be helpful

      Thankyou
      Irfan
      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: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090331/f7b36e6c/attachment.htm>