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
|
function test130
%TEST130 test GrB_apply (hypersparse cases)
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2022, All Rights Reserved.
% SPDX-License-Identifier: Apache-2.0
fprintf ('\ntest130: quick GrB_apply tests\n') ;
rng ('default') ;
m = 8 ;
n = 4 ;
dt = struct ('inp0', 'tran') ;
dr = struct ('outp', 'replace') ;
aclass = 'double' ;
A = GB_spec_random (m, n, 0.3, 100, aclass) ;
Cin = GB_spec_random (m, n, 0.3, 100, aclass) ;
B = GB_spec_random (n, m, 0.3, 100, aclass) ;
cin = GB_mex_cast (0, aclass) ;
Mask = GB_random_mask (m, n, 0.5, true, false) ;
if (isequal (aclass, 'double'))
hrange = [0 1] ;
crange = [0 1] ;
else
hrange = 0 ;
crange = 1 ;
end
for A_is_hyper = hrange
for A_is_csc = crange
for C_is_hyper = hrange
for C_is_csc = crange
for M_is_hyper = hrange
for M_is_csc = crange
A.is_csc = A_is_csc ; A.is_hyper = A_is_hyper ;
Cin.is_csc = C_is_csc ; Cin.is_hyper = C_is_hyper ;
B.is_csc = A_is_csc ; B.is_hyper = A_is_hyper ;
Mask.is_csc = M_is_csc ; Mask.is_hyper = M_is_hyper ;
op.opname = 'ainv' ;
fprintf ('.') ;
op.optype = 'double' ;
% no mask
C1 = GB_spec_apply (Cin, [], [], op, A, []) ;
C2 = GB_mex_apply (Cin, [], [], op, A, []) ;
GB_spec_compare (C1, C2) ;
% no mask, with accum
C1 = GB_spec_apply (Cin, [], 'plus', op, A, []) ;
C2 = GB_mex_apply (Cin, [], 'plus', op, A, []) ;
GB_spec_compare (C1, C2) ;
% with mask
C1 = GB_spec_apply (Cin, Mask, [], op, A, []) ;
C2 = GB_mex_apply (Cin, Mask, [], op, A, []) ;
GB_spec_compare (C1, C2) ;
% with mask and accum
C1 = GB_spec_apply (Cin, Mask, 'plus', op, A, []) ;
C2 = GB_mex_apply (Cin, Mask, 'plus', op, A, []) ;
GB_spec_compare (C1, C2) ;
Cmask = spones (GB_mex_cast (full (Cin.matrix), Cin.class)) ;
% with C == mask, and outp = replace
C1 = GB_spec_apply (Cin, Cmask, [], op, A, dr) ;
C2 = GB_mex_apply_maskalias (Cin, [], op, A, dr) ;
GB_spec_compare (C1, C2) ;
% with C == mask and accum, and outp = replace
C1 = GB_spec_apply (Cin, Cmask, 'plus', op, A, dr) ;
C2 = GB_mex_apply_maskalias (Cin, 'plus', op, A, dr) ;
GB_spec_compare (C1, C2) ;
% no mask, transpose
C1 = GB_spec_apply (Cin, [], [], op, B, dt) ;
C2 = GB_mex_apply (Cin, [], [], op, B, dt) ;
GB_spec_compare (C1, C2) ;
% no mask, with accum, transpose
C1 = GB_spec_apply (Cin, [], 'plus', op, B, dt) ;
C2 = GB_mex_apply (Cin, [], 'plus', op, B, dt) ;
GB_spec_compare (C1, C2) ;
% with mask, transpose
C1 = GB_spec_apply (Cin, Mask, [], op, B, dt) ;
C2 = GB_mex_apply (Cin, Mask, [], op, B, dt) ;
GB_spec_compare (C1, C2) ;
% with mask and accum, transpose
C1 = GB_spec_apply (Cin, Mask, 'plus', op, B, dt) ;
C2 = GB_mex_apply (Cin, Mask, 'plus', op, B, dt) ;
GB_spec_compare (C1, C2) ;
end
end
end
end
end
end
fprintf ('\n') ;
fprintf ('\ntest130: all tests passed\n') ;
|
