关于动态数组的问题

fcode

稀疏矩阵最好使用专门的存储结构和求解器。

如果你真的要自己存储，可以试试这样：

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

program www_fcode_cn
  implicit none
  integer , allocatable :: a(:) , b(:)
  allocate( a(9) )
  a = [1,2,0,0,4,0,0,0,8,0]
  ! 如果你的编译器不支持，请加入以下语句
  !allocate( b( count(a/=0) ) )
  b = pack(a,a/=0)
  write(*,*) b  
end program www_fcode_cn

fcode

不可压缩 发表于 2015-12-7 20:03
我用的就是稀疏矩阵求解器，pardiso。
我要先把稀疏矩阵 A 计算(组装)出来，然后存储成pardiso需要的形式， ...

我并没有安装所有的编译器，所以也不知道哪些需要，哪些不需要。
但这个用法是 Fortran2008 语法，据我所知，并不是所有编译器都支持到。
你应该自己去尝试，如果无法运行，就加上后面的分配语句。

a 可以是任意维度。