honeycomb_kernels/utils/
routines.rs

1use honeycomb_core::{
2    cmap::{CMap2, DartIdType, OrbitPolicy, VertexIdType},
3    geometry::{CoordsFloat, Vertex2},
4    stm::{StmClosureResult, Transaction, retry},
5};
6use smallvec::SmallVec;
7
8/// Check if all faces incident to the vertex have the same orientation.
9///
10/// Note that this function expects the incident faces to be triangles.
11///
12/// # Errors
13///
14/// This method is meant to be called in a context where the returned `Result` is used to
15/// validate the transaction passed as argument. Errors should not be processed manually,
16/// only processed via the `?` operator.
17pub fn is_orbit_orientation_consistent<T: CoordsFloat>(
18    t: &mut Transaction,
19    map: &CMap2<T>,
20    vid: VertexIdType,
21) -> StmClosureResult<bool> {
22    if let Some(new_v) = map.read_vertex(t, vid)? {
23        let mut tmp: SmallVec<DartIdType, 10> = SmallVec::new();
24        for d in map.orbit_transac(t, OrbitPolicy::Vertex, vid) {
25            tmp.push(d?);
26        }
27
28        let ref_sign = {
29            let d = tmp[0];
30            let b1d = map.beta_transac::<1>(t, d)?;
31            let b1b1d = map.beta_transac::<1>(t, b1d)?;
32            let vid1 = map.vertex_id_transac(t, b1d)?;
33            let vid2 = map.vertex_id_transac(t, b1b1d)?;
34            let v1 = if let Some(v) = map.read_vertex(t, vid1)? {
35                v
36            } else {
37                retry()?
38            };
39            let v2 = if let Some(v) = map.read_vertex(t, vid2)? {
40                v
41            } else {
42                retry()?
43            };
44
45            let crossp = Vertex2::cross_product_from_vertices(&new_v, &v1, &v2);
46            crossp.signum()
47        };
48        for &d in &tmp[1..] {
49            let b1d = map.beta_transac::<1>(t, d)?;
50            let b1b1d = map.beta_transac::<1>(t, b1d)?;
51            let vid1 = map.vertex_id_transac(t, b1d)?;
52            let vid2 = map.vertex_id_transac(t, b1b1d)?;
53            let v1 = if let Some(v) = map.read_vertex(t, vid1)? {
54                v
55            } else {
56                retry()?
57            };
58            let v2 = if let Some(v) = map.read_vertex(t, vid2)? {
59                v
60            } else {
61                retry()?
62            };
63
64            let crossp = Vertex2::cross_product_from_vertices(&new_v, &v1, &v2);
65
66            if ref_sign != crossp.signum() {
67                return Ok(false);
68            }
69        }
70    } else {
71        retry()?;
72    }
73
74    Ok(true)
75}
honeycomb_kernels/utils/routines.rs

honeycomb_kernels/utils/
routines.rs
