Comment 4 for bug 488661

If we return the identity matrix with size equal to the size of the incoming quadrature points, the generated code for tabulate_tensor() looks correct. Even if we add a cell integral such that:

L = v*g*dx + v*g*ds

it looks correct.

The problem with this approach is that since the dofs of a quadrature element is associated with the quadrature points it is problematic to define the element and dof_map information consistently.
What should for instance the output be of element.space_dimension(), element.evaluate_dof(),
dof_map.tabulate_coordinates(), dof_map.tabulate_dofs()?