langou

(S/C/D/Z)TGSJA are calling FORTRAN INTRINSIC HUGE Shall we replace with (S/D)LAMCH( 'Overflow' ) instead?

ZGELSY uses ILAENV with ZGERQF and ZUNMRQ for workspace computation https://github.com/Reference-LAPACK/lapack/blob/f40d2208b99956aa7913aacb18ac4a39bacae3ad/SRC/zgelsy.f#L266 I think this should be ZTZRZF and ZUNMRZ instead.

We should add an early stop in QRCP. Given a threshold parameter given by a user, perform QRCP until step ``k`` where the norm of the trailing matrix (maybe the...

maybe a bug in DGESDD

2

Zlatko Drmač just sent me an email letting me know of a suspicious workspace size in the comments of DGESDD. In 3.1.1 we had, If JOBZ = 'O’, LWORK >=...

Bug report from David Hough Subject: dlasy2.f: substitution for small pivots does not always help Date: Fri, 31 Oct 2014 13:50:08 -0700 (PDT) The LAPACK TESTING functions dchkee.f ddrves.f dlatme.f...

This is a bug report from David Hough Subject: LAPACK slaln2.f: more accuracy hurts if not consistent Date: Wed, 29 Oct 2014 15:29:20 -0700 (PDT) Another example turned up by...

xROTMG issue & proposed fix

2

See: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=5332

We have an error during our testing. If we do: `./LIN/xlintstz < ztest.in ` We see: ``` Testing COMPLEX16 Linear-Equation-routines-ztest.out *** XERBLA was called from ZHESV_AA_2STAGE with INFO = 7...

DTBSM, triangular band matrix solve with multiple right-hand sides

2

Hi OpenBLAS community, I received an email from a BLAS/LAPACK user and they were wondering why there exists Level 2 DTBSV, and why there does not exist a Level 3...