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
|
function test201
%TEST201 test iso reduce to vector and reduce to scalar
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2025, All Rights Reserved.
% SPDX-License-Identifier: Apache-2.0
n = 10 ;
A.matrix = pi * sparse (ones (n)) ;
A.class = 'double' ;
A.iso = true ;
win = sparse (n,1) ;
w1 = GB_mex_reduce_to_vector (win, [ ], [ ], 'plus', A, [ ]) ;
w2 = sum (A.matrix, 2) ;
err = norm (w1.matrix - w2, 1) ;
assert (err < 1e-12)
c0 = 0 ;
c2 = GB_spec_reduce_to_scalar (c0, [ ], 'plus', A) ;
c1 = GB_mex_reduce_to_scalar (c0, [ ], 'plus', A) ;
assert (abs (c1-c2) < 1e-12) ;
c3 = sum (sum (A.matrix)) ;
assert (abs (c2-c3) < 1e-12) ;
% reduce a huge iso full matrix to a scalar
m = 2^40 ;
n = 2^48 ;
s = GB_mex_test36 ;
% expected result:
t = pi * m * n ;
relerr = abs (s - t) / t ;
relerr
assert (relerr < 1e-14) ;
bits = ceil (log2 (t))
fprintf ('test201: all tests passed\n') ;
|
