reciprocal_divide: update/correction of the algorithm