Bug in ScaLAPACK 2.0.1 and 2.0.2 PDLAHQR() implementation

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

Bug in ScaLAPACK 2.0.1 and 2.0.2 PDLAHQR() implementation

Postby nvinson2utk » Wed May 02, 2012 3:21 pm

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.
nvinson2utk
 
Posts: 1
Joined: Wed May 02, 2012 3:11 pm

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

Postby admin » Wed May 02, 2012 9:27 pm

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
admin
Site Admin
 
Posts: 501
Joined: Wed Dec 08, 2004 7:07 pm

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

Postby lowe980 » Mon May 07, 2012 8:46 am

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?
lowe980
 
Posts: 1
Joined: Mon May 07, 2012 8:42 am

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

Postby admin » Mon May 07, 2012 11:53 am

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
admin
Site Admin
 
Posts: 501
Joined: Wed Dec 08, 2004 7:07 pm


Return to Bug report

Who is online

Users browsing this forum: No registered users and 1 guest