对 LAPACK 和 BLAS 子例程的未定义引用

我试图了解 Fortran 中的 BLAS 和 LAPACK 如何工作等等，所以我编写了一个生成矩阵并将其反转的代码。

这是代码

program test
Implicit none
external ZGETRF
external ZGETRI
integer ::M  
complex*16,allocatable,dimension(:,:)::A
complex*16,allocatable,dimension(:)::WORK
integer,allocatable,dimension(:)::IPIV
integer i,j,info,error

Print*, 'Enter size of the matrix'
Read*, M 

Print*, 'Enter file of the matrix'
READ(*,*), A

OPEN(UNIT=10,FILE = '(/A/)' ,STATUS='OLD',ACTION='READ')

allocate(A(M,M),WORK(M),IPIV(M),stat=error)
    if (error.ne.0)then
        print *,"error:not enough memory"
        stop
    end if

 !definition of the test matrix A
  do i=1,M
  do j=1,M
    if(j.eq.i)then
       A(i,j)=(1,0)
    else 
       A(i,j)=0
    end if
 end do
end do  

 call ZGETRF(M,M,A,M,IPIV,info)
 if(info .eq. 0) then
   write(*,*)"succeded"
 else
  write(*,*)"failed"
 end if

 call ZGETRI(M,A,M,IPIV,WORK,M,info)
  if(info .eq. 0) then
   write(*,*)"succeded"
 else
  write(*,*)"failed"
 end if
 deallocate(A,IPIV,WORK,stat=error)
 if (error.ne.0)then
  print *,"error:fail to release"
  stop
 end if  
 close (10)    
end program test

矩阵 A 位于我正在调用的文件中，并且我还表示矩阵的大小 (M )。当我用 gfortran 复制它们时，我收到这些消息

/tmp/ccVkb1zY.o：功能中MAIN__': test.f03:(.text+0x751): undefined reference tozgetrf_' test.f03:(.text+0x85d): 对“zgetri_”的未定义引用 collect2：错误：ld 返回 1 退出状态

我已经安装了 BLAS 和 LAPACK，所以我不知道我是否以正确的方式调用库。

有什么建议吗？

看起来您可能没有链接到这些库。尝试：

gfortran -o test test.f03 -llapack -lblas

这会导致链接器（将所有程序部分连接在一起的程序；在 UNIX 上通常称为“ld”）在程序中包含 LAPACK 调用的库代码（或其动态链接）。

如果上面一行的结果是“cannot find -llapack”或类似的结果，则有两个常见问题：

1. 库可以是“共享”（名称以“.so”结尾）或“静态”（名称以“.a”结尾）；链接器将查找共享链接，因此如果您只有静态链接，则应在库链接之前添加“-static”：

   gfortran -o test test.f03 -static -llapack -lblas
   

   这也将使它寻找 BLAS 的静态版本；如果需要共享版本，请在“-lblas”前面添加“-shared”：

   gfortran -o test test.f03 -static -llapack -shared -lblas

   你可能会发现这一页有帮助。

2. 链接器没有在正确的目录中查找库。您需要找到实际的库（称为“liblapack.so”或“liblapack.a”），并确保它所在的目录包含在链接器查找的目录中，例如让它也可以在“/mylibs/maths”中查找：

   gfortran -o test test.f03 -L/mylibs/maths -llapack -lblas