tpose to bench

Author: Will Kimmerer

benchmarks/benchmarks2.jl

@@ -47,6 +47,7 @@ function mxm(A::SparseMatrixCSC, B)
     printstyled(stdout, "\nC = A::SparseMatrixCSC($(size(A))) * B::$(typeof(B))($(size(B)))\n")
     result = @bench A * B
     println(stdout, result, "s")
+    GC.gc()
     flush(stdout)
     return result
 end
@@ -65,10 +66,34 @@ function mxm(A::SuiteSparseGraphBLAS.GBArray, B::SuiteSparseGraphBLAS.GBArray; a
         result = @gbbench mul!(C, A, B; accum=+)
     end
     println(stdout, result, "s")
+    GC.gc()
+    flush(stdout)
+    return result
+end
+
+function tpose(A::SuiteSparseGraphBLAS.GBArray)
+    Ao = storageorder(A) == ColMajor() ? "C" : "R"
+    Bo = storageorder(B) == ColMajor() ? "C" : "R"
+    printstyled(stdout, "\nC::GBArray = transpose(A::GBArray($Ao, $(size(A))))\n")
+    result = @gbbench copy(transpose(A))
+    println(stdout, result, "s")
+    GC.gc()
     flush(stdout)
     return result
 end
 
+function tpose(A::SparseMatrixCSC)
+    printstyled(stdout, "\nC = transpose(A::SparseMatrixCSC($(size(A))))\n")
+    result = @bench copy(transpose(A))
+    println(stdout, result, "s")
+    GC.gc()
+    flush(stdout)
+    return result
+end
+
+function spdbench(A)
+end
+
 function singlebench(pathornum)
     x = tryparse(Int64, pathornum)
     if x !== nothing
@@ -160,6 +185,26 @@ function singlebench(pathornum)
     println(stdout, "A by col (1, 2, 16 thread): $gbresultsC")
     println(stdout, "SparseArrays: $SAresults")
     flush(stdout)
+
+    printstyled(stdout, "\nSparse * Sparse'"; bold=true)
+    println(stdout, "################################")
+    flush(stdout)
+    gbset(A, :format, SuiteSparseGraphBLAS.BYROW)
+    diag(A)
+    gbresultsR = runthreaded(A, transpose(A))
+    gbset(A, :format, SuiteSparseGraphBLAS.BYCOL)
+    diag(A)
+    gbresultsC = runthreaded(A, transpose(A))
+    A2 = SparseMatrixCSC(A) 
+    SAresults = mxm(A2, transpose(A2))
+    println(stdout, )
+    printstyled(stdout, "\nRESULTS, Sparse * Sparse: \n"; bold=true, color=:green)
+    println(stdout, "################################")
+    println(stdout, "A by row (1, 2, 16 thread): $gbresultsR")
+    println(stdout, "A by col (1, 2, 16 thread): $gbresultsC")
+    println(stdout, "SparseArrays: $SAresults")
+    flush(stdout)
+
     return nothing
 end
 

benchmarks/mtx.txt

@@ -1,22 +1,23 @@
-# You may provide paths or SuiteSparse ID numbers here.
-# A hashtag denotes a commented line.
-# /home/faculty/d/davis/Freescale2.mtx
-2662
-# /home/faculty/d/davis/indochina-2004.mtx
-2451
-# /home/faculty/d/davis/nd24k.mtx
-939
-# /home/faculty/d/davis/nd6k.mtx
-937
-# /home/faculty/d/davis/roadNet-CA.mtx
-2317
-# /raid/GAP/GAP-twitter/GAP-twitter.mtx
-2852
-# /raid/GAP/GAP-road/GAP-road.mtx
-2854
-# /raid/GAP/GAP-web/GAP-web.mtx
-# 2853
-# /raid/GAP/GAP-urand/GAP-urand.mtx
-# 2856
-# /raid/GAP/GAP-kron/GAP-kron.mtx
-# 2855
+# # You may provide paths or SuiteSparse ID numbers here.
+# # A hashtag denotes a commented line.
+# # /home/faculty/d/davis/Freescale2.mtx
+# 2662
+# # /home/faculty/d/davis/indochina-2004.mtx
+# 2451
+# # /home/faculty/d/davis/nd24k.mtx
+# 939
+# # /home/faculty/d/davis/nd6k.mtx
+# 937
+# # /home/faculty/d/davis/roadNet-CA.mtx
+# 2317
+# # /raid/GAP/GAP-twitter/GAP-twitter.mtx
+# 2852
+# # /raid/GAP/GAP-road/GAP-road.mtx
+# 2854
+# # /raid/GAP/GAP-web/GAP-web.mtx
+# # 2853
+# # /raid/GAP/GAP-urand/GAP-urand.mtx
+# # 2856
+# # /raid/GAP/GAP-kron/GAP-kron.mtx
+# # 2855
+# C:\Users\wrkim\Freescale2\Freescale2.mtx