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
|
function test147
%TEST147 test C<M>A*B with very sparse M
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2022, All Rights Reserved.
% SPDX-License-Identifier: Apache-2.0
fprintf ('test147 ----------------------------- C<M>A*B with very sparse M\n') ;
rng ('default') ;
n = 1000 ;
A = sprand (n, n, 0.01) ;
A (:,1:300) = 1 ;
S = sparse (n, n) ;
M = spones (speye (n) + sparse (n, n, 1e-5)) ;
semiring.add = 'plus' ;
semiring.multiply = 'times' ;
semiring.class = 'double' ;
% also create the /tmp/GB_mxm_[code].h file, flipxy = false
C1 = GB_mex_mxm (S, M, [ ], semiring, A, A, [ ], 0) ;
C2 = (A*A) .* M ;
assert (norm (C1.matrix - C2, 1) < 1e-12)
fprintf ('test147: all tests passed\n') ;
|
