我有一些 vector vec,我想通过沿 vector 维度复制值来获得新的“表达式”vec2
Eigen::VectorXf vec(5);
vec << 1, 2, 3, 4, 5;
const auto vec2 = vec.someAwesomeEigenMagic<3>();
//vec2 should contains (1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5)^T
//Not (1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5)^T
当然,我可以手动创建这样的 vector ,或者通过 Eigen::Map 使用复制+向量化:
MatrixXf tmp = vec.replicate(1, 3).transpose();
const Map<VectorXf> vec2(tmp.data(), vec.rows() * 3, 1);
但是我希望 vec2 成为某种没有分配的“特征模板表达式”(vec 可能很大,我会经常调用这个例程)并立即计算值。 (vec 包含每个顶点的权重,我想将它用于加权最小二乘)
我考虑过用 1 vector 进行克罗内克乘积技巧,但我不确定它是否针对 1 乘积进行了优化。另外我更喜欢避免不受支持的模块
PS对不起我的英语
最佳答案
使用 devel 分支,您可以使用 LinSpaced 生成索引序列,然后对输入 vector 进行索引:
#include <iostream>
#include <Eigen/Dense>
using namespace Eigen;
using namespace std;
int main()
{
VectorXf vec(5);
vec << 1, 2, 3, 4, 5;
auto vecrep = vec(ArrayXi::LinSpaced(5*3,0,4));
cout << vecrep.transpose() << endl;
}
然后,您可以将关键行包装在返回 auto
的自由函数中,在 C++14 中:
template<typename XprType>
auto magic_rep(const XprType &xpr, Index K) {
return xpr(Eigen::ArrayXi::LinSpaced(xpr.size()*K,0,xpr.size()-1));
}
在主要部分:
cout << magic_rep(vec,3).transpose() << endl;
关于c++ - Eigen:沿一维复制项目,没有无用的分配,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/46110917/