如何将压缩储存的稀疏矩阵转化为满秩矩阵

frankgty · 发表于 2018-6-25 22:18:08

本帖最后由 frankgty 于 2018-6-25 23:46 编辑

现有一个稀疏矩阵以ipodia法压缩储存，数据文件中包含元素值，下标j和ipodia共三列数组。储存方式如图ipodia法所示，ipodia指向每行最后一个非零元素，此矩阵为对称矩阵，只储存了下三角。跪求大神指导如何用fortran实现将此压缩储存矩阵转化为全矩阵。求一个算法或者思路。

frankgty · 发表于 2018-6-28 04:39:59

fcode 发表于 2018-6-26 08:28
[mw_shl_code=fortran,true]program fcode_cn
Implicit None
Real , allocatable :: rMatrix(:,:) , ...

非常感谢

[稀疏矩阵] 如何将压缩储存的稀疏矩阵转化为满秩矩阵