libMesh
cell_inf_hex16.h
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2024 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
3 
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License, or (at your option) any later version.
8 
9 // This library is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // Lesser General Public License for more details.
13 
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 
18 
19 
20 #ifndef LIBMESH_CELL_INF_HEX16_H
21 #define LIBMESH_CELL_INF_HEX16_H
22 
23 #include "libmesh/libmesh_config.h"
24 
25 #ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
26 
27 // Local includes
28 #include "libmesh/cell_inf_hex.h"
29 
30 namespace libMesh
31 {
32 
70 class InfHex16 final : public InfHex
71 {
72 public:
73 
77  explicit
78  InfHex16 (Elem * p=nullptr) :
80  {}
81 
82  InfHex16 (InfHex16 &&) = delete;
83  InfHex16 (const InfHex16 &) = delete;
84  InfHex16 & operator= (const InfHex16 &) = delete;
85  InfHex16 & operator= (InfHex16 &&) = delete;
86  virtual ~InfHex16() = default;
87 
91  virtual unsigned int n_nodes() const override { return num_nodes; }
92 
96  virtual ElemType type () const override { return INFHEX16; }
97 
101  virtual unsigned int n_sub_elem() const override { return 1; }
102 
107  virtual bool is_node_on_side(const unsigned int n,
108  const unsigned int s) const override;
109 
110  virtual std::vector<unsigned int> nodes_on_side(const unsigned int s) const override;
111 
112  virtual std::vector<unsigned int> nodes_on_edge(const unsigned int e) const override;
113 
118  virtual bool is_node_on_edge(const unsigned int n,
119  const unsigned int e) const override;
120 
124  virtual Order default_order() const override;
125 
129  virtual unsigned int local_side_node(unsigned int side,
130  unsigned int side_node) const override;
131 
135  virtual unsigned int local_edge_node(unsigned int edge,
136  unsigned int edge_node) const override;
137 
144  virtual std::unique_ptr<Elem> build_side_ptr (const unsigned int i,
145  bool proxy=false) override;
146 
151  virtual void build_side_ptr (std::unique_ptr<Elem> & elem,
152  const unsigned int i) override;
153 
160  virtual std::unique_ptr<Elem> build_edge_ptr (const unsigned int i) override;
161 
166  virtual void build_edge_ptr (std::unique_ptr<Elem> & edge, const unsigned int i) override;
167 
168  virtual void connectivity(const unsigned int sc,
169  const IOPackage iop,
170  std::vector<dof_id_type> & conn) const override;
171 
172  unsigned int vtk_element_type (const unsigned int) const
173  { return 12; }
174 
178  virtual unsigned int n_second_order_adjacent_vertices (const unsigned int) const override
179  { return 2; }
180 
187  virtual unsigned short int second_order_adjacent_vertex (const unsigned int n,
188  const unsigned int v) const override;
189 
195  virtual std::pair<unsigned short int, unsigned short int>
196  second_order_child_vertex (const unsigned int n) const override;
197 
201  static const int num_nodes = 16;
202  static const int num_sides = 5;
203  static const int num_edges = 8;
204  static const int num_children = 4;
205  static const int nodes_per_side = 8;
206  static const int nodes_per_edge = 3;
207 
212  static const unsigned int side_nodes_map[num_sides][nodes_per_side];
213 
218  static const unsigned int edge_nodes_map[num_edges][nodes_per_edge];
219 
220  virtual void permute(unsigned int perm_num) override final;
221 
222  virtual void flip(BoundaryInfo *) override final;
223 
224  ElemType side_type (const unsigned int s) const override final;
225 
226 protected:
227 
232 
233 
234 #ifdef LIBMESH_ENABLE_AMR
235 
239  virtual Real embedding_matrix (const unsigned int i,
240  const unsigned int j,
241  const unsigned int k) const override
242  { return _embedding_matrix[i][j][k]; }
243 
249 
251 
252 #endif // LIBMESH_ENABLE_AMR
253 
254 };
255 
256 } // namespace libMesh
257 
258 #endif // ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
259 
260 
261 #endif // LIBMESH_CELL_INF_HEX16_H
ElemType
Defines an enum for geometric element types.
Order
defines an enum for polynomial orders.
Definition: enum_order.h:40
A Node is like a Point, but with more information.
Definition: node.h:52
static const int nodes_per_edge
virtual std::pair< unsigned short int, unsigned short int > second_order_child_vertex(const unsigned int n) const override
IOPackage
libMesh interfaces with several different software packages for the purposes of creating, reading, and writing mesh files.
virtual std::vector< unsigned int > nodes_on_side(const unsigned int s) const override
This is the base class from which all geometric element types are derived.
Definition: elem.h:94
unsigned int vtk_element_type(const unsigned int) const
static const int num_edges
static const int num_nodes
Geometric constants for InfHex16.
The libMesh namespace provides an interface to certain functionality in the library.
static const unsigned int edge_nodes_map[num_edges][nodes_per_edge]
This maps the node of the edge to element node numbers.
Node * _nodelinks_data[num_nodes]
Data for links to nodes.
virtual unsigned int local_edge_node(unsigned int edge, unsigned int edge_node) const override
virtual Order default_order() const override
virtual void flip(BoundaryInfo *) override final
Flips the element (by swapping node and neighbor pointers) to have a mapping Jacobian of opposite sig...
virtual std::unique_ptr< Elem > build_edge_ptr(const unsigned int i) override
static const unsigned int side_nodes_map[num_sides][nodes_per_side]
This maps the node of the side to element node numbers.
static const Real _embedding_matrix[num_children][num_nodes][num_nodes]
Matrix that computes new nodal locations/solution values from current nodes/solution.
ElemType side_type(const unsigned int s) const override final
static const int num_children
The InfHex16 is an infinite element in 3D composed of 16 nodes.
The InfHex is an element in 3D with 5 sides.
Definition: cell_inf_hex.h:46
InfHex16 & operator=(const InfHex16 &)=delete
virtual void permute(unsigned int perm_num) override final
Permutes the element (by swapping node and neighbor pointers) according to the specified index...
The BoundaryInfo class contains information relevant to boundary conditions including storing faces...
Definition: boundary_info.h:57
virtual std::vector< unsigned int > nodes_on_edge(const unsigned int e) const override
virtual unsigned int n_sub_elem() const override
virtual unsigned int local_side_node(unsigned int side, unsigned int side_node) const override
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
virtual bool is_node_on_side(const unsigned int n, const unsigned int s) const override
virtual ElemType type() const override
virtual void connectivity(const unsigned int sc, const IOPackage iop, std::vector< dof_id_type > &conn) const override
virtual unsigned int n_nodes() const override
static const int nodes_per_side
virtual unsigned int n_second_order_adjacent_vertices(const unsigned int) const override
virtual std::unique_ptr< Elem > build_side_ptr(const unsigned int i, bool proxy=false) override
static const int num_sides
InfHex16(Elem *p=nullptr)
Constructor.
virtual Real embedding_matrix(const unsigned int i, const unsigned int j, const unsigned int k) const override
Matrix used to create the elements children.
virtual unsigned short int second_order_adjacent_vertex(const unsigned int n, const unsigned int v) const override
virtual bool is_node_on_edge(const unsigned int n, const unsigned int e) const override
virtual ~InfHex16()=default