QuadraticHex.hpp

Go to the documentation of this file.

#ifndef QUADRATIC_HEX_HPP
#define QUADRATIC_HEX_HPP

#include "moab/ElemEvaluator.hpp"

namespace moab 
{
    
class QuadraticHex
{
public:
  static ErrorCode evalFcn(const double *params, const double *field, const int ndim, const int num_tuples, 
                           double *work, double *result);
        
  static ErrorCode reverseEvalFcn(EvalFcn eval, JacobianFcn jacob, InsideFcn ins, 
                                  const double *posn, const double *verts, const int nverts, const int ndim,
                                  const double iter_tol, const double inside_tol, double *work, 
                                  double *params, int *is_inside);
        
  static ErrorCode jacobianFcn(const double *params, const double *verts, const int nverts, const int ndim, 
                               double *work, double *result);
        
  static ErrorCode integrateFcn(const double *field, const double *verts, const int nverts, const int ndim, const int num_tuples, 
                                double *work, double *result);

  static int insideFcn(const double *params, const int ndim, const double tol);
  
  static EvalSet eval_set() 
      {
        return EvalSet(evalFcn, reverseEvalFcn, jacobianFcn, integrateFcn, NULL, insideFcn);
      }
      
  static bool compatible(EntityType tp, int numv, EvalSet &eset) 
      {
        if (tp == MBHEX && numv == 27) {
          eset = eval_set();
          return true;
        }
        else return false;
      }
  
protected:
  static double SH(const int i, const double params);
  static double DSH(const int i, const double params);
  
    /* Preimages of the vertices -- "canonical vertices" -- are known as "corners". */
  static const int corner[27][3];
  static const double gauss[8][2];// TODO fix me
  static const unsigned int corner_count = 27;
  static const unsigned int gauss_count  = 8; // TODO fix me
      
};// class QuadraticHex

} // moab namespace

#endif