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> quadrature_formula(fe.degree + 1); CylindricalManifold<2> cyl_man(0); MappingManifold<2> mapping; auto &data = dynamic_cast< MappingManifold<2>::InternalData &>(*mapping.get_data(update_values | update_gradients | update_JxW_values, quadrature_formula)); data.manifold = &cyl_man; but it does not seem to work. I am afraid that this is not the correct way to associate the Manifold object with the Mapping object, but I can't figure out how it should be done. Can anyone help? Best regards, Andreas -- 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]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/9d0b3e81-6559-4dc9-bc89-0f4a74ec4921%40googlegroups.com.
