Eigen MatrixXdはC++でプッシュバック

Question

EigenはC++のよく知られた行列ライブラリです。私は、単純に行列の最後にアイテムをプッシュする組み込み関数を見つけるのが面倒です。

Eigen::MatrixXd matrix(10, 3); 
long int count = 0; 
long int topCount = 10; 
for (int i = 0; i < listLength; ++i) { 
    matrix(count, 0) = list.x; 
    matrix(count, 1) = list.y; 
    matrix(count, 2) = list.z; 
    count++; 
    if (count == topCount) { 
     topCount *= 2; 
     matrix.conservativeResize(topCount, 3); 
    } 
} 
matrix.conservativeResize(count, 3); 

そして、これは（構文の一部が出てかもしれ）動作します。現在、私はそれはこのように行うことができることを知っています。しかし、それはかなり単純なことのために巻き込まれました。すでに組み込まれている機能はありますか？

Answer 1

固有行列にはそのような関数はありません。この理由は、そのような関数が非常に遅いか、または過剰なメモリを使用するためです。

push_backは非常に高価ではないため、実行したときにスペースが足りなくなったときに、何らかの要因でマトリックスの容量を増やす必要があります。しかし、行列を扱うときには、メモリ使用量が問題になることが多いので、行列の容量を必要以上に大きくすることは問題になる可能性があります。 代わりに、操作がO(n*m)になるたびに、rows()またはcols()のサイズを増やしてください。行列全体を埋めるためにこれを行うことは、適度なサイズの行列でさえかなり遅くなるO(n*n*m*m)です。

さらに、線形代数行列およびベクトルサイズは、ほとんど常に一定であり、事前にわかっています。行列のサイズを変更すると、行列の以前の値は気にしません。これは、std::vectorのresizeとは異なり、Eigenのresize関数が古い値を保持しない理由です。

ファイルから読み込むときに、事前に行列のサイズがわからない場合があります。この場合、std::vectorなどの標準コンテナにデータをロードしてからpush_backを使ってデータをロードし、それを既にサイズの大きい行列にコピーするか、またはメモリをきちんと読み込んでサイズを取得した後にもう一度コピーしますその価値。

Answer 2

そのような機能はありませんが、しかし、あなたはこのように自分で何かを構築することができます。

using Eigen::MatrixXd; 
using Eigen::Vector3d; 

template <typename DynamicEigenMatrix> 
void push_back(DynamicEigenMatrix& m, Vector3d&& values, std::size_t row) 
{ 
    if(row >= m.rows()) { 
     m.conservativeResize(row + 1, Eigen::NoChange); 
    } 
    m.row(row) = values; 
} 


int main() 
{ 
    MatrixXd matrix(10, 3); 
    for (std::size_t i = 0; i < 10; ++i) { 
     push_back(matrix, Vector3d(1,2,3), i); 
    } 
    std::cout << matrix << "\n"; 
    return 0; 
} 

これはかかわらず、あまりにも多くのサイズ変更を実行する必要がある場合、恐ろしく遅いことになるだろう。