LAPACK-3.2.1 errata file
updated on Thu March 18 2010
maintained by J. Langou, U. Colorado Denver
maintained by J. Langou, U. Tennessee
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)
(*) bug0043 :: VL, VU are referenced in lapack-3.2.1 dsyevr even if range.ne.'V'
(*) svn :: r699
o bug report sent by Michael Chuvelev (Intel) on Wed, 14 Oct 2009 to "lapack@cs.utk.edu".
"VL, VU are referenced in lapack-3.2.1 dsyevr even if range.ne.'V'"
o committed by Julien Langou (SVN revision 699, on Tue Oct 20 2009)
(*) bug0044 :: bug in dsyevr.f when N = 1
(*) svn :: r704
o bug report and fix sent by Ashutosh Mahajan (Mathematics and Computer Science Division of Argonne National Labs) on Fri, 11 Dec 2009 to "lapack@cs.utk.edu".
o committed by Julien Langou (SVN revision 704, on Tue Dec 15 2009)
(*) bug0045 :: bug in dstemr.f when N = 2
(*) svn :: r705
o bug report and fix sent by Ashutosh Mahajan (Mathematics and Computer Science Division of Argonne National Labs) on Thu, 10 Dec 2009 to "lapack@cs.utk.edu".
o committed by Julien Langou (SVN revision 705, on Tue Dec 15 2009)
(*) bug0038 :: no LIWORK query in DGELSD
(*) svn :: r698
o bug report from Mathew Cross on Tue 18 Aug 2009
o see forum topic 1601
o committed by Julie Langou (SVN revision 698, on Thu Oct 20 2009)
(*) bug0046 :: Incorrectly documented RWORK workspace in ZGESDD
(*) svn :: r729
o bug report from Zbigniew on Mon 15 Mar 2010
o see forum topic 1779
o committed by Julie Langou (SVN revision 729, on Thu Mar 18 2010)
(*) bug0047 :: incorrect R/C parameter description in *gerfsx
(*) svn :: r726
o bug report sent by Michael Chuvelev (Intel) on Wed, 1 Mar 2010 to "lapack@cs.utk.edu".
o bug fix sent by James Demmel
"R/C are input only by their usage"
"there's no FACT parameter in these routines"
o committed by Julie Langou (SVN revision r726, on Tue Mar 9 2010)
(*) bug0048 :: Hanging could occur in *gebal if a NaN is in input matrix.
(*) svn :: r730
o bug report and fix sent by Alexander Kobotov (Intel) on Wed, 3 Feb 2010 to "lapack@cs.utk.edu".
o committed by Julie Langou (SVN revision r730, on Tue Mar 9 2010)
(*) bug0049 :: output parameter GIVPTR stays uninitialized in case of quick exit (n.eq.0) or if the rank-1 modifier is small enough in *laed8.
(*) svn :: r731
o bug report and fix sent by Alexander Kobotov (Intel) on Wed, 3 Feb 2010 to "lapack@cs.utk.edu".
o committed by Julie Langou (SVN revision r731, on Tue Mar 9 2010)
(*) bug0050 :: An extra argument after matrix A lines 348-249 in STFSM.f causes a crash of SGEMM call.
(*) svn :: r732
o bug report and fix sent by Alexander Kobotov (Intel) on Wed, 3 Feb 2010 to "lapack@cs.utk.edu".
o committed by Julie Langou (SVN revision r732, on Tue Mar 9 2010)
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!)
(*) bug0048 :: should return info=n+3 if xTGSEN fails to swap eigenvalues
o bug report sent by Ondra Kamenik on Mar 2, 2010 to "lapack@cs.utk.edu".
o on line 444 the routine sets info=n+3 upon error in dtgsen.
o However, then it checks whether the selected eigenvalues are only in the upper part.
o If dtgsen has failed, it is almost certain that it is the case that the selected eigenvalues
o are not only in the upper part, and info=n+2 is set. In this way user does not get info=n+3, but only info=n+2.
o confirmed by Julie Langou on Thu Mar 18 2010
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
o Julie : CANNOT REPRODUCE THE PROBLEM
(*) 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.
(*) bug0039 :: faulty workspace 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 15 2009 (LAPACK-3.2.1)
o Sat Oct 10 2009 (LAPACK-3.2.1)
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)