我刚刚下载了MKL,打算试试它的稀疏BLAS,首先尝试的是用MKL_SDNSCSR将全矩阵转换为稀疏矩阵。
这是我的代码。
implicit none
real::A(2,5)
integer::job(8)
INTEGER:: info,n=2,m=5
INTEGER ja(10), ia(6)
REAL acsr(10)
A(1,:)=(/1, 0 ,1,0,2/)
A(2,:)=(/2,0,0,0,1/)
job=0
job(1)=0
job(2)=1
job(3)=1
job(4)=2
job(5)=10
job(6)=1
print*,A
print*,job
call mkl_sdnscsr(job,n,m,A,m,acsr,ja,ia,info)
end
它给出了以下错误。
Intel MKL ERROR: Parameter 1 was incorrect on entry to MKL_SDNSCSR.
我做错了什么?
最好的
我觉得这个例子没有问题。我在win64上尝试了MKL v.2020.这里是我得到的输出。
1.000000 2.000000 0.0000000E+00 0.0000000E+00 1.000000
0.0000000E+00 0.0000000E+00 0.0000000E+00 2.000000 1.000000
0 1 1 2 10 1
0 0