BELFEM 0.9.0
Berkeley Lab Finite Element Framework
Loading...
Searching...
No Matches
fn_intpoints_gauss_quad20.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_QUAD20_HPP
13#define BELFEM_FN_INTPOINTS_GAUSS_QUAD20_HPP
14#include "typedefs.hpp"
15#include "cl_Vector.hpp"
16#include "cl_Matrix.hpp"
17
18namespace belfem
19{
20 namespace integration
21 {
22// ----------------------------------------------------------------------------
23 // doi.org/10.1016/j.camwa.2015.03.017
24 inline void
25 gauss_quad20(
26 Vector <real> & aWeights,
27 Matrix <real> & aPoints )
28 {
29 aPoints.set_size( 2, 20 ) ;
30 aWeights.set_size( 20 ) ;
31
32 aPoints( 0, 0 ) = 0.4889268569743691 ;
33 aPoints( 1, 0 ) = 0.0000000000000000 ;
34
35 aPoints( 0, 1 ) = -0.4889268569743691 ;
36 aPoints( 1, 1 ) = 0.0000000000000000 ;
37
38 aPoints( 0, 2 ) = 0.0000000000000000 ;
39 aPoints( 1, 2 ) = 0.4889268569743691 ;
40
41 aPoints( 0, 3 ) = 0.0000000000000000 ;
42 aPoints( 1, 3 ) = -0.4889268569743691 ;
43
44 aPoints( 0, 4 ) = 0.9396552580968377 ;
45 aPoints( 1, 4 ) = 0.9396552580968377 ;
46
47 aPoints( 0, 5 ) = 0.6908805504863439 ;
48 aPoints( 1, 5 ) = 0.6908805504863439 ;
49
50 aPoints( 0, 6 ) = 0.9396552580968377 ;
51 aPoints( 1, 6 ) = -0.9396552580968377 ;
52
53 aPoints( 0, 7 ) = 0.6908805504863439 ;
54 aPoints( 1, 7 ) = -0.6908805504863439 ;
55
56 aPoints( 0, 8 ) = -0.9396552580968377 ;
57 aPoints( 1, 8 ) = 0.9396552580968377 ;
58
59 aPoints( 0, 9 ) = -0.6908805504863439 ;
60 aPoints( 1, 9 ) = 0.6908805504863439 ;
61
62 aPoints( 0, 10 ) = -0.9396552580968377 ;
63 aPoints( 1, 10 ) = -0.9396552580968377 ;
64
65 aPoints( 0, 11 ) = -0.6908805504863439 ;
66 aPoints( 1, 11 ) = -0.6908805504863439 ;
67
68 aPoints( 0, 12 ) = 0.9186204410567222 ;
69 aPoints( 1, 12 ) = 0.3448720253644036 ;
70
71 aPoints( 0, 13 ) = -0.9186204410567222 ;
72 aPoints( 1, 13 ) = 0.3448720253644036 ;
73
74 aPoints( 0, 14 ) = 0.9186204410567222 ;
75 aPoints( 1, 14 ) = -0.3448720253644036 ;
76
77 aPoints( 0, 15 ) = -0.9186204410567222 ;
78 aPoints( 1, 15 ) = -0.3448720253644036 ;
79
80 aPoints( 0, 16 ) = 0.3448720253644036 ;
81 aPoints( 1, 16 ) = 0.9186204410567222 ;
82
83 aPoints( 0, 17 ) = -0.3448720253644036 ;
84 aPoints( 1, 17 ) = 0.9186204410567222 ;
85
86 aPoints( 0, 18 ) = 0.3448720253644036 ;
87 aPoints( 1, 18 ) = -0.9186204410567222 ;
88
89 aPoints( 0, 19 ) = -0.3448720253644036 ;
90 aPoints( 1, 19 ) = -0.9186204410567222 ;
91
92 aWeights( 0 ) = 0.4541639606867490 ;
93 aWeights( 1 ) = 0.4541639606867490 ;
94 aWeights( 2 ) = 0.4541639606867490 ;
95 aWeights( 3 ) = 0.4541639606867490 ;
96 aWeights( 4 ) = 0.0427312318657758 ;
97 aWeights( 5 ) = 0.2142003609268616 ;
98 aWeights( 6 ) = 0.0427312318657758 ;
99 aWeights( 7 ) = 0.2142003609268616 ;
100 aWeights( 8 ) = 0.0427312318657758 ;
101 aWeights( 9 ) = 0.2142003609268616 ;
102 aWeights( 10 ) = 0.0427312318657758 ;
103 aWeights( 11 ) = 0.2142003609268616 ;
104 aWeights( 12 ) = 0.1444522232603068 ;
105 aWeights( 13 ) = 0.1444522232603068 ;
106 aWeights( 14 ) = 0.1444522232603068 ;
107 aWeights( 15 ) = 0.1444522232603068 ;
108 aWeights( 16 ) = 0.1444522232603068 ;
109 aWeights( 17 ) = 0.1444522232603068 ;
110 aWeights( 18 ) = 0.1444522232603068 ;
111 aWeights( 19 ) = 0.1444522232603068 ;
112 }
113
114// ----------------------------------------------------------------------------
115 } /* namespace integration */
116} /* namespace belfem */
117#endif //BELFEM_FN_INTPOINTS_GAUSS_QUAD20_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_quad20
void gauss_quad20(Vector< real > &aWeights, Matrix< real > &aPoints)
Definition fn_intpoints_gauss_quad20.hpp:25
belfem
USER GUIDES:
Definition cl_Capacitor.cpp:16