t-mitchell

Can you please clarify which algorithm is used when N >= 26 but no vectors are requested? In that case, the singular values seem to be computed correctly. This would...

Martin Köhler and I just answered my own question directly above. When N >= 26 but no vectors are requested, xGESDD does not use D&C. So this behavior seems to...

Martin and I changed the LAPACK tuning parameter so xGESDD would switch to D&C for N >= 10 and the singular vectors requested. What we saw is that the singular...

Another update from Martin and I. We also added gsvd (MATLAB) and DGGSVD to the demos, to compare to the generalized singular values of A,I, which should agree. In this...

As a more cautious alternative to taking the absolute value before returning the singular values, how about just explicitly setting any exactly zero singular value to +0 before returning? The...