Point FEInterface::map(unsigned int dim,
const FEType& fe_t,
const Elem* elem,
const Point& p)
{
fe_with_vec_switch(map(elem, p));
libmesh_error();
return Point();
}
Point FEInterface::map(unsigned int dim,
const FEType & fe_t,
const Elem * elem,
const Point & p)
{
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
if (is_InfFE_elem(elem->type()))
return ifem_map(dim, fe_t, elem, p);
#endif
fe_with_vec_switch(map(elem, p));
libmesh_error_msg("We'll never get here!");
return Point();
}
unsigned int FEInterface::n_dofs_per_elem(const unsigned int dim,
const FEType& fe_t,
const ElemType t)
{
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
if ( is_InfFE_elem(t) )
return ifem_n_dofs_per_elem(dim, fe_t, t);
#endif
const Order o = fe_t.order;
fe_with_vec_switch(n_dofs_per_elem(t, o));
libmesh_error();
return 0;
}
Point FEInterface::inverse_map (const unsigned int dim,
const FEType& fe_t,
const Elem* elem,
const Point& p,
const Real tolerance,
const bool secure)
{
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
if ( is_InfFE_elem(elem->type()) )
return ifem_inverse_map(dim, fe_t, elem, p,tolerance, secure);
#endif
fe_with_vec_switch(inverse_map(elem, p, tolerance, secure));
libmesh_error();
return Point();
}
unsigned int FEInterface::n_dofs_at_node(const unsigned int dim,
const FEType & fe_t,
const ElemType t,
const unsigned int n)
{
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
if ( is_InfFE_elem(t) )
return ifem_n_dofs_at_node(dim, fe_t, t, n);
#endif
const Order o = fe_t.order;
fe_with_vec_switch(n_dofs_at_node(t, o, n));
libmesh_error_msg("We'll never get here!");
return 0;
}
unsigned int FEInterface::n_shape_functions(const unsigned int dim,
const FEType& fe_t,
const ElemType t)
{
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
/*
* Since the FEType, stored in DofMap/(some System child), has to
* be the _same_ for InfFE and FE, we have to catch calls
* to infinite elements through the element type.
*/
if ( is_InfFE_elem(t) )
return ifem_n_shape_functions(dim, fe_t, t);
#endif
const Order o = fe_t.order;
fe_with_vec_switch(n_shape_functions(t, o));
libmesh_error();
return 0;
}
