geomcore 0.34.0

A pure, standalone geometric kernel for CAD-grade curves and surfaces: parametric evaluation, transformations, and curve-on-surface parametrization.
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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89

//! Boxes of surfaces: closed form for spheres and tori, control points for
//! B-splines; planes, cylinders and cones are unbounded.

use super::{Bound3D, BoundingBox3D, coordinates, point};
use crate::Point3D;
use crate::surfaces::{
    BSplineSurface, Cone, Cylinder, ElementarySurface, Plane, Sphere, Surface, Torus,
};

/// The box `center ± half` on each axis.
fn centered(center: Point3D, half: [f64; 3]) -> BoundingBox3D {
    let c = coordinates(center);
    BoundingBox3D {
        min: point(std::array::from_fn(|k| c[k] - half[k])),
        max: point(std::array::from_fn(|k| c[k] + half[k])),
    }
}

impl Bound3D for Plane {
    fn bounding_box(&self) -> Option<BoundingBox3D> {
        None
    }
}

impl Bound3D for Cylinder {
    fn bounding_box(&self) -> Option<BoundingBox3D> {
        None
    }
}

impl Bound3D for Cone {
    fn bounding_box(&self) -> Option<BoundingBox3D> {
        None
    }
}

impl Bound3D for Sphere {
    fn bounding_box(&self) -> Option<BoundingBox3D> {
        Some(centered(self.center(), [self.radius(); 3]))
    }
}

impl Bound3D for Torus {
    /// Along a unit vector `e`, the torus reaches `R h + r` from its
    /// center, `R` and `r` its major and minor radii and `h` the length of
    /// the projection of `e` onto the torus's plane: the center circle
    /// reaches `R h`, and the tube adds `r` around it.
    fn bounding_box(&self) -> Option<BoundingBox3D> {
        let frame = self.frame();
        let (x, y) = (frame.x_direction(), frame.y_direction());
        let (x, y) = ([x.x, x.y, x.z], [y.x, y.y, y.z]);
        let (major, minor) = (self.major_radius(), self.minor_radius());
        Some(centered(
            self.center(),
            std::array::from_fn(|k| major * x[k].hypot(y[k]) + minor),
        ))
    }
}

impl Bound3D for BSplineSurface {
    fn bounding_box(&self) -> Option<BoundingBox3D> {
        BoundingBox3D::from_points(self.poles().iter().copied())
    }
}

impl Bound3D for Surface {
    fn bounding_box(&self) -> Option<BoundingBox3D> {
        match self {
            Surface::Plane(s) => s.bounding_box(),
            Surface::Cylinder(s) => s.bounding_box(),
            Surface::Cone(s) => s.bounding_box(),
            Surface::Sphere(s) => s.bounding_box(),
            Surface::Torus(s) => s.bounding_box(),
            Surface::BSpline(s) => s.bounding_box(),
        }
    }
}

impl Bound3D for ElementarySurface {
    fn bounding_box(&self) -> Option<BoundingBox3D> {
        match self {
            ElementarySurface::Plane(s) => s.bounding_box(),
            ElementarySurface::Cylinder(s) => s.bounding_box(),
            ElementarySurface::Cone(s) => s.bounding_box(),
            ElementarySurface::Sphere(s) => s.bounding_box(),
            ElementarySurface::Torus(s) => s.bounding_box(),
        }
    }
}