c++ - ublas行列式のチュートリアル/例

Question

特定の行列演算を実装しようとしていますが、ublasライブラリの内部で迷子になっています。新しいublas行列式を実装する方法に関するチュートリアルや例などのリソースはありますか？

ありがとう

Answer 1

それが役立つかどうかはわかりませんが、uBlasの拡張に関するwikiページがここにあります。その式テンプレートのものは本当に私の心を吹き飛ばします。:)

Answer 2

私の提案は、行列式やublas内部について心配する必要がないように、新しい関数をテンプレート化することです。たとえば、独自の逆関数を記述したい場合は、それをテンプレートとして記述します。

template<typename MATRIX_IN, typename MATRIX_OUT>
bool invert_matrix(const MATRIX& A, MATRIX_OUT& B)
{
    // now use A(i,j), B(i,j), etc
    ...
}

Answer 3

それらが少しでも役立つならば、ここにいくつかの例もあります。