45 #ifdef LIBMESH_ENABLE_AMR 105 const unsigned int s)
const 107 libmesh_assert_less (s,
n_sides());
113 std::vector<unsigned>
116 libmesh_assert_less(s,
n_sides());
128 libmesh_assert_less (i, this->
n_sides());
131 return libmesh_make_unique<Side<Edge2,Tri3>>(
this,i);
135 std::unique_ptr<Elem> edge = libmesh_make_unique<Edge2>();
139 for (
unsigned n=0; n<edge->n_nodes(); ++n)
140 edge->set_node(n) = this->
node_ptr(Tri3::side_nodes_map[i][n]);
149 const unsigned int i)
151 this->simple_build_side_ptr<Tri3>(
side, i,
EDGE2);
158 std::vector<dof_id_type> & conn)
const 185 libmesh_error_msg(
"Unsupported IO package " << iop);
215 theta0=std::acos(( v10*v20)/len_10/len_20),
216 theta1=std::acos((-v10*v21)/len_10/len_21),
220 libmesh_assert_greater (theta0, 0.);
221 libmesh_assert_greater (theta1, 0.);
222 libmesh_assert_greater (theta2, 0.);
239 Real dot00 = v0 * v0;
240 Real dot01 = v0 * v1;
241 Real dot02 = v0 * v2;
242 Real dot11 = v1 * v1;
243 Real dot12 = v1 * v2;
250 Real invDenom = 1 / (dot00 * dot11 - dot01 * dot01);
251 Real u = (dot11 * dot02 - dot01 * dot12) * invDenom;
252 Real v = (dot00 * dot12 - dot01 * dot02) * invDenom;
255 return (u > -tol) && (v > -tol) && (u + v < 1 + tol);
static const int num_nodes
virtual unsigned int n_sub_elem() const override
static const unsigned int side_nodes_map[num_sides][nodes_per_side]
static const int num_sides
T cross_norm(const TypeVector< T > &b, const TypeVector< T > &c)
virtual bool is_node_on_side(const unsigned int n, const unsigned int s) const override
virtual BoundingBox loose_bounding_box() const
std::pair< Real, Real > min_and_max_angle() const
virtual std::vector< unsigned int > nodes_on_side(const unsigned int s) const override
virtual BoundingBox loose_bounding_box() const override
virtual void connectivity(const unsigned int sf, const IOPackage iop, std::vector< dof_id_type > &conn) const override
long double max(long double a, double b)
virtual bool is_edge(const unsigned int i) const override
virtual unsigned int n_sides() const override final
T triple_product(const TypeVector< T > &a, const TypeVector< T > &b, const TypeVector< T > &c)
virtual bool contains_point(const Point &p, Real tol) const override
static const int nodes_per_side
virtual std::unique_ptr< Elem > build_side_ptr(const unsigned int i, bool proxy) override
virtual bool is_vertex(const unsigned int i) const override
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
subdomain_id_type subdomain_id() const
const Node * node_ptr(const unsigned int i) const
virtual Order default_order() const override
static const int num_children
virtual bool is_face(const unsigned int i) const override
virtual Real volume() const override
long double min(long double a, double b)
A geometric point in (x,y,z) space.
dof_id_type node_id(const unsigned int i) const
const Point & point(const unsigned int i) const
static const float _embedding_matrix[num_children][num_nodes][num_nodes]
std::unique_ptr< Elem > side(const unsigned int i) const