ok. Something like this:
cell = dof_handler.begin_active ();
endc = dof_handler.end ();
for (; cell != endc; ++cell)
{
for (unsigned int vertex = 0;
vertex < GeometryInfo<dim>::vertices_per_cell; ++vertex)
{
cell->vertex(vertex)
}
}
Thank you!
... the only thing that if one has a loop through Q1-dofs, then to get the
corresponding vertex coordinates of the i-th dof in such a way might be
expensive. The solution is to put this part into the preprocessing-phase
though.
On Thursday, November 10, 2016 at 2:45:34 PM UTC+1, Deepak Gupta wrote:
>
> Dear Julian,
>
> This question has been asked in the past and there are already some
> threads which might be of help to you:
>
> https://groups.google.com/forum/#!searchin/dealii/dof$20coordinates%7Csort:relevance
>
> Best,
> Deepak
>
> On Thu, Nov 10, 2016 at 2:41 PM, Julian Dorn <julian...@gmail.com
> <javascript:>> wrote:
>
>> Dear all,
>>
>> if I have
>> FE_Q<dim> fe; // to be Q1 finite element in 2D
>> DoFHandler<dim> dof_handler;
>>
>> how to get coordinates of dofs (for Q1 this will be exactly coordinates
>> of quads' vertices)?
>>
>> Thank you in advance!
>>
>> --
>> 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+un...@googlegroups.com <javascript:>.
>> 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.
