conspire 0.6.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
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
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196

#[cfg(test)]
mod test;

use crate::{
    fem::block::element::{
        ElementNodalEitherCoordinates, FRAC_SQRT_3_5, FiniteElement, ParametricCoordinate,
        ParametricCoordinates, ParametricReference, ShapeFunctions, ShapeFunctionsGradients,
        quadratic::{M, QuadraticElement, QuadraticFiniteElement},
    },
    math::{Scalar, ScalarList},
};

const G: usize = 18;
const N: usize = 15;
const P: usize = N;

pub type Wedge = QuadraticElement<G, N>;

impl FiniteElement<G, M, N, P> for Wedge {
    fn integration_points() -> ParametricCoordinates<G, M> {
        const ONE_SIXTH: Scalar = 1.0 / 12.0;
        const TWO_THIRDS: Scalar = 2.0 / 3.0;
        [
            [ONE_SIXTH, ONE_SIXTH, FRAC_SQRT_3_5],
            [ONE_SIXTH, TWO_THIRDS, FRAC_SQRT_3_5],
            [TWO_THIRDS, ONE_SIXTH, FRAC_SQRT_3_5],
            [ONE_SIXTH, ONE_SIXTH, -FRAC_SQRT_3_5],
            [ONE_SIXTH, TWO_THIRDS, -FRAC_SQRT_3_5],
            [TWO_THIRDS, ONE_SIXTH, -FRAC_SQRT_3_5],
            [ONE_SIXTH, ONE_SIXTH, 0.0],
            [ONE_SIXTH, TWO_THIRDS, 0.0],
            [TWO_THIRDS, ONE_SIXTH, 0.0],
            [0.5, 0.5, FRAC_SQRT_3_5],
            [0.5, 0.0, FRAC_SQRT_3_5],
            [0.0, 0.5, FRAC_SQRT_3_5],
            [0.5, 0.5, -FRAC_SQRT_3_5],
            [0.5, 0.0, -FRAC_SQRT_3_5],
            [0.0, 0.5, -FRAC_SQRT_3_5],
            [0.5, 0.5, 0.0],
            [0.5, 0.0, 0.0],
            [0.0, 0.5, 0.0],
        ]
        .into()
    }
    fn integration_weights(&self) -> &ScalarList<G> {
        &self.integration_weights
    }
    fn parametric_reference() -> ParametricReference<M, N> {
        [
            [0.0, 0.0, -1.0],
            [1.0, 0.0, -1.0],
            [0.0, 1.0, -1.0],
            [0.0, 0.0, 1.0],
            [1.0, 0.0, 1.0],
            [0.0, 1.0, 1.0],
            [0.5, 0.0, -1.0],
            [0.5, 0.5, -1.0],
            [0.0, 0.5, -1.0],
            [0.0, 0.0, 0.0],
            [1.0, 0.0, 0.0],
            [0.0, 1.0, 0.0],
            [0.5, 0.0, 1.0],
            [0.5, 0.5, 1.0],
            [0.0, 0.5, 1.0],
        ]
        .into()
    }
    fn parametric_weights() -> ScalarList<G> {
        const ONE_TWELFTH: Scalar = 1.0 / 12.0;
        const TWO_FIFTEENTHS: Scalar = 2.0 / 15.0;
        const ONE_ONE_HUNDRED_EIGHTH: Scalar = 1.0 / 108.0;
        const TWO_ONE_HUNDRED_THIRTY_FIFTHS: Scalar = 2.0 / 135.0;
        [
            ONE_TWELFTH,
            ONE_TWELFTH,
            ONE_TWELFTH,
            ONE_TWELFTH,
            ONE_TWELFTH,
            ONE_TWELFTH,
            TWO_FIFTEENTHS,
            TWO_FIFTEENTHS,
            TWO_FIFTEENTHS,
            ONE_ONE_HUNDRED_EIGHTH,
            ONE_ONE_HUNDRED_EIGHTH,
            ONE_ONE_HUNDRED_EIGHTH,
            ONE_ONE_HUNDRED_EIGHTH,
            ONE_ONE_HUNDRED_EIGHTH,
            ONE_ONE_HUNDRED_EIGHTH,
            TWO_ONE_HUNDRED_THIRTY_FIFTHS,
            TWO_ONE_HUNDRED_THIRTY_FIFTHS,
            TWO_ONE_HUNDRED_THIRTY_FIFTHS,
        ]
        .into()
    }
    fn scaled_jacobians<const I: usize>(
        _nodal_coordinates: ElementNodalEitherCoordinates<I, N>,
    ) -> ScalarList<P> {
        todo!()
    }
    fn shape_functions(parametric_coordinate: ParametricCoordinate<M>) -> ShapeFunctions<N> {
        let [xi_1, xi_2, xi_3] = parametric_coordinate.into();
        let xi_0 = 1.0 - xi_1 - xi_2;
        [
            -0.5 * xi_0 * (1.0 - xi_3) * (2.0 * xi_1 + 2.0 * xi_2 + xi_3),
            0.5 * xi_1 * (1.0 - xi_3) * (2.0 * xi_1 - 2.0 - xi_3),
            0.5 * xi_2 * (1.0 - xi_3) * (2.0 * xi_2 - 2.0 - xi_3),
            -0.5 * xi_0 * (1.0 + xi_3) * (2.0 * xi_1 + 2.0 * xi_2 - xi_3),
            0.5 * xi_1 * (1.0 + xi_3) * (2.0 * xi_1 - 2.0 + xi_3),
            0.5 * xi_2 * (1.0 + xi_3) * (2.0 * xi_2 - 2.0 + xi_3),
            2.0 * xi_0 * xi_1 * (1.0 - xi_3),
            2.0 * xi_1 * xi_2 * (1.0 - xi_3),
            2.0 * xi_0 * xi_2 * (1.0 - xi_3),
            xi_0 * (1.0 - xi_3 * xi_3),
            xi_1 * (1.0 - xi_3 * xi_3),
            xi_2 * (1.0 - xi_3 * xi_3),
            2.0 * xi_0 * xi_1 * (1.0 + xi_3),
            2.0 * xi_1 * xi_2 * (1.0 + xi_3),
            2.0 * xi_0 * xi_2 * (1.0 + xi_3),
        ]
        .into()
    }
    fn shape_functions_gradients(
        parametric_coordinate: ParametricCoordinate<M>,
    ) -> ShapeFunctionsGradients<M, N> {
        let [xi_1, xi_2, xi_3] = parametric_coordinate.into();
        let xi_0 = 1.0 - xi_1 - xi_2;
        [
            [
                0.5 * (1.0 - xi_3) * (4.0 * xi_1 + 4.0 * xi_2 + xi_3 - 2.0),
                0.5 * (1.0 - xi_3) * (4.0 * xi_1 + 4.0 * xi_2 + xi_3 - 2.0),
                xi_0 * (xi_1 + xi_2 + xi_3 - 0.5),
            ],
            [
                0.5 * (1.0 - xi_3) * (4.0 * xi_1 - 2.0 - xi_3),
                0.0,
                xi_1 * (xi_3 - xi_1 + 0.5),
            ],
            [
                0.0,
                0.5 * (1.0 - xi_3) * (4.0 * xi_2 - xi_3 - 2.0),
                xi_2 * (xi_3 - xi_2 + 0.5),
            ],
            [
                0.5 * (1.0 + xi_3) * (4.0 * xi_1 + 4.0 * xi_2 - xi_3 - 2.0),
                0.5 * (1.0 + xi_3) * (4.0 * xi_1 + 4.0 * xi_2 - xi_3 - 2.0),
                xi_0 * (-xi_1 - xi_2 + xi_3 + 0.5),
            ],
            [
                0.5 * (1.0 + xi_3) * (4.0 * xi_1 - 2.0 + xi_3),
                0.0,
                xi_1 * (xi_3 + xi_1 - 0.5),
            ],
            [
                0.0,
                0.5 * (1.0 + xi_3) * (4.0 * xi_2 + xi_3 - 2.0),
                xi_2 * (xi_3 + xi_2 - 0.5),
            ],
            [
                2.0 * (1.0 - xi_3) * (1.0 - 2.0 * xi_1 - xi_2),
                -2.0 * xi_1 * (1.0 - xi_3),
                -2.0 * xi_0 * xi_1,
            ],
            [
                2.0 * xi_2 * (1.0 - xi_3),
                2.0 * xi_1 * (1.0 - xi_3),
                -2.0 * xi_1 * xi_2,
            ],
            [
                -2.0 * xi_2 * (1.0 - xi_3),
                2.0 * (1.0 - xi_3) * (xi_0 - xi_2),
                -2.0 * xi_0 * xi_2,
            ],
            [xi_3 * xi_3 - 1.0, xi_3 * xi_3 - 1.0, -2.0 * xi_0 * xi_3],
            [1.0 - xi_3 * xi_3, 0.0, -2.0 * xi_1 * xi_3],
            [0.0, 1.0 - xi_3 * xi_3, -2.0 * xi_2 * xi_3],
            [
                2.0 * (1.0 + xi_3) * (1.0 - 2.0 * xi_1 - xi_2),
                -2.0 * xi_1 * (1.0 + xi_3),
                2.0 * xi_0 * xi_1,
            ],
            [
                2.0 * xi_2 * (1.0 + xi_3),
                2.0 * xi_1 * (1.0 + xi_3),
                2.0 * xi_1 * xi_2,
            ],
            [
                -2.0 * xi_2 * (1.0 + xi_3),
                2.0 * (1.0 + xi_3) * (1.0 - xi_1 - 2.0 * xi_2),
                2.0 * xi_0 * xi_2,
            ],
        ]
        .into()
    }
}

impl QuadraticFiniteElement<G, N> for Wedge {}