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 an ill-conditioned equation system using 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

Gators_vs__Sundevil

Beginner

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

05-08-2015
04:39 PM

59 Views

How to solve an ill-conditioned equation system using PARDISO?

Hi,

I'm trying to use PARDISO to solve a linear system with the left hand side being symmetric but very ill-conditioned. I tried to run the solver but with "zero or negative poivting, and A is not SPD matrix" error message when I set mtype = 2. I also tried to change mtype to -2, although the system was solved, the solution didn't make any sense. It's completely wrong.

Can anybody tell me is this system solvable using PARDISO? If so, how to configure the iparm parameters? I have attached the input matrixes. P is the single right hand side. Thanks for any suggestion.

Link Copied

3 Replies

mecej4

Black Belt

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

05-10-2015
09:43 AM

59 Views

You have made some assertions the bases of which may be open to question, and the sheer size of your matrix is a disincentive to try some ideas that might yield a solution.

My questions: How do you know that the matrix is ill-conditioned? Why do you think that the matrix is positive definite? What is the basis for your assertion that the solution with mtype = -2 is "completely wrong"?

Can you prepare a much smaller matrix problem which is beset with the same issues, whose solution may help solve the larger problem?

Gators_vs__Sundevil

Beginner

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

05-11-2015
10:19 AM

59 Views

mecej4 wrote:

You have made some assertions the bases of which may be open to question, and the sheer size of your matrix is a disincentive to try some ideas that might yield a solution.

My questions: How do you know that the matrix is ill-conditioned? Why do you think that the matrix is positive definite? What is the basis for your assertion that the solution with mtype = -2 is "completely wrong"?

Can you prepare a much smaller matrix problem which is beset with the same issues, whose solution may help solve the larger problem?

Hi, Thanks for your reply. I have updated the matrix with a smaller size.

The background information for the equation system is actually a solution of the deformation of a bi-material plate under bending. The elastic properties for the two materials make the solution of this problem totally different. For the case when the difference of the material properties between these two material is small, the resulting K matrix is SPD and solved correctly using PARDISO without any problem. When the difference is very large, say one of the material properties are 100 times smaller than the previous case, the resulting K matrix becomes "ill-conditioned", and PARDISO can not yield a meaningful solution. When mtype = 2, the output information is "zero or negative pivoting, the A matrix is not SPD". Then I changed the mtype to -2. And this is the only option available for my case. For this case, I have a reference solution and that's why I said the solution is completely wrong.

Alemdar__Bulent

Beginner

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

05-15-2015
08:34 AM

59 Views

I believe results would be unpredictable if it is "ill-conditioned". Especially, if conditioned number is close to 16.

-bulent

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