18 #include "libmesh/libmesh_config.h" 20 #ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS 23 #include "libmesh/cell_inf_prism.h" 24 #include "libmesh/cell_inf_prism6.h" 25 #include "libmesh/face_tri3.h" 26 #include "libmesh/face_inf_quad4.h" 27 #include "libmesh/fe_type.h" 28 #include "libmesh/fe_interface.h" 29 #include "libmesh/inf_fe_map.h" 30 #include "libmesh/enum_order.h" 71 libmesh_assert_less (s, this->
n_sides());
89 libmesh_error_msg(
"Invalid side s = " << s);
97 libmesh_assert_less (s, this->
n_sides());
115 libmesh_error_msg(
"Invalid side s = " << s);
122 unsigned int side_node)
const 124 libmesh_assert_less (side, this->
n_sides());
138 unsigned int edge_node)
const 140 libmesh_assert_less(edge, this->
n_edges());
150 libmesh_assert_less (i, this->
n_sides());
152 std::unique_ptr<Elem> face;
158 face = std::make_unique<Tri3>();
166 face = std::make_unique<InfQuad4>();
171 libmesh_error_msg(
"Invalid side i = " << i);
175 for (
auto n : face->node_index_range())
184 const unsigned int i)
186 libmesh_assert_less (i, this->
n_sides());
193 if (!side.get() || side->type() !=
TRI3)
207 if (!side.get() || side->type() !=
INFQUAD4)
216 libmesh_error_msg(
"Invalid side i = " << i);
222 for (
auto n : side->node_index_range())
229 const unsigned int s)
const 232 libmesh_assert_less (s, this->
n_sides());
234 return (s == 0 || c+1 == s || c == s%3);
240 const unsigned int s)
const 242 libmesh_assert_less (e, this->
n_edges());
243 libmesh_assert_less (s, this->
n_sides());
269 const Real tmp_min_distance_sq = std::min(pt0_o.
norm_sq(),
277 const Real min_distance_sq = tmp_min_distance_sq
284 const Real conservative_p_dist_sq = 1.01 * (
Point(p - my_origin).
norm_sq());
286 if (conservative_p_dist_sq < min_distance_sq)
294 Point p_o(p - my_origin);
295 pt0_o /= pt0_o.
norm();
296 pt1_o /= pt1_o.
norm();
297 pt2_o /= pt2_o.
norm();
324 std::vector<unsigned>
327 libmesh_assert_less(e,
n_edges());
345 #endif // ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS class FEType hides (possibly multiple) FEFamily and approximation orders, thereby enabling specialize...
auto norm() const -> decltype(std::norm(T()))
static const Real _master_points[12][3]
Master element node locations.
virtual dof_id_type key() const
virtual unsigned int n_edges() const override final
The libMesh namespace provides an interface to certain functionality in the library.
virtual bool is_flipped() const override final
static Point inverse_map(const unsigned int dim, const Elem *elem, const Point &p, const Real tolerance=TOLERANCE, const bool secure=true)
virtual unsigned int n_children() const override final
std::vector< unsigned int > sides_on_edge(const unsigned int e) const override final
virtual unsigned int local_side_node(unsigned int side, unsigned int side_node) const override
T triple_product(const TypeVector< T > &a, const TypeVector< T > &b, const TypeVector< T > &c)
auto norm_sq() const -> decltype(std::norm(T()))
static const unsigned int side_nodes_map[num_sides][nodes_per_side]
This maps the node of the side to element node numbers.
virtual unsigned int local_edge_node(unsigned int edge, unsigned int edge_node) const override
virtual std::unique_ptr< Elem > side_ptr(const unsigned int i) override final
virtual bool is_child_on_side(const unsigned int c, const unsigned int s) const override final
static bool on_reference_element(const Point &p, const ElemType t, const Real eps=TOLERANCE)
virtual bool is_edge_on_side(const unsigned int e, const unsigned int s) const override final
static const unsigned int edge_nodes_map[num_edges][nodes_per_edge]
This maps the node of the edge to element node numbers.
virtual Point origin() const override
static const unsigned int edge_sides_map[6][2]
This maps each edge to the sides that contain said edge.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
subdomain_id_type subdomain_id() const
static const int nodes_per_side
const Node * node_ptr(const unsigned int i) const
virtual unsigned short dim() const override
virtual bool contains_point(const Point &p, Real tol=TOLERANCE) const override
static const int nodes_per_edge
static dof_id_type compute_key(dof_id_type n0)
virtual Order default_order() const =0
virtual unsigned int n_sides() const override final
virtual ElemType type() const =0
A Point defines a location in LIBMESH_DIM dimensional Real space.
dof_id_type node_id(const unsigned int i) const
const Point & point(const unsigned int i) const
virtual dof_id_type low_order_key(const unsigned int s) const override