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

frankgty

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

fcode

[Fortran] 纯文本查看 复制代码

program fcode_cn
  Implicit None
  Real , allocatable    :: rMatrix(:,:) , value(:)
  Integer , allocatable :: ipodia(:) , indexJ(:)
  Integer :: N , M , i , j
  Open( 12 , File = "Ipodia.txt")
  Open( 13 , File = "Index J.txt")
  Open( 14 , File = "Value.txt")
  N = GetFileN( 12 )
  M = GetFileN( 13 )
  Allocate( rMatrix(N,N) , ipodia(0:N) , indexJ(M) , value(M) )
  rMatrix   = 0
  ipodia(0) = 0
  Read(12,*) ipodia(1:)
  Read(13,*) indexJ
  Read(14,*) value
  Do i = 1 , N
    rMatrix( indexJ( ipodia(i-1)+1 : ipodia(i) ) , i ) = value( ipodia(i-1)+1 : ipodia(i) )
    rMatrix( i , :i-1 ) = rMatrix( :i-1 , i )
  End Do
  Close(12)
  Close(13)
  Close(14)
  
contains

  Integer Function GetFileN( iFileUnit )
    Integer , Intent( IN ) :: iFileUnit
    character( Len = 1 ) :: cDummy
    integer :: ierr
    GetFileN = 0
    Rewind( iFileUnit )
    Do
      Read( iFileUnit , * , ioStat = ierr ) cDummy
      If( ierr /= 0 ) Exit
      GetFileN = GetFileN + 1
    End Do
    Rewind( iFileUnit )
  End Function GetFileN 
  
end program fcode_cn

pasuka

lz用ivf的话，直接查阅MKL帮助文档吧
譬如
https://software.intel.com/en-us ... -fortran-mkl-csrcoo

frankgty

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

非常感谢