> I am trying to merge two triangulations together.
>
> Lets say I want to fit a hyper_cube () triangulation into the slit of a
> hyper_cube_slit () triangulation; which would give me the equivalent
> of a (different) hyper_cube () triangulation. How would I do this?
>
> Is there an existing way to join (merge) to triangulations into a single
> triangulation in the deal.II library? I can not find one...
>
> Can anyone suggest a relatively simple way to create grids that are
> also quite simple using the tools existing in deal.II?
I don't think this functionality exists yet. But it shouldn't be too hard to
do. You'd need the following functions:
1/ A function
GridTools::extract_coarse_cell_information (const Triangulation &, ...)
that returns a std::vector<CellData>, SubCellData, and std::vector<Point> for
all the nodes. This function would essentially be the opposite of
Triangulation::create_triangulation.
2/ You'd do this for your two triangulations separately. Then you'd need a
function, say
GridTools::merge_coarse_cell_information
that would merge these two data sets.
3/ This latter function would probably have to unify vertices that come from
the two different triangulations and that therefore have different numbers
but that are located at the same places. There is already a function that
does that:
GridTools::delete_duplicated_vertices
4/ After this step, you can pass the merged information back to
Triangulation::create_triangulation.
Let me know if you need help with any of these steps!
Best
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
