surge-math 0.2.6-alpha.0

surge synthesizer -- dsp math used in synthesis
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69

crate::ix!();

pub trait LimitRange: Sized {

    fn limit_range(self, low: Self, hi: Self) -> Self;

    /**
      | provides a possibility for taking a
      | non-simd path. if we don't customize it,
      | the ordinary simd code will be called if it
      | is used in "limit_range"
      |
      */
    #[inline] fn limit_range_nosimd(self, low: Self, hi: Self) -> Self {
        self.limit_range(low,hi)
    }
}

pub fn limit_range<T: LimitRange>(x: T, low: T, high: T) -> T {
    x.limit_range(low,high)
}

impl LimitRange for f32 {

    #[cfg(target_arch = "x86_64")] #[inline] 
    fn limit_range(self, low: Self, high: Self) -> Self {

        let mut result: f32 = 0.0;

        unsafe {

            let _s    = _mm_load_ss(&self);
            let _low  = _mm_load_ss(&low);
            let _high = _mm_load_ss(&high);
            let limited = _mm_min_ss(_mm_max_ss(_s, _low), _high);

            _mm_store_ss( &mut result, limited) 
        };

        result
    }
    fn limit_range_nosimd(self, min: f32, max: f32) -> f32 {
        std::cmp::max(std::cmp::min(FloatOrd(self),FloatOrd(max)), FloatOrd(min)).0
    }
}

impl LimitRange for i32 {

    fn limit_range(self, l: i32, h: i32) -> i32 
    {
        std::cmp::max(std::cmp::min(self,h), l)
    }
}

impl LimitRange for u32 {
    fn limit_range(self, l: u32, h: u32) -> u32 
    {
        std::cmp::max(std::cmp::min(self,h), l)
    }
}


impl LimitRange for f64 {
    fn limit_range(self, min:f64, max:f64) -> f64 
    {
        std::cmp::max(std::cmp::min(FloatOrd(self),FloatOrd(max)), FloatOrd(min)).0
    }
}