BELFEM 0.9.0
Berkeley Lab Finite Element Framework
Loading...
Searching...
No Matches
fn_intpoints_gauss_hex34.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_HEX34_HPP
13#define BELFEM_FN_INTPOINTS_GAUSS_HEX34_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
24 // source 10.2307/2004963
25 inline void
26 gauss_hex34(
27 Vector <real> & aWeights,
28 Matrix <real> & aPoints )
29 {
30 aPoints.set_size( 3, 34 ) ;
31 aWeights.set_size( 34 ) ;
32
33 aPoints( 0, 0 ) = 0.9317379999993866 ;
34 aPoints( 1, 0 ) = 0.0000000000000000 ;
35 aPoints( 2, 0 ) = 0.0000000000000000 ;
36
37 aPoints( 0, 1 ) = 0.0000000000000000 ;
38 aPoints( 1, 1 ) = 0.9317379999993866 ;
39 aPoints( 2, 1 ) = 0.0000000000000000 ;
40
41 aPoints( 0, 2 ) = 0.0000000000000000 ;
42 aPoints( 1, 2 ) = 0.0000000000000000 ;
43 aPoints( 2, 2 ) = 0.9317379999993866 ;
44
45 aPoints( 0, 3 ) = -0.9317379999993866 ;
46 aPoints( 1, 3 ) = 0.0000000000000000 ;
47 aPoints( 2, 3 ) = 0.0000000000000000 ;
48
49 aPoints( 0, 4 ) = 0.0000000000000000 ;
50 aPoints( 1, 4 ) = -0.9317379999993866 ;
51 aPoints( 2, 4 ) = 0.0000000000000000 ;
52
53 aPoints( 0, 5 ) = 0.0000000000000000 ;
54 aPoints( 1, 5 ) = 0.0000000000000000 ;
55 aPoints( 2, 5 ) = -0.9317379999993866 ;
56
57 aPoints( 0, 6 ) = 0.9167441779155981 ;
58 aPoints( 1, 6 ) = 0.9167441779155981 ;
59 aPoints( 2, 6 ) = 0.0000000000000000 ;
60
61 aPoints( 0, 7 ) = 0.9167441779155981 ;
62 aPoints( 1, 7 ) = 0.0000000000000000 ;
63 aPoints( 2, 7 ) = 0.9167441779155981 ;
64
65 aPoints( 0, 8 ) = 0.0000000000000000 ;
66 aPoints( 1, 8 ) = 0.9167441779155981 ;
67 aPoints( 2, 8 ) = 0.9167441779155981 ;
68
69 aPoints( 0, 9 ) = 0.9167441779155981 ;
70 aPoints( 1, 9 ) = -0.9167441779155981 ;
71 aPoints( 2, 9 ) = 0.0000000000000000 ;
72
73 aPoints( 0, 10 ) = 0.9167441779155981 ;
74 aPoints( 1, 10 ) = 0.0000000000000000 ;
75 aPoints( 2, 10 ) = -0.9167441779155981 ;
76
77 aPoints( 0, 11 ) = 0.0000000000000000 ;
78 aPoints( 1, 11 ) = 0.9167441779155981 ;
79 aPoints( 2, 11 ) = -0.9167441779155981 ;
80
81 aPoints( 0, 12 ) = -0.9167441779155981 ;
82 aPoints( 1, 12 ) = 0.9167441779155981 ;
83 aPoints( 2, 12 ) = 0.0000000000000000 ;
84
85 aPoints( 0, 13 ) = -0.9167441779155981 ;
86 aPoints( 1, 13 ) = 0.0000000000000000 ;
87 aPoints( 2, 13 ) = 0.9167441779155981 ;
88
89 aPoints( 0, 14 ) = 0.0000000000000000 ;
90 aPoints( 1, 14 ) = -0.9167441779155981 ;
91 aPoints( 2, 14 ) = 0.9167441779155981 ;
92
93 aPoints( 0, 15 ) = -0.9167441779155981 ;
94 aPoints( 1, 15 ) = -0.9167441779155981 ;
95 aPoints( 2, 15 ) = 0.0000000000000000 ;
96
97 aPoints( 0, 16 ) = -0.9167441779155981 ;
98 aPoints( 1, 16 ) = 0.0000000000000000 ;
99 aPoints( 2, 16 ) = -0.9167441779155981 ;
100
101 aPoints( 0, 17 ) = 0.0000000000000000 ;
102 aPoints( 1, 17 ) = -0.9167441779155981 ;
103 aPoints( 2, 17 ) = -0.9167441779155981 ;
104
105 aPoints( 0, 18 ) = 0.4086003799989615 ;
106 aPoints( 1, 18 ) = 0.4086003799989615 ;
107 aPoints( 2, 18 ) = 0.4086003799989615 ;
108
109 aPoints( 0, 19 ) = -0.4086003799989615 ;
110 aPoints( 1, 19 ) = 0.4086003799989615 ;
111 aPoints( 2, 19 ) = 0.4086003799989615 ;
112
113 aPoints( 0, 20 ) = 0.4086003799989615 ;
114 aPoints( 1, 20 ) = -0.4086003799989615 ;
115 aPoints( 2, 20 ) = 0.4086003799989615 ;
116
117 aPoints( 0, 21 ) = -0.4086003799989615 ;
118 aPoints( 1, 21 ) = -0.4086003799989615 ;
119 aPoints( 2, 21 ) = 0.4086003799989615 ;
120
121 aPoints( 0, 22 ) = 0.4086003799989615 ;
122 aPoints( 1, 22 ) = 0.4086003799989615 ;
123 aPoints( 2, 22 ) = -0.4086003799989615 ;
124
125 aPoints( 0, 23 ) = -0.4086003799989615 ;
126 aPoints( 1, 23 ) = 0.4086003799989615 ;
127 aPoints( 2, 23 ) = -0.4086003799989615 ;
128
129 aPoints( 0, 24 ) = 0.4086003799989615 ;
130 aPoints( 1, 24 ) = -0.4086003799989615 ;
131 aPoints( 2, 24 ) = -0.4086003799989615 ;
132
133 aPoints( 0, 25 ) = -0.4086003799989615 ;
134 aPoints( 1, 25 ) = -0.4086003799989615 ;
135 aPoints( 2, 25 ) = -0.4086003799989615 ;
136
137 aPoints( 0, 26 ) = 0.7398529499455402 ;
138 aPoints( 1, 26 ) = 0.7398529499455402 ;
139 aPoints( 2, 26 ) = 0.7398529499455402 ;
140
141 aPoints( 0, 27 ) = -0.7398529499455402 ;
142 aPoints( 1, 27 ) = 0.7398529499455402 ;
143 aPoints( 2, 27 ) = 0.7398529499455402 ;
144
145 aPoints( 0, 28 ) = 0.7398529499455402 ;
146 aPoints( 1, 28 ) = -0.7398529499455402 ;
147 aPoints( 2, 28 ) = 0.7398529499455402 ;
148
149 aPoints( 0, 29 ) = -0.7398529499455402 ;
150 aPoints( 1, 29 ) = -0.7398529499455402 ;
151 aPoints( 2, 29 ) = 0.7398529499455402 ;
152
153 aPoints( 0, 30 ) = 0.7398529499455402 ;
154 aPoints( 1, 30 ) = 0.7398529499455402 ;
155 aPoints( 2, 30 ) = -0.7398529499455402 ;
156
157 aPoints( 0, 31 ) = -0.7398529499455402 ;
158 aPoints( 1, 31 ) = 0.7398529499455402 ;
159 aPoints( 2, 31 ) = -0.7398529499455402 ;
160
161 aPoints( 0, 32 ) = 0.7398529499455402 ;
162 aPoints( 1, 32 ) = -0.7398529499455402 ;
163 aPoints( 2, 32 ) = -0.7398529499455402 ;
164
165 aPoints( 0, 33 ) = -0.7398529499455402 ;
166 aPoints( 1, 33 ) = -0.7398529499455402 ;
167 aPoints( 2, 33 ) = -0.7398529499455402 ;
168
169 aWeights( 0 ) = 0.2846544716914973 ;
170 aWeights( 1 ) = 0.2846544716914973 ;
171 aWeights( 2 ) = 0.2846544716914973 ;
172 aWeights( 3 ) = 0.2846544716914973 ;
173 aWeights( 4 ) = 0.2846544716914973 ;
174 aWeights( 5 ) = 0.2846544716914973 ;
175 aWeights( 6 ) = 0.0998314215798541 ;
176 aWeights( 7 ) = 0.0998314215798541 ;
177 aWeights( 8 ) = 0.0998314215798541 ;
178 aWeights( 9 ) = 0.0998314215798541 ;
179 aWeights( 10 ) = 0.0998314215798541 ;
180 aWeights( 11 ) = 0.0998314215798541 ;
181 aWeights( 12 ) = 0.0998314215798541 ;
182 aWeights( 13 ) = 0.0998314215798541 ;
183 aWeights( 14 ) = 0.0998314215798541 ;
184 aWeights( 15 ) = 0.0998314215798541 ;
185 aWeights( 16 ) = 0.0998314215798541 ;
186 aWeights( 17 ) = 0.0998314215798541 ;
187 aWeights( 18 ) = 0.4229418393139369 ;
188 aWeights( 19 ) = 0.4229418393139369 ;
189 aWeights( 20 ) = 0.4229418393139369 ;
190 aWeights( 21 ) = 0.4229418393139369 ;
191 aWeights( 22 ) = 0.4229418393139369 ;
192 aWeights( 23 ) = 0.4229418393139369 ;
193 aWeights( 24 ) = 0.4229418393139369 ;
194 aWeights( 25 ) = 0.4229418393139369 ;
195 aWeights( 26 ) = 0.2138201745476591 ;
196 aWeights( 27 ) = 0.2138201745476591 ;
197 aWeights( 28 ) = 0.2138201745476591 ;
198 aWeights( 29 ) = 0.2138201745476591 ;
199 aWeights( 30 ) = 0.2138201745476591 ;
200 aWeights( 31 ) = 0.2138201745476591 ;
201 aWeights( 32 ) = 0.2138201745476591 ;
202 aWeights( 33 ) = 0.2138201745476591 ;
203
204 }
205
206// ----------------------------------------------------------------------------
207 } /* namespace integration */
208} /* namespace belfem */
209#endif //BELFEM_FN_INTPOINTS_GAUSS_HEX34_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_hex34
void gauss_hex34(Vector< real > &aWeights, Matrix< real > &aPoints)
Definition fn_intpoints_gauss_hex34.hpp:26
belfem
USER GUIDES:
Definition cl_Capacitor.cpp:16