#include <ElemEvaluator.hpp>

Collaboration diagram for moab::EvalSet:

Public Member Functions
	EvalSet ()
	Bare constructor. More...

	EvalSet (EvalFcn eval, ReverseEvalFcn rev, NormalFcn normal, JacobianFcn jacob, IntegrateFcn integ, InitFcn initf, InsideFcn insidef)
	Constructor. More...

	EvalSet (EvalSet const &eval)
	Copy constructor. More...

EvalSet &	operator= (const EvalSet &eval)
	Operator=. More...

Static Public Member Functions
static ErrorCode	get_eval_set (Interface *mb, EntityHandle eh, EvalSet &eval_set)
	Given an entity handle, get an appropriate eval set, based on type & #vertices. More...

static ErrorCode	get_eval_set (EntityType tp, unsigned int num_vertices, EvalSet &eval_set)
	Given type & #vertices, get an appropriate eval set. More...

static ErrorCode	evaluate_reverse (EvalFcn eval, JacobianFcn jacob, InsideFcn inside_f, 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 *inside)
	Common function to do reverse evaluation based on evaluation and jacobian functions. More...

static int	inside_function (const double *params, const int ndims, const double tol)
	Common function that returns true if params is in [-1,1]^ndims. More...

Public Attributes
EvalFcn	evalFcn
	Forward-evaluation of field at parametric coordinates. More...

ReverseEvalFcn	reverseEvalFcn
	Reverse-evaluation of parametric coordinates at physical space position. More...

NormalFcn	normalFcn
	Evaluate the normal at a local facet (edge/face for 2D/3D) More...

JacobianFcn	jacobianFcn
	Evaluate the jacobian at a specified parametric position. More...

IntegrateFcn	integrateFcn
	Forward-evaluation of field at parametric coordinates. More...

InitFcn	initFcn
	Initialization function for an element. More...

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

Detailed Description

Definition at line 56 of file ElemEvaluator.hpp.

Constructor & Destructor Documentation

◆ EvalSet() [1/3]

moab::EvalSet::EvalSet ( )

inline

Bare constructor.

Definition at line 82 of file ElemEvaluator.hpp.

         : evalFcn( NULL ), reverseEvalFcn( NULL ), normalFcn( NULL ), jacobianFcn( NULL ), integrateFcn( NULL ),
           initFcn( NULL ), insideFcn( NULL )
     {
     }

◆ EvalSet() [2/3]

moab::EvalSet::EvalSet	(	EvalFcn	eval,
		ReverseEvalFcn	rev,
		NormalFcn	normal,
		JacobianFcn	jacob,
		IntegrateFcn	integ,
		InitFcn	initf,
		InsideFcn	insidef
	)

inline

Constructor.

Definition at line 89 of file ElemEvaluator.hpp.

         : evalFcn( eval ), reverseEvalFcn( rev ), normalFcn( normal ), jacobianFcn( jacob ), integrateFcn( integ ),
           initFcn( initf ), insideFcn( insidef )
     {
     }

◆ EvalSet() [3/3]

moab::EvalSet::EvalSet ( EvalSet const & eval )

inline

Copy constructor.

Definition at line 102 of file ElemEvaluator.hpp.

         : evalFcn( eval.evalFcn ), reverseEvalFcn( eval.reverseEvalFcn ), normalFcn( eval.normalFcn ),
           jacobianFcn( eval.jacobianFcn ), integrateFcn( eval.integrateFcn ), initFcn( eval.initFcn ),
           insideFcn( eval.insideFcn )
     {
     }

Member Function Documentation

◆ evaluate_reverse()

ErrorCode moab::EvalSet::evaluate_reverse	(	EvalFcn	eval,
		JacobianFcn	jacob,
		InsideFcn	inside_f,
		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 *	inside
	)

static

Common function to do reverse evaluation based on evaluation and jacobian functions.

