fix: msbfs parents, clarray.copyto bugs

Author: artemiipatov

src/GraphBLAS-sharp.Backend/Algorithms/MSBFS.fs

@@ -119,6 +119,8 @@ module internal MSBFS =
                 let vertexCount = matrix.RowCount
                 let sourceVertexCount = source.Length
 
+                let source = source |> List.sort
+
                 let startMatrix =
                     source |> List.mapi (fun i vertex -> i, vertex, 1)
 
@@ -185,8 +187,6 @@ module internal MSBFS =
 
             let copyIndices = ClArray.copyTo clContext workGroupSize
 
-            // let copyMatrix = Matrix.copy clContext workGroupSize
-
             fun (queue: MailboxProcessor<Msg>) allocationMode (front: ClMatrix.COO<_>) (parents: ClMatrix.COO<_>) ->
 
                 // Find intersection of levels and front indices.
@@ -214,8 +214,8 @@ module internal MSBFS =
 
             let spGeMM =
                 Operations.SpGeMM.COO.expand
-                    (ArithmeticOperations.min -1)
-                    (ArithmeticOperations.fst -1)
+                    (ArithmeticOperations.min)
+                    (ArithmeticOperations.fst)
                     clContext
                     workGroupSize
 
@@ -226,6 +226,8 @@ module internal MSBFS =
                 let vertexCount = inputMatrix.RowCount
                 let sourceVertexCount = source.Length
 
+                let source = source |> List.sort
+
                 let matrix =
                     match inputMatrix with
                     | ClMatrix.CSR m ->

src/GraphBLAS-sharp.Backend/Common/ClArray.fs

@@ -133,7 +133,7 @@ module ClArray =
                 let i = ndRange.GlobalID0
 
                 if i < inputArrayLength then
-                    source.[i] <- destination.[i] @>
+                    destination.[i] <- source.[i] @>
 
         let program = clContext.Compile(copy)
 

src/GraphBLAS-sharp.Backend/Quotes/Arithmetic.fs

@@ -255,17 +255,9 @@ module ArithmeticOperations =
             | None, Some y -> Some y
             | _ -> None @>
 
-    let min zero =
-        <@ fun x y ->
-            let result = min x y
+    let min<'a when 'a: comparison> = <@ fun (x: 'a) (y: 'a) -> Some (min x y) @>
 
-            if result = zero then
-                None
-            else
-                Some result @>
-
-    let fst zero =
-        <@ fun x _ -> if x = zero then None else Some x @>
+    let fst<'a> = <@ fun (x: 'a) (_: 'a) -> Some x @>
 
     //PageRank specific
     let squareOfDifference =