稀疏Blas示例

Question

我需要使用稀疏算法来存储和乘以Fortran 90中的许多矩阵。我已经看到可以使用稀疏的Blas库来完成它。我想知道在哪里可以找到一些编译和使用这个库的例子(我对Fortran非常陌生)。我试着自己实现，我甚至不可能编写第一行。这是我的节目：

program main
use sparse_blas
IMPLICIT NONE
integer::istat
integer::i,j,NonZA,NonZB,nonZmul
double precision::A(4,4),B(4,4),mul(4,4)
!initialization of matrices

call DUSCR_BEGIN(4,4,A,istat)
end program main

我用的makefile：

objects = test.o
f90 = gfortran
fflags = -g -wall


test.x: $(objects)
    $(f90) -o test.x $(objects) -lblas

test.o: test.f90 
$(f90) -c test.f90 -lblas

Answer 1

要编译使用模块(例如sparse_blas)的文件，必须在makefile中指定正确的include path。

将makefile中的最后一行更改为

$(f90) -c test.f90 $(fflags) -Idirectory_which_contains_sparse_blas.mod