20 #include "libmesh/libmesh_config.h" 22 #ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS 25 #include "libmesh/face_inf_quad6.h" 26 #include "libmesh/edge_edge3.h" 27 #include "libmesh/edge_inf_edge2.h" 28 #include "libmesh/enum_io_package.h" 29 #include "libmesh/enum_order.h" 54 const unsigned int s)
const 56 libmesh_assert_less (s,
n_sides());
65 libmesh_assert_less(s,
n_sides());
66 auto trim = (s == 0) ? 0 : 1;
76 #ifdef LIBMESH_ENABLE_AMR 83 { 1.0, 0.0, 0.0, 0.0, 0.0, 0.0 },
84 { 0.0, 0.0, 0.0, 0.0, 1.0, 0.0 },
85 { 0.0, 0.0, 1.0, 0.0, 0.0, 0.0 },
86 { 0.0, 0.0, 0.0, 0.0, 0.0, 1.0 },
87 { 0.375, -0.125, 0.0, 0.0, 0.75, 0.0 },
88 { 0.0, 0.0, 0.375, -0.125, 0.0, 0.75 }
94 { 0.0, 0.0, 0.0, 0.0, 1.0, 0.0 },
95 { 0.0, 1.0, 0.0, 0.0, 0.0, 0.0 },
96 { 0.0, 0.0, 0.0, 0.0, 0.0, 1.0 },
97 { 0.0, 0.0, 0.0, 1.0, 0.0, 0.0 },
98 { -0.125, 0.375, 0.0, 0.0, 0.75, 0.0 },
99 { 0.0, 0.0, -0.125, 0.375, 0.0, 0.75 }
117 libmesh_assert_less (s, this->
n_sides());
131 libmesh_error_msg(
"Invalid side s = " << s);
138 unsigned int side_node)
const 140 libmesh_assert_less (side, this->
n_sides());
152 std::unique_ptr<Elem> edge;
158 edge = std::make_unique<Edge3>();
166 edge = std::make_unique<InfEdge2>();
171 libmesh_error_msg(
"Invalid side i = " << i);
175 for (
auto n : edge->node_index_range())
178 edge->set_interior_parent(
this);
179 edge->inherit_data_from(*
this);
187 const unsigned int i)
189 libmesh_assert_less (i, this->
n_sides());
197 if (!side.get() || side->type() !=
EDGE3)
209 if (!side.get() || side->type() !=
INFEDGE2)
218 libmesh_error_msg(
"Invalid side i = " << i);
223 #ifdef LIBMESH_ENABLE_AMR 224 side->set_p_level(this->
p_level());
227 for (
auto n : side->node_index_range())
235 std::vector<dof_id_type> & conn)
const 267 libmesh_error_msg(
"Invalid sf = " << sf);
272 libmesh_error_msg(
"Unsupported IO package " << iop);
280 const unsigned int v)
const 282 libmesh_assert_greater_equal (n, this->
n_vertices());
283 libmesh_assert_less (n, this->
n_nodes());
284 libmesh_assert_less (v, 2);
298 std::pair<unsigned short int, unsigned short int>
301 libmesh_assert_greater_equal (n, this->
n_vertices());
302 libmesh_assert_less (n, this->
n_nodes());
304 return std::pair<unsigned short int, unsigned short int>
312 libmesh_assert_less (s, 3);
334 #endif // ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS virtual bool is_node_on_side(const unsigned int n, const unsigned int s) const override
ElemType
Defines an enum for geometric element types.
void swap2boundaryedges(unsigned short e1, unsigned short e2, BoundaryInfo *boundary_info) const
Swaps two edges in boundary_info, if it is non-null.
static const int nodes_per_side
static const unsigned int side_nodes_map[num_sides][nodes_per_side]
This maps the node of the side to element node numbers.
Order
defines an enum for polynomial orders.
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 dof_id_type key() const
virtual unsigned short int second_order_adjacent_vertex(const unsigned int n, const unsigned int v) const override
IOPackage
libMesh interfaces with several different software packages for the purposes of creating, reading, and writing mesh files.
void swap2boundarysides(unsigned short s1, unsigned short s2, BoundaryInfo *boundary_info) const
Swaps two sides in boundary_info, if it is non-null.
unsigned int p_level() const
static const int num_sides
Geometric constants for all InfQuads.
The libMesh namespace provides an interface to certain functionality in the library.
virtual unsigned int n_nodes() const override
virtual unsigned int n_vertices() const override final
ElemMappingType mapping_type() const
void swap2nodes(unsigned int n1, unsigned int n2)
Swaps two node_ptrs.
virtual std::vector< unsigned int > nodes_on_side(const unsigned int s) const override
virtual unsigned int n_sub_elem() const override
The BoundaryInfo class contains information relevant to boundary conditions including storing faces...
static const unsigned short int _second_order_adjacent_vertices[2][2]
Matrix that tells which vertices define the location of mid-side (or second-order) nodes...
void swap2neighbors(unsigned int n1, unsigned int n2)
Swaps two neighbor_ptrs.
static const int num_nodes
Geometric constants for InfQuad6.
virtual std::vector< unsigned int > nodes_on_edge(const unsigned int e) const override
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
subdomain_id_type subdomain_id() const
virtual unsigned int n_sides() const override final
virtual std::pair< unsigned short int, unsigned short int > second_order_child_vertex(const unsigned int n) const override
ElemType side_type(const unsigned int s) const override final
virtual unsigned int local_side_node(unsigned int side, unsigned int side_node) const override
static dof_id_type compute_key(dof_id_type n0)
static const int num_children
virtual void connectivity(const unsigned int sf, const IOPackage iop, std::vector< dof_id_type > &conn) const override
static const Real _embedding_matrix[num_children][num_nodes][num_nodes]
Matrix that computes new nodal locations/solution values from current nodes/solution.
dof_id_type node_id(const unsigned int i) const
virtual std::unique_ptr< Elem > build_side_ptr(const unsigned int i) override