Abstract
The linear equations that arise in interior methods for constrained optimization are sparse symmetric indefinite, and they become extremely ill-conditioned as the interior method converges. These linear systems present a challenge for existing solver frameworks based on sparse LU or LDLT decompositions. We benchmark five well known direct linear solver packages on CPU- and GPU-based hardware, using matrices extracted from power grid optimization problems. The achieved solution accuracy varies greatly among the packages. None of the tested packages delivers significant GPU acceleration for our test cases. For completeness of the comparison we include results for MA57, which is one of the most efficient and reliable CPU solvers for this class of problem.
Original language | American English |
---|---|
Article number | 102870 |
Number of pages | 9 |
Journal | Parallel Computing |
Volume | 111 |
DOIs | |
State | Published - Jul 2022 |
Bibliographical note
Publisher Copyright:© 2021
NREL Publication Number
- NREL/JA-2C00-81949
Keywords
- GPU
- Grid optimization
- Solvers
- Sparse linear equations