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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
|
function gbtest70
%GBTEST70 test GrB.random
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2022, All Rights Reserved.
% SPDX-License-Identifier: Apache-2.0
rng ('default') ; A = sprand (4, 5, 0.5) ;
rng ('default') ; C0 = sprand (A) ;
rng ('default') ; C1 = GrB.random (A) ;
assert (isequal (C0, C1)) ;
types = gbtest_types ;
rng ('default') ;
for k = 1:length(types)
type = types {k} ;
rng ('default') ;
G = GrB.random (30, 40, 0.6) ; %#ok<*NASGU>
r = gbtest_cast ([3 40], type) ;
G = GrB.random (300, 400, 0.6, 'range', r) ;
assert (isequal (GrB.type (G), type)) ;
if (~isequal (type, 'logical'))
[i,j,x] = find (G) ; %#ok<*ASGLU>
if (isinteger (r))
assert (min (r) == min (r)) ;
assert (max (r) == max (r)) ;
elseif (isreal (r))
d = min (x) - min (r) ; assert (d > 0 && d < 0.01) ;
d = max (r) - max (x) ; assert (d > 0 && d < 0.01) ;
end
end
G = GrB.random (30, 40, 0.6, 'normal') ;
assert (isequal (GrB.type (G), 'double')) ;
G = GrB.random (30, 40, inf, 'normal') ;
assert (isequal (GrB.type (G), 'double')) ;
assert (nnz (G) == prod (size (G))) ; %#ok<*PSIZE>
G = GrB.random (30, 40, 0.6, 'normal', 'range', r) ;
assert (isequal (GrB.type (G), type)) ;
G = GrB.random (30, 40, 0.6, 'uniform') ;
assert (isequal (GrB.type (G), 'double')) ;
G = GrB.random (30, 40, 0.6, 'uniform', 'range', r) ;
assert (isequal (GrB.type (G), type)) ;
G = GrB.random (30, 0.6, 'symmetric') ;
assert (issymmetric (G)) ;
assert (isequal (GrB.type (G), 'double')) ;
G = GrB.random (30, inf, 'symmetric') ;
assert (issymmetric (G)) ;
assert (isequal (GrB.type (G), 'double')) ;
assert (nnz (G) == prod (size (G))) ;
G = GrB.random (30, 30, 0.6, 'unsymmetric') ;
assert (~issymmetric (G)) ;
assert (isequal (GrB.type (G), 'double')) ;
G = GrB.random (30, 0.6, 'normal', 'symmetric') ;
assert (issymmetric (G)) ;
assert (isequal (GrB.type (G), 'double')) ;
G = GrB.random (30, 0.6, 'normal', 'range', r, 'symmetric') ;
assert (issymmetric (G)) ;
assert (isequal (GrB.type (G), type)) ;
G = GrB.random (30, 0.6, 'normal', 'range', r, 'hermitian') ;
assert (ishermitian (G)) ;
assert (isequal (GrB.type (G), type)) ;
S = sprandsym (30, 0.6) ;
G = sprandsym (GrB (S)) ;
assert (isequal (spones (G), spones (S))) ;
assert (issymmetric (G)) ;
assert (isequal (GrB.type (G), 'double')) ;
S = sprandn (30, 40, 0.6) ;
G = sprandn (GrB (S)) ;
assert (isequal (spones (G), spones (S))) ;
assert (isequal (GrB.type (G), 'double')) ;
S = sprand (30, 40, 0.6) ;
G = sprand (GrB (S)) ;
assert (isequal (spones (G), spones (S))) ;
assert (isequal (GrB.type (G), 'double')) ;
G = sprand (10, 12, GrB (0.5)) ;
assert (isequal (GrB.type (G), 'double')) ;
assert (isa (G, 'GrB')) ;
assert (isequal (size (G), [10 12])) ;
G = sprandn (10, 12, GrB (0.5)) ;
assert (isequal (GrB.type (G), 'double')) ;
assert (isa (G, 'GrB')) ;
assert (isequal (size (G), [10 12])) ;
gnz = nnz (G) ;
% nnz (G) is hard to predict because of duplicates
assert (abs (10*12*0.5 - gnz) < 30) ;
G = sprandn (10, 12, GrB (inf)) ;
assert (isequal (GrB.type (G), 'double')) ;
assert (isa (G, 'GrB')) ;
assert (isequal (size (G), [10 12])) ;
assert (nnz (G) == 120) ;
G = sprandsym (10, GrB (0.5)) ;
assert (isequal (GrB.type (G), 'double')) ;
assert (isa (G, 'GrB')) ;
assert (isequal (size (G), [10 10])) ;
gnz = nnz (G) ;
assert (abs (10*10*0.5 - gnz) < 30) ;
G = sprandsym (10, GrB (inf)) ;
assert (isequal (GrB.type (G), 'double')) ;
assert (isa (G, 'GrB')) ;
assert (isequal (size (G), [10 10])) ;
assert (nnz (G) == 100)
assert (GrB.isfull (G)) ;
assert (GrB.isfull (double (G))) ;
assert (GrB.isfull (full (G))) ;
assert (GrB.isfull (full (double (G)))) ;
end
fprintf ('gbtest70: all tests passed\n') ;
|
