pygeodesy 25.9.9__py2.py3-none-any.whl → 25.11.5__py2.py3-none-any.whl

pygeodesy/units.py

@@ -820,6 +820,11 @@ def _isMeter(obj, iscalar=True):
     return isinstance(obj, _Meters) or (iscalar and _isScalar(obj))
 
 
+def _isRadians(obj, iscalar=True):
+    # Check for valid radian types.
+    return isinstance(obj, _Radians) or (iscalar and _isScalar(obj))
+
+
 def _isRadius(obj, iscalar=True):
     # Check for valid earth radius types.
     return isinstance(obj, _Radii) or (iscalar and _isScalar(obj))

pygeodesy/utily.py

@@ -28,7 +28,7 @@ from math import acos, asin, asinh, atan2 as _atan2, cos, degrees, fabs, \
                  radians, sin, sinh, tan as _tan  # pow
 
 __all__ = _ALL_LAZY.utily
-__version__ = '25.09.09'
+__version__ = '25.10.29'
 
 # sqrt(3) <https://WikiPedia.org/wiki/Square_root_of_3>
 _COS_30,  _SIN_30 = 0.86602540378443864676, _0_5  # sqrt(3) / 2
@@ -807,15 +807,16 @@ def _sin0cos2(q, r, sign):
     return s, c
 
 
-def SinCos2(x):
+def SinCos2(x, unit=Radians):
     '''Get C{sin} and C{cos} of I{typed} angle.
 
        @arg x: Angle (L{Degrees}, L{Radians} or scalar C{radians}).
+       @kwarg unit: The C{B{x}} unit (L{Degrees}, L{Radians}).
 
        @return: 2-Tuple (C{sin(B{x})}, C{cos(B{x})}).
     '''
