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>
