From c37662fb86f062eae25a82779ebd6068763ecc0e Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Sat, 20 Nov 2010 22:34:28 +0000 Subject: [PATCH] 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 --- .../include/deal.II/fe/fe_poly.templates.h | 46 ------------------- tests/fail/interpolate_boundary_values_01.cc | 17 +++++++ 2 files changed, 17 insertions(+), 46 deletions(-) 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()); + } } -- 2.39.5