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
64
65
66
67
68
69
70
71
72
73
74
|
function gbtest38
%GBTEST38 test sqrt, eps, ceil, floor, round, fix, real, conj, ...
% isfinite, isinf, isnan, spfun, eig
% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2022, All Rights Reserved.
% SPDX-License-Identifier: Apache-2.0
rng ('default') ;
for trial = 1:40
A = 1e3 * rand (3) ;
B = single (A) ;
G = GrB (A) ;
H = GrB (B) ;
err = norm (sqrt (A) - sqrt (G), 1) ; assert (err < 8 * eps ('double')) ;
err = norm (sqrt (B) - sqrt (H), 1) ; assert (err < 8 * eps ('single')) ;
assert (gbtest_eq (eps (A), eps (G))) ;
assert (gbtest_eq (eps (B), eps (H))) ;
assert (gbtest_eq (ceil (A), ceil (G))) ;
assert (gbtest_eq (ceil (B), ceil (H))) ;
assert (gbtest_eq (floor (A), floor (G))) ;
assert (gbtest_eq (floor (B), floor (H))) ;
assert (gbtest_eq (round (A), round (G))) ;
assert (gbtest_eq (round (B), round (H))) ;
assert (gbtest_eq (fix (A), fix (G))) ;
assert (gbtest_eq (fix (B), fix (H))) ;
assert (gbtest_eq (real (A), real (G))) ;
assert (gbtest_eq (real (B), real (H))) ;
assert (gbtest_eq (conj (A), conj (G))) ;
assert (gbtest_eq (conj (B), conj (H))) ;
C = A ;
C (1,1) = inf ;
C (2,2) = nan ;
G = GrB (C) ;
assert (gbtest_eq (isfinite (C), isfinite (G))) ;
assert (gbtest_eq (isnan (C), isnan (G))) ;
A = sprand (10, 10, 0.5) ;
G = GrB (A) ;
assert (gbtest_eq (spfun (@exp, A), double (spfun (@exp, G)))) ;
A = rand (10) ;
G = GrB (A) ;
assert (isequal (eig (A), double (eig (G)))) ;
A = sparse (A+A') ;
G = GrB (A) ;
err = norm (eig (A) - double (eig (G)), 2) ;
assert (err < 1e-12) ;
A = rand (10) ;
B = rand (10) ;
G = GrB (A) ;
H = GrB (B) ;
[V1,D1] = eig (A, B) ;
[V2,D2] = eig (G, H) ;
assert (isequal (V1, double (V2))) ;
assert (isequal (D1, double (D2))) ;
end
fprintf ('gbtest38: all tests passed\n') ;
|
