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
|
function test227
%TEST227 test kron
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2025, All Rights Reserved.
% SPDX-License-Identifier: Apache-2.0
[binops, ~, ~, ~, ~, ~] = GB_spec_opsall ;
binops = binops.all ;
fprintf ('-------------- tests of GrB_kronecker:\n') ;
rng ('default') ;
dnn = struct ;
dtn = struct ( 'inp0', 'tran' ) ;
dnt = struct ( 'inp1', 'tran' ) ;
dtt = struct ( 'inp0', 'tran', 'inp1', 'tran' ) ;
types = { 'int32', 'int64', 'single', 'double' } ;
am = 5 ;
an = 3 ;
bm = 4 ;
bn = 2 ;
Ax = sparse (100 * sprandn (am,an, 0.5)) ;
Bx = sparse (100 * sprandn (bm,bn, 0.5)) ;
cm = am * bm ;
cn = an * bn ;
Cx = sparse (cm,cn) ;
AT = Ax' ;
BT = Bx' ;
for k2 = [4 7 45:52 ]
for k1 = 1:4
type = types {k1} ;
binop = binops {k2} ;
op.opname = binop ;
op.optype = type ;
if (GB_spec_is_positional (op.opname))
if (~(isequal (type, 'int32') || isequal (type, 'int64')))
continue
end
end
fprintf ('[ %s %s ] ', binop, type) ;
for A_is_hyper = 0:1
for A_is_csc = 0:1
for B_is_hyper = 0:1
for B_is_csc = 0:1
for C_is_csc = 0:1
fprintf ('.') ;
clear A
A.matrix = Ax ;
A.is_hyper = A_is_hyper ;
A.is_csc = A_is_csc ;
clear B
B.matrix = Bx ;
B.is_hyper = B_is_hyper ;
B.is_csc = B_is_csc ;
clear C
C.matrix = Cx ;
C.is_csc = C_is_csc ;
%---------------------------------------
% kron(A,B)
%---------------------------------------
% C = kron(A,B)
C0 = GB_spec_kron (C, [ ], [ ], op, A, B, dnn) ;
C1 = GB_mex_kron (C, [ ], [ ], op, A, B, dnn) ;
GB_spec_compare (C0, C1) ;
%---------------------------------------
% kron(A',B)
%---------------------------------------
% C = kron(A',B), no Mask
C0 = GB_spec_kron (C, [ ], [ ], op, AT, B, dtn) ;
C1 = GB_mex_kron (C, [ ], [ ], op, AT, B, dtn) ;
GB_spec_compare (C0, C1) ;
%---------------------------------------
% kron(A,B')
%---------------------------------------
% no mask
C0 = GB_spec_kron (C, [ ], [ ], op, A, BT, dnt) ;
C1 = GB_mex_kron (C, [ ], [ ], op, A, BT, dnt) ;
GB_spec_compare (C0, C1) ;
%---------------------------------------
% kron(A',B')
%---------------------------------------
% no Mask
C0 = GB_spec_kron (C, [ ], [ ], op, AT, BT, dtt) ;
C1 = GB_mex_kron (C, [ ], [ ], op, AT, BT, dtt) ;
GB_spec_compare (C0, C1) ;
end
end
end
end
end
fprintf ('\n') ;
end
end
% test the empty case
clear op
op.opname = 'times' ;
op.optype = 'double' ;
clear A C
A.matrix = sparse (10,10) ;
A.pattern = zeros (10) ;
C.matrix = sparse (100,100) ;
C.pattern = zeros (100) ;
C0 = GB_spec_kron (C, [ ], [ ], op, A, A, dnn) ;
C1 = GB_mex_kron (C, [ ], [ ], op, A, A, dnn) ;
GB_spec_compare (C0, C1) ;
fprintf ('\ntest227: all tests passed\n') ;
|
