curvo 0.1.88

NURBS modeling library
Documentation 
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

use crate::misc::FloatingPoint;

mod continuous;
mod curve;

/// The type of discontinuity.
#[derive(Debug, Clone, Copy, PartialEq, Eq, Hash)]
pub enum DiscontinuityType {
    /// C1 (1st derivative)
    C1,
    /// C2 (2nd derivative)
    C2,
    /// G2 (continuous unit tangent & curvature)
    G2,
}

/// Search for a derivative, tangent, or curvature discontinuity.
pub trait Discontinuity<T> {
    /// Get the next discontinuity.
    fn get_next_discontinuity(&self, ty: DiscontinuityType, start: T, end: T) -> Option<T>;

    /// Get an iterator for searching for discontinuities.
    fn discontinuity_iter(
        &self,
        ty: DiscontinuityType,
        start: T,
        end: T,
    ) -> DiscontinuityIterator<'_, Self, T>
    where
        Self: Sized,
    {
        DiscontinuityIterator {
            geometry: self,
            ty,
            next: start,
            end,
        }
    }
}

/// Iterator for searching for discontinuities.
pub struct DiscontinuityIterator<'a, G, T>
where
    G: Discontinuity<T>,
{
    geometry: &'a G,
    ty: DiscontinuityType,
    next: T,
    end: T,
}

impl<G, T> Iterator for DiscontinuityIterator<'_, G, T>
where
    G: Discontinuity<T>,
    T: FloatingPoint,
{
    type Item = T;

    /// Get the next discontinuity.
    fn next(&mut self) -> Option<Self::Item> {
        let item = self
            .geometry
            .get_next_discontinuity(self.ty, self.next, self.end);
        if let Some(t) = item {
            let h = (self.end - self.next) * T::from_f64(1e-12).unwrap();
            self.next = t + h;
            // println!("next: {}, end: {}", self.next, self.end);
            Some(t)
        } else {
            None
        }
    }
}