WIP: PM smoothing; BUG: convergence of binary search