template <int dim>
-void FiniteElement<dim>::get_unit_support_points (vector<Point<dim> > &) const {
+void
+FiniteElement<dim>::get_unit_support_points (vector<Point<dim> > &) const
+{
Assert (false, ExcPureFunctionCalled());
};
template <int dim>
-void FiniteElement<dim>::get_support_points (const DoFHandler<dim>::cell_iterator &,
+void
+FiniteElement<dim>::get_support_points (const DoFHandler<dim>::cell_iterator &,
const Boundary<dim> &,
- vector<Point<dim> > &) const {
+ vector<Point<dim> > &) const
+{
Assert (false, ExcPureFunctionCalled());
};
+template <int dim>
+unsigned int
+FiniteElement<dim>::n_base_elements() const
+{
+ return 1;
+}
-
+template <int dim>
+const FiniteElement<dim>&
+FiniteElement<dim>::base_element(unsigned index) const
+{
+ Assert (index==0, ExcIndexRange(index,0,1));
+ return *this;
+}
/*------------------------------- Explicit Instantiations -------------*/
void FEValuesBase<dim>::get_function_values (const Vector<double> &fe_function,
vector< Vector<double> > &values) const
{
- Assert (fe->n_components == values.size(),
+ Assert (n_quadrature_points == values.size(),
ExcWrongNoOfComponents());
Assert (selected_dataset<shape_values.size(),
ExcInvalidIndex (selected_dataset, shape_values.size()));
for (unsigned i=0;i<values.size();++i)
- Assert (values[i].size() == n_quadrature_points,
+ Assert (values[i].size() == fe->n_components,
ExcWrongVectorSize(values.size(), n_quadrature_points));
// get function values of dofs