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
|
function test14
%TEST14 test cs_droptol
%
% Example:
% test14
% See also: testall
% Copyright 2006-2012, Timothy A. Davis, http://www.suitesparse.com
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
|
