Abstract:The Conjugate Gradient (CG) method has been widely used to solve linear systems of the type when the matrix A is symmetric positive definitive. In the context of the implementation of CG as an iterative algorithm one important bottleneck is in the matrix vector multiplication which has the disadvantages of a low ratio of computation memory accesses. To improve this ratio several strategies have been researched ranging from new paradigm machine architectures to optimizing the cache data usage through the usage of a matrix multivector multiplications Kernel. The latter in the CG context yields the Block Conjugate Gradient (BCG) originally used to solve linear systems with multiple right -hand -side. A main drawback of BCG for the case of a single right hand side is the determination of an optimal block size. Moreover, if an adequate block size is not chosen the effect on the rate of convergence of the CG can be negative and consequently its rate of convergence deteriorates. In this work we compare strategies for determining adaptively the optimal block size in the implementation of BCG considering a single right hand side seeking a trade-off between convergence rate and computational cost. The strategies analyzed arc based on the Rayleigh Values associated to the search subspace, the associated Ritz Values, and the decay rate of the residual measured locally. Numerical experiments on real matrices arc presented. It is observed that the Ritz Values strategy shows better results than the other strategies in terms of decreasing the computational cost. These results encourage further analysis and design of adaptive block algorithms.

PCG Sparse Algorithm for Close-Range Block Bundle Adjustment

A scalable parallel preconditioned conjugate gradient method for bundle adjustment

Bundle block adjustment of large-scale remote sensing data with Block-based Sparse Matrix Compression combined with Preconditioned Conjugate Gradient

Restrictive Preconditioners for Conjugate Gradient Methods for Symmetric Positive Definite Linear Systems

CUDA-based PCG algorithm optimization for a large sparse matrix

On a type of matrix splitting preconditioners for a class of block two-by-two linear systems.

Practical Restrictively Preconditioned Conjugate Gradient Methods for a Class of Block Two-by-Two Linear Systems

A New Projected Variant of the Deflated Block Conjugate Gradient Method

RESTRICTIVELY PRECONDITIONED CONJUGATE GRADIENT METHOD FOR A SERIES OF CONSTANTLY AUGMENTED LEAST SQUARES PROBLEMS

Square Root Bundle Adjustment for Large-Scale Reconstruction

A Non-Krylov subspace Method for Solving Large and Sparse Linear System of Equations

A Method for Solving Sparse Linear Equations of Power Systems Based on GPU

Efficient Spectral Graph Sparsification Via Krylov-Subspace Based Spectral Perturbation Analysis

A Robust Algebraic Domain Decomposition Preconditioner for Sparse Normal Equations

A sketch-and-project method for solving the matrix equation AXB = C

BCCB preconditioners for solving linear systems from delay differential equations

Collective Communication Optimization for Solving Linear Algebraic Equations

Varying The Block Size On Block Conjugate Gradient: Comparison Of Strategies

A Refined Block Conjugate Gradient Method for Symmetric Eigenvalue Problems

A Block Coordinate Descent Approach for Sparse Principal Component Analysis

A power sparse approximate inverse preconditioning procedure for large sparse linear systems