Kyusik,
As you said above, I only touched the points on the boundary as
follows... (I attached the file for detail)
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)
{
if(cell->at_boundary()==true)
{
for(unsigned int face=0;face <
GeometryInfo<2>::faces_per_cell;++face)
{
if(cell->face(face)->at_boundary()==true)
{
for(unsigned int v=0;v <
GeometryInfo<2>::vertices_per_face;++v)
{
if(vertex_touched[cell->vertex_index(v)] == false)
You need to call
cell->face(face)->vertex_index(v)
here.
{
vertex_touched[cell->vertex_index(v)] = true;
And here as well.
I checked num_check value in the code is same as the number of points on
the boundary.
But, nothing happened...
As before, there is no kind of smoothness in the mesh.
There is no change in the mesh.
Could you please help me?
I agree that something looks broken since nodes are not moved even with
my fixed above. But I can't immediately see what exactly it is that's
wrong. There are a number of tests in deal.II/tests/grid/grid_transform*
that use this function and that you could take a look at to see how the
function is used.
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.
