LinearQuad.hpp

Go to the documentation of this file.

    1 #ifndef LINEAR_QUAD_HPP
    2 #define LINEAR_QUAD_HPP
    3 /**\brief Shape function space for trilinear quadahedron, obtained by a pushforward of the canonical
    4  * linear (affine) functions. */
    5  
    6 #include "moab/ElemEvaluator.hpp"
    7 #include "moab/CN.hpp"
    8  
    9 namespace moab
   10 {
   11  
   12 class LinearQuad
   13 {
   14   public:
   15     /** \brief Forward-evaluation of field at parametric coordinates */
   16     static ErrorCode evalFcn( const double* params,
   17                               const double* field,
   18                               const int ndim,
   19                               const int num_tuples,
   20                               double* work,
   21                               double* result );
   22  
   23     /** \brief Reverse-evaluation of parametric coordinates at physical space position */
   24     static ErrorCode reverseEvalFcn( EvalFcn eval,
   25                                      JacobianFcn jacob,
   26                                      InsideFcn ins,
   27                                      const double* posn,
   28                                      const double* verts,
   29                                      const int nverts,
   30                                      const int ndim,
   31                                      const double iter_tol,
   32                                      const double inside_tol,
   33                                      double* work,
   34                                      double* params,
   35                                      int* is_inside );
   36  
   37     /** \brief Evaluate the normal at a specified facet*/
   38     static ErrorCode normalFcn( const int ientDim,
   39                                 const int facet,
   40                                 const int nverts,
   41                                 const double* verts,
   42                                 double normal[3] );
   43  
   44     /** \brief Evaluate the jacobian at a specified parametric position */
   45     static ErrorCode jacobianFcn( const double* params,
   46                                   const double* verts,
   47                                   const int nverts,
   48                                   const int ndim,
   49                                   double* work,
   50                                   double* result );
   51  
   52     /** \brief Forward-evaluation of field at parametric coordinates */
   53     static ErrorCode integrateFcn( const double* field,
   54                                    const double* verts,
   55                                    const int nverts,
   56                                    const int ndim,
   57                                    const int num_tuples,
   58                                    double* work,
   59                                    double* result );
   60  
   61     /** \brief Function that returns whether or not the parameters are inside the natural space of
   62      * the element */
   63     static int insideFcn( const double* params, const int ndim, const double tol );
   64  
   65     static EvalSet eval_set()
   66     {
   67         return EvalSet( evalFcn, reverseEvalFcn, normalFcn, jacobianFcn, integrateFcn, NULL, insideFcn );
   68     }
   69  
   70     static bool compatible( EntityType tp, int numv, EvalSet& eset )
   71     {
   72         if( tp == MBQUAD && numv == 4 )
   73         {
   74             eset = eval_set();
   75             return true;
   76         }
   77         else
   78             return false;
   79     }
   80  
   81   protected:
   82     /* Preimages of the vertices -- "canonical vertices" -- are known as "corners". */
   83     static const double corner[4][2];
   84     static const double gauss[1][2];
   85     static const unsigned int corner_count = 4;
   86     static const unsigned int gauss_count  = 1;
   87  
   88 };  // class LinearQuad
   89  
   90 }  // namespace moab
   91  
   92 #endif