BELFEM 0.9.3
Berkeley Lab Finite Element Framework
Loading...
Searching...
No Matches
fn_intpoints_gauss_pyra8.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_PYRA8_HPP
13#define BELFEM_FN_INTPOINTS_GAUSS_PYRA8_HPP
14
15#include "typedefs.hpp"
16#include "cl_Vector.hpp"
17#include "cl_Matrix.hpp"
18
19
20namespace belfem
21{
22 namespace integration
23 {
24// ----------------------------------------------------------------------------
25
34 inline void
35 gauss_pyra8(
36 Vector <real> & aWeights,
37 Matrix <real> & aPoints )
38 {
39 aPoints.set_size( 3, 8 );
40 aWeights.set_size( 8 );
41
42 aPoints( 0, 0 ) = -0.50661630334978724 ;
43 aPoints( 1, 0 ) = -0.50661630334978724 ;
44 aPoints( 2, 0 ) = 0.1225148226554415 ;
45 aWeights( 0 ) = 0.23254745125350801 ;
46
47 aPoints( 0, 1 ) = 0.50661630334978724 ;
48 aPoints( 1, 1 ) = -0.50661630334978724 ;
49 aPoints( 2, 1 ) = 0.1225148226554415 ;
50 aWeights( 1 ) = 0.23254745125350801 ;
51
52 aPoints( 0, 2 ) = -0.50661630334978724 ;
53 aPoints( 1, 2 ) = 0.50661630334978724 ;
54 aPoints( 2, 2 ) = 0.1225148226554415 ;
55 aWeights( 2 ) = 0.23254745125350801 ;
56
57 aPoints( 0, 3 ) = 0.50661630334978724 ;
58 aPoints( 1, 3 ) = 0.50661630334978724 ;
59 aPoints( 2, 3 ) = 0.1225148226554415 ;
60 aWeights( 3 ) = 0.23254745125350801 ;
61
62 aPoints( 0, 4 ) = -0.26318405556971358 ;
63 aPoints( 1, 4 ) = -0.26318405556971358 ;
64 aPoints( 2, 4 ) = 0.54415184401122529 ;
65 aWeights( 4 ) = 0.10078588207982532 ;
66
67 aPoints( 0, 5 ) = 0.26318405556971358 ;
68 aPoints( 1, 5 ) = -0.26318405556971358 ;
69 aPoints( 2, 5 ) = 0.54415184401122529 ;
70 aWeights( 5 ) = 0.10078588207982532 ;
71
72 aPoints( 0, 6 ) = -0.26318405556971358 ;
73 aPoints( 1, 6 ) = 0.26318405556971358 ;
74 aPoints( 2, 6 ) = 0.54415184401122529 ;
75 aWeights( 6 ) = 0.10078588207982532 ;
76
77 aPoints( 0, 7 ) = 0.26318405556971358 ;
78 aPoints( 1, 7 ) = 0.26318405556971358 ;
79 aPoints( 2, 7 ) = 0.54415184401122529 ;
80 aWeights( 7 ) = 0.10078588207982532 ;
81 }
82
83// ----------------------------------------------------------------------------
84 } /* namespace integration */
85} /* namespace belfem */
86
87#endif //BELFEM_FN_INTPOINTS_GAUSS_PYRA8_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_pyra8
void gauss_pyra8(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_pyra8.hpp:35
belfem
USER GUIDES:
Definition cl_Capacitor.cpp:16