> 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
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);
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