face_inf_quad4.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 #ifndef LIBMESH_FACE_INF_QUAD4_H 00020 #define LIBMESH_FACE_INF_QUAD4_H 00021 00022 00023 #include "libmesh/libmesh_config.h" 00024 00025 #ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS 00026 00027 // Local includes 00028 #include "libmesh/face_inf_quad.h" 00029 00030 // C++ includes 00031 #include <cstddef> 00032 00033 namespace libMesh 00034 { 00035 00036 00037 00053 // ------------------------------------------------------------ 00054 // InfQuad4 class definition 00055 class InfQuad4 : public InfQuad 00056 { 00057 public: 00058 00062 explicit 00063 InfQuad4 (Elem* p=NULL) : 00064 InfQuad(InfQuad4::n_nodes(), p, _nodelinks_data) {} 00065 00069 unsigned int n_nodes() const { return 4; } 00070 00074 ElemType type () const { return INFQUAD4; } 00075 00079 unsigned int n_sub_elem() const { return 1; } 00080 00084 virtual bool is_vertex(const unsigned int i) const; 00085 00089 virtual bool is_edge(const unsigned int i) const; 00090 00094 virtual bool is_face(const unsigned int i) const; 00095 00096 /* 00097 * @returns true iff the specified (local) node number is on the 00098 * specified side 00099 */ 00100 virtual bool is_node_on_side(const unsigned int n, 00101 const unsigned int s) const; 00102 00103 /* 00104 * @returns true iff the specified (local) node number is on the 00105 * specified edge (== is_node_on_side in 2D) 00106 */ 00107 virtual bool is_node_on_edge(const unsigned int n, 00108 const unsigned int e) const 00109 { return this->is_node_on_side(n,e); } 00110 00114 Order default_order() const { return FIRST; } 00115 00120 AutoPtr<Elem> build_side (const unsigned int i, 00121 bool proxy) const; 00122 00123 virtual void connectivity(const unsigned int sf, 00124 const IOPackage iop, 00125 std::vector<dof_id_type>& conn) const; 00126 00127 // void tecplot_connectivity(const unsigned int sf, 00128 // std::vector<unsigned int>& conn) const; 00129 00130 // void vtk_connectivity(const unsigned int sc, 00131 // std::vector<unsigned int>*conn = NULL) const; 00132 00133 // unsigned int vtk_element_type (const unsigned int) const 00134 // { return 9; } 00135 00142 bool contains_point (const Point& p, Real tol=TOLERANCE) const; 00143 00148 static const unsigned int side_nodes_map[3][2]; 00149 00150 00151 protected: 00152 00156 Node* _nodelinks_data[4]; 00157 00158 00159 00160 #ifdef LIBMESH_ENABLE_AMR 00161 00165 float embedding_matrix (const unsigned int i, 00166 const unsigned int j, 00167 const unsigned int k) const 00168 { return _embedding_matrix[i][j][k]; } 00169 00174 static const float _embedding_matrix[2][4][4]; 00175 00176 #endif 00177 00178 }; 00179 00180 00181 } // namespace libMesh 00182 00183 #endif // ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS 00184 00185 #endif // LIBMESH_FACE_INF_QUAD4_H
Site Created By: libMesh Developers
Last modified: February 05 2013 19:54:46 UTC
Hosted By: