Question

想法是将值从std :: vector分配给Eigen :: MatrixXd中的块。

我的问题是，有什么好方法吗？

 1 1 1        1   1   1
 1 1 1    to  1 102 123
 1 1 1        1   1   1

我尝试将std :: vector转换为Eigen :: Map，但没有成功。我有一个有效的代码[以下代码段]，但它吸引力不大。也许有一种更简单的方法？

    void do_work(Eigen::MatrixXd &m, const std::vector<double> &v,
                 const Index i, const Index j, 
                 const Index p, const Index q) {
      auto stop = j + q;
      for (Index start = j, idx = 0; start < stop; ++start, ++idx)
        m.block(i, start, p, 1) << v[idx];
    }

    Eigen::MatrixXd m(3, 3);
    m << 1, 1, 1, 1, 1, 1, 1, 1, 1;
    std::vector<double> v = {102, 123};
    Index change_row = 0;
    Index change_column_from = 1, change_column_to = v.size();
    do_work(m, v, change_row, change_column_from, 1, change_column_to);

预期结果将是以高效（甚至更干净）的方式执行操作。

Answer 1

例如将std::vector<double>转换为Eigen::Map，例如

Eigen::Map<Eigen::VectorXd> map(v.data(), v.size());

用于可读写访问，或者，如果您对矢量具有只读访问权限，则执行以下操作：

Eigen::Map<const Eigen::VectorXd> map(v.data(), v.size());

这可以分配给一个块表达式，假设块的大小与Map的大小匹配：

m.block(i,j, map.size(), 1) = map.transpose();
                               // transpose is actually optional here

使用std :: vector的Eigen :: MatrixXd.block分配

1 个答案: