Mercurial > hg > octave-lyh

comparison scripts/sparse/spones.m @ 6498:2c85044aa63f

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2007-04-05 17:59:47 by jwe]
author jwe
date Thu, 05 Apr 2007 17:59:47 +0000
parents 4c8a2e4e0717
children 93c65f2a5668
comparison
equal deleted inserted replaced
6497:fc8ed0c77e08 6498:2c85044aa63f
19 ## @deftypefn {Function File} {@var{y} =} spones (@var{x}) 19 ## @deftypefn {Function File} {@var{y} =} spones (@var{x})
20 ## Replace the non-zero entries of @var{x} with ones. This creates a 20 ## Replace the non-zero entries of @var{x} with ones. This creates a
21 ## sparse matrix with the same structure as @var{x}. 21 ## sparse matrix with the same structure as @var{x}.
22 ## @end deftypefn 22 ## @end deftypefn
23 23
24 function s = spones(s) 24 function s = spones (s)
25 if issparse(s) 25 if (issparse (s))
26 [i,j,v,m,n] = spfind(s); 26 [i, j, v, m, n] = spfind (s);
27 else 27 else
28 [i,j,v] = find(s); 28 [i, j, v] = find (s);
29 [m,n] = size(s); 29 [m, n] = size (s);
30 end 30 end
31 s = sparse(i,j,1,m,n); 31 s = sparse (i, j, 1, m, n);
32 endfunction 32 endfunction
33 33
34 %!assert(issparse(spones([1,2;3,0]))) 34 %!assert(issparse(spones([1,2;3,0])))
35 %!assert(spones([1,2;3,0]),sparse([1,1;1,0])) 35 %!assert(spones([1,2;3,0]),sparse([1,1;1,0]))
36 %!assert(spones(sparse([1,2;3,0])),sparse([1,1;1,0])) 36 %!assert(spones(sparse([1,2;3,0])),sparse([1,1;1,0]))