Re: [Getfem-users] Get Nedelec DOF orientations

[Torquil Macdonald Sørensen]

Hi Yves!

Thanks, I got the same answer from Ronan Perrussel (off the email list).
After implementing it, it solved my problem, so I managed to construct
the vectors in the kernel of curl that correspond to the divergence of
scalar functions, in both 2d and 3d.

Now it remains for me to find the vector fields in the kernel of curl
that are not the divergence of scalar functions.

Best regards and thanks,
Torquil Sørensen

On 11/12/15 15:21, Yves Renard wrote:
>
> 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
>