BELFEM 0.9.0
Berkeley Lab Finite Element Framework
Loading...
Searching...
No Matches
fn_intpoints_gauss_quad37.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_QUAD37_HPP
13#define BELFEM_FN_INTPOINTS_GAUSS_QUAD37_HPP
14
15
16#include "typedefs.hpp"
17#include "cl_Vector.hpp"
18#include "cl_Matrix.hpp"
19
20namespace belfem
21{
22 namespace integration
23 {
24//----------------------------------------------------------------------------
25
26 // doi.org/10.1016/j.camwa.2015.03.017
27 inline void
28 gauss_quad37(
29 Vector <real> & aWeights,
30 Matrix <real> & aPoints )
31 {
32 aPoints.set_size( 2, 37 ) ;
33 aWeights.set_size( 37 ) ;
34
35 aPoints( 0, 0 ) = 0.0000000000000000 ;
36 aPoints( 1, 0 ) = 0.0000000000000000 ;
37
38 aPoints( 0, 1 ) = 0.9834613326132455 ;
39 aPoints( 1, 1 ) = 0.0000000000000000 ;
40
41 aPoints( 0, 2 ) = 0.6398614183671097 ;
42 aPoints( 1, 2 ) = 0.0000000000000000 ;
43
44 aPoints( 0, 3 ) = -0.9834613326132455 ;
45 aPoints( 1, 3 ) = 0.0000000000000000 ;
46
47 aPoints( 0, 4 ) = -0.6398614183671097 ;
48 aPoints( 1, 4 ) = 0.0000000000000000 ;
49
50 aPoints( 0, 5 ) = 0.0000000000000000 ;
51 aPoints( 1, 5 ) = 0.9834613326132455 ;
52
53 aPoints( 0, 6 ) = 0.0000000000000000 ;
54 aPoints( 1, 6 ) = 0.6398614183671097 ;
55
56 aPoints( 0, 7 ) = 0.0000000000000000 ;
57 aPoints( 1, 7 ) = -0.9834613326132455 ;
58
59 aPoints( 0, 8 ) = 0.0000000000000000 ;
60 aPoints( 1, 8 ) = -0.6398614183671097 ;
61
62 aPoints( 0, 9 ) = 0.9187784880797114 ;
63 aPoints( 1, 9 ) = 0.9187784880797114 ;
64
65 aPoints( 0, 10 ) = 0.3796285051867438 ;
66 aPoints( 1, 10 ) = 0.3796285051867438 ;
67
68 aPoints( 0, 11 ) = 0.6995542165133511 ;
69 aPoints( 1, 11 ) = 0.6995542165133511 ;
70
71 aPoints( 0, 12 ) = 0.9187784880797114 ;
72 aPoints( 1, 12 ) = -0.9187784880797114 ;
73
74 aPoints( 0, 13 ) = 0.3796285051867438 ;
75 aPoints( 1, 13 ) = -0.3796285051867438 ;
76
77 aPoints( 0, 14 ) = 0.6995542165133511 ;
78 aPoints( 1, 14 ) = -0.6995542165133511 ;
79
80 aPoints( 0, 15 ) = -0.9187784880797114 ;
81 aPoints( 1, 15 ) = 0.9187784880797114 ;
82
83 aPoints( 0, 16 ) = -0.3796285051867438 ;
84 aPoints( 1, 16 ) = 0.3796285051867438 ;
85
86 aPoints( 0, 17 ) = -0.6995542165133511 ;
87 aPoints( 1, 17 ) = 0.6995542165133511 ;
88
89 aPoints( 0, 18 ) = -0.9187784880797114 ;
90 aPoints( 1, 18 ) = -0.9187784880797114 ;
91
92 aPoints( 0, 19 ) = -0.3796285051867438 ;
93 aPoints( 1, 19 ) = -0.3796285051867438 ;
94
95 aPoints( 0, 20 ) = -0.6995542165133511 ;
96 aPoints( 1, 20 ) = -0.6995542165133511 ;
97
98 aPoints( 0, 21 ) = 0.6435973749966181 ;
99 aPoints( 1, 21 ) = 0.9750688839059836 ;
100
101 aPoints( 0, 22 ) = 0.3335398811647831 ;
102 aPoints( 1, 22 ) = 0.8644276092670619 ;
103
104 aPoints( 0, 23 ) = -0.6435973749966181 ;
105 aPoints( 1, 23 ) = 0.9750688839059836 ;
106
107 aPoints( 0, 24 ) = -0.3335398811647831 ;
108 aPoints( 1, 24 ) = 0.8644276092670619 ;
109
110 aPoints( 0, 25 ) = 0.6435973749966181 ;
111 aPoints( 1, 25 ) = -0.9750688839059836 ;
112
113 aPoints( 0, 26 ) = 0.3335398811647831 ;
114 aPoints( 1, 26 ) = -0.8644276092670619 ;
115
116 aPoints( 0, 27 ) = -0.6435973749966181 ;
117 aPoints( 1, 27 ) = -0.9750688839059836 ;
118
119 aPoints( 0, 28 ) = -0.3335398811647831 ;
120 aPoints( 1, 28 ) = -0.8644276092670619 ;
121
122 aPoints( 0, 29 ) = 0.9750688839059836 ;
123 aPoints( 1, 29 ) = 0.6435973749966181 ;
124
125 aPoints( 0, 30 ) = 0.8644276092670619 ;
126 aPoints( 1, 30 ) = 0.3335398811647831 ;
127
128 aPoints( 0, 31 ) = -0.9750688839059836 ;
129 aPoints( 1, 31 ) = 0.6435973749966181 ;
130
131 aPoints( 0, 32 ) = -0.8644276092670619 ;
132 aPoints( 1, 32 ) = 0.3335398811647831 ;
133
134 aPoints( 0, 33 ) = 0.9750688839059836 ;
135 aPoints( 1, 33 ) = -0.6435973749966181 ;
136
137 aPoints( 0, 34 ) = 0.8644276092670619 ;
138 aPoints( 1, 34 ) = -0.3335398811647831 ;
139
140 aPoints( 0, 35 ) = -0.9750688839059836 ;
141 aPoints( 1, 35 ) = -0.6435973749966181 ;
142
143 aPoints( 0, 36 ) = -0.8644276092670619 ;
144 aPoints( 1, 36 ) = -0.3335398811647831 ;
145
146 aWeights( 0 ) = 0.2999933443589096 ;
147 aWeights( 1 ) = 0.0381286253498527 ;
148 aWeights( 2 ) = 0.1845354689698072 ;
149 aWeights( 3 ) = 0.0381286253498527 ;
150 aWeights( 4 ) = 0.1845354689698072 ;
151 aWeights( 5 ) = 0.0381286253498527 ;
152 aWeights( 6 ) = 0.1845354689698072 ;
153 aWeights( 7 ) = 0.0381286253498527 ;
154 aWeights( 8 ) = 0.1845354689698072 ;
155 aWeights( 9 ) = 0.0395071406474524 ;
156 aWeights( 10 ) = 0.2313985194006854 ;
157 aWeights( 11 ) = 0.1372420967130035 ;
158 aWeights( 12 ) = 0.0395071406474524 ;
159 aWeights( 13 ) = 0.2313985194006854 ;
160 aWeights( 14 ) = 0.1372420967130035 ;
161 aWeights( 15 ) = 0.0395071406474524 ;
162 aWeights( 16 ) = 0.2313985194006854 ;
163 aWeights( 17 ) = 0.1372420967130035 ;
164 aWeights( 18 ) = 0.0395071406474524 ;
165 aWeights( 19 ) = 0.2313985194006854 ;
166 aWeights( 20 ) = 0.1372420967130035 ;
167 aWeights( 21 ) = 0.0335197800503814 ;
168 aWeights( 22 ) = 0.1135751263643542 ;
169 aWeights( 23 ) = 0.0335197800503814 ;
170 aWeights( 24 ) = 0.1135751263643542 ;
171 aWeights( 25 ) = 0.0335197800503814 ;
172 aWeights( 26 ) = 0.1135751263643542 ;
173 aWeights( 27 ) = 0.0335197800503814 ;
174 aWeights( 28 ) = 0.1135751263643542 ;
175 aWeights( 29 ) = 0.0335197800503814 ;
176 aWeights( 30 ) = 0.1135751263643542 ;
177 aWeights( 31 ) = 0.0335197800503814 ;
178 aWeights( 32 ) = 0.1135751263643542 ;
179 aWeights( 33 ) = 0.0335197800503814 ;
180 aWeights( 34 ) = 0.1135751263643542 ;
181 aWeights( 35 ) = 0.0335197800503814 ;
182 aWeights( 36 ) = 0.1135751263643542 ;
183
184 }
185
186//----------------------------------------------------------------------------
187 } /* namespace integration */
188} /* namespace belfem */
189
190#endif //BELFEM_FN_INTPOINTS_GAUSS_QUAD37_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_quad37
void gauss_quad37(Vector< real > &aWeights, Matrix< real > &aPoints)
Definition fn_intpoints_gauss_quad37.hpp:28
belfem
USER GUIDES:
Definition cl_Capacitor.cpp:16