thank you, John. I need to do in the code. If I get ksp object from the PetscLinearSolver object, it might not help me, because the ::solve() method of the PetscLinearSolver will override my changes. Am I right? Michael.
On 07/20/2016 02:28 PM, John Peterson wrote: > > > On Wed, Jul 20, 2016 at 12:20 PM, Michael Povolotskyi > <mpovo...@purdue.edu <mailto:mpovo...@purdue.edu>> wrote: > > Dear Libmesh developers/users, > I have the following question: > > How to configure the linear solver of a linear system to use MUMPS > package available though PETSc? > In PETSc I would do something like this: > > ierr = KSPSetType(ksp,KSPPREONLY); > ierr = PCSetType(pc,PCLU); > ierr = PCFactorSetMatSolverPackage(pc,MATSOLVERMUMPS); > > > You can do it from the command line: > > -ksp_type preonly -pc_type lu -pc_factor_mat_solver_package mumps > > If you want to do it in code, you can always call ksp() and pc() to > get the corresponding objects from a PetscLinearSolver object, and > make PETSc calls directly on them. > > -- > John -- Michael Povolotskyi, PhD Research Assistant Professor Network for Computational Nanotechnology Hall for Discover and Learning Research, Room 441 207 South Martin Jischke Drive West Lafayette, IN 47907 Phone (765) 4949396 ------------------------------------------------------------------------------ What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic patterns at an interface-level. Reveals which users, apps, and protocols are consuming the most bandwidth. Provides multi-vendor support for NetFlow, J-Flow, sFlow and other flows. Make informed decisions using capacity planning reports.http://sdm.link/zohodev2dev _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
