use crate::traits::{
math::{Normed, Norm}
};
use crate::structure::matrix::Matrix;
use crate::numerical::{
integral,
integral::Integral::GaussLegendre,
eigen,
eigen::{Eigen, EigenMethod::Jacobi},
};
pub trait SimpleNorm: Normed {
fn norm(&self) -> Self::Scalar;
fn normalize(&self) -> Self;
}
pub fn integrate<F: Fn(f64) -> f64>(f: F, (a, b): (f64, f64)) -> f64 {
integral::integrate(f, (a, b), GaussLegendre(15))
}
pub fn eigen(m: &Matrix) -> Eigen {
eigen::eigen(m, Jacobi)
}
impl SimpleNorm for Vec<f64> {
fn norm(&self) -> Self::Scalar {
Normed::norm(self, Norm::L2)
}
fn normalize(&self) -> Self {
Normed::normalize(self, Norm::L2)
}
}
impl SimpleNorm for Matrix {
fn norm(&self) -> Self::Scalar {
Normed::norm(self, Norm::F)
}
fn normalize(&self) -> Self {
unimplemented!()
}
}