LAPACK/BLAS など行列計算系ライブラリのメモ(C/C++ から使う用) #Eigen

背景

C/C++ アプリで行列計算をする必要がある(e.g. SVD など).
- 主に点群レジストレーション
Eigen でもいいが, 計算が速くなるなら専用のライブラリを使うのも検討したい.
- しかしなるべくクロスプラットフォームで, OSS(non-copyleft)がよくて, アーキテクチャや OS 依存は避けたい.
数値計算界隈はよくわからない

なんとなくの理解です. 基本は BLAS と LAPACK を押さえておけばよさそう.

BLAS : 基本となる行列計算ルーチン
LAPACK : SVD とか行列計算アルゴリズムが実装されている. BLAS を使っている. 通常 C/C++ アプリではだたい LAPACK の API を呼び出すのが多い... はず

基本 API は変わらない(fortran 時代)まま, みなさん BLAS/LAPACK を使っているようです(みんな昔から使っているので, あとから変えるのが難しいというのもありそう).

コードベース(source repo)はいろいろ派生があってわかりずらいです.

素(?) の BLAS : fortran で書かれている. 古い.
CBLAS : 基本的には f2c で C に変換したもの. OpenBLAS など, 最適化された BLAS が使えないときに, こちらをビルドして使うのが想定される. アセンブラは使っていないのでポータブルに動く(はず).
OpenBLAS : アーキテクチャごとに最適化された BLAS. アセンブラを使っている. arm もある. https://qiita.com/syoyo/items/6f15ddc8eeb260bda433
ATLAS : http://math-atlas.sourceforge.net/ BLAS API で BLAS を最適化したもの. LAPACK API もいくらかあるらしい. 最近は OpenBLAS が主流?
LAPACK : fortran で書かれている.
clapack : f2c で C に変換したもの. コンパイルは比較的楽. https://qiita.com/syoyo/items/09009ca77393cc0bd986
lapacke : https://www.netlib.org/lapack/lapacke.html f2c ではなく手で書いた LAPACK C interface っぽい. より C から呼び出しやすい API になっている? https://stackoverflow.com/questions/26875415/difference-between-lapacke-and-lapack
- 現在は lapack の repo https://github.com/Reference-LAPACK/lapack (本家?) に取り込まれている
ARPACK : https://www.caam.rice.edu/software/ARPACK/ 疎行列計算. f77 だったり f90 だったりいろいろ派生があるっぽい? LAPACK/BLAS を使っている模様
superlu : https://directory.fsf.org/wiki/SuperLU 疎行列計算. LAPACK/BLAS を使っている.

Eigen : LAPACK/BLAS を使うオプションがある. https://eigen.tuxfamily.org/dox/TopicUsingBlasLapack.html
- MPL or GPL なのでちょっとライセンスが微妙(Eigen のコードをいじらなければ MPL で実質 MIT/BSD 系として使えはする)
arnmadillo : http://arma.sourceforge.net/ ビルドに LAPACK/BLAS は必須. 疎行列の superlu, ARPACK は optional
- Apache 2.0 ライセンスなのがよい.
blaze: https://bitbucket.org/blaze-lib/blaze/src/master/ BSD ライセンスでよい.
- Eigen のように基本テンプレートでヘッダオンリー.
- 自前(?)で疎行列も実装しているっぽい
- C++14 or later にはなるが, 機能一通りあり, 開発結構アクティブである, 今後の本命かも

疎行列ライブラリはさらにいろいろあってよくわからない.

とりあえずは suitesparse と, superlu を押さえておけばよいか?

suitesparse はライセンスが混在していてつらい