为什么 Eigen's Cholesky 分解在我的 Ubuntu 上很慢

Why is Eigen's Cholesky decomposition very slow on my Ubuntu

本文关键字：我的 Ubuntu 分解 Cholesky Eigen 为什么      更新时间：2023-10-16

我在核心i5机器上使用Ubuntu 14.04。以下代码在不同的计算机上进行了测试，运行时间非常快(3258 毫秒(，但在我的系统上需要 112921 毫秒。我使用 g++ 作为编译器。

#include <iostream>
#include <chrono>
#include <Eigen/Core>
#include <Eigen/Cholesky>
using namespace std;
using namespace std::chrono;
using namespace Eigen;

int main()
{
const MatrixXd::Index size = 4200;
MatrixXd m = MatrixXd::Random(size, size);
m = (m + m.transpose()) / 2.0 + 10000 * MatrixXd::Identity(size, size);
LLT<MatrixXd> llt;
auto start = high_resolution_clock::now();
llt.compute(m);
if (llt.info() != Success)
    cout << "Cholesky decomposition failed!" << endl;
auto stop = high_resolution_clock::now();
cout << "Cholesky decomposition in "
     << duration_cast<milliseconds>(stop - start).count()
     << " ms." << endl;
return 0;
}

当您在调试模式下编译程序(例如g++ -std=c++11(时，很多事情都会在后台发生。回想一下发生崩溃时可以看到的信息。所有这些信息和家政服务都需要时间。

因此，当您的程序经过测试并且运行良好时，您不再需要内务管理，您希望性能发挥作用。

使用优化标志，如 O1、O2 或 O3(还有更多(。使用哪一个？很大程度上取决于您的应用程序，因此请至少尝试这三个选项。

注意：请确保使用大写的 Omikron，否则(如果您键入小写的 omikron 或零(，您会认为您在启用优化的情况下编译了代码，但它仍然会在调试模式下编译。

采购 GCC 教程：

生成 gdb 的符号信息和许多警告消息。

g++ -g -Wall myprog.C -o myprog

在带有警告的 Solaris 机器上生成优化的代码。-O 是大写的 o，而不是数字 0！

g++ -Wall -O1 -mv8 myprog.C -o myprog

我假设您是在调试模式下编译而没有优化的。

如果您使用的是 CMake，请使用 cmake -DCMAKE_BUILD_TYPE=Release 或 cmake -DCMAKE_BUILD_TYPE=ReleaseWithDebugInfo 进行构建。否则，请尝试将优化标志-O3添加到 g++ 编译器。

相关文章：