0001    // Auxilliary.swift
0002    //
0003    // Copyright (c) 2014–2015 Mattt Thompson (http://mattt.me)//
0005    // Permission is hereby granted, free of charge, to any person obtaining a copy
0006    // of this software and associated documentation files (the "Software"), to deal
0007    // in the Software without restriction, including without limitation the rights
0008    // to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
0009    // copies of the Software, and to permit persons to whom the Software is
0010    // furnished to do so, subject to the following conditions:
0011    //
0012    // The above copyright notice and this permission notice shall be included in
0013    // all copies or substantial portions of the Software.
0014    //
0015    // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
0016    // IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
0017    // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
0018    // AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
0019    // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
0020    // OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
0021    // THE SOFTWARE.
0022    
0023    import Accelerate
0024    
0025    // MARK: Absolute Value
0026    
0027    public func abs(x: [Double]) -> [Double] {
0028        var results = [Double](count: x.count, repeatedValue: 0.0)
0029        vvfabs(&results, x, [Int32(x.count)])
0030    
0031        return results
0032    }
0033    
0034    public func abs(x: [Float]) -> [Float] {
0035        var results = [Float](count: x.count, repeatedValue: 0.0)
0036        vvfabsf(&results, x, [Int32(x.count)])
0037    
0038        return results
0039    }
0040    
0041    // MARK: Ceiling
0042    
0043    public func ceil(x: [Float]) -> [Float] {
0044        var results = [Float](count: x.count, repeatedValue: 0.0)
0045        vvceilf(&results, x, [Int32(x.count)])
0046    
0047        return results
0048    }
0049    
0050    public func ceil(x: [Double]) -> [Double] {
0051        var results = [Double](count: x.count, repeatedValue: 0.0)
0052        vvceil(&results, x, [Int32(x.count)])
0053    
0054        return results
0055    }
0056    
0057    // MARK: Clip
0058    
0059    public func clip(x: [Float], low: Float, high: Float) -> [Float] {
0060        var results = [Float](count: x.count, repeatedValue: 0.0), y = low, z = high
0061        vDSP_vclip(x, , &y, &z, &results, , vDSP_Length(x.count))
0062    
0063        return results
0064    }
0065    
0066    public func clip(x: [Double], low: Double, high: Double) -> [Double] {
0067        var results = [Double](count: x.count, repeatedValue: 0.0), y = low, z = high
0068        vDSP_vclipD(x, , &y, &z, &results, , vDSP_Length(x.count))
0069    
0070        return results
0071    }
0072    
0073    // MARK: Copy Sign
0074    
0075    public func copysign(sign: [Float], magnitude: [Float]) -> [Float] {
0076        var results = [Float](count: sign.count, repeatedValue: 0.0)
0077        vvcopysignf(&results, magnitude, sign, [Int32(sign.count)])
0078    
0079        return results
0080    }
0081    
0082    public func copysign(sign: [Double], magnitude: [Double]) -> [Double] {
0083        var results = [Double](count: sign.count, repeatedValue: 0.0)
0084        vvcopysign(&results, magnitude, sign, [Int32(sign.count)])
0085    
0086        return results
0087    }
0088    
0089    // MARK: Floor
0090    
0091    public func floor(x: [Float]) -> [Float] {
0092        var results = [Float](count: x.count, repeatedValue: 0.0)
0093        vvfloorf(&results, x, [Int32(x.count)])
0094    
0095        return results
0096    }
0097    
0098    public func floor(x: [Double]) -> [Double] {
0099        var results = [Double](count: x.count, repeatedValue: 0.0)
0100        vvfloor(&results, x, [Int32(x.count)])
0101    
0102        return results
0103    }
0104    
0105    // MARK: Negate
0106    
0107    public func neg(x: [Float]) -> [Float] {
0108        var results = [Float](count: x.count, repeatedValue: 0.0)
0109        vDSP_vneg(x, , &results, , vDSP_Length(x.count))
0110    
0111        return results
0112    }
0113    
0114    public func neg(x: [Double]) -> [Double] {
0115        var results = [Double](count: x.count, repeatedValue: 0.0)
0116        vDSP_vnegD(x, , &results, , vDSP_Length(x.count))
0117    
0118        return results
0119    }
0120    
0121    // MARK: Reciprocal
0122    
0123    public func rec(x: [Float]) -> [Float] {
0124        var results = [Float](count: x.count, repeatedValue: 0.0)
0125        vvrecf(&results, x, [Int32(x.count)])
0126    
0127        return results
0128    }
0129    
0130    public func rec(x: [Double]) -> [Double] {
0131        var results = [Double](count: x.count, repeatedValue: 0.0)
0132        vvrec(&results, x, [Int32(x.count)])
0133    
0134        return results
0135    }
0136    
0137    // MARK: Round
0138    
0139    public func round(x: [Float]) -> [Float] {
0140        var results = [Float](count: x.count, repeatedValue: 0.0)
0141        vvnintf(&results, x, [Int32(x.count)])
0142    
0143        return results
0144    }
0145    
0146    public func round(x: [Double]) -> [Double] {
0147        var results = [Double](count: x.count, repeatedValue: 0.0)
0148        vvnint(&results, x, [Int32(x.count)])
0149    
0150        return results
0151    }
0152    
0153    // MARK: Threshold
0154    
0155    public func threshold(x: [Float], low: Float) -> [Float] {
0156        var results = [Float](count: x.count, repeatedValue: 0.0), y = low
0157        vDSP_vthr(x, , &y, &results, , vDSP_Length(x.count))
0158    
0159        return results
0160    }
0161    
0162    public func threshold(x: [Double], low: Double) -> [Double] {
0163        var results = [Double](count: x.count, repeatedValue: 0.0), y = low
0164        vDSP_vthrD(x, , &y, &results, , vDSP_Length(x.count))
0165    
0166        return results
0167    }
0168    
0169    // MARK: Truncate
0170    
0171    public func trunc(x: [Float]) -> [Float] {
0172        var results = [Float](count: x.count, repeatedValue: 0.0)
0173        vvintf(&results, x, [Int32(x.count)])
0174    
0175        return results
0176    }
0177    
0178    public func trunc(x: [Double]) -> [Double] {
0179        var results = [Double](count: x.count, repeatedValue: 0.0)
0180        vvint(&results, x, [Int32(x.count)])
0181    
0182        return results
0183    }
0184