Workspace allocation problems with pdhseqr

Post here if you have a question about LAPACK or ScaLAPACK algorithm or data format

Workspace allocation problems with pdhseqr

Postby Trevor » Fri Mar 23, 2012 5:51 pm

Hello,

I have been having trouble calling workspace query for the routine PDHSEQR, possibly due to a bug in Scalapack. Regardless, is there a way to determine safe sizes for the WORK and IWORK arrays without doing a workspace query?

Best regards,
Trevor
Trevor
 
Posts: 1
Joined: Fri Mar 23, 2012 5:43 pm

Re: Workspace allocation problems with pdhseqr

Postby janK_FHI » Tue Aug 21, 2012 8:02 am

Hi!

I am having the same problem for the integer workspace. Putting -1 into the size decriptor gives on a two processor parallel 9x9 matrix size test problem on proc 0 a 0 and on proc 1 a 72756096. This does not make sense and consequently this routine produces segfaults, strangely they come from PXERBLA,

+ /lib/x86_64-linux-gnu/libc.so.6(+0x364c0) [0x2adc193444c0]
+ function pxerbla_ (0x1714741)
+ function pdhseqr_ (0x1743FDF)
+ ......

which does handle the error output for wrong parameters. Looks like a bug to me and I suspect it to be coming from the integer workspace array, whos size I cannot safely determine.

Is there any other routine that can determine the size of the integer array needed?

The workspace query to the real array size works fine with -1 tho, it returns reasonable numbers.

Have a nice day, Jan
janK_FHI
 
Posts: 1
Joined: Tue Aug 21, 2012 7:28 am


Return to Algorithm / Data

Who is online

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