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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
|
function gbtest31
%GBTEST31 test GrB and casting
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2022, All Rights Reserved.
% SPDX-License-Identifier: Apache-2.0
types = gbtest_types ;
fprintf ('gbtest31: typecasting\n') ;
rng ('default') ;
for k = 1:length (types)
type = types {k} ;
fprintf ('%s ', type) ;
for m = 0:5
for n = 0:5
A = gbtest_cast (zeros (m, n), type) ;
G = GrB (m, n, type) ;
C = GrB (G, type) ;
assert (gbtest_eq (A, C)) ;
end
end
A = 100 * rand (5, 5) ;
A = gbtest_cast (A, type) ;
G = GrB (A) ;
G2 = sparse (G) ;
assert (gbtest_eq (G, G2)) ;
for k2 = 1:length (types)
type2 = types {k} ;
G2 = GrB (G, type2) ;
A2 = gbtest_cast (A, type2) ;
C = GrB (G2, type2) ;
assert (gbtest_eq (A2, C)) ;
end
F = 100 * ones (5, 5) ;
F = gbtest_cast (F, type) ;
id = F (1,1) ;
A = 100 * sparse (diag (1:5)) ;
G = GrB (A, type) ;
G2 = GrB (F) ;
G2 (logical (speye (5))) = 100:100:500 ;
for k2 = 1:length (types)
type2 = types {k} ;
G3 = full (G, type2, id) ;
H5 = full (G, type2, GrB (id)) ;
G4 = GrB (G2, type2) ;
assert (gbtest_eq (G3, G4)) ;
assert (gbtest_eq (G3, H5)) ;
assert (gbtest_eq (double (G3), double (G4))) ;
assert (gbtest_eq (single (G3), single (G4))) ;
assert (gbtest_eq (uint16 (G3), uint16 (G4))) ;
end
end
fprintf ('\ngbtest31: all tests passed\n') ;
|
