EIGEN地图类:将C阵列映射到VectorXD指针以有效的方式

Question

我需要将普通的双C阵列转换为vectorxd*。我找到了似乎可以完成这项工作的地图课程，但有我不理解的东西。

class A {
  private:
    VectorXd *jpos_;
  public:
    int init( double *v, int len )
    {
      // Here I would like to assign v to jpos_ in an efficient way
      // (something like a pointer assignment without allocation or vector iterations)
      // Note that length of v is len
    }
};

我尝试了以下内容，但它不起作用：

int A::init( double *v, int len )
{
  jpos_ = &VectorXd::Map( v, len );
}

或

int A::init( double *v, int len )
{
  jpos_ = &Map<VectorXd>( v, len );
}

怎么了？谢谢。

emanuele

Answer 1

如果要通过Eigen::Map映射数组，您的成员jpos_应为类型Eigen::Map：

class A
{
private:
    Eigen::Map<Eigen::VectorXd> jpos_;
public:    
    A(double* v, int len) : jpos_ {v, len} {}
};

如果您需要保留VectorXd，则需要将值复制到向量。看到VectorXd*有点奇怪。它类似于std::vector<double>*。我建议不要使用VectorXd*。

class A
{
private:
    Eigen::VectorXd jpos_;
public:    
    A(double* v, int len) : jpos_ {v}
    {
        for (int i=0; i<len; ++i)
            jpos_(i) = v[i];
    }
};

编辑：

如果您确实需要VectorXd*，则可以用new分配一个(并且不要忘记delete(。