del-geo-core 0.1.38

2D/3D geometry utility codes
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

pub const DEGREE2IDX: [usize; 4] = [0, 1, 4, 11];

#[allow(clippy::excessive_precision)]
pub const IDX2POSWEIGHT_F32: [[f32; 3]; 11] = [
    // linear
    [0.3333333333, 0.3333333333, 1.0],
    // quadratic
    [0.1666666667, 0.1666666667, 0.3333333333],
    [0.6666666667, 0.1666666667, 0.3333333333],
    [0.1666666667, 0.6666666667, 0.3333333333],
    // cubic
    [0.1012865073, 0.1012865073, 0.1259391805],
    [0.7974269854, 0.1012865073, 0.1259391805],
    [0.1012865073, 0.7974269854, 0.1259391805],
    [0.4701420641, 0.0597158718, 0.1323941527],
    [0.4701420641, 0.4701420641, 0.1323941527],
    [0.0597158718, 0.4701420641, 0.1323941527],
    [0.3333333333, 0.3333333333, 0.225],
];

#[allow(clippy::excessive_precision)]
pub const IDX2POSWEIGHT_F64: [[f64; 3]; 11] = [
    // linear
    [0.3333333333, 0.3333333333, 1.0],
    // quadratic
    [0.1666666667, 0.1666666667, 0.3333333333],
    [0.6666666667, 0.1666666667, 0.3333333333],
    [0.1666666667, 0.6666666667, 0.3333333333],
    // cubic
    [0.1012865073, 0.1012865073, 0.1259391805],
    [0.7974269854, 0.1012865073, 0.1259391805],
    [0.1012865073, 0.7974269854, 0.1259391805],
    [0.4701420641, 0.0597158718, 0.1323941527],
    [0.4701420641, 0.4701420641, 0.1323941527],
    [0.0597158718, 0.4701420641, 0.1323941527],
    [0.3333333333, 0.3333333333, 0.225],
];

pub struct Quad<Real> {
    _marker: std::marker::PhantomData<fn() -> Real>,
}

pub trait QuadratureTri<Real> {
    fn hoge(i_gauss_degree: usize) -> &'static [[Real; 3]];
}

impl<Real> QuadratureTri<f64> for Quad<Real> {
    fn hoge(i_gauss_degree: usize) -> &'static [[f64; 3]] {
        let i0 = DEGREE2IDX[i_gauss_degree];
        let i1 = DEGREE2IDX[i_gauss_degree + 1];
        &IDX2POSWEIGHT_F64[i0..i1]
    }
}

impl<Real> QuadratureTri<f32> for Quad<Real> {
    fn hoge(i_gauss_degree: usize) -> &'static [[f32; 3]] {
        let i0 = DEGREE2IDX[i_gauss_degree];
        let i1 = DEGREE2IDX[i_gauss_degree + 1];
        &IDX2POSWEIGHT_F32[i0..i1]
    }
}