LAPACK error: info -5

4 次查看（过去 30 天）
Robert
Robert 2016-7-28
评论： Robert 2016-8-15
Hi All,
I'm getting an error when running a Simulink model. I'm not sure what its from. As far as I know I'm not specifically calling out an use of LAPACK. Any insight into possible causes would be great!
Thanks
An error occurred while running the simulation and the simulation was terminated Caused by: The LAPACK call to 'LAPACKE_zgeev' failed with info -5. Please contact MathWorks Technical Support if you can reproduce this error.
  3 个评论
显示 1更早的评论
James Tursa
James Tursa 2016-7-28
Can you reproduce the error and have you contacted Technical Support?
Robert
Robert 2016-7-29
编辑：Robert 2016-7-29
Thanks for the thoughtful response. Yes it can be reproduced and yes technical support was contacted. Wanted to see if anyone out here has seen anything close to this issue before (I've never heard of LAPACK until today) and could give any insight into what could possibly be causing it. Since I couldn't find any help through searching online, thought it would be helpful to others in the future if an answer is found.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

James Tursa
James Tursa 2016-7-29
编辑：James Tursa 2016-7-29
BLAS and LAPACK are the math libraries that MATLAB uses for linear algebra. zgeev is the name of one of the routines in the LAPACK library. E.g., a description can be found on netlib.org (showing Fortran code although the actual libraries that MATLAB uses are likely highly optimized C/assembly using the same interface):
http://www.netlib.org/lapack/explore-3.1.1-html/zgeev.f.html
If you look at the comments for this routine, you see this for "info":
* INFO (output) INTEGER
* = 0: successful exit
* < 0: if INFO = -i, the i-th argument had an illegal value.
* > 0: if INFO = i, the QR algorithm failed to compute all the
* eigenvalues, and no eigenvectors have been computed;
* elements and i+1:N of W contain eigenvalues which have
* converged.
So it appears that the 5th argument had an illegal value. If you look at the signature and the code for generating the error you see this:
SUBROUTINE ZGEEV( JOBVL, JOBVR, N, A, LDA, W, VL, LDVL, VR, LDVR,
$ WORK, LWORK, RWORK, INFO )
:
* ZGEEV computes for an N-by-N complex nonsymmetric matrix A, the
* eigenvalues and, optionally, the left and/or right eigenvectors.
:
ELSE IF( LDA.LT.MAX( 1, N ) ) THEN
INFO = -5
So the 5th input argument, LDA, was less than MAX(1,N), generating the error. Based on nothing more than what you have posted, this looks like a bug in the calling routine. So technical support needs to look at the underlying code that is generating this call.
Is there something unusual about the inputs you are using for your apparent eigenvalue calculations? E.g., if LDA = 0 then this would generate this error. And if you look at what LDA is in the code:
COMPLEX*16 A( LDA, * ), VL( LDVL, * ), VR( LDVR, * ),
$ W( * ), WORK( * )
You see that it is the first dimension of one of the input arrays. Are you passing in (i.e. trying to compute the eigenvalues of) a matrix with this first dimension equal to 0?
  3 个评论
显示 1更早的评论
James Tursa
James Tursa 2016-7-29
Hmmm ... that link shows a different interface with different meanings for the error codes (-5 being a NaN check failure instead of a dimension failure). Does your matrix have NaN's in it?
Robert
Robert 2016-8-15
So in the end, it did end up being NaN's going into the eig() function!
According to mathworks support: "Pre-R2015b, passing NaNs to 'eig' in a MATLAB function block did not error out. A matrix with real NaNs and a diagonal NaN matrix were returned as the eigenvalues and vectors. "
Unfortunately, as they pointed out, I think the resulting error was less helpful than usual because the eig() call was inside of function block. So instead of Simulink pointing to where the error was, it simply gave the LAPACK error above.
Thanks for the help! Hope this helps people in the future.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Simulink Functions 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by