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
|
function testc5
%TESTC5 test complex subref
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2020, All Rights Reserved.
% http://suitesparse.com See GraphBLAS/Doc/License.txt for license.
seed = 1 ;
for m = [1 5 10 100]
for n = [1 5 10 100]
seed = seed + 1 ;
A = GB_mex_random (m, n, 10*(m+n), 1, seed) ;
for trials = 1:10
J = randperm (n, 1+floor(n/2)) ;
I = randperm (m, 1+floor(m/2)) ;
J0 = uint64 (J-1) ;
I0 = uint64 (I-1) ;
C1 = GB_mex_Matrix_subref (A, I0, J0) ;
C2 = GB_mex_complex (A (I,J)) ;
assert (isequal (C1, C2))
end
end
end
fprintf ('testc5: all complex subref C = A(I,J) tests passed\n') ;
|