sparse right hand side reordering problem

[Hailong_X_]

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

[Chao_Y_Intel]

Hailong, 

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

Regerds,
Chao 

[Chao_Y_Intel]

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