Page 1 of 1

Incorrect ZUPMTR call inside ZHPEVX

PostPosted: Fri Jun 22, 2012 9:02 am
by 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

Re: Incorrect ZUPMTR call inside ZHPEVX

PostPosted: Tue Jul 03, 2012 2:19 pm
by rodney
Have you tried using the latest version of LAPACK, version 3.4.1?

Rodney

Re: Incorrect ZUPMTR call inside ZHPEVX

PostPosted: Tue Aug 07, 2012 8:37 am
by 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.

Re: Incorrect ZUPMTR call inside ZHPEVX

PostPosted: Tue Aug 07, 2012 3:29 pm
by 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

Re: Incorrect ZUPMTR call inside ZHPEVX

PostPosted: Mon Aug 13, 2012 7:04 am
by 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.