スパース行列の固有対角更新

Question

Eigenの疎行列の対角を編集する最も速い方法は何ですか？私はすでにトリプレットを使ってそれを埋めましたが、時にはすべての対角値を変更する必要があります（既にゼロ以外に設定されています）。

Answer 1

許容できる場合は、各アクセスはバイナリ検索を必要とするとして、あなたは、単に

sp_mat(i, i) 

のように、それぞれのインデックスを使用することができ、パフォーマンスは満足できない場合があります。この場合でも、いくつかの制限事項でより良いパフォーマンスを達成することができます。

1. マトリックスのレイアウトを変更することはありません。追加/削除nnzs、圧縮;
2. 対角線に繰り返しアクセスします。

これらの制限を使用すると、後でO（1）時間後にデータにアクセスできるように対角コンポーネントのポインタを格納できます。あなたはkangshiyinの答えを完了するために

&sp_mat(i, i) 

Answer 2

でポインタを得ることができる、あなたはまた、各対角要素に1を追加するたとえば、密ベクトル表現を通じて対角線を編集することができます：警告

sp_mat.diagonal().array() += 1; 

：対角要素がすでに存在する場合にのみ動作します。これは元の質問の場合です。