sparse right hand side reordering problem

Hailong_X_ — Mon, 27 May 2013 23:05:13 GMT

Dear All

I am trying to use the sparse right hand feature of pardiso ( iparm(31)=1 ).

I tested with an identity matrix A.

When I set perm = 1 for last several entries except the last one, I got the following error,

*** Error in PARDISO ( reordering_phase) error_num= -180

*** error PARDISO: reordering, symbolic factorization

perm before reordering
0 0 0 0 0 0 0 0 1 1 1 1 1 0
perm after reordering
8 7 4 3 6 2 10 1 5 9 10 11 12 13

Notice: 10 appears twice.

But when I set, the last entry of perm is 1, i.e. perm[last] = 1, then there is no problem.

perm before reordering
0 0 0 0 0 0 0 0 1 1 1 1 1 1
perm after reordering
8 7 6 5 4 3 2 1 9 10 11 12 13 14

The attachment is my C++ test code,

Could you please give me some suggestions?

Hailong

Hailong,

Chao_Y_Intel — Tue, 28 May 2013 02:32:18 GMT

Hailong,

Thanks for your report. We will have some check for the problem.

Regerds,
Chao

Chao_Y_Intel — Wed, 26 Feb 2014 07:45:06 GMT

Hello,

An update on this problem: it was fixed since MKL 11.1.1. Here is list of the MKL bug fix: http://software.intel.com/en-us/articles/intel-mkl-111-bug-fixes

Thanks,
Chao