BELFEM 0.9.0
Berkeley Lab Finite Element Framework
Loading...
Searching...
No Matches
fn_intpoints_gauss_penta50.hpp
Go to the documentation of this file.
1/*
2 * BELFEM -- The Berkeley Lab Finite Element Framework
3 * Copyright (c) 2026, The Regents of the University of California,
4 * through Lawrence Berkeley National Laboratory (subject to receipt of any required
5 * approvals from the U.S. Dept. of Energy). All rights reserved.
6 *
7 * Developers: Christian Messe, Gregory Giard
8 *
9 * See the top-level LICENSE file for the complete license and disclaimer.
10 */
11
12#ifndef BELFEM_FN_INTPOINTS_GAUSS_PENTA50_HPP
13#define BELFEM_FN_INTPOINTS_GAUSS_PENTA50_HPP
14
15#include "typedefs.hpp"
16#include "cl_Vector.hpp"
17#include "cl_Matrix.hpp"
18
19namespace belfem
20{
21 namespace integration
22 {
23// ----------------------------------------------------------------------------
24
25 // source: 10.1016/j.compfluid.2013.01.002
26 void
27 gauss_penta50(
28 Vector< real > & aWeights,
29 Matrix< real > & aPoints )
30 {
31 aPoints.set_size( 3, 50 );
32
33 aPoints( 0, 0 ) = 0.33333333333333;
34 aPoints( 1, 0 ) = 0.33333333333333;
35 aPoints( 2, 0 ) = 0.00000000000000;
36
37 aPoints( 0, 1 ) = 0.45963760180813;
38 aPoints( 1, 1 ) = 0.08072479638374;
39 aPoints( 2, 1 ) = 0.00000000000000;
40
41 aPoints( 0, 2 ) = 0.45963760180813;
42 aPoints( 1, 2 ) = 0.45963760180813;
43 aPoints( 2, 2 ) = 0.00000000000000;
44
45 aPoints( 0, 3 ) = 0.08072479638374;
46 aPoints( 1, 3 ) = 0.45963760180813;
47 aPoints( 2, 3 ) = 0.00000000000000;
48
49 aPoints( 0, 4 ) = 0.33333333333333;
50 aPoints( 1, 4 ) = 0.33333333333333;
51 aPoints( 2, 4 ) = 0.65690685402812;
52
53 aPoints( 0, 5 ) = 0.33333333333333;
54 aPoints( 1, 5 ) = 0.33333333333333;
55 aPoints( 2, 5 ) = -0.65690685402812;
56
57 aPoints( 0, 6 ) = 0.33333333333333;
58 aPoints( 1, 6 ) = 0.33333333333333;
59 aPoints( 2, 6 ) = 0.22289042139052;
60
61 aPoints( 0, 7 ) = 0.33333333333333;
62 aPoints( 1, 7 ) = 0.33333333333333;
63 aPoints( 2, 7 ) = -0.22289042139052;
64
65 aPoints( 0, 8 ) = 0.45891140214883;
66 aPoints( 1, 8 ) = 0.08217719570233;
67 aPoints( 2, 8 ) = 0.87508982900993;
68
69 aPoints( 0, 9 ) = 0.45891140214883;
70 aPoints( 1, 9 ) = 0.08217719570233;
71 aPoints( 2, 9 ) = -0.87508982900993;
72
73 aPoints( 0, 10 ) = 0.45891140214883;
74 aPoints( 1, 10 ) = 0.45891140214883;
75 aPoints( 2, 10 ) = 0.87508982900993;
76
77 aPoints( 0, 11 ) = 0.45891140214883;
78 aPoints( 1, 11 ) = 0.45891140214883;
79 aPoints( 2, 11 ) = -0.87508982900993;
80
81 aPoints( 0, 12 ) = 0.08217719570233;
82 aPoints( 1, 12 ) = 0.45891140214883;
83 aPoints( 2, 12 ) = 0.87508982900993;
84
85 aPoints( 0, 13 ) = 0.08217719570233;
86 aPoints( 1, 13 ) = 0.45891140214883;
87 aPoints( 2, 13 ) = -0.87508982900993;
88
89 aPoints( 0, 14 ) = 0.17213427575152;
90 aPoints( 1, 14 ) = 0.65573144849697;
91 aPoints( 2, 14 ) = 0.39891578047339;
92
93 aPoints( 0, 15 ) = 0.17213427575152;
94 aPoints( 1, 15 ) = 0.65573144849697;
95 aPoints( 2, 15 ) = -0.39891578047339;
96
97 aPoints( 0, 16 ) = 0.17213427575152;
98 aPoints( 1, 16 ) = 0.17213427575152;
99 aPoints( 2, 16 ) = 0.39891578047339;
100
101 aPoints( 0, 17 ) = 0.17213427575152;
102 aPoints( 1, 17 ) = 0.17213427575152;
103 aPoints( 2, 17 ) = -0.39891578047339;
104
105 aPoints( 0, 18 ) = 0.65573144849697;
106 aPoints( 1, 18 ) = 0.17213427575152;
107 aPoints( 2, 18 ) = 0.39891578047339;
108
109 aPoints( 0, 19 ) = 0.65573144849697;
110 aPoints( 1, 19 ) = 0.17213427575152;
111 aPoints( 2, 19 ) = -0.39891578047339;
112
113 aPoints( 0, 20 ) = 0.05182054676067;
114 aPoints( 1, 20 ) = 0.89635890647866;
115 aPoints( 2, 20 ) = 0.07604428020899;
116
117 aPoints( 0, 21 ) = 0.05182054676067;
118 aPoints( 1, 21 ) = 0.89635890647866;
119 aPoints( 2, 21 ) = -0.07604428020899;
120
121 aPoints( 0, 22 ) = 0.05182054676067;
122 aPoints( 1, 22 ) = 0.05182054676067;
123 aPoints( 2, 22 ) = 0.07604428020899;
124
125 aPoints( 0, 23 ) = 0.05182054676067;
126 aPoints( 1, 23 ) = 0.05182054676067;
127 aPoints( 2, 23 ) = -0.07604428020899;
128
129 aPoints( 0, 24 ) = 0.89635890647866;
130 aPoints( 1, 24 ) = 0.05182054676067;
131 aPoints( 2, 24 ) = 0.07604428020899;
132
133 aPoints( 0, 25 ) = 0.89635890647866;
134 aPoints( 1, 25 ) = 0.05182054676067;
135 aPoints( 2, 25 ) = -0.07604428020899;
136
137 aPoints( 0, 26 ) = 0.04905592580168;
138 aPoints( 1, 26 ) = 0.90188814839664;
139 aPoints( 2, 26 ) = 0.83813453082600;
140
141 aPoints( 0, 27 ) = 0.04905592580168;
142 aPoints( 1, 27 ) = 0.90188814839664;
143 aPoints( 2, 27 ) = -0.83813453082600;
144
145 aPoints( 0, 28 ) = 0.04905592580168;
146 aPoints( 1, 28 ) = 0.04905592580168;
147 aPoints( 2, 28 ) = 0.83813453082600;
148
149 aPoints( 0, 29 ) = 0.04905592580168;
150 aPoints( 1, 29 ) = 0.04905592580168;
151 aPoints( 2, 29 ) = -0.83813453082600;
152
153 aPoints( 0, 30 ) = 0.90188814839664;
154 aPoints( 1, 30 ) = 0.04905592580168;
155 aPoints( 2, 30 ) = 0.83813453082600;
156
157 aPoints( 0, 31 ) = 0.90188814839664;
158 aPoints( 1, 31 ) = 0.04905592580168;
159 aPoints( 2, 31 ) = -0.83813453082600;
160
161 aPoints( 0, 32 ) = 0.16556303163007;
162 aPoints( 1, 32 ) = 0.66887393673986;
163 aPoints( 2, 32 ) = 0.96242296786128;
164
165 aPoints( 0, 33 ) = 0.16556303163007;
166 aPoints( 1, 33 ) = 0.66887393673986;
167 aPoints( 2, 33 ) = -0.96242296786128;
168
169 aPoints( 0, 34 ) = 0.16556303163007;
170 aPoints( 1, 34 ) = 0.16556303163007;
171 aPoints( 2, 34 ) = 0.96242296786128;
172
173 aPoints( 0, 35 ) = 0.16556303163007;
174 aPoints( 1, 35 ) = 0.16556303163007;
175 aPoints( 2, 35 ) = -0.96242296786128;
176
177 aPoints( 0, 36 ) = 0.66887393673986;
178 aPoints( 1, 36 ) = 0.16556303163007;
179 aPoints( 2, 36 ) = 0.96242296786128;
180
181 aPoints( 0, 37 ) = 0.66887393673986;
182 aPoints( 1, 37 ) = 0.16556303163007;
183 aPoints( 2, 37 ) = -0.96242296786128;
184
185 aPoints( 0, 38 ) = 0.72851469748201;
186 aPoints( 1, 38 ) = 0.00843554067391;
187 aPoints( 2, 38 ) = 0.57735026919296;
188
189 aPoints( 0, 39 ) = 0.72851469748201;
190 aPoints( 1, 39 ) = 0.00843554067391;
191 aPoints( 2, 39 ) = -0.57735026919296;
192
193 aPoints( 0, 40 ) = 0.26304976184409;
194 aPoints( 1, 40 ) = 0.00843554067391;
195 aPoints( 2, 40 ) = 0.57735026919296;
196
197 aPoints( 0, 41 ) = 0.26304976184409;
198 aPoints( 1, 41 ) = 0.00843554067391;
199 aPoints( 2, 41 ) = -0.57735026919296;
200
201 aPoints( 0, 42 ) = 0.00843554067391;
202 aPoints( 1, 42 ) = 0.26304976184409;
203 aPoints( 2, 42 ) = 0.57735026919296;
204
205 aPoints( 0, 43 ) = 0.00843554067391;
206 aPoints( 1, 43 ) = 0.26304976184409;
207 aPoints( 2, 43 ) = -0.57735026919296;
208
209 aPoints( 0, 44 ) = 0.00843554067391;
210 aPoints( 1, 44 ) = 0.72851469748201;
211 aPoints( 2, 44 ) = 0.57735026919296;
212
213 aPoints( 0, 45 ) = 0.00843554067391;
214 aPoints( 1, 45 ) = 0.72851469748201;
215 aPoints( 2, 45 ) = -0.57735026919296;
216
217 aPoints( 0, 46 ) = 0.26304976184409;
218 aPoints( 1, 46 ) = 0.72851469748201;
219 aPoints( 2, 46 ) = 0.57735026919296;
220
221 aPoints( 0, 47 ) = 0.26304976184409;
222 aPoints( 1, 47 ) = 0.72851469748201;
223 aPoints( 2, 47 ) = -0.57735026919296;
224
225 aPoints( 0, 48 ) = 0.72851469748201;
226 aPoints( 1, 48 ) = 0.26304976184409;
227 aPoints( 2, 48 ) = 0.57735026919296;
228
229 aPoints( 0, 49 ) = 0.72851469748201;
230 aPoints( 1, 49 ) = 0.26304976184409;
231 aPoints( 2, 49 ) = -0.57735026919296;
232
233 aWeights.set_size( 50 );
234 aWeights( 0 ) = 0.01615571966122;
235 aWeights( 1 ) = 0.05320377087767;
236 aWeights( 2 ) = 0.05320377087767;
237 aWeights( 3 ) = 0.05320377087767;
238 aWeights( 4 ) = 0.05571482717252;
239 aWeights( 5 ) = 0.05571482717252;
240 aWeights( 6 ) = 0.00829184780739;
241 aWeights( 7 ) = 0.00829184780739;
242 aWeights( 8 ) = 0.02092087579991;
243 aWeights( 9 ) = 0.02092087579991;
244 aWeights( 10 ) = 0.02092087579991;
245 aWeights( 11 ) = 0.02092087579991;
246 aWeights( 12 ) = 0.02092087579991;
247 aWeights( 13 ) = 0.02092087579991;
248 aWeights( 14 ) = 0.04004467109788;
249 aWeights( 15 ) = 0.04004467109788;
250 aWeights( 16 ) = 0.04004467109788;
251 aWeights( 17 ) = 0.04004467109788;
252 aWeights( 18 ) = 0.04004467109788;
253 aWeights( 19 ) = 0.04004467109788;
254 aWeights( 20 ) = 0.00891956497473;
255 aWeights( 21 ) = 0.00891956497473;
256 aWeights( 22 ) = 0.00891956497473;
257 aWeights( 23 ) = 0.00891956497473;
258 aWeights( 24 ) = 0.00891956497473;
259 aWeights( 25 ) = 0.00891956497473;
260 aWeights( 26 ) = 0.00730764244975;
261 aWeights( 27 ) = 0.00730764244975;
262 aWeights( 28 ) = 0.00730764244975;
263 aWeights( 29 ) = 0.00730764244975;
264 aWeights( 30 ) = 0.00730764244975;
265 aWeights( 31 ) = 0.00730764244975;
266 aWeights( 32 ) = 0.01159992033528;
267 aWeights( 33 ) = 0.01159992033528;
268 aWeights( 34 ) = 0.01159992033528;
269 aWeights( 35 ) = 0.01159992033528;
270 aWeights( 36 ) = 0.01159992033528;
271 aWeights( 37 ) = 0.01159992033528;
272 aWeights( 38 ) = 0.01362196415005;
273 aWeights( 39 ) = 0.01362196415005;
274 aWeights( 40 ) = 0.01362196415005;
275 aWeights( 41 ) = 0.01362196415005;
276 aWeights( 42 ) = 0.01362196415005;
277 aWeights( 43 ) = 0.01362196415005;
278 aWeights( 44 ) = 0.01362196415005;
279 aWeights( 45 ) = 0.01362196415005;
280 aWeights( 46 ) = 0.01362196415005;
281 aWeights( 47 ) = 0.01362196415005;
282 aWeights( 48 ) = 0.01362196415005;
283 aWeights( 49 ) = 0.01362196415005;
284 }
285// ----------------------------------------------------------------------------
286 } /* namespace integration */
287} /* namespace belfem */
288#endif //BELFEM_FN_INTPOINTS_GAUSS_PENTA50_HPP
cl_Matrix.hpp
cl_Vector.hpp
belfem::Matrix::set_size
void set_size(const size_t aNumRows, const size_t aNumCols)
Definition cl_AR_Matrix.hpp:186
belfem::Vector::set_size
void set_size(const size_t aNumRows)
change the size of the vector
Definition cl_AR_Vector.hpp:237
belfem::integration
Definition fn_intpoints.cpp:576
belfem::integration::gauss_penta50
void gauss_penta50(Vector< real > &aWeights, Matrix< real > &aPoints)
Definition fn_intpoints_gauss_penta50.hpp:27
belfem
USER GUIDES:
Definition cl_Capacitor.cpp:16