legato 0.0.10

Legato is a WIP audiograph and DSL for quickly developing audio applications
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
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88

use std::simd::{
    StdFloat,
    num::{SimdFloat, SimdUint},
};

use crate::{
    math::{ONE_VUSIZE, TWO_VUSIZE, cubic_hermite, cubic_hermite_simd, lerp, lerp_simd},
    simd::{LANES, Vf32, Vusize},
};

/// It became clear after a while, that delay lines,
/// samplers, granular, etc. need the same underlying
/// abstraction, a window into a slice, with a few fractional
/// indexing utilities.
///
/// This primative can be used to make delay lines, samplers, etc.
///
/// This also allows us to handle all of these resources as views
/// into one giant continous buffer, which should give better cache locality.
pub struct Window<'a> {
    data: &'a [f32],
}

impl<'a> Window<'a> {
    /// Construct directly from a slice. The slice should already be
    /// positioned so that index 0 is the first logical sample.
    #[inline(always)]
    pub fn new(data: &'a [f32]) -> Self {
        Self { data }
    }

    #[inline(always)]
    pub fn read_linear(&self, offset: f32) -> f32 {
        let i = offset as usize;
        let t = offset - i as f32;
        unsafe {
            let a = *self.data.get_unchecked(i);
            let b = *self.data.get_unchecked(i + 1);
            lerp(a, b, t)
        }
    }

    #[inline(always)]
    pub fn read_cubic(&self, offset: f32) -> f32 {
        let i = offset as usize;
        let t = offset - i as f32;
        unsafe {
            let p0 = *self.data.get_unchecked(i.wrapping_sub(1));
            let p1 = *self.data.get_unchecked(i);
            let p2 = *self.data.get_unchecked(i + 1);
            let p3 = *self.data.get_unchecked(i + 2);
            cubic_hermite(p0, p1, p2, p3, t)
        }
    }

    #[inline(always)]
    pub fn read_linear_simd(&self, offsets: Vf32) -> Vf32 {
        let floor_float = offsets.floor();
        let floor_usize = floor_float.cast::<usize>();
        let t = offsets - floor_float;
        let a = self.gather(floor_usize);
        let b = self.gather(floor_usize + Vusize::splat(1));
        lerp_simd(a, b, t)
    }

    #[inline(always)]
    pub fn read_cubic_simd(&self, offsets: Vf32) -> Vf32 {
        let floor_float = offsets.floor();
        let floor_usize = floor_float.cast::<usize>();
        let t = offsets - floor_float;
        let a = self.gather(floor_usize.saturating_sub(ONE_VUSIZE));
        let b = self.gather(floor_usize);
        let c = self.gather(floor_usize + ONE_VUSIZE);
        let d = self.gather(floor_usize + TWO_VUSIZE);
        cubic_hermite_simd(a, b, c, d, t)
    }

    #[inline(always)]
    fn gather(&self, indices: Vusize) -> Vf32 {
        let mut out = [0.0f32; LANES];
        for lane in 0..LANES {
            unsafe {
                out[lane] = *self.data.get_unchecked(indices[lane]);
            }
        }
        Vf32::from_array(out)
    }
}