#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/grid/grid_out.h>
// Helper function
for(unsigned int i=0; i<dim; ++i)
center[i] = .5;
- HyperBallBoundary<dim,spacedim> boundary(center,center.norm());
+ PolarManifold<spacedim> boundary(center);
Triangulation<dim,spacedim> tria;
GridGenerator::hyper_cube (tria);
typename Triangulation<dim,spacedim>::active_cell_iterator cell;
tria.begin_active()->face(0)->set_manifold_id(1);
- tria.set_boundary(1,boundary);
+ tria.set_manifold(1,boundary);
tria.refine_global(2);
#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/grid/grid_out.h>
// Helper function
double radius=center.norm();
- HyperBallBoundary<dim,spacedim> boundary(center,radius);
+ PolarManifold<spacedim> boundary(center);
Triangulation<dim,spacedim> tria;
GridGenerator::hyper_cube (tria);
typename Triangulation<dim,spacedim>::active_cell_iterator cell;
tria.refine_global(1);
-
for(cell=tria.begin_active(); cell!=tria.end(); ++cell)
if(dim<spacedim && cell->center().distance(center)<radius)
if(cell->face(f)->center().distance(center)< radius)
cell->face(f)->set_all_manifold_ids(1);
- tria.set_boundary(1,boundary);
+ tria.set_manifold(1,boundary);
tria.refine_global(2);
GridOut gridout;
#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/grid/grid_out.h>
// Helper function
center[i] = .5;
Triangulation<dim,spacedim> tria;
- HyperBallBoundary<dim,spacedim> boundary(center,center.norm());
+ PolarManifold<spacedim> boundary(center);
GridGenerator::hyper_cube (tria);
typename Triangulation<dim,spacedim>::active_cell_iterator cell;
tria.begin_active()->face(0)->set_manifold_id(1);
- tria.set_boundary(1,boundary);
+ tria.set_manifold(1,boundary);
tria.refine_global(2);
dim, spacedim);
std::ofstream ofile(fname);
gridout.write_msh(tria,ofile);
- tria.set_boundary(1);
+ tria.set_manifold(1);
}
int main ()