Back to Search Start Over

Parallel Cholesky Factorization for Banded Matrices using OpenMP Tasks

Authors :
Liu, Felix
Fredriksson, Albin
Markidis, Stefano
Publication Year :
2023

Abstract

Cholesky factorization is a widely used method for solving linear systems involving symmetric, positive-definite matrices, and can be an attractive choice in applications where a high degree of numerical stability is needed. One such application is numerical optimization, where direct methods for solving linear systems are widely used and often a significant performance bottleneck. An example where this is the case, and the specific type of optimization problem motivating this work, is radiation therapy treatment planning, where numerical optimization is used to create individual treatment plans for patients. To address this bottleneck, we propose a task-based multi-threaded method for Cholesky factorization of banded matrices with medium-sized bands. We implement our algorithm using OpenMP tasks and compare our performance with state-of-the-art libraries such as Intel MKL. Our performance measurements show a performance that is on par or better than Intel MKL (up to ~26%) for a wide range of matrix bandwidths on two different Intel CPU systems.

Details

Database :
arXiv
Publication Type :
Report
Accession number :
edsarx.2305.04635
Document Type :
Working Paper