BELFEM 0.9.3
Berkeley Lab Finite Element Framework
Loading...
Searching...
No Matches
fn_intpoints_gauss_pyra27.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_PYRA27_HPP
13#define BELFEM_FN_INTPOINTS_GAUSS_PYRA27_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
33 inline void
34 gauss_pyra27(
35 Vector <real> & aWeights,
36 Matrix <real> & aPoints )
37 {
38 aPoints.set_size( 3, 27 );
39 aWeights.set_size( 27 );
40
41 aPoints( 0, 0 ) = -0.71805574131988903 ;
42 aPoints( 1, 0 ) = -0.71805574131988903 ;
43 aPoints( 2, 0 ) = 0.072994024073149699 ;
44 aWeights( 0 ) = 0.04849887687187858 ;
45
46 aPoints( 0, 1 ) = 0 ;
47 aPoints( 1, 1 ) = -0.71805574131988903 ;
48 aPoints( 2, 1 ) = 0.072994024073149699 ;
49 aWeights( 1 ) = 0.077598202995005638 ;
50
51 aPoints( 0, 2 ) = 0.71805574131988903 ;
52 aPoints( 1, 2 ) = -0.71805574131988903 ;
53 aPoints( 2, 2 ) = 0.072994024073149699 ;
54 aWeights( 2 ) = 0.04849887687187858 ;
55
56 aPoints( 0, 3 ) = -0.71805574131988903 ;
57 aPoints( 1, 3 ) = 0 ;
58 aPoints( 2, 3 ) = 0.072994024073149699 ;
59 aWeights( 3 ) = 0.077598202995005638 ;
60
61 aPoints( 0, 4 ) = 0 ;
62 aPoints( 1, 4 ) = 0 ;
63 aPoints( 2, 4 ) = 0.072994024073149699 ;
64 aWeights( 4 ) = 0.1241571247920089 ;
65
66 aPoints( 0, 5 ) = 0.71805574131988903 ;
67 aPoints( 1, 5 ) = 0 ;
68 aPoints( 2, 5 ) = 0.072994024073149699 ;
69 aWeights( 5 ) = 0.077598202995005638 ;
70
71 aPoints( 0, 6 ) = -0.71805574131988903 ;
72 aPoints( 1, 6 ) = 0.71805574131988903 ;
73 aPoints( 2, 6 ) = 0.072994024073149699 ;
74 aWeights( 6 ) = 0.04849887687187858 ;
75
76 aPoints( 0, 7 ) = 0 ;
77 aPoints( 1, 7 ) = 0.71805574131988903 ;
78 aPoints( 2, 7 ) = 0.072994024073149699 ;
79 aWeights( 7 ) = 0.077598202995005638 ;
80
81 aPoints( 0, 8 ) = 0.71805574131988903 ;
82 aPoints( 1, 8 ) = 0.71805574131988903 ;
83 aPoints( 2, 8 ) = 0.072994024073149699 ;
84 aWeights( 8 ) = 0.04849887687187858 ;
85
86 aPoints( 0, 9 ) = -0.50580870785392507 ;
87 aPoints( 1, 9 ) = -0.50580870785392507 ;
88 aPoints( 2, 9 ) = 0.34700376603835181 ;
89 aWeights( 9 ) = 0.045137737425884644 ;
90
91 aPoints( 0, 10 ) = 0 ;
92 aPoints( 1, 10 ) = -0.50580870785392507 ;
93 aPoints( 2, 10 ) = 0.34700376603835181 ;
94 aWeights( 10 ) = 0.072220379881415345 ;
95
96 aPoints( 0, 11 ) = 0.50580870785392507 ;
97 aPoints( 1, 11 ) = -0.50580870785392507 ;
98 aPoints( 2, 11 ) = 0.34700376603835181 ;
99 aWeights( 11 ) = 0.045137737425884644 ;
100
101 aPoints( 0, 12 ) = -0.50580870785392507 ;
102 aPoints( 1, 12 ) = 0 ;
103 aPoints( 2, 12 ) = 0.34700376603835181 ;
104 aWeights( 12 ) = 0.072220379881415345 ;
105
106 aPoints( 0, 13 ) = 0 ;
107 aPoints( 1, 13 ) = 0 ;
108 aPoints( 2, 13 ) = 0.34700376603835181 ;
109 aWeights( 13 ) = 0.11555260781026443 ;
110
111 aPoints( 0, 14 ) = 0.50580870785392507 ;
112 aPoints( 1, 14 ) = 0 ;
113 aPoints( 2, 14 ) = 0.34700376603835181 ;
114 aWeights( 14 ) = 0.072220379881415345 ;
115
116 aPoints( 0, 15 ) = -0.50580870785392507 ;
117 aPoints( 1, 15 ) = 0.50580870785392507 ;
118 aPoints( 2, 15 ) = 0.34700376603835181 ;
119 aWeights( 15 ) = 0.045137737425884644 ;
120
121 aPoints( 0, 16 ) = 0 ;
122 aPoints( 1, 16 ) = 0.50580870785392507 ;
123 aPoints( 2, 16 ) = 0.34700376603835181 ;
124 aWeights( 16 ) = 0.072220379881415345 ;
125
126 aPoints( 0, 17 ) = 0.50580870785392507 ;
127 aPoints( 1, 17 ) = 0.50580870785392507 ;
128 aPoints( 2, 17 ) = 0.34700376603835181 ;
129 aWeights( 17 ) = 0.045137737425884644 ;
130
131 aPoints( 0, 18 ) = -0.22850430565396737 ;
132 aPoints( 1, 18 ) = -0.22850430565396737 ;
133 aPoints( 2, 18 ) = 0.70500220988849838 ;
134 aWeights( 18 ) = 0.0092440441384508461 ;
135
136 aPoints( 0, 19 ) = 0 ;
137 aPoints( 1, 19 ) = -0.22850430565396737 ;
138 aPoints( 2, 19 ) = 0.70500220988849838 ;
139 aWeights( 19 ) = 0.014790470621521336 ;
140
141 aPoints( 0, 20 ) = 0.22850430565396737 ;
142 aPoints( 1, 20 ) = -0.22850430565396737 ;
143 aPoints( 2, 20 ) = 0.70500220988849838 ;
144 aWeights( 20 ) = 0.0092440441384508461 ;
145
146 aPoints( 0, 21 ) = -0.22850430565396737 ;
147 aPoints( 1, 21 ) = 0 ;
148 aPoints( 2, 21 ) = 0.70500220988849838 ;
149 aWeights( 21 ) = 0.014790470621521336 ;
150
151 aPoints( 0, 22 ) = 0 ;
152 aPoints( 1, 22 ) = 0 ;
153 aPoints( 2, 22 ) = 0.70500220988849838 ;
154 aWeights( 22 ) = 0.023664752994434112 ;
155
156 aPoints( 0, 23 ) = 0.22850430565396737 ;
157 aPoints( 1, 23 ) = 0 ;
158 aPoints( 2, 23 ) = 0.70500220988849838 ;
159 aWeights( 23 ) = 0.014790470621521336 ;
160
161 aPoints( 0, 24 ) = -0.22850430565396737 ;
162 aPoints( 1, 24 ) = 0.22850430565396737 ;
163 aPoints( 2, 24 ) = 0.70500220988849838 ;
164 aWeights( 24 ) = 0.0092440441384508461 ;
165
166 aPoints( 0, 25 ) = 0 ;
167 aPoints( 1, 25 ) = 0.22850430565396737 ;
168 aPoints( 2, 25 ) = 0.70500220988849838 ;
169 aWeights( 25 ) = 0.014790470621521336 ;
170
171 aPoints( 0, 26 ) = 0.22850430565396737 ;
172 aPoints( 1, 26 ) = 0.22850430565396737 ;
173 aPoints( 2, 26 ) = 0.70500220988849838 ;
174 aWeights( 26 ) = 0.0092440441384508461 ;
175 }
176
177// ----------------------------------------------------------------------------
178 } /* namespace integration */
179} /* namespace belfem */
180
181#endif //BELFEM_FN_INTPOINTS_GAUSS_PYRA27_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:164
belfem::Vector::set_size
void set_size(const size_t aNumRows)
change the size of the vector
Definition cl_AR_Vector.hpp:215
belfem::integration
Definition fn_intpoints.cpp:576
belfem::integration::gauss_pyra27
void gauss_pyra27(Vector< real > &aWeights, Matrix< real > &aPoints)
Conical-product Gauss rule on the reference pyramid ( base [-1,1]^2 at zeta = 0, apex ( 0,...
Definition fn_intpoints_gauss_pyra27.hpp:34
belfem
USER GUIDES:
Definition cl_Capacitor.cpp:16