LAPACK BUG LIST (original) (raw)

bug 159 DSTEMR - consolidates bugs 113 and 126

bug 158 problem with DC SVD least squares

bug 157 gesdd , gesvd, and bdsdc Workspace computation

bug 154 [LAPACKE] lapacke_?lantr: unnecessary work array allocation

bug 153 Fix out of bounds error in dgeev

bug 151 Bug in xGETC2, quick return for N=0, N=1

bug 150 [LAPACKE] Fix bug on [sd]lantr which don’t report the result

bug 149 workspace computation for [CZ]GESVDX when a workspace query is done

bug 147 INFO parameter value when LVT is not correct.

bug 146 properly avoid setting the install rpath when installing to system directories

bug 145 Fix upper case issue in CMakeLists

bug 144 inconsistency in the lapacke cmake config files install directory location

bug141 LAPACKE dlantr causes memory corruption

bug139 xSTEDC eigenvalues not sorted

bug138 :: LAPACKE_dtr_nancheck not working for trapezoidal matrix

bug 137 :: Dlange gives inconsistent/incorrect results for RowMajor

bug 134 :: Rook pivoting factorization goes into an infinite loop

bug 132 :: dsyevr bug

bug 129 :: LAPACK_?syconv interface problem - Work is used to return values

bug 128 :: lapacke_?laset has limitations - pb if initialize the matrix with NaNs (Matrix assumed to be initialized by this functions)

bug 127 :: LAPACKE interface to ?ORCSD / ?UNCSD is not fully functional - Work is used to return values

bug125 :: xLAHQR v3.5 does not converge on some matrices while xLAHQR v2 converges

bug124 :: dpstf2 returns incorrect value of rank

bug122 :: Wrong Info code when checking input parameter QSIZ in SLASD7 and DLASD7

bug119 :: CBLAS - missing the floating point constant argument in sdsdotsub

bug113 :: DLARRB() infinite loop bug

bug112 :: DSYEVR does not seem to be thread-safe

bug111 :: incorrect RWORK size for ZGESDD

bug110 :: possible infinite loop in ZLARTG when input vector has NaN

bug0109 :: [xSTEGR_WORK] do not check ldz correctly

bug0108 :: [xDRVRFP] Test was failing

bug0107 :: [CSD routines] fixed array index problems

bug0106 :: [xGEBAL] very large backward error for certain Hessenberg matrices (Watkins)

bug0105 :: [LAPACKE] lapacke.h problem

bug0104 :: dgeqp3 fall over when asked to do a factorization of the last two columns.

bug0103 :: [dlasd4/slasd4] fixes bug where TAU2 can be uninitialized.

bug0102 :: convergence problem with some particular pencils in xGGEV

bug0101 :: [S,D, C, Z]hsein: NaN in H can cause an infinite loop

bug0100 :: [S,D]gesdd correct minimum LWORK size

bug0099 :: [S,D]sgesvj array SVA(N2=0) can be out of the array boundary

bug0098 :: DSYEV may return non-orthogonal eigenvectors

bug0097 :: [S,D]LASQ3 could cause NaNs appearance in ?GESVD and ?GESDD

bug0096 :: ?STEMR sorting is not done for N=2.

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

bug0093 :: LAPACKE - LAPACKE_xlange routine - workspace problem

bug0092 :: xLAED6 avoid div by zero

bug0091 :: xLANHF div by zero

bug0090 :: xGGEV and xGGEVX: Need to unscale if necessary when there is an error in DHGEQZ (QZ iteration failed)

bug0089 :: unused external declarations

bug0088 :: Uninitialized variable bug in ZHERFSX

bug0087 :: Wrong call name for XBLAS xLA_RPVGRW, needs to be xLA_GERPVGRW

bug0086 :: ILAxLR:commit r1006 reverts the change made in commit r977

bug0085 :: xGESVD Problem in Workspace computation

bug0084 :: Bug Ming Gu

bug0083 :: Bug in the testing of CHERFSX function

bug0082 :: [CZ]laed2 Memory problems when COMPZ='I' in DSTEDC

bug0081 :: [SD]laed2 WORKSPACE size is wrong when NCVT = NRU = NCC = 0

bug0080 :: [CZ]BDSQR WORKSPACE size is wrong when NCVT = NRU = NCC = 0

bug0079 :: ILA[SDCZ]LR out of bound access in the DO LOOP for array A.

bug0078 :: DBDSDC problem with orthogonality of U in certain cases

bug0076 :: Problem in ?(sy/he)tri2 when nbmax greater than n

