ReLAPACK offers a collection of recursive algorithms for many of LAPACK's compute kernels. Since it preserves LAPACK's established interfaces, ReLAPACK integrates effortlessly into existing application codes. ReLAPACK's routines not only outperform the reference LAPACK but also improve upon the performance of tuned implementations, such as OpenBLAS and MKL.
For a detailed list of covered operations and an overview of operations to which recursion is not efficiently applicable, see coverage.md.
To compile with the default configuration, simply run make to create the
library librelapack.a.
Note that to link with MKL, you currently need to set the flag
COMPLEX_FUNCTIONS_AS_ROUTINES to 1 to avoid problems in ctrsyl and
ztrsyl. For further configuration options see config.md.
ReLAPACK builds on top of BLAS and unblocked kernels from LAPACK. There are many optimized and machine specific implementations of these libraries, which are commonly provided by hardware vendors or available as open source (e.g., OpenBLAS).
ReLAPACK's test suite compares its routines numerically with LAPACK's
counterparts. To set up the tests (located in test/) you need to specify
link flags for BLAS and LAPACK (version 3.5.0 or newer) in make.inc; then
make test runs the tests. For details on the performed tests, see
test/README.md.
Since ReLAPACK replaces parts of LAPACK, any LAPACK example involving the
covered routines applies directly to ReLAPACK. A few separate examples are
given in examples/. For details, see examples/README.md.
When referencing ReLAPACK, please cite the preprint of the paper Recursive Algorithms for Dense Linear Algebra: The ReLAPACK Collection:
@article{relapack,
author = {Elmar Peise and Paolo Bientinesi},
title = {Recursive Algorithms for Dense Linear Algebra: The ReLAPACK Collection},
journal = {CoRR},
volume = {abs/1602.06763},
year = {2016},
url = {http://arxiv.org/abs/1602.06763},
}