本征：复制（广播）由行

Question

我想复制矩阵M的每一行，而不发生任何拷贝（通过创建一个视图IE）：

0 1   0 1 
2 3 -> 0 1 
       2 3 
       2 3 

• M.rowwise().replicate(n)为shorcut M.replicate(1,n)这似乎有点无用。

• 以下代码片段会执行副本，并且如果M是表达式，则无法工作。

 
    Eigen::Index rowFactor = 2; 
    Eigen::MatrixXi M2 = Eigen::Map(M.data(), 1, M.size()).replicate(rowFactor, 1); 
    M2.resize(M.rows()*rowFactor, M.cols()) ; 

• 在某些情况下，我可能会重塑其它操作数使用中间视图Eigen::Map<Eigen::MatrixXi>(M.data(), 1, M.size()).replicate(rowFactor, 1)，但是这不是非常令人满意。

是否有适当的方法来实现此广播视图？

Answer 1

你想要的基本上是一个Kronecker product与矩阵的1。您可以使用的（不支持）KroneckerProduct module：

#include <iostream> 
#include <unsupported/Eigen/KroneckerProduct> 

int main() { 
    Eigen::Matrix2i M; M << 0, 1, 2, 3; 
    std::cout << Eigen::kroneckerProduct(M, Eigen::Vector2i::Ones()) << '\n'; 
} 

作为“不支持”意味着，该模块的API不保证是稳定的（尽管这个模块还没有推出以来改变，我认为）。