A cooperative conjugate gradient method for linear systems permitting efficient multi-thread implementation - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Article Dans Une Revue Computational & Applied Mathematics Année : 2017

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

Résumé

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
Origine : Fichiers produits par l'(les) auteur(s)

Dates et versions

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

Identifiants

Citer

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⟩
195 Consultations
763 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More