libMesh
face_inf_quad4.h
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2025 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 #ifndef LIBMESH_FACE_INF_QUAD4_H
20 #define LIBMESH_FACE_INF_QUAD4_H
21 
22 
23 #include "libmesh/libmesh_config.h"
24 
25 #ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
26 
27 // Local includes
28 #include "libmesh/face_inf_quad.h"
29 
30 namespace libMesh
31 {
32 
52 class InfQuad4 : public InfQuad
53 {
54 public:
55 
59  explicit
60  InfQuad4 (Elem * p=nullptr) :
62 
63  InfQuad4 (InfQuad4 &&) = delete;
64  InfQuad4 (const InfQuad4 &) = delete;
65  InfQuad4 & operator= (const InfQuad4 &) = delete;
66  InfQuad4 & operator= (InfQuad4 &&) = delete;
67  virtual ~InfQuad4() = default;
68 
72  virtual unsigned int n_nodes() const override { return num_nodes; }
73 
77  virtual ElemType type () const override { return INFQUAD4; }
78 
82  virtual unsigned int n_sub_elem() const override { return 1; }
83 
88  virtual bool is_node_on_side(const unsigned int n,
89  const unsigned int s) const override;
90 
91  virtual std::vector<unsigned int> nodes_on_side(const unsigned int s) const override;
92 
93  virtual std::vector<unsigned int> nodes_on_edge(const unsigned int e) const override;
94 
99  virtual bool is_node_on_edge(const unsigned int n,
100  const unsigned int e) const override
101  { return this->is_node_on_side(n,e); }
102 
106  virtual Order default_order() const override;
107 
112  virtual std::unique_ptr<Elem> build_side_ptr (const unsigned int i) override;
113 
117  virtual void build_side_ptr (std::unique_ptr<Elem> & elem,
118  const unsigned int i) override;
119 
120  // Avoid hiding deprecated version with different signature
121  using Elem::build_side_ptr;
122 
123  virtual void connectivity(const unsigned int sf,
124  const IOPackage iop,
125  std::vector<dof_id_type> & conn) const override;
126 
133  virtual bool contains_point (const Point & p, Real tol=TOLERANCE) const override;
134 
138  static const int num_nodes = 4;
139  static const int nodes_per_side = 2;
140 
145  static const unsigned int side_nodes_map[num_sides][nodes_per_side];
146 
147  ElemType side_type (const unsigned int s) const override final;
148 
149  virtual void flip(BoundaryInfo *) override final;
150 
151 protected:
152 
157 
158 
159 
160 #ifdef LIBMESH_ENABLE_AMR
161 
165  virtual Real embedding_matrix (const unsigned int i,
166  const unsigned int j,
167  const unsigned int k) const override
168  { return _embedding_matrix[i][j][k]; }
169 
175 
177 
178 #endif // LIBMESH_ENABLE_AMR
179 
180 };
181 
182 
183 } // namespace libMesh
184 
185 #endif // ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
186 
187 #endif // LIBMESH_FACE_INF_QUAD4_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
virtual bool contains_point(const Point &p, Real tol=TOLERANCE) const override
virtual std::unique_ptr< Elem > build_side_ptr(const unsigned int i)=0
virtual ElemType type() const override
virtual ~InfQuad4()=default
static constexpr Real TOLERANCE
virtual bool is_node_on_edge(const unsigned int n, const unsigned int e) const override
The INFQUAD4 is an infinite element in 2D composed of 4 nodes.
static const unsigned int side_nodes_map[num_sides][nodes_per_side]
This maps the node of the side to element node numbers.
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.
IOPackage
libMesh interfaces with several different software packages for the purposes of creating, reading, and writing mesh files.
This is the base class from which all geometric element types are derived.
Definition: elem.h:94
static const int num_sides
Geometric constants for all InfQuads.
Definition: face_inf_quad.h:95
The libMesh namespace provides an interface to certain functionality in the library.
virtual std::vector< unsigned int > nodes_on_edge(const unsigned int e) const override
virtual bool is_node_on_side(const unsigned int n, const unsigned int s) const override
The BoundaryInfo class contains information relevant to boundary conditions including storing faces...
Definition: boundary_info.h:57
The InfQuad is an abstract element type that lives in two dimensions.
Definition: face_inf_quad.h:56
virtual void connectivity(const unsigned int sf, const IOPackage iop, std::vector< dof_id_type > &conn) const override
virtual unsigned int n_nodes() const override
Node * _nodelinks_data[num_nodes]
Data for links to nodes.
virtual Order default_order() const override
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
ElemType side_type(const unsigned int s) const override final
virtual void flip(BoundaryInfo *) override final
Flips the element (by swapping node and neighbor pointers) to have a mapping Jacobian of opposite sig...
static const int nodes_per_side
virtual std::vector< unsigned int > nodes_on_side(const unsigned int s) const override
static const int num_children
Definition: face_inf_quad.h:96
InfQuad4 & operator=(const InfQuad4 &)=delete
virtual unsigned int n_sub_elem() const override
static const int num_nodes
Geometric constants for InfQuad4.
A Point defines a location in LIBMESH_DIM dimensional Real space.
Definition: point.h:39
virtual std::unique_ptr< Elem > build_side_ptr(const unsigned int i) override
InfQuad4(Elem *p=nullptr)
Constructor.
static const Real _embedding_matrix[num_children][num_nodes][num_nodes]
Matrix that computes new nodal locations/solution values from current nodes/solution.