Eigen::MatrixXd 类型定义的替换

Question

用具有long double 元素的向量和矩阵替换所有Eigen::MatrixXds 和Eigen::VectorXds 的最简单方法是什么？

我的代码中的每个基本浮点变量都是long double 类型。此外，每次我使用矩阵或向量时，我都会使用以下 typedef。

typedef Eigen::VectorXd Vec;
typedef Eigen::MatrixXd Mat;

将这些 typedef 切换为最好的方法是什么？如果我让它们保持原样会发生什么？

Answer 1

只需根据 Eigen 自己的 global matrix typedefs 定义自己的 typedef。

如果您使用Eigen::MatrixXd 并用long double 类型的元素填充它，这些值将缩小以适合矩阵的double 元素，这会导致精度损失，或者在最坏的情况下，溢出错误。然而，在许多架构上，双精度浮点运算是以 80 位扩展精度完成的，因此结果可能相同。你当然不应该依赖这个！如需更多信息，请参阅long double vs double。

#include <Eigen/Core>

typedef Eigen::Matrix< long double, Eigen::Dynamic, 1              > Vec;
typedef Eigen::Matrix< long double, Eigen::Dynamic, Eigen::Dynamic > Mat;

int main()
{
  long double ld = 2;

  Mat m(1,1);
  m(0,0) = ld;
}