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
- Pardiso: error during symbolic factorization

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

Highlighted
##

I'm using Pardiso to solve a very large sparse linear system and I get the following error message:

pf289

Beginner

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

06-27-2012
12:43 AM

121 Views

Pardiso: error during symbolic factorization

ERROR during symbolic factorization: -2

According to the description of the Pardiso errors, this corresponds to a problem with the memory allocation in Pardiso.

I've doubled the memory available for the computation but I still get the same error message.

Does anyone have any clues on how to get around this?

Thanks!

15 Replies

Highlighted
##

Hi,

Alexander_K_Intel2

Employee

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

06-27-2012
04:36 AM

121 Views

What version of MKL do you use? To answer this question we need to know additional detailes, like values of iparm, szie of matrix, pardiso output with msglvl=1. The best variant is to implement testcase and send it to us to reproduce it on our side.

With best regards,

Alexander Kalinkin

Highlighted
##

Dear Alexander

pf289

Beginner

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

06-28-2012
02:07 AM

121 Views

Thank you very much for your reply. Here is the info you requested:

-We have MKL 10.3 which comes with the Intel C Intel 64 Compiler XEfor applications running on Intel 64, Version 12.1.4.319 Build 20120410

- options in Pardiso:

iparm[0] = 1;

iparm[1] = 2;

iparm[2] = 120;

iparm[3] = 0;

iparm[4] = 0;

iparm[5] = 0;

iparm[6] = 0;

iparm[7] = 2;

iparm[8] = 0;

iparm[9] = 13;

iparm[10] = 1;

iparm[11] = 0;

iparm[12] = 1;

iparm[13] = 0;

iparm[14] = 0;

iparm[15] = 0;

iparm[16] = 0;

iparm[17] = -1;

iparm[18] = -1;

iparm[19] = 0;

maxfct = 1;

mnum = 1;

msglvl = 1;

error = 0;

- Pardiso output:

Size of matrix = 406440692

*** Error in PARDISO ( insufficient_memory) error_num= 0

PARDISO Internationalization error. Message 0 is unknown.

total memory wanted here: 1511 kbyte

symbolic (max): 0 symbolic (permanent): 0

real (including 1 factor): 0

=== PARDISO: solving a real nonsymmetric system ===

Summary: ( reordering phase )

================

Times:

======

Time spent in additional calculations : 0.025764 s

Total time spent : 0.025764 s

Statistics:

===========

< Parallel Direct Factorization with number of processors: > 120

< Numerical Factorization with BLAS3 and O(n) synchronization >

< Linear system Ax = b >

number of equations: 64000

number of non-zeros in A: 406440692

number of non-zeros in A (%): 9.922868

number of right-hand sides: 1

< Factors L and U >

number of columns for each panel: 72

number of independent subgraphs: 0

< Preprocessing with state of the art partitioning metis>

number of supernodes: 0

size of largest supernode: 0

number of non-zeros in L: 0

number of non-zeros in U: 0

number of non-zeros in L+U: 0

ERROR during symbolic factorization: -2

Thank you very much

best

Pau

Highlighted
##

Hi,

Alexander_K_Intel2

Employee

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

06-28-2012
10:12 AM

121 Views

You are correct; its look like
there is not enough memory for correct PARDISO call. Could you provide size of
your RAM? And is it possible to send you matrix with program to me to
investigate it on my side?

With best regards,

Alexander Kalinkin

Highlighted
##

in this case, pls try to use out of core version ( set iparm(60) == 2 ) and let us know ho it will work.

Gennady_F_Intel

Moderator

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

06-30-2012
12:15 AM

121 Views

Highlighted
##

Hi,

pf289

Beginner

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

07-03-2012
02:29 AM

121 Views

I'm using 300GB of RAM memory. I've also tried setting iparm(60)=2 but it makes no difference. I tried to send you the matrix so that you can test it on your side but it turns out that the file is 3GB. Any ideas about how I could send it to you?

Many thanks

Pau

Highlighted
##

Hi,

Alexander_K_Intel2

Employee

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

07-03-2012
03:09 AM

121 Views

Maybe I wrong, but it's seemed that you solve some problem on mesh size 40x40x40. If it's true, is it possible to made testcase on lower mesh size, for example 10x10x10? In such case there are no problems to send such matrix by any possible way.

With best regards,

Alexander Kalinkin

Highlighted
##

Dear Alexander,

pf289

Beginner

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

07-05-2012
01:16 AM

121 Views

In fact I'm solving a 2d problem on two patches with 80x80 grid points each and there are 5 functions to be determined. The point is that for lower resolutions, say up to 70x70 grid points on each patch, Pardiso works perfectly fine. However, when I try to increase the resolution, then Pardiso gives me the error I've reported. I can try to set up a website and you can download the matrix from there. Would this be ok for you? Otherwise I can send you the matrix for a coarser grid if that's any helpful.

best

Pau

Highlighted
##

>> I can try to set up a website and you can download the matrix from there. Would this be ok for you?

Gennady_F_Intel

Moderator

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

07-05-2012
08:21 AM

121 Views

<< yes, please do that and give us the link to this site. it would be also helpful to give us the test how to read these data and run pardiso...

Highlighted
##

Sorry for the delay, I've been away for a few days.

pf289

Beginner

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

07-18-2012
02:18 AM

121 Views

I've set up a website where you can download the matrix and the right hand side of the system to be solved:

fileLich.out contains the matrix itself, fileLichRow.out contains the rows' indices, fileLichCol.out contains the columns' indices, and fileRHS.out contains the RHS. fileLichPara.out contains the number of rows and the total length of the matrix.

There's also a simple c code (ForTest.c) which loads the data and solves the linear system. It returns the error that I've previously reported.

Do let me if you have any problems with any of those files.

Many thanks

Pau

Highlighted
##

I have the same problem, I'm running 3D simulations on cluster with 120 Gb RAM. For the small meshes (200 000 elements) everything works perfectly, but if I'm trying to increase number of elements (6 mllion or more) I am getting the same error=-2 durig factorithation. I'm using default values of iparm and solving symmetric indefinite matrix. May be you know what it could be and how to resolve this problem?

Bit_Stekov

Beginner

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

04-14-2013
02:36 AM

121 Views

Highlighted
##

pf289

Beginner

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

04-15-2013
01:49 AM

121 Views

Unfortunately I wasn't able to solve the problem, so I can't really help you with your problem...

Highlighted
##

Roman_A_Intel

Employee

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

04-16-2013
09:02 PM

121 Views

Hi,

I have solved your problem but I made some changes in your test (see Attachments). Also you could try to use out of core version of Pardiso( iparm[59] = 2 ) or your could try to reduce the number of threads because for big number of threads it requires more additional memory.

Regards,

Roman

Highlighted
##

pf289

Beginner

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

04-17-2013
12:18 AM

121 Views

Thanks a lot!

Pau

Highlighted
##

Hi Roman, I am running into the same problem and am wondering what changes you made to solve this problem? It seems there is no good standard solution anywhere on the internet. Thanks a ton! Felix

uzebeckatrente

Beginner

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

08-06-2020
10:06 PM

108 Views

Highlighted
##

Kirill_V_Intel

Employee

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

08-10-2020
12:06 PM

88 Views

Hi,

Please create a separate post about your particular problem. The overall problem defintion is quite general and the observed error can be caused by different reasons.

Also, please, provide more details about your use case: MKL version, size of the matrix, msglvl output, your iparm settings, your hardware, and what you see exactly. Ideally, a standalone reproducer would be the best as we can then check what exactly happens rather than give you some general piece of advice.

Best,

Kirill

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