Effect of attaching SphericalManifold to all entities

[Praveen C]

If I do the following

(dim = 2)

GridGenerator::hyper_shell (triangulation,
Point<dim>(0,0),
1.0,
2.0);
static const SphericalManifold<dim> boundary;
triangulation.set_all_manifold_ids(0);
triangulation.set_manifold (0, boundary);

and use MappingQ<dim>(2), does this make all faces in the mesh to be curved (except the radially oriented faces) ?

Thanks
praveen

[Martin Kronbichler]

Hi Praveen,

no, MappingQ does not unless you set the second boolean flag in the
constructor to enable the mapping on all cells and not just the
boundary. Even worse, when you have attached a curved manifold on the
cells and not just the boundary you will have gaps in your computational
domain. I highly recommend using MappingQGeneric instead.

Best,
Martin

P.S. We had a discussion about removing the class MappingQ and just use
MappingQGeneric always just the other day:
https://github.com/dealii/dealii/issues/3874#issuecomment-305513405 (and
subsequent comments).

[Praveen C]

Thanks a lot Martin. This was very useful as I was not aware of the differences b/w MappingQ and MappingQGeneric. I get expected results with MappingQGeneric now.

Best regards

praveen

--
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 dealii+unsubscribe@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.