19 #include "libmesh/face_tri.h" 20 #include "libmesh/edge_edge2.h" 21 #include "libmesh/face_tri3.h" 22 #include "libmesh/enum_elem_quality.h" 56 libmesh_assert_less (s, this->
n_sides());
66 libmesh_assert_less (s, this->
n_sides());
75 unsigned int side_node)
const 77 libmesh_assert_less (side, this->
n_sides());
86 unsigned int edge_node)
const 104 libmesh_assert_less (i, this->
n_sides());
106 std::unique_ptr<Elem> edge = std::make_unique<Edge2>();
108 for (
auto n : edge->node_index_range())
117 const unsigned int i)
119 this->simple_side_ptr<Tri,Tri3>(side, i,
EDGE2);
125 const unsigned int s)
const 128 libmesh_assert_less (s, this->
n_sides());
130 return (c == s || c == (s+1)%3);
140 !this->
point(2)(2) &&
173 if ((l1 <=0.) || (l2 <= 0.) || (l3 <= 0.))
176 const Real s1 = std::sin(std::acos(v1*v2/l1/l2)/2.);
178 const Real s2 = std::sin(std::acos(v1*v3/l1/l3)/2.);
179 const Real s3 = std::sin(std::acos(v2*v3/l2/l3)/2.);
181 return 8. * s1 * s2 * s3;
200 std::array<Real, 6> A =
210 lambda11 = A[0]*A[0] + A[2]*A[2] + A[4]*A[4],
211 lambda12 = A[0]*A[1] + A[2]*A[3] + A[4]*A[5],
212 lambda22 = A[1]*A[1] + A[3]*A[3] + A[5]*A[5];
216 Real den = lambda11 + lambda22 - lambda12;
221 Real alpha =
std::sqrt(lambda11 * lambda22 - lambda12 * lambda12);
242 std::pair<Real, Real> bounds;
264 bounds.second = 1.155;
279 libMesh::out <<
"Warning: Invalid quality measure chosen." << std::endl;
virtual std::unique_ptr< Elem > side_ptr(const unsigned int i) override final
static const Real _master_points[6][3]
Master element node locations.
auto norm() const -> decltype(std::norm(T()))
virtual unsigned int local_edge_node(unsigned int edge, unsigned int edge_node) const override
Calls local_side_node(edge, edge_node).
static const unsigned int side_nodes_map[num_sides][nodes_per_side]
This maps the node of the side to element node numbers.
virtual bool is_flipped() const override final
ADRealEigenVector< T, D, asd > sqrt(const ADRealEigenVector< T, D, asd > &)
The libMesh namespace provides an interface to certain functionality in the library.
virtual unsigned int n_sides() const override final
virtual bool is_child_on_side(const unsigned int c, const unsigned int s) const override final
virtual dof_id_type key() const override final
virtual dof_id_type low_order_key(const unsigned int s) const override
ElemQuality
Defines an enum for element quality metrics.
virtual unsigned int local_side_node(unsigned int side, unsigned int side_node) const override
static const int nodes_per_side
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const Node * node_ptr(const unsigned int i) const
virtual Real quality(const ElemQuality q) const
virtual std::pair< Real, Real > qual_bounds(const ElemQuality q) const override
virtual unsigned int n_children() const override final
static dof_id_type compute_key(dof_id_type n0)
virtual Real quality(const ElemQuality q) const override
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