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 >
