add method

Author: Wimmerer

benchmarks/benchmarks.jl

@@ -206,14 +206,50 @@ function CaccumAxB_allbyrow(S, G, nthreads, sizerhs)
     end
 end
 
+function CaccumAxB_CRC(S, G, nthreads, sizerhs)
+    printstyled("\nFull_bycol += CSR * Full_bycol\n", color=:green)
+    GC.gc()
+    m = rand(size(S, 2), sizerhs)
+    m2 = GBMatrix(m)
+    println("Size of dense matrix is $(size(m))")
+
+    gbset(m2, :format, :bycol)
+    gbset(G, :format, :byrow)
+
+    printstyled("\nSparseMatrixCSC:\n", bold=true)
+    A = @benchmark $S * $m
+    show(stdout, MIME("text/plain"), A)
+    gbset(:burble, true)
+    C = GBMatrix(size(G, 1), size(m2, 2), 0.0)
+    gbset(C, :sparsity_control, :full)
+    gbset(C, :format, :bycol)
+
+    printstyled("\nGBMatrix:\n", bold=true)
+    #print burble for checking
+    gbset(:burble, false)
+    for n ∈ nthreads
+        printstyled("\nF += S' * F with $n threads: \n", bold=true, color=:green)
+        gbset(:nthreads, n)
+        gbset(:burble, true)
+        mul!(C, G, m2; accum=+)
+        gbset(:burble, false)
+        B = @benchmark mul!($C, $G, $m2; accum=+)
+        show(stdout, MIME("text/plain"), B)
+        tratio = ratio(median(A), median(B))
+        color = tratio.time >= 1.0 ? :green : :red
+        printstyled("\nMedian speedup over SparseArrays using $n threads is: $(string(tratio))\n"; bold=true, color)
+    end
+end
+
 # OPTIONS SET 2:
 # run these functions for benchmarking:
-const functorun = [AxB_allbycol, AxB_ColxRow, CaccumAxB_allbycol, CaccumAxB_allbyrow]
+const functorun = [AxB_allbycol, AxB_ColxRow, CaccumAxB_allbycol, CaccumAxB_allbyrow, CaccumAxB_CRC]
 #= The choices are:
 AxB_allbycol - S * F
 AxB_ColxRow - S' * F
 CaccumAxB_allbycol - F += S * F
 CaccumAxB_allbyrow - F' += S' * F'
+CaccumAxB_CRC - F += S' * F
 Please open an issue or message me for further functions to add here.
 =#