axiolid-reference 0.3.0

Portable scalar reference implementation and certified predicates
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
90
91
92
93
94
95
96

use axiolid_core::Point3;
use axiolid_reference::{segment_triangle_relation, SegmentTriangleRelation};

const TRIANGLE: [Point3; 3] = [
    Point3::new(0.0, 0.0, 0.0),
    Point3::new(2.0, 0.0, 0.0),
    Point3::new(0.0, 2.0, 0.0),
];

#[test]
fn proper_crossing_is_certified() {
    assert_eq!(
        segment_triangle_relation(
            Point3::new(0.5, 0.5, -1.0),
            Point3::new(0.5, 0.5, 1.0),
            TRIANGLE,
        ),
        SegmentTriangleRelation::Proper,
    );
}

#[test]
fn edge_and_vertex_contact_are_touching() {
    assert_eq!(
        segment_triangle_relation(
            Point3::new(0.0, 0.0, -1.0),
            Point3::new(0.0, 0.0, 1.0),
            TRIANGLE,
        ),
        SegmentTriangleRelation::Touching,
    );
    assert_eq!(
        segment_triangle_relation(
            Point3::new(1.0, 1.0, -1.0),
            Point3::new(1.0, 1.0, 1.0),
            TRIANGLE,
        ),
        SegmentTriangleRelation::Touching,
    );
}

#[test]
fn parallel_and_coplanar_segments_are_distinct() {
    assert_eq!(
        segment_triangle_relation(
            Point3::new(0.5, 0.5, 1.0),
            Point3::new(1.0, 0.5, 1.0),
            TRIANGLE,
        ),
        SegmentTriangleRelation::Disjoint,
    );
    assert_eq!(
        segment_triangle_relation(
            Point3::new(0.25, 0.25, 0.0),
            Point3::new(1.0, 0.25, 0.0),
            TRIANGLE,
        ),
        SegmentTriangleRelation::Coplanar,
    );
}

#[test]
fn plane_crossing_outside_the_triangle_is_disjoint() {
    assert_eq!(
        segment_triangle_relation(
            Point3::new(2.0, 2.0, -1.0),
            Point3::new(2.0, 2.0, 1.0),
            TRIANGLE,
        ),
        SegmentTriangleRelation::Disjoint,
    );
    assert_eq!(
        segment_triangle_relation(
            Point3::new(2.0, 2.0, 1.0),
            Point3::new(2.0, 2.0, 0.0),
            TRIANGLE,
        ),
        SegmentTriangleRelation::Disjoint,
    );
}

#[test]
fn zero_length_or_collinear_inputs_are_explicit() {
    assert_eq!(
        segment_triangle_relation(Point3::ZERO, Point3::ZERO, TRIANGLE),
        SegmentTriangleRelation::DegenerateSegment,
    );
    assert_eq!(
        segment_triangle_relation(
            Point3::new(0.0, 0.0, -1.0),
            Point3::new(0.0, 0.0, 1.0),
            [Point3::ZERO, Point3::X, Point3::new(2.0, 0.0, 0.0)],
        ),
        SegmentTriangleRelation::DegenerateTriangle,
    );
}