BELFEM 0.9.0
Berkeley Lab Finite Element Framework
Loading...
Searching...
No Matches
fn_intpoints_gauss_quad28.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_QUAD28_HPP
13#define BELFEM_FN_INTPOINTS_GAUSS_QUAD28_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 // doi.org/10.1016/j.camwa.2015.03.017
25 inline void
26 gauss_quad28(
27 Vector <real> & aWeights,
28 Matrix <real> & aPoints )
29 {
30 aPoints.set_size( 2, 28 ) ;
31 aWeights.set_size( 28 ) ;
32
33 aPoints( 0, 0 ) = 0.7146178296646060 ;
34 aPoints( 1, 0 ) = 0.0000000000000000 ;
35
36 aPoints( 0, 1 ) = -0.7146178296646060 ;
37 aPoints( 1, 1 ) = 0.0000000000000000 ;
38
39 aPoints( 0, 2 ) = 0.0000000000000000 ;
40 aPoints( 1, 2 ) = 0.7146178296646060 ;
41
42 aPoints( 0, 3 ) = 0.0000000000000000 ;
43 aPoints( 1, 3 ) = -0.7146178296646060 ;
44
45 aPoints( 0, 4 ) = 0.2736572101714596 ;
46 aPoints( 1, 4 ) = 0.2736572101714596 ;
47
48 aPoints( 0, 5 ) = 0.6366039322123010 ;
49 aPoints( 1, 5 ) = 0.6366039322123010 ;
50
51 aPoints( 0, 6 ) = 0.2736572101714596 ;
52 aPoints( 1, 6 ) = -0.2736572101714596 ;
53
54 aPoints( 0, 7 ) = 0.6366039322123010 ;
55 aPoints( 1, 7 ) = -0.6366039322123010 ;
56
57 aPoints( 0, 8 ) = -0.2736572101714596 ;
58 aPoints( 1, 8 ) = 0.2736572101714596 ;
59
60 aPoints( 0, 9 ) = -0.6366039322123010 ;
61 aPoints( 1, 9 ) = 0.6366039322123010 ;
62
63 aPoints( 0, 10 ) = -0.2736572101714596 ;
64 aPoints( 1, 10 ) = -0.2736572101714596 ;
65
66 aPoints( 0, 11 ) = -0.6366039322123010 ;
67 aPoints( 1, 11 ) = -0.6366039322123010 ;
68
69 aPoints( 0, 12 ) = 0.9516303887840335 ;
70 aPoints( 1, 12 ) = 0.8155654336896384 ;
71
72 aPoints( 0, 13 ) = 0.3462072000476454 ;
73 aPoints( 1, 13 ) = 0.9355678714875911 ;
74
75 aPoints( 0, 14 ) = -0.9516303887840335 ;
76 aPoints( 1, 14 ) = 0.8155654336896384 ;
77
78 aPoints( 0, 15 ) = -0.3462072000476454 ;
79 aPoints( 1, 15 ) = 0.9355678714875911 ;
80
81 aPoints( 0, 16 ) = 0.9516303887840335 ;
82 aPoints( 1, 16 ) = -0.8155654336896384 ;
83
84 aPoints( 0, 17 ) = 0.3462072000476454 ;
85 aPoints( 1, 17 ) = -0.9355678714875911 ;
86
87 aPoints( 0, 18 ) = -0.9516303887840335 ;
88 aPoints( 1, 18 ) = -0.8155654336896384 ;
89
90 aPoints( 0, 19 ) = -0.3462072000476454 ;
91 aPoints( 1, 19 ) = -0.9355678714875911 ;
92
93 aPoints( 0, 20 ) = 0.8155654336896384 ;
94 aPoints( 1, 20 ) = 0.9516303887840335 ;
95
96 aPoints( 0, 21 ) = 0.9355678714875911 ;
97 aPoints( 1, 21 ) = 0.3462072000476454 ;
98
99 aPoints( 0, 22 ) = -0.8155654336896384 ;
100 aPoints( 1, 22 ) = 0.9516303887840335 ;
101
102 aPoints( 0, 23 ) = -0.9355678714875911 ;
103 aPoints( 1, 23 ) = 0.3462072000476454 ;
104
105 aPoints( 0, 24 ) = 0.8155654336896384 ;
106 aPoints( 1, 24 ) = -0.9516303887840335 ;
107
108 aPoints( 0, 25 ) = 0.9355678714875911 ;
109 aPoints( 1, 25 ) = -0.3462072000476454 ;
110
111 aPoints( 0, 26 ) = -0.8155654336896384 ;
112 aPoints( 1, 26 ) = -0.9516303887840335 ;
113
114 aPoints( 0, 27 ) = -0.9355678714875911 ;
115 aPoints( 1, 27 ) = -0.3462072000476454 ;
116
117 aWeights( 0 ) = 0.2174004398687120 ;
118 aWeights( 1 ) = 0.2174004398687120 ;
119 aWeights( 2 ) = 0.2174004398687120 ;
120 aWeights( 3 ) = 0.2174004398687120 ;
121 aWeights( 4 ) = 0.2772741029838511 ;
122 aWeights( 5 ) = 0.2139336378782481 ;
123 aWeights( 6 ) = 0.2772741029838511 ;
124 aWeights( 7 ) = 0.2139336378782481 ;
125 aWeights( 8 ) = 0.2772741029838511 ;
126 aWeights( 9 ) = 0.2139336378782481 ;
127 aWeights( 10 ) = 0.2772741029838511 ;
128 aWeights( 11 ) = 0.2139336378782481 ;
129 aWeights( 12 ) = 0.0440745691149831 ;
130 aWeights( 13 ) = 0.1016213405196113 ;
131 aWeights( 14 ) = 0.0440745691149831 ;
132 aWeights( 15 ) = 0.1016213405196113 ;
133 aWeights( 16 ) = 0.0440745691149831 ;
134 aWeights( 17 ) = 0.1016213405196113 ;
135 aWeights( 18 ) = 0.0440745691149831 ;
136 aWeights( 19 ) = 0.1016213405196113 ;
137 aWeights( 20 ) = 0.0440745691149831 ;
138 aWeights( 21 ) = 0.1016213405196113 ;
139 aWeights( 22 ) = 0.0440745691149831 ;
140 aWeights( 23 ) = 0.1016213405196113 ;
141 aWeights( 24 ) = 0.0440745691149831 ;
142 aWeights( 25 ) = 0.1016213405196113 ;
143 aWeights( 26 ) = 0.0440745691149831 ;
144 aWeights( 27 ) = 0.1016213405196113 ;
145 }
146
147//----------------------------------------------------------------------------
148 } /* namespace integration */
149} /* namespace belfem */
150
151#endif //BELFEM_FN_INTPOINTS_GAUSS_QUAD28_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_quad28
void gauss_quad28(Vector< real > &aWeights, Matrix< real > &aPoints)
Definition fn_intpoints_gauss_quad28.hpp:26
belfem
USER GUIDES:
Definition cl_Capacitor.cpp:16