#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim < 3 || fe_degree < 2)
tria.refine_global(1);
#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim == 2)
tria.refine_global(1);
#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
tria.refine_global(5 - dim);
tria.begin(tria.n_levels() - 1)->set_refine_flag();
#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
void
test()
{
- SphericalManifold<dim> manifold;
- Triangulation<dim> triangulation;
+ Triangulation<dim> triangulation;
GridGenerator::hyper_shell(triangulation, Point<dim>(), 0.5, 1., 96, true);
- triangulation.set_all_manifold_ids(0);
- triangulation.set_manifold(0, manifold);
triangulation.begin_active()->set_refine_flag();
triangulation.last()->set_refine_flag();
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
// refine first and last cell
tria.begin(tria.n_levels() - 1)->set_refine_flag();
test()
{
using number = double;
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
// refine first and last cell
tria.begin(tria.n_levels() - 1)->set_refine_flag();
test()
{
using number = double;
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
// refine first and last cell
tria.begin(tria.n_levels() - 1)->set_refine_flag();
test()
{
using number = double;
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
// refine first and last cell
tria.begin(tria.n_levels() - 1)->set_refine_flag();
test()
{
using number = double;
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
// refine first and last cell
tria.begin(tria.n_levels() - 1)->set_refine_flag();
test()
{
using number = double;
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
// refine first and last cell
tria.begin(tria.n_levels() - 1)->set_refine_flag();
#include <deal.II/fe/fe_values.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_shell(tria, Point<dim>(), 1., 2., 96, true);
- tria.set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
// refine a few cells
for (unsigned int i = 0; i < 11 - 3 * dim; ++i)
{
- typename Triangulation<dim>::active_cell_iterator cell =
- tria.begin_active(),
- endc = tria.end();
- unsigned int counter = 0;
- for (; cell != endc; ++cell, ++counter)
- if (counter % (7 - i) == 0)
- cell->set_refine_flag();
+ unsigned int counter = 0;
+ for (const auto &cell : tria.active_cell_iterators())
+ {
+ if (counter % (7 - i) == 0)
+ cell->set_refine_flag();
+ ++counter;
+ }
tria.execute_coarsening_and_refinement();
}
#include <deal.II/fe/fe_values.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
test()
{
using number = double;
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
-
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 0.2)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
tria.begin(tria.n_levels() - 1)->set_refine_flag();
tria.last()->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
for (unsigned int i = 0; i < 7 - 2 * dim; ++i)
{
- cell = tria.begin_active();
unsigned int counter = 0;
- for (; cell != endc; ++cell, ++counter)
- if (counter % (7 - i) == 0)
- cell->set_refine_flag();
+ for (const auto &cell : tria.active_cell_iterators())
+ {
+ if (counter % (7 - i) == 0)
+ cell->set_refine_flag();
+ ++counter;
+ }
tria.execute_coarsening_and_refinement();
}
#include <deal.II/fe/fe_values.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
test()
{
using number = double;
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
-
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 0.2)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
tria.begin(tria.n_levels() - 1)->set_refine_flag();
tria.last()->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
for (unsigned int i = 0; i < 7 - 2 * dim; ++i)
{
- cell = tria.begin_active();
unsigned int counter = 0;
- for (; cell != endc; ++cell, ++counter)
- if (counter % (7 - i) == 0)
- cell->set_refine_flag();
+ for (const auto &cell : tria.active_cell_iterators())
+ {
+ if (counter % (7 - i) == 0)
+ cell->set_refine_flag();
+ ++counter;
+ }
tria.execute_coarsening_and_refinement();
}
#include <deal.II/fe/fe_values.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
{
// create hyper ball geometry and refine some
// cells
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
-
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 0.2)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
tria.begin(tria.n_levels() - 1)->set_refine_flag();
tria.last()->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
for (unsigned int i = 0; i < 7 - 2 * dim; ++i)
{
- cell = tria.begin_active();
unsigned int counter = 0;
- for (; cell != endc; ++cell, ++counter)
- if (counter % (7 - i) == 0)
- cell->set_refine_flag();
+ for (const auto &cell : tria.active_cell_iterators())
+ {
+ if (counter % (7 - i) == 0)
+ cell->set_refine_flag();
+ ++counter;
+ }
tria.execute_coarsening_and_refinement();
}
#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/precondition.h>
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim < 3 || fe_degree < 2)
tria.refine_global(1);
tria.begin(tria.n_levels() - 1)->set_refine_flag();
tria.last()->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/precondition.h>
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim < 3 || fe_degree < 2)
tria.refine_global(1);
tria.begin(tria.n_levels() - 1)->set_refine_flag();
tria.last()->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/precondition.h>
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim < 3 || fe_degree < 2)
tria.refine_global(1);
tria.begin(tria.n_levels() - 1)->set_refine_flag();
tria.last()->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/precondition.h>
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim < 3 || fe_degree < 2)
tria.refine_global(1);
tria.begin(tria.n_levels() - 1)->set_refine_flag();
tria.last()->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/precondition.h>
void
test(const unsigned int fe_degree)
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- for (const auto &cell : tria.active_cell_iterators())
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim < 3 || fe_degree < 2)
tria.refine_global(1);
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_shell(tria, Point<dim>(), 0.5, 1., 96, true);
- tria.set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim == 2)
tria.refine_global(2);
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_shell(tria, Point<dim>(), 0.5, 1., 96, true);
- tria.set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim == 2)
tria.refine_global(2);
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 0.2)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
tria.execute_coarsening_and_refinement();
tria.refine_global(4 - dim);
#endif
- cell = tria.begin_active();
for (unsigned int i = 0; i < 10 - 3 * dim; ++i)
{
- cell = tria.begin_active();
unsigned int counter = 0;
- for (; cell != endc; ++cell, ++counter)
- if (counter % (7 - i) == 0)
- cell->set_refine_flag();
+ for (const auto &cell : tria.active_cell_iterators())
+ {
+ if (counter % (7 - i) == 0)
+ cell->set_refine_flag();
+ ++counter;
+ }
tria.execute_coarsening_and_refinement();
}
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim < 3 || fe_degree < 2)
tria.refine_global(1);
tria.begin(tria.n_levels() - 1)->set_refine_flag();
tria.last()->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim < 3 || fe_degree < 2)
tria.refine_global(1);
tria.begin(tria.n_levels() - 1)->set_refine_flag();
tria.last()->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim < 3 || fe_degree < 2)
tria.refine_global(1);
tria.begin(tria.n_levels() - 1)->set_refine_flag();
tria.last()->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
#include <deal.II/fe/fe_values.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim < 3 || fe_degree < 2)
tria.refine_global(1);
tria.begin(tria.n_levels() - 1)->set_refine_flag();
tria.last()->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
#include <deal.II/fe/fe_values.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
-
if (dim < 3 || fe_degree < 2)
tria.refine_global(1);
tria.begin(tria.n_levels() - 1)->set_refine_flag();
tria.last()->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
if (dim < 3 || fe_degree < 2)
tria.refine_global(1);
tria.begin(tria.n_levels() - 1)->set_refine_flag();
tria.last()->set_refine_flag();
tria.execute_coarsening_and_refinement();
- cell = tria.begin_active();
- for (; cell != endc; ++cell)
+ for (const auto &cell : tria.active_cell_iterators())
if (cell->center().norm() < 1e-8)
cell->set_refine_flag();
tria.execute_coarsening_and_refinement();
#include <deal.II/fe/fe_values.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
tria.refine_global(4 - dim);
// refine a few cells
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/hp/fe_values.h>
test()
{
using number = double;
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- for (const auto &cell : tria.active_cell_iterators())
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
tria.refine_global(1);
// refine a few cells
return;
using number = double;
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
tria.refine_global(1);
// refine a few cells
for (unsigned int i = 0; i < 11 - 3 * dim; ++i)
{
- typename Triangulation<dim>::active_cell_iterator cell =
- tria.begin_active(),
- endc = tria.end();
- unsigned int counter = 0;
- for (; cell != endc; ++cell, ++counter)
- if (counter % (7 - i) == 0)
- cell->set_refine_flag();
+ unsigned int counter = 0;
+ for (const auto &cell : tria.active_cell_iterators())
+ {
+ if (counter % (7 - i) == 0)
+ cell->set_refine_flag();
+ ++counter;
+ }
tria.execute_coarsening_and_refinement();
}
#include <deal.II/fe/fe_values.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/hp/fe_values.h>
return;
using number = double;
- const SphericalManifold<dim> manifold;
parallel::distributed::Triangulation<dim> tria(MPI_COMM_WORLD);
GridGenerator::hyper_ball(tria);
- for (const auto &cell : tria.active_cell_iterators())
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
tria.refine_global(1);
// refine a few cells
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria(
+ Triangulation<dim> tria(
Triangulation<dim>::limit_level_difference_at_vertices);
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
tria.refine_global(5 - dim);
FE_Q<dim> fe(fe_degree);
#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
void
test()
{
- SphericalManifold<dim> manifold;
- Triangulation<dim> triangulation;
+ Triangulation<dim> triangulation;
GridGenerator::hyper_shell(triangulation, Point<dim>(), 0.5, 1., 96, true);
- triangulation.set_all_manifold_ids(0);
- triangulation.set_manifold(0, manifold);
triangulation.begin_active()->set_refine_flag();
triangulation.last()->set_refine_flag();
#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
void
test(const unsigned int fe_degree)
{
- SphericalManifold<dim> manifold;
- Triangulation<dim> triangulation;
+ Triangulation<dim> triangulation;
GridGenerator::hyper_shell(triangulation, Point<dim>(), 0.5, 1., 96, true);
- triangulation.set_all_manifold_ids(0);
- triangulation.set_manifold(0, manifold);
triangulation.begin_active()->set_refine_flag();
triangulation.last()->set_refine_flag();
#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
void
test(const FESystem<dim> &fe)
{
- SphericalManifold<dim> manifold;
- Triangulation<dim> triangulation;
+ Triangulation<dim> triangulation;
GridGenerator::hyper_shell(triangulation, Point<dim>(), 0.5, 1., 96, true);
- triangulation.set_all_manifold_ids(0);
- triangulation.set_manifold(0, manifold);
triangulation.begin_active()->set_refine_flag();
triangulation.last()->set_refine_flag();
triangulation.execute_coarsening_and_refinement();
#include <deal.II/fe/fe_values.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
void
test()
{
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
// refine first and last cell
tria.begin(tria.n_levels() - 1)->set_refine_flag();
#include <deal.II/fe/mapping_q.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
test()
{
using number = double;
- const SphericalManifold<dim> manifold;
- Triangulation<dim> tria;
+ Triangulation<dim> tria;
GridGenerator::hyper_ball(tria);
- typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
- endc = tria.end();
- for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
- if (cell->at_boundary(f))
- cell->face(f)->set_all_manifold_ids(0);
- tria.set_manifold(0, manifold);
tria.refine_global(5 - dim);
MappingQ<dim> mapping(4);