#include <deal.II/grid/grid_reordering.h>
#include <deal.II/grid/grid_tools.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
template <int dim>
-class MyBoundary : public Boundary<dim>
+class MyManifold : public Manifold<dim>
{
virtual Point<dim>
get_new_point_on_line (const typename Triangulation<dim>::line_iterator &line) const
template <int dim>
void check ()
{
- MyBoundary<dim> my_boundary;
+ MyManifold<dim> my_manifold;
// create two cubes
Triangulation<dim> coarse_grid (Triangulation<dim>::none, true);
Point<dim> p1 (-1,0,0), p2(1,1,1);
GridGenerator::subdivided_hyper_rectangle(coarse_grid, sub, p1, p2, true);
- // set bottom middle edge to use MyBoundary
+ // set bottom middle edge to use MyManifold
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
for (unsigned int e=0; e<GeometryInfo<dim-1>::faces_per_cell; ++e)
if (coarse_grid.begin_active()->face(f)->line(e)->center()[0] == 0)
if (coarse_grid.begin_active()->face(f)->line(e)->center()[1] == 0.5)
if (coarse_grid.begin_active()->face(f)->line(e)->center()[2] == 0)
- coarse_grid.begin_active()->face(f)->line(e)->set_boundary_id (99);
- coarse_grid.set_manifold (99, my_boundary);
+ coarse_grid.begin_active()->face(f)->line(e)->set_manifold_id (99);
+ coarse_grid.set_manifold (99, my_manifold);
// now try to refine this one
// cell. we should not get an exception, but keep it to make sure the