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
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
pub use std::ops::{Add, Sub, Mul, Div, Rem, Neg, AddAssign, SubAssign, MulAssign, DivAssign};
pub trait Operator: Copy {
fn operator_token() -> Self;
}
pub trait Inverse<O: Operator>: Sized {
fn inverse(&self) -> Self;
#[inline]
fn inverse_mut(&mut self) {
*self = self.inverse()
}
}
#[derive(Clone, Copy)]
pub struct Additive;
#[derive(Clone, Copy)]
pub struct Multiplicative;
#[derive(Clone, Copy)]
pub struct AbstractOperator;
impl Operator for Additive {
#[inline]
fn operator_token() -> Self {
Additive
}
}
impl Operator for Multiplicative {
#[inline]
fn operator_token() -> Self {
Multiplicative
}
}
impl Operator for AbstractOperator {
#[inline]
fn operator_token() -> Self {
AbstractOperator
}
}
macro_rules! impl_additive_inverse(
($($T:ty),* $(,)*) => {$(
impl Inverse<Additive> for $T {
fn inverse(&self) -> Self {
-*self
}
}
)*}
);
impl_additive_inverse!(i8, i16, i32, i64, f32, f64);
impl Inverse<Multiplicative> for f32 {
#[inline]
fn inverse(&self) -> f32 {
1.0 / self
}
}
impl Inverse<Multiplicative> for f64 {
#[inline]
fn inverse(&self) -> f64 {
1.0 / self
}
}
pub trait ClosedAdd<Right = Self>: Sized + Add<Right, Output = Self> + AddAssign<Right> { }
pub trait ClosedSub<Right = Self>: Sized + Sub<Right, Output = Self> + SubAssign<Right> { }
pub trait ClosedMul<Right = Self>: Sized + Mul<Right, Output = Self> + MulAssign<Right> { }
pub trait ClosedDiv<Right = Self>: Sized + Div<Right, Output = Self> + DivAssign<Right> { }
pub trait ClosedNeg: Sized + Neg<Output = Self> { }
impl<T, Right> ClosedAdd<Right> for T where T: Add<Right, Output = T> + AddAssign<Right> { }
impl<T, Right> ClosedSub<Right> for T where T: Sub<Right, Output = T> + SubAssign<Right> { }
impl<T, Right> ClosedMul<Right> for T where T: Mul<Right, Output = T> + MulAssign<Right> { }
impl<T, Right> ClosedDiv<Right> for T where T: Div<Right, Output = T> + DivAssign<Right> { }
impl<T> ClosedNeg for T where T: Neg<Output = T> { }