EllipsoidFactor: use wiggle-contract-dilate-update optimisation cycle

In concert with strict contact conditions, this appears to result in
fewer ellipsoids with surface penetration and closer fitting to the
surface