C ++中的慢矩阵求逆

Question

我目前正在尝试使用犰狳将Matlab代码转换为C ++。 我按照aramdillo文档将某些matlab代码转换为C ++。 但是，与matlab相比，性能令人失望。

在Matlab中，尺寸为625x625的矩阵A的求逆大约需要0.1秒，而在C ++中则需要3秒以上。

在C ++中，我都尝试过

solve()

以及

inv() 

我知道inv产生的结果不太准确的事实，因此我不喜欢使用它。 此外，我确实需要矩阵A的逆，因为我稍后将在算法中使用对角线元素。

产生这些结果的代码：

Matlab的

x=A\b
invA = A\eye(size(A))

C ++

arma::mat x = solve(A,b)
arma::mat invA = solve(A,eye(625,625))

我使用的版本：

C ++：

• Visual Studio 2013

• 犰狳8.300.1

• 英特尔MKL 2018.1.156

Matlab：

• Matlab 2016b

• 版本-blas

适用于Intel（R）64架构应用程序的Intel（R）数学内核库版本11.3.1产品内部版本20151021，CNR分支AVX2

• 版本-lapack

英特尔（R）数学内核库11.3.1版，适用于英特尔（R）64架构应用的产品版本20151021，CNR分支AVX2线性代数软件包版本3.5.0

有谁知道如何使用犰狳克服C ++中这种速度不足的问题？

Answer 1

您是否尝试过：将代码添加到#ARMA_USE_LAPACK中？ 这使您的代码可以使用.inv（）函数的更优化版本。 还要检查文档： http : //arma.sourceforge.net/docs.html