13#ifndef dealii_grid_generator_h
14#define dealii_grid_generator_h
102 template <
int dim,
int spacedim>
105 const double left = 0.,
106 const double right = 1.,
153 template <
int dim,
int spacedim>
188 template <
int dim,
int spacedim>
191 const unsigned int repetitions,
192 const double left = 0.,
193 const double right = 1.,
234 template <
int dim,
int spacedim>
294 template <
int dim,
int spacedim>
297 const std::vector<unsigned int> &repetitions,
324 const std::vector<std::vector<double>> &step_sizes,
346 const std::vector<std::vector<double>> &spacing,
380 template <
int dim,
int spacedim>
383 const std::vector<unsigned int> &holes);
410 template <
int dim,
int spacedim>
414 const double voxel_edge_length);
468 const double inner_radius = 0.4,
470 const double pad_bottom = 2.,
471 const double pad_top = 2.,
472 const double pad_left = 1.,
473 const double pad_right = 1.,
478 const unsigned int n_slices = 2,
576 const double shell_region_width = 0.03,
577 const unsigned int n_shells = 2,
578 const double skewness = 2.0,
670 const std::vector<unsigned int> &lengths_and_heights,
671 const double depth = 1,
672 const unsigned int depth_division = 1,
673 const double shell_region_radius = 0.75,
674 const unsigned int n_shells = 2,
675 const double skewness = 2.0,
676 const bool use_transfinite_region =
false,
749 const std::vector<unsigned int> &lengths_heights_widths,
750 const double inner_radius,
752 const bool use_transfinite_region =
false,
781 const double inner_radius = 0.25,
808 template <
int dim,
int spacedim>
881 const unsigned int n_subdivisions,
901 const unsigned int (&n_subdivisions)[dim],
903 const unsigned int *n_subdivisions,
931 template <
int dim,
int spacedim>
936 const std::vector<unsigned int> &subdivisions = {},
962 const double left = 0.,
963 const double right = 1.,
964 const double thickness = 1.,
1049 template <
int dim,
int spacedim>
1053 const double radius = 1.,
1054 const bool attach_spherical_manifold_on_boundary_cells =
false);
1095 const double radius = 1.);
1117 const unsigned int n_rotate_middle_square);
1142 const bool face_orientation,
1143 const bool face_flip,
1144 const bool face_rotation,
1145 const bool manipulate_left_cube);
1176 template <
int spacedim>
1180 const double radius = 1.);
1222 const double radius = 1.);
1243 const double radius = 1.);
1274 const double radius = 1.,
1275 const double half_length = 1.);
1318 const unsigned int x_subdivisions,
1319 const double radius = 1.,
1320 const double half_length = 1.);
1360 const double radius_0 = 1.0,
1361 const double radius_1 = 0.5,
1362 const double half_length = 1.0);
1487 template <
int dim,
int spacedim>
1520 template <
int dim,
int spacedim>
1523 const std::vector<unsigned int> &sizes,
1524 const bool colorize_cells =
false);
1564 const double left = -1.,
1565 const double right = 1.,
1599 template <
int dim,
int spacedim>
1602 const std::vector<unsigned int> &repetitions,
1605 const std::vector<int> &n_cells_to_remove);
1629 const double left = 0.,
1630 const double right = 1.,
1708 template <
int dim,
int spacedim>
1712 const double inner_radius,
1714 const unsigned int n_cells = 0,
1769 const double inner_radius,
1771 const unsigned int n_cells);
1806 const double inner_radius,
1808 const unsigned int n_cells = 0,
1841 const double inner_radius,
1843 const unsigned int n_cells = 0,
1880 const double length,
1881 const double inner_radius,
1883 const unsigned int n_radial_cells = 0,
1884 const unsigned int n_axial_cells = 0,
1938 template <
int dim,
int spacedim>
1941 const double centerline_radius,
1942 const double inner_radius,
1943 const unsigned int n_cells_toroidal = 6,
1979 template <
int dim,
int spacedim>
1982 const double inner_radius = .25,
1984 const double L = .5,
1985 const unsigned int repetitions = 1,
2062 const double inner_radius = 0.125,
2064 const unsigned int n_shells = 1,
2065 const double skewness = 0.1,
2066 const unsigned int n_cells_per_shell = 0,
2084 const unsigned int n_cells,
2085 const unsigned int n_rotations,
2124 template <
int dim,
int spacedim>
2128 const std::string &grid_generator_function_name,
2129 const std::string &grid_generator_function_arguments);
2201 const double &outer_ball_radius = 1.0);
2290 template <
int dim,
int spacedim>
2295 const double duplicated_vertex_tolerance = 1.0e-12,
2296 const bool copy_manifold_ids =
false,
2297 const bool copy_boundary_ids =
false);
2318 template <
int dim,
int spacedim>
2323 const double duplicated_vertex_tolerance = 1.0e-12,
2324 const bool copy_manifold_ids =
false,
2325 const bool copy_boundary_ids =
false);
2377 template <
int dim,
int spacedim = dim>
2380 const std::vector<unsigned int> &extents,
2420 template <
int dim,
int spacedim>
2463 template <
int dim,
int spacedim>
2544 const unsigned int n_slices,
2547 const bool copy_manifold_ids =
false,
2548 const std::vector<types::manifold_id> &manifold_priorities = {});
2559 const unsigned int n_slices,
2562 const bool copy_manifold_ids =
false,
2563 const std::vector<types::manifold_id> &manifold_priorities = {});
2587 const std::vector<double> &slice_coordinates,
2589 const bool copy_manifold_ids =
false,
2590 const std::vector<types::manifold_id> &manifold_priorities = {});
2601 const std::vector<double> &slice_coordinates,
2603 const bool copy_manifold_ids =
false,
2604 const std::vector<types::manifold_id> &manifold_priorities = {});
2647 template <
int dim,
int spacedim1,
int spacedim2>
2705 template <
int dim,
int spacedim>
2709 const unsigned int n_divisions = (dim == 2 ?
2778 template <
int dim,
int spacedim>
2805 template <
int dim,
int spacedim>
3013 template <
int dim,
int spacedim>
3017 const std::vector<unsigned int> &repetitions,
3039 template <
int dim,
int spacedim>
3042 const unsigned int repetitions,
3043 const double left = 0.0,
3044 const double right = 1.0,
3063 template <
template <
int,
int>
class MeshType,
int dim,
int spacedim>
3066 struct ExtractBoundaryMesh
3069 std::map<
typename MeshType<dim - 1, spacedim>::cell_iterator,
3070 typename MeshType<dim, spacedim>::face_iterator>;
3171 template <
template <
int,
int>
class MeshType,
int dim,
int spacedim>
3178 std::map<
typename MeshType<dim - 1, spacedim>::cell_iterator,
3179 typename MeshType<dim, spacedim>::face_iterator>
3181 typename ExtractBoundaryMesh<MeshType, dim, spacedim>::return_type
3185 MeshType<dim - 1, spacedim> &surface_mesh,
3186 const std::set<types::boundary_id> &boundary_ids =
3187 std::set<types::boundary_id>());
3207 <<
"The number of repetitions " << arg1 <<
" must be >=1.");
3213 <<
"The vector of repetitions must have " << arg1
3220 "The input to this function is oriented in a way that will"
3221 " cause all cells to have negative measure.");
3281 const std::vector<unsigned int> &,
3290 const std::vector<unsigned int> &,
3299 const std::vector<unsigned int> &,
* const unsigned int aspect_ratio
* * const Number outer_radius
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_CXX20_REQUIRES(condition)
#define DEAL_II_NAMESPACE_CLOSE
#define DeclException0(Exception0)
static ::ExceptionBase & ExcInvalidInputOrientation()
static ::ExceptionBase & ExcInvalidRepetitions(int arg1)
#define DeclExceptionMsg(Exception, defaulttext)
static ::ExceptionBase & ExcInvalidRadii()
#define DeclException1(Exception1, type1, outsequence)
static ::ExceptionBase & ExcInvalidRepetitionsDimension(int arg1)
std::vector< index_type > data
void create_triangulation(Triangulation< dim, dim > &tria, const AdditionalData &additional_data=AdditionalData())
void parallelepiped(Triangulation< dim > &tria, const Point< dim >(&corners)[dim], const bool colorize=false)
void hyper_cross(Triangulation< dim, spacedim > &tria, const std::vector< unsigned int > &sizes, const bool colorize_cells=false)
A center cell with stacks of cell protruding from each surface.
void hyper_cube_with_cylindrical_hole(Triangulation< dim, spacedim > &triangulation, const double inner_radius=.25, const double outer_radius=.5, const double L=.5, const unsigned int repetitions=1, const bool colorize=false)
void hyper_ball_balanced(Triangulation< dim > &tria, const Point< dim > ¢er=Point< dim >(), const double radius=1.)
void plate_with_a_hole(Triangulation< dim > &tria, const double inner_radius=0.4, const double outer_radius=1., const double pad_bottom=2., const double pad_top=2., const double pad_left=1., const double pad_right=1., const Point< dim > ¢er=Point< dim >(), const types::manifold_id polar_manifold_id=0, const types::manifold_id tfi_manifold_id=1, const double L=1., const unsigned int n_slices=2, const bool colorize=false)
Rectangular plate with an (offset) cylindrical hole.
void generate_from_name_and_arguments(Triangulation< dim, spacedim > &tria, const std::string &grid_generator_function_name, const std::string &grid_generator_function_arguments)
void enclosed_hyper_cube(Triangulation< dim > &tria, const double left=0., const double right=1., const double thickness=1., const bool colorize=false)
void uniform_channel_with_sphere(Triangulation< dim > &tria, const std::vector< unsigned int > &lengths_heights_widths, const double inner_radius, const double outer_radius, const bool use_transfinite_region=false, const bool colorize=false)
void replicate_triangulation(const Triangulation< dim, spacedim > &input, const std::vector< unsigned int > &extents, Triangulation< dim, spacedim > &result)
Replicate a given triangulation in multiple coordinate axes.
void parallelogram(Triangulation< dim > &tria, const Point< dim >(&corners)[dim], const bool colorize=false)
void general_cell(Triangulation< dim, spacedim > &tria, const std::vector< Point< spacedim > > &vertices, const bool colorize=false)
void subdivided_hyper_cube(Triangulation< dim, spacedim > &tria, const unsigned int repetitions, const double left=0., const double right=1., const bool colorize=false)
void hyper_shell(Triangulation< dim, spacedim > &tria, const Point< spacedim > ¢er, const double inner_radius, const double outer_radius, const unsigned int n_cells=0, bool colorize=false)
void hyper_L(Triangulation< dim > &tria, const double left=-1., const double right=1., const bool colorize=false)
void subdivided_hyper_cube_with_simplices(Triangulation< dim, spacedim > &tria, const unsigned int repetitions, const double left=0.0, const double right=1.0, const bool colorize=false)
void hyper_cube_slit(Triangulation< dim > &tria, const double left=0., const double right=1., const bool colorize=false)
void hyper_cube_with_spherical_hole(Triangulation< dim > &tria, const double inner_radius=0.25, const double outer_radius=0.5, const bool colorize=false)
void eccentric_hyper_shell(Triangulation< dim > &triangulation, const Point< dim > &inner_center, const Point< dim > &outer_center, const double inner_radius, const double outer_radius, const unsigned int n_cells)
void hyper_rectangle(Triangulation< dim, spacedim > &tria, const Point< dim > &p1, const Point< dim > &p2, const bool colorize=false)
void surface_mesh_to_volumetric_mesh(const Triangulation< 2, 3 > &surface_tria, Triangulation< 3 > &vol_tria, const CGALWrappers::AdditionalData< 3 > &data=CGALWrappers::AdditionalData< 3 >{})
void uniform_channel_with_cylinder(Triangulation< dim > &tria, const std::vector< unsigned int > &lengths_and_heights, const double depth=1, const unsigned int depth_division=1, const double shell_region_radius=0.75, const unsigned int n_shells=2, const double skewness=2.0, const bool use_transfinite_region=false, const bool colorize=false)
void cylinder(Triangulation< dim > &tria, const double radius=1., const double half_length=1.)
void moebius(Triangulation< 3, 3 > &tria, const unsigned int n_cells, const unsigned int n_rotations, const double R, const double r)
void pipe_junction(Triangulation< dim, spacedim > &tria, const std::vector< std::pair< Point< spacedim >, double > > &openings, const std::pair< Point< spacedim >, double > &bifurcation, const double aspect_ratio=0.5)
void extrude_triangulation(const Triangulation< 2, 2 > &input, const unsigned int n_slices, const double height, Triangulation< 3, 3 > &result, const bool copy_manifold_ids=false, const std::vector< types::manifold_id > &manifold_priorities={})
void half_hyper_shell(Triangulation< dim > &tria, const Point< dim > ¢er, const double inner_radius, const double outer_radius, const unsigned int n_cells=0, const bool colorize=false)
void quarter_hyper_ball(Triangulation< dim > &tria, const Point< dim > ¢er=Point< dim >(), const double radius=1.)
void cheese(Triangulation< dim, spacedim > &tria, const std::vector< unsigned int > &holes)
Rectangular domain with rectangular pattern of holes.
void implicit_function(Triangulation< dim, 3 > &tria, const Function< 3 > &implicit_function, const CGALWrappers::AdditionalData< dim > &data=CGALWrappers::AdditionalData< dim >{}, const Point< 3 > &interior_point=Point< 3 >(), const double &outer_ball_radius=1.0)
void create_union_triangulation(const Triangulation< dim, spacedim > &triangulation_1, const Triangulation< dim, spacedim > &triangulation_2, Triangulation< dim, spacedim > &result)
void hyper_ball(Triangulation< dim, spacedim > &tria, const Point< spacedim > ¢er={}, const double radius=1., const bool attach_spherical_manifold_on_boundary_cells=false)
void subdivided_hyper_rectangle_with_simplices(Triangulation< dim, spacedim > &tria, const std::vector< unsigned int > &repetitions, const Point< dim > &p1, const Point< dim > &p2, const bool colorize=false)
void non_standard_orientation_mesh(Triangulation< 2 > &tria, const unsigned int n_rotate_middle_square)
return_type extract_boundary_mesh(const MeshType< dim, spacedim > &volume_mesh, MeshType< dim - 1, spacedim > &surface_mesh, const std::set< types::boundary_id > &boundary_ids=std::set< types::boundary_id >())
void alfeld_split_of_simplex_mesh(const Triangulation< dim, spacedim > &in_tria, Triangulation< dim, spacedim > &out_tria)
void subdivided_parallelepiped(Triangulation< dim > &tria, const unsigned int n_subdivisions, const Point< dim >(&corners)[dim], const bool colorize=false)
void subdivided_cylinder(Triangulation< dim > &tria, const unsigned int x_subdivisions, const double radius=1., const double half_length=1.)
void convert_hypercube_to_simplex_mesh(const Triangulation< dim, spacedim > &in_tria, Triangulation< dim, spacedim > &out_tria, const unsigned int n_divisions=(dim==2 ? 8u :24u))
void channel_with_cylinder(Triangulation< dim > &tria, const double shell_region_width=0.03, const unsigned int n_shells=2, const double skewness=2.0, const bool colorize=false)
void subdivided_hyper_L(Triangulation< dim, spacedim > &tria, const std::vector< unsigned int > &repetitions, const Point< dim > &bottom_left, const Point< dim > &top_right, const std::vector< int > &n_cells_to_remove)
void convert_simplex_to_hypercube_mesh(const Triangulation< dim, spacedim > &in_tria, Triangulation< dim, spacedim > &out_tria)
void hyper_sphere(Triangulation< spacedim - 1, spacedim > &tria, const Point< spacedim > ¢er=Point< spacedim >(), const double radius=1.)
void concentric_hyper_shells(Triangulation< dim > &triangulation, const Point< dim > ¢er, const double inner_radius=0.125, const double outer_radius=0.25, const unsigned int n_shells=1, const double skewness=0.1, const unsigned int n_cells_per_shell=0, const bool colorize=false)
void reference_cell(Triangulation< dim, spacedim > &tria, const ReferenceCell< dim > &reference_cell)
void subdivided_hyper_rectangle(Triangulation< dim, spacedim > &tria, const std::vector< unsigned int > &repetitions, const Point< dim > &p1, const Point< dim > &p2, const bool colorize=false)
void quarter_hyper_shell(Triangulation< dim > &tria, const Point< dim > ¢er, const double inner_radius, const double outer_radius, const unsigned int n_cells=0, const bool colorize=false)
void hyper_cube(Triangulation< dim, spacedim > &tria, const double left=0., const double right=1., const bool colorize=false)
void torus(Triangulation< dim, spacedim > &tria, const double centerline_radius, const double inner_radius, const unsigned int n_cells_toroidal=6, const double phi=2.0 *numbers::PI)
void create_triangulation_with_removed_cells(const Triangulation< dim, spacedim > &input_triangulation, const std::set< typename Triangulation< dim, spacedim >::active_cell_iterator > &cells_to_remove, Triangulation< dim, spacedim > &result)
void simplex(Triangulation< dim, dim > &tria, const std::vector< Point< dim > > &vertices)
void truncated_cone(Triangulation< dim > &tria, const double radius_0=1.0, const double radius_1=0.5, const double half_length=1.0)
void merge_triangulations(const Triangulation< dim, spacedim > &triangulation_1, const Triangulation< dim, spacedim > &triangulation_2, Triangulation< dim, spacedim > &result, const double duplicated_vertex_tolerance=1.0e-12, const bool copy_manifold_ids=false, const bool copy_boundary_ids=false)
void half_hyper_ball(Triangulation< dim > &tria, const Point< dim > ¢er=Point< dim >(), const double radius=1.)
void cylinder_shell(Triangulation< dim > &tria, const double length, const double inner_radius, const double outer_radius, const unsigned int n_radial_cells=0, const unsigned int n_axial_cells=0, const bool colorize=false)
void flatten_triangulation(const Triangulation< dim, spacedim1 > &in_tria, Triangulation< dim, spacedim2 > &out_tria)
unsigned int n_subdivision_x_0
unsigned int n_subdivision_y
unsigned int n_subdivision_x_2
void add_parameters(ParameterHandler &prm)
unsigned int airfoil_sampling_factor
unsigned int n_subdivision_x_1
Point< 2, double > joukowski_center