本帖最后由 frankgty 于 2019-9-28 00:19 编辑
[Fortran] 纯文本查看 复制代码 include 'lapack.f90'
Program kontrolle
use lapack95
Implicit None
Integer,parameter :: n=2400 , nrhs=1
Integer :: i
Integer,allocatable :: ipiv(:)
Real,allocatable :: A(:,:) , B(:,:)
Allocate( A(n,n) , B(n,nrhs) , ipiv(n) )
Open(12,File="Vm.txt")
Open(13,File="RHS.txt")
Read(12,*) A
Read(13,*) B
Close(12)
Close(13)
Call sgetrf(n,n,A,n,ipiv,i)
Call sgetrs('N',n,n,A,n,ipiv,B,n,i)
Print*,B
Pause
End Program kontrolle
您好,十分感谢您的指导,我参照您的例子改写了代码,出现了错误(0x00007FF93CC2B811 (mkl_core.dll)处(位于 Kontrolle.exe 中)引发的异常: 0xC0000005: 读取位置 0x000001D26396D000 时发生访问冲突。) |