1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38
|
function test14
%TEST14 test cs_droptol
%
% Example:
% test14
% See also: testall
% Copyright 2006-2007, Timothy A. Davis.
% http://www.cise.ufl.edu/research/sparse
rand ('state', 0) ;
for trial = 1:100
m = fix (100 * rand (1)) ;
n = fix (100 * rand (1)) ;
d = 0.1*rand (1) ;
A = sprandn (m,n,d) ;
[i j x] = find (A) ;
A = sparse (i,j,2*x-1) ;
fprintf ('test14 m %3d n %3d nz %d\n', m, n, nnz (A)) ;
% using CSparse
tol = 0.5 ;
B = cs_droptol (A, tol) ;
% using MATLAB
A = A .* (abs (A) > tol) ;
% [m n] = size (A) ;
% s = abs (A) > tol ;
% [i j] = find (s) ;
% x = A (find (s)) ;
% A = sparse (i, j, x, m, n) ;
if (norm (A-B,1) > 0)
error ('!') ;
end
end
|