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)