如何使用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 ++中也有任何決定性操作。計算行列式