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
39
40
41
42
|
function test2
%TEST2 test sparse2
% Example:
% test2
% See also cholmod_test
% Copyright 2007, Timothy A. Davis, http://www.suitesparse.com
fprintf ('=================================================================\n');
fprintf ('test2: test sparse2\n') ;
i = [ 2 3 ] %#ok
j = [ 3 4 ] %#ok
s = [11.4 9.2] + 1i * [3.4 1.2] %#ok
sparse (i,j,s) %#ok
sparse2 (i,j,s) %#ok
n = 100 ;
nz = 4000 ;
i = fix (n * rand (nz,1)) + 1 ;
j = fix (n * rand (nz,1)) + 1 ;
s = rand (nz,1) + 1i * rand (nz,1) ;
A = sparse (i,j,s,n,n) ;
B = sparse2 (i,j,s,n,n) ;
nnz(A)
if (norm (A-B,1) > 1e-14)
A_minus_B = A-B %#ok
error ('!') ;
end
C = sparse (A) ;
D = sparse2 (B) ;
if (norm (C-D,1) > 1e-14)
C_minus_D = C-D %#ok
error ('!') ;
end
% spy(C)
fprintf ('test2 passed\n') ;
|
