superlu_dist doesn't work in peysc-3.0.0-p1

[Hong Zhang]

Comment out
//PCFactorSetMatSolverPackage(pc,MAT_SOLVER_SUPERLU_DIST)

and use the runtime option
'-pc_type lu -pc_factor_mat_solver_package superlu_dist'
This option is well tested.

Let us know what you get.

Hong

On Tue, 31 Mar 2009, Ravi Kannan wrote:

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