conspire 0.7.1

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

use crate::{
    geometry::{
        Coordinate,
        mesh::{Connectivity, Mesh, Verdict},
    },
    math::Tensor,
};

fn grid() -> (Vec<[usize; 8]>, Vec<[f64; 3]>) {
    let id = |i: usize, j: usize, k: usize| i + 3 * j + 9 * k;
    let mut coordinates = Vec::new();
    for k in 0..3 {
        for j in 0..3 {
            for i in 0..3 {
                coordinates.push([i as f64, j as f64, k as f64]);
            }
        }
    }
    let mut connectivity = Vec::new();
    for k in 0..2 {
        for j in 0..2 {
            for i in 0..2 {
                connectivity.push([
                    id(i, j, k),
                    id(i + 1, j, k),
                    id(i + 1, j + 1, k),
                    id(i, j + 1, k),
                    id(i, j, k + 1),
                    id(i + 1, j, k + 1),
                    id(i + 1, j + 1, k + 1),
                    id(i, j + 1, k + 1),
                ]);
            }
        }
    }
    (connectivity, coordinates)
}

fn mesh(connectivity: Vec<[usize; 8]>, coordinates: Vec<[f64; 3]>) -> Mesh<3> {
    Mesh::from((
        vec![Connectivity::Hexahedral(connectivity.into())],
        coordinates.into(),
    ))
}

fn minimum_scaled(mesh: &Mesh<3>) -> f64 {
    mesh.minimum_scaled_jacobians()
        .into_iter()
        .flatten()
        .fold(f64::INFINITY, f64::min)
}

#[test]
fn interior_node_recenters_without_dropping_quality() {
    let (connectivity, mut coordinates) = grid();
    coordinates[13] = [1.4, 1.3, 0.6];
    let mut mesh = mesh(connectivity, coordinates);
    let center = Coordinate::from([1.0, 1.0, 1.0]);
    let before = minimum_scaled(&mesh);
    let offset = (&mesh.coordinates()[13] - &center).norm();
    mesh.smart_laplace_smooth(10, 1.0);
    assert!(minimum_scaled(&mesh) >= before);
    assert!((&mesh.coordinates()[13] - &center).norm() < offset);
}

#[test]
fn perfect_grid_does_not_move() {
    let (connectivity, coordinates) = grid();
    let mut mesh = mesh(connectivity, coordinates);
    let before: Vec<_> = mesh.coordinates().iter().cloned().collect();
    mesh.smart_laplace_smooth(5, 1.0);
    mesh.coordinates()
        .iter()
        .zip(before)
        .for_each(|(now, was)| assert!((now - &was).norm() < 1.0e-12));
}