In gauss integration factory use new ElementGeometry

Allows us to determine the geometry of an element robustly, might fix
bugs in nodal integration.
In variable_order_integrator_factory use new ElementGeometry