From: Wolfgang Bangerth Date: Sat, 20 Nov 2010 22:34:28 +0000 (+0000) Subject: Make VectorTools::interpolate_boundary_values work for codim=1, at least for 2d trian... X-Git-Tag: v8.0.0~4830 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=c37662fb86f062eae25a82779ebd6068763ecc0e;p=dealii.git Make VectorTools::interpolate_boundary_values work for codim=1, at least for 2d triangulations. git-svn-id: https://svn.dealii.org/trunk@22829 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/include/deal.II/fe/fe_poly.templates.h b/deal.II/include/deal.II/fe/fe_poly.templates.h index f7642932dc..e397274e52 100644 --- a/deal.II/include/deal.II/fe/fe_poly.templates.h +++ b/deal.II/include/deal.II/fe/fe_poly.templates.h @@ -458,52 +458,6 @@ FE_Poly,1,2>::fill_fe_face_values -template <> -inline -void -FE_Poly,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,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,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 void diff --git a/tests/fail/interpolate_boundary_values_01.cc b/tests/fail/interpolate_boundary_values_01.cc index fbf6bef99c..e3f2eabfb2 100644 --- a/tests/fail/interpolate_boundary_values_01.cc +++ b/tests/fail/interpolate_boundary_values_01.cc @@ -55,6 +55,23 @@ void test(std::string filename) { for (std::map::const_iterator i = bv.begin(); i != bv.end(); ++i) deallog << i->first << ' ' << i->second << std::endl; + + for (typename DoFHandler::active_cell_iterator + cell = dof_handler.begin_active(); cell != dof_handler.end(); ++cell) + for (unsigned int f=0; f::faces_per_cell; ++f) + if (cell->at_boundary(f)) + for (unsigned int v=0; v::vertices_per_face; ++v) + for (unsigned int i=0; iface(f)->vertex_dof_index(v,i)) + != bv.end(), + ExcInternalError()); + Assert (bv[cell->face(f)->vertex_dof_index(v,i)] + == + Functions::SquareFunction() + .value(cell->face(f)->vertex(v),i), + ExcInternalError()); + } }