#include <deal.II/base/logstream.h>
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/grid_generator.h>
p1[0] = -1;
p2[0] = 1;
static const ConeBoundary<dim> boundary (1, 0.5, p1, p2);
- triangulation.set_boundary (0, boundary);
+ triangulation.set_manifold (0, boundary);
triangulation.refine_global (2);
-// check the cells generated by the CylinderBoundary the default
+// check the cells generated by the CylinderManifold the default
// axiparallel cylinder
#include "../tests.h"
#include <deal.II/base/logstream.h>
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/grid_generator.h>
{
Triangulation<dim> triangulation;
GridGenerator::cylinder (triangulation);
- static const CylinderBoundary<dim> boundary (1,0);
- triangulation.set_boundary (0, boundary);
+ static const CylinderManifold<dim> boundary (0);
+ triangulation.set_manifold (0, boundary);
triangulation.refine_global (2);
for (typename Triangulation<dim>::active_cell_iterator
#include <deal.II/base/logstream.h>
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/grid_generator.h>
GridTools::transform ((Point<dim> ( *)(const Point<dim> &))&rotate_to_y<dim>, triangulation);
- static const CylinderBoundary<dim> boundary (1,1);
- triangulation.set_boundary (0, boundary);
+ static const CylinderManifold<dim> boundary (1);
+ triangulation.set_manifold (0, boundary);
triangulation.refine_global (2);
for (typename Triangulation<dim>::active_cell_iterator
#include <deal.II/base/logstream.h>
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/grid_generator.h>
GridTools::transform ((Point<dim> ( *)(const Point<dim> &))&rotate_to_y<dim>, triangulation);
- static const CylinderBoundary<dim> boundary (1,2);
- triangulation.set_boundary (0, boundary);
+ static const CylinderManifold<dim> boundary (2);
+ triangulation.set_manifold (0, boundary);
triangulation.refine_global (2);
for (typename Triangulation<dim>::active_cell_iterator
#include <deal.II/base/logstream.h>
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/grid_generator.h>
GridTools::transform ((Point<dim> ( *)(const Point<dim> &))&rotate_to_y<dim>, triangulation);
- static const CylinderBoundary<dim> boundary (1,
- Point<dim>(std::cos(xy_angle), std::sin(xy_angle), 0),
+ static const CylinderManifold<dim> boundary (Point<dim>(std::cos(xy_angle), std::sin(xy_angle), 0),
Point<dim>());
- triangulation.set_boundary (0, boundary);
+ triangulation.set_manifold (0, boundary);
triangulation.refine_global (2);
for (typename Triangulation<dim>::active_cell_iterator
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/grid_reordering.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>
// set bottom face to use MyBoundary
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (coarse_grid.begin_active()->face(f)->center()[dim-1] == -1)
- coarse_grid.begin_active()->face(f)->set_boundary_indicator (1);
- coarse_grid.set_boundary (1, my_boundary);
+ coarse_grid.begin_active()->face(f)->set_manifold_indicator (1);
+ coarse_grid.set_manifold (1, my_boundary);
// now try to refine this one
// cell. we should get an exception
#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>
// set bottom face to use MyBoundary
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (coarse_grid.begin_active()->face(f)->center()[dim-1] == -1)
- coarse_grid.begin_active()->face(f)->set_boundary_indicator (1);
- coarse_grid.set_boundary (1, my_boundary);
+ coarse_grid.begin_active()->face(f)->set_manifold_indicator (1);
+ coarse_grid.set_manifold (1, my_boundary);
// now try to refine this one
// cell. we should get an exception
// set bottom face to use MyBoundary
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (coarse_grid.begin_active()->face(f)->center()[dim-1] == -1)
- coarse_grid.begin_active()->face(f)->set_boundary_indicator (1);
- coarse_grid.set_boundary (1, my_boundary);
+ coarse_grid.begin_active()->face(f)->set_manifold_indicator (1);
+ coarse_grid.set_manifold (1, my_boundary);
// now try to refine this one
// cell. we should get an exception
#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>
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_indicator (99);
- coarse_grid.set_boundary (99, my_boundary);
+ coarse_grid.begin_active()->face(f)->line(e)->set_manifold_indicator (99);
+ coarse_grid.set_manifold (99, my_boundary);
// now try to refine this one
// cell. we should get an exception
#include <deal.II/fe/fe_values.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
void test()
// shell is so narrow that MappingQ(2) distorts the cell
GridGenerator::quarter_hyper_shell (tria, Point<dim>(), 0.95, 1, 1);
- static HyperShellBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
FE_Nothing<dim> dummy;
MappingQ<dim> mapping(2);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/numerics/fe_field_function.h>
#include <deal.II/numerics/vector_tools.h>
template<int dim>
void test()
{
- const HyperBallBoundary<dim> boundary_description;
+ const PolarManifold<dim> boundary_description;
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- triangulation.set_boundary (0, boundary_description);
+ triangulation.set_manifold (0, boundary_description);
triangulation.refine_global (1);
FE_Q<dim> fe(2);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/numerics/fe_field_function.h>
#include <deal.II/numerics/vector_tools.h>
template<int dim>
void test()
{
- const HyperBallBoundary<dim> boundary_description;
+ const PolarManifold<dim> boundary_description;
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- triangulation.set_boundary (0, boundary_description);
+ triangulation.set_manifold (0, boundary_description);
triangulation.refine_global (1);
FESystem<dim> fe(FE_Q<dim> (2),2);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/numerics/fe_field_function.h>
#include <deal.II/numerics/vector_tools.h>
template<int dim>
void test()
{
- const HyperBallBoundary<dim> boundary_description;
+ const PolarManifold<dim> boundary_description;
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- triangulation.set_boundary (0, boundary_description);
+ triangulation.set_manifold (0, boundary_description);
triangulation.refine_global (1);
FE_Q<dim> fe(2);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/numerics/fe_field_function.h>
#include <deal.II/numerics/vector_tools.h>
template<int dim>
void test()
{
- const HyperBallBoundary<dim> boundary_description;
+ const PolarManifold<dim> boundary_description;
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- triangulation.set_boundary (0, boundary_description);
+ triangulation.set_manifold (0, boundary_description);
triangulation.refine_global (1);
FESystem<dim> fe(FE_Q<dim> (2),2);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/numerics/fe_field_function.h>
#include <deal.II/numerics/vector_tools.h>
template<int dim>
void test()
{
- const HyperBallBoundary<dim> boundary_description;
+ const PolarManifold<dim> boundary_description;
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- triangulation.set_boundary (0, boundary_description);
+ triangulation.set_manifold (0, boundary_description);
triangulation.refine_global (1);
FESystem<dim> fe(FE_Q<dim> (2),2);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/numerics/fe_field_function.h>
#include <deal.II/numerics/vector_tools.h>
template<int dim>
void test()
{
- const HyperBallBoundary<dim> boundary_description;
+ const PolarManifold<dim> boundary_description;
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- triangulation.set_boundary (0, boundary_description);
+ triangulation.set_manifold (0, boundary_description);
triangulation.refine_global (1);
FESystem<dim> fe(FE_Q<dim> (2),2);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/numerics/fe_field_function.h>
#include <deal.II/numerics/vector_tools.h>
template<int dim>
void test()
{
- const HyperBallBoundary<dim> boundary_description;
+ const PolarManifold<dim> boundary_description;
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- triangulation.set_boundary (0, boundary_description);
+ triangulation.set_manifold (0, boundary_description);
triangulation.refine_global (1);
FESystem<dim> fe(FE_Q<dim> (4),2);
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_accessor.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/fe/fe_values.h>
{
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
triangulation.refine_global (1);
}
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/hp/dof_handler.h>
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
MappingQ<dim> mapping(degree);
const FE_Q<dim> dummy_fe (1);
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
const MappingQ<dim> mapping (degree);
const FE_Q<dim> fe (1);
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/hp/dof_handler.h>
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
for (unsigned int refinement=0; refinement<2;
++refinement, triangulation.refine_global(1))
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
MappingQ<dim> mapping(degree);
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
const MappingQ<dim> mapping (degree);
const FE_Q<dim> fe (1);
#include <deal.II/lac/precondition.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/dofs/dof_handler.h>
void LaplaceProblem<dim>::run ()
{
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
for (unsigned int cycle=0; cycle<6; ++cycle, triangulation.refine_global(1))
{
#include <deal.II/grid/grid_in.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <fstream>
#include <sstream>
grid_in.read_ucd (input_file);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
}
deallog << " Number of active cells: "
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_accessor.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/hp/fe_values.h>
{
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
triangulation.refine_global (1);
}
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/grid/grid_tools.h>
const int dim = 3;
deallog << "Testing hyper_cube in dim: " << dim << "..."<< endl;
- const HyperBallBoundary<dim> boundary_description;
+ const PolarManifold<dim> boundary_description;
Triangulation<dim> volume_mesh;
GridGenerator::hyper_ball(volume_mesh);
- volume_mesh.set_boundary (0, boundary_description);
+ volume_mesh.set_manifold (0, boundary_description);
// exclude one of the 6 faces
// from the surface mesh
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (cell->at_boundary(f))
{
- cell->face(f)->set_boundary_indicator(1);
+ cell->face(f)->set_manifold_indicator(1);
done = true;
break;
}
const HyperBallBoundary<dim-1,dim> surface_description;
Triangulation<dim-1,dim> boundary_mesh;
- boundary_mesh.set_boundary (0, surface_description);
+ boundary_mesh.set_manifold (0, surface_description);
// now extract a mesh of the 5
// surface faces
cell != boundary_mesh.end(); ++cell)
for (unsigned int f=0; f<GeometryInfo<dim-1>::faces_per_cell; ++f)
if (cell->at_boundary(f))
- cell->face(f)->set_boundary_indicator(1);
+ cell->face(f)->set_manifold_indicator(1);
boundary_mesh.refine_global (2);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/grid/grid_tools.h>
map< Triangulation<dim-1,dim>::cell_iterator,
Triangulation<dim,dim>::face_iterator>
surface_to_volume_mapping;
- const HyperBallBoundary<dim> boundary_description;
+ const PolarManifold<dim> boundary_description;
Triangulation<dim> volume_mesh;
GridGenerator::hyper_ball(volume_mesh);
- volume_mesh.set_boundary (0, boundary_description);
+ volume_mesh.set_manifold (0, boundary_description);
volume_mesh.refine_global (1);
const HyperBallBoundary<dim-1,dim> surface_description;
Triangulation<dim-1,dim> boundary_mesh;
- boundary_mesh.set_boundary (0, surface_description);
+ boundary_mesh.set_manifold (0, surface_description);
surface_to_volume_mapping
= GridTools::extract_boundary_mesh (volume_mesh, boundary_mesh);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/grid/grid_tools.h>
map< Triangulation<dim-1,dim>::cell_iterator,
Triangulation<dim,dim>::face_iterator>
surface_to_volume_mapping;
- const HyperBallBoundary<dim> boundary_description;
+ const PolarManifold<dim> boundary_description;
Triangulation<dim> volume_mesh;
GridGenerator::hyper_ball(volume_mesh);
for (Triangulation<dim>::active_cell_iterator
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (cell->at_boundary(f))
cell->face(f)->set_all_boundary_indicators (1);
- volume_mesh.set_boundary (1, boundary_description);
+ volume_mesh.set_manifold (1, boundary_description);
volume_mesh.refine_global (1);
const HyperBallBoundary<dim-1,dim> surface_description;
Triangulation<dim-1,dim> boundary_mesh;
- boundary_mesh.set_boundary (1, surface_description);
+ boundary_mesh.set_manifold (1, surface_description);
surface_to_volume_mapping
= GridTools::extract_boundary_mesh (volume_mesh, boundary_mesh);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/grid/grid_tools.h>
map< Triangulation<dim-1,dim>::cell_iterator,
Triangulation<dim,dim>::face_iterator>
surface_to_volume_mapping;
- const HyperBallBoundary<dim> boundary_description;
+ const PolarManifold<dim> boundary_description;
Triangulation<dim> volume_mesh;
GridGenerator::hyper_ball(volume_mesh);
for (Triangulation<dim>::active_cell_iterator
if (cell->at_boundary(f) && (cell->center()[dim-1]>0))
cell->face(f)->set_all_boundary_indicators (1);
- volume_mesh.set_boundary (1, boundary_description);
- volume_mesh.set_boundary (0, boundary_description);
+ volume_mesh.set_manifold (1, boundary_description);
+ volume_mesh.set_manifold (0, boundary_description);
volume_mesh.refine_global (1);
const HyperBallBoundary<dim-1,dim> surface_description;
Triangulation<dim-1,dim> boundary_mesh;
- boundary_mesh.set_boundary (1, surface_description);
- boundary_mesh.set_boundary (0, surface_description);
+ boundary_mesh.set_manifold (1, surface_description);
+ boundary_mesh.set_manifold (0, surface_description);
set<types::boundary_id> boundary_ids;
boundary_ids.insert(1);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_in.h>
#include <deal.II/grid/grid_out.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <fstream>
#include <string>
{
HyperBallBoundary<dim, spacedim> boundary;
Triangulation<dim, spacedim> tria;
- tria.set_boundary(1, boundary);
+ tria.set_manifold(1, boundary);
GridIn<dim, spacedim> gi;
gi.attach_triangulation (tria);
std::ifstream in (filename.c_str());
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/grid/grid_tools.h>
#include <deal.II/fe/fe_values.h>
Triangulation<dim> volume_mesh;
GridGenerator::hyper_ball(volume_mesh);
- const HyperBallBoundary<dim-1,dim> surface_description;
+ const PolarManifold<dim> surface_description;
Triangulation<dim-1,dim> boundary_mesh;
- boundary_mesh.set_boundary (0, surface_description);
+ boundary_mesh.set_manifold (0, surface_description);
GridTools::extract_boundary_mesh (volume_mesh, boundary_mesh);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/grid/grid_tools.h>
#include <deal.II/fe/fe_values.h>
const HyperBallBoundary<dim-1,dim> surface_description;
Triangulation<dim-1,dim> boundary_mesh;
- boundary_mesh.set_boundary (0, surface_description);
+ boundary_mesh.set_manifold (0, surface_description);
GridTools::extract_boundary_mesh (volume_mesh, boundary_mesh);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_out.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <fstream>
#include <string>
TorusBoundary<dim, spacedim> boundary (1.5, .5);
Triangulation<dim, spacedim> tria;
- tria.set_boundary (0, boundary);
+ tria.set_manifold (0, boundary);
GridGenerator::torus (tria, 1.5, .5);
tria.refine_global(2);
#include <deal.II/lac/vector.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_accessor.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/dofs/dof_handler.h>
}
else
GridGenerator::hyper_cube(tr, -1./std::sqrt(static_cast<double>(dim)),1./std::sqrt(static_cast<double>(dim)));
- static const HyperBallBoundary<dim> boundary;
+ static const PolarManifold<dim> boundary;
if (dim != 1)
{
- tr.set_boundary (0, boundary);
+ tr.set_manifold (0, boundary);
}
tr.refine_global (1);
tr.begin_active()->set_refine_flag ();
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/fe/fe_q.h>
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/lac/sparse_matrix.h>
template <int dim>
class Ball
:
- public StraightBoundary<dim>
+ public FlatManifold<dim>
{
public:
virtual Point<dim>
get_new_point_on_line (const typename Triangulation<dim>::line_iterator &line) const
{
- Point<dim> middle = StraightBoundary<dim>::get_new_point_on_line(line);
+ Point<dim> middle = FlatManifold<dim>::get_new_point_on_line(line);
for (int i=0; i<dim; ++i)
middle(i) -= .5;
virtual Point<dim>
get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const
{
- Point<dim> middle = StraightBoundary<dim>::get_new_point_on_quad(quad);
+ Point<dim> middle = FlatManifold<dim>::get_new_point_on_quad(quad);
for (int i=0; i<dim; ++i)
middle(i) -= .5;
template <int dim>
class CurvedLine
:
- public StraightBoundary<dim>
+ public FlatManifold<dim>
{
public:
virtual Point<dim>
Point<dim>
CurvedLine<dim>::get_new_point_on_line (const typename Triangulation<dim>::line_iterator &line) const
{
- Point<dim> middle = StraightBoundary<dim>::get_new_point_on_line (line);
+ Point<dim> middle = FlatManifold<dim>::get_new_point_on_line (line);
// if the line is at the top of bottom
// face: do a special treatment on
Point<dim>
CurvedLine<dim>::get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const
{
- Point<dim> middle = StraightBoundary<dim>::get_new_point_on_quad (quad);
+ Point<dim> middle = FlatManifold<dim>::get_new_point_on_quad (quad);
// if the face is at the top of bottom
// face: do not move the midpoint in
{
if (dim==3)
{
- tria->begin_active()->face(4)->set_boundary_indicator(1);
- tria->begin_active()->face(5)->set_boundary_indicator(1);
+ tria->begin_active()->face(4)->set_manifold_indicator(1);
+ tria->begin_active()->face(5)->set_manifold_indicator(1);
};
// set the boundary function
- tria->set_boundary(1, (test_case==2)
+ tria->set_manifold(1, (test_case==2)
? ((Boundary<dim> &)ball) : ((Boundary<dim> &)curved_line));
// refine once
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_Q<dim>(2), 2,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_Q<dim>(2), 2,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_raviart_thomas.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_RaviartThomas<dim>(1), 1,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_RaviartThomas<dim>(1), 1,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_Q<dim>(2), 2,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_Q<dim>(2), 2,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_Q<dim>(2), 2,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_RaviartThomas<dim>(1), 1,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_RaviartThomas<dim>(1), 1,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_RaviartThomas<dim>(1), 1,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FE_RaviartThomas<dim> fe(1);
test(tr, fe);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FE_Nedelec<dim> fe(1);
test(tr, fe);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe(FE_RaviartThomas<dim>(1),1);
test(tr, fe);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_RaviartThomas<dim>(1), 1,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_Q<dim>(2), 2,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_Q<dim>(2), 2,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_Q<dim>(2), 2,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_RaviartThomas<dim>(1), 1,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_RaviartThomas<dim>(1), 1,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_RaviartThomas<dim>(1), 1,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_Q<dim>(2), 2,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_RaviartThomas<dim>(1), 1,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_Q<dim>(2), 2,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1), 1,
FE_RaviartThomas<dim>(1), 1,
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_system.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1),
SymmetricTensor<2,dim>::n_independent_components);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_system.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1),
Tensor<2,dim>::n_independent_components);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_system.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1),
SymmetricTensor<2,dim>::n_independent_components);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_system.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
FESystem<dim> fe (FE_Q<dim>(1),
Tensor<2,dim>::n_independent_components);
// refinement
#include "../tests.h"
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/tria.h>
Point<dim> center;
Triangulation<dim> tria (Triangulation<dim>::none, true);
GridGenerator::hyper_shell (tria, center, r1, r2, n, true);
- static const HyperShellBoundary<dim> boundary(center);
- tria.set_boundary(0, boundary);
+ static const PolarManifold<dim> boundary(center);
+ tria.set_manifold(0, boundary);
for (unsigned int i=0; i<2; ++i)
{
// this series test the other possibilities for the number of cells
#include "../tests.h"
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/tria.h>
Point<dim> center;
Triangulation<dim> tria (Triangulation<dim>::none);
GridGenerator::hyper_shell (tria, center, r1, r2, n);
- static const HyperShellBoundary<dim> boundary(center);
- tria.set_boundary(0, boundary);
+ static const PolarManifold<dim> boundary(center);
+ tria.set_manifold(0, boundary);
tria.refine_global(1);
GridOut grid_out;
// this series test the other possibilities for the number of cells
#include "../tests.h"
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/tria.h>
Point<dim> center;
Triangulation<dim> tria (Triangulation<dim>::none);
GridGenerator::hyper_shell (tria, center, r1, r2, n);
- static const HyperShellBoundary<dim> boundary(center);
- tria.set_boundary(0, boundary);
+ static const PolarManifold<dim> boundary(center);
+ tria.set_manifold(0, boundary);
tria.refine_global(1);
GridOut grid_out;
// this series test the other possibilities for the number of cells
#include "../tests.h"
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/tria.h>
Point<dim> center;
Triangulation<dim> tria (Triangulation<dim>::none);
GridGenerator::hyper_shell (tria, center, r1, r2, n);
- static const HyperShellBoundary<dim> boundary(center);
- tria.set_boundary(0, boundary);
+ static const PolarManifold<dim> boundary(center);
+ tria.set_manifold(0, boundary);
GridOut grid_out;
grid_out.write_gnuplot (tria, logfile);
#include "../tests.h"
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/base/logstream.h>
void test ()
{
Triangulation<dim> tria;
- static const HyperBallBoundary<dim> x;
+ static const PolarManifold<dim> x;
if (dim == 2)
{
- tria.set_boundary (0, x);
+ tria.set_manifold (0, x);
GridGenerator::hyper_ball (tria);
}
else
#include <deal.II/base/logstream.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/grid/grid_in.h>
Triangulation<dim> tr;
GridGenerator::hyper_cube_with_cylindrical_hole(tr, .3, .4, 1, 1, false);
- CylinderBoundary<dim> boundary(.3, 2);
- tr.set_boundary(1, boundary);
+ CylinderManifold<dim> boundary(.3, 2);
+ tr.set_manifold(1, boundary);
{
std::ofstream grid_file("coarse_grid.inp");
go.write_ucd(tr, grid_file);
deallog << "Writing refined from constructor" << std::endl;
go.write_ucd(tr, out);
- tr.set_boundary(1);
+ tr.set_manifold(1);
tr.clear();
GridIn<dim> gi;
grid_file.close();
}
- tr.set_boundary(1, boundary);
+ tr.set_manifold(1, boundary);
tr.refine_global(1);
deallog << "Writing refined from file" << std::endl;
go.write_ucd(tr, out);
grid_file.close();
}
- tr.set_boundary(1);
+ tr.set_manifold(1);
tr.clear();
}
#include "../tests.h"
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/tria.h>
template <int dim>
class Ball :
- public StraightBoundary<dim>
+ public FlatManifold<dim>
{
public:
virtual Point<dim>
get_new_point_on_line (const typename Triangulation<dim>::line_iterator &line) const
{
- Point<dim> middle = StraightBoundary<dim>::get_new_point_on_line(line);
+ Point<dim> middle = FlatManifold<dim>::get_new_point_on_line(line);
for (int i=0; i<dim; ++i)
middle(i) -= .5;
virtual Point<dim>
get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const
{
- Point<dim> middle = StraightBoundary<dim>::get_new_point_on_quad(quad);
+ Point<dim> middle = FlatManifold<dim>::get_new_point_on_quad(quad);
for (int i=0; i<dim; ++i)
middle(i) -= .5;
template <int dim>
class CurvedLine :
- public StraightBoundary<dim>
+ public FlatManifold<dim>
{
public:
virtual Point<dim>
Point<dim>
CurvedLine<dim>::get_new_point_on_line (const typename Triangulation<dim>::line_iterator &line) const
{
- Point<dim> middle = StraightBoundary<dim>::get_new_point_on_line (line);
+ Point<dim> middle = FlatManifold<dim>::get_new_point_on_line (line);
// if the line is at the top of bottom
// face: do a special treatment on
Point<dim>
CurvedLine<dim>::get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const
{
- Point<dim> middle = StraightBoundary<dim>::get_new_point_on_quad (quad);
+ Point<dim> middle = FlatManifold<dim>::get_new_point_on_quad (quad);
// if the face is at the top of bottom
// face: do not move the midpoint in
{
if (dim==3)
{
- tria.begin_active()->face(4)->set_boundary_indicator(1);
- tria.begin_active()->face(5)->set_boundary_indicator(1);
+ tria.begin_active()->face(4)->set_manifold_indicator(1);
+ tria.begin_active()->face(5)->set_manifold_indicator(1);
};
Ball<dim> ball;
CurvedLine<dim> curved_line;
if (test_case==2)
- tria.set_boundary (1, ball);
+ tria.set_manifold (1, ball);
else
- tria.set_boundary (1, curved_line);
+ tria.set_manifold (1, curved_line);
// refine once
tria.begin_active()->set_refine_flag();
tria.execute_coarsening_and_refinement();
};
- tria.set_boundary (1);
+ tria.set_manifold (1);
break;
}
}
#include "../tests.h"
#include <deal.II/base/logstream.h>
#include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_tools.h>
#include <deal.II/grid/grid_out.h>
Triangulation<dim> tria;
GridGenerator::hyper_ball(tria, Point<dim>(), 1);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
for (unsigned int i=0; i<4; ++i)
{
#include <deal.II/base/logstream.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/tria_accessor.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/fe/mapping_q.h>
{
const unsigned int dim=2;
Point<dim> origin;
- HyperShellBoundary<dim> boundary(origin);
+ PolarManifold<dim> boundary(origin);
MappingQ<dim> mapping(2);
Triangulation<dim> tria;
- tria.set_boundary(0, boundary);
+ tria.set_manifold(0, boundary);
const double inner_radius=1.;
const double outer_radius=5.;
GridGenerator::hyper_shell(tria, origin, inner_radius, outer_radius, 8);
// circle.
new_points.insert(std::pair<types::global_dof_index, Point<dim> > (
face->vertex_index(vertex_no), n_radius/inner_radius*v+n_center));
- face->set_boundary_indicator(1);
+ face->set_manifold_indicator(1);
}
else
Assert(false, ExcInternalError());
}
GridGenerator::laplace_transformation (tria, new_points);
- HyperBallBoundary<dim> inner_ball(n_center, n_radius);
- tria.set_boundary(1, inner_ball);
+ PolarManifold<dim> inner_ball(n_center, n_radius);
+ tria.set_manifold(1, inner_ball);
GridOut grid_out;
std::ofstream eps_stream2("output");
#include <deal.II/base/logstream.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/tria_accessor.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_out.h>
{
const unsigned int dim=3;
Point<dim> origin;
- CylinderBoundary<dim> boundary;
+ CylinderManifold<dim> boundary;
Triangulation<dim> tria;
- tria.set_boundary(0, boundary);
+ tria.set_manifold(0, boundary);
GridGenerator::cylinder(tria, 1, .7);
tria.refine_global(2);
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_accessor.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/fe/fe_values.h>
{
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
triangulation.refine_global (1);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/constraint_matrix.h>
#include <deal.II/fe/fe_q.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
tr.refine_global(2);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/constraint_matrix.h>
#include <deal.II/fe/fe_q.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
tr.refine_global(2);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/dofs/dof_renumbering.h>
#include <deal.II/lac/constraint_matrix.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
tr.refine_global(1);
#include <grid/tria.h>
#include <grid/grid_generator.h>
-#include <grid/tria_boundary_lib.h>
+#include <grid/manifold_lib.h>
#include <grid/grid_out.h>
#include <dofs/dof_handler.h>
#include <fe/fe_system.h>
double radius = cell->face(f)->center().norm() - center.norm();
if (std::fabs(cell->face(f)->center()(2) - sqrt(3.)*cell->face(f)->center()(0)) < eps ) // z = sqrt(3)x set boundary 2
{
- cell->face(f)->set_boundary_indicator(2);
+ cell->face(f)->set_manifold_indicator(2);
for (unsigned int j=0; j<GeometryInfo<3>::lines_per_face; ++j)
if (cell->face(f)->line(j)->at_boundary())
if (std::fabs(cell->face(f)->line(j)->vertex(0).norm() - cell->face(f)->line(j)->vertex(1).norm()) > eps)
- cell->face(f)->line(j)->set_boundary_indicator(2);
+ cell->face(f)->line(j)->set_manifold_indicator(2);
}
else if (std::fabs(cell->face(f)->center()(2) + sqrt(3.)*cell->face(f)->center()(0)) < eps) // z = -sqrt(3)x set boundary 3
{
- cell->face(f)->set_boundary_indicator(3);
+ cell->face(f)->set_manifold_indicator(3);
for (unsigned int j=0; j<GeometryInfo<3>::lines_per_face; ++j)
if (cell->face(f)->line(j)->at_boundary())
if (std::fabs(cell->face(f)->line(j)->vertex(0).norm() - cell->face(f)->line(j)->vertex(1).norm()) > eps)
- cell->face(f)->line(j)->set_boundary_indicator(3);
+ cell->face(f)->line(j)->set_manifold_indicator(3);
}
else if (std::fabs(cell->face(f)->center()(2) - sqrt(3.)*cell->face(f)->center()(1)) < eps ) // z = sqrt(3)y set boundary 4
{
- cell->face(f)->set_boundary_indicator(4);
+ cell->face(f)->set_manifold_indicator(4);
for (unsigned int j=0; j<GeometryInfo<3>::lines_per_face; ++j)
if (cell->face(f)->line(j)->at_boundary())
if (std::fabs(cell->face(f)->line(j)->vertex(0).norm() - cell->face(f)->line(j)->vertex(1).norm()) > eps)
- cell->face(f)->line(j)->set_boundary_indicator(4);
+ cell->face(f)->line(j)->set_manifold_indicator(4);
}
else if (std::fabs(cell->face(f)->center()(2) + sqrt(3.)*cell->face(f)->center()(1)) < eps ) // z = -sqrt(3)y set boundary 5
{
- cell->face(f)->set_boundary_indicator(5);
+ cell->face(f)->set_manifold_indicator(5);
for (unsigned int j=0; j<GeometryInfo<3>::lines_per_face; ++j)
if (cell->face(f)->line(j)->at_boundary())
if (std::fabs(cell->face(f)->line(j)->vertex(0).norm() - cell->face(f)->line(j)->vertex(1).norm()) > eps)
- cell->face(f)->line(j)->set_boundary_indicator(5);
+ cell->face(f)->line(j)->set_manifold_indicator(5);
}
else if (radius < middle) // inner radius set boundary 0
{
- cell->face(f)->set_boundary_indicator(0);
+ cell->face(f)->set_manifold_indicator(0);
for (unsigned int j=0; j<GeometryInfo<3>::lines_per_face; ++j)
if (cell->face(f)->line(j)->at_boundary())
if (std::fabs(cell->face(f)->line(j)->vertex(0).norm() - cell->face(f)->line(j)->vertex(1).norm()) < eps)
- cell->face(f)->line(j)->set_boundary_indicator(0);
+ cell->face(f)->line(j)->set_manifold_indicator(0);
}
else if (radius > middle) // outer radius set boundary 1
{
- cell->face(f)->set_boundary_indicator(1);
+ cell->face(f)->set_manifold_indicator(1);
for (unsigned int j=0; j<GeometryInfo<3>::lines_per_face; ++j)
if (cell->face(f)->line(j)->at_boundary())
if (std::fabs(cell->face(f)->line(j)->vertex(0).norm() - cell->face(f)->line(j)->vertex(1).norm()) < eps)
- cell->face(f)->line(j)->set_boundary_indicator(1);
+ cell->face(f)->line(j)->set_manifold_indicator(1);
}
else
AssertThrow (false, ExcInternalError());
0.5,
1.0);
- static HyperShellBoundary<dim> boundary((Point<dim>()));
- triangulation.set_boundary (0, boundary);
+ static PolarManifold<dim> boundary((Point<dim>()));
+ triangulation.set_manifold (0, boundary);
// write out the mesh. may not be
// strictly necessary here, but is
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/hp/dof_handler.h>
#include <deal.II/hp/fe_collection.h>
#include <deal.II/lac/constraint_matrix.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
tr.refine_global(2);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/hp/dof_handler.h>
#include <deal.II/hp/fe_collection.h>
#include <deal.II/lac/constraint_matrix.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
tr.refine_global(2);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/hp/dof_handler.h>
#include <deal.II/dofs/dof_renumbering.h>
#include <deal.II/lac/constraint_matrix.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
tr.refine_global(1);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/constraint_matrix.h>
#include <deal.II/fe/fe_q.h>
Triangulation<dim> tr;
GridGenerator::hyper_ball(tr);
- static const HyperBallBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
tr.refine_global(2);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_tools.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/hp/dof_handler.h>
#include <deal.II/lac/constraint_matrix.h>
template<int dim>
void test()
{
- HyperBallBoundary<dim> boundary;
+ PolarManifold<dim> boundary;
// create 2 triangulations with the
// same coarse grid, and refine
// them differently
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- tria.set_boundary (0, boundary);
+ tria.set_manifold (0, boundary);
FE_Q<dim> fe(1);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/base/logstream.h>
#include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/grid_generator.h>
deallog << "dim=" << dim << std::endl;
Triangulation<dim> tria;
- StraightBoundary<dim> boundary;
+ FlatManifold<dim> boundary;
for (unsigned int case_no=0; case_no<2; ++case_no)
{
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/base/logstream.h>
#include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/grid_generator.h>
deallog << "dim=" << dim << std::endl;
Triangulation<dim> tria;
- StraightBoundary<dim> boundary;
+ FlatManifold<dim> boundary;
for (unsigned int case_no=0; case_no<2; ++case_no)
{
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/base/logstream.h>
#include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/grid_generator.h>
deallog << "dim=" << dim << std::endl;
Triangulation<dim> tria;
- StraightBoundary<dim> boundary;
+ FlatManifold<dim> boundary;
GridGenerator::hyper_cube(tria, 0, 1);
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgq.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_monomial.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_monomial.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_monomial.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_monomial.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_monomial.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_monomial.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_monomial.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_monomial.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_monomial.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_monomial.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_nonparametric.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_nonparametric.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_nonparametric.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_nonparametric.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_nonparametric.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_nonparametric.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_nonparametric.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_nonparametric.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_nonparametric.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/lac/vector.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_dgp_nonparametric.h>
subdivisions[0] = 2;
GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
- static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
- tr.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary(tr.begin_active()->center());
+ tr.set_manifold (1, boundary);
// set boundary id on cell 1
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (tr.begin_active()->at_boundary(f))
- tr.begin_active()->face(f)->set_boundary_indicator (1);
+ tr.begin_active()->face(f)->set_manifold_indicator (1);
test(tr);
}
#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
template<int dim>
void test()
{
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
MappingQ<dim> mapping(3);
FE_Nothing<dim> dummy;
#include <deal.II/dofs/dof_accessor.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_out.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/mapping_cartesian.h>
#include <deal.II/fe/mapping_q1.h>
#include <deal.II/fe/mapping_q.h>
v2(1) = 3.;
v3(0) = 3.;
v3(1) = 3.;
- tria->set_boundary(1,*boundary1);
- tria->set_boundary(2,*boundary2);
- tria->begin_active()->face(1)->set_boundary_indicator(1);
- tria->begin_active()->face(3)->set_boundary_indicator(2);
+ tria->set_manifold(1,*boundary1);
+ tria->set_manifold(2,*boundary2);
+ tria->begin_active()->face(1)->set_manifold_indicator(1);
+ tria->begin_active()->face(3)->set_manifold_indicator(2);
double pi=std::acos(-1.);
double alpha=2*std::atan(0.5);
exact_areas.push_back(4+pi-2.5*(alpha-std::sin(alpha)));
v2(1) = 1.;
v3(0) = 0.5;
v3(1) = 1.5;
- tria->set_boundary(1,*boundary1);
- tria->begin_active()->face(1)->set_boundary_indicator(1);
+ tria->set_manifold(1,*boundary1);
+ tria->begin_active()->face(1)->set_manifold_indicator(1);
exact_areas.push_back(0.);
for (unsigned int i=0; i<=4; ++i)
show[4][i]=1;
tria=new Triangulation<3>();
tria_ptr.push_back(tria);
GridGenerator::hyper_cube(*tria, 1., 3.);
- tria->set_boundary(1,*boundary1);
- tria->begin_active()->face(1)->set_boundary_indicator(1);
+ tria->set_manifold(1,*boundary1);
+ tria->begin_active()->face(1)->set_manifold_indicator(1);
exact_areas.push_back(8.+pi/3*h*h*(3*r-h));
}
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/numerics/fe_field_function.h>
template<int dim>
void test2()
{
- const HyperBallBoundary<dim> boundary_description;
+ const PolarManifold<dim> boundary_description;
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- triangulation.set_boundary (0, boundary_description);
+ triangulation.set_manifold (0, boundary_description);
triangulation.refine_global (1);
MappingQ1< dim > mapping;
points[95]=p;
fe_function.value_list (points, m); // <<<< this fails at point[95] but only if the other points are filled in?!
- triangulation.set_boundary (0);
+ triangulation.set_manifold (0);
deallog << "OK" << std::endl;
}
#include <deal.II/base/utilities.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/mapping_q.h>
// set the boundary indicator for
// one face of the single cell
- triangulation.set_boundary (1, boundary);
- triangulation.begin_active()->face(0)->set_boundary_indicator (1);
+ triangulation.set_manifold (1, boundary);
+ triangulation.begin_active()->face(0)->set_manifold_indicator (1);
const unsigned int n_points = 5;
std::vector< Point<dim> > unit_points(Utilities::fixed_power<dim>(n_points));
#include <deal.II/base/utilities.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/mapping_q.h>
// set the boundary indicator for
// one face of the single cell
- triangulation.set_boundary (1, boundary);
- triangulation.begin_active()->face(0)->set_boundary_indicator (1);
+ triangulation.set_manifold (1, boundary);
+ triangulation.begin_active()->face(0)->set_manifold_indicator (1);
const unsigned int n_points = 5;
std::vector< Point<dim> > unit_points(Utilities::fixed_power<dim>(n_points));
#include <deal.II/base/utilities.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/mapping_q.h>
// define a spherical cap boundary
// to be used for one of the faces
const double radius = Point<dim>(1.43757e-10, 4.48023e+06, 4.48023e+06).norm();
- HyperBallBoundary<dim> boundary (Point<dim>(), radius);
+ PolarManifold<dim> boundary (Point<dim>(), radius);
// create the mesh: a single cell
// with the following coordinates:
// set the boundary indicator for
// one face and adjacent edges of
// the single cell
- triangulation.set_boundary (1, boundary);
+ triangulation.set_manifold (1, boundary);
triangulation.begin_active()->face(5)->set_all_boundary_indicators (1);
// now try to find the coordinates
#include <deal.II/base/utilities.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/mapping_q.h>
// define a spherical cap boundary
// to be used for one of the faces
const double radius = Point<dim>(1.43757e-10, 4.48023e+06, 4.48023e+06).norm();
- HyperBallBoundary<dim> boundary (Point<dim>(), radius);
+ PolarManifold<dim> boundary (Point<dim>(), radius);
// create the mesh: a single cell
// with the following coordinates:
// set the boundary indicator for
// one face and adjacent edges of
// the single cell
- triangulation.set_boundary (1, boundary);
+ triangulation.set_manifold (1, boundary);
triangulation.begin_active()->face(5)->set_all_boundary_indicators (1);
// now try to find the coordinates
#include <deal.II/base/utilities.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/mapping_q.h>
// define a spherical cap boundary
// to be used for one of the faces
const double radius = Point<dim>(1.43757e-10, 4.48023e+06, 4.48023e+06).norm();
- HyperBallBoundary<dim> boundary (Point<dim>(), radius);
+ PolarManifold<dim> boundary (Point<dim>(), radius);
// create the mesh: a single cell
// with the following coordinates:
// set the boundary indicator for
// one face and adjacent edges of
// the single cell
- triangulation.set_boundary (1, boundary);
+ triangulation.set_manifold (1, boundary);
triangulation.begin_active()->face(5)->set_all_boundary_indicators (1);
// now try to find the coordinates
#include <deal.II/base/utilities.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/fe/mapping_q.h>
// define a spherical cap boundary
// to be used for one of the faces
const double radius = Point<dim>(1.43757e-10, 4.48023e+06, 4.48023e+06).norm();
- HyperBallBoundary<dim> boundary (Point<dim>(), radius);
+ PolarManifold<dim> boundary (Point<dim>(), radius);
// create the mesh: a single cell
// with the following coordinates:
// set the boundary indicator for
// one face and adjacent edges of
// the single cell
- triangulation.set_boundary (1, boundary);
+ triangulation.set_manifold (1, boundary);
triangulation.begin_active()->face(5)->set_all_boundary_indicators (1);
// now see if the point is inside
#include <deal.II/base/tensor.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/grid_out.h>
#include <fstream>
if (dim>1) p3[1] = 1.;
if (dim>2) p3[2] = 4.;
- HalfHyperBallBoundary<dim> boundary_description (p1, 3);
+ HalfPolarManifold<dim> boundary_description (p1, 3);
GridOut go;
GridOut::OutputFormat format = GridOut::gnuplot;
deallog << "half_hyper_ball" << std::endl;
Triangulation<dim> tr;
GridGenerator::half_hyper_ball(tr, p1, 3.);
- tr.set_boundary (0, boundary_description);
+ tr.set_manifold (0, boundary_description);
tr.refine_global (2);
go.write(tr, out, format);
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_tools.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/hp/dof_handler.h>
#include <deal.II/hp/q_collection.h>
template<int dim>
void test()
{
- HyperBallBoundary<dim> boundary;
+ PolarManifold<dim> boundary;
// create 2 triangulations with the
// same coarse grid, and refine
// them differently
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- tria.set_boundary (0, boundary);
+ tria.set_manifold (0, boundary);
hp::FECollection<dim> fe;
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/hp/dof_handler.h>
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
for (unsigned int refinement=0; refinement<2;
++refinement, triangulation.refine_global(1))
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
MappingQ<dim> m(degree);
const hp::MappingCollection<dim> mapping (m);
Triangulation<dim> triangulation;
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
const MappingQ<dim> m (degree);
const hp::MappingCollection<dim> mapping(m);
#include <deal.II/lac/precondition.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/dofs/dof_handler.h>
void LaplaceProblem<dim>::run ()
{
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
for (unsigned int cycle=0; cycle<6; ++cycle, triangulation.refine_global(1))
{
#include <deal.II/lac/precondition.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/dofs/dof_handler.h>
void LaplaceProblem<dim>::run ()
{
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
for (unsigned int cycle=0; cycle<6; ++cycle, triangulation.refine_global(1))
{
#include <deal.II/grid/grid_in.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <fstream>
#include <sstream>
grid_in.read_ucd (input_file);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
}
deallog << " Number of active cells: "
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_accessor.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/hp/fe_values.h>
{
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
triangulation.refine_global (1);
}
{
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
// refine first and last cell
tria.begin(tria.n_levels()-1)->set_refine_flag();
tria.last()->set_refine_flag();
typedef double number;
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
// refine first and last cell
tria.begin(tria.n_levels()-1)->set_refine_flag();
tria.last()->set_refine_flag();
typedef double number;
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
// refine first and last cell
tria.begin(tria.n_levels()-1)->set_refine_flag();
tria.last()->set_refine_flag();
typedef double number;
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
// refine first and last cell
tria.begin(tria.n_levels()-1)->set_refine_flag();
tria.last()->set_refine_flag();
typedef double number;
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
// refine first and last cell
tria.begin(tria.n_levels()-1)->set_refine_flag();
tria.last()->set_refine_flag();
#include <deal.II/lac/vector.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/constraint_matrix.h>
{
Triangulation<dim> tria;
GridGenerator::hyper_shell (tria, Point<dim>(), 1., 2., 96, true);
- static const HyperShellBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
- tria.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
+ tria.set_manifold (1, boundary);
// refine a few cells
for (unsigned int i=0; i<11-3*dim; ++i)
#include <deal.II/lac/vector.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/constraint_matrix.h>
typedef double number;
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
typename Triangulation<dim>::active_cell_iterator
cell = tria.begin_active (),
endc = tria.end();
#include <deal.II/lac/vector.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/constraint_matrix.h>
// cells
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
typename Triangulation<dim>::active_cell_iterator
cell = tria.begin_active (),
endc = tria.end();
Triangulation<dim> tria;
GridGenerator::hyper_shell (tria, Point<dim>(),
0.5, 1., 96, true);
- static const HyperShellBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
- tria.set_boundary (1, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
+ tria.set_manifold (1, boundary);
if (dim == 2)
tria.refine_global (2);
{
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
typename Triangulation<dim>::active_cell_iterator
cell = tria.begin_active (),
endc = tria.end();
#include <deal.II/lac/block_vector.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/dofs/dof_renumbering.h>
{
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
tria.refine_global(4-dim);
// refine a few cells
typedef double number;
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
tria.refine_global(1);
// refine a few cells
{
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
tria.refine_global(5-dim);
FE_Q<dim> fe (fe_degree);
#include <deal.II/lac/block_vector.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/dofs/dof_renumbering.h>
Triangulation<dim> triangulation;
GridGenerator::hyper_shell (triangulation, Point<dim>(),
0.5, 1., 96, true);
- static HyperShellBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
- triangulation.set_boundary (1, boundary);
+ static PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
+ triangulation.set_manifold (1, boundary);
triangulation.begin_active()->set_refine_flag();
triangulation.last()->set_refine_flag();
triangulation.execute_coarsening_and_refinement();
#include <deal.II/lac/vector.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/constraint_matrix.h>
{
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
// refine first and last cell
tria.begin(tria.n_levels()-1)->set_refine_flag();
tria.last()->set_refine_flag();
#include <deal.II/lac/vector.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/constraint_matrix.h>
typedef double number;
Triangulation<dim> tria;
GridGenerator::hyper_ball (tria);
- static const HyperBallBoundary<dim> boundary;
- tria.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ tria.set_manifold (0, boundary);
tria.refine_global(5-dim);
MappingQ<dim> mapping (4);
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/grid/grid_refinement.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_accessor.h>
#include <deal.II/dofs/dof_tools.h>
{
GridGenerator::hyper_cube (triangulation);
- // static const HyperBallBoundary<dim> boundary;
- // triangulation.set_boundary (0, boundary);
+ // static const PolarManifold<dim> boundary;
+ // triangulation.set_manifold (0, boundary);
triangulation.refine_global (2);
}
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/grid/grid_refinement.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_accessor.h>
#include <deal.II/dofs/dof_tools.h>
{
GridGenerator::hyper_cube (triangulation);
- // static const HyperBallBoundary<dim> boundary;
- // triangulation.set_boundary (0, boundary);
+ // static const PolarManifold<dim> boundary;
+ // triangulation.set_manifold (0, boundary);
triangulation.refine_global (3);
}
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/dofs/dof_accessor.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/numerics/data_out.h>
#include <deal.II/lac/trilinos_vector.h>
R1,
12,
true);
- static HyperShellBoundary<dim> boundary;
- tr.set_boundary (0, boundary);
- tr.set_boundary (1, boundary);
+ static PolarManifold<dim> boundary;
+ tr.set_manifold (0, boundary);
+ tr.set_manifold (1, boundary);
tr.refine_global (1);
for (unsigned int step=0; step<20; ++step)
cat_file((std::string("dat.") + Utilities::int_to_string(0)).c_str());
}
- tr.set_boundary (0);
- tr.set_boundary (1);
+ tr.set_manifold (0);
+ tr.set_manifold (1);
}
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/dofs/dof_accessor.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/numerics/data_out.h>
#include <deal.II/lac/trilinos_vector.h>
R1,
12,
true);
- static HyperShellBoundary<dim> boundary;
-// tr.set_boundary (0, boundary);
-// tr.set_boundary (1, boundary);
+ static PolarManifold<dim> boundary;
+// tr.set_manifold (0, boundary);
+// tr.set_manifold (1, boundary);
tr.refine_global (3);
if (1)
std::ofstream output (filename.c_str());
data_out.write_deal_II_intermediate (output);
- tr.set_boundary (0);
- tr.set_boundary (1);
+ tr.set_manifold (0);
+ tr.set_manifold (1);
if (Utilities::MPI::this_mpi_process (MPI_COMM_WORLD) == 0)
deallog << "OK" << std::endl;
}
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/dofs/dof_accessor.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/numerics/data_out.h>
#include <deal.II/lac/trilinos_block_vector.h>
R1,
12,
true);
- static HyperShellBoundary<dim> boundary;
-// tr.set_boundary (0, boundary);
-//tr.set_boundary (1, boundary);
+ static PolarManifold<dim> boundary;
+// tr.set_manifold (0, boundary);
+//tr.set_manifold (1, boundary);
tr.refine_global (1);
if (1)
std::ofstream output (filename.c_str());
data_out.write_vtu (output);
- tr.set_boundary (0);
- tr.set_boundary (1);
+ tr.set_manifold (0);
+ tr.set_manifold (1);
if (Utilities::MPI::this_mpi_process (MPI_COMM_WORLD) == 0)
deallog << "OK" << std::endl;
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_refinement.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/integrators/laplace.h>
{
GridGenerator::hyper_cube (triangulation, -1, 1);
-// static const HyperBallBoundary<dim> boundary;
-// triangulation.set_boundary (0, boundary);
+// static const PolarManifold<dim> boundary;
+// triangulation.set_manifold (0, boundary);
triangulation.refine_global (1);
}
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_refinement.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
#include <deal.II/dofs/dof_accessor.h>
#include <deal.II/dofs/dof_tools.h>
{
GridGenerator::hyper_ball (triangulation);
- static const HyperBallBoundary<dim> boundary;
- triangulation.set_boundary (0, boundary);
+ static const PolarManifold<dim> boundary;
+ triangulation.set_manifold (0, boundary);
triangulation.refine_global (1);
}