signal_processing 0.3.0

A signal processing library.
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
50
51
52
53
54
55

use core::ops::{AddAssign, BitAnd, Shr};

use num::{complex::ComplexFloat, pow::Pow, traits::Inv, Integer, One};

use crate::{quantities::MaybeList, operations::Simplify, systems::Tf, transforms::system::ToTf};

impl<T, B, A, I> Pow<I> for Tf<T, B, A>
where
    T: ComplexFloat + AddAssign,
    B: MaybeList<T>,
    A: MaybeList<T>,
    I: Integer + BitAnd<I, Output = I> + Shr<usize, Output = I> + Copy,
    Self: ToTf<T, Vec<T>, Vec<T>, (), ()> + Inv<Output: ToTf<T, Vec<T>, Vec<T>, (), ()>>,
    Tf<T, Vec<T>, Vec<T>>: Simplify<Output = Tf<T, Vec<T>, Vec<T>>>
{
    type Output = Tf<T, Vec<T>, Vec<T>>;

    fn pow(self, mut n: I) -> Self::Output
    {
        let mut x = if n < I::zero()
        {
            self.inv()
                .to_tf((), ())
        }
        else
        {
            self.to_tf((), ())
        };
        let mut r = if (n & I::one()) == I::one()
        {
            x.clone()
        }
        else
        {
            One::one()
        };
    
        let two = I::one() + I::one();
        loop
        {
            n = n/two;
            if n == I::zero()
            {
                break;
            }
            x = x.as_view()*x.as_view();
            if (n & I::one()) == I::one()
            {
                r = r*x.as_view();
            }
        }
    
        r
    }
}