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
|
function test92
%TEST92 test GB_subref (symbolic case)
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2020, All Rights Reserved.
% http://suitesparse.com See GraphBLAS/Doc/License.txt for license.
fprintf ('test92: test GB_subref (symbolic case)\n') ;
rng ('default') ;
n = 100 ;
A = sprand (n, n, 0.1) ;
% S = A (I,J) ;
% No column of A can be dense
% length(I) > 1
% I is not ":", not lo:hi
% sort not needed so I must be sorted
% 32*length(I) > entries in columns of A
I = [ 1 2 2 4:8 20:30 32 32 32 32] ;
I0 = uint64 (I) - 1 ;
J = 1:n ;
J0 = uint64 (J) - 1 ;
nI = length (I) ;
nJ = n ;
S = GB_mex_subref_symbolic (A, I0, J0) ;
[si,sj,sx] = find (S.matrix) ;
[i,j,x] = find (A) ;
nz = length (x) ;
x = 1:nz ;
X = sparse (i,j,x,n,n) ;
S2 = X (I,J) ;
[si2,sj2,sx2] = find (S2) ;
assert (isequal (si, si2)) ;
assert (isequal (sj, sj2)) ;
assert (isequal (sx+1, sx2)) ;
fprintf ('test92: all tests passed\n') ;
|
