[Octave-bug-tracker] [bug #55564] GESDD: vastly different singular values when vectors are also requested

[Mike Miller]

Update of bug #55564 (project octave):

                  Status:                    None => Confirmed              

    _______________________________________________________

Follow-up Comment #2:

Here is what I've distilled out of the original problem statement on
https://github.com/Reference-LAPACK/lapack/issues/316:


## N ≤ 25 works ok
N = 25;
A = compan (fliplr ([1, 1 ./ cumprod(1:N)]));
s1 = svd (A);
[U, S, V] = svd (A);
s2 = diag (S);
s1
s2

## N ≥ 26 produces different singular values
N = 25;
A = compan (fliplr ([1, 1 ./ cumprod(1:N)]));
s1 = svd (A);
[U, S, V] = svd (A);
s2 = diag (S);
s1
s2



    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?55564>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/