Mercurial > forge

changeset 12649:9c65ab6b3f52 octave-forge

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
introduced a 'demo' section in sparsersb.cc.
author michelemartone
date Sun, 21 Jun 2015 22:26:38 +0000
parents bdd6ce97fa37
children 8d2bdadecce4
files main/sparsersb/src/sparsersb.cc
diffstat 1 files changed, 38 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- a/main/sparsersb/src/sparsersb.cc	Sun Jun 21 20:28:26 2015 +0000
+++ b/main/sparsersb/src/sparsersb.cc	Sun Jun 21 22:26:38 2015 +0000
@@ -2473,3 +2473,41 @@
 %!test
 %! sparsersb(sparsersb(sprand(100,100,0.4)),"autotune","n",20,4,1,1,1)
 */
+
+/*
+%!demo
+%! s=sparsersb([2])
+%! s=sparsersb([1,2],[1,1],[11,21],2,2)
+%! 
+%! M=10000;N=10000;P=1 / M;
+%! s=sparse(sprand(M,N,P));
+%! r=sparsersb(s);
+%! x=ones(M,1);
+%! bdt=1.0;
+%! nnz(s)
+%! nnz(r)
+%! 
+%! tic()
+%! nnz(s);
+%! sc=0;
+%! for i=1,10
+%!   s*x;
+%!   sc=sc+1;
+%! end
+%! toc()
+
+%! tic()
+%! nnz(s);
+%! sc=0;
+%! for i=1,10
+%!   r*x;
+%!   sc=sc+1;
+%! end
+%! sc
+%! toc()
+%! 
+%! # FIXME: this is incomplete.
+*/
+
+
+