Re: [Getfem-users] Enumerer tout les DOF d'une face à l'ordre 2

[Renard Yves]

Dear Cedric,

The langage used on this mailing list is preferably english.

In order to enumerate the dof on a face you should use the method
mf.ind_basic_dof_of_face_of_element(cv, f)
(and mf.nb_basic_dof_of_face_of_element(cv, f) to have the number of dof).

the method cvs->nb_points_of_face(face_id) return the number of points  
of the reference element. It is independent of the finite element  
method.
The construction
m_shape_mesh->ind_dof_of_element(element_id)[cvs->ind_points_of_face(face_id)[i]];
is not valid because it mixes two different concepts : the vertex of  
the element and the degrees of freedom.

Yves.




Cédric Venet <address@hidden> a écrit :

> Bonjour,
>
> Je voudrais énumérer tout les DOF d'une face d'un element. J'utilise les
> fonctions:
>
> bgeot::pconvex_structure cvs = m_mesh->structure_of_convex(element_id);
> const int np = cvs->nb_points_of_face(face_id);
> getfem::mesh_fem::ind_dof_ct l2g =
> m_shape_mesh->ind_dof_of_element(element_id);
> l2g[cvs->ind_points_of_face(face_id)[i]];
>
> Ceci marche bien pour des elements d'ordre 1, mais quand j'utiliser des
> IM_GAUSS_PARALLELEPIPED(2,2) ou IM_HEXAHEDRON(2), ca ne marche plus.
> D'une part, nb_points_of_face renvoit 2 dans le cas 2d et 4 dans le cas
> 3D au lieu de 3 et 9 respectivement, d'autre part, la numérotation
> renvoyer par ind_dof_of_element semble bizarre...
> Est-ce que j'ai raté quelque chose ou est ce normal? Y a t il un moyen
> de récuperer tout les dof de la face?
> (je travail sur la version stable de getfem)
>
> cordialement,
> Cédric Venet
>
> ps: vaut il mieux poster en francais ou en anglais sur cette ml?
>
> _______________________________________________
> Getfem-users mailing list
> address@hidden
> https://mail.gna.org/listinfo/getfem-users