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 Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
- cluster_sparse_solver Schur complement matrix - can't get it to work

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

Laura_S_3

Beginner

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

04-06-2018
08:37 AM

57 Views

cluster_sparse_solver Schur complement matrix - can't get it to work

Hi all,

Has anyone had luck getting the new Schur Complement matrix option (iparm(36)) to work with MKL 2018 release 2? I am not having luck, and can't find any example programs that use this option. If anyone has an example, I'd love to see it. Otherwise I'll post my best attempt at making it work to see if anyone can figure out what is wrong.

Laura

Link Copied

4 Replies

Gennady_F_Intel

Moderator

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

04-06-2018
10:07 AM

57 Views

You may take a look at the mklroot\examples\solverc\source\ pardiso_schur_c.c example

Laura_S_3

Beginner

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

04-11-2018
11:30 AM

57 Views

OK, that helped with the parameters. Thank you!

A couple of notes for future readers:

1) The interface is touchy - some parameters that look innocuous can cause the program not to run. Specifically, initialize the "ignored" parameters maxfct and mnum to 1 and make sure iparm(13) = 0. (At least, that is how the current release is working.)

2) The pardiso_schur_c.c example included with this release uses an asymmetrical matrix but tells Pardiso it is symmetric. So the reported Schur matrix won't match the Schur matrix computed other ways unless you account for this.

Gennady_F_Intel

Moderator

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

04-14-2018
10:38 PM

57 Views

wrt #2: mkl 2018 u2 contains symmetric input:

MKL_INT ia[9] = { 1, 5, 8, 10, 12, 15, 17, 18, 19};

MKL_INT ja[18] = { 1, 3, 6, 7,

2, 3, 5,

3, 8,

4, 7,

5, 6, 7,

6, 8,

7,

8

};

double a[18] = { 7.0, 1.0, 2.0, 7.0,

-4.0, 8.0, 2.0,

1.0, 5.0,

7.0, 9.0,

5.0, 1.0, 5.0,

-1.0, 5.0,

11.0,

5.0

};

Laura_S_3

Beginner

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

04-16-2018
09:29 AM

57 Views

Let me clarify: The matrix A as specified by the input ia, ja, and a above is neither symmetric nor structurally symmetric. However, mtype is set to -2 (real symmetric). This is not a bug in the program, merely a confusing detail.

Edit: Only the upper right triangle is specified, as that is all that is needed, and the lower triangle is the reflection.

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