const InVector &fe_function,
const Point<spacedim> &point)
{
-#ifdef DEBUG
- const FiniteElement<dim> &fe = dof.get_fe();
- Assert(fe.n_components() == 1,
+ Assert(dof.get_fe().n_components() == 1,
ExcMessage ("Finite element is not scalar as is necessary for this function"));
-#endif
Vector<double> value(1);
point_value(mapping, dof, fe_function, point, value);
const InVector &fe_function,
const Point<spacedim> &point)
{
-#ifdef DEBUG
- const hp::FECollection<dim, spacedim> &fe = dof.get_fe();
- Assert(fe.n_components() == 1,
+ Assert(dof.get_fe().n_components() == 1,
ExcMessage ("Finite element is not scalar as is necessary for this function"));
-#endif
Vector<double> value(1);
point_value(mapping, dof, fe_function, point, value);
const InVector &fe_function,
const Point<spacedim> &point)
{
-#ifdef DEBUG
- const FiniteElement<dim> &fe = dof.get_fe();
- Assert(fe.n_components() == 1,
+ Assert(dof.get_fe().n_components() == 1,
ExcMessage ("Finite element is not scalar as is necessary for this function"));
-#endif
std::vector<Tensor<1, dim> > gradient(1);
point_gradient (mapping, dof, fe_function, point, gradient);
const InVector &fe_function,
const Point<spacedim> &point)
{
-#ifdef DEBUG
- const hp::FECollection<dim, spacedim> &fe = dof.get_fe();
- Assert(fe.n_components() == 1,
+ Assert(dof.get_fe().n_components() == 1,
ExcMessage ("Finite element is not scalar as is necessary for this function"));
-#endif
std::vector<Tensor<1, dim> > gradient(1);
point_gradient (mapping, dof, fe_function, point, gradient);