conspire 0.7.2

The Rust interface to conspire.
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

use super::super::geometry::star_volume;
use super::super::test::{hexahedron, signed_volumes, sphere};
use crate::{
    geometry::{
        Coordinates,
        mesh::{Connectivity, Mesh},
    },
    math::CrossProduct,
};
use std::collections::HashSet;

#[test]
fn assemble_single_hexahedron() {
    let tessellation = sphere(3);
    let mesh = hexahedron([0.9, -0.1, -0.1], [1.1, 0.1, 0.1]);
    let classes = tessellation.classify(&mesh);
    let tables = tessellation
        .tables(&mesh, &classes, &HashSet::new())
        .unwrap();
    let result = tessellation.assemble(&mesh, &classes, &tables).unwrap();
    assert_eq!(result.number_of_element_blocks(), 1);
    assert_eq!(result.number_of_nodes(), 8);
    match &result.connectivities()[0] {
        Connectivity::Hexahedral(hexes) => {
            assert_eq!(hexes.iter().count(), 1);
            let element: Vec<usize> = hexes.iter().flatten().copied().collect();
            let coordinates = result.coordinates();
            let base = &(&coordinates[element[1]] - &coordinates[element[0]])
                .cross(&(&coordinates[element[3]] - &coordinates[element[0]]))
                * &(&coordinates[element[4]] - &coordinates[element[0]]);
            assert!(base > 0.0)
        }
        _ => panic!(),
    }
}

#[test]
fn agglomerate_sliver() {
    let coordinates = Coordinates::from(vec![
        [-1.0, -1.0, -2.0],
        [3.0, -1.0, -2.0],
        [3.0, 2.0, -2.0],
        [-1.0, 2.0, -2.0],
        [-1.0, -1.0, -0.04],
        [3.0, -1.0, 0.28],
        [3.0, 2.0, 0.28],
        [-1.0, 2.0, -0.04],
    ]);
    let triangles = vec![
        [0, 3, 2],
        [0, 2, 1],
        [4, 5, 6],
        [4, 6, 7],
        [0, 1, 5],
        [0, 5, 4],
        [1, 2, 6],
        [1, 6, 5],
        [3, 7, 6],
        [3, 6, 2],
        [3, 0, 4],
        [3, 4, 7],
    ];
    let tessellation = crate::geometry::mesh::tessellation::Tessellation::from(Mesh::from((
        vec![Connectivity::Triangular(triangles.into())],
        coordinates,
    )));
    let mesh = Mesh::from((
        vec![Connectivity::Hexahedral(
            vec![[0, 1, 2, 3, 4, 5, 6, 7], [1, 8, 9, 2, 5, 10, 11, 6]].into(),
        )],
        Coordinates::from(vec![
            [0.0, 0.0, 0.0],
            [1.0, 0.0, 0.0],
            [1.0, 1.0, 0.0],
            [0.0, 1.0, 0.0],
            [0.0, 0.0, 1.0],
            [1.0, 0.0, 1.0],
            [1.0, 1.0, 1.0],
            [0.0, 1.0, 1.0],
            [2.0, 0.0, 0.0],
            [2.0, 1.0, 0.0],
            [2.0, 0.0, 1.0],
            [2.0, 1.0, 1.0],
        ]),
    ));
    let classes = tessellation.classify(&mesh);
    assert_eq!(
        classes,
        vec![super::super::Class::Cut, super::super::Class::Cut]
    );
    let tables = tessellation
        .tables(&mesh, &classes, &HashSet::new())
        .unwrap();
    let result = tessellation.assemble(&mesh, &classes, &tables).unwrap();
    assert_eq!(result.number_of_element_blocks(), 1);
    match &result.connectivities()[0] {
        Connectivity::Polyhedral(polyhedra) => {
            assert_eq!(polyhedra.elements_faces().len(), 1);
            assert_eq!(polyhedra.elements_faces()[0].len(), 9);
            let faces: Vec<Vec<usize>> = polyhedra.elements_faces()[0]
                .iter()
                .map(|&face| polyhedra.faces_nodes()[face].clone())
                .collect();
            let volume = star_volume(&faces, result.coordinates());
            assert!((volume - 0.220095389507154).abs() < 1e-12, "{volume}");
            let signed = signed_volumes(polyhedra, result.coordinates())[0];
            assert!((signed - 0.220095389507154).abs() < 1e-12, "{signed}")
        }
        _ => panic!(),
    }
}