本徵：複製（廣播）由行

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不保證是穩定的（儘管這個模塊還沒有推出以來改變，我認爲）。