特征库，将矩阵的列乘以方阵

【问题标题】：Eigen library, multiply columns of matrix by square matrix特征库，将矩阵的列乘以方阵
【发布时间】：2012-11-21 11:30:51
【问题描述】：

我想通过算子 N 来转换矩阵 M 的每一列。Eigen 允许用预乘法来表达这一点：

M.colwise() *= N;

但乘法 M_j * N 在数学上是未定义的。

有没有办法避免写循环？

【问题讨论】：

【解决方案1】：

如果你想将M的每一列从左边乘以N，只需执行一个普通的matrix-matrix-multiplication：

M = N * M;

这会将N*M 评估为一个临时对象，然后将其移至M。如果您经常这样做并希望为此重新使用分配的内存，请在之前的某处声明一个临时矩阵 M_temp 并写入

M_temp.noalias() = N * M;
M.swap(M_temp);  // M_temp has the old memory of M; Swapping is O(1)

如果怕内存消耗太多，可以写类似

for(long i=0; i<M.cols()-3; i+=4)
    M.middleCols<4>(i) = N * M.middleCols<4>(i);
M.rightCols(M.cols()%4) = N * M.rightCols(M.cols()%4);

【讨论】：