Thanks a lot Martin,
I had been waiting for this type of solution in dealii for a long time!
I also read in another thread about using MappingQGeneric instead of
MappingQ, so I updated my test code as below, and everything worked like a
charm!
PolarManifold<dim> boundary;
TransfiniteInterpolationManifold<dim> inner_manifold;
Triangulation<dim> triangulation;
Geom_parameters gp;
unsigned int curved_faces_label=10, not_curved = 1;
cell_rod ( triangulation, 1.0, R, gp, false );
triangulation.set_manifold ( curved_faces_label,
boundary ); // labels defined with the grid
inner_manifold.initialize(triangulation);
triangulation.set_manifold (not_curved, inner_manifold);
const MappingQGeneric<dim> mapping (degree);
const QGauss<dim> quadrature(degree);
// const FE_Q<dim> dummy_fe (1);
const FE_Q<dim> dummy_fe (QGaussLobatto<1>(degree + 1));
DoFHandler<dim> dof_handler (triangulation);
FEValues<dim> fe_values (mapping, dummy_fe, quadrature,
update_JxW_values);
Juan Carlos Araújo,
Umeå Universitet
El sábado, 3 de junio de 2017, 10:36:58 (UTC+2), Martin Kronbichler
escribió:
>
> Dear Juan Carlos,
>
> We have now finalized the TransfiniteInterpolationManifold implementation
> in deal.II,
>
> http://dealii.org/developer/doxygen/deal.II/classTransfiniteInterpolationManifold.html
>
> If you are interested, please try it out and let us know. As its
> implementation is completely new and I have only applied it to two
> applications in my group, there might still be issues, so let us know.
>
> By the way, if you are working on a 2D disk I recommend using
> PolarManifold rather than SphericalManifold. I have not yet looked in your
> code, but I recommend to try the transfinite interpolation manifold first.
>
> Best,
> Martin
>
--
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].
For more options, visit https://groups.google.com/d/optout.
