-template <>
-inline
-void
-FE_Poly<TensorProductPolynomials<2>,2,3>::fill_fe_face_values (const Mapping<2,3> &,
- const Triangulation<2,3>::cell_iterator &,
- const unsigned int,
- const Quadrature<1> &,
- Mapping<2,3>::InternalDataBase &,
- Mapping<2,3>::InternalDataBase &,
- FEValuesData<2,3> &) const
-{
- AssertThrow(false, ExcNotImplemented());
-}
-
-
-
-template <>
-inline
-void
-FE_Poly<PolynomialSpace<1>,1,2>::fill_fe_face_values (const Mapping<1,2> &,
- const Triangulation<1,2>::cell_iterator &,
- const unsigned int,
- const Quadrature<0> &,
- Mapping<1,2>::InternalDataBase &,
- Mapping<1,2>::InternalDataBase &,
- FEValuesData<1,2> &) const
-{
- AssertThrow(false, ExcNotImplemented());
-}
-
-
-
-template <>
-inline
-void
-FE_Poly<PolynomialSpace<2>,2,3>::fill_fe_face_values (const Mapping<2,3> &,
- const Triangulation<2,3>::cell_iterator &,
- const unsigned int,
- const Quadrature<1> &,
- Mapping<2,3>::InternalDataBase &,
- Mapping<2,3>::InternalDataBase &,
- FEValuesData<2,3> &) const
-{
- AssertThrow(false, ExcNotImplemented());
-}
-
template <class POLY, int dim, int spacedim>
void
for (std::map<unsigned int, double>::const_iterator i = bv.begin();
i != bv.end(); ++i)
deallog << i->first << ' ' << i->second << std::endl;
+
+ for (typename DoFHandler<dim,spacedim>::active_cell_iterator
+ cell = dof_handler.begin_active(); cell != dof_handler.end(); ++cell)
+ for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
+ if (cell->at_boundary(f))
+ for (unsigned int v=0; v<GeometryInfo<dim>::vertices_per_face; ++v)
+ for (unsigned int i=0; i<fe.dofs_per_vertex; ++i)
+ {
+ Assert (bv.find(cell->face(f)->vertex_dof_index(v,i))
+ != bv.end(),
+ ExcInternalError());
+ Assert (bv[cell->face(f)->vertex_dof_index(v,i)]
+ ==
+ Functions::SquareFunction<spacedim>()
+ .value(cell->face(f)->vertex(v),i),
+ ExcInternalError());
+ }
}