MOAB: Mesh Oriented datABase
(version 5.2.1)
|
00001 /* ***************************************************************** 00002 MESQUITE -- The Mesh Quality Improvement Toolkit 00003 00004 Copyright 2006 Lawrence Livermore National Laboratory. Under 00005 the terms of Contract B545069 with the University of Wisconsin -- 00006 Madison, Lawrence Livermore National Laboratory retains certain 00007 rights in this software. 00008 00009 This library is free software; you can redistribute it and/or 00010 modify it under the terms of the GNU Lesser General Public 00011 License as published by the Free Software Foundation; either 00012 version 2.1 of the License, or (at your option) any later version. 00013 00014 This library is distributed in the hope that it will be useful, 00015 but WITHOUT ANY WARRANTY; without even the implied warranty of 00016 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00017 Lesser General Public License for more details. 00018 00019 You should have received a copy of the GNU Lesser General Public License 00020 (lgpl.txt) along with this library; if not, write to the Free Software 00021 Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA 00022 00023 (2006) kraftche@cae.wisc.edu 00024 00025 ***************************************************************** */ 00026 00027 #ifndef MSQ_LINEAR_QUADRILATERAL_HPP 00028 #define MSQ_LINEAR_QUADRILATERAL_HPP 00029 00030 /** \file LinearQuadrilateral.hpp 00031 * \brief Linear mapping funtion for Quadrilateral elements. 00032 * \author Jason Kraftcheck 00033 */ 00034 00035 #include "Mesquite.hpp" 00036 #include "MappingFunction.hpp" 00037 00038 namespace MBMesquite 00039 { 00040 00041 /**\brief Linear shape function for quadrilateral elements 00042 * 00043 * This class implements the MappingFunction interface, providing 00044 * a Linear shape function for quadrilateral elements. 00045 * 00046 * \f$\vec{x}(\xi,\eta) = \sum_{i=0}^{3} N_i(\xi,\eta) \vec{x_i}\f$ 00047 * 00048 * \f$N_0(\xi,\eta) = (1-\xi)(1-\eta)\f$ 00049 * 00050 * \f$N_1(\xi,\eta) = \xi (1-\eta)\f$ 00051 * 00052 * \f$N_2(\xi,\eta) = \xi \eta \f$ 00053 * 00054 * \f$N_3(\xi,\eta) = (1-\xi) \eta\f$ 00055 */ 00056 class MESQUITE_EXPORT LinearQuadrilateral : public MappingFunction2D 00057 { 00058 public: 00059 virtual EntityTopology element_topology() const; 00060 00061 virtual int num_nodes() const; 00062 00063 virtual void coefficients( Sample location, NodeSet nodeset, double* coeff_out, size_t* indices_out, 00064 size_t& num_coeff_out, MsqError& err ) const; 00065 00066 static void coefficients( Sample location, NodeSet nodeset, double* coeff_out, size_t* indices_out, 00067 size_t& num_coeff_out ); 00068 00069 virtual void derivatives( Sample location, NodeSet nodeset, size_t* vertex_indices_out, 00070 MsqVector< 2 >* d_coeff_d_xi_out, size_t& num_vtx, MsqError& err ) const; 00071 00072 static void derivatives( Sample location, NodeSet nodeset, size_t* vertex_indices_out, 00073 MsqVector< 2 >* d_coeff_d_xi_out, size_t& num_vtx ); 00074 00075 virtual void ideal( Sample location, MsqMatrix< 3, 2 >& jacobian_out, MsqError& err ) const; 00076 }; 00077 00078 } // namespace MBMesquite 00079 00080 #endif