use crate::algebra::linear::matrix::{General, UpperTriangular};
#[cfg(feature = "serde")]
use serde::{Deserialize, Serialize};
use std::clone::Clone;
pub trait QRDecomposition<T> {
fn dec_qr(&self) -> Result<QRDec<T>, ()>;
}
#[cfg_attr(feature = "serde", derive(Serialize, Deserialize))]
#[derive(Debug, Clone)]
pub struct QRDec<T> {
q: General<T>,
r: UpperTriangular<T>,
}
impl<T> QRDec<T> {
pub(super) fn new(q: General<T>, r: UpperTriangular<T>) -> QRDec<T> {
QRDec { q, r }
}
pub fn q(self) -> General<T> {
self.q
}
pub fn r(self) -> UpperTriangular<T> {
self.r
}
pub fn qr(self) -> (General<T>, UpperTriangular<T>) {
(self.q, self.r)
}
}