Eigen: Replizieren (broadcast) durch Reihen

Question

Ich mag würde jede Zeile einer Matrix M ohne Kopie replizieren auftritt (dh durch Erstellen einer Ansicht):

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

• M.rowwise().replicate(n) ist ein shorcut für M.replicate(1,n) das scheint irgendwie nutzlos.

• Das folgende Snippet erstellt eine Kopie und funktioniert nicht, wenn M ein Ausdruck ist.

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

• In einigen Situation kann ich die Zwischenansicht Eigen::Map<Eigen::MatrixXi>(M.data(), 1, M.size()).replicate(rowFactor, 1) durch Umformen die anderen Operanden verwenden, aber das ist nicht sehr befriedigend.

Gibt es einen richtigen Weg, um diese Broadcast-Ansicht zu erreichen?

Answer 1

Was Sie wollen, ist im Wesentlichen ein Kronecker product mit einer Matrix von Einsen. Sie können die (nicht unterstützt) KroneckerProduct module dafür verwenden:

#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'; 
} 

Sein ‚nicht unterstützt‘ bedeutet, dass die API des Moduls nicht stabil sein wird garantiert (auch wenn dieses Modul seit seiner Einführung hat sich nicht geändert, glaube ich).