Page 1 of 1

Bug in ScaLAPACK 2.0.1 and 2.0.2 PDLAHQR() implementation

PostPosted: Wed May 02, 2012 3:21 pm
by nvinson2utk
According to http://www.netlib.org/scalapack/html/src/pdlahqr.f when LWORK is set to -1 and passed into PDLAQHR() the function is supposed to return MAX( 2*MAX(DESCZ(LLD_),DESCA(LLD_)) + 2*LOCq(N), 7*Ceil(N/HBL)/LCM(NPROW,NPCOL)) + MAX( 2*N, (8*LCM(NPROW,NPCOL)+2)**2 ) via the WORK variable. However, current implementations just error with a message indicating that argument 15 has an invalid value.

The documentation should be updated to indicate that querying PDLAHQR() is no longer allowed or (preferrably) the definition of PDLAHQR() should be updated to allow querying again.

Re: Bug in ScaLAPACK 2.0.1 and 2.0.2 PDLAHQR() implementatio

PostPosted: Wed May 02, 2012 9:27 pm
by admin
Yes you're right, PDLAHQR does not allow anymore workspace query, but actually the comment has already been modified.
Your link pointed to ScaLAPACK version 1.5. Those files are going to be deleted as they are no longer relevant.

The current html documentation for ScaLAPACK 2.0.2 is here:
http://www.netlib.org/scalapack/explore ... qr_8f.html for interface
http://www.netlib.org/scalapack/explore ... ource.html for source code

We are planning to further integrate Doxygen format into ScaLAPACK like we did for LAPACK.
Let me know of you have any other concern,
Julie

Re: Bug in ScaLAPACK 2.0.1 and 2.0.2 PDLAHQR() implementatio

PostPosted: Mon May 07, 2012 8:46 am
by lowe980
Good day,

Does the ScaLAPACK 2.0.1 include bug fixes to the PxHSEQR Nonsymmetric Eigenvalue Problem routines? If so, where could i look for bug patches?

Re: Bug in ScaLAPACK 2.0.1 and 2.0.2 PDLAHQR() implementatio

PostPosted: Mon May 07, 2012 11:53 am
by admin
ScaLAPACK 2.0.1 and 2.0.2 include bug fixes on PxHSEQR.
Please look at http://netlib.org/scalapack/scalapack-2.0.2.html
and http://netlib.org/scalapack/scalapack-2.0.1.html for mor info