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