Mesh Oriented datABase  (version 5.5.1)
An array-based unstructured mesh library
SpectralHex.hpp
Go to the documentation of this file.
1 #ifndef SPECTRAL_HEX_HPP
2 #define SPECTRAL_HEX_HPP
3 /**\brief Shape function space for spectral hexahedron
4  */
5 
6 #include "moab/ElemEvaluator.hpp"
7 #include "SpectralFuncs.hpp"
8 
9 namespace moab
10 {
11 
12 class SpectralHex
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 jacobian at a specified parametric position */
38  static ErrorCode jacobianFcn( const double* params,
39  const double* verts,
40  const int nverts,
41  const int ndim,
42  double* work,
43  double* result );
44 
45  /** \brief Forward-evaluation of field at parametric coordinates */
46  static ErrorCode integrateFcn( const double* field,
47  const double* verts,
48  const int nverts,
49  const int ndim,
50  const int num_tuples,
51  double* work,
52  double* result );
53 
54  /** \brief Initialize this EvalSet */
55  static ErrorCode initFcn( const double* verts, const int nverts, double*& work );
56 
57  /** \brief Function that returns whether or not the parameters are inside the natural space of
58  * the element */
59  static int insideFcn( const double* params, const int ndim, const double tol );
60 
61  static EvalSet eval_set()
62  {
63  return EvalSet( evalFcn, reverseEvalFcn, jacobianFcn, integrateFcn, initFcn );
64  }
65 
66  static bool compatible( EntityType tp, int numv, EvalSet& eset )
67  {
68  if( tp != MBHEX ) return false;
69  int i;
70  for( i = 3; i * i * i == numv || i * i * i > numv; i++ )
71  ;
72  if( i * i * i != numv ) return false;
73  eset = eval_set();
74  return true;
75  }
76 
77  protected:
78  static int _n;
79  static double* _z[3];
80  static lagrange_data _ld[3];
81  static opt_data_3 _data;
82  static double* _odwork; // work area
83  static bool init_;
84 }; // class SpectralHex
85 
86 } // namespace moab
87 
88 #endif