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

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

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

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

解决方法:

只需根据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;
}