use lapack95编译提示错误,请教大家。
我的一个子程序中用到了use lapack95,如下所示,主要功能就是矩阵求逆。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里的属性设置,但不知为何就是出错。
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
InvV 如何定义的试试把
Real(kind=4), Allocatable :: ipiv(:)
改为
Integer, Allocatable :: ipiv(:)
fcode 发表于 2022-12-2 08:20
InvV 如何定义的试试把
Real(kind=4), Allocatable :: ipiv(:)
改为
按您的方法已经解决。佩服,真的佩服,一眼就看出来问题,我找了半天,……。
页:
[1]