[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Getfem-users] Enumerer tout les DOF d'une face à l'ordre 2
From: |
Cédric Venet |
Subject: |
[Getfem-users] Enumerer tout les DOF d'une face à l'ordre 2 |
Date: |
Sun, 11 Oct 2009 22:54:56 +0200 |
User-agent: |
Mozilla/5.0 (Windows; U; Windows NT 6.1; fr; rv:1.9.1.4pre) Gecko/20090915 Lightning/1.0pre Thunderbird/3.0b4 |
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] Enumerer tout les DOF d'une face à l'ordre 2,
Cédric Venet <=