Add reader for int matrices

Author: dpanfilyonok

src/GraphBLAS-sharp/IO/MtxReader.fs

@@ -159,10 +159,68 @@ type MtxReader(pathToFile: string) =
         | Coordinate -> matrixFromCoordinateFormat ()
         | Array -> failwith "Unsupported matrix format"
 
+    member this.ReadMatrixInteger(converter: string -> 'a) : Matrix<'a> =
+        if object <> MtxMatrix then failwith "Object is not matrix"
+        if field <> Integer then failwith "Field is not real"
+
+        use streamReader = new StreamReader(pathToFile)
+        while streamReader.Peek() = int '%' do
+            streamReader.ReadLine() |> ignore
+
+        let matrixFromCoordinateFormat () =
+            let size =
+                streamReader.ReadLine().Split(' ')
+                |> Array.map int
+
+            let n = size.[0]
+            let m = size.[1]
+            let nnz = size.[2]
+
+            let pack x y = (uint64 x <<< 32) ||| (uint64 y)
+            let unpack x = int ((x &&& 0xFFFFFFFF0000000UL) >>> 32), int (x &&& 0xFFFFFFFUL)
+
+            let sortedData =
+                [0 .. nnz - 1]
+                |> List.map (fun _ -> streamReader.ReadLine().Split(' '))
+                |> Array.ofList
+                |> Array.Parallel.map
+                    (fun line ->
+                        let i = int line.[0]
+                        let j = int line.[1]
+                        let v = converter line.[2]
+                        struct(pack i j, v)
+                    )
+                |> Array.sortBy (fun struct(packedIndex, _) -> packedIndex)
+
+            let rows = Array.zeroCreate sortedData.Length
+            let cols = Array.zeroCreate sortedData.Length
+            let values = Array.zeroCreate sortedData.Length
+
+            Array.Parallel.iteri (fun i struct(packedIndex, value) ->
+                let (rowIdx, columnIdx) = unpack packedIndex
+                // in mtx indecies start at 1
+                rows.[i] <- rowIdx - 1
+                cols.[i] <- columnIdx - 1
+                values.[i] <- value
+            ) sortedData
+
+            MatrixCOO {
+                Rows = rows
+                Columns = cols
+                Values = values
+                RowCount = n
+                ColumnCount = m
+            }
+
+        match format with
+        | Coordinate -> matrixFromCoordinateFormat ()
+        | Array -> failwith "Unsupported matrix format"
+
     member this.ReadMatrix(converter: string -> 'a) : Matrix<'a> =
         match field with
         | Real -> this.ReadMatrixReal(converter)
         | Pattern -> this.ReadMatrixBoolean(converter)
+        | Integer -> this.ReadMatrixInteger(converter)
         | _ -> failwith "Not implemented"
 
 and MtxObject =

tests/GraphBLAS-sharp.Tests/MatrixOperationsTests/EWiseAddTests.fs

@@ -144,7 +144,7 @@ let testFixtures case = [
 
     case
     |> correctnessGenericTest<bool> AnyAll.bool (=)
-    |> testPropertyWithConfig config (getCorrectnessTestName "bool")
+    |> ptestPropertyWithConfig config (getCorrectnessTestName "bool")
 ]
 
 let tests =

tests/GraphBLAS-sharp.Tests/MatrixOperationsTests/GetTuplesTests.fs

@@ -126,7 +126,7 @@ let testFixtures case = [
 
     case
     |> correctnessGenericTest<bool> (=) false
-    |> testPropertyWithConfig config (getCorrectnessTestName "bool")
+    |> ptestPropertyWithConfig config (getCorrectnessTestName "bool")
 ]
 
 let tests =

tests/GraphBLAS-sharp.Tests/MatrixOperationsTests/MxmTests.fs

@@ -159,7 +159,7 @@ let testFixtures case = [
 
     case
     |> correctnessGenericTest<bool> AnyAll.bool (=)
-    |> testPropertyWithConfig config (getCorrectnessTestName "bool")
+    |> ptestPropertyWithConfig config (getCorrectnessTestName "bool")
 ]
 
 let tests =

tests/GraphBLAS-sharp.Tests/MatrixOperationsTests/MxvTests.fs

@@ -179,7 +179,7 @@ let testFixtures case = [
 
     case
     |> correctnessGenericTest<bool> AnyAll.bool (=)
-    |> testPropertyWithConfig config (getCorrectnessTestName "bool")
+    |> ptestPropertyWithConfig config (getCorrectnessTestName "bool")
 
     testCase (sprintf "Explicit zero test on %A" case) <| fun () ->
         let matrix = array2D [

tests/GraphBLAS-sharp.Tests/MatrixOperationsTests/TransposeTests.fs

@@ -131,7 +131,7 @@ let testFixtures case = [
 
     case
     |> correctnessGenericTest<bool> (=) false
-    |> testPropertyWithConfig config (getCorrectnessTestName "bool")
+    |> ptestPropertyWithConfig config (getCorrectnessTestName "bool")
 ]
 
 let tests =

tests/GraphBLAS-sharp.Tests/MatrixOperationsTests/VxmTests.fs

@@ -164,7 +164,7 @@ let testFixtures case = [
 
     case
     |> correctnessGenericTest<bool> AnyAll.bool (=)
-    |> testPropertyWithConfig config (getCorrectnessTestName "bool")
+    |> ptestPropertyWithConfig config (getCorrectnessTestName "bool")
 ]
 
 let tests =