ScaLAPACK Archives

[Scalapack] Bugs in integer workspace queries for P{S,D}LAQR1

Dear ScaLAPACK developers,

I have been experimenting with ScaLAPACK's AED implementations and my
drivers were crashing for small problem sizes. I subsequently found and
fixed the problem: P{S,D}LAQR1 is not setting IWORK(1)=3 when the
integer workspace is queries, but it is accessing the first three
entries of IWORK during normal operation.

The reason that this causes crashes for tiny matrices is that large
matrices will result in P{S,D}HSEQR also querying workspaces from
P{S,D}LAQR0, and the integer workspace query will almost certainly be
larger than 3.

You can find my two-line changesets here:

Please let me know if it is better to submit this to the forum instead.


-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 901 bytes
Desc: OpenPGP digital signature

<Prev in Thread] Current Thread [Next in Thread>

For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or