mkl_sdnscsr 问题

frankgty · 发表于 2019-10-16 05:27:59

请教一个问题，用mkl中mkl_sdnscsr将dense矩阵转化为CSR格式，出现错误如下图，这里的parameter1指的是什么？

[Fortran] 纯文本查看 复制代码

01program converters
02    implicit none
03     
04    ! Definition arrays for sparse matrix formats
05    integer,parameter :: m=2400, n=2400, lda=2400, nzmax=11820
06    integer :: ia, ja
07    real :: adns(m,n) , acsr(nzmax)
08     
09    ! Declairation of local variables:
10    integer :: job(8), i, j, info, ibase1, ibase2, locat
11    info = 0
12    locat = 2
13    ibase1 = 1
14    ibase2 = 2
15    job(2) = ibase1
16    job(3) = ibase2
17    job(4) = locat
18    job(5) = nzmax
19     
20    ! Task1: densly matrix to crs
21    open(11,file="1.txt")
22    read(11,*) adns
23    job(1) = 0
24    job(6) = 1
25    call mkl_sdnscsr(job, m, n, adns, lda, acsr, ja, ia ,info)
26    print*, acsr, ia, ja
27    pause
28    end program converters

li913 · 发表于 2019-10-16 14:10:53

本帖最后由 li913 于 2019-10-16 14:12 编辑

查询mkl帮助，了解每个参数的信息。提示是指第一个参数的值有问题，也就是job有问题。
1、ja, ia应该是数组；
2、job(3)只能为0或1.

frankgty · 发表于 2019-10-16 18:51:45

感谢，是粗心造成的问题，将job(3) 改为1后，出现:0x00A422DC 处(位于 Formatconverters.exe 中)引发的异常: 0xC0000005: 写入位置 0x02502000 时发生访问冲突。
0x00000123 处有未经处理的异常(在 Formatconverters.exe 中): 0xC00001A5: 检测到无效的异常处理程序例程。 (参数: 0x00000003)。
在call mkl_sdnscsr处引发的问题。
程序窗口里是空白，没有信息，这种情况是哪里出了问题呢?

frankgty · 发表于 2019-10-16 19:26:02

li913 发表于 2019-10-16 14:10
查询mkl帮助，了解每个参数的信息。提示是指第一个参数的值有问题，也就是job有问题。
1、ja, ia应该是数组 ...

感谢，是粗心造成的问题，将job(3) 改为1后，出现:0x00A422DC 处(位于 Formatconverters.exe 中)引发的异常: 0xC0000005: 写入位置 0x02502000 时发生访问冲突。
0x00000123 处有未经处理的异常(在 Formatconverters.exe 中): 0xC00001A5: 检测到无效的异常处理程序例程。 (参数: 0x00000003)。
在call mkl_sdnscsr处引发的问题。
程序窗口里是空白，没有信息，这种情况是哪里出了问题呢?

li913 · 发表于 2019-10-17 14:30:43

ja ia 改了没？给出最新的代码，和输入文件。

frankgty · 发表于 2019-10-17 19:55:27

li913 发表于 2019-10-17 14:30
ja ia 改了没？给出最新的代码，和输入文件。

read文件时候出了问题，已解决得到正确结果，非常感谢您！

[数学库] mkl_sdnscsr 问题