topic User error or bug in p?geevx ? in Intel® oneAPI Math Kernel Library

User error or bug in p?geevx ?

MOZ — Sat, 25 Jul 2020 01:43:32 GMT

Dear MKL Community,

I have encountered a problem while converting C code that uses LAPACK into code that will use ScaLAPACK. I need to find the right eigenvectors of a non-symmetric real matrix and both sets of eigenvectors (along with eigenvalues) of a non-symmetric, non-hermitian complex matrix.

It is tempting to jump to the conclusion that there is a bug in the relatively newly released p?geevx routines. But I realize that it probably is more likely that I misunderstood the documentation and am sending one or more incorrect parameters to routines or mis-initializing arrays.

I would greatly appreciate it if a more experienced user could take a look at my code to help me find the problem (or problems).

As part of the process of debugging, I stripped the code to more-or-less the bare minimum needed to produce the undesirable behavior. I am posting the real version as pdtest.c and the complex version as pztest.c. The compiler switches I used to compile the code are provided in comments at the beginning of these files.

I am compiling using the gnu C compiler on Ubuntu 16.04 LTS. I have tried both 2020 releases of Intel MKL. I'm using the 2019.7.217 version of the Intel MPI.

The pdtest.c code appears to work with matrix size=12. Example output is provided in the file output01.txt. However, with size=14, the application aborts with calls to MPI_Abort (as seen in output02.txt).

The pztest.c code appears to get stuck in an endless loop, never returning from the call to pzgeevx. (By never I mean about an hour.) Example output is provided in output03.txt. It may be helpful to know that top shows 100% CPU usage by the processes (see top-output.txt).

Thank you in advance for your help with this problem!

Yours sincerely,

Matthew

Re:User error or bug in p?geevx ?

Khang_N_Intel — Tue, 15 Jun 2021 21:48:49 GMT

Error confirmed with oneMKL 2021.2. Escalated!

Re:User error or bug in p?geevx ?

Khang_N_Intel — Wed, 14 Jul 2021 00:32:42 GMT

Hi Mathew,

The problem was in the compiler.

This issue is fixed in the next release of oneMKL, 2021.4.

Best,

Khang

Re: Re:User error or bug in p?geevx ?

MOZ — Sat, 24 Jul 2021 00:08:08 GMT

Dear Khang,

This is good news. I'm looking forward to the release of 2021.4, then!

Yours sincerely,
Matthew

Re:User error or bug in p?geevx ?

Khang_N_Intel — Mon, 26 Jul 2021 18:51:37 GMT

Hi Mathew,

oneMKL 2021.4 will be released sometime in September this year, unless there are some kinds of delay.

Khang

Re:User error or bug in p?geevx ?

Khang_N_Intel — Thu, 30 Sep 2021 20:02:04 GMT

Hi Matthew,

oneMKL 2021.4.will be released tomorrow.

Please install the 2021.4 version when it becomes available and let me know if the issue is fixed from your side.

Best regards,

Khang

Re: Re:User error or bug in p?geevx ?

Gennady_F_Intel — Fri, 01 Oct 2021 03:17:22 GMT

Intel oneAPI Math Kernel Library (oneMKL) version 2021.4 is now available

Re:User error or bug in p?geevx ?

Khang_N_Intel — Mon, 11 Oct 2021 22:43:13 GMT

Hi Matthew,

We haven't heard from you for a while.

Since the solution has been provided in oneMKL 2021.4, I would like to close this thread.

There will be no more communication on this thread.

Best regards,

Khang

Re:User error or bug in p?geevx ?

Khang_N_Intel — Fri, 15 Oct 2021 05:14:17 GMT

Hi Matthew,

I apologize for the previous post.

Could you install oneMKL 2021.4 and check to see if the issue is fixed?

Please let me know as soon as possible so that we can close this thread?

Thanks,

Khang

Re:User error or bug in p?geevx ?

Khang_N_Intel — Thu, 28 Oct 2021 17:46:56 GMT

Hi Matthew,

I haven't heard from you for a while.

Since the solution has been provided, I would like to close this thread.

There will be no more communication on this thread.

Best regards,

Khang