Mercurial > hg > octave-nkf

comparison scripts/sparse/svds.m @ 12329:5f203b5bbf98

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use testif to only run some sparse tests when necessary libraries are installed.
author Rik <octave@nomad.inbox5.com>
date Wed, 02 Feb 2011 21:31:33 -0800
parents d25dfa9ed18b
children c2e8c9d9e284
comparison
equal deleted inserted replaced
12328:ad2ee858180b 12329:5f203b5bbf98
256 %!test 256 %!test
257 %! [u2,s2,v2,flag] = svds(A,k); 257 %! [u2,s2,v2,flag] = svds(A,k);
258 %! s2 = diag(s2); 258 %! s2 = diag(s2);
259 %! assert(flag,!1); 259 %! assert(flag,!1);
260 %! assert(s2, s(end:-1:end-k+1), 1e-10); 260 %! assert(s2, s(end:-1:end-k+1), 1e-10);
261 %!test 261 %!testif HAVE_UMFPACK
262 %! [u2,s2,v2,flag] = svds(A,k,0,opts); 262 %! [u2,s2,v2,flag] = svds(A,k,0,opts);
263 %! s2 = diag(s2); 263 %! s2 = diag(s2);
264 %! assert(flag,!1); 264 %! assert(flag,!1);
265 %! assert(s2, s(k:-1:1), 1e-10); 265 %! assert(s2, s(k:-1:1), 1e-10);
266 %!test 266 %!testif HAVE_UMFPACK
267 %! idx = floor(n/2); 267 %! idx = floor(n/2);
268 %! % Don't put sigma right on a singular value or there are convergence issues 268 %! % Don't put sigma right on a singular value or there are convergence issues
269 %! sigma = 0.99*s(idx) + 0.01*s(idx+1); 269 %! sigma = 0.99*s(idx) + 0.01*s(idx+1);
270 %! [u2,s2,v2,flag] = svds(A,k,sigma,opts); 270 %! [u2,s2,v2,flag] = svds(A,k,sigma,opts);
271 %! s2 = diag(s2); 271 %! s2 = diag(s2);