替代Eigen :: MatrixXd typedef

Question

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

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

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

将这些typedef切换到什么是最好的？如果我将它们保留原样会怎样？

Answer 1

只需根据Eigen自己的全局矩阵typedef定义自己的typedef。

如果您使用Eigen::MatrixXdtype元素填充它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;
}