Dear Julian,
In general its not easy to provide the functionality that you're looking
for because not every degree-of-freedom is associated with a support point
(e.g. the FE_DGPMonomial element) and let alone a vertex (e.g. an FE_Q of
polynomial order 2 has DoFs with support points at face and
ok. Something like this:
cell = dof_handler.begin_active ();
endc = dof_handler.end ();
for (; cell != endc; ++cell)
{
for (unsigned int vertex = 0;
vertex < GeometryInfo::vertices_per_cell; ++vertex)
{
cell->vertex(vertex)
}
}
Thank you!
...
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
wrote:
> Dear all,
>
> if I hav
Dear all,
if I have
FE_Q fe; // to be Q1 finite element in 2D
DoFHandler 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
