conspire 0.7.0

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

#[cfg(test)]
mod test;

use crate::{
    geometry::Coordinates,
    math::{Scalar, Tensor},
};

const CORNERS: [[usize; 2]; 4] = [[1, 3], [2, 0], [3, 1], [0, 2]];

const EDGES: [[usize; 2]; 4] = [[0, 1], [1, 2], [2, 3], [3, 0]];

pub(super) fn maximum_edge_ratio<const D: usize>(
    element: &[usize],
    coordinates: &Coordinates<D>,
) -> Scalar {
    super::maximum_edge_ratio(&EDGES, element, coordinates)
}

pub(super) fn minimum_jacobian<const D: usize>(
    element: &[usize],
    coordinates: &Coordinates<D>,
) -> Scalar {
    super::min_jacobian(&CORNERS, element, coordinates)
}

pub(super) fn minimum_scaled_jacobian<const D: usize>(
    element: &[usize],
    coordinates: &Coordinates<D>,
) -> Scalar {
    super::min_scaled_jacobian(&CORNERS, element, coordinates, 1.0)
}

pub(super) fn maximum_skew<const D: usize>(
    element: &[usize],
    coordinates: &Coordinates<D>,
) -> Scalar {
    let p = |i: usize| &coordinates[element[i]];
    let x1 = (p(1) - p(0)) + (p(2) - p(3));
    let x2 = (p(3) - p(0)) + (p(2) - p(1));
    let (n1, n2) = (x1.norm(), x2.norm());
    if n1 > 0.0 && n2 > 0.0 {
        ((&x1 * &x2) / (n1 * n2)).abs()
    } else {
        0.0
    }
}

pub(super) fn volume<const D: usize>(element: &[usize], coordinates: &Coordinates<D>) -> Scalar {
    super::triangle_area(&[element[0], element[1], element[2]], coordinates)
        + super::triangle_area(&[element[0], element[2], element[3]], coordinates)
}