LAPACK-3.2.1 errata file
updated on Thu October 01 2009
maintained by J. Langou, U. Colorado Denver

LAPACK-3.2.1 errata file
CLAPACK-3.2.1 errata file
BLACS-1.1.patch3 errata file
ScaLAPACK-1.8.0 errata file

Wishlist

Corrected in SVN -- to be released in 3.2.2
===========================================

(*) bug0014 :: CSYEQUB/ZSYEQUB/SSYEQUB/DSYEQUB/CHEEQUB/CHEEQUB: EXTERNAL and INTRINSIC list are not accurate (or not even present) 
(*) svn :: r660

    o bug report sent by Alexander V. Kobotov (Intel) on Mon, 6 Apr 2009 to "lapack@cs.utk.edu".
      "c(he/sy)equb: no slamch in external lists and lsame, intrinsics not described"
      Same problem for D and S versions and for the WORK array.
    o committed by Julien Langou (SVN revision 660 on Sun May 10 2009)

(*) bug0015 :: referencing WORK(1) and IWORK(1) when IJOB=0 in CTGSEN/DTGSEN/STGSEN/ZTGSEN contradicts header description
(*) svn :: r658,659

    o bug report sent by Alexander V. Kobotov (Intel) on Mon, 6 Apr 2009 to "lapack@cs.utk.edu".
      "(d/s)tgsen.f: iwork(1) always referenced: line 455: iwork( 1 ) = LIWMIN, documentation says that
      if IJOB=0 it shouldn't, so NULL pointer causes a segfault."
      Same problem for C and Z versions and for the WORK array.
    o committed by Julien Langou (SVN revision 658 and 659, on Sun May 10 2009)

(*) bug0016 :: problem in workspace query for CGEESX/DGEESX/SGEESX/ZGEESX
(*) svn :: r657

    o bug report sent by Alexander V. Kobotov (Intel) on Mon, 6 Apr 2009 to "lapack@cs.utk.edu".
      "(c/z)geesx: no lwork=-1 branch at all, info = -15 is returned while doing lquery."
      The S and D versions were not exiting after WORKSPACE query as well.
    o committed by Julien Langou (SVN revision 657, on Sun May 10 2009)

(*) bug0017 :: in CHBGST/ZHBGST, reference to uninitialized elements of RWORK in the DO 240 and DO 470 loops
(*) svn :: r656

    o patch given by Lawrence Mulholland (NAG) on Wed May 06 2009 
    o see forum topic 1408
    o committed by Julien Langou (SVN revision 656, on Sun May 10 2009)

(*) bug0019 :: 'make install' generates a "*** Error: Time for operations was zero" too easily
(*) bug0019 :: negative flop counts can be generated by 'make install'
(*) svn :: r693
    o reported by Sergei Steshenko on Thu Nov 20 2008
    o see: emails
    o committed by Julie Langou (SVN revision 693, on Fri Sep 11 2009)
    
(*) bug0021 :: in DGESDD, workspace query gives a value smaller than the minimal value given in the header to run the routine
(*) svn :: r649
    o reported by Guy Bencteux on Sat Dec 06 2008
    o see forum topic 846
    o committed by Julie Langou (SVN revision 649, on Tue Mar 17 2009)

(*) bug0023 :: SROTMG and DROTMG uses deprecated Fortran ASSIGN statement and assigned GOTO statement
(*) svn :: r691

    o reported by Marco Atzeri on Fri Mar 06 2009
    o see forum topic 1420
    o committed by Julie Langou (SVN revision 691, on Thu Sep 10 2009)
 
(*) bug0031 :: whitespace fix in comments
(*) svn :: r692

    o reported by "Jan" on Tue Jun 02 2009
    o see forum topic 1449
    o committed by Julie Langou (SVN revision 692, on Fri Sep 11 2009)

(*) bug0038 :: out of range element access in DLASET when using DGELSS
(*) svn :: r695

    o bug report sent by Tim Hopkins on Fri 31 Jul 2009
    o see forum topic 1582
    o committed by Julie Langou (SVN revision 695, on Fri Sep 25 2009)

