Fortran Coder»论坛 › 研讨区 › 编程工具交流 › Intel Visual Fortran MKL DSS 求解大型稀疏矩阵方程出 ...

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

查看数: 20853 | 评论数: 11 | 收藏 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

duty

发布时间: 2014-3-1 15:12

正文摘要:

IVF入门，现需要调用DSS求解大型稀疏矩阵方程，修改了mkl自带的example dss_sym_f.f，仅修改了部分输入数据，可以正常编译，但debug时执行到 error = dss_reorder( handle, MKL_DSS_DEFAULTS, idum) 出现中断，显示 ...

aliouying 发表于 2014-3-2 18:14:04

duty 发表于 2014-3-2 15:36
想问一下superLU和MUMPS有源码吗？

网上开源可以下载

duty 发表于 2014-3-2 15:36:11

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

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

aliouying 发表于 2014-3-2 08:39:08

排序后果然可以求解，结果一致

aliouying 发表于 2014-3-2 08:29:41

本帖最后由 aliouying 于 2014-3-2 08:30 编辑

duty 发表于 2014-3-2 02:16
谢谢斑竹大神，我在Intel的论坛发帖问了，下面是回复
HI,

哦，怪不得，superLU和MUMPS对数据的顺序没要求

只需要排序下即可，我到算法板块发下排序算法

aliouying 发表于 2014-3-2 08:29:41

duty 发表于 2014-3-2 02:16
谢谢斑竹大神，我在Intel的论坛发帖问了，下面是回复
HI,

哦，怪不得，superLU和MUMPS对数据的顺序没要求

duty 发表于 2014-3-2 02:16:55

谢谢斑竹大神，我在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发现输入数据确实是乱序的~正在改

aliouying 发表于 2014-3-1 16:44:11

retval = HeapFree(_crtheap, 0, pBlock); 定在这一行

aliouying 发表于 2014-3-1 16:43:42

chuxf 发表于 2014-3-1 16:40
我想这个问题恐怕不是你的错，你的算例和代码都没有问题。

我在 intel 的论坛发现了去年的一个帖子提及 ds ...

DSS有错误提示，UMFPack错误提示是-8，pardiso直接跳出C的子程序free.c

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

void __cdecl _free_base (void * pBlock)
{

        int retval = 0;


        if (pBlock == NULL)
            return;

        RTCCALLBACK(_RTC_Free_hook, (pBlock, 0));

        retval = HeapFree(_crtheap, 0, pBlock);
        if (retval == 0)
        {
            errno = _get_errno_from_oserr(GetLastError());
        }
}

楚香饭 发表于 2014-3-1 16:40:18

我想这个问题恐怕不是你的错，你的算例和代码都没有问题。

我在 intel 的论坛发现了去年的一个帖子提及 dss_reorder 出错，后来发布了新版本解决这一问题。但是是 C 接口的。

我不确定是不是 Fotran 的接口没有改。

MKL 是一个商业化的函数库，就算是矩阵问题，也会有良性的错误提示或返回值。而且 reorder 只是一个简单操作，对数据应该不会有多大依赖性。