superlu_dist doesn't work in peysc-3.0.0-p1
Ravi Kannan
rxk at cfdrc.com
Wed Mar 25 11:26:45 CDT 2009
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/20090325/d6fba4c8/attachment.htm>
More information about the petsc-users
mailing list