Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Intel Community
- Software
- Software Development SDKs and Libraries
- Intel® oneAPI 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

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

72 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)

Link Copied

9 Replies

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

72 Views

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

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

72 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.

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

72 Views

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

72 Views

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

72 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.

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

72 Views

Can you reproduce this issue? Thank you.

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

72 Views

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

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

72 Views

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

72 Views

Topic Options

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

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