(*) bug0041 :: printing format misuse in ./TESTING/EIG/[s,c,d,z]drvgbx.f
(*) svn :: r696

    o bug report sent by Michael Chuvelev (Intel) on Tue 29 Jul 2009 to "lapack@cs.utk.edu".
      "Bug in lapack-3.2.1/TESTING/EIG/zdrvgbx.f"
    o printing format misuse, this would possibly lead to a segmentation fault in case of errors in zgbsvxx. 
    o committed by Julien Langou (SVN revision 696, on Tue Sep 29 2009)
    
(*) bug0042 :: printing format misuse in ./TESTING/EIG/[s,c,d,z]drvgbx.f
(*) svn :: r697

    o bug report sent by Kevin Wadleigh on Tue 29 Sep 2009 to "lapack@cs.utk.edu".
      "LAPACK routine ieeeck.f and compiler optimizations"
    o The LAPACK function IEEECK should be modified to keep compilers from optimizing away the line
      NAN6 = NAN5*0.0
      and replacing it by NAN6 = 0.0.
      I've seen this happen using two different compilers at their default level of optimization. However, using
      NAN6 = NAN5*ZERO
      as in the attached file keeps this from happening.
    o committed by Julie Langou (SVN revision 697, on Wed Sep 30 2009)
    
Confirmed problems -- bug needs to be fixed!
============================================

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

    o reported by Stephan Ripke (Broad Institute, Boston) on Thu May 21 2009
    o confirmed by Julien Langou on Thu May 21 2009
    o see forum topic 1418

(*) bug0022 :: workspace query bug with DGEJSV (when only left or right singular vectors desired)

    o reported by Allin Cottrell on Tue Jan 15 2009
    o confirmed by Julien Langou on Fri Jan 30 2009
    o see forum topic 866
    o all information sent to Zlatko Drmac

(*) bug0024 :: the notation A' is confusing for complex matrices. Does A' stand for A transpose of A Hermitian transpose?

    o reported by Clint Whaley early July 2009
    o LAPACK is not consistent in it use of the notation for transpose and Hermitian transpose.
    o The use of LAPACK A' is not consistent with the BLAS useage.
    o We note that the Matlab/Scilab/Octave convention is A' stands for
      "conjugate transpose" and A.' stands for "transpose".
    o In the Z routines, LAPACK uses the notation A' 669 times, it means "conjugate
      transpose" most of the time but it means "transpose" a few times (maybe
      10 times max, for BLAS-like routines). The notation A**H, A**T, A**N is
      used 549 times, The notation A^T, A^N, A^H 19 times. (This is just for
      the Z routines.)
    o see Julien's emails
    o vote on the lapackers mailing list:
          A**T (or A^T) and A**H (or A^H): Clint Whaley, Jim Demmel, Sven Hammarling, Antoine Petitet, Daniel Kressner
          A' (for A**H) and A.' for (A**T): Julien Langou
      (I love to be that supported by the lapackers!)


Bug reports -- bug needs to be confirmed!
=========================================

(*) bug0025 :: serious problem in DGESDD: DGESDD returns negative singular values with INFO = 0

    o reported by Guy Bencteux on Sat Dec 06 2008
    o see forum topic 847
    o Julie : CANNOT REPRODUCE THE PROBLEM

(*) bug0026 :: problem when playing with DGEEV and gfortran optimization flags

    o reported by Pauli Virtanen (scipy team) on Sat Dec 06 2008 (email to lapack@cs.utk.edu)
    o see: email from Pauli Virtanen
      see: http://scipy.org/scipy/scipy/ticket/709

(*) bug0027 :: misleading comment in the header of DGHEQZ, the 2x2 block of the real Schur form are not as claimed

    o reported by Vasile Sima on Tue Jan 27 2009
    o see: email from Vasile
    o more information:
      more specifically Vasile is reporting that the comments:

*          [...] H contains the upper quasi-triangular
*          matrix S from the generalized Schur factorization;
*          2-by-2 diagonal blocks (corresponding to complex conjugate
*          pairs of eigenvalues) are returned in standard form, with
*          H(i,i) = H(i+1,i+1) and H(i+1,i)*H(i,i+1) < 0.

     are not correct

