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
- cluster_sparse_solver and distributed format

- 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

09-01-2014
12:47 AM

62 Views

cluster_sparse_solver and distributed format

Hello,

I'm trying to use cluster_sparse_solver and solve a system in-place (iparm(6) = 1), with a distributed format (iparm(40) = 2). I adapted the example cl_solver_unsym_distr_c.c as you can see attached, and at runtime, on two MPI processes, I get the following output:

The solution out-of-place of the system is:

on zero process x [0] = 0.263109 rhs [0] = 1.000000

on zero process x [1] = 0.305243 rhs [1] = 1.000000

on zero process x [2] = -0.347378 rhs [2] = 0.250000

The solution out-of-place of the system is:

on first process x [0] = -0.347378 rhs [0] = 0.750000

on first process x [1] = 0.205993 rhs [1] = 1.000000

on first process x [2] = 0.288390 rhs [2] = 1.000000

Solving system in-place...

The solution in-place of the system is:

on zero process x [0] = 1.000000

on zero process x [1] = 1.000000

on zero process x [2] = 0.250000

The solution in-place of the system is:

on first process x [0] = 0.750000

on first process x [1] = 1.000000

on first process x [2] = 1.000000

Can you reproduce this behavior ? The solution in-place is obviously wrong. Do you see how to fix that ? Thank you in advance.

Link Copied

6 Replies

Chao_Y_Intel

Employee

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

09-02-2014
01:49 AM

62 Views

Hello,

Thanks for your report. Our engineer owner could verify this is problem, and they are planning to fix it in our near future release.

Regards,

Chao

Gennady_F_Intel

Moderator

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

11-06-2014
01:09 AM

62 Views

asd__asdqwe

Beginner

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

11-16-2014
10:10 AM

62 Views

With the new update, I'm now getting errors that I can't reproduce on a simple program such as:

*** Error in PARDISO ( insufficient_memory) error_num= 10

*** Error in PARDISO memory allocation: SOLVING_ITERREF_WORK_DATA, allocation of 1 bytes failed

total memory wanted here: 159 kbyte

The thing is that the matrix is really small (36 x 36 on 8 processes), and I'm turning iterative refinement off, so I don't understand why is PARDISO trying to allocate memory, and why it is failing (I have 6 GB of RAM on the machine I'm getting these errors on).

asd__asdqwe

Beginner

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

11-16-2014
10:50 AM

62 Views

Here is the exact same test case as above, except this time, the RHS and the solution are centralized while the matrix is distributed. Clearly the problem is still here. Can you reproduce this problem on your side ?

Thank you.

$ icc -v

icc version 15.0.1 (gcc version 4.9.0 compatibility)

$ mpicc -cc=icc cl_solver_unsym_distr_c.c -lmkl_intel_thread -lmkl_core -lmkl_intel_lp64 -liomp5

$ mpirun -np 2 ./a.out

The solution out-of-place of the system is:

on zero process x [0] = 0.149579 rhs [0] = 1.000000

on zero process x [1] = 0.259831 rhs [1] = 1.000000

on zero process x [2] = -0.370084 rhs [2] = 0.250000

on zero process x [3] = 0.011236 rhs [3] = 1.000000

on zero process x [4] = 0.415730 rhs [4] = 1.000000

Solving system in-place...

The solution in-place of the system is:

on zero process x [0] = 0.149579

on zero process x [1] = 0.259831

on zero process x [2] = -0.370084

on zero process x [3] = 1.000000

on zero process x [4] = 1.000000

asd__asdqwe

Beginner

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

12-01-2014
12:46 AM

62 Views

asd__asdqwe

Beginner

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

12-08-2014
12:21 AM

62 Views

It has been three weeks since I posted this new issue, but I still don't have a clear answer, could you please let me know if the problem is on my side ?

Thank you.

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.