Towards a Reproducible Solution of Linear Systems

Abstract : Solving a linear system of equations is an important underlying part of numerous scientific applications. In this article, we address the issue of non-deterministic and, therefore, non-reproducible solution of linear systems and propose an approach to ensure its reproducibility. Our approach is based on the hierarchical and modular structure of linear algebra algorithms. Consequently, we divide computations into smaller logical blocks – such as a blocked LU factorization, triangular system solve, and matrix-matrix multiplication – and ensure their reproducible results. In this manner, we also split the blocked LU factorization into the unblocked LU and the BLAS-3 routines; the former is built on top of scaling a vector and outer product of two vectors routines from BLAS-1 and-2, accordingly. In this work, our focus is on constructing these building blocks that eventually lead, as we will prove, to the reproducible solution of linear systems.
Complete list of metadatas

Cited literature [4 references]  Display  Hide  Download

https://hal.archives-ouvertes.fr/hal-01633980
Contributor : Roman Iakymchuk <>
Submitted on : Monday, November 13, 2017 - 3:47:22 PM
Last modification on : Thursday, March 21, 2019 - 2:43:10 PM
Long-term archiving on : Wednesday, February 14, 2018 - 4:14:14 PM

File

reproducible-solution-linear.p...
Files produced by the author(s)

Identifiers

  • HAL Id : hal-01633980, version 1

Citation

Roman Iakymchuk, Stef Graillat, David Defour, Erwin Laure, Enrique Quintana-Ortí. Towards a Reproducible Solution of Linear Systems. Supercomputing Conference 2017-Computational Reproducibility at Exascale Workshop, Nov 2017, Denver, United States. ⟨hal-01633980⟩

Share

Metrics

Record views

387

Files downloads

83