possible bug in dirichletconstraints.hh
The values
vector that contains the Dirichlet data has the dimension equal to numdofs
of the polygon that the entity under consideration belongs to, while logically the entity should only touch the values corresponding to its own vertices. As the vector values
is being set to zero each time an entity is visited, wLocal
will contain wrong values at the end of the assembly process.