From ea1982f5943075646ea3029d491882fcbd5a4252 Mon Sep 17 00:00:00 2001
From: Don Syme <dsyme@github.com>
Date: Sun, 12 Oct 2025 13:28:01 +0100
Subject: [PATCH 1/2] Added comprehensive matrix operation benchmarks (14
 operations)

---
 benchmarks/FsMath.Benchmarks/FsMath.Benchmarks.fsproj | 1 +
 benchmarks/FsMath.Benchmarks/Program.fs               | 3 ++-
 2 files changed, 3 insertions(+), 1 deletion(-)

diff --git a/benchmarks/FsMath.Benchmarks/FsMath.Benchmarks.fsproj b/benchmarks/FsMath.Benchmarks/FsMath.Benchmarks.fsproj
index 6b13aa6..e32e1c0 100644
--- a/benchmarks/FsMath.Benchmarks/FsMath.Benchmarks.fsproj
+++ b/benchmarks/FsMath.Benchmarks/FsMath.Benchmarks.fsproj
@@ -7,6 +7,7 @@
 
   <ItemGroup>
     <Compile Include="Vector.fs" />
+    <Compile Include="Matrix.fs" />
     <Compile Include="Program.fs" />
   </ItemGroup>
 
diff --git a/benchmarks/FsMath.Benchmarks/Program.fs b/benchmarks/FsMath.Benchmarks/Program.fs
index d3c9b13..3c8a441 100644
--- a/benchmarks/FsMath.Benchmarks/Program.fs
+++ b/benchmarks/FsMath.Benchmarks/Program.fs
@@ -5,8 +5,9 @@ open FsMath.Benchmarks
 [<EntryPoint>]
 let Main args =
     // Register multiple benchmark classes
-    let switcher = BenchmarkSwitcher [| 
+    let switcher = BenchmarkSwitcher [|
         typeof<VectorBenchmarks>
+        typeof<MatrixBenchmarks>
     |]
     switcher.Run args |> ignore
     0

From 12868c58a363cdc9387fd6ce643297a6da99a0cb Mon Sep 17 00:00:00 2001
From: Don Syme <dsyme@github.com>
Date: Sun, 12 Oct 2025 13:29:17 +0100
Subject: [PATCH 2/2] Added comprehensive matrix operation benchmarks (14
 operations)

---
 benchmarks/FsMath.Benchmarks/Matrix.fs | 107 +++++++++++++++++++++++++
 1 file changed, 107 insertions(+)
 create mode 100644 benchmarks/FsMath.Benchmarks/Matrix.fs

diff --git a/benchmarks/FsMath.Benchmarks/Matrix.fs b/benchmarks/FsMath.Benchmarks/Matrix.fs
new file mode 100644
index 0000000..30d9d15
--- /dev/null
+++ b/benchmarks/FsMath.Benchmarks/Matrix.fs
@@ -0,0 +1,107 @@
+namespace FsMath.Benchmarks
+
+open System
+open BenchmarkDotNet.Attributes
+open FsMath
+
+[<MemoryDiagnoser>]
+type MatrixBenchmarks() =
+
+    let mutable matrixA = Unchecked.defaultof<Matrix<float>>
+    let mutable matrixB = Unchecked.defaultof<Matrix<float>>
+    let mutable vector = [||]
+
+    // Parameterize matrix sizes: small, medium
+    [<Params(10, 50, 100)>]
+    member val Size = 0 with get, set
+
+    [<GlobalSetup>]
+    member this.Setup() =
+        // Initialize square matrices with some values
+        matrixA <- Matrix.init<float> this.Size this.Size (fun i j -> float (i + j))
+        matrixB <- Matrix.init<float> this.Size this.Size (fun i j -> float (i * 2 + j))
+        vector <- Array.init this.Size (fun i -> float i)
+
+    // Element-wise operations
+
+    [<Benchmark>]
+    member _.ElementWiseAdd() =
+        let result = Matrix.add matrixA matrixB
+        GC.KeepAlive(result)
+
+    [<Benchmark>]
+    member _.ElementWiseSubtract() =
+        let result = Matrix.subtract matrixA matrixB
+        GC.KeepAlive(result)
+
+    [<Benchmark>]
+    member _.ElementWiseMultiply() =
+        let result = Matrix.multiply matrixA matrixB
+        GC.KeepAlive(result)
+
+    [<Benchmark>]
+    member _.ElementWiseDivide() =
+        let result = Matrix.divide matrixA matrixB
+        GC.KeepAlive(result)
+
+    // Scalar operations
+
+    [<Benchmark>]
+    member _.ScalarAdd() =
+        let result = Matrix.addScalar matrixA 5.0
+        GC.KeepAlive(result)
+
+    [<Benchmark>]
+    member _.ScalarMultiply() =
+        let result = Matrix.multiplyScalar matrixA 2.5
+        GC.KeepAlive(result)
+
+    // Matrix multiplication
+
+    [<Benchmark>]
+    member _.MatrixMultiply() =
+        let result = Matrix.matmul matrixA matrixB
+        GC.KeepAlive(result)
+
+    // Matrix-vector operations
+
+    [<Benchmark>]
+    member _.MatrixVectorMultiply() =
+        let result = Matrix.muliplyVector matrixA vector
+        GC.KeepAlive(result)
+
+    [<Benchmark>]
+    member _.VectorMatrixMultiply() =
+        let result = Matrix.multiplyRowVector vector matrixA
+        GC.KeepAlive(result)
+
+    // Transpose
+
+    [<Benchmark>]
+    member _.Transpose() =
+        let result = matrixA.Transpose()
+        GC.KeepAlive(result)
+
+    // Row/column access patterns
+
+    [<Benchmark>]
+    member this.GetRow() =
+        let result = Matrix.getRow (this.Size / 2) matrixA
+        GC.KeepAlive(result)
+
+    [<Benchmark>]
+    member this.GetCol() =
+        let result = Matrix.getCol (this.Size / 2) matrixA
+        GC.KeepAlive(result)
+
+    // Broadcast operations
+
+    [<Benchmark>]
+    member _.AddRowVector() =
+        let result = Matrix.addRowVector matrixA vector
+        GC.KeepAlive(result)
+
+    [<Benchmark>]
+    member _.AddColVector() =
+        let result = Matrix.addColVector matrixA vector
+        GC.KeepAlive(result)