2016-11-19 32 views
2

我有尺寸3329×3329有很多零個的字段的矩陣m,我想計算m^9乘稀疏矩陣。 使用matrix包(Data.Matrix易於使用)嘗試此操作後,我發現稀疏矩陣可以更好地表示內存使用情況和計算速度。所以我想弄清楚如何使用hmatrix包。我已經成功地創建一個稀疏矩陣:如何使用HMATRIX

module Example where 

import Numeric.LinearAlgebra as LA 

assocExample :: AssocMatrix 
assocExample = [((0,0), 1),((3329,5),1)] 

sparseExample :: GMatrix 
sparseExample = LA.mkSparse assocExample 

我在這一點上的問題似乎是,我有一個GMatrix,但乘法運算(<>)我需要一個Matrix t來代替。 通過查找凹槽文件hmatrix我沒有設法弄清楚如何在這裏獲得Matrix t

我也有一個在introduction to hmatrix快速凝視,但沒有提到它的術語sparse

我的預感是,這應該很容易做到,但我缺少一些簡單的東西。

回答

1

稀疏矩陣對我來說是相當年輕的hmatrix。通過文檔查看,似乎沒有稀疏矩陣的產品。你必須自己實現它。

編輯:如果你這樣做了,在這裏評論:https://github.com/albertoruiz/hmatrix/issues/162(也證實了以上我的發言)

+0

嘿感謝 - 我會研究這一點,並很可能接受你的答案。這種可能性甚至不會發生在我身上 –