-    return sincos2d(x) if isinstanceof(x, Degrees, Degrees_) else (
-#          sincos2(x)  if isinstanceof(x, Radians, Radians_) else
+    return sincos2d(x) if unit is Degrees or isinstanceof(x, Degrees, Degrees_) else (
+#          sincos2(x)  if unit is Radians or isinstanceof(x, Radians, Radians_) else
            sincos2(Radians(x)))  # assume C{radians}
 
 

pygeodesy/vector2d.py

@@ -641,7 +641,7 @@ def triaxum5(points, useZ=True):
     A = []
     for i, p in enumerate(ps):
         v = _otherV3d(useZ=useZ, i=i, points=p)
-        A.append(v.x2y2z2)
+        A.append(v.x2y2z23)
 
     with _numpy(triaxum5, n=n) as _np:
         A = _np.array(A)

pygeodesy/vector3dBase.py

@@ -10,9 +10,9 @@ A pure Python implementation of vector-based functions by I{(C) Chris Veness
 
 from pygeodesy.basics import _copysign, islistuple, isscalar, \
                               map1, map2, _signOf, _zip
-from pygeodesy.constants import EPS, EPS0, INT0, PI, PI2, \
-                               _1_0, isnear0, isnear1, isneg0, \
-                               _copysignINF, _float0, _pos_self
+from pygeodesy.constants import EPS, EPS0, INT0, PI, PI2, _1_0, \
+                               _copysignINF, isnear0, isnear1, \
+                               _flipsign, _float0, _pos_self
 from pygeodesy.errors import CrossError, VectorError, _xcallable, _xError
 from pygeodesy.fmath import euclid_, fdot, fdot_, hypot_, hypot2_  # _MODS.fmath.fma
 from pygeodesy.interns import _coincident_, _colinear_, _COMMASPACE_, _xyz_
@@ -29,7 +29,7 @@ from pygeodesy.utily import atan2, sincos2,  fabs
 from math import ceil as _ceil, floor as _floor, trunc as _trunc
 
 __all__ = _ALL_LAZY.vector3dBase
-__version__ = '25.08.31'
+__version__ = '25.10.25'
 
 
 class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
@@ -112,9 +112,20 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
 
     cmp = __cmp__
 
-    def __divmod__(self, other):  # PYCHOK no cover
-        '''Not implemented.'''
-        return _NotImplemented(self, other)
+    def __divmod__(self, scalar):  # PYCHOK no cover
+        '''Apply C{scalar} divmod to this vector's components.
+
+           @arg scalar: Divisor (C{scalar}).
+
+           @return: 2-Tuple C{(d, m)} each a L{Vector3d}.
+        '''
+        s = Scalar(divisor=scalar)
+        d, m = [], []
+        for x in self.xyz3:
+            q, r = divmod(x, s)
+            d.append(q)
+            m.append(r)
+        return self.classof(d), self.classof(m)
 
     def __eq__(self, other):
         '''Is this vector equal to an other vector?
@@ -132,15 +143,22 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         return _NotImplemented(self)  # must return C{float}
 
     def __floor__(self):  # PYCHOK no cover
-        '''Return a vector with the C{floor} of these components.
+        '''Return a vector with the C{floor} of each components.
 
            @return: Floor-ed (L{Vector3d}).
         '''
         return self._mapped(_floor)
 
-    def __floordiv__(self, other):  # PYCHOK no cover
-        '''Not implemented.'''
-        return _NotImplemented(self, other)
+    def __floordiv__(self, scalar):  # PYCHOK no cover
+        '''Floor-divide this vector by a scalar, C{this // B{scalar}}.
+
+           @arg scalar: The divisor (C{scalar}).
+
+           @return: Floored quotient (L{Vector3d}).
+
+           @raise TypeError: Non-scalar B{C{scalar}}.
+        '''
+        return self.floorDividedBy(scalar)
 
     def __ge__(self, other):
         '''Is this vector longer than or equal to an other vector?
@@ -183,9 +201,14 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self._xyz(self.plus(other))
 
-    def __ifloordiv__(self, other):  # PYCHOK no cover
-        '''Not implemented.'''
-        return _NotImplemented(self, other)
+    def __ifloordiv__(self, scalar):  # PYCHOK no cover
+        '''Floor-divide this vector by a scalar I{in-place}, C{this //= B{scalar}}.
+
+           @arg scalar: The divisor (C{scalar}).
+
+           @raise TypeError: Non-scalar B{C{scalar}}.
+        '''
+        return self._xyz(self.floorDividedBy(scalar))
 
     def __imatmul__(self, other):  # PYCHOK Python 3.5+
         '''Cross multiply this and an other vector I{in-place}, C{this @= B{other}}.
@@ -215,9 +238,16 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''Not implemented, see method C{ints}.'''
         return _NotImplemented(self)  # must return C{int}
 
-    def __ipow__(self, other, *mod):  # PYCHOK no cover
-        '''Not implemented.'''
-        return _NotImplemented(self, other, *mod)
+    def __ipow__(self, scalar, *mod):  # PYCHOK no cover
+        '''Raise each component I{in-place} as C{pow(C, B{scalar})}.
+
+           @arg scalar: Exponent (C{scalar}).
+
+           @return: Power (L{Vector3d}).
+
+           @raise TypeError: Non-scalar B{C{scalar}}.
+        '''
+        return self._xyz(self.pow(scalar, *mod))
 
     def __isub__(self, other):
         '''Subtract an other vector from this one I{in-place}, C{this -= B{other}}.
@@ -280,9 +310,19 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self.cross(other)
 
-    def __mod__(self, other):  # PYCHOK no cover
-        '''Not implemented.'''
-        return _NotImplemented(self, other)
+    def __mod__(self, modulus):  # PYCHOK no cover
+        '''Apply C{scalar} modulo to this vector's components.
+
+           @arg modulus: Modulus (C{scalar}).
+
+           @return: Modulo (L{Vector3d}).
+        '''
+        m = Scalar(modulus=modulus)
+
+        def _mod(x):
+            return x % m
+
+        return self._mapped(_mod)
 
     def __mul__(self, scalar):
         '''Multiply this vector by a scalar, C{this * B{scalar}}.
@@ -318,13 +358,29 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self if _pos_self else self.copy()
 
-    def __pow__(self, other, *mod):  # PYCHOK no cover
-        '''Not implemented.'''
-        return _NotImplemented(self, other, *mod)
+    def __pow__(self, scalar, *mod):  # PYCHOK no cover
+        '''Return a vector with components as C{pow(C, B{scalar})}.
+
+           @arg scalar: Exponent (C{scalar}).
+
+           @return: Power (L{Vector3d}).
+
+           @raise TypeError: Non-scalar B{C{scalar}}.
+        '''
+        return self.pow(scalar, *mod)
 
-    __radd__ = __add__  # PYCHOK no cover
+    def __radd__(self, other):  # PYCHOK no cover
+        '''Add this vector to an other vector, C{B{other} + this}.
+
+           @arg other: The other vector (L{Vector3d}).
+
+           @return: Sum (L{Vector3d}).
+
+           @raise TypeError: Incompatible B{C{other}} C{type}.
+        '''
+        return self.others(other).plus(self)
 
-    def __rdivmod__ (self, other):  # PYCHOK no cover
+    def __rdivmod__(self, other):
         '''Not implemented.'''
         return _NotImplemented(self, other)
 
@@ -333,7 +389,7 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
 #       '''
 #       return self.toRepr()
 
-    def __rfloordiv__(self, other):  # PYCHOK no cover
+    def __rfloordiv__(self, other):
         '''Not implemented.'''
         return _NotImplemented(self, other)
 
@@ -348,11 +404,11 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self.others(other).cross(self)
 
-    def __rmod__(self, other):  # PYCHOK no cover
+    def __rmod__(self, other):
         '''Not implemented.'''
         return _NotImplemented(self, other)
 
-    __rmul__ = __mul__
+    __rmul__ = __mul__  # scalar other
 
     def __round__(self, *ndigits):  # PYCHOK no cover
         '''Return a vector with these components C{rounded}.
@@ -361,8 +417,11 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
 
            @return: Rounded (L{Vector3d}).
         '''
+        def _rnd(x):
+            return round(x, *ndigits)
+
         # <https://docs.Python.org/3.12/reference/datamodel.html?#object.__round__>
-        return self.classof(*(round(_, *ndigits) for _ in self.xyz3))
+        return self._mapped(_rnd)
 
     def __rpow__(self, other, *mod):  # PYCHOK no cover
         '''Not implemented.'''
@@ -379,9 +438,9 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self.others(other).minus(self)
 
-    def __rtruediv__(self, scalar):  # PYCHOK no cover
+    def __rtruediv__(self, other):
         '''Not implemented.'''
-        return _NotImplemented(self, scalar)
+        return _NotImplemented(self, other)
 
 #   def __str__(self):
 #       '''Return the default C{str(self)}.
@@ -527,22 +586,39 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         self._crosserrors = bool(raiser)
 
-    def dividedBy(self, divisor):
+    def dividedBy(self, scalar):
         '''Divide this vector by a scalar.
 
-           @arg divisor: The divisor (C{scalar}).
+           @arg scalar: The divisor (C{scalar}).
 
            @return: New, scaled vector (L{Vector3d}).
 
-           @raise TypeError: Non-scalar B{C{divisor}}.
+           @raise TypeError: Non-scalar B{C{scalar}}.
 
-           @raise VectorError: Invalid or zero B{C{divisor}}.
+           @raise VectorError: Invalid or zero B{C{scalar}}.
         '''
-        d = Scalar(divisor=divisor)
+        d = Scalar(divisor=scalar)
         try:
             return self._times(_1_0 / d)
         except (ValueError, ZeroDivisionError) as x:
-            raise VectorError(divisor=divisor, cause=x)
+            raise VectorError(divisor=scalar, cause=x)
+
+    def dividedBy_(self, scalar_x, *y_z):
+        '''Divide this vector by separate X, Y and Z divisors.
+
+           @arg scalar_x: X divisor (C{scalar}) or a vector's
+                          X, Y, and Z components as divisors
+                          (C{Cartesian}, L{Ecef9Tuple}, C{Nvector},
+                          L{Vector3d}, L{Vector3Tuple}, L{Vector4Tuple}).
+           @arg y_z: Y and Z divisors (C{scalar}, C{scalar}), ignored
+                     if B{C{scalar_x}} is not C{scalar}.
+
+           @return: New, scaled vector (L{Vector3d}).
+
+           @raise VectorError: Invalid B{C{scalar_x}} or B{C{y_z}}.
+        '''
+        x, y, z = _xyz3(self.dividedBy_, scalar_x, *y_z)
+        return self.classof(self.x / x, self.y / y, self.z / z)
 
     def dot(self, other):
         '''Compute the dot (scalar) product of this and an other vector.
@@ -595,6 +671,43 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self._mapped(_float0)
 
+    def floorDividedBy(self, scalar):
+        '''Floor-divide this vector by a scalar.
+
+           @arg scalar: The divisor (C{scalar}).
+
+           @return: New, scaled vector (L{Vector3d}).
+
+           @raise TypeError: Non-scalar B{C{scalar}}.
+
+           @raise VectorError: Invalid or zero B{C{scalar}}.
+        '''
+        d = Scalar(divisor=scalar)
+        try:
+            def _floor_d(x):
+                return x // d
+
+            return self._mapped(_floor_d)
+        except (ValueError, ZeroDivisionError) as x:
+            raise VectorError(divisor=scalar, cause=x)
+
+    def floorDividedBy_(self, scalar_x, *y_z):
+        '''Floor-divide this vector by separate X, Y and Z divisors.
+
+           @arg scalar_x: X divisor (C{scalar}) or a vector's
+                          X, Y, and Z components as divisors
+                          (C{Cartesian}, L{Ecef9Tuple}, C{Nvector},
+                          L{Vector3d}, L{Vector3Tuple}, L{Vector4Tuple}).
+           @arg y_z: Y and Z divisors (C{scalar}, C{scalar}), ignored
+                     if B{C{scalar_x}} is not C{scalar}.
+
+           @return: New, scaled vector (L{Vector3d}).
+
+           @raise VectorError: Invalid B{C{scalar_x}} or B{C{y_z}}.
+        '''
+        x, y, z = _xyz3(self.floorDividedBy_, scalar_x, *y_z)
+        return self.classof(self.x // x, self.y // y, self.z // z)
+
     @Property
     def _fromll(self):
         '''(INTERNAL) Get the latlon reference (C{LatLon}) or C{None}.
@@ -617,11 +730,8 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
             y =  y / z
 #           z = _1_0
         else:
-            if isneg0(z):
-                x = -x
-                y = -y
-            x = _copysignINF(x)
-            y = _copysignINF(y)
+            x = _copysignINF(_flipsign(x, z))
+            y = _copysignINF(_flipsign(y, z))
 #           z =  NAN
         return self.classof(x, y, _1_0)
 
@@ -705,7 +815,7 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
 
            @see: Properties L{length2} and L{euclid}.
         '''
-        return Float(length=hypot_(self.x, self.y, self.z))
+        return Float(length=hypot_(*self.xyz3))
 
     @Property_RO
     def length2(self):  # __dict__ value overwritten by Property_RO C{_united}
@@ -713,7 +823,7 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
 
            @see: Property L{length} and method C{equirectangular}.
         '''
-        return Float(length2=hypot2_(self.x, self.y, self.z))
+        return Float(length2=hypot2_(*self.xyz3))
 
     def _mapped(self, func):
         '''(INTERNAL) Apply C{func} to all components.
@@ -815,12 +925,28 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
            @arg y_z: Y and Z components (C{scalar}, C{scalar}),
                      ignored if B{C{other_x}} is not C{scalar}.
 
-           @return: New, vectiorial vector (L{Vector3d}).
+           @return: New, vectorial vector (L{Vector3d}).
 
            @raise ValueError: Invalid B{C{other_x}} or B{C{y_z}}.
         '''
         return self._plus(*_xyz3(self.plus_, other_x, *y_z))
 
+    def pow(self, scalar, *mod):
+        '''Raise each X, Y and Z to C{scalar} power.
+
+           @arg scalar: Exponent (C{scalar}).
+
+           @return: Power (L{Vector3d}).
+
+           @raise TypeError: Non-scalar B{C{scalar}}.
+        '''
+        p = Scalar(power=scalar)
+
+        def _pow(x):
+            return pow(x, p, *mod)
+
+        return self._mapped(_pow)
+
     def rotate(self, axis, theta, fma=False):
         '''Rotate this vector around an axis by a specified angle.
 
@@ -863,6 +989,16 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''DEPRECATED, use method C{rotate}.'''
         return self.rotate(axis, theta)  # PYCHOK no cover
 
+    def _roty(self, *pos, **name):
+        '''(INTERNAL) Prolate rotation, for C{+1 if pos else -1}.
+        '''
+        v = self.copy(**name)
+        if pos:
+            x, _, z = v.xyz3
+            _update_all(v, needed=3)
+            v._x, v._z = (-z, x) if pos[0] else (z, -x)
+        return v
+
     def times(self, factor):
         '''Multiply this vector by a scalar.
 
@@ -989,8 +1125,10 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
     def _xyz(self, x_xyz, *y_z):
         '''(INTERNAL) Set the C{_x}, C{_y} and C{_z} attributes.
         '''
-        _update_all(self, needed=3)
-        self._x, self._y, self._z = _xyz3(_xyz_, x_xyz, *y_z)
+        xyz = _xyz3(_xyz_, x_xyz, *y_z)
+        if self.xyz3 != xyz:
+            _update_all(self, needed=3)
+            self._x, self._y, self._z = xyz
         return self
 
     @property_RO
@@ -999,9 +1137,15 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self.x, self.y, self.z
 
-    @property_RO
+    @Property_RO
     def x2y2z2(self):
-        '''Get the X, Y and Z components I{squared} (3-tuple C{(x**2, y**2, z**2)}).
+        '''Get the X, Y and Z components, I{squared} (L{Vector3Tuple}).
+        '''
+        return _MODS.namedTuples.Vector3Tuple(*self.x2y2z23, name=self.name)
+
+    @property_RO
+    def x2y2z23(self):
+        '''Get the X, Y and Z components, I{squared} (3-tuple C{(x**2, y**2, z**2)}).
         '''
         return self.x**2, self.y**2, self.z**2
 
@@ -1037,7 +1181,7 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
 
 
 def _xyz3(where, x_xyz, *y_z):  # in .cartesianBase._rtp3
-    '''(INTERNAL) Get , Y and Z as 3-tuple C{(x, y, z)}.
+    '''(INTERNAL) Get X, Y and Z as 3-tuple C{(x, y, z)}.
     '''
     try:
         xyz3 = map1(_float0, x_xyz, *y_z) if y_z else (  # islistuple for Vector*Tuple

{pygeodesy-25.9.9.dist-info → pygeodesy-25.11.5.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: pygeodesy
-Version: 25.9.9
+Version: 25.11.5
 Summary: Pure Python geodesy tools
 Home-page: https://GitHub.com/mrJean1/PyGeodesy
 Author: Jean M. Brouwers
@@ -8,7 +8,7 @@ Author-email: mrJean1@Gmail.com
 Maintainer: Jean M. Brouwers
 Maintainer-email: mrJean1@Gmail.com
 License: MIT
-Keywords: AER Albers altitude Andoyer annulus antipode area attitude Authalic auxiliary azimuth azimuthal azimuth-elevation-range bearing bank Barsky Barth beta bi-quadratic boolean cached Cagnoli cartesian Cassini Cassini-Soldner chord circle-intersections circumcenter circumcircle circumradius clip Cohen Cohen-Sutherland Collins composite conformal conic constants contact-triangle Cook Correia cosines-law coverage curvature cylindrical datum deprecation deficit development discrete distance Douglas earth east-north-up eccentricity ECEF elevation ellipsoid ellipsoidal-latitude-beta ellipsoidal-longitude-omega elliptic ENU EPSG equal-area equidistant equirectangular ETM ETRF Euclidean even-odd-rule ExactTM excess Farrell Farrell-Barth Ferrari-solution Field-Of-View flattening fma fmath footpoint footprint Forster Forster-Hormann-Popa Forsythe FOV fractional Frechet Fréchet frustum Fsum fused-multiply-add GARS geocentric GeoConvert GeodesicExact geodesy geodetic GeodSolve GeodTest geographiclib Geohash geoid geoidHeight GeoidHeights georef Girard gnomonic gons grades gradians Greiner Greiner-Hormann Hartzell Hausdorff Haversine heading hectare height Heikkinen Heron Hodgman horizon Hormann Hubeny IDW incenter incirle infix_@_operator inradius intermediate interpolate intersect intersection intersection3d intersections IntersectTool Inverse-Distance-Weighting Isometric ITRF Jacobi Jacobi-Conformal Jarque-Bera Jekel Karney Krueger Krüger kurtosis Lambert latitude law-of-cosines least-squares Lesh L_Huilier LHuilier Liang Liang-Barsky linearize Line-Of-Sight LocalCartesian local-tangent-plane local-x-y-z longitude LOS loxodrome lstsq LTP lune LV03 LV95 mean memoize memoized Mercator Meeus MGRS nearest NED Niemeyer non-finite normalize Norrdine north-east-down numpy n-vector Nvector oblate omega orthographic orthometric-height OSGB OSGR overlap parallel parallel-of-latitude Parametric path-intersection perimeter Peucker Pierlot pitch plumb Point-Of-View polar Popa POV precision-cubic-root precision-hypotenuse precision-powers precision-running-summation precision-square-root precision-summation prolate Pseudo-Mercator pygeodesy PyInstaller PyPy quartic radical radii radius Ramer Ramer-Douglas-Peucker Rectifying Reduced resect resection Rey-Jer Reumann Reumann-Witkam rho-theta-phi rhumb RhumbSolve running-linear-regression running-statistics running-stats running-summation scipy secant semi-perimeter sexagecimal simplify skewness Snellius Snellius-Pothenot Snyder Soddy Soddy-circles Soldner sphere sphere-intersections spherical-deficit spherical-excess spherical-polar spherical-triangle squared-quartic standard-deviation stereographic Sudano surface-area Sutherland Sutherland-Hodgman tangent-circles Terrestrial-Reference-Frame Thomas Tienstra tilt TMcoords TMExact toise transverse TransverseMercatorExact TRF triangle triangulate triaxial triaxial-ellipsoid trigonometry trilaterate trilaterate-2d trilaterate-3d TwoProduct TwoSum umbilic-point unit unroll UPS UTM UTM/UPS variance velocities Veness Vermeille viewing-frustum Vincenty Visvalingam Visvalingam-Whyatt volume  volumetric Web-Mercator Welford WGRS WGS Whyatt Wildberger Witkam winding-number XYZ yaw You zenzi-cubic zenzi-quartic
+Keywords: AER Albers altitude Andoyer annulus antipode area attitude Authalic auxiliary azimuth azimuthal azimuth-elevation-range bearing bank Barsky Barth beta bi-quadratic boolean cached Cagnoli cartesian Cassini Cassini-Soldner chord circle-intersections circumcenter circumcircle circumradius clip Cohen Cohen-Sutherland Collins composite conformal conformal-sphere conic constants contact-triangle Cook Correia cosines-law coverage curvature cylindrical datum deprecation deficit development discrete distance Douglas earth east-north-up eccentricity ECEF elevation ellipsoid ellipsoidal-latitude-beta ellipsoidal-longitude-omega elliptic ENU EPSG equal-area equidistant equirectangular ETM ETRF Euclidean even-odd-rule ExactTM excess Farrell Farrell-Barth Ferrari-solution Field-Of-View flattening fma fmath footpoint footprint Forster Forster-Hormann-Popa Forsythe FOV fractional Frechet Fréchet frustum Fsum fused-multiply-add GARS geocentric GeoConvert GeodesicExact geodesy geodetic GeodSolve GeodTest geographiclib Geohash geoid geoidHeight GeoidHeights georef Girard gnomonic gons grades gradians Greiner Greiner-Hormann Hartzell Hausdorff Haversine heading hectare height Heikkinen Heron Hodgman horizon Hormann Hubeny IDW incenter incirle infix_@_operator inradius intermediate interpolate intersect intersection intersection3d intersections IntersectTool Inverse-Distance-Weighting Isometric ITRF Jacobi Jacobi-Conformal Jarque-Bera Jekel Karney Krueger Krüger kurtosis Lambert latitude law-of-cosines least-squares Lesh L_Huilier LHuilier Liang Liang-Barsky linearize Line-Of-Sight LocalCartesian local-tangent-plane local-x-y-z longitude LOS loxodrome lstsq LTP lune LV03 LV95 mean memoize memoized Mercator Meeus MGRS nearest NED Niemeyer non-finite normalize Norrdine north-east-down numpy n-vector Nvector oblate omega orthographic orthometric-height OSGB OSGR overlap parallel parallel-of-latitude Parametric path-intersection perimeter Peucker Pierlot pitch plumb Point-Of-View polar Popa POV precision-cubic-root precision-hypotenuse precision-powers precision-running-summation precision-square-root precision-summation prolate Pseudo-Mercator pygeodesy PyInstaller PyPy quartic radical radii radius Ramer Ramer-Douglas-Peucker Rectifying Reduced resect resection Rey-Jer Reumann Reumann-Witkam rho-theta-phi rhumb RhumbSolve running-linear-regression running-statistics running-stats running-summation scipy secant semi-perimeter sexagecimal simplify skewness Snellius Snellius-Pothenot Snyder Soddy Soddy-circles Soldner sphere sphere-intersections spherical-deficit spherical-excess spherical-polar spherical-triangle squared-quartic standard-deviation stereographic Sudano surface-area Sutherland Sutherland-Hodgman tangent-circles Terrestrial-Reference-Frame Thomas Tienstra tilt TMcoords TMExact toise transverse TransverseMercatorExact TRF triangle triangulate triaxial triaxial-ellipsoid trigonometry trilaterate trilaterate-2d trilaterate-3d TwoProduct TwoSum umbilic-point unit unroll UPS UTM UTM/UPS variance velocities Veness Vermeille viewing-frustum Vincenty Visvalingam Visvalingam-Whyatt volume  volumetric Web-Mercator Welford WGRS WGS Whyatt Wildberger Witkam winding-number XYZ yaw You zenzi-cubic zenzi-quartic
 Classifier: Development Status :: 5 - Production/Stable
 Classifier: Environment :: Console
 Classifier: Intended Audience :: Developers
@@ -16,7 +16,6 @@ Classifier: License :: OSI Approved :: MIT License
 Classifier: Operating System :: OS Independent
 Classifier: Programming Language :: Python
 Classifier: Programming Language :: Python :: 2.7
-Classifier: Programming Language :: Python :: 3.10
 Classifier: Programming Language :: Python :: 3.11
 Classifier: Programming Language :: Python :: 3.12
 Classifier: Programming Language :: Python :: 3.13
@@ -43,7 +42,7 @@ provides a geodetic LatLon_ and a geocentric Cartesian_ class with methods and f
 distance, surface area, perimeter, forward and reverse azimuth, initial and final bearing, intermediate
 and nearest points, intersections of geodesic, great circle and rhumb lines, circle intersections and
 secants, `3-point resections`_, triangulation, trilateration (by intersection, by overlap and in
-3d), conversions and unrolling, among other things.
+3d), among other things.
 
 Also included are modules for conversions to and from `Cassini-Soldner`_, ECEF_ (Earth-Centered,
 Earth-Fixed cartesian), UTM_ (Universal Transverse Mercator and Exact_), UPS_ (Universal Polar
@@ -60,8 +59,9 @@ operations between (composite) polygons, functions to simplify_ or linearize a p
 points (or a `numpy array`_), including implementations of the `Ramer-Douglas-Peucker`_,
 `Visvalingam-Whyatt`_ and `Reumann-Witkam`_ algorithms and modified versions of the former.
 
-Plus modules and classes to interpolate_ the Height_ of *LatLon* points and Geoid_ models, compute
-various Fréchet_ or Hausdorff_ distances or perform *boolean* operations between (composite) polygons.
+Plus modules and classes to interpolate_ the Height_ of *LatLon* points and Geoid_ models, to compute
+various Fréchet_ or Hausdorff_ distances or to perform *boolean* operations between (composite)
+polygons of *LatLon* points.
 
 For further details see the documentation_, the descriptions of `Latitude/Longitude`_, Vincenty_ and
 `Vector-based`_ geodesy, the original `JavaScript source`_ or docs_ and *Karney*\'s Python geographiclib_
@@ -120,22 +120,22 @@ test results (on macOS only) and the complete documentation_ generated by Epydoc
 Tests
 =====
 
-The tests ran with Python 3.13.7 (with geographiclib_ 2.1), 3.12.7 (with geographiclib_ 2.0, numpy_ 2.1.0,
-scipy_ 1.14.1, GeodSolve_ 2.5, IntersectTool_ 2.5 and RhumbSolve_ 2.5), 3.11.5 (with geographiclib_ 2.0,
-numpy_ 1.24.2 and scipy_ 1.10.1), Python 3.10.8 (with geographiclib_ 2.0, numpy_ 1.23.3, scipy_ 1.9.1,
-GeoConvert_ 2.5, GeodSolve_ 2.5), Python 3.9.6 and Python 2.7.18 (with geographiclib_ 1.50, numpy_ 1.16.6,
-scipy_ 1.2.2, GeoConvert_ 2.5, GeodSolve_ 2.5, IntersectTool_ 2.5 and RhumbSolve_ 2.5), all on macOS 15.6.1
-Sequoia in 64-bit.
+The tests ran with Python 3.14 (with geographiclib_ 2.1), Python 3.13.7 (with geographiclib_ 2.1,
+numpy_ 2.3.3, scipy_ 1.16.2, GeoConvert_ 2.5 and GeodSolve_ 2.5), Python 3.12.7 (with geographiclib_ 2.0,
+numpy_ 2.1.0, scipy_ 1.14.1, GeodSolve_ 2.5, IntersectTool_ 2.5 and RhumbSolve_ 2.5), Python 3.11.5 (with
+geographiclib_ 2.0, numpy_ 1.24.2 and scipy_ 1.10.1) and with Python 2.7.18 (with geographiclib_ 1.50,
+numpy_ 1.16.6, scipy_ 1.2.2, GeoConvert_ 2.5, GeodSolve_ 2.5, IntersectTool_ 2.5 and RhumbSolve_ 2.5),
+all on macOS 26.0.1 Tahoe in 64-bit.
 
 All tests ran with and without ``lazy import`` for Python 3 and with command line option ``-W default``
 and env variable ``PYGEODESY_WARNINGS=on`` for all Python versions.  The results of those tests are
 included in the distribution files.
 
-Python 3.13.7, 3.12.7, 3.11.5 and 3.10.8 run on Apple M4 Si (``arm64``), *natively*.  Python 2.7.18 runs
+Python 3.14, 3.13.7, 3.12.7 and 3.11.5 run on Apple Si M4 (``arm64``), *natively*.  Python 2.7.18 runs
 on Intel (``x86_64``) or Intel *emulation* (\"``arm64_x86_64``\", see function `pygeodesy.machine`_).
 
-Test coverage has been measured with coverage_ 7.6.1 using Python 3.13.7, 3.12.7, 3.11.5 and 3.10.8.  The
-complete coverage report in HTML and a PDF summary are included in the distribution files.
+Test coverage has been measured with coverage_ 7.10.7 using Python 3.14, 3.13.7 and 3.12.7.  The complete
+coverage report in HTML and a PDF summary are included in the distribution files.
 
 The tests also ran with Python 3.13.7 (and geographiclib_ 2.1) on `Debian 12`_ in 64-bit only and with
 Python 3.12.8 (and geographiclib_ 2.0) on `Windows 2019Server`_ in 64-bit only and with Python 2.7.18
@@ -144,13 +144,13 @@ Python 3.12.8 (and geographiclib_ 2.0) on `Windows 2019Server`_ in 64-bit only a
 A single-File and single-Directory application with ``pygeodesy`` has been bundled using PyInstaller_ 3.4
 and 64-bit Python 3.7.4 and 3.7.3 on macOS 10.13.6 High Sierra.
 
-Previously, the tests were run with Python 3.13.0-5, 3.12.0-6, 3.11.2-4, 3.10.1-7, 3.9.1, 3.8.7, 3.7.1, 2.7.15,
+Previously, the tests were run with Python 3.13.0-6, 3.12.0-6, 3.11.2-4, 3.10.1-7, 3.9.1, 3.8.7, 3.7.1, 2.7.15,
 PyPy_ 7.3.12 (Python 3.10.12), 7.3.1 (Python 3.6.9) and PyPy_ 7.1.1 (Python 2.7.13) (and geographiclib_ 1.52,
 numpy_ 1.16.3, 1.16.4, 1.16.6, 1.19.0, 1.19.4, 1.19.5 or 1.22.4 and scipy_ 1.2.1, 1.4.1, 1.5.2 or 1.8.1) on
 Ubuntu 16.04, with Python 3.10.0-1, 3.9.0-5, 3.8.0-6, 3.7.2-6, 3.7.0, 3.6.2-5, 3.5.3, 2.7.13-17, 2.7.10
 and 2.6.9 (and numpy_ 1.19.0, 1.16.5, 1.16.2, 1.15.2, 1.14.0, 1.13.1, 1.8.0rc1 or 1.6.2 and scipy_ 1.5.0),
 PyPy_ 7.3.0 (Python 2.7.13 and 3.6.9), PyPy_ 6.0.0 (Python 2.7.13 and 3.5.3) and `Intel-Python`_ 3.5.3 (and
-numpy_ 1.11.3) on macOS 15.0-5 Sequoia, 14.0-6.1 Sonoma, 13.0-5.2 Ventura, 12.1-6 Monterey, 11.0-5.2-6.1 Big
+numpy_ 1.11.3) on macOS 15.0-6 Sequoia, 14.0-6.1 Sonoma, 13.0-5.2 Ventura, 12.1-6 Monterey, 11.0-5.2-6.1 Big
 Sur (aka 10.16), 10.15.3, 10.15.5-7 Catalina, 10.14 Mojave, 10.13.6 High Sierra and 10.12 Sierra, MacOS X
 10.11 El Capitan and/or MacOS X 10.10 Yosemite, with Pythonista_ 3.2 (with geographiclib 1.50 or 1.49 and
 numpy 1.8.0) on iOS 14.4.2, 11.4.1, 12.0-3 on iPad4, iPhone6, iPhone10 and/or iPhone12, with Pythonista_ 3.1
@@ -160,12 +160,12 @@ Server 2012R2, Windows 10 Pro and 32-bit Python 2.6.6 on Windows XP SP3.
 Notes
 =====
 
-All Python source code has been statically checked_ with Ruff_ using Python 3.13.5 and with PyChecker_, PyFlakes_,
-PyCodeStyle_ (formerly Pep8) and McCabe_ using Python 2.7.18, both in 64-bit on macOS 15.6.1 Sequoia only.
+All Python source code has been statically checked_ with Ruff_ using Python 3.13.7 and with PyChecker_, PyFlakes_,
+PyCodeStyle_ (formerly Pep8) and McCabe_ using Python 2.7.18, both in 64-bit on macOS 26.0.1 Tahoe only.
 
 For a summary of all *Karney*-based functionality in ``pygeodesy``, see module karney_.
 
-*Last updated: Sep 09, 2025.*
+*Last updated: Nov 05, 2025.*
 
 License
 =======