Parallel Selected Inversion of Block-Tridiagonal with Arrowhead Matrices

The inversion of structured sparse matrices is a fundamental yet computationally and memory-intensive task in many scientific applications, such as Bayesian statistical modeling and material science. In certain cases, only particular entries of the full inverse are required. This has motivated the d...

Celý popis

Uloženo v:
Podrobná bibliografie
Vydáno v:Proceedings / IEEE International Conference on Cluster Computing s. 1 - 12
Hlavní autoři: Maillou, Vincent, Gaedke-Merzhauser, Lisa, Ziogas, Alexandros Nikolaos, Schenk, Olaf, Luisier, Mathieu
Médium: Konferenční příspěvek
Jazyk:angličtina
Vydáno: IEEE 02.09.2025
Témata:
Computational modeling
Distributed algorithms
Graphics processing units
Libraries
Linear algebra
Materials science and technology
Parallel algorithms
Predictive models
Selected linear algebra
Sparse matrices
Supercomputers
Temperature distribution
ISSN:2168-9253
On-line přístup:Získat plný text
Tagy: Přidat tag
Žádné tagy, Buďte první, kdo vytvoří štítek k tomuto záznamu!
Popis
Shrnutí:The inversion of structured sparse matrices is a fundamental yet computationally and memory-intensive task in many scientific applications, such as Bayesian statistical modeling and material science. In certain cases, only particular entries of the full inverse are required. This has motivated the development of so-called selected inversion algorithms (SIA), capable of computing only specific elements of the full inverse. Currently, most SIA implementations are restricted to shared-/distributed-memory CPU architectures or to single GPUs. Here, we introduce novel numerical methods to perform the parallel selected inversion and Cholesky decomposition of positive-definite, block-tridiagonal with arrowhead matrices. A distributed memory, GPU-accelerated implementation of our approach is presented and integrated into the structured solver library Serinv. We demonstrate its performance on synthetic and real datasets from statistical air temperature prediction models and achieve CPU (GPU) speedups of up to 2.6 \times(71.4 \times) over the SIA of the PARDISO library and up to 14 \times(380.9 \times) over the MUMPS library, when scaling to 16 processes.
ISSN:2168-9253
DOI:10.1109/CLUSTER59342.2025.11186484