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
77
78
79
80
81
82
83

use crate::{
    geometry::{
        Coordinate, Coordinates,
        mesh::{Connectivity, Mesh, differential::laplace::Weighting},
    },
    math::test::{TestError, assert_eq_within_tols},
};

fn triangle(coordinates: [Coordinate<3>; 3]) -> Mesh<3> {
    let connectivities = vec![Connectivity::Triangular(vec![[0_usize, 1, 2]].into())];
    Mesh::from((connectivities, Coordinates::from(coordinates)))
}

#[test]
fn single_triangle() -> Result<(), TestError> {
    let laplacian = triangle([
        Coordinate::const_from([0.0, 0.0, 0.0]),
        Coordinate::const_from([2.0, 0.0, 0.0]),
        Coordinate::const_from([0.0, 2.0, 0.0]),
    ])
    .laplacian(Weighting::Uniform);
    assert_eq_within_tols(&laplacian[0], &[-1.0, -1.0, 0.0].into())?;
    assert_eq_within_tols(&laplacian[1], &[2.0, -1.0, 0.0].into())?;
    assert_eq_within_tols(&laplacian[2], &[-1.0, 2.0, 0.0].into())
}

#[test]
fn vertex_at_neighbor_centroid_is_fixed() -> Result<(), TestError> {
    let laplacian = triangle([
        Coordinate::const_from([0.0, 0.0, 0.0]),
        Coordinate::const_from([2.0, 0.0, 0.0]),
        Coordinate::const_from([1.0, 0.0, 0.0]),
    ])
    .laplacian(Weighting::Uniform);
    assert_eq_within_tols(&laplacian[2], &[0.0, 0.0, 0.0].into())
}

#[test]
fn translation_invariant() -> Result<(), TestError> {
    let base = triangle([
        Coordinate::const_from([0.0, 0.0, 0.0]),
        Coordinate::const_from([2.0, 0.0, 0.0]),
        Coordinate::const_from([0.0, 2.0, 0.0]),
    ])
    .laplacian(Weighting::Uniform);
    let shifted = triangle([
        Coordinate::const_from([5.0, -3.0, 7.0]),
        Coordinate::const_from([7.0, -3.0, 7.0]),
        Coordinate::const_from([5.0, -1.0, 7.0]),
    ])
    .laplacian(Weighting::Uniform);
    (0..3).try_for_each(|n| assert_eq_within_tols(&base[n], &shifted[n]))
}

#[test]
fn cotangent_single_right_triangle() -> Result<(), TestError> {
    let laplacian = triangle([
        Coordinate::const_from([0.0, 0.0, 0.0]),
        Coordinate::const_from([2.0, 0.0, 0.0]),
        Coordinate::const_from([0.0, 2.0, 0.0]),
    ])
    .laplacian(Weighting::Cotangent);
    assert_eq_within_tols(&laplacian[0], &[-1.0, -1.0, 0.0].into())?;
    assert_eq_within_tols(&laplacian[1], &[2.0, 0.0, 0.0].into())?;
    assert_eq_within_tols(&laplacian[2], &[0.0, 2.0, 0.0].into())
}

#[test]
fn cotangent_translation_invariant() -> Result<(), TestError> {
    let base = triangle([
        Coordinate::const_from([0.0, 0.0, 0.0]),
        Coordinate::const_from([2.0, 0.0, 0.0]),
        Coordinate::const_from([0.0, 2.0, 0.0]),
    ])
    .laplacian(Weighting::Cotangent);
    let shifted = triangle([
        Coordinate::const_from([5.0, -3.0, 7.0]),
        Coordinate::const_from([7.0, -3.0, 7.0]),
        Coordinate::const_from([5.0, -1.0, 7.0]),
    ])
    .laplacian(Weighting::Cotangent);
    (0..3).try_for_each(|n| assert_eq_within_tols(&base[n], &shifted[n]))
}