Real FEInterface::shape(const unsigned int dim,
const FEType& fe_t,
const Elem* elem,
const unsigned int i,
const Point& p)
{
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
if ( is_InfFE_elem(elem->type()) )
return ifem_shape(dim, fe_t, elem, i, p);
#endif
const Order o = fe_t.order;
fe_switch(shape(elem,o,i,p));
libmesh_error();
return 0.;
}
Real FEInterface::shape(const unsigned int dim,
const FEType & fe_t,
const ElemType t,
const unsigned int i,
const Point & p)
{
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
if ( is_InfFE_elem(t) )
return ifem_shape(dim, fe_t, t, i, p);
#endif
const Order o = fe_t.order;
fe_switch(shape(t,o,i,p));
libmesh_error_msg("We'll never get here!");
return 0.;
}
void FEInterface::shape<Real>(const unsigned int dim,
const FEType& fe_t,
const ElemType t,
const unsigned int i,
const Point& p,
Real& phi)
{
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
if ( is_InfFE_elem(t) )
phi = ifem_shape(dim, fe_t, t, i, p);
#endif
const Order o = fe_t.order;
switch(dim)
{
case 0:
fe_scalar_vec_error_switch(0, shape(t,o,i,p), phi = , ; break;);
break;
case 1:
fe_scalar_vec_error_switch(1, shape(t,o,i,p), phi = , ; break;);