Re: [Getfem-users] Get Nedelec DOF orientations

[Yves Renard]

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
> address@hidden
> https://mail.gna.org/listinfo/getfem-users


-- 

  Yves Renard (address@hidden)       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

---------