#include <QuadraticHex.hpp>

Static Public Member Functions
static ErrorCode	evalFcn (const double params, const double field, const int ndim, const int num_tuples, double work, double result)
	Forward-evaluation of field at parametric coordinates. More...

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)
	Reverse-evaluation of parametric coordinates at physical space position. More...

static ErrorCode	normalFcn (const int ientDim, const int facet, const int nverts, const double *verts, double normal[])
	Evaluate the normal at a specified facet. More...

static ErrorCode	jacobianFcn (const double params, const double verts, const int nverts, const int ndim, double work, double result)
	Evaluate the jacobian at a specified parametric position. More...

static ErrorCode	integrateFcn (const double field, const double verts, const int nverts, const int ndim, const int num_tuples, double work, double result)
	Forward-evaluation of field at parametric coordinates. More...

static int	insideFcn (const double *params, const int ndim, const double tol)
	Function that returns whether or not the parameters are inside the natural space of the element. More...

static EvalSet	eval_set ()

static bool	compatible (EntityType tp, int numv, EvalSet &eset)

Static Protected Member Functions
static double	SH (const int i, const double params)

static double	DSH (const int i, const double params)

Static Protected Attributes
static const int	corner [27][3]

static const double	gauss [8][2]

static const unsigned int	corner_count = 27

static const unsigned int	gauss_count = 8

Detailed Description

Definition at line 11 of file QuadraticHex.hpp.

Member Function Documentation

◆ compatible()

static bool moab::QuadraticHex::compatible	(	EntityType	tp,
		int	numv,
		EvalSet &	eset
	)

inlinestatic

Definition at line 69 of file QuadraticHex.hpp.

     {
         if( tp == MBHEX && numv == 27 )
         {
             eset = eval_set();
             return true;
         }
         else
             return false;
     }

References eval_set(), and MBHEX.

Referenced by moab::EvalSet::get_eval_set().

◆ DSH()

double moab::QuadraticHex::DSH	(	const int	i,
		const double	params
	)

staticprotected

Definition at line 45 of file QuadraticHex.cpp.

 {
     switch( i )
     {
         case -1:
             return params - 0.5;
         case 0:
             return -2 * params;
         case 1:
             return params + 0.5;
         default:
             return 0.;
     }
 }

References moab::params.

Referenced by jacobianFcn().

◆ eval_set()

static EvalSet moab::QuadraticHex::eval_set ( )

inlinestatic

Definition at line 64 of file QuadraticHex.hpp.

     {
         return EvalSet( evalFcn, reverseEvalFcn, normalFcn, jacobianFcn, integrateFcn, NULL, insideFcn );
     }

References evalFcn(), insideFcn(), integrateFcn(), jacobianFcn(), normalFcn(), and reverseEvalFcn().

Referenced by compatible().

◆ evalFcn()

ErrorCode moab::QuadraticHex::evalFcn	(	const double *	params,
		const double *	field,
		const int	ndim,
		const int	num_tuples,
		double *	work,
		double *	result
	)

static

Forward-evaluation of field at parametric coordinates.

Definition at line 60 of file QuadraticHex.cpp.

 {
     assert( params && field && num_tuples > 0 );
     std::fill( result, result + num_tuples, 0.0 );
     for( int i = 0; i < 27; i++ )
     {
         const double sh = SH( corner[i][0], params[0] ) * SH( corner[i][1], params[1] ) * SH( corner[i][2], params[2] );
         for( int j = 0; j < num_tuples; j++ )
             result[j] += sh * field[num_tuples * i + j];
     }
  
     return MB_SUCCESS;
 }

References corner, MB_SUCCESS, moab::params, and SH().

Referenced by eval_set().

◆ insideFcn()

int moab::QuadraticHex::insideFcn	(	const double *	params,
		const int	ndim,
		const double	tol
	)

static

Function that returns whether or not the parameters are inside the natural space of the element.

Definition at line 136 of file QuadraticHex.cpp.

 {
     return EvalSet::inside_function( params, ndim, tol );
 }

References moab::EvalSet::inside_function(), and moab::params.

Referenced by eval_set().

◆ integrateFcn()

ErrorCode moab::QuadraticHex::integrateFcn	(	const double *	field,
		const double *	verts,
		const int	nverts,
		const int	ndim,
		const int	num_tuples,
		double *	work,
		double *	result
	)

static

Forward-evaluation of field at parametric coordinates.

Definition at line 107 of file QuadraticHex.cpp.

 {
     return MB_NOT_IMPLEMENTED;
 }

References MB_NOT_IMPLEMENTED.

Referenced by eval_set().

◆ jacobianFcn()

ErrorCode moab::QuadraticHex::jacobianFcn	(	const double *	params,
		const double *	verts,
		const int	nverts,
		const int	ndim,
		double *	work,
		double *	result
	)

static

Evaluate the jacobian at a specified parametric position.

Definition at line 79 of file QuadraticHex.cpp.

 {
     assert( 27 == nverts && params && verts );
     if( 27 != nverts ) return MB_FAILURE;
     Matrix3* J = reinterpret_cast< Matrix3* >( result );
     for( int i = 0; i < 27; i++ )
     {
         const double sh[3]  = { SH( corner[i][0], params[0] ), SH( corner[i][1], params[1] ),
                                 SH( corner[i][2], params[2] ) };
         const double dsh[3] = { DSH( corner[i][0], params[0] ), DSH( corner[i][1], params[1] ),
                                 DSH( corner[i][2], params[2] ) };
  
         for( int j = 0; j < 3; j++ )
         {
             ( *J )( j, 0 ) += dsh[0] * sh[1] * sh[2] * verts[ndim * i + j];  // dxj/dr first column
             ( *J )( j, 1 ) += sh[0] * dsh[1] * sh[2] * verts[ndim * i + j];  // dxj/ds
             ( *J )( j, 2 ) += sh[0] * sh[1] * dsh[2] * verts[ndim * i + j];  // dxj/dt
         }
     }
  
     return MB_SUCCESS;
 }

References corner, DSH(), MB_SUCCESS, moab::params, and SH().

Referenced by eval_set().

◆ normalFcn()

ErrorCode moab::QuadraticHex::normalFcn	(	const int	ientDim,
		const int	facet,
		const int	nverts,
		const double *	verts,
		double	normal[]
	)

static

Evaluate the normal at a specified facet.

Definition at line 141 of file QuadraticHex.cpp.

 {
     return MB_NOT_IMPLEMENTED;
 }

References MB_NOT_IMPLEMENTED.

Referenced by eval_set().

◆ reverseEvalFcn()

ErrorCode moab::QuadraticHex::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

Reverse-evaluation of parametric coordinates at physical space position.

Definition at line 118 of file QuadraticHex.cpp.

 {
     assert( posn && verts );
     return EvalSet::evaluate_reverse( eval, jacob, ins, posn, verts, nverts, ndim, iter_tol, inside_tol, work, params,
                                       is_inside );
 }

References moab::EvalSet::evaluate_reverse(), and moab::params.

Referenced by eval_set().

◆ SH()

double moab::QuadraticHex::SH	(	const int	i,
		const double	params
	)

staticprotected

Definition at line 31 of file QuadraticHex.cpp.

 {
     switch( i )
     {
         case -1:
             return ( params * params - params ) / 2;
         case 0:
             return 1 - params * params;
         case 1:
             return ( params * params + params ) / 2;
         default:
             return 0.;
     }
 }

References moab::params.

Referenced by evalFcn(), and jacobianFcn().

Member Data Documentation

◆ corner

const int moab::QuadraticHex::corner

staticprotected

Initial value:

= { { -1, -1, -1 }, { 1, -1, -1 }, { 1, 1, -1 },  
                                          { -1, 1, -1 },                                
                                          { -1, -1, 1 },                 
                                          { 1, -1, 1 },                  
                                          { 1, 1, 1 },    { -1, 1, 1 },  
                                          { 0, -1, -1 },                 
                                          { 1, 0, -1 },                  
                                          { 0, 1, -1 },                  
                                          { -1, 0, -1 },                 
                                          { -1, -1, 0 },                 
                                          { 1, -1, 0 },                  
                                          { 1, 1, 0 },                   
                                          { -1, 1, 0 },                  
                                          { 0, -1, 1 },                  
                                          { 1, 0, 1 },                   
                                          { 0, 1, 1 },                   
                                          { -1, 0, 1 },                  
                                          { 0, -1, 0 },                  
                                          { 1, 0, 0 },                   
                                          { 0, 1, 0 },                   
                                          { -1, 0, 0 },   { 0, 0, -1 },  { 0, 0, 1 },  { 0, 0, 0 } }

Definition at line 85 of file QuadraticHex.hpp.

Referenced by evalFcn(), and jacobianFcn().

◆ corner_count

const unsigned int moab::QuadraticHex::corner_count = 27

staticprotected

Definition at line 87 of file QuadraticHex.hpp.

◆ gauss

const double moab::QuadraticHex::gauss[8][2]

staticprotected

Definition at line 86 of file QuadraticHex.hpp.

◆ gauss_count

const unsigned int moab::QuadraticHex::gauss_count = 8

staticprotected

Definition at line 88 of file QuadraticHex.hpp.

The documentation for this class was generated from the following files:

Static Public Member Functions

Static Protected Member Functions

Static Protected Attributes

Detailed Description

Member Function Documentation

◆ compatible()

◆ DSH()

◆ eval_set()

◆ evalFcn()

◆ insideFcn()

◆ integrateFcn()

◆ jacobianFcn()

◆ normalFcn()

◆ reverseEvalFcn()

◆ SH()

Member Data Documentation

◆ corner

◆ corner_count

◆ gauss

◆ gauss_count