[Octave-bug-tracker] [bug #41426] OOM for issymmetric on sparse matrices

[David Bateman]

Follow-up Comment #2, bug #41426 (project octave):



retval = any ((x != x.')(:));


seems slightly faster for full matrices and about the same speed as the
original code, but the "(:)" creates an out of memory error for large sparse
matrices as well. The difference in speed for the two versions for full
matrices only seems to be about 5% so maybe we should just take your patch as
is otherwise something like


if (issparse (x))
  retval = nnz (x != x.') == 0;
else
  retval = all ((x == x.')(:));
endif


might be appropriate 

D.

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?41426>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/