主题:这段语句怎么写成稀疏矩阵?
大家好,我对稀疏矩阵很模糊,在网上和论坛上搜了下一直没找到该怎么去构建稀疏矩阵。我这里有个图形,是六角格子,沿着x方向具有平移对称性,也就是说x方向是一样的,y方向我取8个点为一个原包,也就是y方向8个点为一个最小单元,我现在把这8个点的哈密顿矩阵写出来:
0 t 0 0 0 0 0 t
t 0 t 0 0 0 0 0
0 t 0 t 0 0 0 0
0 0 t 0 t 0 0 0
0 0 0 t 0 t 0 0
0 0 0 0 t 0 t 0
0 0 0 0 0 t 0 t
t 0 0 0 0 0 t 0
t的近邻两点之间相互作用,从哈密顿形式上看,是稀疏矩阵,我有一个问题:怎么用fortran把它写成稀疏矩阵?比如我把上面的哈密顿写出来之后,怎么把它转成稀疏矩阵?
这个小程序就是上面这个哈密顿矩阵的写法。
program xishu
implicit none
real(8) :: t=1.0
integer :: i1,ny
real(8),allocatable :: h(:,:)
ny=8
allocate(h(ny,ny))
h=0.0
do i1=1,ny-1
h(i1,i1+1)=t
h(i1+1,i1)=t
enddo
h(1,ny)=t
h(ny,1)=t
end
我想把它转成稀疏矩阵,该怎么办呢?sparse命令到底是怎么用的呢?纠结啊!!