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
use std::cmp::Ordering;
use std::f32;
pub struct BoundedFunction {
pub func: fn(f32) -> f32,
pub lower: f32,
pub higher: f32,
}
pub struct PartialFunction {
funcs: Vec<BoundedFunction>,
}
impl PartialFunction {
pub fn new() -> PartialFunctionBuilder {
PartialFunctionBuilder::new()
}
pub fn eval(&self, x: f32) -> f32 {
let iter = self.funcs.iter().enumerate();
for (i, bounded) in iter {
let next = self.funcs.get(i + 1);
if (x >= bounded.lower && x < bounded.higher) ||
(next.is_none() && x == bounded.higher) ||
(next.is_some() && next.unwrap().lower != bounded.higher)
{
let f = bounded.func;
return f(x);
}
}
f32::NAN
}
}
pub struct PartialFunctionBuilder {
funcs: Vec<BoundedFunction>,
}
impl PartialFunctionBuilder {
pub fn new() -> Self {
PartialFunctionBuilder { funcs: vec![] }
}
pub fn with(mut self, lower: f32, higher: f32, func: fn(f32) -> f32) -> Self {
let f = BoundedFunction {
func: func,
lower: lower,
higher: higher,
};
debug_assert!(self.can_insert(lower, higher));
self.funcs.push(f);
self
}
pub fn can_insert(&self, lower: f32, higher: f32) -> bool {
!self.funcs.iter().any(|b| {
(lower >= b.lower && lower < b.higher) || (higher > b.lower && higher <= b.higher) ||
(lower <= b.lower && higher >= b.higher)
})
}
pub fn build(mut self) -> PartialFunction {
self.funcs.sort_by(|a, b| {
a.lower
.partial_cmp(&b.lower)
.unwrap_or(a.higher.partial_cmp(&b.higher).unwrap_or(Ordering::Equal))
});
PartialFunction { funcs: self.funcs }
}
}