Add vDSP_vclrD, vDSP_vfillD

Author: ptoffy

README.md

@@ -54,8 +54,8 @@ The package is structured as follows:
         - [x] `vDSP_vsubD` - Subtract two vectors
         - [x] `vDSP_vmulD` - Multiply two vectors
         - [x] `vvpow` - Raise a vector to a power
-        - [ ] `vDSP_vclrD` - Clear a vector
-        - [ ] `vDSP_vfillD` - Fill a vector with a value
+        - [x] `vDSP_vclrD` - Clear a vector
+        - [x] `vDSP_vfillD` - Fill a vector with a value
         - [ ] `vDSP_vabsD` - Absolute value of a vector
         - [ ] `vDSP_vnegD` - Negate a vector
         - [ ] `vDSP_vsqD` - Square a vector

Sources/AccelerateLinux/VectorOps/VectorBasicOps.swift

@@ -151,4 +151,32 @@ public func vvpow(
         i += 1
     }
 }
+
+/// Populates a double-precision vector with zeros.
+/// - Parameters:
+///   - __C: Double-precision real output vector.
+///   - __IC: Address stride for C.
+///   - __N: The number of elements to process.
+public func vDSP_vclrD(
+    _ __C: UnsafeMutablePointer<Double>,
+    _ __IC: vDSP_Stride,
+    _ __N: vDSP_Length
+) {
+    var zero: Double = 0
+    vDSP_vfillD(&zero, __C, __IC, __N)
+}
+
+/// Populates a double-precision vector with a specified scalar value.
+public func vDSP_vfillD(
+    _ __A: UnsafePointer<Double>,
+    _ __C: UnsafeMutablePointer<Double>,
+    _ __IC: vDSP_Stride,
+    _ __N: vDSP_Length
+) {
+    var i = 0
+    while i < __N {
+        __C[i * __IC] = __A.pointee
+        i += 1
+    }
+}
 #endif

Tests/AccelerateLinuxTests/VectorTests/VectorBasicOpsTests.swift

@@ -1,4 +1,5 @@
 import AccelerateLinux
+import FoundationEssentials
 import Testing
 
 @Suite("Vector Basic Ops Tests")
@@ -130,4 +131,44 @@ struct VectorBasicOpsTests {
         vvpow(&z, &y, &x, &n)
         #expect(z == [9.0, 16.0, 1000.0, 36.0])
     }
+
+    // https://developer.apple.com/documentation/accelerate/1450639-vdsp_vclrd
+    @Test("vDSP_vclrD")
+    func vDSP_vclrDTest() {
+        let n = vDSP_Length(10)
+        let stride = vDSP_Stride(1)
+
+        var c = [Double](
+            repeating: .nan,
+            count: Int(n))
+
+        vDSP_vclrD(&c, stride, n)
+
+        #expect(c == [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0])
+    }
+
+    // https://developer.apple.com/documentation/accelerate/1450171-vdsp_vfilld
+    @Test("vDSP_vfillD")
+    func vDSP_vfillDTest() {
+        let n = vDSP_Length(10)
+        let stride = vDSP_Stride(1)
+
+        var a = Double.pi
+
+        var c = [Double](
+            repeating: .nan,
+            count: Int(n))
+
+        vDSP_vfillD(&a, &c, stride, n)
+
+        #expect(
+            c.map {
+                ($0 * 100).rounded() / 100
+            } == [
+                3.14, 3.14, 3.14,
+                3.14, 3.14, 3.14,
+                3.14, 3.14, 3.14,
+                3.14,
+            ])
+    }
 }