稀疏矩陣的特徵 - 對角線更新

Question

在特徵中編輯稀疏矩陣對角線的最快方法是什麼？我已經使用三元組填充它，但我偶爾需要更改所有對角線值（它們已被設置爲非零）。

Answer 1

如果接受，你可以簡單地使用他們的指標像

sp_mat(i, i) 

由於每個訪問包括二進制搜索，其性能可能無法滿足。在這種情況下，您仍然可以通過一些限制獲得更好的性能。

1. 您不會更改矩陣的基礎佈局，例如，添加/刪除nnzs，壓縮;
2. 您將重複訪問對角線。

有了這些限制，您可以存儲對角線組件的指針，以便您可以稍後在O（1）時間訪問數據。你可以得到的指針與

&sp_mat(i, i) 

Answer 2

要完成kangshiyin的回答，您還可以編輯通過密集的矢量表達式對角線，例如到1添加到每個對角線元素：

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

警告：這隻有在對角元素已經存在的情況下才有效，原始問題就是這種情況。