Fortran Coder

标题: use lapack95编译提示错误，请教大家。 [打印本页]

作者: 静待花开 时间: 2022-12-1 22:08
标题: use lapack95编译提示错误，请教大家。
我的一个子程序中用到了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

网上找的例子.rar (1.85 KB, 下载次数: 0)

作者: fcode 时间: 2022-12-2 08:20
InvV 如何定义的试试把
Real(kind=4), Allocatable :: ipiv(:)
改为
Integer, Allocatable :: ipiv(:)

作者: 静待花开 时间: 2022-12-2 11:38

fcode 发表于 2022-12-2 08:20
InvV 如何定义的试试把
Real(kind=4), Allocatable :: ipiv(:)
改为

按您的方法已经解决。佩服，真的佩服，一眼就看出来问题，我找了半天，……。

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)