(*) bug0028 :: failure of convergence of DHGEQZ for some specific matrix pencils

    o reported by Vasile Sima on Thu Jan 22 2009
    o see: email from Vasile
    o more information at: lapack known issues
     
(*) bug0029 :: problem of WORKSPACE query with DSTEDC

    o reported by Tiago Requeijo on Mon Mar 30 2009
    o see forum topic 1389

(*) bug0030 :: line 641 of DTGSY2 (present as well in STGSY2, ZTGSY2, CTGSY2)

    o reported by Zbigniew Leyk (College Station, TX) on Mon Jun 01, 2009
    o see forum topic 1445

(*) bug0032 :: matrix splitting problem with dstebz and sstebz

    o bug report sent by Phil DeMier (IBM) on Mon 29 Jun 2009 to "lapack@cs.utk.edu".
      "matrix splitting problem with dstebz and sstebz"
    o see Phil DeMier's email and Phil DeMier's test code

(*) bug0033 :: failure of zheevd on a specific matrix

    o bug report sent by Cezary Sliwa (Instytut Fizyki PAN, Poland) on Thu 2 Jul 2009 to "lapack@cs.utk.edu".
    o see Cezary Sliwa's test code, Cezary Sliwa's matrix, and
      Cezary Sliwa's email.

(*) bug0034 :: comments in SSPTRF not correct 

    o clarification asked by "sarnath" on Thu Jul 02 2009
    o confirmation of bug by "sarnath" and Julien Langou on Fri Jul 03 2009
    o The statement "JMAX is the column-index of the largest off-diagonal element in row IMAX" is not correct
    o see forum topic 1553

(*) bug0035 :: dlarfb.f and likes

    o bug report sent by Michael Chuvelev (Intel) on Tue 21 Jul 2009 to "lapack@cs.utk.edu".
      "Bug in LAPACK 3.2 dlarfb.f and likes"
    o see Michael Chuvelev and Julie's emails, q.dat,
      tau.dat, and test_dorgql.f

(*) bug0036 :: scaling in xSTEQR, xSTERF

    o bug report sent by Pat Quillen on Fri 24 Jul 2009 to "lapack@cs.utk.edu".
      "Scaling in xSTEQR, xSTERF"
    o see Pat Quillen's email, quillen dsteqr.f.

(*) bug0037 :: xLARFP and scaling

    o bug report sent by Pat Quillen on Wed 29 Jul 2009 to "lapack@cs.utk.edu".
      "xLARFP and scaling"
    o see Pat Quillen's email, and quillen dlarfp2.f,
    o see Sven Hammarling's email
    o see Pat Quillen's email, quillen dlarfp3.f, and quillen dlarf2.f.

(*) bug0038 :: no LIWORK query in DGELSD

    o bug report from Mathew Cross on Tue 18 Aug 2009 
    o see forum topic 1601

(*) bug0039 :: faulty wkspce queries in GEESX

    o bug report from Mathew Cross on Tue 18 Aug 2009 
    o see forum topic 1602

(*) bug0040 :: Q referenced in real TGSENs when WANTQ .FALSE.

    o bug report from Mathew Cross on Tue 18 Aug 2009 
    o see forum topic 1604

Deprecated routines
===================

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

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

Previous Errata (LAPACK 3.2.x)
==============================

   o Thu Oct 01 2009 (LAPACK-3.2.1)
   o Tue Sep 29 2009 (LAPACK-3.2.1)
   o Fri Sep 25 2009 (LAPACK-3.2.1)
   o Mon Sep 11 2009 (LAPACK-3.2.1)
   o Mon Sep 07 2009 (LAPACK-3.2.1)
   o Sun May 10 2009 (LAPACK-3.2.1)
   o Fri Apr 17 2009 (LAPACK-3.2.1)
   o Tue Mar 10 2009 (LAPACK-3.2.0)
   o Tue Mar 03 2009 (LAPACK-3.2.0)
   o Mon Mar 02 2009 (LAPACK-3.2.0)
   o Mon Feb 01 2009 (LAPACK-3.2.0)