Thanks Wolfgang for your insight, and a quick note:
> > I am trying to merge two triangulations together.
> >
> > 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.
>
> 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.
I like this option. Say, after creating sub_triangulations:
GridTools::merge_coarse_cell_information
(std::vector<dealii::Triangulation> coarse_grids,
bool destroy_coarse_grids = false);
where the boolean could free up memory if the old triangualtions are not
needed. This should be false by default because the old triangulations can
(in principle) be used for postprocessing initiatives.
> 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
Ouch! I missed that function somehow; great that function exist
already.... There are a few others that may help here that I missed.
GridTools::shift(), GridTools::transform()... and my favourite,
GridTools::scale().
Thanks very much. I will see what I can do.
Best,
Toby
-----
Toby D. Young
Assistant Professor
Philosophy-Physics
Polish Academy of Sciences
Warszawa, Polska
www: http://www.ippt.gov.pl/~tyoung
skype: stenografia
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii