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
- Pardiso : OOC mode + iparm[1]!=0 + phase=13 = error after factorisation

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

Thierry_LE_SOMMER__E

Beginner

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

09-16-2010
01:29 AM

136 Views

Pardiso : OOC mode + iparm[1]!=0 + phase=13 = error after factorisation

I have a problem (Access violation reading) when I use Pardiso in OOC mode and I try to do phases 1 to 3 at the same time (phase parameter = 13). But when I do phase 3 separately, it works.

Here is a summary of my code :

[cpp]iparm[0] = 1; /* No solver default */ iparm[1] = 2; /* Fill-in reordering from METIS */ iparm[59] = 2; /*In-Core or Out-Of-Core */ maxfct = 1; /* Maximum number of numerical factorizations. */ mnum = 1; /* Which factorization to use. */ msglvl = 1; /* Print statistical information to the screen */ error = 0; /* Initialize error flag */ phase = 12; PARDISO (_pt, &maxfct, &mnum, &mtype, &phase, &n, a, ia, ja, &idum, &nrhs, iparm, &msglvl, B, X, &error); phase = 33; PARDISO (_pt, &maxfct, &mnum, &mtype, &phase, &n, a, ia, ja, &idum, &nrhs, iparm, &msglvl, B, X, &error);[/cpp]This code works fine in all cases.

But if I want to do phases 1-3 in a single call (phase=13), I get a "Access violation reading" in Pardiso.

This occurs only on large scale matrix when OOC is really used (files write/read) and if I set iparm[1]=0, it works fine.

Thierry

Link Copied

7 Replies

Gennady_F_Intel

Moderator

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

09-16-2010
02:32 AM

136 Views

Thierry,

is it win64?

what is thenumbers of non-zero elements?

--Gennady

Thierry_LE_SOMMER__E

Beginner

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

09-16-2010
03:20 AM

136 Views

For this problem I tested 2 matrices :

(1) 54243x54243 / 2053206 non-zeros

(2) 408483*408483 / 16082406 non-zeros

To force the OOC mode, I set the environment variable MKL_PARDISO_OOC_MAX_CORE_SIZE to 100 for the first matrix and 3000 for the second matrix.

I use MKL 10.2.6 but I did some more tests with MKL 10.2.5 and 10.3.0 beta. I can't reproduce the problem with MKL 10.2.5 and 10.3.0 beta. So,

I don't see any difference in my compiler/linker parameters between MKL 10.2.5 / MKL 10.2.6 / MKL 10.3.0

Thierry

Gennady_F_Intel

Moderator

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

09-16-2010
03:34 AM

136 Views

I'm sure that probably does not matter, but nevertheless, is this SPD matrixes?

and

is this static or dynamic linking?

--Gennady

Thierry_LE_SOMMER__E

Beginner

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

09-16-2010
03:57 AM

136 Views

For the matrix type, I set mtype=-2. Can I set mtype=2 but I am not sure that my matrices are really SPD ?

Gennady_F_Intel

Moderator

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

09-16-2010
04:08 AM

136 Views

yes, you can try

Gennady_F_Intel

Moderator

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

09-17-2010
09:11 AM

136 Views

Thierry,I am interested in the result of this attempt?

Thierry_LE_SOMMER__E

Beginner

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

09-19-2010
11:57 PM

136 Views

I tried to put mtype = 2. It works for a few tens of iterations, but then I get an error indicating that my MKL matrix is not SPD.

So, I use again mtype=-2

Thierry

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.