Just to follow up this point: >> I would like to use a parallel direct solver from petsc which I have >> compiled with full set of solvers, umfpack, superlu etc.. but the >> command >> line argument -mat_type superlu etc. all throw an error in petsc >> 3.1.4 as >> unrecognized, so what is the command line way to use a parallel sparse >> direct solver in libmesh? >
I regularly use the command line options: -ksp_type preonly -pc_type lu -pc_factor_mat_solver_package umfpack to use a direct solver with libMesh/PETSc (this works for me with petsc 3.0.0-p12). To change the solver, you can change umfpack to mumps or superlu or whatever. Dave ------------------------------------------------------------------------------ The ultimate all-in-one performance toolkit: Intel(R) Parallel Studio XE: Pinpoint memory and threading errors before they happen. Find and fix more than 250 security defects in the development cycle. Locate bottlenecks in serial and parallel code that limit performance. http://p.sf.net/sfu/intel-dev2devfeb _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
