What is uint* facet_dofs = new uint[10*dof_map->num_facet_dofs()]; doing? Why is it multiplied by 10? /Anders