21 #ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS 91 const unsigned int s)
const 93 libmesh_assert_less (s,
n_sides());
102 libmesh_assert_less(s,
n_sides());
107 const unsigned int e)
const 109 libmesh_assert_less (e,
n_edges());
127 libmesh_assert_less (i, this->
n_sides());
135 return libmesh_make_unique<Side<Quad4,InfHex8>>(
this,i);
142 return libmesh_make_unique<Side<InfQuad4,InfHex8>>(
this,i);
145 libmesh_error_msg(
"Invalid side i = " << i);
152 std::unique_ptr<Elem> face;
159 face = libmesh_make_unique<Quad4>();
169 face = libmesh_make_unique<InfQuad4>();
174 libmesh_error_msg(
"Invalid side i = " << i);
180 for (
unsigned n=0; n<face->n_nodes(); ++n)
181 face->set_node(n) = this->
node_ptr(InfHex8::side_nodes_map[i][n]);
189 const unsigned int i)
198 libmesh_assert_less (i, this->
n_edges());
201 return libmesh_make_unique<SideEdge<Edge2,InfHex8>>(
this,i);
204 return libmesh_make_unique<SideEdge<InfEdge2,InfHex8>>(
this,i);
210 std::vector<dof_id_type> & conn)
const 233 libmesh_error_msg(
"Unsupported IO package " << iop);
239 #ifdef LIBMESH_ENABLE_AMR 246 { 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
247 { 0.5, 0.5, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
248 { 0.25, 0.25, 0.25, 0.25, 0.0, 0.0, 0.0, 0.0},
249 { 0.5, 0.0, 0.0, 0.5, 0.0, 0.0, 0.0, 0.0},
250 { 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0},
251 { 0.0, 0.0, 0.0, 0.0, 0.5, 0.5, 0.0, 0.0},
252 { 0.0, 0.0, 0.0, 0.0, 0.25, 0.25, 0.25, 0.25},
253 { 0.0, 0.0, 0.0, 0.0, 0.5, 0.0, 0.0, 0.5}
259 { 0.5, 0.5, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
260 { 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
261 { 0.0, 0.5, 0.5, 0.0, 0.0, 0.0, 0.0, 0.0},
262 { 0.25, 0.25, 0.25, 0.25, 0.0, 0.0, 0.0, 0.0},
263 { 0.0, 0.0, 0.0, 0.0, 0.5, 0.5, 0.0, 0.0},
264 { 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0},
265 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.5, 0.5, 0.0},
266 { 0.0, 0.0, 0.0, 0.0, 0.25, 0.25, 0.25, 0.25}
272 { 0.5, 0.0, 0.0, 0.5, 0.0, 0.0, 0.0, 0.0},
273 { 0.25, 0.25, 0.25, 0.25, 0.0, 0.0, 0.0, 0.0},
274 { 0.0, 0.0, 0.5, 0.5, 0.0, 0.0, 0.0, 0.0},
275 { 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0},
276 { 0.0, 0.0, 0.0, 0.0, 0.5, 0.0, 0.0, 0.5},
277 { 0.0, 0.0, 0.0, 0.0, 0.25, 0.25, 0.25, 0.25},
278 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.5, 0.5},
279 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0}
285 { 0.25, 0.25, 0.25, 0.25, 0.0, 0.0, 0.0, 0.0},
286 { 0.0, 0.5, 0.5, 0.0, 0.0, 0.0, 0.0, 0.0},
287 { 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0},
288 { 0.0, 0.0, 0.5, 0.5, 0.0, 0.0, 0.0, 0.0},
289 { 0.0, 0.0, 0.0, 0.0, 0.25, 0.25, 0.25, 0.25},
290 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.5, 0.5, 0.0},
291 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0},
292 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.5, 0.5}
302 #endif // ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS virtual bool is_node_on_edge(const unsigned int n, const unsigned int e) const override
static const float _embedding_matrix[num_children][num_nodes][num_nodes]
virtual unsigned int n_sub_elem() const override
static const int nodes_per_side
virtual bool is_edge(const unsigned int i) const override
virtual std::unique_ptr< Elem > build_edge_ptr(const unsigned int i) override
virtual std::unique_ptr< Elem > side_ptr(const unsigned int i) override final
virtual Order default_order() const override
virtual bool is_face(const unsigned int i) const override
virtual void connectivity(const unsigned int sc, const IOPackage iop, std::vector< dof_id_type > &conn) const override
static const int num_children
virtual std::unique_ptr< Elem > build_side_ptr(const unsigned int i, bool proxy) override
virtual std::vector< unsigned int > nodes_on_side(const unsigned int s) const override
static const unsigned int side_nodes_map[num_sides][nodes_per_side]
static const int num_sides
subdomain_id_type subdomain_id() const
const Node * node_ptr(const unsigned int i) const
static const unsigned int edge_nodes_map[num_edges][nodes_per_edge]
virtual unsigned int n_edges() const override final
static const int num_edges
virtual bool is_vertex(const unsigned int i) const override
static const int num_nodes
dof_id_type node_id(const unsigned int i) const
static const int nodes_per_edge
virtual unsigned int n_sides() const override final
virtual bool is_node_on_side(const unsigned int n, const unsigned int s) const override