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

[Mike Miller]

Follow-up Comment #3, bug #55564 (project octave):

Copy/paste error, this should be a little clearer:


function testit (N)
  A = compan (fliplr ([1, 1 ./ cumprod(1:N)]));
  s1 = svd (A);
  [U, S, V] = svd (A);
  s2 = diag (S);
  assert (s1, s2, -10*eps);
endfunction


testit (25)  ## ok
testit (26)  ## not ok



    _______________________________________________________

Reply to this item at:

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

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