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

use core::ops::Add;

use num::{complex::ComplexFloat, Float, Zero};

use crate::quantities::{List, ListOrSingle, Lists};

pub trait Rssq<'a, T>: Lists<T>
where
    T: ComplexFloat
{
    fn rssq(&'a self) -> Self::RowsMapped<T::Real>;
}

impl<'a, T, L> Rssq<'a, T> for L
where
    L: Lists<T, RowView<'a>: List<T>> + 'a,
    T: ComplexFloat + 'a
{
    fn rssq(&'a self) -> Self::RowsMapped<T::Real>
    {
        self.map_rows_to_owned(|x| {
            let x = x.to_vec();
            x.into_iter()
                .map(|x| (x.conj()*x).re())
                .reduce(Add::add)
                .map(|x| Float::sqrt(x))
                .unwrap_or(T::Real::zero())
        })
    }
}