use lapack95编译提示错误，请教大家。

静待花开 · 发表于 2022-12-1 22:08:16

我的一个子程序中用到了use lapack95，如下所示，主要功能就是矩阵求逆。

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

Subroutine InverseV2D
  use variables, Only : Np, InvV, V2D
  Use lapack95
  Implicit none
    Real(kind=4),    Allocatable :: ipiv(:)
  Allocate(ipiv(Np))
  ipiv = 0.0
  InvV = V2D
  Call getrf(InvV, ipiv)
  Call getri(InvV, ipiv)
End Subroutine InverseV2D

我按照下面做了设置，我用的是VS2022+Intel OneAPI Base toolkit+Intel OneAPI hpc toolkit,

但是编译时，总是提示下面错误。因为我在其它子程序中也用到了矩阵求逆，这些子程序中都出现了这种编译错误，有没有知道这是什么原因呢？

上面的代码是从网上参考了一个例子，并且我编译这个例子，也没有提示有错，下面是代码。我也把这个代码，以及.sln和.vfproj也上传上来。上面自己的项目属性设置，就是参考了这个.sln和.vfproj里的属性设置，但不知为何就是出错。

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

program MKL_getrfANDgetri
        use lapack95
        implicit none
        integer, parameter :: n = 3
        integer :: i, j, ipiv(n)
        real(kind=8) :: a(n,n), aa(n,n)

        call random_seed()
        call random_number(a)
        aa = a

        write(*,'(1x,a)') "a = "
        do i = 1, n
                write(*,'(*(f12.6,3x))') a(i,:)
        end do

        !// 使用lapack求解逆矩阵
        call getrf( a, ipiv )
        !write(6,*)ipiv
        call getri( a, ipiv )
        !write(6,*)ipiv

        write(*,'(1x,a)') 'inv(a) = '
        do i = 1, n
                write(*,'(*(f12.6,3x))') a(i,:)
        end do

        write(*,'(1x,a)') "checking..."
        aa = matmul(aa,a)  !// 原矩阵与其逆矩阵的结果为单位矩阵
        do i = 1, n
                write(*,'(*(f12.6,3x))') aa(i,:)
        end do

end program MKL_getrfANDgetri