module MatrixMult
open System
open System.Threading

let calcEl(l, r , max, str1, stol1, long, m1:int[,], m2:int[,], resMatrix:int[,]) =
 if (r<max) then do
  let mutable res = 0
  for x=l to r do
   res <- 0
   let str = x / stol1
   let stol = x%stol1
   for i=0 to long-1 do
    res <- res + (m1.[str,i] * m2.[i,stol])
    resMatrix.[str,stol] <- res
  resMatrix


let multiply (m1: int[,] , m2: int[,] , threadNumber:int) =
 let l11 = m1.GetLength 0
 let l12 = m1.GetLength 1

 let l21 = m2.GetLength 0
 let l22 = m2.GetLength 1
 let resLength = l11*l22
 let resMatrix = (Array2D.zeroCreate l11 l22)
 let step = resLength / threadNumber
 let threadArray = Array.init (threadNumber) (fun i -> new Thread(ThreadStart(fun _ ->
    calcEl (i*step , (i+1)*step, resLength, l11, l22, l12, m1, m2, resMatrix)
    )))
 calcEl ((threadNumber-1)*step, resLength-1, resLength, l11,l22,l12,m1,m2,resMatrix)

 for t in threadArray do
  t.Start()
 for t in threadArray do
  t.Join()
  resMatrix



[<EntryPoint>]
let main argv =
 0