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
|
function gbtest65
%GBTEST65 test GrB.mis
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2022, All Rights Reserved.
% SPDX-License-Identifier: Apache-2.0
rng ('default') ;
% load west0479 ; %#ok<LOAD>
load west0479_correct ; %#ok<*LOAD>
west0479 = Problem.A ;
A = GrB.offdiag (spones (west0479)) ;
A = A+A' ;
maxisize = 0 ;
n = size (A, 1) ;
for trial = 1:100
if (mod (trial, 4) == 1)
iset = GrB.mis (A, 'check') ;
else
iset = GrB.mis (A) ;
end
% assert that iset is an independent set
p = find (iset) ;
assert (nnz (A (p,p)) == 0) ;
isize = length (p) ;
if (isize > maxisize)
fprintf ('trial %3d: iset size: %d\n', trial, isize) ;
maxisize = isize ;
end
% assert that iset is maximal
q = find (~iset) ;
d = GrB.entries (A (p, q), 'col', 'degree') ; %#ok<FNDSB>
assert (all (d > 0)) ;
end
fprintf ('max independent set found: %d of %d nodes\n', maxisize, n) ;
fprintf ('gbtest65: all tests passed\n') ;
|
