c++ - lapack++ (2.5+) で行列式を計算する方法

Question

lapack ++ライブラリを使用して(非対称、二乗) LaMatGenDouble行列の行列式を計算する最良の(最速の)方法は何ですか?

Answer 1

行列式を計算する 1 つの方法は、LU 分解を使用することです。

  LaVectorLongInt pivots(A.cols());

  LUFactorizeIP(A, pivots);

  double detA = 1;
  for (int i = 0; i < A.cols(); ++i)
    detA *= A(i, i);

警告、Aは変更されるため、コピーを作成することをお勧めします。

Answer 2

lapack ++についてはわかりませんが、標準のlapackにはないはずです。確認してください。私の知る限り、lapack ++はマトリシアル操作自体を実装していませんが、他のものを使用しています。実際には、それらのいくつか（atlas、mkl（インテルマスカーネルライブラリ）など）を切り替えることができます。したがって、私の仮定は、lapack++にも行列式の操作があるということです。