mesh-sieve 4.0.1

Modular, high-performance Rust library for mesh and data management, designed for scientific computing and PDE codes.
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

// cargo run --example meshgen_basic
use mesh_sieve::algs::meshgen::{
    MeshGenOptions, StructuredCellType, cylinder_shell, sphere_shell, structured_box_2d,
    structured_box_3d,
};
use mesh_sieve::data::section::Section;
use mesh_sieve::data::storage::VecStorage;
use mesh_sieve::mesh_error::MeshSieveError;
use mesh_sieve::topology::cell_type::CellType;

fn count_type(section: &Section<CellType, VecStorage<CellType>>, cell_type: CellType) -> usize {
    section
        .iter()
        .filter(|(_, slice)| slice[0] == cell_type)
        .count()
}

fn main() -> Result<(), MeshSieveError> {
    let tri_mesh = structured_box_2d(
        1,
        1,
        [0.0, 0.0],
        [1.0, 1.0],
        StructuredCellType::Triangle,
        MeshGenOptions::default(),
    )?;
    let tri_types = tri_mesh.cell_types.as_ref().unwrap();
    assert_eq!(count_type(tri_types, CellType::Vertex), 4);
    assert_eq!(count_type(tri_types, CellType::Triangle), 2);

    let hex_mesh = structured_box_3d(
        1,
        1,
        1,
        [0.0, 0.0, 0.0],
        [1.0, 1.0, 1.0],
        StructuredCellType::Hexahedron,
        MeshGenOptions::default(),
    )?;
    let hex_types = hex_mesh.cell_types.as_ref().unwrap();
    assert_eq!(count_type(hex_types, CellType::Vertex), 8);
    assert_eq!(count_type(hex_types, CellType::Hexahedron), 1);

    let sphere = sphere_shell(1.0, 2, 4, MeshGenOptions::default())?;
    let sphere_types = sphere.cell_types.as_ref().unwrap();
    assert_eq!(count_type(sphere_types, CellType::Vertex), 6);
    assert_eq!(count_type(sphere_types, CellType::Triangle), 8);

    let cylinder = cylinder_shell(1.0, 1.0, 4, 1, MeshGenOptions::default())?;
    let cylinder_types = cylinder.cell_types.as_ref().unwrap();
    assert_eq!(count_type(cylinder_types, CellType::Vertex), 8);
    assert_eq!(count_type(cylinder_types, CellType::Quadrilateral), 4);

    Ok(())
}