c++ - Eigen 3/C++ : MatrixXd multiply one row with another

Question

使用 Eigen3/C++ 库，给定一个 MatrixXd

    / x0   ...  y0  \
    | x1   ...  y1  |
M = | ...  ...  ... |
    |               |
    \ xN   ...  yN  /

快速实现如下所示修改版本的方法是什么？

     / x0 * y0   ...  y0  \
     | x1 * y1   ...  y1  |
M' = | ...  ...  ...      |
     |                    |
     \ xN * yN   ...  yN  /

即一列(带x-s的那一列)被自己替换与另一列(与 y-s 相乘)相乘。

Answer 1

你的意思是如何对第一个和最后一个列 vector 进行系数赋值和乘法？有很多方法可以做到这一点，但最简单/最快的可能是

Eigen::MatrixXd M2 = M;

M2.leftCols<1>().array() *= M2.rightCols<1>().array();

另一种方法可能是构造一个具有给定行数/列数的未初始化矩阵，然后像这样进行 block 分配

Eigen::MatrixXd M2{ M.rows(), M.cols() };

M2.rightCols( M.cols() - 1 ) = M.rightCols( M.cols() - 1 );
M2.leftCols<1>() = M.leftCols<1>().cwiseProduct( M.rightCols<1>() );

我不知道哪个更快(但您首选的分析器确实如此)。

以后的问题，这里是官方Eigen quick reference ;)