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
|
% run this output of this program as a script:
% GraphBLAS matrix A: nrows: 3 ncols 5 entries: 10
A = sparse (3,5) ;
A (1,2) = ( 0) + ( 0.8253972948485638)*1i ;
A (1,3) = ( 0.8158771738902956) + ( 0.6199473944886393)*1i ;
A (1,4) = ( 0.230335927869842) + ( 0)*1i ;
A (1,5) = ( 0.6834455780600209) + ( 0)*1i ;
A (2,1) = ( 0) + ( 0.8145545554415048)*1i ;
A (2,2) = ( 0) + ( 0.8106511441383925)*1i ;
A (2,3) = ( 0) + ( 0.6296551667948371)*1i ;
A (2,4) = ( 0.6749993971837149) + ( 0.1028046574093258)*1i ;
A (2,5) = ( 0.1781985185630041) + ( 0)*1i ;
A (3,1) = ( 0.6668132827629559) + ( 0)*1i ;
A
% GraphBLAS matrix B: nrows: 5 ncols 4 entries: 13
B = sparse (5,4) ;
B (1,1) = ( 0.03470661760412896) + ( 0)*1i ;
B (1,2) = ( 0.2320846810101108) + ( 0)*1i ;
B (1,4) = ( 0.4037695773165162) + ( 0)*1i ;
B (2,1) = ( 0) + ( 0.70501889367129)*1i ;
B (2,3) = ( 0.2114568817022959) + ( 0)*1i ;
B (2,4) = ( 0) + ( 0.4423346124115619)*1i ;
B (3,1) = ( 0.929676618909714) + ( 0)*1i ;
B (3,2) = ( 0.03154825424997802) + ( 0)*1i ;
B (3,3) = ( 0) + ( 0.1598994553484671)*1i ;
B (3,4) = ( 0.9813322475732387) + ( 0.9048904118202962)*1i ;
B (4,1) = ( 0) + ( 0.9376008381429164)*1i ;
B (4,4) = ( 0) + ( 0.8269171297987302)*1i ;
B (5,4) = ( 0.1371733231185191) + ( 0)*1i ;
B
% GraphBLAS matrix C: nrows: 3 ncols 4 entries: 11
C = sparse (3,4) ;
C (1,1) = ( 0.1765812448145327) + ( 0.7923137566352751)*1i ;
C (1,2) = ( 0.02573950051864457) + ( 0.01955825802293901)*1i ;
C (1,3) = (-0.09912925072343469) + ( 0.3049942538704926)*1i ;
C (1,4) = (-0.03168916369751931) + ( 1.537122526250736)*1i ;
C (2,1) = ( -0.6679141057457941) + ( 1.24652612056368)*1i ;
C (2,2) = ( 0) + ( 0.2089101554568301)*1i ;
C (2,3) = ( -0.1006815182278427) + ( 0.1714177630879029)*1i ;
C (2,4) = ( -0.9889148320979797) + ( 1.50496183271377)*1i ;
C (3,1) = ( 0.02314283361820783) + ( 0)*1i ;
C (3,2) = ( 0.1547571480233455) + ( 0)*1i ;
C (3,4) = ( 0.2692389173302374) + ( 0)*1i ;
C
E = A*B
err = norm (C-E,1)
assert (err < 1e-12)
|
