19 #include "libmesh/side.h" 20 #include "libmesh/edge_edge3.h" 21 #include "libmesh/face_quad8.h" 22 #include "libmesh/enum_io_package.h" 23 #include "libmesh/enum_order.h" 47 #ifdef LIBMESH_ENABLE_AMR 54 { 1.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000 },
55 { 0.00000, 0.00000, 0.00000, 0.00000, 1.00000, 0.00000, 0.00000, 0.00000 },
56 { -0.250000, -0.250000, -0.250000, -0.250000, 0.500000, 0.500000, 0.500000, 0.500000 },
57 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 1.00000 },
58 { 0.375000, -0.125000, 0.00000, 0.00000, 0.750000, 0.00000, 0.00000, 0.00000 },
59 { -0.187500, -0.187500, -0.187500, -0.187500, 0.750000, 0.375000, 0.250000, 0.375000 },
60 { -0.187500, -0.187500, -0.187500, -0.187500, 0.375000, 0.250000, 0.375000, 0.750000 },
61 { 0.375000, 0.00000, 0.00000, -0.125000, 0.00000, 0.00000, 0.00000, 0.750000 }
67 { 0.00000, 0.00000, 0.00000, 0.00000, 1.00000, 0.00000, 0.00000, 0.00000 },
68 { 0.00000, 1.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000 },
69 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 1.00000, 0.00000, 0.00000 },
70 { -0.250000, -0.250000, -0.250000, -0.250000, 0.500000, 0.500000, 0.500000, 0.500000 },
71 { -0.125000, 0.375000, 0.00000, 0.00000, 0.750000, 0.00000, 0.00000, 0.00000 },
72 { 0.00000, 0.375000, -0.125000, 0.00000, 0.00000, 0.750000, 0.00000, 0.00000 },
73 { -0.187500, -0.187500, -0.187500, -0.187500, 0.375000, 0.750000, 0.375000, 0.250000 },
74 { -0.187500, -0.187500, -0.187500, -0.187500, 0.750000, 0.375000, 0.250000, 0.375000 }
80 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 1.00000 },
81 { -0.250000, -0.250000, -0.250000, -0.250000, 0.500000, 0.500000, 0.500000, 0.500000 },
82 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 1.00000, 0.00000 },
83 { 0.00000, 0.00000, 0.00000, 1.00000, 0.00000, 0.00000, 0.00000, 0.00000 },
84 { -0.187500, -0.187500, -0.187500, -0.187500, 0.375000, 0.250000, 0.375000, 0.750000 },
85 { -0.187500, -0.187500, -0.187500, -0.187500, 0.250000, 0.375000, 0.750000, 0.375000 },
86 { 0.00000, 0.00000, -0.125000, 0.375000, 0.00000, 0.00000, 0.750000, 0.00000 },
87 { -0.125000, 0.00000, 0.00000, 0.375000, 0.00000, 0.00000, 0.00000, 0.750000 }
93 { -0.250000, -0.250000, -0.250000, -0.250000, 0.500000, 0.500000, 0.500000, 0.500000 },
94 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 1.00000, 0.00000, 0.00000 },
95 { 0.00000, 0.00000, 1.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000 },
96 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 1.00000, 0.00000 },
97 { -0.187500, -0.187500, -0.187500, -0.187500, 0.375000, 0.750000, 0.375000, 0.250000 },
98 { 0.00000, -0.125000, 0.375000, 0.00000, 0.00000, 0.750000, 0.00000, 0.00000 },
99 { 0.00000, 0.00000, 0.375000, -0.125000, 0.00000, 0.00000, 0.750000, 0.00000 },
100 { -0.187500, -0.187500, -0.187500, -0.187500, 0.250000, 0.375000, 0.750000, 0.375000 }
131 const unsigned int s)
const 133 libmesh_assert_less (s,
n_sides());
139 std::vector<unsigned>
142 libmesh_assert_less(s,
n_sides());
146 std::vector<unsigned>
181 libmesh_assert_less (s, this->
n_sides());
206 libmesh_error_msg(
"Invalid side s = " << s);
213 unsigned int side_node)
const 215 libmesh_assert_less (side, this->
n_sides());
226 return this->simple_build_side_ptr<Edge3, Quad8>(i, proxy);
232 const unsigned int i)
234 this->simple_build_side_ptr<Quad8>(side, i,
EDGE3);
244 std::vector<dof_id_type> & conn)
const 307 libmesh_error_msg(
"Invalid sf = " << sf);
323 libmesh_error_msg(
"Unsupported IO package " << iop);
340 for (
unsigned d=0; d<LIBMESH_DIM; ++d)
343 for (
unsigned int p=1; p != 8; ++p)
344 center += this->
point(p)(d);
348 for (
unsigned int p=0; p != 8; ++p)
351 pmin(d) = center - hd;
352 pmax(d) = center + hd;
382 a1 = -x0/4 + x1/4 + x2/4 - x3/4 - x5/2 + x7/2,
383 b1 = -x0/2 - x1/2 + x2/2 + x3/2 + x4 - x6,
384 c1 = x0/2 + x1/2 + x2/2 + x3/2 - x4 - x6,
385 d1 = x0/4 - x1/4 + x2/4 - x3/4,
387 a2 = -x0/4 - x1/4 + x2/4 + x3/4 + x4/2 - x6/2,
388 b2 = -x0/2 + x1/2 + x2/2 - x3/2 - x5 + x7,
389 c2 = x0/4 - x1/4 + x2/4 - x3/4,
390 d2 = x0/2 + x1/2 + x2/2 + x3/2 - x5 - x7,
394 const unsigned int N = 3;
396 const Real w[N] = {5./9, 8./9, 5./9};
399 for (
unsigned int i=0; i<N; ++i)
400 for (
unsigned int j=0; j<N; ++j)
401 vol += w[i] * w[j] *
cross_norm(q[j]*q[j]*a1 + q[i]*q[j]*b1 + q[i]*c1 + q[j]*d1 + e1,
402 q[i]*q[i]*a2 + q[i]*q[j]*b2 + q[i]*c2 + q[j]*d2 + e2);
410 const unsigned int v)
const 412 libmesh_assert_greater_equal (n, this->
n_vertices());
413 libmesh_assert_less (n, this->
n_nodes());
414 libmesh_assert_less (v, 2);
421 std::pair<unsigned short int, unsigned short int>
424 libmesh_assert_greater_equal (n, this->
n_vertices());
425 libmesh_assert_less (n, this->
n_nodes());
431 return std::pair<unsigned short int, unsigned short int>
439 libmesh_assert_less (perm_num, 4);
441 for (
unsigned int i = 0; i != perm_num; ++i)
471 libmesh_assert_less (s, 4);
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.
virtual unsigned int n_vertices() const override final
Order
defines an enum for polynomial orders.
static const int num_sides
static const int nodes_per_side
T cross_norm(const TypeVector< T > &b, const TypeVector< T > &c)
Calls cross_norm_sq() and takes the square root of the result.
virtual BoundingBox loose_bounding_box() const override
IOPackage
libMesh interfaces with several different software packages for the purposes of creating, reading, and writing mesh files.
virtual std::unique_ptr< Elem > build_side_ptr(const unsigned int i, bool proxy=false) override
virtual std::vector< unsigned int > nodes_on_side(const unsigned int s) const override
void swap2boundarysides(unsigned short s1, unsigned short s2, BoundaryInfo *boundary_info) const
Swaps two sides in boundary_info, if it is non-null.
virtual std::pair< unsigned short int, unsigned short int > second_order_child_vertex(const unsigned int n) const override
virtual dof_id_type key() const override
virtual unsigned int n_nodes() const override
ADRealEigenVector< T, D, asd > sqrt(const ADRealEigenVector< T, D, asd > &)
The libMesh namespace provides an interface to certain functionality in the library.
virtual void flip(BoundaryInfo *) override final
Flips the element (by swapping node and neighbor pointers) to have a mapping Jacobian of opposite sig...
virtual bool is_vertex(const unsigned int i) const override
ADRealEigenVector< T, D, asd > abs(const ADRealEigenVector< T, D, asd > &)
virtual Order default_order() const override
void swap4nodes(unsigned int n1, unsigned int n2, unsigned int n3, unsigned int n4)
Swaps four node_ptrs, "rotating" them.
virtual bool is_face(const unsigned int i) const override
ElemMappingType mapping_type() const
void swap2nodes(unsigned int n1, unsigned int n2)
Swaps two node_ptrs.
static const int num_nodes
Geometric constants for Quad8.
virtual void connectivity(const unsigned int sf, const IOPackage iop, std::vector< dof_id_type > &conn) const override
static const unsigned short int _second_order_vertex_child_index[9]
Vector that names the child vertex index for each second order node.
The BoundaryInfo class contains information relevant to boundary conditions including storing faces...
virtual std::vector< unsigned int > nodes_on_edge(const unsigned int e) const override
static const unsigned short int _second_order_vertex_child_number[9]
Vector that names a child sharing each second order node.
static constexpr Real affine_tol
Default tolerance to use in has_affine_map().
static const unsigned short int _second_order_adjacent_vertices[4][2]
Matrix that tells which vertices define the location of mid-side (or second-order) nodes...
virtual unsigned int n_sub_elem() const override
virtual unsigned int n_sides() const override final
void swap2neighbors(unsigned int n1, unsigned int n2)
Swaps two neighbor_ptrs.
virtual bool is_edge(const unsigned int i) const override
Defines a Cartesian bounding box by the two corner extremum.
virtual void permute(unsigned int perm_num) override final
Permutes the element (by swapping node and neighbor pointers) according to the specified index...
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
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_node_on_side(const unsigned int n, const unsigned int s) const override
virtual Real volume() const
void swap4neighbors(unsigned int n1, unsigned int n2, unsigned int n3, unsigned int n4)
Swaps four neighbor_ptrs, "rotating" them.
virtual bool has_affine_map() const override
virtual unsigned short int second_order_adjacent_vertex(const unsigned int n, const unsigned int v) const override
static const int num_children
unsigned int center_node_on_side(const unsigned short side) const override final
static dof_id_type compute_key(dof_id_type n0)
virtual Real volume() const override
An optimized method for approximating the area of a QUAD8 using quadrature.
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
bool relative_fuzzy_equals(const TypeVector< T > &rhs, Real tol=TOLERANCE) const
auto index_range(const T &sizable)
Helper function that returns an IntRange<std::size_t> representing all the indices of the passed-in v...
ElemType side_type(const unsigned int s) const override final
static const Real _embedding_matrix[num_children][num_nodes][num_nodes]
Matrix that computes new nodal locations/solution values from current nodes/solution.
virtual unsigned int local_side_node(unsigned int side, unsigned int side_node) const override