Hi Christian, > Now I need a mapping between the index of a vertex in the vector > vertices and the dof index.
A degree of freedom does not need to correspond to a vertex of the cell, it doesn't even have to correspond to a point value at all (depending on the finite element space). Maybe you can elaborate on what you want to achieve? In case you are doing an integration loop, look at the quadrature points instead. They always have a position associated with them. If you really want to stick to what you are doing (and you are sure you never want to support higher order elements or non-lagrange element types), take a look at the following functions: DoFAccessor::vertex_dof_index, TriaAccessor::vertex_index Best, Timo -- Timo Heister http://num.math.uni-goettingen.de/~heister _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
