Interfaces are added.

Author: gsvgit

2021-03-23-Note-10-09.xopp

@@ -21,6 +21,7 @@
     <Compile Include="AssemblyInfo.fs" />
     <Compile Include="ImgHelpers.fs" />
     <Compile Include="Streaming.fs" />
+    <Compile Include="Matrices.fs" />
     <Compile Include="Main.fs" />
     <None Include="YC.OpenCL.NET.dll.config">
       <CopyToOutputDirectory>Always</CopyToOutputDirectory>

src/ImgProcessing/ImgProcessing.fsproj

@@ -19,8 +19,7 @@ module Main =
         let final = (System.DateTime.Now - start).TotalMilliseconds / 1000.0
         r,final
 
-    [<EntryPoint>]
-    let main (argv: string array) =
+    let imgPerf () =
         let img1 = "../../../Examples/sara-budhwani-h_P71-B8BPw-unsplash.jpg"
         let img2 = "../../../Examples/nick-gavrilov-F-rvSJl6qI0-unsplash.jpg"
         let inDir = "../../../Examples/input"
@@ -47,9 +46,22 @@ module Main =
             |> filter ImgHelpers.edgesKernel
 *)
         //let r,t = time r
-        //let _do () = ImgPorcessing.Streaming.processAllFiles inDir outDir filters
-        let _do () = ImgPorcessing.Streaming.processAllFilesAsync inDir outDir filters
+        let _do () = ImgPorcessing.Streaming.processAllFiles inDir outDir filters
+        //let _do () = ImgPorcessing.Streaming.processAllFilesAsync inDir outDir filters
         let t = time _do
         printfn "Execution time: %A seconds" t
         //ImgHelpers.save2DByteArrayAsImage r (img2 + "_out_blur_edg_10_gpu_2_br.jpg")
+
+    [<EntryPoint>]
+    let main (argv: string array) =
+        (*let platform, queue = Matrices.getCommandQueue "*Intel*"
+
+        let m1 = Matrices.getRandomMatrix 4096 4096
+        let m2 = Matrices.getRandomMatrix 4096 4096
+
+        let res, time = time (fun () -> Matrices.multiply platform queue 16 m1 m2)
+
+        printfn "Time: %A" time
+*)
+        imgPerf ()
         0

src/ImgProcessing/Main.fs

