OK, I see.
Anyway, as always, Thanks for educating me free of charge :-)
Kyusik.
--
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 Googl
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);
Dr. Bangerth,
Thanks for your reply again.
> Yes, but that's because you have a small bug in your code: you fix *all*
> vertices before you call laplace_transform. You should only fix the ones
> that
> are are at the boundary of the domain.
>
>
As you said above, I only touched the points on
Kyusik,
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?
I'm sorry, I should have been more specific.
What I mean that the result has no difference compared to when I didn't use
GridTools::lapace_transform(
Thank you very much for your quick reply.
>
> 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?
>
>
> I'm sorry, I should have been more specific.
What I mean that the result has no difference compared to when I didn
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 = triangula
