Mercurial > hg > octave-lyh

diff scripts/sparse/svds.m @ 12649:8579e7a6e6e0

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Improve performance (2X) of svds by avoiding issymmetric calculation (bug #33073) * svds.m: Pass issym option to eigs() call. * eigs.cc: Avoid issymetric calculation if 'issym' option given.
author Rik <octave@nomad.inbox5.com>
date Thu, 05 May 2011 16:58:43 -0700
parents c2e8c9d9e284
children 9529392f020e
line wrap: on
line diff
--- a/scripts/sparse/svds.m
+++ b/scripts/sparse/svds.m
@@ -150,6 +150,8 @@
     ## Scale everything by the 1-norm to make things more stable.
     b = A / max_a;
     b_opts = opts;
+    ## Call to eigs is always a symmetric matrix by construction
+    b_opts.issym = true;
     b_opts.tol = opts.tol / max_a;
     b_sigma = sigma;
     if (!ischar (b_sigma))