Showing results for

- Intel Community
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
- Info 512 returned by 'dsyevr' in the Intel MKL library

- 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
##

Hi,

We use Intel MKL 10.1.1.019 for numerical calculation in C++. The 'dsyevr' function is used for eigenvalue decomposition. It works as expected in most cases. But for some cases, it returns info = 512. According to the doc, it means internal error in the Intel MKL library. I wonder what info 512 means exactly for the error.

It is a single threaded calculation. The matrix size is4000 x 4000 at most. Any input will be appreciated.

dengzj

Beginner

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

12-15-2010
07:59 PM

6 Views

Info 512 returned by 'dsyevr' in the Intel MKL library

We use Intel MKL 10.1.1.019 for numerical calculation in C++. The 'dsyevr' function is used for eigenvalue decomposition. It works as expected in most cases. But for some cases, it returns info = 512. According to the doc, it means internal error in the Intel MKL library. I wonder what info 512 means exactly for the error.

It is a single threaded calculation. The matrix size is4000 x 4000 at most. Any input will be appreciated.

5 Replies

Highlighted
##

Hello dengzj,

Gennady_F_Intel

Moderator

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

12-15-2010
10:16 PM

6 Views

we had the similar problems with dsyevr which has been fixed in 10.2 update 3 versions.

You can check if this issue persists into thelatest10.3 version. ( please try the evaluation version for checking the problem and let us know if anyfurtherproblems).

--Gennady

Highlighted
##

Hello Gennady,

Thanks for the reply. We will follow your advice and see if the problem will go away.

Meanwhile, I wonder if some hints could be given onwhat kind of error it is. Any chances that it could be caused by user errors? Orany precautionsthe users cantakein their code tolower the probability of stepping into the same situation again?

Anyhow, thank you very much for the suggestion.

dengzj

Beginner

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

12-16-2010
11:47 AM

6 Views

Thanks for the reply. We will follow your advice and see if the problem will go away.

Meanwhile, I wonder if some hints could be given onwhat kind of error it is. Any chances that it could be caused by user errors? Orany precautionsthe users cantakein their code tolower the probability of stepping into the same situation again?

Anyhow, thank you very much for the suggestion.

Highlighted
##

Hi Gennady,

Would you mind one more question? If info=511 is returned, is it some kind of error similar to info=512? Will it go away by trying the 10.2 update 3 or later versions? Thanks again.

dengzj

Beginner

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

12-16-2010
12:13 PM

6 Views

Would you mind one more question? If info=511 is returned, is it some kind of error similar to info=512? Will it go away by trying the 10.2 update 3 or later versions? Thanks again.

Highlighted
##

Hello dengzj,

Alexander_K_Intel3

Employee

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

12-17-2010
08:31 AM

6 Views

Behavior of DSYEVR in Intel MKL corresponds to the function available in NETLIB LAPACK. And there isdefinitelycould be a better description for values of INFO. The positive values represents following:

1. For JOBZ='N' this is the INFO value returned by underlying DSTEBZ (0<5, see description of DSTEBZ for more details)

2. For JOBZ='V' this is the INFO value returned by underlyingDSTEIN and means number ofeigenvectors failed to converge. IWORK(2*N:2*N+N-1)stores the indices of eigenvectorsthat corresponding to eigenvectors that fail to converge.

So the reported INFO could relate to how the MRRR algorithm handles your input data.You could try other symmetrical eigensolvers which in general more stable in terms of finding solution to see if they could solve you task: DSYEVD (similar by performance) and DSYEV (slower due higher operations count).

W.B.R.

Alexander

Highlighted
##

dengzj

Beginner

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

12-17-2010
09:29 AM

6 Views

Thank you all for the help.

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