1
2% run this output of this program as a script in MATLAB:
3
4% GraphBLAS matrix A: nrows: 3 ncols 5 entries: 10
5A = sparse (3,5) ;
6    A (1,2) =  (                   0) + (  0.8253972948485638)*1i ;
7    A (1,3) =  (  0.8158771738902956) + (  0.6199473944886393)*1i ;
8    A (1,4) =  (   0.230335927869842) + (                   0)*1i ;
9    A (1,5) =  (  0.6834455780600209) + (                   0)*1i ;
10    A (2,1) =  (                   0) + (  0.8145545554415048)*1i ;
11    A (2,2) =  (                   0) + (  0.8106511441383925)*1i ;
12    A (2,3) =  (                   0) + (  0.6296551667948371)*1i ;
13    A (2,4) =  (  0.6749993971837149) + (  0.1028046574093258)*1i ;
14    A (2,5) =  (  0.1781985185630041) + (                   0)*1i ;
15    A (3,1) =  (  0.6668132827629559) + (                   0)*1i ;
16A
17
18% GraphBLAS matrix B: nrows: 5 ncols 4 entries: 13
19B = sparse (5,4) ;
20    B (1,1) =  ( 0.03470661760412896) + (                   0)*1i ;
21    B (1,2) =  (  0.2320846810101108) + (                   0)*1i ;
22    B (1,4) =  (  0.4037695773165162) + (                   0)*1i ;
23    B (2,1) =  (                   0) + (    0.70501889367129)*1i ;
24    B (2,3) =  (  0.2114568817022959) + (                   0)*1i ;
25    B (2,4) =  (                   0) + (  0.4423346124115619)*1i ;
26    B (3,1) =  (   0.929676618909714) + (                   0)*1i ;
27    B (3,2) =  ( 0.03154825424997802) + (                   0)*1i ;
28    B (3,3) =  (                   0) + (  0.1598994553484671)*1i ;
29    B (3,4) =  (  0.9813322475732387) + (  0.9048904118202962)*1i ;
30    B (4,1) =  (                   0) + (  0.9376008381429164)*1i ;
31    B (4,4) =  (                   0) + (  0.8269171297987302)*1i ;
32    B (5,4) =  (  0.1371733231185191) + (                   0)*1i ;
33B
34
35% GraphBLAS matrix C: nrows: 3 ncols 4 entries: 11
36C = sparse (3,4) ;
37    C (1,1) =  (  0.1765812448145327) + (  0.7923137566352751)*1i ;
38    C (1,2) =  ( 0.02573950051864457) + ( 0.01955825802293901)*1i ;
39    C (1,3) =  (-0.09912925072343469) + (  0.3049942538704926)*1i ;
40    C (1,4) =  (-0.03168916369751931) + (   1.537122526250736)*1i ;
41    C (2,1) =  ( -0.6679141057457941) + (    1.24652612056368)*1i ;
42    C (2,2) =  (                   0) + (  0.2089101554568301)*1i ;
43    C (2,3) =  ( -0.1006815182278427) + (  0.1714177630879029)*1i ;
44    C (2,4) =  ( -0.9889148320979797) + (    1.50496183271377)*1i ;
45    C (3,1) =  ( 0.02314283361820783) + (                   0)*1i ;
46    C (3,2) =  (  0.1547571480233455) + (                   0)*1i ;
47    C (3,4) =  (  0.2692389173302374) + (                   0)*1i ;
48C
49E = A*B
50err = norm (C-E,1)
51assert (err < 1e-12)
52