Add quickgraph benchmark

Author: dpanfilyonok

.vscode/settings.json

@@ -4,5 +4,5 @@
     "FSharp.analyzersPath": [
         "./packages/analyzers"
     ],
-    "FSharp.suggestSdkScripts": false
+    "FSharp.suggestSdkScripts": false,
 }

.vscode/tasks.json

@@ -21,4 +21,4 @@
             "problemMatcher": "$msCompile"
         }
     ]
-}
+}

benchmarks/GraphBLAS-sharp.Benchmarks/BenchmarksQG.fs

@@ -0,0 +1,75 @@
+namespace GraphBLAS.FSharp.Benchmarks
+
+open GraphBLAS.FSharp
+open GraphBLAS.FSharp.Algorithms
+open BenchmarkDotNet.Attributes
+open BenchmarkDotNet.Configs
+open BenchmarkDotNet.Columns
+open System.IO
+open System
+open System.IO
+open System.Text.RegularExpressions
+open GraphBLAS.FSharp
+open GraphBLAS.FSharp.Predefined
+open BenchmarkDotNet.Attributes
+open BenchmarkDotNet.Configs
+open BenchmarkDotNet.Columns
+open BenchmarkDotNet.Filters
+open Brahma.FSharp.OpenCL.WorkflowBuilder.Basic
+open Brahma.FSharp.OpenCL.WorkflowBuilder.Evaluation
+open OpenCL.Net
+open GraphBLAS.FSharp.IO
+open QuickGraph
+
+type Config2() =
+    inherit ManualConfig()
+
+    do
+        base.AddColumn(
+            MatrixShapeColumn("RowCount", fun mtxReader -> mtxReader.ReadMatrixShape().RowCount) :> IColumn,
+            MatrixShapeColumn("ColumnCount", fun mtxReader -> mtxReader.ReadMatrixShape().ColumnCount) :> IColumn,
+            MatrixShapeColumn("NNZ", fun mtxReader -> mtxReader.ReadMatrixShape().Nnz) :> IColumn,
+            TEPSColumn() :> IColumn,
+            StatisticColumn.Min,
+            StatisticColumn.Max
+        ) |> ignore
+
+[<IterationCount(5)>]
+[<WarmupCount(3)>]
+[<Config(typeof<Config2>)>]
+type QGBenchmarks() =
+    let random = Random()
+
+    let graph = AdjacencyGraph<int, Edge<int>>(false)
+    let mutable source = 0
+    let mutable bfs = Unchecked.defaultof<Algorithms.Search.BreadthFirstSearchAlgorithm<int, Edge<int>>>
+
+    [<ParamsSource("InputMatricesProvider")>]
+    member val InputMatrixReader = Unchecked.defaultof<MtxReader> with get, set
+
+    [<GlobalSetup>]
+    member this.BuildMatrixAndSetSource() =
+        let matrix = this.InputMatrixReader.ReadMatrixBoolean(fun _ -> true)
+
+        match matrix with
+        | MatrixCSR csr -> failwith "Not implemented"
+        | MatrixCOO coo ->
+            for i = 0 to coo.Values.Length - 1 do
+                graph.AddVerticesAndEdge(Edge(coo.Rows.[i], coo.Columns.[i])) |> ignore
+
+        source <- random.Next <| Matrix.rowCount matrix
+        bfs <- Algorithms.Search.BreadthFirstSearchAlgorithm(graph)
+
+    [<Benchmark>]
+    member this.LevelBFS() =
+        bfs.Compute(source)
+
+    static member InputMatricesProvider =
+        "BFSBenchmarks.txt"
+        |> Utils.getMatricesFilenames
+        |> Seq.map
+            (fun matrixFilename ->
+                match Path.GetExtension matrixFilename with
+                | ".mtx" -> MtxReader(Utils.getFullPathToMatrix "BFSDatasets" matrixFilename)
+                | _ -> failwith "Unsupported matrix format"
+            )

benchmarks/GraphBLAS-sharp.Benchmarks/GraphBLAS-sharp.Benchmarks.fsproj

@@ -17,7 +17,8 @@
     <Compile Include="Helpers.fs" />
     <Compile Include="BenchmarksBFS.fs" />
     <Compile Include="BenchmarksEWiseAdd.fs" />
