updated on April 20 2012

maintained by Julien Langou, U. Colorado Denver and Julie Langou, U. Tennessee

CLAPACK-3.2.1 errata file

BLACS and ScaLAPACK errata file

Wishlist

RED means BUG WAS CONFIRMED but IT IS STILL HERE

ORANGE means BUG WAS NOT CONFIRMED

GREEN means BUG WAS CORRECTED :

bug0094 :: [S,D]SYEVR and [C,Z]SYEVR routine - workspace indices problem

  • CORRECTED

  • Present in LAPACK 3.4.1

  • bug report by Keita Teranishi (CRAY) on May 18 2012

  • INDIWO should be INDIFL+N. Otherwise, subsequent {S,D}STEIN call takes the same address for IWORK(INDIWO) and INWROK(INDIFL).

  • Sent to LAPACK team

  • corrected by Julie on May 21 2012 in rev 1302

bug0093 :: LAPACKE - LAPACKE_xlange routine - workspace problem

  • CORRECTED

  • Present in LAPACK 3.4.1

  • bug report by Mathieu (PLASMA)

  • Was allocating workspace for Frobenius norm, but not inf norm - This is wrong, only the inf norm needs a workspace.

  • Sent directly yo Julie

  • corrected by Julie on May 11 2012 in rev 1299

bug0061 :: zgehrd.f is overflowing

  • CONFIRMED

  • bug report by Victor Liu on Tuesday August 03rd, 2010

  • zgehrd.f is overflowing the maximum size in bytes of the largest array that

  • can be put on the stack with gfortran (⇒ leads ,e.g., to non-reproducible and uncorrect execution

  • of zgeev in multithreaded context

  • lapackers involved in the bug fix process: Julie and Julien Langou and Mathieu Faverge

  • bug fix by Victor Liu on Wednesday August 04th, 2010

  • see forum topic 1930

  • two outstanding issues: (1) how do we go to fix this?, (2) more LAPACK subroutines might be concerned.

bug0057 :: Non-eigenvectors from DGEEV, problem in balancing

  • CONFIRMED

  • bug report sent by Numpy team on April 22nd, 2010

  • see forum topic 1818

  • see follow-up on lapackers on April 22nd from Jim "bug in DGEHRD" "The balancing is choosing to scale some rows/columns by factors like 1e18, which rather amplifies tiny rounding errors in the (presumably correctly) computed eigenvectors of the balanced matrix, when balancing is undone at the end. The expert driver interface dgeevx would let the user disable balancing, until we figure out whether this is a bug or a feature, I’m not sure which."

bug0056 :: Non-orthogonal eigenvectors returned from DSYEVR

  • CONFIRMED

  • bug report sent by Scipy team on April 22nd

  • see LAPACK Mailing list msg 874

  • Confirmed by Rodney on Tue, 17 Feb 2011

  • " I can reproduce similar problems that the user encountered with non-orthogonal eigenvectors when only

  • a few eigenvalues are requested. This also happens with SSYEVR. If I set ABSTOL to SLAMCH(Safe minimum),

  • the problem goes away in SSYEVR, and improves with DSYEVR but does not completely go away."

bug0033 :: failure of zheevd on a specific matrix

bug0032 :: matrix splitting problem with dstebz and sstebz

  • CONFIRMED

  • bug report sent by Phil DeMier (IBM) on Mon 29 Jun 2009 to "lapack@cs.utk.edu". "matrix splitting problem with dstebz and sstebz"

  • see Phil DeMier’s email and Phil DeMier’s test code

  • Confirmed by Rodney on Tue, 17 Feb 2011

bug0020 :: in xGESDD, workspace query can overflow for large matrices

  • CONFIRMED

  • reported by Stephan Ripke (Broad Institute, Boston) on Thu May 21 2009

  • confirmed by Julien Langou on Thu May 21 2009

  • see forum topic 1418

dead code* :: xLAHRD is now replaced by xLAHR2, xLAHRD is a dead code since LAPACK 3.1

  • xLAHRD will be removed at some point in the future from the distribution