Serinv: A Scalable Library for the Selected Inversion of Block-Tridiagonal with Arrowhead Matrices
Fuente:
arXiv
Saved in:
| Main Authors: | , , , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866917964934021120 |
|---|---|
| author | Maillou, Vincent Gaedke-Merzhaeuser, Lisa Ziogas, Alexandros Nikolaos Schenk, Olaf Luisier, Mathieu |
| author_facet | Maillou, Vincent Gaedke-Merzhaeuser, Lisa Ziogas, Alexandros Nikolaos Schenk, Olaf Luisier, Mathieu |
| contents | The inversion of structured sparse matrices is a key but computationally and memory-intensive operation in many scientific applications. There are cases, however, where only particular entries of the full inverse are required. This has motivated the development of so-called selected-inversion algorithms, capable of computing only specific elements of the full inverse. Currently, most of them are either shared-memory codes or limited to CPU implementations. Here, we introduce Serinv, a scalable library providing distributed, GPU-based algorithms for the selected inversion and Cholesky decomposition of positive-definite, block-tridiagonal arrowhead matrices. This matrix class is highly relevant in statistical climate modeling and materials science applications. The performance of Serinv is demonstrated on synthetic and real datasets from statistical air temperature prediction models. In our numerical tests, Serinv achieves 32.3% strong and 47.2% weak scaling efficiency and up to two orders of magnitude speedup over the sparse direct solvers PARDISO and MUMPS on 16 GPUs. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2503_17528 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Serinv: A Scalable Library for the Selected Inversion of Block-Tridiagonal with Arrowhead Matrices Maillou, Vincent Gaedke-Merzhaeuser, Lisa Ziogas, Alexandros Nikolaos Schenk, Olaf Luisier, Mathieu Distributed, Parallel, and Cluster Computing Numerical Analysis Performance G.4; F.2 The inversion of structured sparse matrices is a key but computationally and memory-intensive operation in many scientific applications. There are cases, however, where only particular entries of the full inverse are required. This has motivated the development of so-called selected-inversion algorithms, capable of computing only specific elements of the full inverse. Currently, most of them are either shared-memory codes or limited to CPU implementations. Here, we introduce Serinv, a scalable library providing distributed, GPU-based algorithms for the selected inversion and Cholesky decomposition of positive-definite, block-tridiagonal arrowhead matrices. This matrix class is highly relevant in statistical climate modeling and materials science applications. The performance of Serinv is demonstrated on synthetic and real datasets from statistical air temperature prediction models. In our numerical tests, Serinv achieves 32.3% strong and 47.2% weak scaling efficiency and up to two orders of magnitude speedup over the sparse direct solvers PARDISO and MUMPS on 16 GPUs. |
| title | Serinv: A Scalable Library for the Selected Inversion of Block-Tridiagonal with Arrowhead Matrices |
| topic | Distributed, Parallel, and Cluster Computing Numerical Analysis Performance G.4; F.2 |
| url | https://arxiv.org/abs/2503.17528 |