Interface Mesh3d

splashsurf_lib/src/mesh.rs:269–372

Basic functionality that is provided by all meshes of the library Meshes consist of vertices and cells. Cells identify their associated vertices using indices into the mesh's slice of vertices.

/// Meshes consist of vertices and cells. Cells identify their associated vertices using indices
/// into the mesh's slice of vertices.
pub trait Mesh3d<R: Real>
where
    Self: Sized,
{
    /// The cell connectivity type of the mesh
    type Cell: CellConnectivity + Clone;

    /// Returns a slice of all vertices of the mesh
    fn vertices(&self) -> &[Vector3<R>];
    /// Returns a mutable slice of all vertices of the mesh
    fn vertices_mut(&mut self) -> &mut [Vector3<R>];
    /// Returns a slice of all cells of the mesh
    fn cells(&self) -> &[Self::Cell];

    /// Constructs a mesh from the given vertices and connectivity (does not check inputs for validity)
    fn from_vertices_and_connectivity(
        vertices: Vec<Vector3<R>>,
        connectivity: Vec<Self::Cell>,
    ) -> Self;

    /// Returns a mapping of all mesh vertices to the set of their connected neighbor vertices
    fn vertex_vertex_connectivity(&self) -> Vec<Vec<usize>> {
        profile!("vertex_vertex_connectivity");

        let mut connectivity_map: Vec<Vec<usize>> =
            vec![Vec::with_capacity(4); self.vertices().len()];
        for cell in self.cells() {
            for &i in cell.vertices() {
                for &j in cell.vertices() {
                    if i != j && !connectivity_map[i].contains(&j) {
                        connectivity_map[i].push(j);
                    }
                }
            }
        }

        connectivity_map
    }

    /// Returns a mapping of all mesh vertices to the set of the cells they belong to
    fn vertex_cell_connectivity(&self) -> Vec<Vec<usize>> {
        profile!("vertex_cell_connectivity");
        let mut connectivity_map: Vec<Vec<usize>> = vec![Vec::new(); self.vertices().len()];
        for (cell_idx, cell) in self.cells().iter().enumerate() {
            for &v_i in cell.vertices() {
                if !connectivity_map[v_i].contains(&cell_idx) {
                    connectivity_map[v_i].push(cell_idx);
                }
            }
        }

        connectivity_map
    }

    /// Returns a new mesh containing only the specified cells and removes all unreferenced vertices
    fn keep_cells(&self, cell_indices: &[usize], keep_vertices: bool) -> Self {
        if keep_vertices {
            keep_cells_impl(self, cell_indices, &[])