BELFEM 0.9.0
Berkeley Lab Finite Element Framework
Loading...
Searching...
No Matches
fn_intpoints_gauss_penta11.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_PENTA11_HPP
13#define BELFEM_FN_INTPOINTS_GAUSS_PENTA11_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 inline void
27 gauss_penta11(
28 Vector< real > & aWeights,
29 Matrix< real > & aPoints )
30 {
31 aPoints.set_size( 3, 11 );
32
33 aPoints( 0, 0 ) = 0.06268838027601;
34 aPoints( 1, 0 ) = 0.46865580986200;
35 aPoints( 2, 0 ) = 0.00000000000000;
36
37 aPoints( 0, 1 ) = 0.46865580986200;
38 aPoints( 1, 1 ) = 0.46865580986200;
39 aPoints( 2, 1 ) = 0.00000000000000;
40
41 aPoints( 0, 2 ) = 0.46865580986200;
42 aPoints( 1, 2 ) = 0.06268838027601;
43 aPoints( 2, 2 ) = 0.00000000000000;
44
45 aPoints( 0, 3 ) = 0.33333333333333;
46 aPoints( 1, 3 ) = 0.33333333333333;
47 aPoints( 2, 3 ) = -0.86686197400903;
48
49 aPoints( 0, 4 ) = 0.33333333333333;
50 aPoints( 1, 4 ) = 0.33333333333333;
51 aPoints( 2, 4 ) = 0.86686197400903;
52
53 aPoints( 0, 5 ) = 0.79851918840218;
54 aPoints( 1, 5 ) = 0.10074040579891;
55 aPoints( 2, 5 ) = -0.67563982368227;
56
57 aPoints( 0, 6 ) = 0.79851918840218;
58 aPoints( 1, 6 ) = 0.10074040579891;
59 aPoints( 2, 6 ) = 0.67563982368227;
60
61 aPoints( 0, 7 ) = 0.10074040579891;
62 aPoints( 1, 7 ) = 0.10074040579891;
63 aPoints( 2, 7 ) = -0.67563982368227;
64
65 aPoints( 0, 8 ) = 0.10074040579891;
66 aPoints( 1, 8 ) = 0.10074040579891;
67 aPoints( 2, 8 ) = 0.67563982368227;
68
69 aPoints( 0, 9 ) = 0.10074040579891;
70 aPoints( 1, 9 ) = 0.79851918840218;
71 aPoints( 2, 9 ) = -0.67563982368227;
72
73 aPoints( 0, 10 ) = 0.10074040579891;
74 aPoints( 1, 10 ) = 0.79851918840218;
75 aPoints( 2, 10 ) = 0.67563982368227;
76
77 aWeights.set_size( 11 );
78 aWeights( 0 ) = 0.13641461260548;
79 aWeights( 1 ) = 0.13641461260548;
80 aWeights( 2 ) = 0.13641461260548;
81 aWeights( 3 ) = 0.10791197481553;
82 aWeights( 4 ) = 0.10791197481553;
83 aWeights( 5 ) = 0.06248870209208;
84 aWeights( 6 ) = 0.06248870209208;
85 aWeights( 7 ) = 0.06248870209208;
86 aWeights( 8 ) = 0.06248870209208;
87 aWeights( 9 ) = 0.06248870209208;
88 aWeights( 10 ) = 0.06248870209208;
89 }
90
91// ----------------------------------------------------------------------------
92 } /* namespace integration */
93} /* namespace belfem */
94
95#endif //BELFEM_FN_INTPOINTS_GAUSS_PENTA11_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_penta11
void gauss_penta11(Vector< real > &aWeights, Matrix< real > &aPoints)
Definition fn_intpoints_gauss_penta11.hpp:27
belfem
USER GUIDES:
Definition cl_Capacitor.cpp:16