Skip to Main content Skip to Navigation
Journal articles

Resolution of a large number of small random symmetric linear systems in single precision arithmetic on GPUs

Abstract : This paper focuses on the resolution of a large number of small symmetric linear systems and its parallel implementation on single precision on GPUs. The computations involved by each linear system are independent from the others and the number of unknowns does not exceed 64. For this purpose, we present the adaptation to our context of largely used methods that include: LDLt, House-holder reduction to a tridiagonal matrix, parallel cyclic reduction that is not a power of two and the divide and conquer algorithm for tridiagonal eigenprob-lems. We not only detail the implementation and optimization of each method but we also compare the sustainability of each solution and its performance which include both parallel complexity and cache memory occupation. In the context of solving a large number of small random linear systems on GPU with no information about their conditioning, we show that the best strategy seems to be the use of Householder tridiagonalization + PCR followed if necessary by a divide & conquer diagonalization.
Complete list of metadatas

Cited literature [24 references]  Display  Hide  Download

https://hal.archives-ouvertes.fr/hal-01295549
Contributor : Lokman Abbas-Turki <>
Submitted on : Thursday, March 31, 2016 - 11:43:06 AM
Last modification on : Friday, March 27, 2020 - 3:56:17 AM
Document(s) archivé(s) le : Friday, July 1, 2016 - 12:20:32 PM

File

LokmanStef.pdf
Files produced by the author(s)

Identifiers

Citation

Lokman A. Abbas-Turki, Stef Graillat. Resolution of a large number of small random symmetric linear systems in single precision arithmetic on GPUs. Journal of Supercomputing, Springer Verlag, 2017, 73 (4), pp.1360-1386. ⟨10.1007/s11227-016-1813-9⟩. ⟨hal-01295549⟩

Share

Metrics

Record views

447

Files downloads

530