petsc-3.8.4 2018-03-24
Report Typos and Errors


Newton basis GMRES implementation with adaptive augmented eigenvectors The techniques used are best described in [1]. The contribution of this work is that it combines many of the previous work to reduce the amount of MPI messages and improve the robustness of the global approach by using deflation techniques. It has been successfully applied to a class of real and industrial problems. Please see [1] for numerical experiments and [2] for a description of these problems. There are many ongoing work that aim at avoiding (or minimizing) the communication in Krylov subspace methods. This code can be used as an experimental framework to combine several techniques in the particular case of GMRES. For instance, the computation of the shifts can be improved with techniques described in [3]. The orthogonalization technique can be replaced by TSQR [4]. The generation of the basis can be done using s-steps approaches[5].

Options Database Keys

+ -ksp_gmres_restart <restart> - the number of Krylov directions . -ksp_gmres_krylov_monitor - plot the Krylov space generated . -ksp_agmres_eigen <neig> - Number of eigenvalues to deflate (Number of vectors to augment) . -ksp_agmres_maxeigen <max_neig> - Maximum number of eigenvalues to deflate . -ksp_agmres_MinRatio <1> - Relaxation parameter in the adaptive strategy; smallest multiple of the remaining number of steps allowed . -ksp_agmres_MaxRatio <1> - Relaxation parameter in the adaptive strategy; Largest multiple of the remaining number of steps allowed . --ksp_agmres_DeflPrecond Apply deflation as a preconditioner, this is similar to DGMRES but it rather builds a Newton basis. This is an experimental option. . -ksp_dgmres_force <0, 1> - Force the deflation at each restart. . - There are many experimental parameters. Run with -help option to see the whole list

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

Developer Notes: This object is subclassed off of KSPDGMRES

Contributed by Desire NUENTSA WAKAM, INRIA <[email protected]> Inputs from Guy Atenekeng <[email protected]> and R.B. Sidje <[email protected]>


+ [1] D. Nuentsa Wakam and J. Erhel, Parallelism and robustness in GMRES with the Newton basis and the deflated restarting. Research report INRIA RR-7787, November 2011,, in revision for ETNA. . [2] D. NUENTSA WAKAM and F. PACULL, Memory Efficient Hybrid Algebraic Solvers for Linear Systems Arising from Compressible Flows, Computers and Fluids, In Press, . [3] B. Philippe and L. Reichel, On the generation of Krylov subspace bases, Applied Numerical Mathematics, 62(9), pp. 1171-1186, 2012 . [4] J. Demmel, L. Grigori, M. F. Hoemmen, and J. Langou, Communication-optimal parallel and sequential QR and LU factorizations, SIAM journal on Scientific Computing, 34(1), A206-A239, 2012 . [5] M. Mohiyuddin, M. Hoemmen, J. Demmel, and K. Yelick, Minimizing communication in sparse matrix solvers, in SC '09: Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis, New York, NY, USA, 2009, ACM, pp. 1154-1171. . Sidje, Roger B. Alternatives for parallel Krylov subspace basis computation. Numer. Linear Algebra Appl. 4 (1997), no. 4, 305-331

.seealso: KSPCreate(), KSPSetType(), KSPType (for list of available types), KSP, KSPDGMRES, KSPPGMRES, KSPGMRESSetRestart(), KSPGMRESSetHapTol(), KSPGMRESSetPreAllocateVectors(), KSPGMRESSetOrthogonalization(), KSPGMRESGetOrthogonalization(), KSPGMRESClassicalGramSchmidtOrthogonalization(), KSPGMRESModifiedGramSchmidtOrthogonalization(), KSPGMRESCGSRefinementType, KSPGMRESSetCGSRefinementType(), KSPGMRESGetCGSRefinementType(), KSPGMRESMonitorKrylov(), KSPSetPCSide()

Index of all KSP routines
Table of Contents for all manual pages
Index of all manual pages