bug0075 :: not orhtogonal: eigenvectors of a symmetric matrix with DSYEV

bug0074 :: ?sysv and ?hesv

bug0073 :: secondtst.f and dsecnd.f

bug0072 :: cchkhs.f and zchkhs.f

bug0071 :: slarre and dlarre

bug0070 :: IWORK dimension is incorrect

bug0069 :: sytri2x is writing on other part of A

bug0068 :: Need to add LDA in *sytrs2

bug0067 :: Need to add LDA in *syswapr

bug0066 :: bugs in {s,d}gejsv

bug0065 :: workspace queries not being used

bug0064 :: got a bug in divide and conquer (non convergence)

bug0063 :: Division by 0 in *GESVJ & *GEVSJ with zero columns

bug0062 :: Array out-of-bounds reference in xLAQR5

bug0061 :: zgehrd.f is overflowing

bug0060 :: Bug in dtgsy2, missing IWORK call

bug0060 :: Typo in TESTING/EIG/dchkee.f and TESTING/EIG/zchkee.f

bug0059 :: ZGELSD crashes on large matrices (eq M=N=NHRS>=88), if real workspace allocated exactly as reported on workspace query.

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

bug0056 :: Non-orthogonal eigenvectors returned from DSYEVR

bug0055 :: Missing arguments descriptions of work and uplo in [csz]syequb

bug0053 :: Fix external declaration issue in testing suite file: cebchvxx.f

bug0052 :: Incorrect NRHS parameter in call to ?SYTRS from ?la_{po,sy,he}rfsx_extended.f.

bug0051 :: INFO=N+2 in xDGGES due to too stringent threshold in xTGEX2.

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.

bug0048 :: Hanging could occur in *gebal if a NaN is in input matrix.

bug0047 :: incorrect R/C parameter description in *gerfsx

bug0046 :: Incorrectly documented RWORK workspace in ZGESDD

bug0045 :: bug in dstemr.f when N = 2

bug0044 :: bug in dsyevr.f when N = 1

bug0043 :: VL, VU are referenced in lapack-3.2.1 dsyevr even if range.ne.V

bug0042 :: printing format misuse in ./TESTING/EIG/[s,c,d,z]drvgbx.f

bug0041 :: printing format misuse in ./TESTING/EIG/[s,c,d,z]drvgbx.f

bug0040 :: Q referenced in real TGSENs when WANTQ .FALSE.

bug0039 :: faulty workspace queries in GEESX

bug0038 :: no LIWORK query in DGELSD

bug0038 :: out of range element access in DLASET when using DGELSS

bug0037 :: xLARFP and scaling

bug0036 :: scaling in xSTEQR, xSTERF

bug0035 :: dlarfb.f and likes

bug0033 :: failure of zheevd on a specific matrix

bug0032 :: matrix splitting problem with dstebz and sstebz

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

bug0029 :: problem of WORKSPACE query with DSTEDC

bug0028 :: failure of convergence of DHGEQZ for some specific matrix pencils

*          [...] 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

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

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

bug0023 :: SROTMG and DROTMG uses deprecated Fortran ASSIGN statement and assigned GOTO statement

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

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

bug0019 :: make install generates a "** Error: Time for operations was zero" too easily

bug0017 :: in CHBGST/ZHBGST, reference to uninitialized elements of RWORK in the DO 240 and DO 470 loops

bug0016 :: problem in workspace query for CGEESX/DGEESX/SGEESX/ZGEESX

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

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

bugXXXX Wrong type for workspace in cpot01 and zpot01

bugXXXX Revert DLARFG to 3.1.0, move DLARFG to DLARFP and improved DLARFP

1) Fix as much as possible the xLARFP routines. This was corrected by  
   Julien (see (r706, r707, r708, r710, r712, r713)) based on Jim's idea.  
2) Back out the changes in LAPACK to go back to using the old reflector  
   (xlarfg) everywhere (e.g. in xGEQRF, xGEQR2, etc.).  
3) Keep the change that saves flops when u has lots of trailing zeros,  
   since this is independently a good idea.  
4) Change the name from xLARFP to xLARFGP and introduce new QR routines  
   (xGEQRFP and xGEQR2P) that uses xlarfp to guarantee a new positive  
   diagonal of R, for who needs it.  
5) This requires to adapt the test suites  
6) Note: only the QR routines have an xGEQRFP. While it seems that quite  
   a few users care for QR factorization with nonnegative diagonal R. Noone  
   seemed to care for this feature on LQ, QL or RQ. So we have left out  
   these routines.  
7) Steps 2 to 6 were done during r755 by Julie.

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