]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Make VectorTools::interpolate_boundary_values work for codim=1, at least for 2d trian...
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Sat, 20 Nov 2010 22:34:28 +0000 (22:34 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Sat, 20 Nov 2010 22:34:28 +0000 (22:34 +0000)
git-svn-id: https://svn.dealii.org/trunk@22829 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/include/deal.II/fe/fe_poly.templates.h
tests/fail/interpolate_boundary_values_01.cc

index f7642932dc7087067355347f417221407e585069..e397274e52b0d61279fa2a708dde7674bd969cee 100644 (file)
@@ -458,52 +458,6 @@ FE_Poly<TensorProductPolynomials<1>,1,2>::fill_fe_face_values
 
 
 
-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
index fbf6bef99c87eb470a44ce28c99ad7f3736f98e8..e3f2eabfb2540908eac93e5f6ac51fee85c01ccc 100644 (file)
@@ -55,6 +55,23 @@ void test(std::string filename) {
     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());
+             }
 }
 
 

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.