face_tri6.h
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2014 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
3 
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License, or (at your option) any later version.
8 
9 // This library is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // Lesser General Public License for more details.
13 
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 
18 
19 
20 #ifndef LIBMESH_FACE_TRI6_H
21 #define LIBMESH_FACE_TRI6_H
22 
23 
24 // Local includes
25 #include "libmesh/libmesh_common.h"
26 #include "libmesh/face_tri.h"
27 
28 // C++ includes
29 #include <cstddef>
30 
31 namespace libMesh
32 {
33 
34 
35 
36 // Forward declarations
37 
38 
39 
40 
57 // ------------------------------------------------------------
58 // Tri6 class definition
59 class Tri6 : public Tri
60 {
61 public:
62 
66  explicit
67  Tri6 (Elem* p=NULL) :
68  Tri(Tri6::n_nodes(), p, _nodelinks_data) {}
69 
73  ElemType type () const { return TRI6; }
74 
78  unsigned int n_nodes() const { return 6; }
79 
83  unsigned int n_sub_elem() const { return 4; }
84 
88  virtual bool is_vertex(const unsigned int i) const;
89 
93  virtual bool is_edge(const unsigned int i) const;
94 
98  virtual bool is_face(const unsigned int i) const;
99 
100  /*
101  * @returns true iff the specified (local) node number is on the
102  * specified side
103  */
104  virtual bool is_node_on_side(const unsigned int n,
105  const unsigned int s) const;
106 
107  /*
108  * @returns true iff the specified (local) node number is on the
109  * specified edge (== is_node_on_side in 2D)
110  */
111  virtual bool is_node_on_edge(const unsigned int n,
112  const unsigned int e) const
113  { return this->is_node_on_side(n,e); }
114 
115  /*
116  * @returns true iff the element map is definitely affine within
117  * numerical tolerances
118  */
119  virtual bool has_affine_map () const;
120 
124  Order default_order() const { return SECOND; }
125 
135  dof_id_type key (const unsigned int s) const;
136 
137  AutoPtr<Elem> build_side (const unsigned int i,
138  bool proxy) const;
139 
140  virtual void connectivity(const unsigned int sf,
141  const IOPackage iop,
142  std::vector<dof_id_type>& conn) const;
143 
147  unsigned int n_second_order_adjacent_vertices (const unsigned int) const
148  { return 2; }
149 
155  unsigned short int second_order_adjacent_vertex (const unsigned int n,
156  const unsigned int v) const;
157 
167  virtual std::pair<unsigned short int, unsigned short int>
168  second_order_child_vertex (const unsigned int n) const;
169 
174  static const unsigned int side_nodes_map[3][3];
175 
176 
177 protected:
178 
183 
184 
185 
186 #ifdef LIBMESH_ENABLE_AMR
187 
191  float embedding_matrix (const unsigned int i,
192  const unsigned int j,
193  const unsigned int k) const
194  { return _embedding_matrix[i][j][k]; }
195 
200  static const float _embedding_matrix[4][6][6];
201 
202 #endif
203 
204 
205 private:
206 
211  static const unsigned short int _second_order_adjacent_vertices[3][2];
212 
216  static const unsigned short int _second_order_vertex_child_number[6];
217 
221  static const unsigned short int _second_order_vertex_child_index[6];
222 };
223 
224 
225 } // namespace libMesh
226 
227 #endif // LIBMESH_FACE_TRI6_H

Site Created By: libMesh Developers
Last modified: February 07 2014 16:57:04 UTC

Hosted By:
SourceForge.net Logo