yagi 0.1.0

Batteries-included DSP 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
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
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117

use num_complex::{Complex, ComplexFloat};
use num_traits::{Float, Zero};

// Extension trait for complex number support
pub trait IsComplex {
    fn is_complex() -> bool;
}

impl IsComplex for f32 {
    fn is_complex() -> bool {
        false
    }
}

impl IsComplex for f64 {
    fn is_complex() -> bool {
        false
    }
}

impl IsComplex for Complex<f32> {
    fn is_complex() -> bool {
        true
    }
}

impl IsComplex for Complex<f64> {
    fn is_complex() -> bool {
        true
    }
}

pub trait FloatComplex:
    ComplexFloat<Real: Float + Zero + From<f32> + std::fmt::LowerExp>
    + std::fmt::Display
    + std::ops::AddAssign
    + std::ops::SubAssign
    + std::ops::MulAssign
    + Default
    + std::iter::Sum
    + IsComplex
{
}

impl<
        T: ComplexFloat<Real: Float + Zero + From<f32> + std::fmt::LowerExp>
            + std::fmt::Display
            + std::ops::AddAssign
            + std::ops::SubAssign
            + std::ops::MulAssign
            + Default
            + std::iter::Sum
            + IsComplex,
    > FloatComplex for T
{
}

/// Print matrix to stdout
pub fn matrix_print<T>(x: &[T], rows: usize, cols: usize)
where
    T: FloatComplex,
{
    println!("matrix [{} x {}] :", rows, cols);
    for r in 0..rows {
        for c in 0..cols {
            print!("{:12}", x[r * cols + c]);
        }
        println!();
    }
}

/// Initialize square matrix to the identity matrix
pub fn matrix_eye<T>(x: &mut [T], n: usize)
where
    T: FloatComplex,
{
    for r in 0..n {
        for c in 0..n {
            x[r * n + c] = if r == c { T::one() } else { T::zero() };
        }
    }
}

/// Initialize matrix to ones
pub fn matrix_ones<T>(x: &mut [T], rows: usize, cols: usize)
where
    T: FloatComplex,
{
    for i in 0..(rows * cols) {
        x[i] = T::one();
    }
}

/// Initialize matrix to zeros
pub fn matrix_zeros<T>(x: &mut [T], rows: usize, cols: usize)
where
    T: FloatComplex,
{
    for i in 0..(rows * cols) {
        x[i] = T::zero();
    }
}

// Helper function to access matrix elements
#[inline]
pub fn matrix_access<T>(matrix: &[T], _rows: usize, cols: usize, r: usize, c: usize) -> T
where
    T: FloatComplex,
{
    matrix[r * cols + c]
}

// Helper function to mutably access matrix elements
#[inline]
pub fn matrix_access_mut<T>(matrix: &mut [T], _rows: usize, cols: usize, r: usize, c: usize, value: T) {
    matrix[r * cols + c] = value;
}