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

ix!();

pub trait WindowInput = 
std::cmp::PartialOrd 
+ std::ops::Mul<Output=f64> 
+ From<i32> 
+ std::ops::SubAssign 
+ Into<f64>;

pub fn blackman< T: WindowInput>(i: T, n: i32) 
    -> f64 
{
    let i: f64 = i.into();
    let n: f64 = n.into();
    0.42 - 
        0.5 * (2.0 * PI  * i / (n - 1.0)).cos() + 
        0.08 * (4.0 * PI * i / (n - 1.0)).cos()
}

pub fn symmetric_blackman< T: WindowInput >(i: T, n: i32) 
    -> f64 
{
    let mut i: f64 = i.into();
    let n: f64 = n.into();
    i -= n / 2.0;
    0.42 - 
        0.5 * (2.0 * PI  * i / (n)).cos() + 
        0.08 * (4.0 * PI * i / (n)).cos()
}

pub fn blackman_harris< T: WindowInput >(i: T, n: i32) 
    -> f64 
{
    let i: f64 = i.into();
    let n: f64 = n.into();
    0.35875 - 
        0.48829 * (2.0 * PI * i / (n - 1.0)).cos() + 
        0.14128 * (4.0 * PI * i / (n - 1.0)).cos() -
        0.01168 * (6.0 * PI * i / (n - 1.0)).cos()
}

pub fn symmetric_blackman_harris< T: WindowInput >(i: T, n: i32) 
    -> f64 
{
    let mut i: f64 = i.into();
    let n: f64 = n.into();
    i -= n / 2.0;
    0.35875 - 
        0.48829 * (2.0 * PI * i / (n)).cos() + 
        0.14128 * (4.0 * PI * i / (n - 1.0)).cos() -
        0.01168 * (6.0 * PI * i / (n)).cos()
}

pub fn hamming< T: WindowInput>(i: T, n: i32) 
    -> f64 
{
    let i: f64 = i.into();
    let n: f64 = n.into();

    if i >= n {
        return 0.0;
    }

    0.54 - 
        0.46 * (2.0 * PI * i / (n - 1.0)).cos()
}

pub fn hanning< T: WindowInput >( i: T, n: i32) 
    -> f64 
{
    let i: f64 = i.into();
    let n: f64 = n.into();

    if i >= n {
        return 0.0;
    }

    0.5 * 
        (1.0 - (2.0 * PI * i / (n - 1.0)).cos())
}