Bugfix: the method geometryInFather() should not return a Geometry<dim,dimworld>.

Instead, it now returns LocalGeometry (Geometry<dim,dim>), because the return
value of geometryInFather is the element expressed in local coordinates in the
father, and the father is dim-dimensional.  Of course this only makes a
difference when dim!=dimworld.

[[Imported from SVN: r4185]]

grid/common/entity.hh

@@ -179,6 +179,11 @@ namespace Dune
     /** \brief The geometry type of this entity */
     typedef typename GridImp::template Codim<0>::Geometry Geometry;
 
+    /** \brief The geometry type of this entity when the geometry is expressed
+       embedded in the father element.  This differs from Geometry only when
+       dim != dimworld.*/
+    typedef typename GridImp::template Codim<0>::LocalGeometry LocalGeometry;
+
     /** \brief EntityPointer types of the different codimensions */
     template <int cd>
     struct Codim
@@ -275,7 +280,7 @@ namespace Dune
        implementation of numerical algorithms is only done for simple discretizations.
        Assumes that meshes are nested.
      */
-    const Geometry& geometryInFather () const
+    const LocalGeometry& geometryInFather () const
     {
       return realEntity.geometryInFather();
     }