Forwarded in name of Lauren who is having issues to get her emails
accepted by the list.
---------- Forwarded message ----------
Hi everyone,
I want to access the global dof index of a solution vector associated with
a Discontinuous Galerkin element, but when I call vertex_dof_index, it
gives me an error which seems to indicate that the number of dofs per
vertex is zero.
In particular, I have a solution vector, rstress, associated with
corresponding dof handler and fe, where the finite element was defined:
FESystem<dim> rstress_fe (FE_DGQ<dim>(1), dim*dim)
Then later on, inside a loop over the cells, I want to call:
rstress(cell->vertex_dof_index(quad_point, degree))
But I get the following error message:
An error occurred in line <1209> of file
</w/ferguson/deal.II/include/deal.II/dofs/dof_accessor.templates.h> in
function
static unsigned int
dealii::internal::DoFAccessor::Implementation::get_vertex_dof_index(const
dealii::DoFHandler<dim, spacedim>&, unsigned int, unsigned int, unsigned
int) [with int dim = 2, int spacedim = 2]
The violated condition was:
local_index < dof_handler.selected_fe->dofs_per_vertex
The name and call sequence of the exception was:
ExcIndexRange(local_index, 0,
dof_handler.selected_fe->dofs_per_vertex)
Additional Information:
Index 0 is not in [0,0[
The problem seems to be that using DGQ(1)
dof_handler.selected_fe->dofs_per_vertex returns 0,
which is not what one would expect.
Any ideas why or how to make it work as expected? Thanks!
Lauren Ferguson
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
