To learn it, I changed grid_3() in step-49 a little bit as follows
void grid_3 ()
{
Triangulation<2> triangulation;
GridGenerator::hyper_cube_with_cylindrical_hole (triangulation, 0.25, 1.0);
Triangulation<2>::active_cell_iterator
cell = triangulation.begin_active(),
endc = triangulation.end();
for (; cell!=endc; ++cell)
{
for (unsigned int i=0; i<GeometryInfo<2>::vertices_per_cell; ++i)
{
Point<2> &v = cell->vertex(i);
if (std::abs(v(1)-1.0)<1e-5)
{
v(1) -= 0.5;
}
}
}
GridTools::copy_boundary_to_manifold_id(triangulation);
const SphericalManifold<2> boundary_description(Point<2>(0,0));
triangulation.set_manifold (1, boundary_description);
triangulation.refine_global(2);
std::vector<bool> vertex_touched (triangulation.n_vertices(),false);
std::map< unsigned int, Point< 2 > > new_points;
cell = triangulation.begin_active(),
endc = triangulation.end();
for (; cell!=endc; ++cell)
{
for(unsigned int v=0;v < GeometryInfo<2>::vertices_per_cell;++v)
{
if(vertex_touched[cell->vertex_index(v)] == false)
{
vertex_touched[cell->vertex_index(v)] = true;
unsigned int vertex_number = cell->vertex_index(v);
new_points[vertex_number]= cell->vertex(v);
}
}
}
GridTools::laplace_transform (new_points,triangulation);
But, as Denis said, there was nothing special compared to original result of
step-49 (FYI I'm using 8.3.0 version)
I'm not quite sure what you're getting at here. The mesh you show looks
correct to me. What is it that you're surprised about?
you said that it is added to the testsuite.
But, what is the testsuite...?
It's here:
https://www.dealii.org/developer/developers/testsuite.html
It is, in essence, a collection of several thousand small programs that we run
for each change we apply to the library to make sure that what used to work in
the past continues to work in the future.
where can I get those things and learn how to apply those changes in my code...?
That would be here:
https://www.dealii.org/participate.html
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
