Recursive Algorithms for Dense Linear Algebra: The ReLAPACK Collection

02/22/2016
by   Elmar Peise, et al.
0

To exploit both memory locality and the full performance potential of highly tuned kernels, dense linear algebra libraries such as LAPACK commonly implement operations as blocked algorithms. However, to achieve next-to-optimal performance with such algorithms, significant tuning is required. On the other hand, recursive algorithms are virtually tuning free, and yet attain similar performance. In this paper, we first analyze and compare blocked and recursive algorithms in terms of performance, and then introduce ReLAPACK, an open-source library of recursive algorithms to seamlessly replace most of LAPACK's blocked algorithms. In many scenarios, ReLAPACK clearly outperforms reference LAPACK, and even improves upon the performance of optimizes libraries.

READ FULL TEXT

page 1

page 4

page 6

page 7

page 14

page 15

research
03/22/2021

Kokkos Kernels: Performance Portable Sparse/Dense Linear Algebra and Graph Kernels

As hardware architectures are evolving in the push towards exascale, dev...
research
04/08/2017

BLASFEO: basic linear algebra subroutines for embedded optimization

BLASFEO is a dense linear algebra library providing high-performance imp...
research
04/25/2021

Efficient Binary Decision Diagram Manipulation in External Memory

We follow up on the idea of Lars Arge to rephrase the Reduce and Apply a...
research
08/20/2021

On the Parallel I/O Optimality of Linear Algebra Kernels: Near-Optimal Matrix Factorizations

Matrix factorizations are among the most important building blocks of sc...
research
10/11/2017

Deriving Correct High-Performance Algorithms

Dijkstra observed that verifying correctness of a program is difficult a...
research
07/05/2022

FLOPs as a Discriminant for Dense Linear Algebra Algorithms

Expressions that involve matrices and vectors, known as linear algebra e...
research
05/02/2020

On infinite guarded recursive specifications in process algebra

In most presentations of ACP with guarded recursion, recursive specifica...

Please sign up or login with your details

Forgot password? Click here to reset