lapack
lapack copied to clipboard
Reference-LAPACK
LAPACK development repository
DLATB4 is called from DLATTB with IMAT negative for lower triangular matrices: ``` IF( UPPER ) THEN CALL DLATB4( PATH, IMAT, N, N, TYPE, KL, KU, ANORM, MODE, $ CNDNUM,...
For routines [DGEJSV](http://www.netlib.org/lapack/explore-html/d1/d7e/group__double_g_esing_ga8767bfcf983f8dc6ef2842029ab25599.html#ga8767bfcf983f8dc6ef2842029ab25599) and [SGEJSV](http://www.netlib.org/lapack/explore-html/d4/dca/group__real_g_esing_ga4b629535d69b782d6421518962e639f1.html#ga4b629535d69b782d6421518962e639f1), there is no way to query the optimal size of workspace (lwork). The lack of this query create significant inconvenience when using these routines. e.g....
The argument qmax is labeled as an 'in' parameter in the doc, but it is potentially changed (line 320), and used in xLASQ2 after the call to xLAQS3. I believe...
In SLAEBZ and DLAEBZ, the argument 'e' (off-diagonal elements) is unused. I believe it can safely be removed, provided the corresponding arguments in xLARRD and xSTEBZ are also removed.
If exactly zero is returned as a singular value (unlikely but still happens sometimes), it may be that it is returned as +0 or -0. While allowed in IEEE floating...
Trying to build lapack with cmake and these options ```shell CMAKE_ARGS+= -DCBLAS=ON \ -DLAPACKE=ON \ -DBUILD_SHARED_LIBS=ON ``` The build succeeds but linking fails because it cannot find the fortran compiler....
Last night I found what appears to be a bad bug when computing a full SVD using xGESDD. On the particular matrix, companion_demo(n) from EigTool with n>=26, the computed singular...
In the quick return code path for N=0, the RCONDE/RCONDV is not set if SENSE is not 'N'. We may need a short piece of code to set them in...
In lapacke_dgeesx.c, line 82-88, it only allocates iwork if sense = 'b' or 'v'. But in LAPACK version of dgeesx.f, it always set IWORK(1) = LIWRK at line 403 &...
I expected that a workspace query would always return a valid LWORK. This is not the case, a workspace query can yield LWORK=0. Should a warning be added to the...
