runmat-meshing-core 0.6.0

Deterministic meshing preparation contracts for RunMat
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
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127

use super::*;
use crate::quality::tolerance::MeshingTolerance;

#[test]
fn tetrahedron_signed_volume_reports_orientation() {
    let tetrahedron = [
        [0.0, 0.0, 0.0],
        [1.0, 0.0, 0.0],
        [0.0, 1.0, 0.0],
        [0.0, 0.0, 1.0],
    ];

    assert!((tetrahedron_signed_volume(tetrahedron) - 1.0 / 6.0).abs() < 1.0e-12);
    let (node_ids, volume) = orient_tetrahedron_node_ids(
        [1, 3, 2, 4],
        [
            tetrahedron[0],
            tetrahedron[2],
            tetrahedron[1],
            tetrahedron[3],
        ],
    );
    assert_eq!(node_ids, [1, 2, 3, 4]);
    assert!((volume - 1.0 / 6.0).abs() < 1.0e-12);
}

#[test]
fn tetrahedron_scaled_jacobian_reports_shape_quality() {
    let orthogonal = [
        [0.0, 0.0, 0.0],
        [1.0, 0.0, 0.0],
        [0.0, 1.0, 0.0],
        [0.0, 0.0, 1.0],
    ];
    let sliver = [
        [0.0, 0.0, 0.0],
        [1.0, 0.0, 0.0],
        [0.0, 1.0, 0.0],
        [1.0e-4, 1.0e-4, 1.0e-5],
    ];

    assert!(tetrahedron_scaled_jacobian(orthogonal) > 0.5);
    assert!(tetrahedron_scaled_jacobian(sliver) < 0.01);
}

#[test]
fn ray_triangle_intersection_returns_distance_and_barycentric_coordinates() {
    let hit = ray_triangle_intersection(
        [-1.0, 0.25, 0.25],
        [1.0, 0.0, 0.0],
        [[0.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0]],
        MeshingTolerance::default(),
    )
    .expect("ray should hit triangle");

    assert!((hit.distance - 1.0).abs() < 1.0e-12);
    assert!((hit.barycentric_u - 0.25).abs() < 1.0e-12);
    assert!((hit.barycentric_v - 0.25).abs() < 1.0e-12);
}

#[test]
fn point_in_closed_surface_classifies_cube_points() {
    let triangles = cube_triangles();
    let tolerance = MeshingTolerance::from_bounds([0.0, 0.0, 0.0], [1.0, 1.0, 1.0]);

    assert_eq!(
        point_in_closed_triangle_surface([0.5, 0.5, 0.5], &triangles, tolerance),
        PointInClosedSurface::Inside
    );
    assert_eq!(
        point_in_closed_triangle_surface([1.5, 0.5, 0.5], &triangles, tolerance),
        PointInClosedSurface::Outside
    );
    assert_eq!(
        point_in_closed_triangle_surface([1.0, 0.5, 0.5], &triangles, tolerance),
        PointInClosedSurface::OnBoundary
    );
}

#[test]
fn circumsphere_contains_regular_tetrahedron_center() {
    let tetrahedron = [
        [0.0, 0.0, 0.0],
        [1.0, 0.0, 0.0],
        [0.0, 1.0, 0.0],
        [0.0, 0.0, 1.0],
    ];

    assert!(tetrahedron_circumsphere_contains_point(
        tetrahedron,
        [0.5, 0.5, 0.5],
        MeshingTolerance::default()
    ));
    assert!(!tetrahedron_circumsphere_contains_point(
        tetrahedron,
        [2.0, 2.0, 2.0],
        MeshingTolerance::default()
    ));
}

fn cube_triangles() -> Vec<Triangle3> {
    let p = [
        [0.0, 0.0, 0.0],
        [1.0, 0.0, 0.0],
        [1.0, 1.0, 0.0],
        [0.0, 1.0, 0.0],
        [0.0, 0.0, 1.0],
        [1.0, 0.0, 1.0],
        [1.0, 1.0, 1.0],
        [0.0, 1.0, 1.0],
    ];
    [
        [p[0], p[2], p[1]],
        [p[0], p[3], p[2]],
        [p[4], p[5], p[6]],
        [p[4], p[6], p[7]],
        [p[0], p[1], p[5]],
        [p[0], p[5], p[4]],
        [p[1], p[2], p[6]],
        [p[1], p[6], p[5]],
        [p[2], p[3], p[7]],
        [p[2], p[7], p[6]],
        [p[3], p[0], p[4]],
        [p[3], p[4], p[7]],
    ]
    .to_vec()
}