如何用lapack++ (2.5+)计算行列式矩阵

Question

用lapack++库计算(非对称的，平方的) LaMatGenDouble矩阵的行列式的最佳(最快)方法是什么？

Answer 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中没有，check。据我所知，lapack++本身没有实现矩阵运算，而是使用了别人的，实际上你可以在其中的几个之间切换(atlas，mkl (英特尔数学内核库)等等)。因此，我的假设是lapack++中也有行列式运算。