dsp 0.11.0

Digital Signal Processing
Documentation 
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

/// Helper functions for vector operations based on block of real data
/// All operations here should be alloc free 
use std::ops::Mul;
use std::f32;


/// Multiply vectors element wise
pub fn multiply<T>(xs: &[T], ys: &[T], output: &mut [T]) 
where
    T: Mul<Output = T> + Copy
{
    assert_eq!(xs.len(), output.len());
    assert_eq!(ys.len(), output.len());
    for i in 0..output.len() {
        output[i] = xs[i] * ys[i];
    }
}

/// Calculate arg max for complex numbers
pub fn argmax(xs: &[f32]) -> usize {
    let mut max_value = f32::MIN;
    let mut arg_max = 0;
    for i in 0..xs.len() {
        let v = xs[i];
        if max_value < v {
            max_value = v;
            arg_max = i;
        }
    }
    arg_max
}


/// ------------------------------------------------------------------------------------------------
/// Module unit tests
/// ------------------------------------------------------------------------------------------------
#[cfg(test)]
mod tests {
    use super::*;

    #[test]
    fn test_multiply() {
        let v = vec![3.0; 5];
        let u = vec![2.0; 5];
        let mut output = vec![0.0; 5];
        multiply(&v, &u, &mut output);
        assert_eq!(output, vec![6.0; 5]);
    }
}