-
Notifications
You must be signed in to change notification settings - Fork 75
Expand file tree
/
Copy pathmod.rs
More file actions
196 lines (181 loc) · 5.67 KB
/
mod.rs
File metadata and controls
196 lines (181 loc) · 5.67 KB
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
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
//! Items to ease the application of windowing functions to signals.
use crate::{ConstHz, FromIterator, Phase, Signal};
use core::marker::PhantomData;
use dasp_frame::Frame;
use dasp_sample::Sample;
use dasp_window::Window as WindowType;
#[cfg(feature = "window-hann")]
pub use hann::hann;
#[cfg(feature = "window-rectangle")]
pub use rectangle::rectangle;
#[cfg(feature = "window-hann")]
mod hann;
#[cfg(feature = "window-rectangle")]
mod rectangle;
/// A `Signal` type that for every yielded `phase`, yields the amplitude across the `window::Type`
/// for that phase.
///
/// ### Required Features
///
/// - When using `dasp_signal`, this item requires the **window** feature to be enabled.
/// - When using `dasp`, this item requires the **signal-window** feature to be enabled.
#[derive(Clone)]
pub struct Window<F, W>
where
F: Frame,
W: WindowType<f64, Output = f64>,
{
/// Yields phase stepped at a constant rate to be passed to the window function `W`.
pub phase: Phase<ConstHz>,
/// Size of the window
pub len: usize,
marker: PhantomData<(F, W)>,
}
/// Takes a long slice of frames and yields `Windowed` chunks of size `bin` once every `hop` frames.
///
/// ### Required Features
///
/// - When using `dasp_signal`, this item requires the **window** feature to be enabled.
/// - When using `dasp`, this item requires the **signal-window** feature to be enabled.
#[derive(Clone)]
pub struct Windower<'a, F, W>
where
F: 'a + Frame,
W: WindowType<f64, Output = f64>,
{
/// The size of each `Windowed` chunk to be yielded.
pub bin: usize,
/// The step size over `frames` for the start of each new `Windowed` chunk.
pub hop: usize,
/// The beginning of the remaining slice to be yielded by the `Windower`.
pub frames: &'a [F],
wttype: PhantomData<W>,
}
/// An Iterator that multiplies a Signal with a Window.
///
/// Returns `None` once the `Window` has been exhausted.
///
/// ### Required Features
///
/// - When using `dasp_signal`, this item requires the **window** feature to be enabled.
/// - When using `dasp`, this item requires the **signal-window** feature to be enabled.
#[derive(Clone)]
pub struct Windowed<S, W>
where
S: Signal,
W: WindowType<f64, Output = f64>,
{
signal: S,
window: Window<<S::Frame as Frame>::Float, W>,
}
impl<F, W> Window<F, W>
where
F: Frame,
W: WindowType<f64, Output = f64>,
{
/// Construct a new `Window` with the given length as a number of frames.
///
/// ### Required Features
///
/// - When using `dasp_signal`, this item requires the **window** feature to be enabled.
/// - When using `dasp`, this item requires the **signal-window** feature to be enabled.
pub fn new(len: usize) -> Self {
let step = crate::rate(len as f64 - 1.0).const_hz(1.0);
Window {
len,
phase: crate::phase(step),
marker: PhantomData,
}
}
}
impl<'a, F, W> Windower<'a, F, W>
where
F: 'a + Frame,
W: WindowType<f64, Output = f64>,
{
/// Constructor for a new `Windower` iterator.
///
/// ### Required Features
///
/// - When using `dasp_signal`, this item requires the **window** feature to be enabled.
/// - When using `dasp`, this item requires the **signal-window** feature to be enabled.
pub fn new(frames: &'a [F], bin: usize, hop: usize) -> Self {
Windower {
bin: bin,
hop: hop,
frames: frames,
wttype: PhantomData,
}
}
}
impl<F, W> Iterator for Window<F, W>
where
F: Frame,
W: WindowType<f64, Output = f64>,
{
type Item = F;
fn next(&mut self) -> Option<Self::Item> {
// make sure we didn't produce items indefinitely
if self.len == 0 {
return None;
}
self.len -= 1;
let v = W::window(self.phase.next_phase());
let v_f: <F::Sample as Sample>::Float = v.to_sample();
Some(F::from_fn(|_| v_f.to_sample::<F::Sample>()))
}
}
impl<'a, F, W> Iterator for Windower<'a, F, W>
where
F: 'a + Frame,
W: WindowType<f64, Output = f64>,
{
type Item = Windowed<FromIterator<core::iter::Cloned<core::slice::Iter<'a, F>>>, W>;
fn next(&mut self) -> Option<Self::Item> {
let num_frames = self.frames.len();
if self.bin <= num_frames {
let frames = &self.frames[..self.bin];
let window = Window::new(self.bin);
self.frames = if self.hop < num_frames {
&self.frames[self.hop..]
} else {
&[]
};
Some(Windowed {
signal: crate::from_iter(frames.iter().cloned()),
window: window,
})
} else {
None
}
}
fn size_hint(&self) -> (usize, Option<usize>) {
let num_frames = self.frames.len();
// Must have at least `bin` number of frames left to iterate at all.
if self.bin < num_frames {
// If the hop size is 0, we'll iterate forever.
if self.hop == 0 {
return (core::usize::MAX, None);
}
// Otherwise we can determine exactly how many iterations remain.
let remaining_hop_frames = self.frames.len() - self.bin;
let remaining_iterations = remaining_hop_frames / self.hop;
(remaining_iterations, Some(remaining_iterations))
} else {
(0, Some(0))
}
}
}
impl<S, W> Iterator for Windowed<S, W>
where
S: Signal,
W: WindowType<f64, Output = f64>,
{
type Item = S::Frame;
fn next(&mut self) -> Option<Self::Item> {
self.window.next().map(|w_f| {
let s_f = self.signal.next();
s_f.mul_amp(w_f)
})
}
}