A cooperative conjugate gradient method for linear systems permitting efficient multi-thread implementation - Archive ouverte HAL Access content directly
Journal Articles Computational & Applied Mathematics Year : 2017

A cooperative conjugate gradient method for linear systems permitting efficient multi-thread implementation

(1) , (2, 3) , (4) , (5)
1
2
3
4
5

Abstract

This paper revisits, in a multi-thread context, the so-called mul-tiparameter or block conjugate gradient (BCG) methods, first proposed as sequential algorithms by O'Leary and Brezinski, for the solution of the linear system Ax = b, for an n-dimensional symmetric positive definite matrix A. Instead of the scalar parameters of the classical CG algorithm, which minimizes a scalar functional at each iteration, multiple descent and conjugate directions are updated simultaneously. Implementation involves the use of multiple threads and the algorithm is referred to as cooperative CG (CCG) in order to emphasize that each thread now uses information that comes from the other threads. It is shown that for a sufficiently large matrix dimension n, the use of an optimal number of threads results in a worst case flop count of O(n 7/3) in exact arithmetic. Numerical experiments on a multicore, multi-thread computer , for synthetic and real matrices, illustrate the theoretical results.
Fichier principal
Vignette du fichier
bbnp16b.pdf (382.13 Ko) Télécharger le fichier
Origin : Files produced by the author(s)

Dates and versions

hal-01558765 , version 1 (10-07-2017)

Identifiers

Cite

Amit A Bhaya, Pierre-Alexandre A Bliman, Guilherme Niedu, Fernando A Pazos. A cooperative conjugate gradient method for linear systems permitting efficient multi-thread implementation. Computational & Applied Mathematics, 2017, pp.1-28. ⟨10.1007/s40314-016-0416-7⟩. ⟨hal-01558765⟩
188 View
684 Download

Altmetric

Share

Gmail Facebook Twitter LinkedIn More