@@ -0,0 +1,75 @@
+module ImgPorcessing.Matrices
+
+open OpenCL.Net
+open Brahma.OpenCL
+open Brahma.FSharp.OpenCL.Core
+open Brahma.FSharp.OpenCL.Extensions
+
+
+let getCommandQueue platformName =
+    let deviceType = DeviceType.Default
+
+    let provider =
+        try  ComputeProvider.Create(platformName, deviceType)
+        with
+        | ex -> failwith ex.Message
+
+    let mutable commandQueue = new CommandQueue(provider, provider.Devices |> Seq.head)
+
+    provider, commandQueue
+
+let getRandomMatrix n m =
+    let rnd = System.Random()
+    let res = Array2D.zeroCreate n m
+    for i in 0 .. n - 1 do
+        for j in 0 .. m - 1 do
+            res.[i,j] <- rnd.NextDouble()
+    res
+let multiply (provider:ComputeProvider) (commandQueue:CommandQueue) localWorkSize (m1: float[,]) (m2: _[,]) =
+    let m1Rows = m1.GetLength 0
+    let m1Cols = m1.GetLength 1
+
+    let m2Rows = m2.GetLength 0
+    let m2Cols = m2.GetLength 1
+
+    let resRows = m1Rows
+    let resCols = m2Cols
+
+    let mToArr m =
+            [| for x in 0 .. Array2D.length1 m - 1 do
+               yield! [| for y in 0 .. Array2D.length2 m - 1 -> m.[x, y] |]
+            |]
+
+    let m1Arr = mToArr m1
+    let m2Arr = mToArr m2
+
+    let command =
+        <@
+            fun (r:_2D) (m1:array<_>) (m2:array<_>) (res:array<_>) ->
+                let row = r.GlobalID0
+                let col = r.GlobalID1
+                let mutable buf = res.[row * resCols + col]
+                for k in 0 .. m2Rows - 1 do
+                    //res.[row * resCols + col] <- res.[row * resCols + col] + (m1.[row * m1Cols + k] * m2.[k * m2Cols + col])
+                    buf <- buf + (m1.[row * m1Cols + k] * m2.[k * m2Cols + col])
+                res.[row * resCols + col] <- buf
+        @>
+
+    let kernel, kernelPrepare, kernelRun =
+        provider.Compile command
+
+    let d = _2D(resRows, resCols, localWorkSize)
+    let result' = Array.zeroCreate (resRows * resCols)
+
+    kernelPrepare d m1Arr m2Arr result'
+
+    commandQueue.Add(kernelRun()) |> ignore
+    let _ = commandQueue.Add(result'.ToHost provider).Finish()
+
+    let result = Array2D.zeroCreate resRows resCols
+
+    Array.Parallel.iteri (fun x v -> result.[x / resCols, x % resCols] <-  v) result'
+
+    provider.CloseAllBuffers()
+
+    result

src/ImgProcessing/Matrices.fs

@@ -0,0 +1,82 @@
+module Lecture10.Interfaces
+
+type IMatrix<'elem> =
+    abstract Item: int * int -> 'elem with get, set
+    abstract Size: int * int with get
+
+type MyMatrix<'elem> (n,m) =
+    let store = Array2D.zeroCreate n m
+    interface IMatrix<'elem> with
+        member this.Item
+            with get (i, j) = store.[i, j]
+            and set (i, j) v = store.[i, j] <- v
+
+        member this.Size = n,m
+
+type ISmb<'t> =
+    abstract PrettyString: string
+
+type Smb<'t> =
+    | Eps
+    | Smb of 't
+    interface ISmb<'t> with
+        member this.PrettyString =
+            match this with
+            | Eps -> "Eps"
+            | Smb t -> t.ToString()
+
+[<AbstractClass>]
+type FA<'smb> () =
+    abstract StartState: int with get//, set
+    abstract FinalState: int with get//, set
+    abstract TransitionsList : array<int * int * ISmb<'smb>>
+    member this.ToDot (outFile) =
+        let header =
+            [
+                "digraph nfa"
+                "{"
+                "rankdir = LR"
+                "node [shape = circle];"
+                sprintf "%A[shape = circle, label = \"%A_Start\"]" this.StartState this.StartState
+            ]
+
+        let footer =
+            [
+                sprintf "%A[shape = doublecircle]" this.FinalState
+                "}"
+            ]
+
+        let content =
+            [
+                for (s,f,t) in this.TransitionsList do
+                    sprintf
+                        "%A -> %A [label = \"%s\"]"
+                        s
+                        f
+                        t.PrettyString
+            ]
+
+        System.IO.File.WriteAllLines (outFile, header @ content @ footer)
+
+
+type MyDFA<'smb> (transitions, startSate, finalState) =
+    inherit FA<'smb>()
+    override this.StartState = startSate
+    override this.FinalState = finalState
+    override this.TransitionsList =
+        let n,m = this.Transitions.Size
+        [|
+         for i in 0 .. n - 1 do
+            for j in 0 .. m - 1 do
+                yield i, j, this.Transitions.[i,j]
+        |]
+
+    member this.Transitions: IMatrix<ISmb<'smb>> = transitions
+
+
+
+let m = new MyMatrix<_>(2,3) :> IMatrix<_>
+m.[1,2] <- Smb(0) :> ISmb<_>
+
+let fa = new MyDFA<_>(m,1,0)
+fa.ToDot "file"

src/Lecture10/Interfaces.fs

@@ -22,6 +22,7 @@
     <Compile Include="Automata.fs" />
     <Compile Include="Regexp.fs" />
     <Compile Include="Main.fs" />
+    <Compile Include="Interfaces.fs" />
   </ItemGroup>
   <Import Project="..\..\.paket\Paket.Restore.targets" />
 </Project>