use std::collections::VecDeque;
use num_traits::{Num, Float};
use crate::{Norm, Distance};
use crate::desc::{Abs, Sup, PNorm, PNormReal};
use crate::utility::{supnorm_iterable, pnorm_iterable, pnorm_real_iterable};
impl<T: Norm<Abs, Output = R>, R: Num + PartialOrd> Norm<Sup> for Vec<T> {
type Output = <T as Norm<Abs>>::Output;
fn norm(&self, _desc: Sup) -> <Self as Norm<Sup>>::Output {
supnorm_iterable(self.iter().map(|a| a.norm(Abs::new())))
}
}
impl<T: Norm<Abs, Output = R>, R: Float + From<f32>> Norm<PNorm> for Vec<T> {
type Output = <T as Norm<Abs>>::Output;
fn norm(&self, desc: PNorm) -> <Self as Norm<PNorm>>::Output {
pnorm_iterable(self.iter().map(|a| a.norm(Abs::new())), desc)
}
}
impl<T: Norm<Abs, Output = R>, R: Float + From<f32>> Norm<PNormReal> for Vec<T> {
type Output = <T as Norm<Abs>>::Output;
fn norm(&self, desc: PNormReal) -> <Self as Norm<PNormReal>>::Output {
pnorm_real_iterable(self.iter().map(|a| a.norm(Abs::new())), desc)
}
}
impl<T: Distance<Abs, Output = R>, R: Num + PartialOrd> Distance<Sup> for Vec<T> {
type Output = <T as Distance<Abs>>::Output;
fn distance(&self, other: &Self, _desc: Sup) -> <Self as Distance<Sup>>::Output {
assert_eq!(
self.len(), other.len(),
"Sizes of Vecs do not match. Left: {}; Right: {}.",
self.len(), other.len()
);
supnorm_iterable(self.iter().zip(other).map(|(a, b)| a.distance(b, Abs::new())))
}
}
impl<T: Distance<Abs, Output = R>, R: Float + From<f32>> Distance<PNorm> for Vec<T> {
type Output = <T as Distance<Abs>>::Output;
fn distance(&self, other: &Self, desc: PNorm) -> <Self as Distance<PNorm>>::Output {
assert_eq!(
self.len(), other.len(),
"Sizes of Vecs do not match. Left: {}; Right: {}.",
self.len(), other.len()
);
pnorm_iterable(self.iter().zip(other).map(|(a, b)| a.distance(b, Abs::new())), desc)
}
}
impl<T: Distance<Abs, Output = R>, R: Float + From<f32>> Distance<PNormReal> for Vec<T> {
type Output = <T as Distance<Abs>>::Output;
fn distance(&self, other: &Self, desc: PNormReal) -> <Self as Distance<PNormReal>>::Output {
assert_eq!(
self.len(), other.len(),
"Sizes of Vecs do not match. Left: {}; Right: {}.",
self.len(), other.len()
);
pnorm_real_iterable(self.iter().zip(other).map(|(a, b)| a.distance(b, Abs::new())), desc)
}
}
impl<T: Norm<Abs, Output = R>, R: Num + PartialOrd> Norm<Sup> for VecDeque<T> {
type Output = <T as Norm<Abs>>::Output;
fn norm(&self, _desc: Sup) -> <Self as Norm<Sup>>::Output {
supnorm_iterable(self.iter().map(|a| a.norm(Abs::new())))
}
}
impl<T: Norm<Abs, Output = R>, R: Float + From<f32>> Norm<PNorm> for VecDeque<T> {
type Output = <T as Norm<Abs>>::Output;
fn norm(&self, desc: PNorm) -> <Self as Norm<PNorm>>::Output {
pnorm_iterable(self.iter().map(|a| a.norm(Abs::new())), desc)
}
}
impl<T: Norm<Abs, Output = R>, R: Float + From<f32>> Norm<PNormReal> for VecDeque<T> {
type Output = <T as Norm<Abs>>::Output;
fn norm(&self, desc: PNormReal) -> <Self as Norm<PNormReal>>::Output {
pnorm_real_iterable(self.iter().map(|a| a.norm(Abs::new())), desc)
}
}
impl<T: Distance<Abs, Output = R>, R: Num + PartialOrd> Distance<Sup> for VecDeque<T> {
type Output = <T as Distance<Abs>>::Output;
fn distance(&self, other: &Self, _desc: Sup) -> <Self as Distance<Sup>>::Output {
assert_eq!(
self.len(), other.len(),
"Sizes of VecDeques do not match. Left: {}; Right: {}.",
self.len(), other.len()
);
supnorm_iterable(self.iter().zip(other).map(|(a, b)| a.distance(b, Abs::new())))
}
}
impl<T: Distance<Abs, Output = R>, R: Float + From<f32>> Distance<PNorm> for VecDeque<T> {
type Output = <T as Distance<Abs>>::Output;
fn distance(&self, other: &Self, desc: PNorm) -> <Self as Distance<PNorm>>::Output {
assert_eq!(
self.len(), other.len(),
"Sizes of VecDeques do not match. Left: {}; Right: {}.",
self.len(), other.len()
);
pnorm_iterable(self.iter().zip(other).map(|(a, b)| a.distance(b, Abs::new())), desc)
}
}
impl<T: Distance<Abs, Output = R>, R: Float + From<f32>> Distance<PNormReal> for VecDeque<T> {
type Output = <T as Distance<Abs>>::Output;
fn distance(&self, other: &Self, desc: PNormReal) -> <Self as Distance<PNormReal>>::Output {
assert_eq!(
self.len(), other.len(),
"Sizes of VecDeques do not match. Left: {}; Right: {}.",
self.len(), other.len()
);
pnorm_real_iterable(self.iter().zip(other).map(|(a, b)| a.distance(b, Abs::new())), desc)
}
}
#[cfg(test)]
mod tests {
use std::collections::VecDeque;
use crate::{Norm, Distance};
use crate::desc::{Sup, PNorm, PNormReal};
#[test]
fn supnorm_vec() {
let a = vec![3.0f32, -4.0, 2.0];
assert_eq!(a.norm(Sup::new()), 4.0);
let b = vec![2.0f32, 2.0, 2.0];
assert_eq!(a.distance(&b, Sup::new()), 6.0);
}
#[test]
fn pnorm_vec() {
let a = vec![3.0f32, -4.0, 2.0];
assert_eq!(a.norm(PNorm::new(2)), 29.0f32.sqrt());
let b = vec![2.0f32, 2.0, 2.0];
assert_eq!(a.distance(&b, PNorm::new(2)), 37.0f32.sqrt());
}
#[test]
fn pnorm_real_vec() {
let a = vec![3.0f32, -4.0, 2.0];
assert_eq!(a.norm(PNormReal::from_f32(2.0)), 29.0f32.sqrt());
let b = vec![2.0f32, 2.0, 2.0];
assert_eq!(a.distance(&b, PNormReal::from_f32(2.0)), 37.0f32.sqrt());
}
#[test]
fn supnorm_vec_deque() {
let a = vec![3.0f32, -4.0, 2.0];
assert_eq!(a.norm(Sup::new()), 4.0);
let b = vec![2.0f32, 2.0, 2.0];
assert_eq!(a.distance(&b, Sup::new()), 6.0);
}
#[test]
fn pnorm_vec_deque() {
let a = VecDeque::from(vec![3.0f32, -4.0, 2.0]);
assert_eq!(a.norm(PNorm::new(2)), 29.0f32.sqrt());
let b = VecDeque::from(vec![2.0f32, 2.0, 2.0]);
assert_eq!(a.distance(&b, PNorm::new(2)), 37.0f32.sqrt());
}
#[test]
fn pnorm_real_vec_deque() {
let a = VecDeque::from(vec![3.0f32, -4.0, 2.0]);
assert_eq!(a.norm(PNormReal::from_f32(2.0)), 29.0f32.sqrt());
let b = VecDeque::from(vec![2.0f32, 2.0, 2.0]);
assert_eq!(a.distance(&b, PNormReal::from_f32(2.0)), 37.0f32.sqrt());
}
}