Dear Torquil,
The determination of the orientation is determined by a particular choice made by the method P1_nedelec_::mat_trans(...) in getfem_fem.cc. The choice is made by a comparison with an a priori direction "norient" : the dof represent always a vector having a positive scalar product with norient ... Best regards, Yves. Le 09/12/2015 12:09, Torquil Macdonald Sørensen a écrit : > Hi! > > Given a tetrahedron "cv" with Nedelec DOFs, and a local DOF number d, is > there a simple way in GetFEM to get a vector that points in the same > direction as the vectorial DOF d? > > I'm trying to solve an electromagnetic problem, and I need to construct > the kernel of the curl operator. For each node, I must have a list of > all Nedelec DOFs on edges connected to the node, and I must know if they > points towards or away from the node. Therefore I need the direction of > the DOF. > > I have tried to obtain the value of the DOF at its "home location", by > using a fem_interpolation_context, but I'm not completely sure I'm doing > it correctly. By this method I am able to contruct vectors in the kernel > of curl in 2d, but in 3d I get vectors not in the kernel. > > Best regards and thanks, > Torquil Sørensen > > > _______________________________________________ > Getfem-users mailing list > [email protected] > https://mail.gna.org/listinfo/getfem-users -- Yves Renard ([email protected]) tel : (33) 04.72.43.87.08 Pole de Mathematiques, INSA-Lyon fax : (33) 04.72.43.85.29 20, rue Albert Einstein 69621 Villeurbanne Cedex, FRANCE http://math.univ-lyon1.fr/~renard --------- _______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
