Reverted to Todd's implementation of position_estimation().

After comparing the two, his turned out to be way faster. I suspect
NumPy is doing something parallel and cool under the hood ...