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 settings to better avoid errmsg -4?

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

dahlberg

Beginner

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

07-06-2011
02:14 AM

68 Views

Pardiso settings to better avoid errmsg -4?

I'm using the Pardiso solver in a research code (Fortran90). It has been working quite well, but since I stepped it up a notch recently it has been giving me a few failures of the type "Zero pivot, numerical factorization or iterative refinement problem", i.e. errmsg = -4.

I previously had some problems with this which was mostly sorted out on this forum in the following thread: http://software.intel.com/en-us/forums/showthread.php?t=81123&o=a&s=lr But now the problem is back - and it is nasty!

The problem I wish to solve is a structurally symmetric (real)linear system with a single rhs. However, I use mtype = 11 due to recommendations from Konstantin Arturov in previous thread. I have Lagrange multipliers in my problem so there are a few zeros on the diagonal, whichI store explicitly. The problem size differs, but is on the order of 20k-100k number of equations, where about 0.25%-2.5% of these are Lagrange constraint(multiplier) equations. The resulting number of non-zerosis typically about 1-5 million, so it is not a crazy huge problem, but not in anyway tiny either.

I suspect that my system matrix may be quite ill-conditioned. (Is there any good way to determine that?) I have been using the solver defaults, i.e.iparm(1) = 0. Since my linear algebra skills are not as sharp as I would like, and my knowledge of sparse matrix methods areclose tonon-existant I was hoping to get some pointers as to what I may try to change in the settings to make this problem go away?

Best regards,

Carl

Link Copied

2 Replies

Konstantin_A_Intel

Employee

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

07-06-2011
03:41 AM

68 Views

Thank you for reporting the issue. I will look at it as soon as possible. I'll keep you posted!

Regards,

Konstantin

dahlberg

Beginner

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

07-15-2011
11:43 AM

68 Views

Hi Konstantin

Any news? My problem persists - unfortunately...

Best regards,

Carl

Any news? My problem persists - unfortunately...

Best regards,

Carl

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