Definition at line 19 of file ElemEvaluator.cpp.

 {
     // TODO: should differentiate between epsilons used for
     // Newton Raphson iteration, and epsilons used for curved boundary geometry errors
     // right now, fix the tolerance used for NR
     const double error_tol_sqr = iter_tol * iter_tol;
     CartVect* cvparams         = reinterpret_cast< CartVect* >( params );
     const CartVect* cvposn     = reinterpret_cast< const CartVect* >( posn );
  
     // initialize to center of element
     *cvparams = CartVect( -.4 );
  
     CartVect new_pos;
     // evaluate that first guess to get a new position
     ErrorCode rval = ( *eval )( cvparams->array(), verts, ndim,
                                 3,  // hardwire to num_tuples to 3 since the field is coords
                                 work, new_pos.array() );
     if( MB_SUCCESS != rval ) return rval;
  
     // residual is diff between old and new pos; need to minimize that
     CartVect res = new_pos - *cvposn;
     Matrix3 J;
     int dum, *tmp_inside = ( inside ? inside : &dum );
  
     int iters = 0;
     // while |res| larger than tol
     while( res % res > error_tol_sqr )
     {
         if( ++iters > 10 )
         {
             // if we haven't converged but we're outside, that's defined as success
             *tmp_inside = ( *inside_f )( params, ndim, inside_tol );
             if( !( *tmp_inside ) )
                 return MB_SUCCESS;
             else
                 return MB_FAILURE;
         }
  
         // get jacobian at current params
         rval       = ( *jacob )( cvparams->array(), verts, nverts, ndim, work, J.array() );
         double det = J.determinant();
         if( det < std::numeric_limits< double >::epsilon() )
         {
             *tmp_inside = ( *inside_f )( params, ndim, inside_tol );
             if( !( *tmp_inside ) )
                 return MB_SUCCESS;
             else
                 return MB_INDEX_OUT_OF_RANGE;
         }
  
         // new params tries to eliminate residual
         *cvparams -= J.inverse() * res;
  
         // get the new forward-evaluated position, and its difference from the target pt
         rval = ( *eval )( params, verts, ndim,
                           3,  // hardwire to num_tuples to 3 since the field is coords
                           work, new_pos.array() );
         if( MB_SUCCESS != rval ) return rval;
         res = new_pos - *cvposn;
     }
  
     if( inside ) *inside = ( *inside_f )( params, ndim, inside_tol );
  
     return MB_SUCCESS;
 }  // Map::evaluate_reverse()

References moab::CartVect::array(), moab::Matrix3::array(), moab::Matrix3::determinant(), moab::dum, ErrorCode, moab::Matrix3::inverse(), MB_INDEX_OUT_OF_RANGE, MB_SUCCESS, and moab::params.

Referenced by moab::LinearHex::reverseEvalFcn(), moab::LinearQuad::reverseEvalFcn(), and moab::QuadraticHex::reverseEvalFcn().

◆ get_eval_set() [1/2]

ErrorCode moab::EvalSet::get_eval_set	(	EntityType	tp,
		unsigned int	num_vertices,
		EvalSet &	eval_set
	)

static

Given type & #vertices, get an appropriate eval set.

Definition at line 107 of file ElemEvaluator.cpp.

 {
     switch( tp )
     {
         case MBEDGE:
             break;
         case MBTRI:
             if( LinearTri::compatible( tp, num_vertices, eval_set ) ) return MB_SUCCESS;
             break;
         case MBQUAD:
             if( LinearQuad::compatible( tp, num_vertices, eval_set ) ) return MB_SUCCESS;
             //            if (SpectralQuad::compatible(tp, num_vertices, eval_set)) return
             //            MB_SUCCESS;
             break;
         case MBTET:
             if( LinearTet::compatible( tp, num_vertices, eval_set ) ) return MB_SUCCESS;
             break;
         case MBHEX:
             if( LinearHex::compatible( tp, num_vertices, eval_set ) ) return MB_SUCCESS;
             if( QuadraticHex::compatible( tp, num_vertices, eval_set ) ) return MB_SUCCESS;
             //            if (SpectralHex::compatible(tp, num_vertices, eval_set)) return
             //            MB_SUCCESS;
             break;
         default:
             break;
     }
  
     return MB_NOT_IMPLEMENTED;
 }

References moab::LinearHex::compatible(), moab::LinearQuad::compatible(), moab::LinearTet::compatible(), moab::LinearTri::compatible(), moab::QuadraticHex::compatible(), MB_NOT_IMPLEMENTED, MB_SUCCESS, MBEDGE, MBHEX, MBQUAD, MBTET, and MBTRI.

◆ get_eval_set() [2/2]

ErrorCode moab::EvalSet::get_eval_set	(	Interface *	mb,
		EntityHandle	eh,
		EvalSet &	eval_set
	)

inlinestatic

