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
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
use num_complex::Complex;
use num_traits::Float;
use crate::Vector;
pub trait VAbs: Vector
{
type Output;
fn abs_sqr(&self) -> Self::Output;
fn abs(&self) -> Self::Output;
}
impl<const N: usize> VAbs for [f32; N]
{
type Output = f32;
fn abs_sqr(&self) -> Self::Output
{
self.iter().map(|x| x.abs().powi(2)).reduce(|a, b| a + b).unwrap_or(0.0)
}
fn abs(&self) -> Self::Output
{
self.abs_sqr().sqrt()
}
}
impl<const N: usize> VAbs for [f64; N]
{
type Output = f64;
fn abs_sqr(&self) -> Self::Output
{
self.iter().map(|x| x.abs().powi(2)).reduce(|a, b| a + b).unwrap_or(0.0)
}
fn abs(&self) -> Self::Output
{
self.abs_sqr().sqrt()
}
}
impl<F: Float, const N: usize> VAbs for [Complex<F>; N]
{
type Output = F;
fn abs_sqr(&self) -> Self::Output
{
self.iter().map(|x| x.norm_sqr()).reduce(|a, b| a + b).unwrap_or(F::zero())
}
fn abs(&self) -> Self::Output
{
self.abs_sqr().sqrt()
}
}