Incorrect ZUPMTR call inside ZHPEVX

Post here if you want to report a bug to the LAPACK team

Incorrect ZUPMTR call inside ZHPEVX

Postby MaksimZh » Fri Jun 22, 2012 9:02 am

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
MaksimZh
 
Posts: 3
Joined: Fri Jun 22, 2012 8:53 am

Re: Incorrect ZUPMTR call inside ZHPEVX

Postby rodney » Tue Jul 03, 2012 2:19 pm

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

Rodney
rodney
 
Posts: 48
Joined: Thu Feb 10, 2011 8:20 pm
Location: Colorado College

Re: Incorrect ZUPMTR call inside ZHPEVX

Postby MaksimZh » Tue Aug 07, 2012 8:37 am

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.
MaksimZh
 
Posts: 3
Joined: Fri Jun 22, 2012 8:53 am

Re: Incorrect ZUPMTR call inside ZHPEVX

Postby rodney » Tue Aug 07, 2012 3:29 pm

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
rodney
 
Posts: 48
Joined: Thu Feb 10, 2011 8:20 pm
Location: Colorado College

Re: Incorrect ZUPMTR call inside ZHPEVX

Postby MaksimZh » Mon Aug 13, 2012 7:04 am

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.
MaksimZh
 
Posts: 3
Joined: Fri Jun 22, 2012 8:53 am


Return to Bug report

Who is online

Users browsing this forum: Bing [Bot] and 1 guest