利用MKL求解矩阵方程出错
利用一个简单的方程验证自己用的对不对的，但是发现出错了program mainuse mkl95_precision
use mkl95_lapack
implicit none
integer i,info
integer ipiv(3)
real*8 ::a(3,3),b(3,1),X(3,1)
a=0
b=1
forall( i=1:3)
a(i,i)=1
end forall
print *,'Matrix A is:'
do i=1,3
print *,a(i,:)
end do
print *,'Matrix B is:'
print *,b
call getrf(a,ipiv,info)
call getrs(a,ipiv,b,'N',info)
print *,'Solution X is:'
print *,b
end program
运行的结果，不对，不知道哪里出错了，请教一下，谢谢~
奇怪，我这里是正常的。得到的结果也是正确的
fcode 发表于 2017-12-6 17:17
奇怪，我这里是正常的。得到的结果也是正确的
难道是我的库的链接设置错了吗还是？ 这个是include的位置：C:\Program Files (x86)\Intel\Composer XE 2013 SP1\mkl\include
linker的library位置：C:\Program Files (x86)\Intel\Composer XE 2013 SP1\mkl\lib\intel64
library那里改成了parallel
页:
[1]