def distance(self, p1, p2): '''Return the L{equirectangular_} distance in C{radians squared}. ''' d, _ = unrollPI(p1.b, p2.b, wrap=self._wrap) if self._adjust: d *= _scaler(p1.a, p2.a) return d**2 + (p2.a - p1.a)**2 # like equirectangular_ d2
def distance(self, p1, p2): '''Return the L{equirectangular_} distance in C{radians squared}. ''' d, _ = unrollPI(p1.lam, p2.lam, wrap=self._wrap) if self._adjust: d *= _scaler(p1.phi, p2.phi) return hypot2(d, p2.phi - p1.phi) # like equirectangular_ d2
def _distances(self, x, y): # (x, y) radians for xk, yk in zip(self._xs, self._ys): d, _ = unrollPI(xk, x, wrap=self._wrap) if self._adjust: d *= _scaler(yk, y) yield d**2 + (yk - y)**2 # like equirectangular_ distance2