如何使用lapack ++（2.5+）计算行列式矩阵++（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); 

警告，一个会发生变化，所以制作一份拷贝大概建议。

Answer 2

我不知道lapack ++，但我确定没有标准的lapack，check。据我所知，lapack ++没有实现本身的操作，而是使用别人的操作，实际上你可以在其中的几个（atlas，mkl（intel数学内核库）等等）之间切换。因此，我的假设是，在lapack ++中也有任何决定性操作。计算行列式