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
- PARDISO Schur Complement BUG

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

mokwinski__yvan

Beginner

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

01-11-2019
02:15 PM

61 Views

PARDISO Schur Complement BUG

Here's the following sparse CSR matrix.

n = 10

ia = {0, 3, 6, 9, 12, 15, 18, 19, 20, 21, 22 };

ja = { 4, 7, 8, 0, 5, 7, 1, 6, 8, 0, 2, 7, 1, 2, 8, 2, 3, 4, 3, 5, 6, 4 };

a = { -1, 1, 1, 1, 1, 5.5, 1, 1, 4, 1, -1, -2, 1, -1, -4, 1, -1, -6, 1, 1, 1, 1 };

perm = {0, 0,0, 0,0,0 ,0,1,1,1}

The exact Schur complement is {-13.5,-6,0,-6,-14,0,1,1,0}.

When calculating the Schur complement, Pardiso is handling a pivoting issue.

A/ If Pardiso perturbs the pivot elements with 1E-13, then I get {-7.5, 1.72432057769409E-10, 0,0,-8,0,0,1,0

B/ If Pardiso perturbs the pivot elements with 1E-6, then I get {-13.5009434324193, -6.0009554327987, 0,-6,-14.0002128945854,0,1,1.00000000046566,0}

If I include all the zeros in the sparse matrix (full matrix with a sparse matrix representation of 100 coefficients) then I get the exact right answer.

If I am doing the factorization of the upper left corner (7x7) of the matrix and calculating myself the Schur complement

then I get the exact right answer.

What's happening? Is it a bug?

Yvan

Link Copied

1 Reply

Kirill_V_Intel

Employee

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

01-15-2019
01:40 PM

61 Views

Hello Yvan,

Could you please provide us also the values of iparm and mtype which you have used? And, also, what was your configuration on which you ran the code? I managed to get different results for the Schur complement with your test matrix but for me only one entry #4 varied significantly with different pivoting.

We'll have a look at it.

Best,

Kirill

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.