Intel Visual Fortran MKL DSS 求解大型稀疏矩阵方程出错

duty

IVF入门，现需要调用DSS求解大型稀疏矩阵方程，修改了mkl自带的example dss_sym_f.f，仅修改了部分输入数据，可以正常编译，但debug时执行到
error = dss_reorder( handle, MKL_DSS_DEFAULTS, idum)
出现中断，显示如下。跪求大神指点。小弟用的是VS 2010 和 Intel Parallel Studio XE 2013。代码和输入数据文件见附件。

duty

谢谢斑竹大神，我在Intel的论坛发帖问了，下面是回复
HI,

I've put your data into PARDISO interface and obtain following error:

*** Error in PARDISO  (incorrect input matrix  ) error_num= 24
*** Input check: j=1410, ja(j)=72, ja(j+1)=67 are incompatible

The reason of this failure: in 72th row elements store in incorrect format - column of 1411th element is less than column of element with 1410 number in spite of they lie in one row!

Thanks,

Alex发现输入数据确实是乱序的~正在改

duty

aliouying 发表于 2014-3-2 08:39
排序后果然可以求解，结果一致

想问一下superLU和MUMPS有源码吗？