Speed up the error estimation by a factor 4!

The thing was to call outside() on the intersection iterator only once and
store the EntityPointer. Now it is only slightly more expensive than assembly on tetrahedral
mesh and the same cost as assembly on hex mesh.

[[Imported from SVN: r3296]]