This array is declared with a size of N, the order of the matrix to be decomposed.
However, in the documentation https://software.intel.com/en-us/node/521126
'isuppz' is an array of size at least 2 *max(1, m), where m is the total number of eigenvalues found.
Would it be better to declared 'isuppz' with a size of 2 * N? Thank you.
thanks, Patric. Yes, it probably makes sense but we align with the original lapack's (http://www.netlib.org/lapack/explore-html/d2/d8a/group__double_s_yeigen_gaeed8a131adf56eaa2a9e5b1e0cce5718.html) and then you may submit the similar query against netlib community.