Fortran Coder

查看: 426|回复: 2

[数学库] use lapack95编译提示错误,请教大家。

[复制链接]

36

帖子

12

主题

0

精华

熟手

F 币
142 元
贡献
80 点
发表于 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,
1.PNG 2.PNG

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

上面的代码是从网上参考了一个例子,并且我编译这个例子,也没有提示有错,下面是代码。我也把这个代码, 以及.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)

1820

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
760 元
贡献
305 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2022-12-2 08:20:24 | 显示全部楼层
InvV 如何定义的试试把
Real(kind=4),    Allocatable :: ipiv(:)
改为
Integer,    Allocatable :: ipiv(:)

36

帖子

12

主题

0

精华

熟手

F 币
142 元
贡献
80 点
 楼主| 发表于 2022-12-2 11:38:01 | 显示全部楼层
fcode 发表于 2022-12-2 08:20
InvV 如何定义的试试把
Real(kind=4),    Allocatable :: ipiv(:)
改为

按您的方法已经解决。佩服,真的佩服,一眼就看出来问题,我找了半天,……。
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2023-2-4 16:14

Powered by Tencent X3.4

© 2013-2023 Tencent

快速回复 返回顶部 返回列表