> CylindricalManifold<2> cylman(0);
> triangulation.set_manifold(numbers::flat_manifold_id, cylman);
The flat_manifold_id is usually reserved for the flat manifold, and yours is
not a flat manifold. Setting the flat manifold to the cylindrical manifold
implies
You should set the manifold id
Dear Luca,
Thank you very much for your help. I tried to use your recommendation, but
it does not seem to work correctly. My code is:
CylindricalManifold<2> cylman(0);
triangulation.set_manifold(numbers::flat_manifold_id, cylman);
QGauss<2> quadrature_formula(fe.degree + 1);
MappingManifold<2
You don’t need to do anything special. Just create a MappingManifold, and pass
it where required. The manifold that is used is the one associated to the
triangulation. You don’t need to associate it to the Mapping, since the mapping
is used on *all* objects of the triangulation, and it will know
Hi,
I want to use a real-to-reference-cell mapping which corresponds to
cylindrical symmetry in 2D. In other words, the integration jacobian and
the shape function gradients should contain the scale factor of the
cylindrical coordinates. I tried to do it with the following code:
QGauss<2> q
