Incorrect ZUPMTR call inside ZHPEVX

[MaksimZh]

I'm trying to get few eigenvalues and corresponding eigenvectors of some 488x488 matrix with ZHPEVX routine.
Everything works fine if I ask for eingenvalues only.
If I try to find eigenvectors the following error message appears:

Code: Select all

 ** On entry to ZUPMTR parameter number  9 had an illegal value


This happens only if the number of requested eigenvectors is larger then 8.
This does not happen for every matrix, but quite easy to reproduce.

PS:
LAPACK version: 3.3.1
Platform: Linux 64-bit

[rodney]

Have you tried using the latest version of LAPACK, version 3.4.1?

Rodney

[MaksimZh]

The same error for LAPACK-3.4.1

It looks like only size of the matrix matters (488 x 488).
I see the same error for zero matrix, identity matrix and for a matrix where each element is equal to 1.

[rodney]

I have been unable to reproduce the problem you describe, even using n=488. Could you provide a small example program that demonstrates the problem?

--Rodney

[MaksimZh]

There is a bug in the wrapping function I used.
The dimension of IFAIL parameter was always 1 instead of N.
Now I call ZHPEVX directly and it works well.