Fix a numerical issue in the function 'angular_separation'

For identical input coordinates, the function might have
sometimes returned _NaN. An error is now thrown if this
happens. Moreover, for identical input coordinates, it just
returns zero.