Add vDSP_mmulD

Author: ptoffy

README.md

@@ -35,7 +35,7 @@ The package is structured as follows:
 - `Sources/AccelerateLinux/MatrixOps`: Contains matrix operations:
     - `BasicOps` - Basic matrix operations:
         - [x] `vDSP_mtransD` - Transpose a matrix
-        - [ ] `vDSP_mmulD` - Matrix multiplication
+        - [x] `vDSP_mmulD` - Matrix multiplication
     - `LAPACK` - LAPACK functions:
         - [x] `dgesv_` - Solve a system of linear equations
         - [x] `dgesvd_` - Singular Value Decomposition

Sources/AccelerateLinux/MatrixOps/MatrixBasicOps.swift

@@ -53,6 +53,18 @@ public func vDSP_mmulD(
     _ __N: vDSP_Length,
     _ __P: vDSP_Length
 ) {
-
+    var a = __A
+    for i in 0..<__M {
+        var b = __B
+        for _ in 0..<__P {
+            var c = __C.advanced(by: Int(i) * __IC * Int(__N))
+            for _ in 0..<__N {
+                c.pointee += a.pointee * b.pointee
+                b = b.advanced(by: __IB)
+                c = c.advanced(by: __IC)
+            }
+            a = a.advanced(by: __IA)
+        }
+    }
 }
 #endif  // canImport(Accelerate)