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 <
ngockien.lea...@gmail.com> đã 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 <
> bruno.turck...@gmail.com> đã viết:
>
>> Pham Ngoc Kien,
>>
>> Le lun. 18 févr. 2019 à 19:33, Phạm Ngọc Kiên <ngockien.lea...@gmail.com>
>> 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 dealii+unsubscr...@googlegroups.com.
>> 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 dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
