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
|
function test179
%TEST179 bitmap select
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2022, All Rights Reserved.
% SPDX-License-Identifier: Apache-2.0
fprintf ('test179: --------------------------------- bitmap select\n') ;
ops = { 'nonzero', 'eq_zero', 'ne_thunk', 'eq_thunk' } ;
n = 20 ;
rng ('default') ;
Cin = GB_spec_random (n, n, 0.5, 1, 'double complex') ;
A = GB_spec_random (n, n, 0.5, 1, 'double complex') ;
scalar = 3+1i ;
A.matrix (2:3,2:4) = scalar ;
A.pattern (2:3,2:4) = true ;
for builtin = 0:1
GB_builtin_complex_set (builtin) ;
for sparsity_control = [1 2 4]
Cin.sparsity = sparsity_control ;
A.sparsity = sparsity_control ;
for k = 1:length (ops)
op = ops {k} ;
C1 = GB_spec_select (Cin, [], [], op, A, scalar, []) ;
C2 = GB_mex_select (Cin, [], [], op, A, scalar, []) ;
GB_spec_compare (C1, C2) ;
end
end
end
% try a user-defined selectop
Cin = GB_spec_random (n, n, 0.5, 1, 'double') ;
A = GB_spec_random (n, n, 0.5, 1, 'double') ;
A.matrix (2:3,2:4) = nan ;
A.pattern (2:3,2:4) = true ;
for sparsity_control = [1 2 4]
Cin.sparsity = sparsity_control ;
A.sparsity = sparsity_control ;
C1 = GB_spec_select (Cin, [], [], 'isnan', A, 0, []) ;
C2 = GB_mex_select (Cin, [], [], 'isnan', A) ;
GB_spec_compare (C1, C2) ;
end
GB_builtin_complex_set (true) ;
fprintf ('\ntest179: all tests passed\n') ;
|
