KSPPIPECG2#
Pipelined conjugate gradient method with a single non-blocking reduction per two iterations [TV20]. Pipelined Krylov Methods
Notes#
This method has only a single non-blocking reduction per two iterations, compared to 2 blocking for standard KSPCG. The
non-blocking reduction is overlapped by two matrix-vector products and two preconditioner applications.
The solver has a two-step inner iteration, each of which computes the solution and updates the residual norm.
Hence the values from KSPGetResidualHistory() and KSPGetIterationNumber() will differ.
MPI configuration may be necessary for reductions to make asynchronous progress, which is important for performance of pipelined methods. See What steps are necessary to make the pipelined solvers execute efficiently?
Developer Note#
The implementation code contains a good amount of hand-tuned fusion of multiple inner products and similar computations on multiple vectors
Contributed by#
Manasi Tiwari, Computational and Data Sciences, Indian Institute of Science, Bangalore
References#
Manasi Tiwari and Sathish Vadhiyar. Pipelined preconditioned conjugate gradient methods for distributed memory systems. In 2020 IEEE 27th International Conference on High Performance Computing, Data, and Analytics (HiPC), 151–160. IEEE, 2020.
See Also#
KSP: Linear System Solvers, What steps are necessary to make the pipelined solvers execute efficiently?, Pipelined Krylov Methods, KSPCreate(), KSPSetType(), KSPCG, KSPPIPECG, KSPGROPPCG
Level#
intermediate
Location#
src/ksp/ksp/impls/cg/pipecg2/pipecg2.c
Index of all KSP routines
Table of Contents for all manual pages
Index of all manual pages