Question

在例如数据之间复制数据很容易。 Eigen::VectorXd和std::vector<double>或std::vector<Eigen::Vector3d>，例如

std::vector<Eigen::Vector3> vec1(10, {0,0,0});
Eigen::VectorXd vec2(30);
VectorXd::Map(&vec2[0], vec1.size()) = vec1;

（请参阅例如https://stackoverflow.com/a/26094708/4069571或https://stackoverflow.com/a/21560121/4069571）

此外，可以从Matrix块/列/ ...创建Eigen::Ref<VectorXd>，例如

MatrixXd mat(10,10);
Eigen::Ref<VectorXd> vec = mat.col(0);

问题

是否可以在不首先复制数据的情况下从Eigen::Ref<VectorXd>甚至std::vector<double>创建std::vector<Eigen::Vector3d>？

Answer 1

我试过，它实际上正如我在评论中描述的那样，通过首先映射然后将其包装为Eigen::Ref对象。这里通过谷歌测试显示。

void processVector(Eigen::Ref<Eigen::VectorXd> refVec) {
  size_t size = refVec.size();
  ASSERT_TRUE(10 == size);
  std::cout << "Sum before change: " << refVec.sum(); // output is 50 = 10 * 5.0
  refVec(0) = 10.0; // for a sum of 55
  std::cout << "Sum after change: " << refVec.sum() << std::endl;
}

TEST(testEigenRef, onStdVector) {
  std::vector<double> v10(10, 5.0);
  Eigen::Map<Eigen::VectorXd> mPtr(&v10[0], 10);
  processVector(mPtr);
  // confirm that no copy is made and std::vector is changed as well
  std::cout << "Std vec[0]: " << v10[0] << std::endl; // output is 10.0
}

在第二次编辑后更精细。现在我对Eigen::Ref进行了谷歌单元测试（谢谢）。希望这会有所帮助。

从std :: vector创建Eigen :: Ref

1 个答案: