runmat-geometry-core 0.5.6

Canonical geometry domain model 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

use serde::{Deserialize, Serialize};

#[derive(Debug, Clone, Copy, PartialEq, Eq, Serialize, Deserialize)]
#[serde(rename_all = "snake_case")]
pub enum MeshKind {
    Surface,
    Volume,
}

#[derive(Debug, Clone, PartialEq, Eq, Serialize, Deserialize)]
pub struct MeshDescriptor {
    pub mesh_id: String,
    pub kind: MeshKind,
    pub vertex_count: u64,
    pub element_count: u64,
}

#[derive(Debug, Clone, PartialEq, Serialize, Deserialize)]
pub struct SurfaceMesh {
    pub mesh_id: String,
    pub vertices: Vec<[f64; 3]>,
    pub triangles: Vec<[u32; 3]>,
}

impl SurfaceMesh {
    pub fn new(
        mesh_id: impl Into<String>,
        vertices: Vec<[f64; 3]>,
        triangles: Vec<[u32; 3]>,
    ) -> Self {
        Self {
            mesh_id: mesh_id.into(),
            vertices,
            triangles,
        }
    }

    pub fn validate(&self) -> Result<(), &'static str> {
        if self.mesh_id.is_empty() {
            return Err("surface mesh id must not be empty");
        }
        if self
            .vertices
            .iter()
            .flatten()
            .any(|coordinate| !coordinate.is_finite())
        {
            return Err("surface mesh vertices must be finite");
        }
        let vertex_count = self.vertices.len();
        if self
            .triangles
            .iter()
            .flatten()
            .any(|index| *index as usize >= vertex_count)
        {
            return Err("surface mesh triangle index out of bounds");
        }
        Ok(())
    }
}