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
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
function gbtest39
%GBTEST39 test amd, colamd, symamd, symrcm, dmperm, etree
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2022, All Rights Reserved.
% SPDX-License-Identifier: Apache-2.0
for trial = 1:40
fprintf ('.') ;
n = 20 ;
A = sprand (n, n, 0.1) ;
S = A + A' ;
G = GrB (A) ;
H = GrB (S) ;
assert (isequal (amd (A), amd (G))) ;
assert (isequal (amd (S), amd (H))) ;
assert (isequal (colamd (A), colamd (G))) ;
assert (isequal (colamd (S), colamd (H))) ;
assert (isequal (symamd (A), symamd (G))) ;
assert (isequal (symamd (S), symamd (H))) ;
assert (isequal (symrcm (A), symrcm (G))) ;
assert (isequal (symrcm (S), symrcm (H))) ;
assert (isequal (etree (A), etree (G))) ;
assert (isequal (etree (S), etree (H))) ;
[p1, post1] = dmperm (A) ;
[p2, post2] = dmperm (G) ;
assert (isequal (p1, p2)) ;
assert (isequal (post1, post2)) ;
[p1, post1] = dmperm (S) ;
[p2, post2] = dmperm (H) ;
assert (isequal (p1, p2)) ;
assert (isequal (post1, post2)) ;
assert (isequal (dmperm (A), dmperm (G))) ;
assert (isequal (dmperm (S), dmperm (H))) ;
[p1, q1, r1, s1, cc1, rr1] = dmperm (A) ;
[p2, q2, r2, s2, cc2, rr2] = dmperm (G) ;
assert (isequal (p1, p2)) ;
assert (isequal (q1, q2)) ;
assert (isequal (r1, r2)) ;
assert (isequal (s1, s2)) ;
assert (isequal (cc1, cc2)) ;
assert (isequal (rr1, rr2)) ;
[p1, q1, r1, s1, cc1, rr1] = dmperm (S) ;
[p2, q2, r2, s2, cc2, rr2] = dmperm (H) ;
assert (isequal (p1, p2)) ;
assert (isequal (q1, q2)) ;
assert (isequal (r1, r2)) ;
assert (isequal (s1, s2)) ;
assert (isequal (cc1, cc2)) ;
assert (isequal (rr1, rr2)) ;
end
fprintf ('\ngbtest39: all tests passed\n') ;
|
