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

use crate::algebra::linear::matrix::{General, UpperTriangular};
#[cfg(feature = "serde")]
use serde::{Deserialize, Serialize};
use std::clone::Clone;

/// Result of Schur decomposition
#[cfg_attr(feature = "serde", derive(Serialize, Deserialize))]
#[derive(Debug, Clone)]
pub struct SchurDec<T> {
    q: General<T>,
    u: UpperTriangular<T>,
}

///
/// ```math
/// A = Q * U Q ^-1
/// ```
impl<T> SchurDec<T> {
    pub(super) fn new(q: General<T>, u: UpperTriangular<T>) -> SchurDec<T> {
        SchurDec { q, u }
    }

    /// Return the unitary matrix q
    ///
    /// # Arguments
    ///
    /// * `self`
    pub fn q(self) -> General<T> {
        self.q
    }

    /// Return the upper triangular matrix u
    pub fn u(self) -> UpperTriangular<T> {
        self.u
    }

    pub fn qu(self) -> (General<T>, UpperTriangular<T>) {
        (self.q, self.u)
    }
}