This article is cited in 1 paper

Some ways of parallel implementation of the conjugate gradient method with an implicit factorized preconditioner

O. Yu. Milyukova

Keldysh Institute of Applied Mathematics of RAS

Abstract: The paper considers two ways of using MPI and MPI+OpenMP technologies for constructing and inverting the preconditioner of an incomplete triangular Cholesky decomposition without filling IC(0) for solving systems of linear algebraic equations with an arbitrary symmetric positive definite matrix. They differ in the way in which the preconditioning matrix IC(0) is computed. Methods of using MPI and MPI+OpenMP technologies are based on the use of grid node orderings consistent with the division of the calculation area. Comparative timing results for the MPI+OpenMP and MPI implementations of the proposed preconditioning used with the conjugate gradient method for a model problems and the sparse matrix collections SuiteSpars e as well as comparing the time of solving these problems using two methods of using MPI and MPI + OpenMP technology are presented.

Keywords: incomplete Cholesky factorization, Domain Decomposition ordering, parallel preconditioning, conjugate gradient method.

Received: 22.08.2023
Revised: 22.08.2023
Accepted: 16.10.2023

DOI: 10.20948/mm-2024-02-10

 English version: Mathematical Models and Computer Simulations, 2024, 16:4, 638–653