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

use crate::circulant_matrix::CirculantMatrix;
use crate::number::Number;
pub struct ToeplitzMatrix<U>
where
    U: Number,
{
    dim: usize,
    row: Vec<U>,
    column: Vec<U>,
}

impl<U> ToeplitzMatrix<U>
where
    U: Number,
{
    pub fn new(row: Vec<U>, column: Vec<U>) -> Self {
        let dim = row.len();

        if column.len() != dim {
            panic!("dimension mismatch")
        }
        if dim < 2 || row[0] != column[0] {
            panic!("")
        }

        Self { dim, row, column }
    }

    pub fn get_dim(&self) -> usize {
        self.dim
    }

    pub fn get_row(&self) -> &[U] {
        &self.row
    }

    pub fn get_column(&self) -> &[U] {
        &self.column
    }

    pub fn embedded_circulant(&self) -> CirculantMatrix<U> {
        let row = (0..self.dim)
            .into_iter()
            .chain((1..self.dim - 1).rev().into_iter())
            .map(|i| self.row[i])
            .collect();

        CirculantMatrix::<U>::new(row)
    }
}