Branch data     Line data    Source code

       1                 :            : /* *****************************************************************
       2                 :            :     MESQUITE -- The Mesh Quality Improvement Toolkit
       3                 :            : 
       4                 :            :     Copyright 2006 Lawrence Livermore National Laboratory.  Under
       5                 :            :     the terms of Contract B545069 with the University of Wisconsin --
       6                 :            :     Madison, Lawrence Livermore National Laboratory retains certain
       7                 :            :     rights in this software.
       8                 :            : 
       9                 :            :     This library is free software; you can redistribute it and/or
      10                 :            :     modify it under the terms of the GNU Lesser General Public
      11                 :            :     License as published by the Free Software Foundation; either
      12                 :            :     version 2.1 of the License, or (at your option) any later version.
      13                 :            : 
      14                 :            :     This library is distributed in the hope that it will be useful,
      15                 :            :     but WITHOUT ANY WARRANTY; without even the implied warranty of
      16                 :            :     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
      17                 :            :     Lesser General Public License for more details.
      18                 :            : 
      19                 :            :     You should have received a copy of the GNU Lesser General Public License
      20                 :            :     (lgpl.txt) along with this library; if not, write to the Free Software
      21                 :            :     Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
      22                 :            : 
      23                 :            :     (2006) [email protected]
      24                 :            : 
      25                 :            :   ***************************************************************** */
      26                 :            : 
      27                 :            : /** \file LinearPrism.hpp
      28                 :            :  *  \brief mapping function for linear prism
      29                 :            :  *  \author Jason Kraftcheck
      30                 :            :  */
      31                 :            : 
      32                 :            : #ifndef MSQ_LINEAR_PRISM_HPP
      33                 :            : #define MSQ_LINEAR_PRISM_HPP
      34                 :            : 
      35                 :            : #include "Mesquite.hpp"
      36                 :            : #include "MappingFunction.hpp"
      37                 :            : 
      38                 :            : namespace MBMesquite
      39                 :            : {
      40                 :            : 
      41                 :            : /**\brief Linear mapping function for a prism element
      42                 :            :  *
      43                 :            :  * \f$ \vec{x}(\vec{\xi})=\sum_{i=1}^6 N_i(\vec{\xi})\vec{x_i} \f$
      44                 :            :  * - \f$ N_1(\vec{\xi})=(1-\xi)(1-\eta-\zeta) \f$
      45                 :            :  * - \f$ N_2(\vec{\xi})=(1-\xi)(\eta)         \f$
      46                 :            :  * - \f$ N_3(\vec{\xi})=(1-\xi)(\zeta)        \f$
      47                 :            :  * - \f$ N_4(\vec{\xi})=(  \xi)(1-\eta-\zeta) \f$
      48                 :            :  * - \f$ N_5(\vec{\xi})=(  \xi)(\eta)         \f$
      49                 :            :  * - \f$ N_6(\vec{\xi})=(  \xi)(\zeta)        \f$
      50                 :            :  * \f$ \xi,\eta,\zeta \in [0,1] \f$
      51                 :            :  * \f$ \eta + \zeta \leq 1 \f$
      52                 :            :  *
      53                 :            :  * The logical center of the element is at:
      54                 :            :  * \f$ (\xi,\eta,\zeta) = (0,\frac{1}{3},\frac{1}{3}) \f$
      55                 :            :  */
      56         [ -  + ]:        392 : class MESQUITE_EXPORT LinearPrism : public MappingFunction3D
      57                 :            : {
      58                 :            :   public:
      59                 :            :     virtual EntityTopology element_topology() const;
      60                 :            : 
      61                 :            :     virtual int num_nodes() const;
      62                 :            : 
      63                 :            :     virtual void coefficients( Sample location, NodeSet nodeset, double* coeff_out, size_t* indices_out,
      64                 :            :                                size_t& num_coeff_out, MsqError& err ) const;
      65                 :            : 
      66                 :            :     virtual void derivatives( Sample location, NodeSet nodeset, size_t* vertex_indices_out,
      67                 :            :                               MsqVector< 3 >* d_coeff_d_xi_out, size_t& num_vtx, MsqError& err ) const;
      68                 :            : 
      69                 :            :     virtual void ideal( Sample location, MsqMatrix< 3, 3 >& jacobian_out, MsqError& err ) const;
      70                 :            : };
      71                 :            : 
      72                 :            : }  // namespace MBMesquite
      73                 :            : 
      74                 :            : #endif