Osni Marques
As one can see, the script calls DSYEVR with JOBZ='N', RANGE='I' , and UPLO='L'. This indeed leads to a problem in DSTEBZ (bisection), a situation that I have not observed...
I am unable to reproduce the problem for this matrix of dimension 11, with gfortran and LAPACK 3.9.0. I would have to know the version of LAPACK that is being...
We (Beresford Parlett and I) have implemented a new (and simpler) version of DLARRF, which is called by DLARRV, and which passed all our exhaustive tests. However, we would like...
We do have a (Fortran) standalone test for the *lasq* family. We are mostly interested in accuracy and speed (in comparison with QR, for example). Is that what you mean...
Please send me a case (matrix) for which you have found the problem and I will try to reproduce it.
Thank you, Dmitry. This week is bad. I will take a look at this next week. Osni
Ooops. Clicked on the wrong button. Issue reopened.
I confirm there is a bug in the code. It is not propagating an error message (INFO>0) from DBDSVDX. The 6x4 matrix has two very small singular values, one of...
The problem has been identified and potentially fixed. More tests need to be performed before committing the changes.
This bug is related to the splitting criterion used in DSTEMR, starting at line 606 in the corresponding source file (as in LAPACK 3.7.0). Based on the input variable TRYRAC,...