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 <[email protected]> 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 > [email protected] > https://mail.gna.org/listinfo/getfem-users > _______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