Given an entity handle, get an appropriate eval set, based on type & #vertices.

Definition at line 147 of file ElemEvaluator.hpp.

 {
     int nv;
     EntityType tp = mb->type_from_handle( eh );
     const EntityHandle* connect;
     std::vector< EntityHandle > dum_vec;
     ErrorCode rval = mb->get_connectivity( eh, connect, nv, false, &dum_vec );
     if( MB_SUCCESS != rval ) return rval;
  
     return get_eval_set( tp, nv, eval_set );
 }

References ErrorCode, mb, and MB_SUCCESS.

Referenced by moab::ElemEvaluator::set_ent_handle(), and moab::ElemEvaluator::set_eval_set().

◆ inside_function()

int moab::EvalSet::inside_function	(	const double *	params,
		const int	ndims,
		const double	tol
	)

static

Common function that returns true if params is in [-1,1]^ndims.

Definition at line 96 of file ElemEvaluator.cpp.

 {
     if( params[0] >= -1 - tol && params[0] <= 1 + tol &&
         ( ndims < 2 || ( params[1] >= -1 - tol && params[1] <= 1 + tol ) ) &&
         ( ndims < 3 || ( params[2] >= -1 - tol && params[2] <= 1 + tol ) ) )
         return true;
     else
         return false;
 }

References moab::params.

Referenced by moab::LinearHex::insideFcn(), moab::LinearQuad::insideFcn(), and moab::QuadraticHex::insideFcn().

◆ operator=()

EvalSet& moab::EvalSet::operator= ( const EvalSet & eval )

inline

Operator=.

Definition at line 116 of file ElemEvaluator.hpp.

     {
         evalFcn        = eval.evalFcn;
         reverseEvalFcn = eval.reverseEvalFcn;
         normalFcn      = eval.normalFcn;
         jacobianFcn    = eval.jacobianFcn;
         integrateFcn   = eval.integrateFcn;
         initFcn        = eval.initFcn;
         insideFcn      = eval.insideFcn;
         return *this;
     }

References evalFcn, initFcn, insideFcn, integrateFcn, jacobianFcn, normalFcn, and reverseEvalFcn.

Member Data Documentation

◆ evalFcn

EvalFcn moab::EvalSet::evalFcn

Forward-evaluation of field at parametric coordinates.

Definition at line 60 of file ElemEvaluator.hpp.

Referenced by operator=(), and moab::ElemEvaluator::reverse_eval().

◆ initFcn

InitFcn moab::EvalSet::initFcn

Initialization function for an element.

Definition at line 75 of file ElemEvaluator.hpp.

Referenced by operator=(), and moab::ElemEvaluator::set_eval_set().

◆ insideFcn

InsideFcn moab::EvalSet::insideFcn

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

Definition at line 79 of file ElemEvaluator.hpp.

Referenced by operator=(), and moab::ElemEvaluator::reverse_eval().

◆ integrateFcn

IntegrateFcn moab::EvalSet::integrateFcn

Forward-evaluation of field at parametric coordinates.

Definition at line 72 of file ElemEvaluator.hpp.

Referenced by operator=().

◆ jacobianFcn

JacobianFcn moab::EvalSet::jacobianFcn

Evaluate the jacobian at a specified parametric position.

Definition at line 69 of file ElemEvaluator.hpp.

Referenced by operator=(), and moab::ElemEvaluator::reverse_eval().

◆ normalFcn

NormalFcn moab::EvalSet::normalFcn

Evaluate the normal at a local facet (edge/face for 2D/3D)

Definition at line 66 of file ElemEvaluator.hpp.

Referenced by operator=().

◆ reverseEvalFcn

ReverseEvalFcn moab::EvalSet::reverseEvalFcn

Reverse-evaluation of parametric coordinates at physical space position.

Definition at line 63 of file ElemEvaluator.hpp.

Referenced by operator=().

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

Public Member Functions

Static Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ EvalSet() [1/3]

◆ EvalSet() [2/3]

◆ EvalSet() [3/3]

Member Function Documentation

◆ evaluate_reverse()

◆ get_eval_set() [1/2]

◆ get_eval_set() [2/2]

◆ inside_function()

◆ operator=()

Member Data Documentation

◆ evalFcn

◆ initFcn

◆ insideFcn

◆ integrateFcn

◆ jacobianFcn

◆ normalFcn

◆ reverseEvalFcn