[Libmesh-users] PetscDMWrapper functionality

[Salazar De T. M. <sal...@ll...>]

Hello,

I just wanted to ask if this capability (geometric multigrid in libMesh) is working properly before I dig deeper in the possible parameters for the PETSc mg. I’m trying to solve a linear elasticity problem with the settings:

-ksp_type cg
-pc_type mg
-pc_mg_levels 3
-pc_mg_galerkin
-pc_mg_type full

But I get a lof of iterations (>1000) in the linear solver. I do pass the options –use_petsc_dm and –node_major_dofs. I noticed this capability is still very new and that the number of levels seems to be limited to 1. Does that mean that I should not use it at all?

Thanks
Miguel


Miguel A. Salazar de Troya
Graduate Scholar, Lawrence Livermore National Laboratory
Ph.D. Candidate, University of Illinois at Urbana-Champaign
B141
Rm: 1085-5
Ph: 1(925) 422-6411