docs: def get_interp_val(m, theta, phi, nest=False, lonlat=False):
      """
      Drop-in replacement for healpy `~healpy.pixelfunc.get_interp_val`.
      """
      lon, lat = _healpy_to_lonlat(theta, phi, lonlat=lonlat)
      return interpolate_bilinear_lonlat(lon, lat, m, order='nested' if nest else 'ring')