Estimating the parallelism in the solution of sparse triangular linear systems

Saved in:
Bibliographic Details
Title: Estimating the parallelism in the solution of sparse triangular linear systems
Authors: Eduardo Martínez González, Pablo Ezzatti, Ernesto Dufrechou
Source: 2020 39th International Conference of the Chilean Computer Science Society (SCCC). :1-8
Publisher Information: IEEE, 2020.
Publication Year: 2020
Subject Terms: 0202 electrical engineering, electronic engineering, information engineering, 02 engineering and technology
Description: Most direct methods to solve sparse linear systems, as well as preconditioners for iterative methods, factorize the problem to operate on sparse triangular matrices. The solution of sparse triangular linear systems is then one of the fundamental building blocks of numerical methods, which has motivated extensive research dedicated to achieve efficient algorithms for this operation in parallel hardware platforms. However, as it often occurs with sparse problems, the parallel performance of each method depends heavily on the nonzero pattern of the matrix. In this sense, observing these characteristics can allow predicting which method is better suited for each problem.In the case of sparse triangular matrices, one of the most important constrains to parallelism, is the amount of data dependencies during forward or backward substitution. This is related to the number of level-sets, i.e. groups of independent rows in the matrix, but the high cost of computing this number makes its utilization impractical.In this work, we propose different strategies to approximate the number of level-sets through inexpensive procedures, and provide implementations of the heuristics for CPU and GPU.
Document Type: Article
DOI: 10.1109/sccc51225.2020.9281169
Access URL: http://dblp.uni-trier.de/db/conf/sccc/sccc2020.html#GonzalezDE20
https://doi.org/10.1109/SCCC51225.2020.9281169
http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber=9281169
http://xplorestaging.ieee.org/ielx7/9281144/9281146/09281169.pdf?arnumber=9281169
Rights: IEEE Copyright
Accession Number: edsair.doi.dedup.....ff4dfd82ed8fa6c9e4ffac4800a5ba49
Database: OpenAIRE
Description
Abstract:Most direct methods to solve sparse linear systems, as well as preconditioners for iterative methods, factorize the problem to operate on sparse triangular matrices. The solution of sparse triangular linear systems is then one of the fundamental building blocks of numerical methods, which has motivated extensive research dedicated to achieve efficient algorithms for this operation in parallel hardware platforms. However, as it often occurs with sparse problems, the parallel performance of each method depends heavily on the nonzero pattern of the matrix. In this sense, observing these characteristics can allow predicting which method is better suited for each problem.In the case of sparse triangular matrices, one of the most important constrains to parallelism, is the amount of data dependencies during forward or backward substitution. This is related to the number of level-sets, i.e. groups of independent rows in the matrix, but the high cost of computing this number makes its utilization impractical.In this work, we propose different strategies to approximate the number of level-sets through inexpensive procedures, and provide implementations of the heuristics for CPU and GPU.
DOI:10.1109/sccc51225.2020.9281169