BELFEM 0.9.0
Berkeley Lab Finite Element Framework
Loading...
Searching...
No Matches
fn_intpoints_gauss_tet165.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_TET165_HPP
13#define BELFEM_FN_INTPOINTS_GAUSS_TET165_HPP
14
15#include "typedefs.hpp"
16#include "cl_Vector.hpp"
17#include "cl_Matrix.hpp"
18#include "fn_sum.hpp"
19
20namespace belfem
21{
22 namespace integration
23 {
24// ----------------------------------------------------------------------------
25
34 inline void
35 gauss_tet165(
36 Vector< real > & aWeights,
37 Matrix< real > & aPoints )
38 {
39 aPoints.set_size( 4, 165 );
40
41 aPoints( 0, 0 ) = 0.25;
42 aPoints( 1, 0 ) = 0.25;
43 aPoints( 2, 0 ) = 0.25;
44
45 aPoints( 0, 1 ) = 0.02617103277295219;
46 aPoints( 1, 1 ) = 0.02617103277295219;
47 aPoints( 2, 1 ) = 0.02617103277295219;
48
49 aPoints( 0, 2 ) = 0.1006862685057424;
50 aPoints( 1, 2 ) = 0.1006862685057424;
51 aPoints( 2, 2 ) = 0.1006862685057424;
52
53 aPoints( 0, 3 ) = 0.02617103277295219;
54 aPoints( 1, 3 ) = 0.02617103277295219;
55 aPoints( 2, 3 ) = 0.9214869016811434;
56
57 aPoints( 0, 4 ) = 0.1006862685057424;
58 aPoints( 1, 4 ) = 0.1006862685057424;
59 aPoints( 2, 4 ) = 0.6979411944827728;
60
61 aPoints( 0, 5 ) = 0.02617103277295219;
62 aPoints( 1, 5 ) = 0.9214869016811434;
63 aPoints( 2, 5 ) = 0.02617103277295219;
64
65 aPoints( 0, 6 ) = 0.1006862685057424;
66 aPoints( 1, 6 ) = 0.6979411944827728;
67 aPoints( 2, 6 ) = 0.1006862685057424;
68
69 aPoints( 0, 7 ) = 0.9214869016811434;
70 aPoints( 1, 7 ) = 0.02617103277295219;
71 aPoints( 2, 7 ) = 0.02617103277295219;
72
73 aPoints( 0, 8 ) = 0.6979411944827728;
74 aPoints( 1, 8 ) = 0.1006862685057424;
75 aPoints( 2, 8 ) = 0.1006862685057424;
76
77 aPoints( 0, 9 ) = 0.01375236581886419;
78 aPoints( 1, 9 ) = 0.01375236581886419;
79 aPoints( 2, 9 ) = 0.4862476341811358;
80
81 aPoints( 0, 10 ) = 0.3881547512699796;
82 aPoints( 1, 10 ) = 0.3881547512699796;
83 aPoints( 2, 10 ) = 0.11184524873002039;
84
85 aPoints( 0, 11 ) = 0.01375236581886419;
86 aPoints( 1, 11 ) = 0.4862476341811358;
87 aPoints( 2, 11 ) = 0.01375236581886419;
88
89 aPoints( 0, 12 ) = 0.3881547512699796;
90 aPoints( 1, 12 ) = 0.11184524873002039;
91 aPoints( 2, 12 ) = 0.3881547512699796;
92
93 aPoints( 0, 13 ) = 0.4862476341811358;
94 aPoints( 1, 13 ) = 0.01375236581886419;
95 aPoints( 2, 13 ) = 0.01375236581886419;
96
97 aPoints( 0, 14 ) = 0.11184524873002039;
98 aPoints( 1, 14 ) = 0.3881547512699796;
99 aPoints( 2, 14 ) = 0.3881547512699796;
100
101 aPoints( 0, 15 ) = 0.01375236581886419;
102 aPoints( 1, 15 ) = 0.4862476341811358;
103 aPoints( 2, 15 ) = 0.4862476341811358;
104
105 aPoints( 0, 16 ) = 0.3881547512699796;
106 aPoints( 1, 16 ) = 0.11184524873002039;
107 aPoints( 2, 16 ) = 0.11184524873002039;
108
109 aPoints( 0, 17 ) = 0.4862476341811358;
110 aPoints( 1, 17 ) = 0.01375236581886419;
111 aPoints( 2, 17 ) = 0.4862476341811358;
112
113 aPoints( 0, 18 ) = 0.11184524873002039;
114 aPoints( 1, 18 ) = 0.3881547512699796;
115 aPoints( 2, 18 ) = 0.11184524873002039;
116
117 aPoints( 0, 19 ) = 0.4862476341811358;
118 aPoints( 1, 19 ) = 0.4862476341811358;
119 aPoints( 2, 19 ) = 0.01375236581886419;
120
121 aPoints( 0, 20 ) = 0.11184524873002039;
122 aPoints( 1, 20 ) = 0.11184524873002039;
123 aPoints( 2, 20 ) = 0.3881547512699796;
124
125 aPoints( 0, 21 ) = 0.02388787886656957;
126 aPoints( 1, 21 ) = 0.02388787886656957;
127 aPoints( 2, 21 ) = 0.8705443116154074;
128
129 aPoints( 0, 22 ) = 0.01903349443211699;
130 aPoints( 1, 22 ) = 0.01903349443211699;
131 aPoints( 2, 22 ) = 0.7970471485218408;
132
133 aPoints( 0, 23 ) = 0.100859667399371;
134 aPoints( 1, 23 ) = 0.100859667399371;
135 aPoints( 2, 23 ) = 0.01436337595093029;
136
137 aPoints( 0, 24 ) = 0.02434906451212993;
138 aPoints( 1, 24 ) = 0.02434906451212993;
139 aPoints( 2, 24 ) = 0.6322506459263664;
140
141 aPoints( 0, 25 ) = 0.3718666283348195;
142 aPoints( 1, 25 ) = 0.3718666283348195;
143 aPoints( 2, 25 ) = 0.2352758046590782;
144
145 aPoints( 0, 26 ) = 0.2247178848749781;
146 aPoints( 1, 26 ) = 0.2247178848749781;
147 aPoints( 2, 26 ) = 0.02248198827546399;
148
149 aPoints( 0, 27 ) = 0.1171336251346634;
150 aPoints( 1, 27 ) = 0.1171336251346634;
151 aPoints( 2, 27 ) = 0.5212126864871611;
152
153 aPoints( 0, 28 ) = 0.2516437633635237;
154 aPoints( 1, 28 ) = 0.2516437633635237;
155 aPoints( 2, 28 ) = 0.1073345337250223;
156
157 aPoints( 0, 29 ) = 0.02388787886656957;
158 aPoints( 1, 29 ) = 0.8705443116154074;
159 aPoints( 2, 29 ) = 0.02388787886656957;
160
161 aPoints( 0, 30 ) = 0.01903349443211699;
162 aPoints( 1, 30 ) = 0.7970471485218408;
163 aPoints( 2, 30 ) = 0.01903349443211699;
164
165 aPoints( 0, 31 ) = 0.100859667399371;
166 aPoints( 1, 31 ) = 0.01436337595093029;
167 aPoints( 2, 31 ) = 0.100859667399371;
168
169 aPoints( 0, 32 ) = 0.02434906451212993;
170 aPoints( 1, 32 ) = 0.6322506459263664;
171 aPoints( 2, 32 ) = 0.02434906451212993;
172
173 aPoints( 0, 33 ) = 0.3718666283348195;
174 aPoints( 1, 33 ) = 0.2352758046590782;
175 aPoints( 2, 33 ) = 0.3718666283348195;
176
177 aPoints( 0, 34 ) = 0.2247178848749781;
178 aPoints( 1, 34 ) = 0.02248198827546399;
179 aPoints( 2, 34 ) = 0.2247178848749781;
180
181 aPoints( 0, 35 ) = 0.1171336251346634;
182 aPoints( 1, 35 ) = 0.5212126864871611;
183 aPoints( 2, 35 ) = 0.1171336251346634;
184
185 aPoints( 0, 36 ) = 0.2516437633635237;
186 aPoints( 1, 36 ) = 0.1073345337250223;
187 aPoints( 2, 36 ) = 0.2516437633635237;
188
189 aPoints( 0, 37 ) = 0.8705443116154074;
190 aPoints( 1, 37 ) = 0.02388787886656957;
191 aPoints( 2, 37 ) = 0.02388787886656957;
192
193 aPoints( 0, 38 ) = 0.7970471485218408;
194 aPoints( 1, 38 ) = 0.01903349443211699;
195 aPoints( 2, 38 ) = 0.01903349443211699;
196
197 aPoints( 0, 39 ) = 0.01436337595093029;
198 aPoints( 1, 39 ) = 0.100859667399371;
199 aPoints( 2, 39 ) = 0.100859667399371;
200
201 aPoints( 0, 40 ) = 0.6322506459263664;
202 aPoints( 1, 40 ) = 0.02434906451212993;
203 aPoints( 2, 40 ) = 0.02434906451212993;
204
205 aPoints( 0, 41 ) = 0.2352758046590782;
206 aPoints( 1, 41 ) = 0.3718666283348195;
207 aPoints( 2, 41 ) = 0.3718666283348195;
208
209 aPoints( 0, 42 ) = 0.02248198827546399;
210 aPoints( 1, 42 ) = 0.2247178848749781;
211 aPoints( 2, 42 ) = 0.2247178848749781;
212
213 aPoints( 0, 43 ) = 0.5212126864871611;
214 aPoints( 1, 43 ) = 0.1171336251346634;
215 aPoints( 2, 43 ) = 0.1171336251346634;
216
217 aPoints( 0, 44 ) = 0.1073345337250223;
218 aPoints( 1, 44 ) = 0.2516437633635237;
219 aPoints( 2, 44 ) = 0.2516437633635237;
220
221 aPoints( 0, 45 ) = 0.02388787886656957;
222 aPoints( 1, 45 ) = 0.8705443116154074;
223 aPoints( 2, 45 ) = 0.08167993065145351;
224
225 aPoints( 0, 46 ) = 0.01903349443211699;
226 aPoints( 1, 46 ) = 0.7970471485218408;
227 aPoints( 2, 46 ) = 0.16488586261392524;
228
229 aPoints( 0, 47 ) = 0.100859667399371;
230 aPoints( 1, 47 ) = 0.01436337595093029;
231 aPoints( 2, 47 ) = 0.7839172892503277;
232
233 aPoints( 0, 48 ) = 0.02434906451212993;
234 aPoints( 1, 48 ) = 0.6322506459263664;
235 aPoints( 2, 48 ) = 0.3190512250493738;
236
237 aPoints( 0, 49 ) = 0.3718666283348195;
238 aPoints( 1, 49 ) = 0.2352758046590782;
239 aPoints( 2, 49 ) = 0.02099093867128282;
240
241 aPoints( 0, 50 ) = 0.2247178848749781;
242 aPoints( 1, 50 ) = 0.02248198827546399;
243 aPoints( 2, 50 ) = 0.5280822419745798;
244
245 aPoints( 0, 51 ) = 0.1171336251346634;
246 aPoints( 1, 51 ) = 0.5212126864871611;
247 aPoints( 2, 51 ) = 0.24452006324351205;
248
249 aPoints( 0, 52 ) = 0.2516437633635237;
250 aPoints( 1, 52 ) = 0.1073345337250223;
251 aPoints( 2, 52 ) = 0.3893779395479303;
252
253 aPoints( 0, 53 ) = 0.8705443116154074;
254 aPoints( 1, 53 ) = 0.02388787886656957;
255 aPoints( 2, 53 ) = 0.08167993065145351;
256
257 aPoints( 0, 54 ) = 0.7970471485218408;
258 aPoints( 1, 54 ) = 0.01903349443211699;
259 aPoints( 2, 54 ) = 0.16488586261392524;
260
261 aPoints( 0, 55 ) = 0.01436337595093029;
262 aPoints( 1, 55 ) = 0.100859667399371;
263 aPoints( 2, 55 ) = 0.7839172892503277;
264
265 aPoints( 0, 56 ) = 0.6322506459263664;
266 aPoints( 1, 56 ) = 0.02434906451212993;
267 aPoints( 2, 56 ) = 0.3190512250493738;
268
269 aPoints( 0, 57 ) = 0.2352758046590782;
270 aPoints( 1, 57 ) = 0.3718666283348195;
271 aPoints( 2, 57 ) = 0.02099093867128282;
272
273 aPoints( 0, 58 ) = 0.02248198827546399;
274 aPoints( 1, 58 ) = 0.2247178848749781;
275 aPoints( 2, 58 ) = 0.5280822419745798;
276
277 aPoints( 0, 59 ) = 0.5212126864871611;
278 aPoints( 1, 59 ) = 0.1171336251346634;
279 aPoints( 2, 59 ) = 0.24452006324351205;
280
281 aPoints( 0, 60 ) = 0.1073345337250223;
282 aPoints( 1, 60 ) = 0.2516437633635237;
283 aPoints( 2, 60 ) = 0.3893779395479303;
284
285 aPoints( 0, 61 ) = 0.8705443116154074;
286 aPoints( 1, 61 ) = 0.08167993065145351;
287 aPoints( 2, 61 ) = 0.02388787886656957;
288
289 aPoints( 0, 62 ) = 0.7970471485218408;
290 aPoints( 1, 62 ) = 0.16488586261392524;
291 aPoints( 2, 62 ) = 0.01903349443211699;
292
293 aPoints( 0, 63 ) = 0.01436337595093029;
294 aPoints( 1, 63 ) = 0.7839172892503277;
295 aPoints( 2, 63 ) = 0.100859667399371;
296
297 aPoints( 0, 64 ) = 0.6322506459263664;
298 aPoints( 1, 64 ) = 0.3190512250493738;
299 aPoints( 2, 64 ) = 0.02434906451212993;
300
301 aPoints( 0, 65 ) = 0.2352758046590782;
302 aPoints( 1, 65 ) = 0.02099093867128282;
303 aPoints( 2, 65 ) = 0.3718666283348195;
304
305 aPoints( 0, 66 ) = 0.02248198827546399;
306 aPoints( 1, 66 ) = 0.5280822419745798;
307 aPoints( 2, 66 ) = 0.2247178848749781;
308
309 aPoints( 0, 67 ) = 0.5212126864871611;
310 aPoints( 1, 67 ) = 0.24452006324351205;
311 aPoints( 2, 67 ) = 0.1171336251346634;
312
313 aPoints( 0, 68 ) = 0.1073345337250223;
314 aPoints( 1, 68 ) = 0.3893779395479303;
315 aPoints( 2, 68 ) = 0.2516437633635237;
316
317 aPoints( 0, 69 ) = 0.02388787886656957;
318 aPoints( 1, 69 ) = 0.02388787886656957;
319 aPoints( 2, 69 ) = 0.08167993065145351;
320
321 aPoints( 0, 70 ) = 0.01903349443211699;
322 aPoints( 1, 70 ) = 0.01903349443211699;
323 aPoints( 2, 70 ) = 0.16488586261392524;
324
325 aPoints( 0, 71 ) = 0.100859667399371;
326 aPoints( 1, 71 ) = 0.100859667399371;
327 aPoints( 2, 71 ) = 0.7839172892503277;
328
329 aPoints( 0, 72 ) = 0.02434906451212993;
330 aPoints( 1, 72 ) = 0.02434906451212993;
331 aPoints( 2, 72 ) = 0.3190512250493738;
332
333 aPoints( 0, 73 ) = 0.3718666283348195;
334 aPoints( 1, 73 ) = 0.3718666283348195;
335 aPoints( 2, 73 ) = 0.02099093867128282;
336
337 aPoints( 0, 74 ) = 0.2247178848749781;
338 aPoints( 1, 74 ) = 0.2247178848749781;
339 aPoints( 2, 74 ) = 0.5280822419745798;
340
341 aPoints( 0, 75 ) = 0.1171336251346634;
342 aPoints( 1, 75 ) = 0.1171336251346634;
343 aPoints( 2, 75 ) = 0.24452006324351205;
344
345 aPoints( 0, 76 ) = 0.2516437633635237;
346 aPoints( 1, 76 ) = 0.2516437633635237;
347 aPoints( 2, 76 ) = 0.3893779395479303;
348
349 aPoints( 0, 77 ) = 0.02388787886656957;
350 aPoints( 1, 77 ) = 0.08167993065145351;
351 aPoints( 2, 77 ) = 0.02388787886656957;
352
353 aPoints( 0, 78 ) = 0.01903349443211699;
354 aPoints( 1, 78 ) = 0.16488586261392524;
355 aPoints( 2, 78 ) = 0.01903349443211699;
356
357 aPoints( 0, 79 ) = 0.100859667399371;
358 aPoints( 1, 79 ) = 0.7839172892503277;
359 aPoints( 2, 79 ) = 0.100859667399371;
360
361 aPoints( 0, 80 ) = 0.02434906451212993;
362 aPoints( 1, 80 ) = 0.3190512250493738;
363 aPoints( 2, 80 ) = 0.02434906451212993;
364
365 aPoints( 0, 81 ) = 0.3718666283348195;
366 aPoints( 1, 81 ) = 0.02099093867128282;
367 aPoints( 2, 81 ) = 0.3718666283348195;
368
369 aPoints( 0, 82 ) = 0.2247178848749781;
370 aPoints( 1, 82 ) = 0.5280822419745798;
371 aPoints( 2, 82 ) = 0.2247178848749781;
372
373 aPoints( 0, 83 ) = 0.1171336251346634;
374 aPoints( 1, 83 ) = 0.24452006324351205;
375 aPoints( 2, 83 ) = 0.1171336251346634;
376
377 aPoints( 0, 84 ) = 0.2516437633635237;
378 aPoints( 1, 84 ) = 0.3893779395479303;
379 aPoints( 2, 84 ) = 0.2516437633635237;
380
381 aPoints( 0, 85 ) = 0.08167993065145351;
382 aPoints( 1, 85 ) = 0.02388787886656957;
383 aPoints( 2, 85 ) = 0.02388787886656957;
384
385 aPoints( 0, 86 ) = 0.16488586261392524;
386 aPoints( 1, 86 ) = 0.01903349443211699;
387 aPoints( 2, 86 ) = 0.01903349443211699;
388
389 aPoints( 0, 87 ) = 0.7839172892503277;
390 aPoints( 1, 87 ) = 0.100859667399371;
391 aPoints( 2, 87 ) = 0.100859667399371;
392
393 aPoints( 0, 88 ) = 0.3190512250493738;
394 aPoints( 1, 88 ) = 0.02434906451212993;
395 aPoints( 2, 88 ) = 0.02434906451212993;
396
397 aPoints( 0, 89 ) = 0.02099093867128282;
398 aPoints( 1, 89 ) = 0.3718666283348195;
399 aPoints( 2, 89 ) = 0.3718666283348195;
400
401 aPoints( 0, 90 ) = 0.5280822419745798;
402 aPoints( 1, 90 ) = 0.2247178848749781;
403 aPoints( 2, 90 ) = 0.2247178848749781;
404
405 aPoints( 0, 91 ) = 0.24452006324351205;
406 aPoints( 1, 91 ) = 0.1171336251346634;
407 aPoints( 2, 91 ) = 0.1171336251346634;
408
409 aPoints( 0, 92 ) = 0.3893779395479303;
410 aPoints( 1, 92 ) = 0.2516437633635237;
411 aPoints( 2, 92 ) = 0.2516437633635237;
412
413 aPoints( 0, 93 ) = 0.02388787886656957;
414 aPoints( 1, 93 ) = 0.08167993065145351;
415 aPoints( 2, 93 ) = 0.8705443116154074;
416
417 aPoints( 0, 94 ) = 0.01903349443211699;
418 aPoints( 1, 94 ) = 0.16488586261392524;
419 aPoints( 2, 94 ) = 0.7970471485218408;
420
421 aPoints( 0, 95 ) = 0.100859667399371;
422 aPoints( 1, 95 ) = 0.7839172892503277;
423 aPoints( 2, 95 ) = 0.01436337595093029;
424
425 aPoints( 0, 96 ) = 0.02434906451212993;
426 aPoints( 1, 96 ) = 0.3190512250493738;
427 aPoints( 2, 96 ) = 0.6322506459263664;
428
429 aPoints( 0, 97 ) = 0.3718666283348195;
430 aPoints( 1, 97 ) = 0.02099093867128282;
431 aPoints( 2, 97 ) = 0.2352758046590782;
432
433 aPoints( 0, 98 ) = 0.2247178848749781;
434 aPoints( 1, 98 ) = 0.5280822419745798;
435 aPoints( 2, 98 ) = 0.02248198827546399;
436
437 aPoints( 0, 99 ) = 0.1171336251346634;
438 aPoints( 1, 99 ) = 0.24452006324351205;
439 aPoints( 2, 99 ) = 0.5212126864871611;
440
441 aPoints( 0, 100 ) = 0.2516437633635237;
442 aPoints( 1, 100 ) = 0.3893779395479303;
443 aPoints( 2, 100 ) = 0.1073345337250223;
444
445 aPoints( 0, 101 ) = 0.08167993065145351;
446 aPoints( 1, 101 ) = 0.02388787886656957;
447 aPoints( 2, 101 ) = 0.8705443116154074;
448
449 aPoints( 0, 102 ) = 0.16488586261392524;
450 aPoints( 1, 102 ) = 0.01903349443211699;
451 aPoints( 2, 102 ) = 0.7970471485218408;
452
453 aPoints( 0, 103 ) = 0.7839172892503277;
454 aPoints( 1, 103 ) = 0.100859667399371;
455 aPoints( 2, 103 ) = 0.01436337595093029;
456
457 aPoints( 0, 104 ) = 0.3190512250493738;
458 aPoints( 1, 104 ) = 0.02434906451212993;
459 aPoints( 2, 104 ) = 0.6322506459263664;
460
461 aPoints( 0, 105 ) = 0.02099093867128282;
462 aPoints( 1, 105 ) = 0.3718666283348195;
463 aPoints( 2, 105 ) = 0.2352758046590782;
464
465 aPoints( 0, 106 ) = 0.5280822419745798;
466 aPoints( 1, 106 ) = 0.2247178848749781;
467 aPoints( 2, 106 ) = 0.02248198827546399;
468
469 aPoints( 0, 107 ) = 0.24452006324351205;
470 aPoints( 1, 107 ) = 0.1171336251346634;
471 aPoints( 2, 107 ) = 0.5212126864871611;
472
473 aPoints( 0, 108 ) = 0.3893779395479303;
474 aPoints( 1, 108 ) = 0.2516437633635237;
475 aPoints( 2, 108 ) = 0.1073345337250223;
476
477 aPoints( 0, 109 ) = 0.08167993065145351;
478 aPoints( 1, 109 ) = 0.8705443116154074;
479 aPoints( 2, 109 ) = 0.02388787886656957;
480
481 aPoints( 0, 110 ) = 0.16488586261392524;
482 aPoints( 1, 110 ) = 0.7970471485218408;
483 aPoints( 2, 110 ) = 0.01903349443211699;
484
485 aPoints( 0, 111 ) = 0.7839172892503277;
486 aPoints( 1, 111 ) = 0.01436337595093029;
487 aPoints( 2, 111 ) = 0.100859667399371;
488
489 aPoints( 0, 112 ) = 0.3190512250493738;
490 aPoints( 1, 112 ) = 0.6322506459263664;
491 aPoints( 2, 112 ) = 0.02434906451212993;
492
493 aPoints( 0, 113 ) = 0.02099093867128282;
494 aPoints( 1, 113 ) = 0.2352758046590782;
495 aPoints( 2, 113 ) = 0.3718666283348195;
496
497 aPoints( 0, 114 ) = 0.5280822419745798;
498 aPoints( 1, 114 ) = 0.02248198827546399;
499 aPoints( 2, 114 ) = 0.2247178848749781;
500
501 aPoints( 0, 115 ) = 0.24452006324351205;
502 aPoints( 1, 115 ) = 0.5212126864871611;
503 aPoints( 2, 115 ) = 0.1171336251346634;
504
505 aPoints( 0, 116 ) = 0.3893779395479303;
506 aPoints( 1, 116 ) = 0.1073345337250223;
507 aPoints( 2, 116 ) = 0.2516437633635237;
508
509 aPoints( 0, 117 ) = 0.6716079793211197;
510 aPoints( 1, 117 ) = 0.2028653638555222;
511 aPoints( 2, 117 ) = 0.1012015873047946;
512
513 aPoints( 0, 118 ) = 0.5181755993267751;
514 aPoints( 1, 118 ) = 0.3567323365469498;
515 aPoints( 2, 118 ) = 0.1028384582772451;
516
517 aPoints( 0, 119 ) = 0.6716079793211197;
518 aPoints( 1, 119 ) = 0.2028653638555222;
519 aPoints( 2, 119 ) = 0.02432506951856349;
520
521 aPoints( 0, 120 ) = 0.5181755993267751;
522 aPoints( 1, 120 ) = 0.3567323365469498;
523 aPoints( 2, 120 ) = 0.022253605849029953;
524
525 aPoints( 0, 121 ) = 0.6716079793211197;
526 aPoints( 1, 121 ) = 0.1012015873047946;
527 aPoints( 2, 121 ) = 0.2028653638555222;
528
529 aPoints( 0, 122 ) = 0.5181755993267751;
530 aPoints( 1, 122 ) = 0.1028384582772451;
531 aPoints( 2, 122 ) = 0.3567323365469498;
532
533 aPoints( 0, 123 ) = 0.6716079793211197;
534 aPoints( 1, 123 ) = 0.1012015873047946;
535 aPoints( 2, 123 ) = 0.02432506951856349;
536
537 aPoints( 0, 124 ) = 0.5181755993267751;
538 aPoints( 1, 124 ) = 0.1028384582772451;
539 aPoints( 2, 124 ) = 0.022253605849029953;
540
541 aPoints( 0, 125 ) = 0.6716079793211197;
542 aPoints( 1, 125 ) = 0.02432506951856349;
543 aPoints( 2, 125 ) = 0.2028653638555222;
544
545 aPoints( 0, 126 ) = 0.5181755993267751;
546 aPoints( 1, 126 ) = 0.022253605849029953;
547 aPoints( 2, 126 ) = 0.3567323365469498;
548
549 aPoints( 0, 127 ) = 0.6716079793211197;
550 aPoints( 1, 127 ) = 0.02432506951856349;
551 aPoints( 2, 127 ) = 0.1012015873047946;
552
553 aPoints( 0, 128 ) = 0.5181755993267751;
554 aPoints( 1, 128 ) = 0.022253605849029953;
555 aPoints( 2, 128 ) = 0.1028384582772451;
556
557 aPoints( 0, 129 ) = 0.2028653638555222;
558 aPoints( 1, 129 ) = 0.6716079793211197;
559 aPoints( 2, 129 ) = 0.1012015873047946;
560
561 aPoints( 0, 130 ) = 0.3567323365469498;
562 aPoints( 1, 130 ) = 0.5181755993267751;
563 aPoints( 2, 130 ) = 0.1028384582772451;
564
565 aPoints( 0, 131 ) = 0.2028653638555222;
566 aPoints( 1, 131 ) = 0.6716079793211197;
567 aPoints( 2, 131 ) = 0.02432506951856349;
568
569 aPoints( 0, 132 ) = 0.3567323365469498;
570 aPoints( 1, 132 ) = 0.5181755993267751;
571 aPoints( 2, 132 ) = 0.022253605849029953;
572
573 aPoints( 0, 133 ) = 0.2028653638555222;
574 aPoints( 1, 133 ) = 0.1012015873047946;
575 aPoints( 2, 133 ) = 0.6716079793211197;
576
577 aPoints( 0, 134 ) = 0.3567323365469498;
578 aPoints( 1, 134 ) = 0.1028384582772451;
579 aPoints( 2, 134 ) = 0.5181755993267751;
580
581 aPoints( 0, 135 ) = 0.2028653638555222;
582 aPoints( 1, 135 ) = 0.1012015873047946;
583 aPoints( 2, 135 ) = 0.02432506951856349;
584
585 aPoints( 0, 136 ) = 0.3567323365469498;
586 aPoints( 1, 136 ) = 0.1028384582772451;
587 aPoints( 2, 136 ) = 0.022253605849029953;
588
589 aPoints( 0, 137 ) = 0.2028653638555222;
590 aPoints( 1, 137 ) = 0.02432506951856349;
591 aPoints( 2, 137 ) = 0.6716079793211197;
592
593 aPoints( 0, 138 ) = 0.3567323365469498;
594 aPoints( 1, 138 ) = 0.022253605849029953;
595 aPoints( 2, 138 ) = 0.5181755993267751;
596
597 aPoints( 0, 139 ) = 0.2028653638555222;
598 aPoints( 1, 139 ) = 0.02432506951856349;
599 aPoints( 2, 139 ) = 0.1012015873047946;
600
601 aPoints( 0, 140 ) = 0.3567323365469498;
602 aPoints( 1, 140 ) = 0.022253605849029953;
603 aPoints( 2, 140 ) = 0.1028384582772451;
604
605 aPoints( 0, 141 ) = 0.1012015873047946;
606 aPoints( 1, 141 ) = 0.6716079793211197;
607 aPoints( 2, 141 ) = 0.2028653638555222;
608
609 aPoints( 0, 142 ) = 0.1028384582772451;
610 aPoints( 1, 142 ) = 0.5181755993267751;
611 aPoints( 2, 142 ) = 0.3567323365469498;
612
613 aPoints( 0, 143 ) = 0.1012015873047946;
614 aPoints( 1, 143 ) = 0.6716079793211197;
615 aPoints( 2, 143 ) = 0.02432506951856349;
616
617 aPoints( 0, 144 ) = 0.1028384582772451;
618 aPoints( 1, 144 ) = 0.5181755993267751;
619 aPoints( 2, 144 ) = 0.022253605849029953;
620
621 aPoints( 0, 145 ) = 0.1012015873047946;
622 aPoints( 1, 145 ) = 0.2028653638555222;
623 aPoints( 2, 145 ) = 0.6716079793211197;
624
625 aPoints( 0, 146 ) = 0.1028384582772451;
626 aPoints( 1, 146 ) = 0.3567323365469498;
627 aPoints( 2, 146 ) = 0.5181755993267751;
628
629 aPoints( 0, 147 ) = 0.1012015873047946;
630 aPoints( 1, 147 ) = 0.2028653638555222;
631 aPoints( 2, 147 ) = 0.02432506951856349;
632
633 aPoints( 0, 148 ) = 0.1028384582772451;
634 aPoints( 1, 148 ) = 0.3567323365469498;
635 aPoints( 2, 148 ) = 0.022253605849029953;
636
637 aPoints( 0, 149 ) = 0.1012015873047946;
638 aPoints( 1, 149 ) = 0.02432506951856349;
639 aPoints( 2, 149 ) = 0.6716079793211197;
640
641 aPoints( 0, 150 ) = 0.1028384582772451;
642 aPoints( 1, 150 ) = 0.022253605849029953;
643 aPoints( 2, 150 ) = 0.5181755993267751;
644
645 aPoints( 0, 151 ) = 0.1012015873047946;
646 aPoints( 1, 151 ) = 0.02432506951856349;
647 aPoints( 2, 151 ) = 0.2028653638555222;
648
649 aPoints( 0, 152 ) = 0.1028384582772451;
650 aPoints( 1, 152 ) = 0.022253605849029953;
651 aPoints( 2, 152 ) = 0.3567323365469498;
652
653 aPoints( 0, 153 ) = 0.02432506951856349;
654 aPoints( 1, 153 ) = 0.6716079793211197;
655 aPoints( 2, 153 ) = 0.2028653638555222;
656
657 aPoints( 0, 154 ) = 0.022253605849029953;
658 aPoints( 1, 154 ) = 0.5181755993267751;
659 aPoints( 2, 154 ) = 0.3567323365469498;
660
661 aPoints( 0, 155 ) = 0.02432506951856349;
662 aPoints( 1, 155 ) = 0.6716079793211197;
663 aPoints( 2, 155 ) = 0.1012015873047946;
664
665 aPoints( 0, 156 ) = 0.022253605849029953;
666 aPoints( 1, 156 ) = 0.5181755993267751;
667 aPoints( 2, 156 ) = 0.1028384582772451;
668
669 aPoints( 0, 157 ) = 0.02432506951856349;
670 aPoints( 1, 157 ) = 0.2028653638555222;
671 aPoints( 2, 157 ) = 0.6716079793211197;
672
673 aPoints( 0, 158 ) = 0.022253605849029953;
674 aPoints( 1, 158 ) = 0.3567323365469498;
675 aPoints( 2, 158 ) = 0.5181755993267751;
676
677 aPoints( 0, 159 ) = 0.02432506951856349;
678 aPoints( 1, 159 ) = 0.2028653638555222;
679 aPoints( 2, 159 ) = 0.1012015873047946;
680
681 aPoints( 0, 160 ) = 0.022253605849029953;
682 aPoints( 1, 160 ) = 0.3567323365469498;
683 aPoints( 2, 160 ) = 0.1028384582772451;
684
685 aPoints( 0, 161 ) = 0.02432506951856349;
686 aPoints( 1, 161 ) = 0.1012015873047946;
687 aPoints( 2, 161 ) = 0.6716079793211197;
688
689 aPoints( 0, 162 ) = 0.022253605849029953;
690 aPoints( 1, 162 ) = 0.1028384582772451;
691 aPoints( 2, 162 ) = 0.5181755993267751;
692
693 aPoints( 0, 163 ) = 0.02432506951856349;
694 aPoints( 1, 163 ) = 0.1012015873047946;
695 aPoints( 2, 163 ) = 0.2028653638555222;
696
697 aPoints( 0, 164 ) = 0.022253605849029953;
698 aPoints( 1, 164 ) = 0.1028384582772451;
699 aPoints( 2, 164 ) = 0.3567323365469498;
700
701 for( uint k=0; k<165; ++k )
702 {
703 aPoints( 3, k ) = 1.0
704 - aPoints( 0, k )
705 - aPoints( 1, k )
706 - aPoints( 2, k );
707 }
708
709
710 aWeights.set_size( 165 );
711
712 aWeights( 0 ) = 0.004128263926563686;
713 aWeights( 1 ) = 0.0002449879785283559;
714 aWeights( 2 ) = 0.0016552570096826388;
715 aWeights( 3 ) = 0.0002449879785283559;
716 aWeights( 4 ) = 0.0016552570096826388;
717 aWeights( 5 ) = 0.0002449879785283559;
718 aWeights( 6 ) = 0.0016552570096826388;
719 aWeights( 7 ) = 0.0002449879785283559;
720 aWeights( 8 ) = 0.0016552570096826388;
721 aWeights( 9 ) = 0.00028659306196642;
722 aWeights( 10 ) = 0.0018861224110710562;
723 aWeights( 11 ) = 0.00028659306196642;
724 aWeights( 12 ) = 0.0018861224110710562;
725 aWeights( 13 ) = 0.00028659306196642;
726 aWeights( 14 ) = 0.0018861224110710562;
727 aWeights( 15 ) = 0.00028659306196642;
728 aWeights( 16 ) = 0.0018861224110710562;
729 aWeights( 17 ) = 0.00028659306196642;
730 aWeights( 18 ) = 0.0018861224110710562;
731 aWeights( 19 ) = 0.00028659306196642;
732 aWeights( 20 ) = 0.0018861224110710562;
733 aWeights( 21 ) = 0.00012119309860710957;
734 aWeights( 22 ) = 0.0003503365407403157;
735 aWeights( 23 ) = 0.00038979850792971176;
736 aWeights( 24 ) = 0.0005456570060766684;
737 aWeights( 25 ) = 0.0010806694342139658;
738 aWeights( 26 ) = 0.0011851247050948092;
739 aWeights( 27 ) = 0.0022576163852431725;
740 aWeights( 28 ) = 0.002303218527587786;
741 aWeights( 29 ) = 0.00012119309860710957;
742 aWeights( 30 ) = 0.0003503365407403157;
743 aWeights( 31 ) = 0.00038979850792971176;
744 aWeights( 32 ) = 0.0005456570060766684;
745 aWeights( 33 ) = 0.0010806694342139658;
746 aWeights( 34 ) = 0.0011851247050948092;
747 aWeights( 35 ) = 0.0022576163852431725;
748 aWeights( 36 ) = 0.002303218527587786;
749 aWeights( 37 ) = 0.00012119309860710957;
750 aWeights( 38 ) = 0.0003503365407403157;
751 aWeights( 39 ) = 0.00038979850792971176;
752 aWeights( 40 ) = 0.0005456570060766684;
753 aWeights( 41 ) = 0.0010806694342139658;
754 aWeights( 42 ) = 0.0011851247050948092;
755 aWeights( 43 ) = 0.0022576163852431725;
756 aWeights( 44 ) = 0.002303218527587786;
757 aWeights( 45 ) = 0.00012119309860710957;
758 aWeights( 46 ) = 0.0003503365407403157;
759 aWeights( 47 ) = 0.00038979850792971176;
760 aWeights( 48 ) = 0.0005456570060766684;
761 aWeights( 49 ) = 0.0010806694342139658;
762 aWeights( 50 ) = 0.0011851247050948092;
763 aWeights( 51 ) = 0.0022576163852431725;
764 aWeights( 52 ) = 0.002303218527587786;
765 aWeights( 53 ) = 0.00012119309860710957;
766 aWeights( 54 ) = 0.0003503365407403157;
767 aWeights( 55 ) = 0.00038979850792971176;
768 aWeights( 56 ) = 0.0005456570060766684;
769 aWeights( 57 ) = 0.0010806694342139658;
770 aWeights( 58 ) = 0.0011851247050948092;
771 aWeights( 59 ) = 0.0022576163852431725;
772 aWeights( 60 ) = 0.002303218527587786;
773 aWeights( 61 ) = 0.00012119309860710957;
774 aWeights( 62 ) = 0.0003503365407403157;
775 aWeights( 63 ) = 0.00038979850792971176;
776 aWeights( 64 ) = 0.0005456570060766684;
777 aWeights( 65 ) = 0.0010806694342139658;
778 aWeights( 66 ) = 0.0011851247050948092;
779 aWeights( 67 ) = 0.0022576163852431725;
780 aWeights( 68 ) = 0.002303218527587786;
781 aWeights( 69 ) = 0.00012119309860710957;
782 aWeights( 70 ) = 0.0003503365407403157;
783 aWeights( 71 ) = 0.00038979850792971176;
784 aWeights( 72 ) = 0.0005456570060766684;
785 aWeights( 73 ) = 0.0010806694342139658;
786 aWeights( 74 ) = 0.0011851247050948092;
787 aWeights( 75 ) = 0.0022576163852431725;
788 aWeights( 76 ) = 0.002303218527587786;
789 aWeights( 77 ) = 0.00012119309860710957;
790 aWeights( 78 ) = 0.0003503365407403157;
791 aWeights( 79 ) = 0.00038979850792971176;
792 aWeights( 80 ) = 0.0005456570060766684;
793 aWeights( 81 ) = 0.0010806694342139658;
794 aWeights( 82 ) = 0.0011851247050948092;
795 aWeights( 83 ) = 0.0022576163852431725;
796 aWeights( 84 ) = 0.002303218527587786;
797 aWeights( 85 ) = 0.00012119309860710957;
798 aWeights( 86 ) = 0.0003503365407403157;
799 aWeights( 87 ) = 0.00038979850792971176;
800 aWeights( 88 ) = 0.0005456570060766684;
801 aWeights( 89 ) = 0.0010806694342139658;
802 aWeights( 90 ) = 0.0011851247050948092;
803 aWeights( 91 ) = 0.0022576163852431725;
804 aWeights( 92 ) = 0.002303218527587786;
805 aWeights( 93 ) = 0.00012119309860710957;
806 aWeights( 94 ) = 0.0003503365407403157;
807 aWeights( 95 ) = 0.00038979850792971176;
808 aWeights( 96 ) = 0.0005456570060766684;
809 aWeights( 97 ) = 0.0010806694342139658;
810 aWeights( 98 ) = 0.0011851247050948092;
811 aWeights( 99 ) = 0.0022576163852431725;
812 aWeights( 100 ) = 0.002303218527587786;
813 aWeights( 101 ) = 0.00012119309860710957;
814 aWeights( 102 ) = 0.0003503365407403157;
815 aWeights( 103 ) = 0.00038979850792971176;
816 aWeights( 104 ) = 0.0005456570060766684;
817 aWeights( 105 ) = 0.0010806694342139658;
818 aWeights( 106 ) = 0.0011851247050948092;
819 aWeights( 107 ) = 0.0022576163852431725;
820 aWeights( 108 ) = 0.002303218527587786;
821 aWeights( 109 ) = 0.00012119309860710957;
822 aWeights( 110 ) = 0.0003503365407403157;
823 aWeights( 111 ) = 0.00038979850792971176;
824 aWeights( 112 ) = 0.0005456570060766684;
825 aWeights( 113 ) = 0.0010806694342139658;
826 aWeights( 114 ) = 0.0011851247050948092;
827 aWeights( 115 ) = 0.0022576163852431725;
828 aWeights( 116 ) = 0.002303218527587786;
829 aWeights( 117 ) = 0.000743794899989901;
830 aWeights( 118 ) = 0.0010519450784730495;
831 aWeights( 119 ) = 0.000743794899989901;
832 aWeights( 120 ) = 0.0010519450784730495;
833 aWeights( 121 ) = 0.000743794899989901;
834 aWeights( 122 ) = 0.0010519450784730495;
835 aWeights( 123 ) = 0.000743794899989901;
836 aWeights( 124 ) = 0.0010519450784730495;
837 aWeights( 125 ) = 0.000743794899989901;
838 aWeights( 126 ) = 0.0010519450784730495;
839 aWeights( 127 ) = 0.000743794899989901;
840 aWeights( 128 ) = 0.0010519450784730495;
841 aWeights( 129 ) = 0.000743794899989901;
842 aWeights( 130 ) = 0.0010519450784730495;
843 aWeights( 131 ) = 0.000743794899989901;
844 aWeights( 132 ) = 0.0010519450784730495;
845 aWeights( 133 ) = 0.000743794899989901;
846 aWeights( 134 ) = 0.0010519450784730495;
847 aWeights( 135 ) = 0.000743794899989901;
848 aWeights( 136 ) = 0.0010519450784730495;
849 aWeights( 137 ) = 0.000743794899989901;
850 aWeights( 138 ) = 0.0010519450784730495;
851 aWeights( 139 ) = 0.000743794899989901;
852 aWeights( 140 ) = 0.0010519450784730495;
853 aWeights( 141 ) = 0.000743794899989901;
854 aWeights( 142 ) = 0.0010519450784730495;
855 aWeights( 143 ) = 0.000743794899989901;
856 aWeights( 144 ) = 0.0010519450784730495;
857 aWeights( 145 ) = 0.000743794899989901;
858 aWeights( 146 ) = 0.0010519450784730495;
859 aWeights( 147 ) = 0.000743794899989901;
860 aWeights( 148 ) = 0.0010519450784730495;
861 aWeights( 149 ) = 0.000743794899989901;
862 aWeights( 150 ) = 0.0010519450784730495;
863 aWeights( 151 ) = 0.000743794899989901;
864 aWeights( 152 ) = 0.0010519450784730495;
865 aWeights( 153 ) = 0.000743794899989901;
866 aWeights( 154 ) = 0.0010519450784730495;
867 aWeights( 155 ) = 0.000743794899989901;
868 aWeights( 156 ) = 0.0010519450784730495;
869 aWeights( 157 ) = 0.000743794899989901;
870 aWeights( 158 ) = 0.0010519450784730495;
871 aWeights( 159 ) = 0.000743794899989901;
872 aWeights( 160 ) = 0.0010519450784730495;
873 aWeights( 161 ) = 0.000743794899989901;
874 aWeights( 162 ) = 0.0010519450784730495;
875 aWeights( 163 ) = 0.000743794899989901;
876 aWeights( 164 ) = 0.0010519450784730495;
877
878 aWeights( 0 ) = 0.0 ;
879 aWeights( 0 ) = 1./6. - sum( aWeights );
880 }
881
882// ----------------------------------------------------------------------------
883 } /* namespace integration */
884} /* end namespace belfem */
885
886#endif // BELFEM_FN_INTPOINTS_GAUSS_TET165_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
fn_sum.hpp
Sum of the entries of a vector.
belfem::integration
Definition fn_intpoints.cpp:576
belfem::integration::gauss_tet165
void gauss_tet165(Vector< real > &aWeights, Matrix< real > &aPoints)
11th order interpolation
Definition fn_intpoints_gauss_tet165.hpp:35
belfem
USER GUIDES:
Definition cl_Capacitor.cpp:16
belfem::uint
unsigned int uint
Definition typedefs.hpp:30
belfem::sum
auto sum(const Vector< T > &aA) -> decltype(sum(aA.vector_data()))
Definition fn_sum.hpp:49