added interp_polar function

Function for angular interpolation accepting as input any range of angles.
If I am reading the rest of your code correctly, you only need the two smallest distances from the sorting: using `partition` or `argpartition` should result in a significant speed-up.