lapack
lapack copied to clipboard
LAPACK development repository
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...
**Description** Under certain conditions the new *LARTG implementations can essentially increase the 2-norm of initial vector. It can lead to significant loss of accuracy of singular vectors or eigenvectors in...
Closes #629 New algorithms for computing Givens rotations --------------------------------------------- @sergey-v-kuznetsov highlighted in #629 that the new Givens rotations operations may have lower accuracy than the ones that were in LAPACK...
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 >=...
**Description, and suggested fix from Professor Ren-Cang Li, University of Texas at Arlington** Consider the case of computing the last updated eigenvalue using `DLAED4`. If `C=0`, `ETA` is currently set...
Solves #523. This PR adds a Github Action with a test for Windows (latest) using the MinGW generator. Since this new code is not in the master branch, please look...
**Description** This PR adds an implementation of an iterative algorithm for the reduction of a pencil to Hessenberg-upper triangular form. It leverages the connection between generalized and standard eigenvalue problems....
**Description** fixes recent CI failure "InvalidVersionSpecError: Invalid version spec =2.7" (see also conda/conda issue 10618) **Checklist** N/A
Prompted by https://github.com/xianyi/OpenBLAS/pull/3478 I had a look at the significance of the LAPACK_NOOPT macro in the build system. (Setting compiler options to "-O0" or equivalent, presumably to work around cases...
**Description** The reference BLAS implementation of [s/d/c/z]trsm seems to have inconsistent behavior between its varying branches when the right-hand side is all zero and the triangular matrix contains some zeros...