Showing results for

- Intel Community
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
- Wrong solution computed by distributed PARDISO with 0-based indexed BSR matrix

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

Highlighted
##

asd__asdqwe

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-09-2016
08:07 AM

9 Views

Wrong solution computed by distributed PARDISO with 0-based indexed BSR matrix

Hello,

I'm attaching two .c files that should compute the solution of the linear system Ax = b, with

A = [1, -3, 5, 1;

-1, -2, 6, 10;

4, -4, 7, -5;

4, 2, 8, -1];

and

b = [1;1;0.75;1]

The solution given by MATLAB is the same as with the file cl_solver_unsym_distr_1.c, i.e.,

x =

-1.0000

0.1500

0.5500

-0.3000

With cl_solver_unsym_distr_0.c, I get a completly wrong solution: x = [-0.031169, -0.036084, 0.139033, 0.026384]'

Moreover, at the beginning of the execution, I get the following output:

1-based array indexing is turned ON

which is wrong, since I'm setting iparm[34] to 1.

Can you reproduce this issue? How can I fix it? I'm compiling the examples with the following command: make sointel64 mpi=mpich2 compiler=intel (inside the folder /opt/intel/mkl/examples/cluster_sparse_solverc of MKL 11.3.3)

9 Replies

Highlighted
##

asd__asdqwe

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-13-2016
02:53 AM

9 Views

Is this fixed in the upcoming update? Thank you in advance for an answer.

Highlighted
##

asd__asdqwe

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

09-19-2016
06:47 AM

9 Views

With the latest MKL from m_ccompxe_2017.0.036.dmg, I now have the following error with these scripts:

Fatal error in MPI_Reduce: Invalid count, error stack:

MPI_Reduce(1258): MPI_Reduce(sbuf=0x7fbae3700500, rbuf=0x7fbae3607500, count=-16380, MPI_DOUBLE, MPI_SUM, root=0, comm=0x84000004) failed

MPI_Reduce(1167): Negative count, value is -16380

Do you have a fix? Thanks.

Highlighted
##

Actually, this issue is targeted to be fixed in MKL 2017 update 1. We will keep you updated with the status of this issue.

Gennady_F_Intel

Moderator

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

09-19-2016
08:07 AM

9 Views

Highlighted
##

MKL 2017 update 1 has been recently released. Could you please check if the problem is still there and let us know! regards.

Gennady_F_Intel

Moderator

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

11-06-2016
09:04 PM

9 Views

Highlighted
##

asd__asdqwe

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

11-13-2016
02:05 PM

9 Views

The problem is still there, I still get a completely wrong solution with the 0-based approach:

The solution of the system is:

on zero process x [0] = -0.031169

on zero process x [1] = -0.036084

The solution of the system is:

on first process x [0] = 0.139033

on first process x [1] = 0.026384

With the 1-based approach, I still get the correct solution:

The solution of the system is:

on zero process x [0] = -1.000000

on zero process x [1] = 0.150000

The solution of the system is:

on first process x [0] = 0.550000

on first process x [1] = -0.300000

Can you fix this please? Thanks.

Highlighted
##

asd__asdqwe

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

11-16-2016
10:38 AM

9 Views

Can you reproduce this issue? Thank you.

Highlighted
##

Gennady_F_Intel

Moderator

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

11-17-2016
10:42 PM

9 Views

yes, we have to confirm the issue still exists with version 2017 u1.

Highlighted
##

Dear users, please accept our apologize for the wrong message regard to the previous message about the fix of this issue in MKL 2017 u1. Actually the fix of this issue will be available into MKL 2017 update 2. Regards, Gennady

Gennady_F_Intel

Moderator

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

11-23-2016
06:15 PM

9 Views

Highlighted
##

MKL 2017 update 2 is released a few days ago. Could you please check the problem with this update and keep us inform the status.

Gennady_F_Intel

Moderator

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

02-27-2017
06:40 AM

9 Views

For more complete information about compiler optimizations, see our Optimization Notice.