I am sorry as my question is not clear enough.
What I mean is when we loop over all degree of freedom in a cell, for
example,
for (unsigned int i = 0; i < dofs_per_cell; ++i) { .......... }
Do we have any method to get the position of the i_th degree of
freedom inside this loop?
I would like to thank you very much!
Best,
Pham Ngoc Kien
Vào Th 3, 19 thg 2, 2019 vào lúc 14:47 Phạm Ngọc Kiên <
[email protected]> đã viết:
> Hi,
> I have tried with the methods described in the FAQ and figured out that:
> 1. As I used edge-based element (FE_NedelecSZ), I queried if my element
> had support points by *fe.has_support_points()*. As a consequence, it
> returned zero.
> The *FiniteElement::get_unit_support_points() *also returned an empty
> vector.
> As there was no support point, DoFTools::map_dofs_to_support_points()
> could not work as you mentioned above.
>
> 2. I have tried some thing with the third method. And here-below is my
> code:
>
> Point<dim> p{0.5,0.5,0.5}; //position in reference cell
> Quadrature<dim> q(p);
> FEValues<dim> fe_values_q(fe, q, update_quadrature_points);
> fe_values_q.reinit(cell);
>
> //position in the real cell
> std::vector<Point<dim>> dof_pos = fe_values_q.get_quadrature_points();
> for (int k = 0; k <dof_pos.size() ; ++k) {
> for (int j = 0; j <dim ; ++j) {
> std::cout<<dof_pos[k][j];
> std::cout<<"\n";
> }
> }
>
> From this codes, I can get the real position of arbitrary points in my
> current reference cell.
>
> However, my purpose is to get the real position of a given degree of freedom.
>
> Could you please tell me how to get the position of a degree of freedom in
> the reference cell?
>
> Thank you very much for your help.
>
> Best,
>
> Pham Ngoc Kien
>
>
> Vào Th 3, 19 thg 2, 2019 vào lúc 11:47 Phạm Ngọc Kiên <
> [email protected]> đã viết:
>
>> Yes, thank you very much.
>> I will dig deeper into this.
>> Best,
>> Pham Ngoc Kien
>>
>> Vào Th 3, 19 thg 2, 2019 vào lúc 11:23 Bruno Turcksin <
>> [email protected]> đã viết:
>>
>>> Pham Ngoc Kien,
>>>
>>> Le lun. 18 févr. 2019 à 19:33, Phạm Ngọc Kiên <[email protected]>
>>> a écrit :
>>>
>>>> With the DoFTools::map_dofs_to_support_points() , I saw in the library
>>>> manual that it can be used with no edge elements or the like.
>>>> As my finite element is edge-based one, I wonder if this function work?
>>>>
>>> No the function won't work but there are two other methods described in
>>> the FAQ that you could use. In particular, you can create a special
>>> quadrature at the points of interest.
>>>
>>> Best,
>>>
>>> Bruno
>>>
>>> --
>>> The deal.II project is located at http://www.dealii.org/
>>> For mailing list/forum options, see
>>> https://groups.google.com/d/forum/dealii?hl=en
>>> ---
>>> You received this message because you are subscribed to the Google
>>> Groups "deal.II User Group" group.
>>> To unsubscribe from this group and stop receiving emails from it, send
>>> an email to [email protected].
>>> For more options, visit https://groups.google.com/d/optout.
>>>
>>
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
