Community
cancel
Showing results for 
Search instead for 
Did you mean: 
dahlberg
Beginner
68 Views

Pardiso settings to better avoid errmsg -4?

Hi,

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
0 Kudos
2 Replies
68 Views

Hi Carl,
Thank you for reporting the issue. I will look at it as soon as possible. I'll keep you posted!
Regards,
Konstantin
dahlberg
Beginner
68 Views

Hi Konstantin

Any news? My problem persists - unfortunately...

Best regards,
Carl
Reply