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
43
44
45
46
47
48
49
|
function cspy_test
%CSPY_TEST test cspy and cs_dmspy
% Example
% cspy_test
% See also: testall
% Copyright 2006-2012, Timothy A. Davis, http://www.suitesparse.com
index = ssget ;
[ignore f] = sort (max (index.nrows, index.ncols)) ;
clf
% f = f (523:end) ;
% f = f ((find (f == 938)):end) ;
for i = f
Prob = ssget (i,index) ;
disp (Prob) ;
A = Prob.A ;
try
subplot (1,4,1) ; cspy (A) ;
drawnow
subplot (1,4,2) ; cspy (A,64) ;
drawnow
subplot (1,4,3) ; cs_dmspy (A) ;
drawnow
subplot (1,4,4) ; cs_dmspy (A,0) ;
drawnow
catch
fprintf ('failed...\n') ;
end
[m n] = size (A) ;
if (m == n & nnz (diag (A)) == n) %#ok
p = cs_dmperm (A) ;
if (any (p ~= 1:n))
error ('!') ;
end
[p q r s cc rr] = cs_dmperm (A) ; %#ok
if (any (p ~= q))
error ('not sym!') ;
end
end
drawnow
end
|
