* @ingroup boundary
*
* @deprecated The boundary classes have been deprecated in favor of the more
- * general Manifold classes. For this class CylindricalManifold is an appropriate
+ * general Manifold classes. For this class SphericalManifold is an appropriate
* replacement which is expected to produce better results under mesh refinement.
*
* @author Wolfgang Bangerth, 1999
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.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>
center, inner_radius, outer_radius,
10);
- static const HyperShellBoundary<2> boundary_description(center);
- triangulation.set_boundary (0, boundary_description);
+ static const SphericalManifold<2> boundary_description(center);
+ triangulation.set_manifold (0, boundary_description);
for (unsigned int step=0; step<5; ++step)
{
// 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>
for (unsigned int e=0; e<GeometryInfo<dim>::lines_per_face; ++e)
c->face(f)->line(e)->set_boundary_id(0);
- static const HyperShellBoundary<dim> boundary(center);
- tria.set_boundary(0, boundary);
+ static const SphericalManifold<dim> boundary(center);
+ tria.set_manifold(0, boundary);
for (unsigned int i=0; i<2; ++i)
{
3.64058 -2.64503
3.80423 1.23607
2.35114 3.23607
- 1.16736e-16 4.00000
+ 8.88178e-16 4.00000
-2.35114 3.23607
-3.80423 1.23607
-3.80423 -1.23607
-2.35114 -3.23607
- -7.00415e-16 -4.00000
+ -4.44089e-16 -4.00000
2.35114 -3.23607
3.80423 -1.23607
4.16338 1.35277
2.57311 3.54159
- 1.69390e-16 4.37764
+ 5.55112e-16 4.37764
-2.57311 3.54159
-4.16338 1.35277
-4.16338 -1.35277
-2.57311 -3.54159
- -7.38785e-16 -4.37764
+ -6.10623e-16 -4.37764
2.57311 -3.54159
4.16338 -1.35277
4.76127 0.734732
2.46213 3.38883
3.10684 3.09331
1.98184 3.91067
- 1.95717e-16 4.56646
- 1.43063e-16 4.18882
+ 3.88578e-16 4.56646
+ 7.21645e-16 4.18882
0.695288 4.32870
-0.695288 4.32870
-2.68410 3.69435
-2.46213 -3.38883
-3.10684 -3.09331
-1.98184 -3.91067
- -7.57971e-16 -4.56646
- -7.19600e-16 -4.18882
+ -6.93889e-16 -4.56646
+ -5.27356e-16 -4.18882
-0.695288 -4.32870
0.695288 -4.32870
2.68410 -3.69435
2.82843 0.00000 2.82843
0.00000 2.82843 2.82843
0.00000 0.00000 -3.59216
- 3.59216 0.00000 -4.44089e-16
- 0.00000 -4.44089e-16 3.59216
+ 3.59216 0.00000 -2.22045e-16
+ 0.00000 -2.22045e-16 3.59216
-3.59216 0.00000 0.00000
- 0.00000 -3.59216 -4.44089e-16
+ 0.00000 -3.59216 -2.22045e-16
0.00000 3.59216 0.00000
-0.908716 -2.02166 -2.02166
0.908716 -2.02166 -2.02166
-1.41421 3.21029 0.00000
0.00000 -3.21029 1.41421
0.00000 -3.21029 -1.41421
- 1.41421 -3.21029 -2.22045e-16
- -1.41421 -3.21029 -2.22045e-16
- 0.00000 -3.29608 -2.22045e-16
- 0.00000 -3.88824 -2.22045e-16
+ 1.41421 -3.21029 -1.11022e-16
+ -1.41421 -3.21029 -1.11022e-16
+ 0.00000 -3.29608 -1.11022e-16
+ 0.00000 -3.88824 -1.11022e-16
-3.21029 0.00000 1.41421
-3.21029 0.00000 -1.41421
-3.29608 0.00000 0.00000
-3.21029 -1.41421 0.00000
0.00000 1.41421 3.21029
0.00000 -1.41421 3.21029
- 1.41421 -2.22045e-16 3.21029
- -1.41421 -2.22045e-16 3.21029
- 0.00000 -2.22045e-16 3.29608
- 0.00000 -2.22045e-16 3.88824
+ 1.41421 -1.11022e-16 3.21029
+ -1.41421 -1.11022e-16 3.21029
+ 0.00000 -1.11022e-16 3.29608
+ 0.00000 -1.11022e-16 3.88824
3.21029 0.00000 1.41421
3.21029 0.00000 -1.41421
- 3.21029 1.41421 -2.22045e-16
- 3.21029 -1.41421 -2.22045e-16
- 3.29608 0.00000 -2.22045e-16
- 3.88824 0.00000 -2.22045e-16
+ 3.21029 1.41421 -1.11022e-16
+ 3.21029 -1.41421 -1.11022e-16
+ 3.29608 0.00000 -1.11022e-16
+ 3.88824 0.00000 -1.11022e-16
0.00000 0.00000 -3.29608
0.00000 0.00000 -3.88824
1.41421 0.00000 -3.21029
-1.41421 0.00000 -3.21029
0.00000 1.41421 -3.21029
0.00000 -1.41421 -3.21029
- -1.09719 -2.56755 -1.09719
- -1.09719 -2.56755 1.09719
- 1.09719 -2.56755 -1.09719
- 1.09719 -2.56755 1.09719
- -1.09719 -1.09719 -2.56755
- 1.09719 -1.09719 -2.56755
- -1.09719 1.09719 -2.56755
- 1.09719 1.09719 -2.56755
- -2.56755 -1.09719 -1.09719
- -2.56755 1.09719 -1.09719
- -2.56755 -1.09719 1.09719
- -2.56755 1.09719 1.09719
- 2.56755 -1.09719 -1.09719
- 2.56755 -1.09719 1.09719
- 2.56755 1.09719 -1.09719
- 2.56755 1.09719 1.09719
- -1.09719 2.56755 -1.09719
- 1.09719 2.56755 -1.09719
- -1.09719 2.56755 1.09719
- 1.09719 2.56755 1.09719
- -1.09719 -1.09719 2.56755
- -1.09719 1.09719 2.56755
- 1.09719 -1.09719 2.56755
- 1.09719 1.09719 2.56755
+ -1.09051 -2.57324 -1.09051
+ -1.09051 -2.57324 1.09051
+ 1.09051 -2.57324 -1.09051
+ 1.09051 -2.57324 1.09051
+ -1.09051 -1.09051 -2.57324
+ 1.09051 -1.09051 -2.57324
+ -1.09051 1.09051 -2.57324
+ 1.09051 1.09051 -2.57324
+ -2.57324 -1.09051 -1.09051
+ -2.57324 1.09051 -1.09051
+ -2.57324 -1.09051 1.09051
+ -2.57324 1.09051 1.09051
+ 2.57324 -1.09051 -1.09051
+ 2.57324 -1.09051 1.09051
+ 2.57324 1.09051 -1.09051
+ 2.57324 1.09051 1.09051
+ -1.09051 2.57324 -1.09051
+ 1.09051 2.57324 -1.09051
+ -1.09051 2.57324 1.09051
+ 1.09051 2.57324 1.09051
+ -1.09051 -1.09051 2.57324
+ -1.09051 1.09051 2.57324
+ 1.09051 -1.09051 2.57324
+ 1.09051 1.09051 2.57324
-2.96918 -1.33461 -2.96918
-2.96918 1.33461 -2.96918
-2.29529 -1.03171 -2.29529
-1.28446 2.88960 1.28446
1.28446 2.88960 -1.28446
-1.28446 2.88960 -1.28446
- 1.28113 -2.99097 -1.11022e-16
- -1.28113 -2.99097 -1.11022e-16
- 1.59099 -3.53511 -1.11022e-16
- -1.59099 -3.53511 -1.11022e-16
+ 1.28113 -2.99097 -5.55112e-17
+ -1.28113 -2.99097 -5.55112e-17
+ 1.59099 -3.53511 -5.55112e-17
+ -1.59099 -3.53511 -5.55112e-17
1.28446 -2.88960 1.28446
1.28446 -2.88960 -1.28446
-1.28446 -2.88960 1.28446
-2.88960 -1.28446 1.28446
-2.88960 1.28446 -1.28446
-2.88960 -1.28446 -1.28446
- 1.28113 -1.11022e-16 2.99097
- -1.28113 -1.11022e-16 2.99097
- 1.59099 -1.11022e-16 3.53511
- -1.59099 -1.11022e-16 3.53511
+ 1.28113 -5.55112e-17 2.99097
+ -1.28113 -5.55112e-17 2.99097
+ 1.59099 -5.55112e-17 3.53511
+ -1.59099 -5.55112e-17 3.53511
1.28446 1.28446 2.88960
1.28446 -1.28446 2.88960
-1.28446 1.28446 2.88960
0.00000 1.59099 3.53511
0.00000 -1.28113 2.99097
0.00000 -1.59099 3.53511
- 2.99097 1.28113 -1.11022e-16
- 2.99097 -1.28113 -1.11022e-16
- 3.53511 1.59099 -1.11022e-16
- 3.53511 -1.59099 -1.11022e-16
+ 2.99097 1.28113 -5.55112e-17
+ 2.99097 -1.28113 -5.55112e-17
+ 3.53511 1.59099 -5.55112e-17
+ 3.53511 -1.59099 -5.55112e-17
2.88960 1.28446 1.28446
2.88960 1.28446 -1.28446
2.88960 -1.28446 1.28446
1.50238 -1.50238 -3.29172
-1.50238 1.50238 -3.29172
1.50238 1.50238 -3.29172
- -1.19083 -1.19083 -2.72858
- 1.19083 -1.19083 -2.72858
- -1.19083 1.19083 -2.72858
- 1.19083 1.19083 -2.72858
+ -1.18748 -1.18748 -2.73142
+ 1.18748 -1.18748 -2.73142
+ -1.18748 1.18748 -2.73142
+ 1.18748 1.18748 -2.73142
3.29172 -1.50238 -1.50238
3.29172 1.50238 -1.50238
- 2.72858 -1.19083 -1.19083
- 2.72858 1.19083 -1.19083
+ 2.73142 -1.18748 -1.18748
+ 2.73142 1.18748 -1.18748
3.29172 -1.50238 1.50238
3.29172 1.50238 1.50238
- 2.72858 -1.19083 1.19083
- 2.72858 1.19083 1.19083
+ 2.73142 -1.18748 1.18748
+ 2.73142 1.18748 1.18748
-1.50238 -1.50238 3.29172
1.50238 -1.50238 3.29172
- -1.19083 -1.19083 2.72858
- 1.19083 -1.19083 2.72858
+ -1.18748 -1.18748 2.73142
+ 1.18748 -1.18748 2.73142
-1.50238 1.50238 3.29172
1.50238 1.50238 3.29172
- -1.19083 1.19083 2.72858
- 1.19083 1.19083 2.72858
+ -1.18748 1.18748 2.73142
+ 1.18748 1.18748 2.73142
-3.29172 -1.50238 -1.50238
- -2.72858 -1.19083 -1.19083
+ -2.73142 -1.18748 -1.18748
-3.29172 1.50238 -1.50238
- -2.72858 1.19083 -1.19083
+ -2.73142 1.18748 -1.18748
-3.29172 -1.50238 1.50238
- -2.72858 -1.19083 1.19083
+ -2.73142 -1.18748 1.18748
-3.29172 1.50238 1.50238
- -2.72858 1.19083 1.19083
+ -2.73142 1.18748 1.18748
-1.50238 -3.29172 -1.50238
1.50238 -3.29172 -1.50238
- -1.19083 -2.72858 -1.19083
- 1.19083 -2.72858 -1.19083
+ -1.18748 -2.73142 -1.18748
+ 1.18748 -2.73142 -1.18748
-1.50238 -3.29172 1.50238
1.50238 -3.29172 1.50238
- -1.19083 -2.72858 1.19083
- 1.19083 -2.72858 1.19083
+ -1.18748 -2.73142 1.18748
+ 1.18748 -2.73142 1.18748
-1.50238 3.29172 -1.50238
- -1.19083 2.72858 -1.19083
+ -1.18748 2.73142 -1.18748
1.50238 3.29172 -1.50238
- 1.19083 2.72858 -1.19083
+ 1.18748 2.73142 -1.18748
-1.50238 3.29172 1.50238
- -1.19083 2.72858 1.19083
+ -1.18748 2.73142 1.18748
1.50238 3.29172 1.50238
- 1.19083 2.72858 1.19083
+ 1.18748 2.73142 1.18748
object "cells" class array type int rank 1 shape 8 items 384 data follows
8 102 106 298 104 306 314 442
104 306 314 442 29 176 182 441
// 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 SphericalManifold<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 SphericalManifold<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 SphericalManifold<dim> boundary(center);
+ tria.set_manifold(0, boundary);
GridOut grid_out;
grid_out.write_gnuplot (tria, logfile);
// of the boundary
#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, true);
- static const HyperShellBoundary<dim> boundary(center);
- tria.set_boundary(0, boundary);
+ static const SphericalManifold<dim> boundary(center);
+ tria.set_manifold(0, boundary);
for (typename Triangulation<dim>::cell_iterator cell=tria.begin();
cell != tria.end(); ++cell)
#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_tools.h>
#include <deal.II/grid/grid_out.h>
{
const unsigned int dim=2;
Point<dim> origin;
- HyperShellBoundary<dim> boundary(origin);
+ SphericalManifold<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);
#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_tools.h>
#include <deal.II/grid/grid_out.h>
{
const unsigned int dim=2;
Point<dim> origin;
- HyperShellBoundary<dim> boundary(origin);
+ SphericalManifold<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);
Triangulation<d> tr;
GridGenerator::hyper_cube_with_cylindrical_hole(tr,inner,outer,outer,true);
- static HyperShellBoundary<d> boundary;
+ static SphericalManifold<d> boundary;
tr.set_all_manifold_ids_on_boundary(1,1);
tr.set_manifold(1,boundary);
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.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>
center, inner_radius, outer_radius,
10);
- static const HyperShellBoundary<2> boundary_description(center);
- triangulation.set_boundary (0, boundary_description);
+ static const SphericalManifold<2> boundary_description(center);
+ triangulation.set_manifold (0, boundary_description);
for (unsigned int step=0; step<5; ++step)
{
#include <deal.II/grid/grid_generator.h>
#include <deal.II/distributed/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>
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if (cell->at_boundary(f))
cell->face(f)->set_all_manifold_ids (numbers::invalid_manifold_id);
- // static const HyperShellBoundary<dim> boundary_shell;
- // tr.set_boundary (0, boundary_shell);
- // tr.set_boundary (1, boundary_shell);
+ // static const SphericalManifold<dim> boundary_shell;
+ // tr.set_manifold (0, boundary_shell);
+ // tr.set_manifold (1, boundary_shell);
tr.refine_global ((dim==2)?3:1);
#include <deal.II/lac/block_sparsity_pattern.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_tools.h>
#include <deal.II/dofs/dof_renumbering.h>
MPI_Comm mpi_communicator;
- HyperShellBoundary<dim> boundary;
+ SphericalManifold<dim> boundary;
parallel::distributed::Triangulation<dim> triangulation;
FESystem<dim> fe;
DoFHandler<dim> dof_handler;
#endif
- triangulation.set_boundary(0, boundary);
- triangulation.set_boundary(1, boundary);
+ triangulation.set_manifold(0, boundary);
+ triangulation.set_manifold(1, boundary);
triangulation.refine_global (4-dim);
#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 SphericalManifold<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_parallel_block_vector.h>
R1,
12,
true);
- static HyperShellBoundary<dim> boundary;
-// tr.set_boundary (0, boundary);
-//tr.set_boundary (1, boundary);
+ static SphericalManifold<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/lac/block_sparsity_pattern.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_tools.h>
#include <deal.II/dofs/dof_renumbering.h>
MPI_Comm mpi_communicator;
- HyperShellBoundary<dim> boundary;
+ SphericalManifold<dim> boundary;
parallel::distributed::Triangulation<dim> triangulation;
FESystem<dim> fe;
DoFHandler<dim> dof_handler;
#endif
- triangulation.set_boundary(0, boundary);
- triangulation.set_boundary(1, boundary);
+ triangulation.set_manifold(0, boundary);
+ triangulation.set_manifold(1, boundary);
triangulation.refine_global (4-dim);
#include <deal.II/distributed/grid_refinement.h>
#include <deal.II/grid/grid_refinement.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/grid/grid_tools.h>
EquationData::R1,
12,
true);
- static HyperShellBoundary<2> boundary;
- tr.set_boundary (0, boundary);
- tr.set_boundary (1, boundary);
+ static SphericalManifold<2> boundary;
+ tr.set_manifold (0, boundary);
+ tr.set_manifold (1, boundary);
tr.refine_global (initial_refinement);
#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/dofs/dof_handler.h>
#include <deal.II/fe/fe_system.h>
0.5,
1.0);
- static HyperShellBoundary<dim> boundary((Point<dim>()));
- triangulation.set_boundary (0, boundary);
+ static SphericalManifold<dim> boundary((Point<dim>()));
+ triangulation.set_manifold (0, boundary);
// write out the mesh. may not be
// strictly necessary here, but is
-0.2236 -0.2236 0.3873 0 0
-0.2364 -0.1626 0.4095 0 0
--0.1723 -0.1732 0.4362 0 0
+-0.1729 -0.1729 0.4361 0 0
-0.1626 -0.2364 0.4095 0 0
-0.2236 -0.2236 0.3873 0 0
-0.1626 -0.2364 0.4095 0 0
--0.1723 -0.1732 0.4362 0 0
--0.1047 -0.1804 0.4544 0 0
+-0.1729 -0.1729 0.4361 0 0
+-0.1054 -0.1799 0.4544 0 0
-0.0987 -0.2451 0.4245 0 0
-0.1626 -0.2364 0.4095 0 0
-0.0987 -0.2451 0.4245 0 0
--0.1047 -0.1804 0.4544 0 0
--0.0351 -0.1841 0.4636 0 0
+-0.1054 -0.1799 0.4544 0 0
+-0.0354 -0.1835 0.4638 0 0
-0.0331 -0.2495 0.4321 0 0
-0.0987 -0.2451 0.4245 0 0
-0.0331 -0.2495 0.4321 0 0
--0.0351 -0.1841 0.4636 0 0
-0.0351 -0.1841 0.4636 0 0
+-0.0354 -0.1835 0.4638 0 0
+0.0354 -0.1835 0.4638 0 0
0.0331 -0.2495 0.4321 0 0
-0.0331 -0.2495 0.4321 0 0
0.0331 -0.2495 0.4321 0 0
-0.0351 -0.1841 0.4636 0 0
-0.1047 -0.1804 0.4544 0 0
+0.0354 -0.1835 0.4638 0 0
+0.1054 -0.1799 0.4544 0 0
0.0987 -0.2451 0.4245 0 0
0.0331 -0.2495 0.4321 0 0
0.0987 -0.2451 0.4245 0 0
-0.1047 -0.1804 0.4544 0 0
-0.1723 -0.1732 0.4362 0 0
+0.1054 -0.1799 0.4544 0 0
+0.1729 -0.1729 0.4361 0 0
0.1626 -0.2364 0.4095 0 0
0.0987 -0.2451 0.4245 0 0
0.1626 -0.2364 0.4095 0 0
-0.1723 -0.1732 0.4362 0 0
+0.1729 -0.1729 0.4361 0 0
0.2364 -0.1626 0.4095 0 0
0.2236 -0.2236 0.3873 0 0
0.1626 -0.2364 0.4095 0 0
-0.2364 -0.1626 0.4095 0 0
-0.2451 -0.0987 0.4245 0 0
--0.1789 -0.1057 0.4548 0 0
--0.1723 -0.1732 0.4362 0 0
+-0.1799 -0.1054 0.4544 0 0
+-0.1729 -0.1729 0.4361 0 0
-0.2364 -0.1626 0.4095 0 0
--0.1723 -0.1732 0.4362 0 0
--0.1789 -0.1057 0.4548 0 0
--0.1089 -0.1104 0.4753 0 0
--0.1047 -0.1804 0.4544 0 0
--0.1723 -0.1732 0.4362 0 0
+-0.1729 -0.1729 0.4361 0 0
+-0.1799 -0.1054 0.4544 0 0
+-0.1099 -0.1099 0.4752 0 0
+-0.1054 -0.1799 0.4544 0 0
+-0.1729 -0.1729 0.4361 0 0
--0.1047 -0.1804 0.4544 0 0
--0.1089 -0.1104 0.4753 0 0
--0.0366 -0.1128 0.4857 0 0
--0.0351 -0.1841 0.4636 0 0
--0.1047 -0.1804 0.4544 0 0
+-0.1054 -0.1799 0.4544 0 0
+-0.1099 -0.1099 0.4752 0 0
+-0.0370 -0.1122 0.4858 0 0
+-0.0354 -0.1835 0.4638 0 0
+-0.1054 -0.1799 0.4544 0 0
--0.0351 -0.1841 0.4636 0 0
--0.0366 -0.1128 0.4857 0 0
-0.0366 -0.1128 0.4857 0 0
-0.0351 -0.1841 0.4636 0 0
--0.0351 -0.1841 0.4636 0 0
+-0.0354 -0.1835 0.4638 0 0
+-0.0370 -0.1122 0.4858 0 0
+0.0370 -0.1122 0.4858 0 0
+0.0354 -0.1835 0.4638 0 0
+-0.0354 -0.1835 0.4638 0 0
-0.0351 -0.1841 0.4636 0 0
-0.0366 -0.1128 0.4857 0 0
-0.1089 -0.1104 0.4753 0 0
-0.1047 -0.1804 0.4544 0 0
-0.0351 -0.1841 0.4636 0 0
+0.0354 -0.1835 0.4638 0 0
+0.0370 -0.1122 0.4858 0 0
+0.1099 -0.1099 0.4752 0 0
+0.1054 -0.1799 0.4544 0 0
+0.0354 -0.1835 0.4638 0 0
-0.1047 -0.1804 0.4544 0 0
-0.1089 -0.1104 0.4753 0 0
-0.1789 -0.1057 0.4548 0 0
-0.1723 -0.1732 0.4362 0 0
-0.1047 -0.1804 0.4544 0 0
+0.1054 -0.1799 0.4544 0 0
+0.1099 -0.1099 0.4752 0 0
+0.1799 -0.1054 0.4544 0 0
+0.1729 -0.1729 0.4361 0 0
+0.1054 -0.1799 0.4544 0 0
-0.1723 -0.1732 0.4362 0 0
-0.1789 -0.1057 0.4548 0 0
+0.1729 -0.1729 0.4361 0 0
+0.1799 -0.1054 0.4544 0 0
0.2451 -0.0987 0.4245 0 0
0.2364 -0.1626 0.4095 0 0
-0.1723 -0.1732 0.4362 0 0
+0.1729 -0.1729 0.4361 0 0
-0.2451 -0.0987 0.4245 0 0
-0.2495 -0.0331 0.4321 0 0
--0.1823 -0.0355 0.4642 0 0
--0.1789 -0.1057 0.4548 0 0
+-0.1835 -0.0354 0.4638 0 0
+-0.1799 -0.1054 0.4544 0 0
-0.2451 -0.0987 0.4245 0 0
--0.1789 -0.1057 0.4548 0 0
--0.1823 -0.0355 0.4642 0 0
--0.1110 -0.0372 0.4861 0 0
--0.1089 -0.1104 0.4753 0 0
--0.1789 -0.1057 0.4548 0 0
+-0.1799 -0.1054 0.4544 0 0
+-0.1835 -0.0354 0.4638 0 0
+-0.1122 -0.0370 0.4858 0 0
+-0.1099 -0.1099 0.4752 0 0
+-0.1799 -0.1054 0.4544 0 0
--0.1089 -0.1104 0.4753 0 0
--0.1110 -0.0372 0.4861 0 0
--0.0373 -0.0380 0.4972 0 0
--0.0366 -0.1128 0.4857 0 0
--0.1089 -0.1104 0.4753 0 0
+-0.1099 -0.1099 0.4752 0 0
+-0.1122 -0.0370 0.4858 0 0
+-0.0378 -0.0378 0.4971 0 0
+-0.0370 -0.1122 0.4858 0 0
+-0.1099 -0.1099 0.4752 0 0
--0.0366 -0.1128 0.4857 0 0
--0.0373 -0.0380 0.4972 0 0
-0.0373 -0.0380 0.4972 0 0
-0.0366 -0.1128 0.4857 0 0
--0.0366 -0.1128 0.4857 0 0
+-0.0370 -0.1122 0.4858 0 0
+-0.0378 -0.0378 0.4971 0 0
+0.0378 -0.0378 0.4971 0 0
+0.0370 -0.1122 0.4858 0 0
+-0.0370 -0.1122 0.4858 0 0
-0.0366 -0.1128 0.4857 0 0
-0.0373 -0.0380 0.4972 0 0
-0.1110 -0.0372 0.4861 0 0
-0.1089 -0.1104 0.4753 0 0
-0.0366 -0.1128 0.4857 0 0
+0.0370 -0.1122 0.4858 0 0
+0.0378 -0.0378 0.4971 0 0
+0.1122 -0.0370 0.4858 0 0
+0.1099 -0.1099 0.4752 0 0
+0.0370 -0.1122 0.4858 0 0
-0.1089 -0.1104 0.4753 0 0
-0.1110 -0.0372 0.4861 0 0
-0.1823 -0.0355 0.4642 0 0
-0.1789 -0.1057 0.4548 0 0
-0.1089 -0.1104 0.4753 0 0
+0.1099 -0.1099 0.4752 0 0
+0.1122 -0.0370 0.4858 0 0
+0.1835 -0.0354 0.4638 0 0
+0.1799 -0.1054 0.4544 0 0
+0.1099 -0.1099 0.4752 0 0
-0.1789 -0.1057 0.4548 0 0
-0.1823 -0.0355 0.4642 0 0
+0.1799 -0.1054 0.4544 0 0
+0.1835 -0.0354 0.4638 0 0
0.2495 -0.0331 0.4321 0 0
0.2451 -0.0987 0.4245 0 0
-0.1789 -0.1057 0.4548 0 0
+0.1799 -0.1054 0.4544 0 0
-0.2495 -0.0331 0.4321 0 0
-0.2495 0.0331 0.4321 0 0
--0.1823 0.0355 0.4642 0 0
--0.1823 -0.0355 0.4642 0 0
+-0.1835 0.0354 0.4638 0 0
+-0.1835 -0.0354 0.4638 0 0
-0.2495 -0.0331 0.4321 0 0
--0.1823 -0.0355 0.4642 0 0
--0.1823 0.0355 0.4642 0 0
--0.1110 0.0372 0.4861 0 0
--0.1110 -0.0372 0.4861 0 0
--0.1823 -0.0355 0.4642 0 0
+-0.1835 -0.0354 0.4638 0 0
+-0.1835 0.0354 0.4638 0 0
+-0.1122 0.0370 0.4858 0 0
+-0.1122 -0.0370 0.4858 0 0
+-0.1835 -0.0354 0.4638 0 0
--0.1110 -0.0372 0.4861 0 0
--0.1110 0.0372 0.4861 0 0
--0.0373 0.0380 0.4972 0 0
--0.0373 -0.0380 0.4972 0 0
--0.1110 -0.0372 0.4861 0 0
+-0.1122 -0.0370 0.4858 0 0
+-0.1122 0.0370 0.4858 0 0
+-0.0378 0.0378 0.4971 0 0
+-0.0378 -0.0378 0.4971 0 0
+-0.1122 -0.0370 0.4858 0 0
--0.0373 -0.0380 0.4972 0 0
--0.0373 0.0380 0.4972 0 0
-0.0373 0.0380 0.4972 0 0
-0.0373 -0.0380 0.4972 0 0
--0.0373 -0.0380 0.4972 0 0
+-0.0378 -0.0378 0.4971 0 0
+-0.0378 0.0378 0.4971 0 0
+0.0378 0.0378 0.4971 0 0
+0.0378 -0.0378 0.4971 0 0
+-0.0378 -0.0378 0.4971 0 0
-0.0373 -0.0380 0.4972 0 0
-0.0373 0.0380 0.4972 0 0
-0.1110 0.0372 0.4861 0 0
-0.1110 -0.0372 0.4861 0 0
-0.0373 -0.0380 0.4972 0 0
+0.0378 -0.0378 0.4971 0 0
+0.0378 0.0378 0.4971 0 0
+0.1122 0.0370 0.4858 0 0
+0.1122 -0.0370 0.4858 0 0
+0.0378 -0.0378 0.4971 0 0
-0.1110 -0.0372 0.4861 0 0
-0.1110 0.0372 0.4861 0 0
-0.1823 0.0355 0.4642 0 0
-0.1823 -0.0355 0.4642 0 0
-0.1110 -0.0372 0.4861 0 0
+0.1122 -0.0370 0.4858 0 0
+0.1122 0.0370 0.4858 0 0
+0.1835 0.0354 0.4638 0 0
+0.1835 -0.0354 0.4638 0 0
+0.1122 -0.0370 0.4858 0 0
-0.1823 -0.0355 0.4642 0 0
-0.1823 0.0355 0.4642 0 0
+0.1835 -0.0354 0.4638 0 0
+0.1835 0.0354 0.4638 0 0
0.2495 0.0331 0.4321 0 0
0.2495 -0.0331 0.4321 0 0
-0.1823 -0.0355 0.4642 0 0
+0.1835 -0.0354 0.4638 0 0
-0.2495 0.0331 0.4321 0 0
-0.2451 0.0987 0.4245 0 0
--0.1789 0.1057 0.4548 0 0
--0.1823 0.0355 0.4642 0 0
+-0.1799 0.1054 0.4544 0 0
+-0.1835 0.0354 0.4638 0 0
-0.2495 0.0331 0.4321 0 0
--0.1823 0.0355 0.4642 0 0
--0.1789 0.1057 0.4548 0 0
--0.1089 0.1104 0.4753 0 0
--0.1110 0.0372 0.4861 0 0
--0.1823 0.0355 0.4642 0 0
+-0.1835 0.0354 0.4638 0 0
+-0.1799 0.1054 0.4544 0 0
+-0.1099 0.1099 0.4752 0 0
+-0.1122 0.0370 0.4858 0 0
+-0.1835 0.0354 0.4638 0 0
--0.1110 0.0372 0.4861 0 0
--0.1089 0.1104 0.4753 0 0
--0.0366 0.1128 0.4857 0 0
--0.0373 0.0380 0.4972 0 0
--0.1110 0.0372 0.4861 0 0
+-0.1122 0.0370 0.4858 0 0
+-0.1099 0.1099 0.4752 0 0
+-0.0370 0.1122 0.4858 0 0
+-0.0378 0.0378 0.4971 0 0
+-0.1122 0.0370 0.4858 0 0
--0.0373 0.0380 0.4972 0 0
--0.0366 0.1128 0.4857 0 0
-0.0366 0.1128 0.4857 0 0
-0.0373 0.0380 0.4972 0 0
--0.0373 0.0380 0.4972 0 0
+-0.0378 0.0378 0.4971 0 0
+-0.0370 0.1122 0.4858 0 0
+0.0370 0.1122 0.4858 0 0
+0.0378 0.0378 0.4971 0 0
+-0.0378 0.0378 0.4971 0 0
-0.0373 0.0380 0.4972 0 0
-0.0366 0.1128 0.4857 0 0
-0.1089 0.1104 0.4753 0 0
-0.1110 0.0372 0.4861 0 0
-0.0373 0.0380 0.4972 0 0
+0.0378 0.0378 0.4971 0 0
+0.0370 0.1122 0.4858 0 0
+0.1099 0.1099 0.4752 0 0
+0.1122 0.0370 0.4858 0 0
+0.0378 0.0378 0.4971 0 0
-0.1110 0.0372 0.4861 0 0
-0.1089 0.1104 0.4753 0 0
-0.1789 0.1057 0.4548 0 0
-0.1823 0.0355 0.4642 0 0
-0.1110 0.0372 0.4861 0 0
+0.1122 0.0370 0.4858 0 0
+0.1099 0.1099 0.4752 0 0
+0.1799 0.1054 0.4544 0 0
+0.1835 0.0354 0.4638 0 0
+0.1122 0.0370 0.4858 0 0
-0.1823 0.0355 0.4642 0 0
-0.1789 0.1057 0.4548 0 0
+0.1835 0.0354 0.4638 0 0
+0.1799 0.1054 0.4544 0 0
0.2451 0.0987 0.4245 0 0
0.2495 0.0331 0.4321 0 0
-0.1823 0.0355 0.4642 0 0
+0.1835 0.0354 0.4638 0 0
-0.2451 0.0987 0.4245 0 0
-0.2364 0.1626 0.4095 0 0
--0.1723 0.1732 0.4362 0 0
--0.1789 0.1057 0.4548 0 0
+-0.1729 0.1729 0.4361 0 0
+-0.1799 0.1054 0.4544 0 0
-0.2451 0.0987 0.4245 0 0
--0.1789 0.1057 0.4548 0 0
--0.1723 0.1732 0.4362 0 0
--0.1047 0.1804 0.4544 0 0
--0.1089 0.1104 0.4753 0 0
--0.1789 0.1057 0.4548 0 0
+-0.1799 0.1054 0.4544 0 0
+-0.1729 0.1729 0.4361 0 0
+-0.1054 0.1799 0.4544 0 0
+-0.1099 0.1099 0.4752 0 0
+-0.1799 0.1054 0.4544 0 0
--0.1089 0.1104 0.4753 0 0
--0.1047 0.1804 0.4544 0 0
--0.0351 0.1841 0.4636 0 0
--0.0366 0.1128 0.4857 0 0
--0.1089 0.1104 0.4753 0 0
+-0.1099 0.1099 0.4752 0 0
+-0.1054 0.1799 0.4544 0 0
+-0.0354 0.1835 0.4638 0 0
+-0.0370 0.1122 0.4858 0 0
+-0.1099 0.1099 0.4752 0 0
--0.0366 0.1128 0.4857 0 0
--0.0351 0.1841 0.4636 0 0
-0.0351 0.1841 0.4636 0 0
-0.0366 0.1128 0.4857 0 0
--0.0366 0.1128 0.4857 0 0
+-0.0370 0.1122 0.4858 0 0
+-0.0354 0.1835 0.4638 0 0
+0.0354 0.1835 0.4638 0 0
+0.0370 0.1122 0.4858 0 0
+-0.0370 0.1122 0.4858 0 0
-0.0366 0.1128 0.4857 0 0
-0.0351 0.1841 0.4636 0 0
-0.1047 0.1804 0.4544 0 0
-0.1089 0.1104 0.4753 0 0
-0.0366 0.1128 0.4857 0 0
+0.0370 0.1122 0.4858 0 0
+0.0354 0.1835 0.4638 0 0
+0.1054 0.1799 0.4544 0 0
+0.1099 0.1099 0.4752 0 0
+0.0370 0.1122 0.4858 0 0
-0.1089 0.1104 0.4753 0 0
-0.1047 0.1804 0.4544 0 0
-0.1723 0.1732 0.4362 0 0
-0.1789 0.1057 0.4548 0 0
-0.1089 0.1104 0.4753 0 0
+0.1099 0.1099 0.4752 0 0
+0.1054 0.1799 0.4544 0 0
+0.1729 0.1729 0.4361 0 0
+0.1799 0.1054 0.4544 0 0
+0.1099 0.1099 0.4752 0 0
-0.1789 0.1057 0.4548 0 0
-0.1723 0.1732 0.4362 0 0
+0.1799 0.1054 0.4544 0 0
+0.1729 0.1729 0.4361 0 0
0.2364 0.1626 0.4095 0 0
0.2451 0.0987 0.4245 0 0
-0.1789 0.1057 0.4548 0 0
+0.1799 0.1054 0.4544 0 0
-0.2364 0.1626 0.4095 0 0
-0.2236 0.2236 0.3873 0 0
-0.1626 0.2364 0.4095 0 0
--0.1723 0.1732 0.4362 0 0
+-0.1729 0.1729 0.4361 0 0
-0.2364 0.1626 0.4095 0 0
--0.1723 0.1732 0.4362 0 0
+-0.1729 0.1729 0.4361 0 0
-0.1626 0.2364 0.4095 0 0
-0.0987 0.2451 0.4245 0 0
--0.1047 0.1804 0.4544 0 0
--0.1723 0.1732 0.4362 0 0
+-0.1054 0.1799 0.4544 0 0
+-0.1729 0.1729 0.4361 0 0
--0.1047 0.1804 0.4544 0 0
+-0.1054 0.1799 0.4544 0 0
-0.0987 0.2451 0.4245 0 0
-0.0331 0.2495 0.4321 0 0
--0.0351 0.1841 0.4636 0 0
--0.1047 0.1804 0.4544 0 0
+-0.0354 0.1835 0.4638 0 0
+-0.1054 0.1799 0.4544 0 0
--0.0351 0.1841 0.4636 0 0
+-0.0354 0.1835 0.4638 0 0
-0.0331 0.2495 0.4321 0 0
0.0331 0.2495 0.4321 0 0
-0.0351 0.1841 0.4636 0 0
--0.0351 0.1841 0.4636 0 0
+0.0354 0.1835 0.4638 0 0
+-0.0354 0.1835 0.4638 0 0
-0.0351 0.1841 0.4636 0 0
+0.0354 0.1835 0.4638 0 0
0.0331 0.2495 0.4321 0 0
0.0987 0.2451 0.4245 0 0
-0.1047 0.1804 0.4544 0 0
-0.0351 0.1841 0.4636 0 0
+0.1054 0.1799 0.4544 0 0
+0.0354 0.1835 0.4638 0 0
-0.1047 0.1804 0.4544 0 0
+0.1054 0.1799 0.4544 0 0
0.0987 0.2451 0.4245 0 0
0.1626 0.2364 0.4095 0 0
-0.1723 0.1732 0.4362 0 0
-0.1047 0.1804 0.4544 0 0
+0.1729 0.1729 0.4361 0 0
+0.1054 0.1799 0.4544 0 0
-0.1723 0.1732 0.4362 0 0
+0.1729 0.1729 0.4361 0 0
0.1626 0.2364 0.4095 0 0
0.2236 0.2236 0.3873 0 0
0.2364 0.1626 0.4095 0 0
-0.1723 0.1732 0.4362 0 0
+0.1729 0.1729 0.4361 0 0
-0.4472 -0.4472 0.7746 0 0
#include <deal.II/lac/constraint_matrix.h>
#include <deal.II/lac/block_sparsity_pattern.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/distributed/tria.h>
#include <deal.II/dofs/dof_renumbering.h>
#include <deal.II/dofs/dof_tools.h>
MPI_Comm mpi_communicator;
- HyperShellBoundary<dim> boundary;
+ SphericalManifold<dim> boundary;
parallel::distributed::Triangulation<dim> triangulation;
FESystem<dim> fe;
DoFHandler<dim> dof_handler;
0,
true);
- triangulation.set_boundary(0, boundary);
- triangulation.set_boundary(1, boundary);
+ triangulation.set_manifold(0, boundary);
+ triangulation.set_manifold(1, boundary);
setup_dofs ();