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 test96
%TEST96 test dot product
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2022, All Rights Reserved.
% SPDX-License-Identifier: Apache-2.0
n = 1000 ;
A = sprandn (n, n, 0.5) ;
B = sprandn (n, n, 0.5) ;
C1 = A'*B ;
M = spones (sprandn (n, n, 0.5)) ;
C2 = GB_mex_AdotB (A,B) ;
err = norm (C1-C2, 1) / norm (C1, 1)
assert (err < 1e-12)
C4 = C1 .* M ;
C3 = GB_mex_AdotB (A,B,M) ;
err = norm (C3-C4, 1) / norm (C3, 1)
assert (err < 1e-12)
fprintf ('test96: all tests passed\n') ;
|