caexfer 0.2.0

Transfer finite-element meshes, voxel occupancy, and results across CAE formats
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

//! Shared association of Nastran GRID displacement rows with a companion mesh.
use crate::core::{Dataset, Error, Field, FieldLocation, Mesh, Result};
use std::collections::BTreeMap;

/// Build one normalized displacement field after a format-specific decoder.
pub(crate) fn displacement_dataset(
    mesh: &Mesh,
    rows: &BTreeMap<u64, [f64; 6]>,
    width: usize,
    subcase: i64,
    step: usize,
    time: Option<f64>,
    code: &'static str,
) -> Result<Dataset> {
    // Check the companion mesh and result dimensions before allocating values.
    mesh.validate()?;
    if !matches!(width, 3 | 6) || rows.len() != mesh.points.len() {
        return Err(Error::new(
            code,
            "displacement nodes do not match mesh nodes",
        ));
    }
    // Bound the flattened field size even if the point count is extreme.
    let capacity = mesh
        .points
        .len()
        .checked_mul(width)
        .ok_or_else(|| Error::new(code, "displacement value count overflows"))?;
    let mut values = Vec::with_capacity(capacity);
    // Use mesh order for field values; match rows by original GRID ID.
    for point in &mesh.points {
        let row = rows
            .get(&point.id)
            .ok_or_else(|| Error::new(code, format!("no displacement for GRID {}", point.id)))?;
        values.extend_from_slice(&row[..width]);
    }
    // Store the selected result step in the field's signed metadata type.
    let step = i64::try_from(step).map_err(|_| Error::new(code, "result step exceeds Int64"))?;
    // Omit rotation labels when the decoder supplied translations only.
    let components = ["T1", "T2", "T3", "R1", "R2", "R3"]
        .into_iter()
        .take(width)
        .map(str::to_owned)
        .collect();
    let dataset = Dataset {
        mesh: mesh.clone(),
        fields: vec![Field {
            name: format!("DISPLACEMENT_SUBCASE_{subcase}"),
            location: FieldLocation::Point,
            components,
            values,
            step: Some(step),
            time,
        }],
    };
    // Check that the field layout still matches the companion mesh.
    dataset.validate()?;
    Ok(dataset)
}