Fri, 13 Dec 2024 22:37:12 -0500
Add direct products (Pair, RowOp, ColOp, DiagOp)
5 | 1 | /*! |
2 | Some useful (numerical) types and traits. | |
3 | ||
4 | The traits are based on corresponding ones in [`num_traits`], but try to fill some gaps in the | |
5 | super-traits and available constants. | |
0 | 6 | |
5 | 7 | As [`nalgebra`] unnecessarily provides many of the same methods as [`num_traits`], to avoid having |
8 | to refer to the methods with the full path, it is often necesary to use [`ToNalgebraRealField`][crate::nalgebra_support::ToNalgebraRealField] to hide the nalgebra implementations until | |
9 | absolutely necessary to use nalgebra. | |
10 | */ | |
11 | ||
12 | //use trait_set::trait_set; | |
0 | 13 | pub use num_traits::Float as NumTraitsFloat; // needed to re-export functions. |
14 | pub use num_traits::cast::AsPrimitive; | |
15 | ||
16 | /// Typical integer type | |
17 | #[allow(non_camel_case_types)] | |
18 | pub type int = i64; | |
19 | ||
20 | /// Typical unsigned integer type | |
21 | #[allow(non_camel_case_types)] | |
22 | pub type uint = u64; | |
23 | ||
24 | /// Typical floating point number type | |
25 | #[allow(non_camel_case_types)] | |
26 | pub type float = f64; | |
27 | ||
28 | /// Casts of abstract numerical types to others via the standard `as` keyword. | |
29 | pub trait CastFrom<T : 'static + Copy> : num_traits::cast::AsPrimitive<T> { | |
30 | fn cast_from(other : T) -> Self; | |
31 | } | |
32 | ||
33 | macro_rules! impl_casts { | |
34 | ($($type:ty)*) => { $( | |
35 | impl_casts!(@phase2, $type, | |
36 | u8 u16 u32 u64 u128 usize | |
37 | i8 i16 i32 i64 i128 isize | |
38 | f32 f64); | |
39 | )* }; | |
40 | (@phase2, $type:ty, $($type2:ty)*) => { $( | |
41 | impl CastFrom<$type2> for $type { | |
42 | #[inline] | |
43 | fn cast_from(other : $type2) -> Self { other as $type } | |
44 | } | |
45 | )* }; | |
46 | } | |
47 | ||
48 | impl_casts!(u8 u16 u32 u64 u128 usize | |
49 | i8 i16 i32 i64 i128 isize | |
50 | f32 f64); | |
51 | ||
5 | 52 | /// Trait for general numeric types |
8
4e09b7829b51
Multithreaded bisection tree operations
Tuomo Valkonen <tuomov@iki.fi>
parents:
5
diff
changeset
|
53 | pub trait Num : 'static + Copy + Sync + Send + num::Num + num_traits::NumAssign |
0 | 54 | + std::iter::Sum + std::iter::Product |
55 | + std::fmt::Debug + std::fmt::Display + serde::Serialize | |
56 | + CastFrom<u8> + CastFrom<u16> + CastFrom<u32> + CastFrom<u64> | |
57 | + CastFrom<u128> + CastFrom<usize> | |
58 | + CastFrom<i8> + CastFrom<i16> + CastFrom<i32> + CastFrom<i64> | |
59 | + CastFrom<i128> + CastFrom<isize> | |
60 | + CastFrom<f32> + CastFrom<f64> { | |
61 | ||
62 | const ZERO : Self; | |
63 | const ONE : Self; | |
64 | const TWO : Self; | |
8
4e09b7829b51
Multithreaded bisection tree operations
Tuomo Valkonen <tuomov@iki.fi>
parents:
5
diff
changeset
|
65 | /// Generic version of `Self::MAX` |
0 | 66 | const RANGE_MAX : Self; |
8
4e09b7829b51
Multithreaded bisection tree operations
Tuomo Valkonen <tuomov@iki.fi>
parents:
5
diff
changeset
|
67 | /// Generic version of `Self::MIN` |
0 | 68 | const RANGE_MIN : Self; |
69 | } | |
70 | ||
71 | /// Trait for signed numeric types | |
72 | pub trait SignedNum : Num + num::Signed + std::ops::Neg<Output=Self> {} | |
73 | impl<U : Num + num::Signed + std::ops::Neg<Output=Self>> SignedNum for U { } | |
74 | ||
75 | /// Trait for floating point numbers | |
5 | 76 | pub trait Float : SignedNum + num::Float /*+ From<Self::CompatibleSize>*/ { |
77 | // An unsigned integer that can be used for indexing operations and | |
78 | // converted to F without loss. | |
79 | //type CompatibleSize : CompatibleUnsigned<Self>; | |
0 | 80 | |
81 | const PI : Self; | |
82 | const E : Self; | |
83 | const EPSILON : Self; | |
84 | const SQRT_2 : Self; | |
85 | const INFINITY : Self; | |
86 | const NEG_INFINITY : Self; | |
32 | 87 | const NAN : Self; |
0 | 88 | const FRAC_2_SQRT_PI : Self; |
89 | } | |
90 | ||
91 | /// Trait for integers | |
92 | pub trait Integer : Num + num::Integer {} | |
93 | ||
94 | /// Trait for unsigned integers | |
95 | pub trait Unsigned : Num + Integer + num::Unsigned {} | |
96 | ||
97 | /// Trait for signed integers | |
98 | pub trait Signed : SignedNum + Integer {} | |
99 | ||
100 | macro_rules! impl_num_consts { | |
101 | ($($type:ty)*) => { $( | |
102 | impl Num for $type { | |
103 | const ZERO : Self = 0 as $type; | |
104 | const ONE : Self = 1 as $type; | |
105 | const TWO : Self = 2 as $type; | |
106 | const RANGE_MAX : Self = <$type>::MAX; | |
107 | const RANGE_MIN : Self = <$type>::MIN; | |
108 | } | |
109 | )* } | |
110 | } | |
111 | ||
112 | macro_rules! impl_integers { | |
113 | ($signed:ty : $($type:ty)*) => { $( | |
114 | impl_num_consts!($type); | |
115 | impl Integer for $type {} | |
116 | impl $signed for $type {} | |
117 | )* } | |
118 | } | |
119 | ||
120 | impl_integers!(Signed: i8 i16 i32 i64 i128 isize); | |
121 | impl_integers!(Unsigned: u8 u16 u32 u64 u128 usize); | |
122 | ||
123 | impl_num_consts!(f32 f64); | |
124 | ||
125 | impl Float for f64 { | |
5 | 126 | /*#[cfg(any(target_pointer_width = "128", target_pointer_width = "64"))] |
0 | 127 | type CompatibleSize = u32; |
128 | #[cfg(any(target_pointer_width = "32", target_pointer_width = "16"))] | |
5 | 129 | type CompatibleSize = usize;*/ |
0 | 130 | |
131 | const PI : Self = std::f64::consts::PI; | |
132 | const E : Self = std::f64::consts::E; | |
133 | const EPSILON : Self = std::f64::EPSILON; | |
134 | const SQRT_2 : Self = std::f64::consts::SQRT_2; | |
135 | const INFINITY : Self = std::f64::INFINITY; | |
136 | const NEG_INFINITY : Self = std::f64::NEG_INFINITY; | |
32 | 137 | const NAN : Self = std::f64::NAN; |
0 | 138 | const FRAC_2_SQRT_PI : Self = std::f64::consts::FRAC_2_SQRT_PI; |
139 | } | |
140 | ||
141 | impl Float for f32 { | |
5 | 142 | /* |
0 | 143 | #[cfg(any(target_pointer_width = "128", target_pointer_width = "64", target_pointer_width = "32"))] |
144 | type CompatibleSize = u16; | |
145 | #[cfg(any(target_pointer_width = "16"))] | |
146 | type CompatibleSize = usize; | |
5 | 147 | */ |
0 | 148 | |
149 | const PI : Self = std::f32::consts::PI; | |
150 | const E : Self = std::f32::consts::E; | |
151 | const EPSILON : Self = std::f32::EPSILON; | |
152 | const SQRT_2 : Self = std::f32::consts::SQRT_2; | |
153 | const INFINITY : Self = std::f32::INFINITY; | |
154 | const NEG_INFINITY : Self = std::f32::NEG_INFINITY; | |
32 | 155 | const NAN : Self = std::f32::NAN; |
0 | 156 | const FRAC_2_SQRT_PI : Self = std::f32::consts::FRAC_2_SQRT_PI; |
157 | } | |
158 | ||
5 | 159 | /* |
0 | 160 | trait_set! { |
161 | pub trait CompatibleUnsigned<F : Float> = Unsigned + Into<F>; | |
162 | pub trait CompatibleSigned<F : Float> = Signed + Into<F>; | |
163 | } | |
5 | 164 | */ |
13
465fa2121ccb
Better Linear and Mapping structure that can provide consuming and reference `apply`.
Tuomo Valkonen <tuomov@iki.fi>
parents:
8
diff
changeset
|
165 |