+    <Compile Include="BenchmarksQG.fs" />
     <Compile Include="Program.fs" />
   </ItemGroup>
   <Import Project="..\..\.paket\Paket.Restore.targets" />
-</Project>
+</Project>

benchmarks/GraphBLAS-sharp.Benchmarks/Program.fs

@@ -6,7 +6,8 @@ let main argv =
     let benchmarks = BenchmarkSwitcher [|
         // typeof<EWiseAddBenchmarks4Float32>
         // typeof<EWiseAddBenchmarks4Bool>
-        typeof<BFSBenchmarks>
+        // typeof<BFSBenchmarks>
+        typeof<QGBenchmarks>
     |]
 
     benchmarks.Run argv |> ignore

benchmarks/GraphBLAS-sharp.Benchmarks/paket.references

@@ -3,3 +3,4 @@ Microsoft.NETFramework.ReferenceAssemblies
 BenchmarkDotNet
 MathNet.Numerics.FSharp
 MathNet.Numerics.MKL.Win-x64
+YC.QuickGraph

paket.dependencies

@@ -23,6 +23,7 @@ nuget BenchmarkDotNet
 nuget MathNet.Numerics.FSharp
 nuget MathNet.Numerics.MKL.Win-x64 2.5.0
 nuget TypeShape 9.0.0
+nuget YC.QuickGraph 3.7.4
 
 // [ FAKE GROUP ]
 group Build

paket.lock

@@ -65,6 +65,9 @@ NUGET
       Microsoft.NETCore.App (>= 2.0) - restriction: >= netcoreapp2.0
       System.Runtime.Loader (>= 4.3) - restriction: && (>= netcoreapp1.0) (< netcoreapp2.0)
       System.ValueTuple (>= 4.4) - restriction: && (>= netcoreapp1.0) (< netcoreapp2.0)
+    DotNet.Contracts (1.10.20606.1)
+    DotParser (1.1)
+      System.ValueTuple (>= 4.5) - restriction: >= net461
     Expecto (8.13.1)
       FSharp.Core (>= 4.3.4) - restriction: || (>= net461) (>= netstandard2.0)
       Mono.Cecil (>= 0.11) - restriction: || (>= net461) (>= netstandard2.0)
@@ -930,6 +933,10 @@ NUGET
       ExtraConstraints.Fody (>= 1.14) - restriction: >= netstandard2.0
       Microsoft.Build.Framework (>= 16.6) - restriction: >= netstandard2.0
       System.CodeDom (>= 4.7) - restriction: >= netstandard2.0
+    YC.QuickGraph (3.7.4)
+      DotNet.Contracts
+      DotParser (>= 1.0.6)
+      FSharp.Core
     YoloDev.Expecto.TestSdk (0.8)
       Expecto (>= 8.10 < 9.0) - restriction: || (>= net461) (>= netcoreapp2.0)
       FSharp.Core (>= 4.3.4) - restriction: || (>= net461) (>= netcoreapp2.0)

src/GraphBLAS-sharp/IO/MtxReader.fs

@@ -23,6 +23,9 @@ type MtxReader(pathToFile: string) =
     member this.Field = field
     member this.Symmetry = symmetry
 
+    override this.ToString() =
+        Path.GetFileName pathToFile
+
     member this.ReadMatrixShape() =
         use streamReader = new StreamReader(pathToFile)
         while streamReader.Peek() = int '%' do

src/GraphBLAS-sharp/Operations/Matrix.fs

@@ -37,8 +37,9 @@ module Matrix =
         methods
     *)
 
-    let rowCount (matrix: Matrix<'a>) : int = failwith "Not Implemented yet"
-    let columnCount (matrix: Matrix<'a>) : int = failwith "Not Implemented yet"
+    let rowCount (matrix: Matrix<'a>) : int = matrix.RowCount
+    let columnCount (matrix: Matrix<'a>) : int = matrix.ColumnCount
+
     let copy (matrix: Matrix<'a>) : GraphblasEvaluation<Matrix<'a>> = failwith "Not Implemented yet"
     let resize (rowCount: int) (columnCount: int) (matrix: Matrix<'a>) : GraphblasEvaluation<Matrix<'a>> = failwith "Not Implemented yet"