function test14 %TEST14 test cs_droptol % % Example: % test14 % See also: testall % Copyright 2006, 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