[Octave-bug-tracker] [bug #57185] svds returns inconsistent singular values

[Rik]

Follow-up Comment #8, bug #57185 (project octave):

For the second question, I think we need to check the code.  In svds.m there
is


    if (nargout > 1)
      [V, s, flag] = eigs ([sparse(m,m), b; b', sparse(n,n)],
                           b_k, b_sigma, b_opts);
      s = diag (s);
    else
      s = eigs ([sparse(m,m), b; b', sparse(n,n)], b_k, b_sigma, b_opts);
    endif


so we are calling eigs with a different function prototype in each case.  I
seem to recall that if only the singular values are requested we may perform a
sort on them to ensure uniform ordering.  But, we can't do that arbitrarily
when someone has requested the full decomposition because it would also
require changing the U and V matrices to ensure that the decomposition


A = U * S * V'


still holds.  


    _______________________________________________________

Reply to this item at:

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

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