行でEigen :: MatrixXdへのスレッドセーフ書き込み

Question

私の質問は非常に簡単で、うまくいけばうまくいけばいいです：構築済みのEigen::MatrixXd行列があるとき、複数のスレッドを使って同じ行列時間がない場合（つまり、行が並行して書き込まれていないことが保証されている場合）、または各スレッドで一時行オブジェクトを作成して、それらを行列にコピーする必要があります。

Answer 1

Eigen::MatrixXdはカラムのメジャーストレージであるため、異なるスレッドから同じアドレスに書き込まないという点ではスレッドセーフである可能性がありますが、キャッシュ上に混乱が生じる可能性があります（基本的にはfalse sharingです）。一時的な行のメジャー・マトリックスを作成し、それをコラム・メジャー・マトリックスにコピーする方が速いかもしれません。

また、既存のマトリックスの列を行として扱い（寸法が変更/一致していることを確認してから）、m.transposeInPlace()を実行することができます。マトリックスの形状と配置によっては、m = m.transpose().eval()より効率的です。行列が十分に大きいとIDがゼロである場合も

はスレッドのIDを使用することも可能であるベースと例えばOMPとしたり、独自に異なるIDを追跡することなく、std::thread例えば、類似していない（連続しました）。 これはまた、行数がキャッシュラインサイズの倍数になるように行列をパディングする必要があり、各列は整列したメモリブロックで開始されます。 キャッシュラインが64バイトであると仮定します。整数倍のブロックを扱う場合、各スレッドは「独自の」キャッシュ行にしか触れないので、誤った共有を避けることができます。あなたがこれを行うことができるなら、余分な一時的なものやコピー/スワップはないはずです。