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
use super::super::{
Buffer, BufferBorrow, ComplexNumberSpace, DataDomain, DspVec, FrequencyDomain, MetaData,
RededicateForceOps, ToSliceMut, ToTimeResult, Vector,
};
use numbers::*;
pub trait FrequencyDomainOperations<S, T>
where
S: ToSliceMut<T>,
T: RealNumber,
{
fn mirror<B>(&mut self, buffer: &mut B)
where
B: for<'a> Buffer<'a, S, T>;
fn fft_shift(&mut self);
fn ifft_shift(&mut self);
}
impl<S, T, N, D> FrequencyDomainOperations<S, T> for DspVec<S, T, N, D>
where
DspVec<S, T, N, D>: ToTimeResult,
<DspVec<S, T, N, D> as ToTimeResult>::TimeResult: RededicateForceOps<DspVec<S, T, N, D>>,
S: ToSliceMut<T>,
T: RealNumber,
N: ComplexNumberSpace,
D: FrequencyDomain,
{
fn mirror<B>(&mut self, buffer: &mut B)
where
B: for<'a> Buffer<'a, S, T>,
{
if self.domain() != DataDomain::Frequency && !self.is_complex() {
self.valid_len = 0;
return;
}
let len = self.len();
let step = 2;
let temp_len = 2 * len - step;
let mut temp = buffer.borrow(temp_len);
{
let data = self.data.to_slice();
let temp = temp.to_slice_mut();
temp[step..len].clone_from_slice(&data[step..len]);
temp[0..step].clone_from_slice(&data[0..step]);
let mut j = step + 1;
let mut i = temp_len - 1;
while i >= len {
temp[i] = -data[j];
temp[i - 1] = data[j - 1];
i -= 2;
j += 2;
}
self.valid_len = temp_len;
}
temp.trade(&mut self.data);
}
fn fft_shift(&mut self) {
self.swap_halves_priv(true)
}
fn ifft_shift(&mut self) {
self.swap_halves_priv(false)
}
}