quadraturerules 0.9.0

quadrature rules library
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

//! Hammer--Marlowe--Stroud quadrature rule.

use crate::Domain;

/// Get a Hammer--Marlowe--Stroud quadrature rule.
pub fn hammer_marlowe_stroud(
    domain: Domain,
    order: usize,
) -> Result<(Vec<f64>, Vec<f64>), &'static str> {
    match domain {
        Domain::Interval => match order {
            _ => Err("Invalid quadrature order"),
        },
        Domain::Quadrilateral => match order {
            _ => Err("Invalid quadrature order"),
        },
        Domain::EdgeAdjacentQuadrilaterals => match order {
            _ => Err("Invalid quadrature order"),
        },
        Domain::VertexAdjacentQuadrilaterals => match order {
            _ => Err("Invalid quadrature order"),
        },
        Domain::Triangle => match order {
            4 => Ok((
                vec![
                    0.13005607921683443,
                    0.009703785126946106,
                    0.8602401356562195,
                    0.009703785126946106,
                    0.13005607921683443,
                    0.8602401356562195,
                    0.0936377844373285,
                    0.04612207990645205,
                    0.8602401356562195,
                    0.04612207990645205,
                    0.0936377844373285,
                    0.8602401356562195,
                    0.38749748340669415,
                    0.028912084224389012,
                    0.5835904323689168,
                    0.028912084224389012,
                    0.38749748340669415,
                    0.5835904323689168,
                    0.2789904634965088,
                    0.13741910413457437,
                    0.5835904323689168,
                    0.13741910413457437,
                    0.2789904634965088,
                    0.5835904323689168,
                    0.6729468631505064,
                    0.050210123211369806,
                    0.2768430136381237,
                    0.050210123211369806,
                    0.6729468631505064,
                    0.2768430136381237,
                    0.4845083266304333,
                    0.23864865973144297,
                    0.2768430136381237,
                    0.23864865973144297,
                    0.4845083266304333,
                    0.2768430136381237,
                    0.8774288093304679,
                    0.06546699455501448,
                    0.05710419611451767,
                    0.06546699455501448,
                    0.8774288093304679,
                    0.05710419611451767,
                    0.6317312516411253,
                    0.311164552244357,
                    0.05710419611451767,
                    0.311164552244357,
                    0.6317312516411253,
                    0.05710419611451767,
                ],
                vec![
                    0.010846451821050509,
                    0.010846451821050509,
                    0.020334519128957573,
                    0.020334519128957573,
                    0.04516809856473986,
                    0.04516809856473986,
                    0.08467944904349257,
                    0.08467944904349257,
                    0.07077613579617188,
                    0.07077613579617188,
                    0.13268843221409946,
                    0.13268843221409946,
                    0.047136736386764674,
                    0.047136736386764674,
                    0.08837017704472347,
                    0.08837017704472347,
                ],
            )),
            _ => Err("Invalid quadrature order"),
        },
        Domain::EdgeAdjacentTriangleAndQuadrilateral => match order {
            _ => Err("Invalid quadrature order"),
        },
        Domain::VertexAdjacentTriangleAndQuadrilateral => match order {
            _ => Err("Invalid quadrature order"),
        },
        Domain::EdgeAdjacentTriangles => match order {
            _ => Err("Invalid quadrature order"),
        },
        Domain::VertexAdjacentTriangles => match order {
            _ => Err("Invalid quadrature order"),
        },
        Domain::Hexahedron => match order {
            _ => Err("Invalid quadrature order"),
        },
        Domain::Tetrahedron => match order {
            _ => Err("Invalid quadrature order"),
        },
        Domain::SquareBasedPyramid => match order {
            _ => Err("Invalid quadrature order"),
        },
        Domain::TriangularPrism => match order {
            _ => Err("Invalid quadrature order"),
        },
    }
}