Question

我有一些矢量vec，我想获得一个新的＆＃34;表达＆＃34; 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

当然我可以手动创建这样的矢量，或者使用Eigen :: Map：

的复制+矢量化

MatrixXf tmp = vec.replicate(1, 3).transpose();
const Map<VectorXf> vec2(tmp.data(), vec.rows() * 3, 1);

但我希望vec2成为某种特殊的模板表达式＆＃34;没有说明（vec可能很大，我会经常调用这个例程）并立即计算值。（vec包含每个顶点权重，我想将它用于加权最小二乘）

我想到kronecker产品技巧与矢量的一个，但我不确定它是否针对产品优化。此外，我更愿意避免使用不受支持的模块

PS抱歉我的英文

Answer 1

使用devel分支，您可以使用LinSpaced生成索引序列，然后索引输入向量：

#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;

Eigen：沿着一个维度复制项目而没有无用的分配

1 个答案: