Hi, you would not expect a dg element to have indices at vertices.the indices or dofs are located inside the cell.hope this helps. Bärbel -- Diese Nachricht wurde von meinem Android-Mobiltelefon mit K-9 Mail gesendet.
"M. Sebastian Pauletti" <[email protected]> schrieb: 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
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
