Question

我需要在Fortran 90中使用稀疏算法进行存储并乘以许多矩阵。我已经看到它可以使用Sparse 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

要编译使用MODULE的文件（例如sparse_blas），您必须在makefile中指定正确的include path。

将makefile中的最后一行更改为

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

稀疏Blas的例子

1 个答案: