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
- mkl error in p?geevx - fortran

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

sottile__francesco

Beginner

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

04-10-2020
06:52 AM

267 Views

mkl error in p?geevx - fortran

I'm trying the newly introduced pcgeevx (complex non-hermitian problem - single precision), and I obtain "Intel MKL ERROR: Parameter 14 was incorrect on entry to PCGEEVX." I do not understand how come (the 14th should be an integer and a reference an integer).

Here is (in my opinion) the concerned part

complex :: eigenval(n)

complex :: eigenvec(lda,ldb), matrix(lda,ldb)

integer :: desca(9)

integer :: lwork, lrwork, liwork, idum

real :: dum

real, allocatable :: rwork(:)

complex, allocatable :: work(:)

complex :: cdum

integer, allocatable :: iwork(:)

lwork = -1 ! this is just the first call to obtain the minimum work size requirement

allocate(work(1))

allocate(rwork(n))

allocate(iwork(n))

call pcgeevx('B', 'N', 'V', 'N', n, matrix, desca, eigenval, cdum, idum, eigenvec, desca, &

idum, idum, rwork, dum, rwork, dum, work, lwork, info)

I would like to stress that the if I use the hermitian equivalent pcheevx or the non parallel version cgeev everything works like a charm.

Is it possible that I didn't get right the arguments in pcgeevx?

thanks in advance

Link Copied

15 Replies

sottile__francesco

Beginner

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

04-14-2020
12:39 AM

267 Views

I have also tried, instead of complex eigenval, to use real and imaginary part arrays, as in

call pcgeevx('B', 'N', 'V', 'N', n, matrix, desca, wr, wi, cdum, idum, eigenvec, desca, &

idum, idum, rwork, dum, rwork, dum, work, lwork, info)

in which case I obtain "Intel MKL ERROR: Parameter 15 was incorrect on entry to PCGEEVX."

Gennady_F_Intel

Moderator

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

04-14-2020
02:06 AM

267 Views

it might be the issue. What version of mkl did you try?

sottile__francesco

Beginner

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

04-14-2020
02:10 AM

267 Views

mkl 2020, update 1 (last month release). only this one has the p?geevx available in the scalapack part.

thanks

Gennady_F_Intel

Moderator

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

04-14-2020
03:31 AM

267 Views

just to keep the time - could you give us the simplest example which we could compile and run?

sottile__francesco

Beginner

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

04-14-2020
07:19 AM

267 Views

Gennady_F_Intel

Moderator

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

04-14-2020
07:26 AM

267 Views

thanks, we will play with the example

Gennady_F_Intel

Moderator

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

04-14-2020
10:34 PM

267 Views

We reproduced the problem. The issue is escalated. The thread will be updated asap.

thanks

sottile__francesco

Beginner

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

05-11-2020
11:19 PM

267 Views

Hello, are there any progress on this subject? thanks.

Gennady_F_Intel

Moderator

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

05-12-2020
01:34 AM

267 Views

yes, the fix of the issue is planning to release the next update of mkl.

Gennady_F_Intel

Moderator

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

07-22-2020
12:54 AM

240 Views

sottile__francesco

Beginner

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

08-25-2020
08:09 AM

191 Views

Hello,

I have checked update 2, but the problem still seems to be there. I have still doubts about the call pgeevx for the error message says "Intel MKL ERROR: Parameter 15 was incorrect on entry to PCGEEVX." as if I make a mistake in the arguments. But I do not understand where (if it is).

thanks, f.

Gennady_F_Intel

Moderator

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

08-25-2020
09:29 AM

184 Views

ok, we will re-checked the case once again and get back to this thread soon.

Gennady_F_Intel

Moderator

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

08-27-2020
12:30 AM

178 Views

The original test contains the error in input parameters. Namely, using nblk to 2 is not correct. According to the documentation (https://software.intel.com/content/www/us/en/develop/documentation/mkl-developer-reference-c/top/sca...) the minimal value of nblk is 6 (see the section "Current Notes and Restrictions"):

The current implementation of p?lahqr requires the distributed block size to be square and at least six (6); unlike simpler codes like LU, this algorithm is extremely sensitive to block size."

The modified example min_test_modv1.F90 ( attached) shows no errors.

Here is the output:

mpiexec -n 1 ./a.out

|- BLACS and DESC initialized

|- Block cyclic mat with dim (and proc) 12 12 0

|- Reference matrix diagonalized sequentially, using lapack

Intel(R) Math Kernel Library Version 2020.0.2 Product Build 20200624 for Intel(R) 64 architecture applications

liwork 184 384

info 0 0

info 0 0

|- diagonalised using SCALAPACK ...

Gennady_F_Intel

Moderator

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

08-27-2020
12:35 AM

176 Views

Gennady_F_Intel

Moderator

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

08-27-2020
12:50 AM

171 Views

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