cell_tet10.h
Go to the documentation of this file.00001 // The libMesh Finite Element Library. 00002 // Copyright (C) 2002-2012 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner 00003 00004 // This library is free software; you can redistribute it and/or 00005 // modify it under the terms of the GNU Lesser General Public 00006 // License as published by the Free Software Foundation; either 00007 // version 2.1 of the License, or (at your option) any later version. 00008 00009 // This library is distributed in the hope that it will be useful, 00010 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00011 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00012 // Lesser General Public License for more details. 00013 00014 // You should have received a copy of the GNU Lesser General Public 00015 // License along with this library; if not, write to the Free Software 00016 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA 00017 00018 00019 00020 #ifndef LIBMESH_CELL_TET10_H 00021 #define LIBMESH_CELL_TET10_H 00022 00023 // Local includes 00024 #include "libmesh/cell_tet.h" 00025 00026 // C++ includes 00027 #include <cstddef> 00028 00029 namespace libMesh 00030 { 00031 00032 00033 00034 00061 // ------------------------------------------------------------ 00062 // Tet10 class definition 00063 class Tet10 : public Tet 00064 { 00065 public: 00066 00070 explicit 00071 Tet10 (Elem* p=NULL); 00072 00076 ElemType type () const { return TET10; } 00077 00081 unsigned int n_nodes() const { return 10; } 00082 00086 unsigned int n_sub_elem() const { return 8; } 00087 00091 virtual bool is_vertex(const unsigned int i) const; 00092 00096 virtual bool is_edge(const unsigned int i) const; 00097 00101 virtual bool is_face(const unsigned int i) const; 00102 00103 /* 00104 * @returns true iff the specified (local) node number is on the 00105 * specified side 00106 */ 00107 virtual bool is_node_on_side(const unsigned int n, 00108 const unsigned int s) const; 00109 00110 /* 00111 * @returns true iff the specified (local) node number is on the 00112 * specified edge 00113 */ 00114 virtual bool is_node_on_edge(const unsigned int n, 00115 const unsigned int e) const; 00116 00117 /* 00118 * @returns true iff the specified child is on the 00119 * specified side 00120 */ 00121 virtual bool is_child_on_side(const unsigned int c, 00122 const unsigned int s) const; 00123 00124 /* 00125 * @returns true iff the element map is definitely affine within 00126 * numerical tolerances 00127 */ 00128 virtual bool has_affine_map () const; 00129 00133 Order default_order() const { return SECOND; } 00134 00139 AutoPtr<Elem> build_side (const unsigned int i, 00140 bool proxy) const; 00141 00146 AutoPtr<Elem> build_edge (const unsigned int i) const; 00147 00148 virtual void connectivity(const unsigned int sc, 00149 const IOPackage iop, 00150 std::vector<dof_id_type>& conn) const; 00151 00155 unsigned int n_second_order_adjacent_vertices (const unsigned int) const 00156 { return 2; } 00157 00163 unsigned short int second_order_adjacent_vertex (const unsigned int n, 00164 const unsigned int v) const; 00165 00175 virtual std::pair<unsigned short int, unsigned short int> 00176 second_order_child_vertex (const unsigned int n) const; 00177 00182 static const unsigned int side_nodes_map[4][6]; 00183 00188 static const unsigned int edge_nodes_map[6][3]; 00189 00190 protected: 00191 00195 Node* _nodelinks_data[10]; 00196 00197 00198 00199 #ifdef LIBMESH_ENABLE_AMR 00200 00204 float embedding_matrix (const unsigned int i, 00205 const unsigned int j, 00206 const unsigned int k) const; 00207 00212 static const float _embedding_matrix[8][10][10]; 00213 00214 00215 #endif 00216 00217 00218 private: 00219 00224 static const unsigned short int _second_order_adjacent_vertices[6][2]; 00225 00229 static const unsigned short int _second_order_vertex_child_number[10]; 00230 00234 static const unsigned short int _second_order_vertex_child_index[10]; 00235 }; 00236 00237 00238 00239 // ------------------------------------------------------------ 00240 // Tet10 class member functions 00241 inline 00242 Tet10::Tet10(Elem* p) : 00243 Tet(Tet10::n_nodes(), p, _nodelinks_data) 00244 { 00245 } 00246 00247 } // namespace libMesh 00248 00249 00250 #endif // LIBMESH_CELL_TET10_H
Site Created By: libMesh Developers
Last modified: February 05 2013 19:54:45 UTC
Hosted By: