Neanderthal

行列と線形代数の超速い計算を可能にするclojureライブラリ

BLAS と LAPACK をベースにしている

チュートリアルなどはこちら

CPU(native engine)を使う場合でも Intel の MKL で提供されている so なり dll なりを所定の場所に配置する必要がある

code:clj

(dge 2 2 2 2 3 3 {:layout :row})

ベクトル/行列の各要素に alpha 倍する: ax

ベクトル/行列同士の各要素の積: axpy

逆行列の求め方？

code:clj

(def mat (dge 3 3 1 0 -1 3 -2 3 4 1 1 {:layout :row}))

(def lu (trf mat))

(mm (tri lu) mat)