我需要解決電阻網絡研究中出現的一些大的(N〜1e6)拉普拉斯矩陣。網絡分析的其餘部分正在使用boost圖處理,如果可能的話,我想留在C++中。我知道有很多很多C++矩陣庫,但似乎沒有人在速度或可用性方面明顯領先。另外,關於這個主題的許多問題,這裏和其他地方似乎都迅速轉化爲效用有限的洗衣清單。爲了幫助自己和他人,我會盡量保持問題的簡明性和可回答性:什麼是大拉普拉斯矩陣的快速簡單求解器?
什麼是可以有效處理以下要求的最佳庫?
- 矩陣類型:對稱角佔優/拉普拉斯
- 尺寸:非常大(N〜1E6),無動態調整所需
- 稀疏:極限(每行的最大5非零項/列)
- 需要的操作:求解A * x = b和mat/vec中的x
- 語言:C++(C ok)
- 優先級:代碼的速度和簡單性。我真的寧願避免爲這個問題學習一個全新的框架,或者不得不手動編寫太多的幫助代碼。
超愛用最少的工作示例答案...
我相信KLU非常適合您的需求:http://www.cise.ufl.edu/research/sparse/klu/ –