petsc-3.7.1 2016-05-15
Report Typos and Errors

KSPDGMRES

Implements the deflated GMRES as defined in [1,2]. In this implementation, the adaptive strategy allows to switch to the deflated GMRES when the stagnation occurs.

Options Database Keys

GMRES Options (inherited)

-ksp_gmres_restart <restart> - the number of Krylov directions to orthogonalize against Many br
-ksp_gmres_haptol <tol> - sets the tolerance for "happy ending" (exact convergence) Many br
-ksp_gmres_preallocate - preallocate all the Krylov search directions initially (otherwise groups of Many brvectors are allocated as needed) Many br
-ksp_gmres_classicalgramschmidt - use classical (unmodified) Gram-Schmidt to orthogonalize against the Krylov space (fast) (the default) Many br
-ksp_gmres_modifiedgramschmidt - use modified Gram-Schmidt in the orthogonalization (more stable, but slower) Many br
-ksp_gmres_cgs_refinement_type <never,ifneeded,always> - determine if iterative refinement is used to increase the Many brstability of the classical Gram-Schmidt orthogonalization. Many br
-ksp_gmres_krylov_monitor - plot the Krylov space generated Many br

DGMRES Options Database Keys

-ksp_dgmres_eigen <neig> - number of smallest eigenvalues to extract at each restart Many br
-ksp_dgmres_max_eigen <max_neig> - maximum number of eigenvalues that can be extracted during the iterative Many brprocess Many br
-ksp_dgmres_force - use the deflation at each restart; switch off the adaptive strategy. Many br
-ksp_dgmres_view_deflation_vecs <viewerspec> - View the deflation vectors, where viewerspec is a key that can be Many brparsed by PetscOptionsGetViewer(). If neig > 1, viewerspec should Many brend with ":append". No vectors will be viewed if the adaptive Many brstrategy chooses not to deflate, so -ksp_dgmres_force should also Many brbe given. Many brThe deflation vectors span a subspace that may be a good Many brapproximation of the subspace of smallest eigenvectors of the Many brpreconditioned operator, so this option can aid in understanding Many brthe performance of a preconditioner. Many br

Many br

Notes: Left and right preconditioning are supported, but not symmetric preconditioning. Complex arithmetic is not yet supported Many br

References

1. - J. Erhel, K. Burrage and B. Pohl, Restarted GMRES preconditioned by deflation,J. Computational and Applied Mathematics, 69(1996). Many br
2. - D. NUENTSA WAKAM and F. PACULL, Memory Efficient Hybrid Algebraic Solvers for Linear Systems Arising from Compressible Flows, Computers and Fluids, Many brIn Press, http://dx.doi.org/10.1016/j.compfluid.2012.03.023 Many br

Contributed by: Desire NUENTSA WAKAM,INRIA Many br

.seealso: KSPCreate(), KSPSetType(), KSPType (for list of available types), KSP, KSPFGMRES, KSPLGMRES, Many brKSPGMRESSetRestart(), KSPGMRESSetHapTol(), KSPGMRESSetPreAllocateVectors(), KSPGMRESSetOrthogonalization(), KSPGMRESGetOrthogonalization(), Many brKSPGMRESClassicalGramSchmidtOrthogonalization(), KSPGMRESModifiedGramSchmidtOrthogonalization(), Many brKSPGMRESCGSRefinementType, KSPGMRESSetCGSRefinementType(), KSPGMRESGetCGSRefinementType(), KSPGMRESMonitorKrylov(), KSPSetPCSide() Many br

Level:beginner
Location:
src/ksp/ksp/impls/gmres/dgmres/dgmres.c
Index of all KSP routines
Table of Contents for all manual pages
Index of all manual pages