cell_prism6.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_PRISM6_H 00021 #define LIBMESH_CELL_PRISM6_H 00022 00023 // Local includes 00024 #include "libmesh/cell_prism.h" 00025 00026 // C++ includes 00027 #include <cstddef> 00028 00029 namespace libMesh 00030 { 00031 00032 00033 00034 00053 // ------------------------------------------------------------ 00054 // Prism class definition 00055 class Prism6 : public Prism 00056 { 00057 public: 00058 00062 explicit 00063 Prism6 (Elem* p=NULL); 00064 00068 ElemType type () const { return PRISM6; } 00069 00073 unsigned int n_sub_elem() const { return 1; } 00074 00078 virtual bool is_vertex(const unsigned int i) const; 00079 00083 virtual bool is_edge(const unsigned int i) const; 00084 00088 virtual bool is_face(const unsigned int i) const; 00089 00090 /* 00091 * @returns true iff the specified (local) node number is on the 00092 * specified side 00093 */ 00094 virtual bool is_node_on_side(const unsigned int n, 00095 const unsigned int s) const; 00096 00097 /* 00098 * @returns true iff the specified (local) node number is on the 00099 * specified edge 00100 */ 00101 virtual bool is_node_on_edge(const unsigned int n, 00102 const unsigned int e) const; 00103 00104 /* 00105 * @returns true iff the element map is definitely affine within 00106 * numerical tolerances 00107 */ 00108 virtual bool has_affine_map () const; 00109 00113 Order default_order() const { return FIRST; } 00114 00119 AutoPtr<Elem> build_side (const unsigned int i, 00120 bool proxy) const; 00121 00126 AutoPtr<Elem> build_edge (const unsigned int i) const; 00127 00128 virtual void connectivity(const unsigned int sc, 00129 const IOPackage iop, 00130 std::vector<dof_id_type>& conn) const; 00131 00136 static const unsigned int side_nodes_map[5][4]; 00137 00141 static const unsigned int side_elems_map[5][4]; 00142 00147 static const unsigned int edge_nodes_map[9][2]; 00148 00152 virtual Real volume () const; 00153 00154 protected: 00155 00159 Node* _nodelinks_data[6]; 00160 00161 00162 00163 #ifdef LIBMESH_ENABLE_AMR 00164 00168 float embedding_matrix (const unsigned int i, 00169 const unsigned int j, 00170 const unsigned int k) const 00171 { return _embedding_matrix[i][j][k]; } 00172 00177 static const float _embedding_matrix[8][6][6]; 00178 00179 #endif 00180 00181 }; 00182 00183 00184 00185 // ------------------------------------------------------------ 00186 // Prism6 class member functions 00187 inline 00188 Prism6::Prism6(Elem* p) : 00189 Prism(Prism6::n_nodes(), p, _nodelinks_data) 00190 { 00191 } 00192 00193 00194 } // namespace libMesh 00195 00196 #endif // LIBMESH_CELL_PRISM6_H
Site Created By: libMesh Developers
Last modified: February 05 2013 19:54:45 UTC
Hosted By: