1%N2 script to test CHOLMOD septree function
2% Example:
3%   n2
4% See also cholmod_test
5
6% Copyright 2007, Timothy A. Davis, http://www.suitesparse.com
7
8index = ssget ;
9f = find ((index.amd_lnz > 0) & (index.nrows > 200)) ;
10[ignore i] = sort (index.amd_lnz (f)) ;
11f = f (i) ;
12nmat = length (f) ;
13
14for i = f
15
16    Prob = ssget (i, index) ;
17    disp (Prob) ;
18    A = spones (Prob.A) ;
19    [m n] = size (A) ;
20    name = Prob.name ;
21    clear Prob
22
23    if (m == n)
24	mode = 'sym' ;
25	A = A + A' ;
26	len = n ;
27    elseif (m < n)
28	mode = 'row' ;
29	len = m ;
30    else
31	mode = 'col' ;
32	len = n ;
33    end
34
35    [p cp cmem] = nesdis (A, mode) ;
36
37    subplot (2,4,1) ;
38    treeplot (cp) ;
39
40    [cp2 cmem2] = septree (cp, cmem, 0.5, 200) ;	    %#ok
41    subplot (2,4,2) ;
42    treeplot (cp2) ;
43
44    [cp3 cmem3] = septree (cp, cmem, 0.2, 300) ;	    %#ok
45    subplot (2,4,3) ;
46    treeplot (cp3) ;
47
48    [cp4 cmem4] = septree (cp, cmem, 0.12, 500) ;	    %#ok
49    subplot (2,4,4) ;
50    treeplot (cp4) ;
51
52
53    [p cp cmem] = nesdis (A, mode, [200 1]) ;
54
55    subplot (2,4,5) ;
56    treeplot (cp) ;
57
58    [cp2 cmem2] = septree (cp, cmem, 0.5, 200) ;	    %#ok
59    subplot (2,4,6) ;
60    treeplot (cp2) ;
61
62    [cp3 cmem3] = septree (cp, cmem, 0.2, 300) ;	    %#ok
63    subplot (2,4,7) ;
64    treeplot (cp3) ;
65
66    [cp4 cmem4] = septree (cp, cmem, 0.12, 500) ;	    %#ok
67    subplot (2,4,8) ;
68    treeplot (cp4) ;
69
70    drawnow
71    % pause
72
73end
74
75