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
|
function testcb
%TESTCB test complex reduce
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2020, All Rights Reserved.
% http://suitesparse.com See GraphBLAS/Doc/License.txt for license.
tol = 1e-13 ;
seed = 1 ;
for m = [1 5 10 100]
for n = [1 5 10 100]
for trial = 1:100
A = GB_mex_random (m, n, 10*(m+n), 1, seed) ; seed = seed + 1 ;
s = complex (0) ;
c1 = complex (full (sum (A (:)))) ;
c2 = GB_mex_reduce_to_scalar (s, [], 'plus', A) ;
assert (abs (c1 - c2) <= tol * abs (c1) + tol) ;
% assert (isequal (c1,c2))
[i j x] = find (A) ;
c1 = complex (prod (x)) ;
c2 = GB_mex_reduce_to_scalar (s, [], 'times', A) ;
assert (abs (c1 - c2) <= tol * abs (c1) + tol) ;
% assert (isequal (c1,c2))
s = complex (pi,2) ;
c1 = s * complex (full (sum (A (:)))) ;
c2 = GB_mex_reduce_to_scalar (s, 'times', 'plus', A) ;
assert (abs (c1 - c2) <= tol * abs (c1) + tol) ;
% assert (isequal (c1,c2))
s = complex (pi,2) ;
c1 = s + complex (full (sum (A (:)))) ;
c2 = GB_mex_reduce_to_scalar (s, 'plus', 'plus', A) ;
assert (abs (c1 - c2) <= tol * abs (c1) + tol) ;
% assert (isequal (c1,c2))
if (n > 1)
w = GB_mex_complex (sparse (m, 1)) ;
c1 = sum (A.').' ;
c2 = GB_mex_reduce_to_vector (w, [], [], 'plus', A, []) ;
assert (norm (c1 - c2.matrix, 1) <= tol * norm (c1, 1) + tol) ;
% assert (isequal (c1,c2.matrix))
w = GB_mex_random (m, 1, 10, 1, seed) ; seed = seed + 1 ;
c1 = w + sum (A.').' ;
c2 = GB_mex_reduce_to_vector (w, [], 'plus', 'plus', A, []) ;
assert (norm (c1 - c2.matrix, 1) <= tol * norm (c1, 1) + tol) ;
% assert (isequal (c1,c2.matrix))
end
end
end
end
fprintf ('testcb: all complex reduce tests passed\n') ;
|
