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
- How to solve this symmetric indefinite matrix with pardiso?

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

Muthulingam_S_1

Beginner

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

05-23-2013
03:04 AM

78 Views

How to solve this symmetric indefinite matrix with pardiso?

Hello,

Greetings! I would like to know how to effectively solve this(attached file) symmetric indefinite system using PARDISO. This particular matrix has super diagonal dominancy and almost zero and negative offdiagonal terms. I tried with differente pardiso parameters, could not get the solution promised by other softwares. Please suggest and list pardiso parameters/options for solving this matrix. The expected(promised) solution is ~0.7 all.

Thank you

Stoka

Link Copied

5 Replies

mecej4

Black Belt

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

05-23-2013
06:13 AM

78 Views

Included among the MKL source code examples are source files (in C, Fortran) to solve symmetric indefinite systems. Adapt one of them to read your problem data and run the program.

I find your phrase "promised solution" incomprehensible. Please explain.

Muthulingam_S_1

Beginner

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

05-23-2013
09:27 AM

78 Views

Yes, I shall explain more, the matrix given in the text file is generated by a finite element program (ANSYS).The finite element program gives output as ~0.70 (promised solution) for all which is correct practically. When I use the same matrix and solve it through MKL PARDISO it gives values ranging from ~0.6 to ~2.0 which is not practical. I just used phase 11,22,33 to solve this matrix. As I explained earlier, the matrix is super dominant along diagonal but have very less values offdiagonal (zero and negative). I would like to know different parameter tweaks in MKL PARDISO to get result similar to the finite element program.

Any help in this direction is greately appreciated.

mecej4

Black Belt

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

05-23-2013
11:05 AM

78 Views

The example code pardiso_sym_f90.f90 that is distributed with MKL, when run on your data, gives a solution vector whose components range from 0.700 to 0.703. No changes to the algorithmic parameters (IPARM array) were needed other than those already present in the example code.

Rather than asking for tweaks, you should pay attention to providing your data correctly to MKL-Pardiso.

Zhang_Z_Intel

Employee

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

05-23-2013
11:27 AM

78 Views

How do you set PARDISO parameters? Do you just use the default settings? There are a few settings in iparm you can try to tweak to get better results (assuming mtype = -2):

iparm[9] - Pivort perturbation

iparm[10] - Scaling

iparm[12] - Maximum weighted matching

See here for detailed discussions on how to set these parameters.

See the code example $MKLROOT/examples/solverc/source/pardiso_sym_getdiag_c.c for hints on how to handle symmetric indefinie matrix.

Muthulingam_S_1

Beginner

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

05-24-2013
08:13 AM

78 Views

Thank you for both of your help. Sorted it by tweaking iparm values as indicated above.

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.