如何将Eigen::Matrix映射到std::V载体&;lt;Eigen::Vectoramp;gt;?
例如,如果我有一个大小为10列3行的
Eigen::MatrixXd
,我如何将其别名为Eigen::Vector3d
的10个元素的std::vector
?
当我说别名时,我是指使用相同的内存块而不进行复制。
我知道我可以通过如下方式进行反向映射:
std::vector<Vector3d> v(10);
...
Map<Matrix<double,3,Dynamic> > m(v.data().data(), 3, 10);
相反,如果我有一个Eigen::Matrix,我试图将它转换为Eigen::VECTOR的向量,但以下代码行编译失败
Eigen::Matrix2Xd points;
...
std::vector<Eigen::Vector2d> v_points(points.data()
Eigen::aligned_allocator<Eigen::vector2d>(points_data()))
解决方案
可以使用如下所述的自定义分配器:Is it possible to initialize std::vector over already allocated memory?
为此,您必须将其重新解释为Vector3d*
:
Vector3d* buf = reinterpret_cast<Vector3d*>(mat.data());
然后将其传递给您的自定义分配器:
std::vector<Vector3d, PreAllocator<Vector3d>> my_vec(10, PreAllocator<Vector3d>(buf, 10));
另一种选择是将其包装在一些类似std::vector
的视图中,例如gsl::span
:
gsl::span<Vector3d> view(buf,mat.cols());
在这两种情况下,mat
对象必须在std::vector
或span
的整个生存期内保持活动状态。
为了完整性,您还可以在std::vector
:
std::vector<Vector3d> v(10);
Matrix<double,3,Dynamic>::Map(v.data().data(),3,10) = mat;
相关文章