honeycomb-kernels 0.11.0

Meshing kernel implementations using combinatorial maps
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

//! clipping operation routines

use std::collections::VecDeque;

use honeycomb_core::cmap::{CMap2, DartIdType, FaceIdType, NULL_DART_ID, OrbitPolicy};
use honeycomb_core::geometry::{CoordsFloat, Vertex2};
use rustc_hash::FxHashSet as HashSet;

use crate::grisubal::GrisubalError;
use crate::grisubal::model::Boundary;
use crate::utils::VertexAnchor;

/// Clip content on the left side of the boundary.
pub fn clip_left<T: CoordsFloat>(cmap: &mut CMap2<T>) -> Result<(), GrisubalError> {
    // color faces using a bfs starting on multiple nodes
    let marked = mark_faces(cmap, Boundary::Left, Boundary::Right)?;

    // save vertices & split boundary
    delete_darts(cmap, marked, Boundary::Right);

    Ok(())
}

/// Clip content on the right side of the boundary.
pub fn clip_right<T: CoordsFloat>(cmap: &mut CMap2<T>) -> Result<(), GrisubalError> {
    // color faces using a bfs starting on multiple nodes
    let marked = mark_faces(cmap, Boundary::Right, Boundary::Left)?;

    // save vertices & split boundary
    delete_darts(cmap, marked, Boundary::Left);

    Ok(())
}

// --- internals

#[allow(clippy::cast_possible_truncation)]
fn mark_faces<T: CoordsFloat>(
    cmap: &CMap2<T>,
    mark: Boundary,
    other: Boundary,
) -> Result<HashSet<FaceIdType>, GrisubalError> {
    let mut marked: HashSet<FaceIdType> = HashSet::default();
    marked.insert(0);
    let mut queue: VecDeque<FaceIdType> = (1..cmap.n_darts() as DartIdType)
        .filter_map(|dart_id| {
            // use darts on the left side of the boundary as starting points to walk through faces
            if cmap.read_attribute::<Boundary>(dart_id) == Some(mark) && !cmap.is_free(dart_id) {
                return Some(cmap.face_id(dart_id));
            }
            None
        })
        .collect();
    while let Some(face_id) = queue.pop_front() {
        // mark faces
        if marked.insert(face_id) {
            // detect orientation issues / open geometries
            let mut darts = cmap.orbit(OrbitPolicy::Face, face_id as DartIdType);
            if darts.any(|did| cmap.read_attribute::<Boundary>(did) == Some(other)) {
                return Err(GrisubalError::InconsistentOrientation(
                    "between-boundary inconsistency",
                ));
            }
            // find neighbor faces where entry darts aren't tagged
            let darts = cmap.orbit(OrbitPolicy::Face, face_id as DartIdType);
            queue.extend(darts.filter_map(|dart_id| {
                if matches!(
                    cmap.read_attribute::<Boundary>(cmap.beta::<2>(dart_id)),
                    Some(Boundary::None) | None
                ) {
                    return Some(cmap.face_id(cmap.beta::<2>(dart_id)));
                }
                None
            }));
        }
    }
    marked.remove(&0);
    Ok(marked)
}

#[allow(clippy::cast_possible_truncation)]
fn delete_darts<T: CoordsFloat>(
    cmap: &mut CMap2<T>,
    marked: HashSet<FaceIdType>,
    kept_boundary: Boundary,
) {
    let kept_boundary_components: Vec<(DartIdType, Vertex2<T>, Option<VertexAnchor>)> =
        (1..cmap.n_darts() as DartIdType)
            .filter_map(|dart_id| {
                if cmap.read_attribute::<Boundary>(dart_id) == Some(kept_boundary) {
                    return Some((
                        dart_id,
                        cmap.read_vertex(cmap.vertex_id(dart_id))
                            .expect("E: found a topological vertex with no associated coordinates"),
                        if cmap.contains_attribute::<VertexAnchor>() {
                            cmap.read_attribute(cmap.vertex_id(dart_id)) // may be Some or None
                        } else {
                            None
                        },
                    ));
                }
                None
            })
            .collect();

    for face_id in marked {
        let darts: Vec<DartIdType> = cmap
            .orbit(OrbitPolicy::Face, face_id as DartIdType)
            .collect();
        for &dart in &darts {
            let _ = cmap.remove_vertex(cmap.vertex_id(dart));
            cmap.set_betas(dart, [NULL_DART_ID; 3]);
            cmap.release_dart(dart).expect("E: unreachable");
        }
    }

    for (dart, vertex, anchor) in kept_boundary_components {
        cmap.set_beta::<2>(dart, NULL_DART_ID); // set beta2(dart) to 0
        let vid = cmap.vertex_id(dart);
        cmap.write_vertex(vid, vertex);
        // if the map does not contain anchors, this branch is never taken
        if let Some(a) = anchor {
            cmap.write_attribute(vid, a);
        }
    }
}