Hi
In a code for 3D linear elasticity, I'm using an FESystem made up of 3 FE_DGQ
elements. I'm trying to access, in a given cell, the degrees of freedom
associated with a given vertex. Is there a way of doing this, for discontinuous
elements?
The call "cell->vertex_dof_index(vertex_no, component)" ( with, for example,
vertex_no = 1, component = 0) is giving me an error, with the following
description:
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[
(It seems that the dofs on the cell are not directly associated with particular
vertices: if I call fe.dofs_per_vertex, it's returning the value 0, although
fe.dofs_per_cell is returning the value 24.)
Thanks,
Beverley
###
UNIVERSITY OF CAPE TOWN
This e-mail is subject to the UCT ICT policies and e-mail disclaimer published
on our website at http://www.uct.ac.za/about/policies/emaildisclaimer/ or
obtainable from +27 21 650 4500. This e-mail is intended only for the person(s)
to whom it is addressed. If the e-mail has reached you in error, please notify
the author. If you are not the intended recipient of the e-mail you may not
use, disclose, copy, redirect or print the content. If this e-mail is not
related to the business of UCT it is sent by the sender in the sender's
individual capacity.
###
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
