EllipsoidFactor: rotate ellipsoid aligned to contact points

Build covariance matrix from contact points' coordinates and use an
eigenvalue decomposition as the new ellipsoid orientation. By default,
seems that the short axis is aligned with contact points.