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

pygeodesy/__init__.py

@@ -17,7 +17,7 @@ and a geocentric B{C{Cartesian}} class with methods and functions to compute dis
 forward and reverse azimuth, initial and final bearing, intermediate and nearest points, intersections of geodesic,
 great circle and rhumb lines, circle intersections and secants, U{3-point resections
 <https://WikiPedia.org/wiki/Position_resection_and_intersection>}, triangulation, trilateration (by intersection,
-by overlap and in 3-D), conversions and unrolling, among other things.
+by overlap and in 3-D), among other things.
 
 Also included are modules for conversions to and from U{Cassini-Soldner
 <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1CassiniSoldner.html>},
@@ -51,9 +51,9 @@ resources/hull:8338>} and the U{Reumann-Witkam<https://psimpl.SourceForge.net/re
 algorithms and modified versions of the former.
 
 Plus modules and classes to U{interpolate<https://docs.SciPy.org/doc/scipy/reference/interpolate.html>} the
-L{height<pygeodesy.heights>} of C{LatLon} points and L{Geoid<pygeodesy.geoids>} models, compute various U{Fréchet
+L{height<pygeodesy.heights>} of C{LatLon} points and L{Geoid<pygeodesy.geoids>} models, to compute various U{Fréchet
 <https://WikiPedia.org/wiki/Frechet_distance>} or U{Hausdorff<https://WikiPedia.org/wiki/Hausdorff_distance>}
-distances or perform I{boolean} operations between (composite) polygons.
+distances or to perform I{boolean} operations between (composite) polygons of C{LatLon} points.
 
 For further details see the U{documentation<https://mrJean1.GitHub.io/PyGeodesy>}, the descriptions of
 U{Latitude/Longitude<https://www.Movable-Type.co.UK/scripts/latlong.html>}, U{Vincenty
@@ -125,7 +125,11 @@ C{epydoc --html --no-private --no-source --name=pygeodesy --url=... -v pygeodesy
 Tests
 =====
 
-The tests ran with Python 3.13.7 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.1),
+The tests ran with Python 3.14 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.1),
+Python 3.13.7 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.1),
+U{numpy<https://PyPI.org/project/numpy>} 2.3.3, U{scipy<https://PyPI.org/project/scipy>} 1.16.2,
+U{GeoConvert<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5 and
+U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5),
 Python 3.12.7 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0,
 U{numpy<https://PyPI.org/project/numpy>} 2.1.0, U{scipy<https://PyPI.org/project/scipy>} 1.14.1,
 U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5,
@@ -133,27 +137,23 @@ U{IntersectTool<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.
 U{RhumbSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5),
 Python 3.11.5 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0,
 U{numpy<https://PyPI.org/project/numpy>} 1.24.2 and U{scipy<https://PyPI.org/project/scipy>} 1.10.1),
-Python 3.10.8 (with U{geographiclib <https://PyPI.org/project/geographiclib>} 2.0,
-U{numpy<https://PyPI.org/project/numpy>} 1.23.3, U{scipy<https://PyPI.org/project/scipy>} 1.9.1,
-U{GeoConvert<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5,
-U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5) and
-Python 2.7.18 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 1.50,
+and with Python 2.7.18 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 1.50,
 U{numpy<https://PyPI.org/project/numpy>} 1.16.6, U{scipy<https://PyPI.org/project/scipy>} 1.2.2,
 U{GeoConvert<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5,
 U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5,
 U{IntersectTool<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5 and
 U{RhumbSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5), all in 64-bit on
-macOS 15.6.1 Sequoia.
+macOS 26.0.1 Tahoe.
 
 All tests ran with and without C{lazy import} for Python 3 and with command line option C{-W default} and
 env variable C{PYGEODESY_WARNINGS=on} for all Python versions.  The results of those tests are included in
 the distribution files.
 
-Test coverage has been measured with U{coverage<https://PyPI.org/project/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 U{coverage<https://PyPI.org/project/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.
 
-Python 3.13.7, 3.12.7, 3.11.5 and 3.10.8 run on Apple M4 Si (C{arm64}), I{natively}.  Python 2.7.18 runs
+Python 3.14, 3.13.7, 3.12.7 and 3.11.5 run on Apple M4 Si (C{arm64}), I{natively}.  Python 2.7.18 runs
 on Intel (C{x86_64}) or Intel I{emulation} ("C{arm64_x86_64}", see function L{machine<pygeodesy.machine>}).
 
 The tests also ran with Python 3.13.7 (and U{geographiclib<https://PyPI.org/project/geographiclib>} 2.1) on
@@ -166,7 +166,7 @@ in 64- and 32-bit.
 A single-File and single-Directory application with C{pygeodesy} has been bundled using U{PyInstaller
 <https://PyPI.org/project/pyinstaller>} 3.4 and 64-bit Python 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.6, 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.6, 3.9.1, 3.8.7, 3.7.1, 2.7.15,
 U{PyPy<https://PyPy.org>} 7.3.12 (Python 3.10.12), 7.3.1 (Python 3.6.9) and U{PyPy<https://PyPy.org>} 7.1.1 (Python
 2.7.13) (and U{geographiclib <https://PyPI.org/project/geographiclib>} 1.52, U{numpy<https://PyPI.org/project/numpy>}
 1.16.3, 1.16.4, 1.16.6, 1.19.0, 1.19.4, 1.19.5 or 1.22.4 and U{scipy<https://PyPI.org/project/scipy>} 1.2.1, 1.4.1,
@@ -175,7 +175,7 @@ U{PyPy<https://PyPy.org>} 7.3.12 (Python 3.10.12), 7.3.1 (Python 3.6.9) and U{Py
 1.16.5, 1.16.2, 1.15.2, 1.14.0, 1.13.1, 1.8.0rc1 or 1.6.2 and U{scipy<https://PyPI.org/project/scipy>} 1.5.0), U{PyPy
 <https://PyPy.org>} 7.3.0 (Python 2.7.13 and 3.6.9), U{PyPy<https://PyPy.org>} 6.0.0 (Python 2.7.13 and 3.5.3)
 and U{Intel-Python<https://software.Intel.com/en-us/distribution-for-python>} 3.5.3 (and U{numpy
-<https://PyPI.org/project/numpy>} 1.11.3) on macOS 15.0-5 Sequoia, 14.0-6.1 Sonoma, 13.0-5.2 Ventura, 12.1-6 Monterey,
+<https://PyPI.org/project/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 U{Pythonista<https://OMZ-Software.com/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
@@ -187,10 +187,10 @@ Notes
 =====
 
 All Python source code has been statically U{checked<https://GitHub.com/ActiveState/code/tree/master/recipes/Python/
-546532_PyChecker_postprocessor>} with U{Ruff<https://GitHub.com/astral-sh/ruff>} using Python 3.13.5 and with
+546532_PyChecker_postprocessor>} with U{Ruff<https://GitHub.com/astral-sh/ruff>} using Python 3.13.7 and with
 U{PyChecker<https://PyPI.org/project/pychecker>}, U{PyFlakes<https://PyPI.org/project/pyflakes>}, U{PyCodeStyle
 <https://PyPI.org/project/pycodestyle>} (formerly Pep8) and U{McCabe<https://PyPI.org/project/mccabe>} using Python
-2.7.18, both in 64-bit on macOS 15.6.1 Sequoia.
+2.7.18, both in 64-bit on macOS 26.0.1 Tahoe.
 
 For a summary of all I{Karney}-based functionality in C{pygeodesy}, see module U{karney
 <https://mrJean1.GitHub.io/PyGeodesy/docs/pygeodesy.karney-module.html>}.
@@ -605,7 +605,7 @@ else:
 
 from pygeodesy.internals import _version2,  _DOT_  # noqa: E402
 # from pygeodesy.interns import _DOT_  # from .internals
-__version__ = '25.09.09'
+__version__ = '25.10.10'
 # see setup.py for similar logic
 version     = _DOT_(*_version2(__version__, n=3))
 

pygeodesy/azimuthal.py

@@ -49,7 +49,7 @@ from pygeodesy.ellipsoidalBase import LatLonEllipsoidalBase as _LLEB, \
                                      _isin, _xinstanceof
 from pygeodesy.datums import _spherical_datum, _WGS84
 from pygeodesy.errors import _ValueError, _xdatum, _xkwds
-from pygeodesy.fmath import euclid, fdot_, hypot as _hypot,  Fsum
+from pygeodesy.fmath import euclid, hypot as _hypot,  Fsum
 # from pygeodesy.fsums import Fsum  # from .fmath
 # from pygeodesy.formy import antipode  # _MODS
 # from pygeodesy.internals import typename  # from .karney
@@ -71,7 +71,7 @@ from pygeodesy.utily import asin1, atan1, atan2, atan2b, atan2d, \
 from math import acos, degrees, fabs, sin, sqrt
 
 __all__ = _ALL_LAZY.azimuthal
-__version__ = '25.05.12'
+__version__ = '25.09.29'
 
 _EPS_K         = _EPStol * _0_1  # Karney's eps_ or _EPSmin * _0_1?
 _over_horizon_ = 'over horizon'
@@ -147,7 +147,7 @@ class _AzimuthalBase(_NamedBase):
         '''I{Must be overloaded}.'''
         self._notOverloaded(lat, lon, **name)
 
-    def _forward(self, lat, lon, name, _k_t_2):
+    def _forward(self, lat, lon, name, _k_t):
         '''(INTERNAL) Azimuthal (spherical) forward C{lat, lon} to C{x, y}.
         '''
         lat, lon = Lat_(lat), Lon_(lon)
@@ -155,16 +155,16 @@ class _AzimuthalBase(_NamedBase):
         s0, c0 = self._sc0
 
         cb  *= ca
-        k, t = _k_t_2(fdot_(s0, sa, c0, cb))
+        k, t = _k_t(c0 * cb + s0 * sa)
         if t:
             r = k * self.radius
-            y = r * fdot_(c0, sa, -s0, cb)
+            y = r * (c0 * sa - s0 * cb)
             e, n, z, _ = _enzh4(r * sb * ca, y, None)
         else:  # 0 or 180
             e = n = z = _0_0
 
         t = Azimuthal7Tuple(e, n, lat, lon, z, k, self.datum,
-                                       name=self._name__(name))
+                                                  name=self._name__(name))
         return t
 
     def _forwards(self, *lls):
@@ -258,7 +258,7 @@ class _AzimuthalBase(_NamedBase):
             if lea or fabs(c0) > EPS:
                 d = atan2d(e * sc, r * c0 * cc - n * s0 * sc)
             else:
-                d = atan2d(e, (n if s0 < 0 else -n))
+                d = atan2d(e,     (n if s0 < 0 else -n))
             lon = _norm180(self.lon0 + d)
 
         if LatLon is None:
@@ -745,8 +745,8 @@ class _GnomonicBase(_AzimuthalGeodesic):
 
            @arg lat: Latitude of the location (C{degrees90}).
            @arg lon: Longitude of the location (C{degrees180}).
-           @kwarg raiser: Do or don't throw an error (C{bool}) if
-                          the location lies over the horizon.
+           @kwarg raiser: If C{False}, do not throw an error if the location lies
+                          over the horizon (C{bool}).
            @kwarg name: Optional C{B{name}=NN} for the location (C{str}).
 
            @return: An L{Azimuthal7Tuple}C{(x, y, lat, lon, azimuth, scale, datum)}
@@ -760,8 +760,6 @@ class _GnomonicBase(_AzimuthalGeodesic):
            @raise AzimuthalError: Invalid B{C{lat}}, B{C{lon}} or the location lies
                                   over the horizon and C{B{raiser}=True}.
         '''
-        self._iteration = 0
-
         r = self.geodesic.Inverse(self.lat0, self.lon0,
                                   Lat_(lat), Lon_(lon), outmask=self._mask)
         M = r.M21
@@ -776,7 +774,7 @@ class _GnomonicBase(_AzimuthalGeodesic):
             e = n = NAN
 
         t = self._7Tuple(e, n, r, _name__(name), M=M)
-        t._iteraton = 0
+        t._iteration = self._iteration = 0
         return t
 
     def reverse(self, x, y, LatLon=None, **name_LatLon_kwds):  # PYCHOK signature

pygeodesy/constants.py

@@ -26,7 +26,7 @@ except ImportError:  # Python 2-
     _inf, _nan = float(_INF_), float(_NAN_)
 
 __all__ = _ALL_LAZY.constants
-__version__ = '25.09.09'
+__version__ = '25.09.11'
 
 
 def _copysign_0_0(y):
@@ -186,6 +186,7 @@ _0_1     = _float(   0.1)     # PYCHOK expected
 _0_125   = _float(   0.125)   # PYCHOK expected
 _0_25    = _float(   0.25)    # PYCHOK expected
 _0_5     = _float(   0.5)     # PYCHOK expected
+_0_75    = _float(   0.75)     # PYCHOK expected
 _1_0     = _float(   1)       # PYCHOK expected
 _1_0_1T  = _1_0,              # PYCHOK 1-tuple
 _1_5     = _float(   1.5)     # PYCHOK expected
@@ -236,40 +237,35 @@ try:
 #   RADIX    =  Int(  RADIX   =_f_i.radix)     # PYTHON system's float base
     del _f_i
 except ImportError:  # PYCHOK no cover
-    DIG      =  Int(  DIG     =15)          # PYCHOK system's 64-bit float decimal digits
+    DIG      =  Int(  DIG     =15)           # PYCHOK system's 64-bit float decimal digits
     EPS      = _Float(EPS     =2.220446049250313e-16)  # PYCHOK EPSilon 2**-52, M{EPS +/- 1 != 1}
-    MANT_DIG =  Int(  MANT_DIG=53)          # PYCHOK float mantissa bits ≈ 53 (C{int})
+    MANT_DIG =  Int(  MANT_DIG=53)           # PYCHOK float mantissa bits ≈ 53 (C{int})
     MAX      = _Float(MAX     =pow(_2_0,  1023) * (_2_0 - EPS))  # PYCHOK ≈ 10**308
-    MAX_EXP  =  Int(  MAX_ESP =_log2(MAX))  # 308 base 10
+    MAX_EXP  =  Int(  MAX_ESP =_log2(MAX))   # 308 base 10
     MIN      = _Float(MIN     =pow(_2_0, -1021))  # PYCHOK ≈ 10**-308
-    MIN_EXP  =  Int(MIN_EXP   =_log2(MIN))  # -307 base 10
-#   RADIX    =  Int(Radix     =2)           # base
-
-EPS0     = _Float( EPS0  = EPS**2)          # PYCHOK near-/non-zero comparison 4.930381e-32, or EPS or EPS_2
-EPS02    = _Float( EPS02 = EPS**4)          # PYCHOK near-zero-squared comparison 2.430865e-63
-EPS_2    = _Float( EPS_2 = EPS / _2_0)      # PYCHOK ≈ 1.110223024625e-16
-EPS1     = _Float( EPS1  =_1_0 - EPS)       # PYCHOK ≈ 0.9999999999999998
-EPS2     = _Float( EPS2  = EPS * _2_0)      # PYCHOK ≈ 4.440892098501e-16
-EPS4     = _Float( EPS4  = EPS * _4_0)      # PYCHOK ≈ 8.881784197001e-16
-# _1EPS  = _Float(_1EPS  =_1_0 + EPS)       # PYCHOK ≈ 1.0000000000000002
-_1_EPS   = _Float(_1_EPS =_1_0 / EPS)       # PYCHOK = 4503599627370496.0
-# _2_EPS = _Float(_2_EPS =_2_0 / EPS)       # PYCHOK = 9007199254740992.0
-_EPS2e4  = _Float(_EPS2e4= EPS2 * 1.e4)     # PYCHOK ≈ 4.440892098501e-12
-_EPS4e8  = _Float(_EPS4e8= EPS4 * 1.e8)     # PYCHOK ≈ 8.881784197001e-08
-_EPSjam  = _Float(_EPSjam= pow(EPS, 0.75))  # PYCHOK = 1.818989403546e-12
-_EPSmin  = _Float(_EPSmin= sqrt(MIN))       # PYCHOK = 1.49166814624e-154
-_EPSqrt  = _Float(_EPSqrt= sqrt(EPS))       # PYCHOK = 1.49011611938e5-08
-_EPStol  = _Float(_EPStol=_EPSqrt * _0_1)   # PYCHOK = 1.49011611938e5-09 == sqrt(EPS * _0_01)
+    MIN_EXP  =  Int(MIN_EXP   =_log2(MIN))   # -307 base 10
+#   RADIX    =  Int(Radix     =2)            # base
+
+EPS0     = _Float( EPS0  = EPS**2)           # PYCHOK near-/non-zero comparison 4.930381e-32, or EPS or EPS_2
+EPS02    = _Float( EPS02 = EPS**4)           # PYCHOK near-zero-squared comparison 2.430865e-63
+EPS_2    = _Float( EPS_2 = EPS / _2_0)       # PYCHOK ≈ 1.110223024625e-16
+EPS1     = _Float( EPS1  =_1_0 - EPS)        # PYCHOK ≈ 0.9999999999999998
+EPS2     = _Float( EPS2  = EPS * _2_0)       # PYCHOK ≈ 4.440892098501e-16
+EPS4     = _Float( EPS4  = EPS * _4_0)       # PYCHOK ≈ 8.881784197001e-16
+# _1EPS  = _Float(_1EPS  =_1_0 + EPS)        # PYCHOK ≈ 1.0000000000000002
+_1_EPS   = _Float(_1_EPS =_1_0 / EPS)        # PYCHOK = 4503599627370496.0
+# _2_EPS = _Float(_2_EPS =_2_0 / EPS)        # PYCHOK = 9007199254740992.0
+_EPS2e4  = _Float(_EPS2e4= EPS2 * 1.e4)      # PYCHOK ≈ 4.440892098501e-12
+_EPS4e8  = _Float(_EPS4e8= EPS4 * 1.e8)      # PYCHOK ≈ 8.881784197001e-08
+_EPSjam  = _Float(_EPSjam= pow(EPS, _0_75))  # PYCHOK = 1.818989403546e-12
+_EPSmin  = _Float(_EPSmin= sqrt(MIN))        # PYCHOK = 1.49166814624e-154
+_EPSqrt  = _Float(_EPSqrt= sqrt(EPS))        # PYCHOK = 1.49011611938e5-08
+_EPStol  = _Float(_EPStol=_EPSqrt * _0_1)    # PYCHOK = 1.49011611938e5-09 == sqrt(EPS * _0_01)
 
 _89_999  = _Float(_89_999=_90_0 * EPS1)  # just below 90.0
 # <https://Numbers.Computation.Free.FR/Constants/Miscellaneous/digits.html>
 # _1__90 = _Float(_1__90 =_1_0 / _90_0)  # PYCHOK = 0.011_111_111_111_111_111_111_111_111_111_111_111_111_111_111_11111
 _2__PI   = _Float(_2__PI =_2_0 / _pi)    # PYCHOK = 0.636_619_772_367_581_343_075_535_053_490_057_448_137_838_582_96182
-
-_1_16th  = _Float(_1_16th=_1_0 / _16_0)  # PYCHOK in .ellipsoids, .karney
-_1_3rd   = _Float(_1_3rd =_1_0 /  _3_0)  # PYCHOK in .fmath
-_1_6th   = _Float(_1_6th =_1_0 /  _6_0)  # PYCHOK in .fmath
-
 _K0_UTM  = _Float(_K0_UTM = 0.9996)  # PYCHOK in .etm, .ktm, .utm, UTM scale at central meridian
 # sqrt(2) <https://WikiPedia.org/wiki/Square_root_of_2>
 # 1.414213562373095_048_801_688_724_209_698_078_569_671_875_376_948_073_176_679_737_99

pygeodesy/elliptic.py

@@ -100,7 +100,7 @@ from math import asin, asinh, atan, ceil, cosh, fabs, floor, radians, \
                  sin, sinh, sqrt, tan, tanh  # tan as _tan
 
 __all__ = _ALL_LAZY.elliptic
-__version__ = '25.09.04'
+__version__ = '25.10.10'
 
 _TolRD  =  zqrt(EPS * 0.002)
 _TolRF  =  zqrt(EPS * 0.030)
@@ -932,7 +932,7 @@ class Elliptic(_Named):
             raise _ellipticError(Elliptic._RFRD, x, y, z, m, cause=X)
         return float(R)
 
-_allPropertiesOf_n(16, Elliptic)  #  # PYCHOK assert, see Elliptic.reset
+_allPropertiesOf_n(16, Elliptic)  # PYCHOK assert, see Elliptic.reset
 
 
 class EllipticError(_ValueError):

pygeodesy/etm.py

@@ -67,7 +67,7 @@ from pygeodesy.basics import _isin, map1, neg, neg_, _xinstanceof
 from pygeodesy.constants import EPS, EPS02, PI_2, PI_4, _K0_UTM, \
                              _1_EPS, _0_0, _0_1, _0_5, _1_0, _2_0, \
                              _3_0, _90_0, isnear0, isnear90
-from pygeodesy.constants import _4_0  # PYCHOK used!
+from pygeodesy.constants import _0_75, _4_0  # PYCHOK used!
 from pygeodesy.datums import _ellipsoidal_datum, _WGS84,  _EWGS84
 # from pygeodesy.ellipsoids import _EWGS84  # from .datums
 # from pygeodesy.elliptic import Elliptic  # _MODS
@@ -93,7 +93,7 @@ from pygeodesy.utm import _cmlon, _LLEB, _parseUTM5, _toBand, _toXtm8, \
 from math import asinh, degrees, radians, sinh, sqrt
 
 __all__ = _ALL_LAZY.etm
-__version__ = '25.08.31'
+__version__ = '25.09.15'
 
 _OVERFLOW = _1_EPS**2  # ~2e+31
 _TAYTOL   =  pow(EPS,  0.6)
@@ -410,7 +410,7 @@ class ExactTransverseMercator(_NamedBase):
     def _Ev_3cKE_4(self):
         '''(INTERNAL) Get and cache C{_Ev.cKE * 3 / 4}.
         '''
-        return self._Ev_cKE * 0.75  # _0_75
+        return self._Ev_cKE * _0_75
 
     @Property_RO
     def _Ev_5cKE_4(self):

pygeodesy/fmath.py

@@ -10,8 +10,8 @@ from __future__ import division as _; del _  # noqa: E702 ;
 from pygeodesy.basics import _copysign, copysign0, isbool, isint, isodd, \
                               isscalar, len2, map1, _xiterable,  typename
 from pygeodesy.constants import EPS0, EPS02, EPS1, NAN, PI, PI_2, PI_4, \
-                               _0_0, _0_125, _1_6th, _0_25, _1_3rd, _0_5, _1_0, \
-                               _1_5, _copysign_0_0, isfinite, remainder
+                               _0_0, _0_125, _0_25, _0_5, _1_0, _1_5, \
+                               _copysign_0_0, isfinite, remainder
 from pygeodesy.errors import _IsnotError, LenError, _TypeError, _ValueError, \
                              _xError, _xkwds, _xkwds_pop2, _xsError
 from pygeodesy.fsums import _2float, Fsum, fsum, _isFsum_2Tuple,  Fmt, unstr
@@ -25,10 +25,12 @@ from math import fabs, sqrt  # pow
 import operator as _operator  # in .datums, .trf, .utm
 
 __all__ = _ALL_LAZY.fmath
-__version__ = '25.08.31'
+__version__ = '25.09.15'
 
 # sqrt(2) - 1 <https://WikiPedia.org/wiki/Square_root_of_2>
 _0_4142  =  0.41421356237309504880  # ~ 3_730_904_090_310_553 / 9_007_199_254_740_992
+_1_3rd   = _1_0 / 3
+_1_6th   = _1_0 / 6
 _2_3rd   = _1_3rd * 2
 _h_lt_b_ = 'abs(h) < abs(b)'
 

pygeodesy/geodesicx/gx.py

@@ -65,7 +65,7 @@ from pygeodesy.utily import atan2, atan2d as _atan2d_reverse, _unrollon, \
 from math import copysign, cos, degrees, fabs, radians, sqrt
 
 __all__ = ()
-__version__ = '25.09.02'
+__version__ = '25.09.16'
 
 _MAXIT1 = 20
 _MAXIT2 = 10 + _MAXIT1 + MANT_DIG  # MANT_DIG == C++ digits
@@ -188,7 +188,7 @@ class GeodesicExact(_GeodesicBase):
                  <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1GeodesicExact.html>}
                  and Python U{Geodesic.ArcDirect<https://GeographicLib.SourceForge.io/Python/doc/code.html>}.
         '''
-        return self._GDictDirect(lat1, lon1, azi1, True, a12, outmask)
+        return self._GDictDirect(lat1, lon1, azi1, True, a12, outmask=outmask)
 
     def ArcDirectLine(self, lat1, lon1, azi1, a12, caps=Caps.ALL, **name):
         '''Define a L{GeodesicLineExact} in terms of the I{direct} geodesic problem and as arc length.
@@ -348,7 +348,7 @@ class GeodesicExact(_GeodesicBase):
                  <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1GeodesicExact.html>}
                  and Python U{Geodesic.Direct<https://GeographicLib.SourceForge.io/Python/doc/code.html>}.
         '''
-        return self._GDictDirect(lat1, lon1, azi1, False, s12, outmask)
+        return self._GDictDirect(lat1, lon1, azi1, False, s12, outmask=outmask)
 
     def Direct3(self, lat1, lon1, azi1, s12):  # PYCHOK outmask
         '''Return the destination lat, lon and reverse azimuth
@@ -356,7 +356,7 @@ class GeodesicExact(_GeodesicBase):
 
            @return: L{Destination3Tuple}C{(lat, lon, final)}.
         '''
-        r = self._GDictDirect(lat1, lon1, azi1, False, s12, Caps._AZIMUTH_LATITUDE_LONGITUDE)
+        r = self._GDictDirect(lat1, lon1, azi1, False, s12, outmask=Caps._AZIMUTH_LATITUDE_LONGITUDE)
         return Destination3Tuple(r.lat2, r.lon2, r.azi2)  # no iteration
 
     def _DirectLine(self, ll1, azi12, s12=0, **caps_name):
@@ -710,7 +710,7 @@ class GeodesicExact(_GeodesicBase):
            @return: L{Direct9Tuple}C{(a12, lat2, lon2, azi2,
                                       s12, m12,  M12,  M21, S12)}.
         '''
-        r = self._GDictDirect(lat1, lon1, azi1, arcmode, s12_a12, outmask)
+        r = self._GDictDirect(lat1, lon1, azi1, arcmode, s12_a12, outmask=outmask)
         return r.toDirect9Tuple()
 
     def _GenInverse(self, lat1, lon1, lat2, lon2, outmask=Caps.STANDARD):
@@ -719,7 +719,7 @@ class GeodesicExact(_GeodesicBase):
            @return: L{Inverse10Tuple}C{(a12, s12, salp1, calp1, salp2, calp2,
                                              m12,   M12,   M21,   S12)}.
         '''
-        r = self._GDictInverse(lat1, lon1, lat2, lon2, outmask | Caps._SALP_CALPs_)
+        r = self._GDictInverse(lat1, lon1, lat2, lon2, outmask=outmask | Caps._SALP_CALPs_)
         return r.toInverse10Tuple()
 
     def _Inverse(self, ll1, ll2, wrap, **outmask):
@@ -752,7 +752,7 @@ class GeodesicExact(_GeodesicBase):
                  <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1GeodesicExact.html>} and
                  Python U{Geodesic.InverseLine<https://GeographicLib.SourceForge.io/Python/doc/code.html>}.
         '''
-        return self._GDictInverse(lat1, lon1, lat2, lon2, outmask)
+        return self._GDictInverse(lat1, lon1, lat2, lon2, outmask=outmask)
 
     def Inverse1(self, lat1, lon1, lat2, lon2, wrap=False):
         '''Return the non-negative, I{angular} distance in C{degrees}.
@@ -764,7 +764,8 @@ class GeodesicExact(_GeodesicBase):
         # and .HeightIDWkarney._distances
         if wrap:
             _, lat2, lon2 = _Wrap.latlon3(lat1, lat2, lon2, True)  # _Geodesic.LONG_UNROLL
-        return fabs(self._GDictInverse(lat1, lon1, lat2, lon2, Caps.EMPTY).a12)  # a12 always
+        r = self._GDictInverse(lat1, lon1, lat2, lon2, outmask=Caps.EMPTY)
+        return fabs(r.a12)  # a12 always
 
     def Inverse3(self, lat1, lon1, lat2, lon2):  # PYCHOK outmask
         '''Return the distance in C{meter} and the forward and
@@ -772,7 +773,7 @@ class GeodesicExact(_GeodesicBase):
 
            @return: L{Distance3Tuple}C{(distance, initial, final)}.
         '''
-        r = self._GDictInverse(lat1, lon1, lat2, lon2, Caps.AZIMUTH_DISTANCE)
+        r = self._GDictInverse(lat1, lon1, lat2, lon2, outmask=Caps.AZIMUTH_DISTANCE)
         return Distance3Tuple(r.s12, wrap360(r.azi1), wrap360(r.azi2),
                               iteration=r.iteration)
 
@@ -804,7 +805,7 @@ class GeodesicExact(_GeodesicBase):
                  <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1GeodesicExact.html>} and
                  Python U{Geodesic.InverseLine<https://GeographicLib.SourceForge.io/Python/doc/code.html>}.
         '''
-        r = self._GDictInverse(lat1, lon1, lat2, lon2, caps | Caps._SALP_CALPs_)
+        r = self._GDictInverse(lat1, lon1, lat2, lon2, outmask=caps | Caps._SALP_CALPs_)
         return GeodesicLineExact(self, lat1, lon1, None, caps=caps, _s_calp1=(r.salp1, r.calp1),
                                                             **name)._GenSet(self._debug, **r)
 

pygeodesy/geodesicx/gxarea.py

@@ -31,7 +31,7 @@ from pygeodesy.props import Property, Property_RO, property_RO
 from math import fmod as _fmod
 
 __all__ = ()
-__version__ = '25.06.04'
+__version__ = '25.09.16'
 
 
 class GeodesicAreaExact(_NamedBase):
@@ -74,7 +74,7 @@ class GeodesicAreaExact(_NamedBase):
         self._g_gX = g = geodesic
         # use the class-level Caps since the values
         # differ between GeodesicExact and Geodesic
-        self._mask = g.DISTANCE | g.LATITUDE | g.LONGITUDE
+        self._mask =  g.DISTANCE | g.LATITUDE | g.LONGITUDE
         self._Peri = _Accumulator(name='_Peri')
         if not polyline:  # perimeter and area
             self._mask |=  g.AREA | g.LONG_UNROLL
@@ -188,7 +188,7 @@ class GeodesicAreaExact(_NamedBase):
         '''(INTERNAL) Edge helper.
         '''
         lon1 = self.lon1
-        r = self._g_gX._GDictDirect(self.lat1, lon1, azi, False, s, self._mask)
+        r = self._g_gX._GDictDirect(self.lat1, lon1, azi, False, s, outmask=self._mask)
         if self._Area:  # aka transitDirect
             # Count crossings of prime meridian exactly as
             # int(ceil(lon2 / 360)) - int(ceil(lon1 / 360))
@@ -219,7 +219,7 @@ class GeodesicAreaExact(_NamedBase):
     def _Inverse(self, lat1, lon1, lat2, lon2):
         '''(INTERNAL) Point helper.
         '''
-        r = self._g_gX._GDictInverse(lat1, lon1, lat2, lon2, self._mask)
+        r = self._g_gX._GDictInverse(lat1, lon1, lat2, lon2, outmask=self._mask)
         if self._Area:  # aka transit
             # count crossings of prime meridian as +1 or -1
             # if in east or west direction, otherwise 0

pygeodesy/geohash.py

@@ -37,7 +37,7 @@ from pygeodesy.units import Degrees_, Int, Lat_, Lon_, Meter, Precision_, Str
 from math import fabs, ldexp, log10, radians
 
 __all__ = _ALL_LAZY.geohash
-__version__ = '25.04.21'
+__version__ = '25.09.16'
 
 _formy   = _MODS.into(formy=__name__)
 _MASK5   =  16, 8, 4, 2, 1  # PYCHOK used!
@@ -150,10 +150,10 @@ class _GH(object):
         '''
         def _encodes(lat, lon, prec, eps=0):
             s, w, n, e = self.SWNE4
-            E, d, _mid = self.EncodeB32, True, _2mid
+            d, _mid = True, _2mid
             for _ in range(prec):
                 i = 0
-                for _ in range(5):  # len(_MASK5)
+                for _ in _MASK5:
                     i += i
                     if d:  # bisect longitude
                         a = _mid(e, w)
@@ -170,7 +170,7 @@ class _GH(object):
                             s  = a
                             i += 1
                     d = not d
-                yield E[i]
+                yield self.EncodeB32[i]
                 if eps > 0:  # infer prec
                     if _2dab(lon, e, w) < eps and \
                        _2dab(lat, n, s) < eps:
@@ -233,10 +233,10 @@ class _GH(object):
         if not (0 < nc <= _MaxPrec):  # or geohash.startswith(_INV_)
             raise GeohashError(geohash=geohash, len=nc)
         s, w, n, e = self.SWNE4
-        D, d, _mid = self.DecodeB32, True, _2mid
+        d, _mid = True, _2mid
         try:
             for j, c in enumerate(geohash.lower()):
-                i = D[c]
+                i = self.DecodeB32[c]
                 for m in mask5:
                     if d:  # longitude
                         a = _mid(e, w)

pygeodesy/geoids.py

@@ -28,7 +28,7 @@ or L{GeoidPGM} and download a C{geoid} model file, containing locations with
 known heights also referred to as the C{grid knots}.  See the documentation of
 the interpolator class for references to available C{grid} models.
 
-C{>>> from pygeodesy import GeoidEGM96  # or -G2012B, -Karney or -PGM as GeoidXyz}
+C{>>> from pygeodesy import GeoidEGM96 as GeoidXyz  # or GeoidG2012B, -Karney or -PGM}
 
 2. Instantiate an interpolator with the C{geoid} model file and use keyword
 arguments to select different interpolation options
@@ -126,7 +126,7 @@ except ImportError:  # Python 3+
     from io import BytesIO as _BytesIO  # PYCHOK expected
 
 __all__ = _ALL_LAZY.geoids
-__version__ = '25.05.21'
+__version__ = '25.09.26'
 
 _assert_    = 'assert'
 _bHASH_     = b'#'

pygeodesy/heights.py

@@ -21,7 +21,7 @@ C{>>> ...}
 
 2. Select one of the C{Height} classes for height interpolation
 
-C{>>> from pygeodesy import HeightCubic  # or other Height... as HeightXyz}
+C{>>> from pygeodesy import HeightCubic as HeightXyz  # or an other Height... class}
 
 3. Instantiate a height interpolator with the C{knots} and use keyword
 arguments to select different interpolation options
@@ -92,7 +92,7 @@ from pygeodesy.units import _isDegrees, Float_, Int_
 # from math import radians  # from .points
 
 __all__ = _ALL_LAZY.heights
-__version__ = '25.05.26'
+__version__ = '25.09.29'
 
 _error_  = 'error'
 _formy   = _MODS.into(formy=__name__)

pygeodesy/internals.py

@@ -202,7 +202,7 @@ class _MODS_Base(object):
                      (_macOS_,   p.mac_ver),
                      (_Windows_, p.win32_ver),
                      (_Nix,  _MODS.nix_ver),
-                     ('Java',    p.java_ver),
+# removed Py 3.15    ('Java',    p.java_ver),
                      ('uname',   p.uname)):
             v = v()[0]
             if v and n:
@@ -659,11 +659,10 @@ def _usage_argv(argv0, *args):
     '''(INTERNAL) Return 3-tuple C{(python, '-m', module, *args)}.
     '''
     o = _MODS.os
-    m =  o.path.dirname(argv0)
-    m =  m.replace(o.getcwd(), _ELLIPSIS_) \
-          .replace(o.sep,      _DOT_).strip()
-    b =  o.path.basename(argv0)
-    b, x = o.path.splitext(b)
+    p =  o.path
+    m =  p.dirname(argv0).replace(o.getcwd(), _ELLIPSIS_) \
+                         .replace(o.sep,      _DOT_).strip()
+    b, x = p.splitext(p.basename(argv0))
     if x == '.py' and b != _DMAIN_:
         m = _DOT_(m or _pygeodesy_, b)
     p = NN(_python_, _MODS.sys_version_info2[0])
@@ -710,7 +709,7 @@ def _versions(sep=_SPACE_):
 
 
 __all__ = tuple(map(typename, (machine, print_, printf, typename)))
-__version__ = '25.09.05'
+__version__ = '25.10.06'
 
 if __name__ == _DMAIN_:
 

pygeodesy/karney.py

@@ -145,8 +145,8 @@ from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import _copysign, _isin, isint, neg, unsigned0, \
                              _xgeographiclib, _zip
-from pygeodesy.constants import NAN, _isfinite as _math_isfinite, _0_0, \
-                               _1_16th, _1_0, _2_0, _180_0, _N_180_0, _360_0
+from pygeodesy.constants import NAN, _isfinite as _math_isfinite, \
+                               _0_0, _1_0, _2_0, _180_0, _N_180_0, _360_0
 from pygeodesy.errors import GeodesicError, _ValueError, _xkwds
 from pygeodesy.fmath import cbrt, fremainder, norm2  # Fhorner, Fsum
 from pygeodesy.internals import _getenv, _popen2, _PYGEODESY_ENV, typename, \
@@ -165,8 +165,9 @@ from pygeodesy.utily import atan2d, sincos2d, tand, _unrollon,  fabs
 # from math import fabs  # from .utily
 
 __all__ = _ALL_LAZY.karney
-__version__ = '25.08.31'
+__version__ = '25.09.13'
 
+_1_16th     = _1_0 / 16
 _2_4_       = '2.4'
 _K_2_0      = _getenv(_PYGEODESY_ENV(typename(_xgeographiclib)[2:]), _2_)
 _K_2_4      = _K_2_0 ==  _2_4_
@@ -723,11 +724,10 @@ def _around(x):  # in .utily.sincos2d
     try:
         return _wrapped.Math.AngRound(x)
     except AttributeError:
-        b, a = _1_16th, fabs(x)
-        if a < b:
-            a -= b
-            a += b
-            x = _copysign(a, x)
+        z = _1_16th
+        w =  z - fabs(x)
+        if w > 0:  # don't "simplify" z - (z - x) to x
+            x = _copysign(z - w, x)
         return x
 
 
@@ -772,13 +772,13 @@ def _diff182(deg0, deg, K_2_0=False):
     try:
         return _wrapped.Math.AngDiff(deg0, deg)
     except AttributeError:
-        if K_2_0 or _K_2_0:  # geographiclib 2.0
-            _r, _360 = fremainder, _360_0
-            d, t = _sum2(_r(-deg0, _360),
-                         _r( deg,  _360))
-            d, t = _sum2(_r( d,    _360), t)
+        if K_2_0 or _K_2_0:  # geographiclib 2.0+
+            _r = fremainder
+            d, t = _sum2(_r(-deg0, _360_0),
+                         _r( deg,  _360_0))
+            d, t = _sum2(_r( d,    _360_0), t)
             if _isin(d, _0_0, _180_0, _N_180_0):
-                d = _copysign(d, -t if t else (deg - deg0))
+                d = _copysign(d, (-t) if t else (deg - deg0))
         else:
             _n = _norm180
             d, t = _sum2(_n(-deg0), _n(deg))

pygeodesy/named.py

@@ -1247,7 +1247,7 @@ def modulename(clas, prefixed=None):  # in .basics._xversion
     '''Return the class name optionally prefixed with the
        module name.
 
-       @arg clas: The class (any C{class}).
+       @arg clas: The class (any C{class} or C{str}).
        @kwarg prefixed: Include the module name (C{bool}), see
                         function C{classnaming}.
 

pygeodesy/streprs.py

@@ -22,7 +22,7 @@ from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS
 from math import fabs, log10 as _log10
 
 __all__ = _ALL_LAZY.streprs
-__version__ = '25.04.14'
+__version__ = '25.10.10'
 
 _at_        = 'at'         # PYCHOK used!
 _EN_PREC    =  6           # max MGRS/OSGR precision, 1 micrometer
@@ -545,22 +545,22 @@ def unstr(where, *args, **kwds_):
        @arg where: Class, function, method (C{type}) or name (C{str}).
        @arg args: Optional positional arguments.
        @kwarg kwds_: Optional keyword arguments, except C{B{_Cdot}=None},
-                     C{B{_ELLIPSIS}=False} and C{B{_fmt}=Fmt.g}.
+                     C{B{_ELLIPSIS}=False}, C{B{_fmt}=Fmt.g} and C{B{_prec}=6}.
 
        @return: Representation (C{str}).
     '''
-    def _C_e_g_kwds3(_Cdot=None, _ELLIPSIS=0, _fmt=Fmt.g, **kwds):
-        return _Cdot, _ELLIPSIS, _fmt, kwds
+    def _C_e_fmt_prec_kwds(_Cdot=None, _ELLIPSIS=0, _fmt=Fmt.g, _prec=6, **kwds):
+        return _Cdot, _ELLIPSIS, dict(fmt=_fmt, prec=_prec), kwds
 
-    C, e, g, kwds = _C_e_g_kwds3(**kwds_)
+    C, e, fmt_prec, kwds = _C_e_fmt_prec_kwds(**kwds_)
     if e and len(args) > (e + 1):
-        t  =  reprs(args[:e], fmt=g)
+        t  =  reprs(args[:e], **fmt_prec)
         t += _ELLIPSIS_,
-        t +=  reprs(args[-1:], fmt=g)
+        t +=  reprs(args[-1:], **fmt_prec)
     else:
-        t  =  reprs(args, fmt=g) if args else ()
+        t  =  reprs(args, **fmt_prec) if args else ()
     if kwds:
-        t += pairs(itemsorted(kwds), fmt=g)
+        t += pairs(itemsorted(kwds), **fmt_prec)
     n = where if isstr(where) else typename(where)  # _NN_
     if C and hasattr(C, n):
         try:  # bound method of class C?

pygeodesy/trf.py

@@ -71,7 +71,7 @@ en/how-to-deal-with-etrs89-datum-and-time-dependent-transformation-parameters-45
 
 from pygeodesy.basics import _isin, map1, neg, isidentifier, isstr, _xinstanceof, \
                              _xscalar,  typename
-from pygeodesy.constants import _float as _F, _0_0s, _0_0, _0_001, _0_5, _1_0
+from pygeodesy.constants import _float as _F, _0_0s, _0_0, _0_001, _0_5, _0_75, _1_0
 from pygeodesy.datums import Datums, _earth_datum, _equall, _GDA2020_, _Names7, \
                             _negastr, Transform, _WGS84,  _EWGS84, _operator
 # from pygeodesy.ellipsoids import _EWGS84  # from .datums
@@ -93,7 +93,7 @@ from math import ceil as _ceil, fabs
 # import operator as _operator  # from .datums
 
 __all__ = _ALL_LAZY.trf
-__version__ = '25.05.12'
+__version__ = '25.09.11'
 
 _EP0CH    =  Epoch(0, low=0)
 _Es       = {_EP0CH: _EP0CH}  # L{Epoch}s, deleted below
@@ -1336,7 +1336,7 @@ _trfX(_ITRF2020_, _ITRF94_,   epoch=_E(2015),
                               xform=_X(   6.5,     -3.9,   -77.9,    3.98,     _0_0,      _0_0,      0.36),
                               rates=_R(   0.1,     -0.6,    -3.1,    0.12,     _0_0,      _0_0,      0.02))
 _trfX(_ITRF2020_, _ITRF93_,   epoch=_E(2015),
-                              xform=_X( -65.8,      1.9,   -71.3,    4.47,     -3.36,     -4.33,     0.75),
+                              xform=_X( -65.8,      1.9,   -71.3,    4.47,     -3.36,     -4.33,    _0_75),
                               rates=_R(  -2.8,     -0.2,    -2.3,    0.12,     -0.11,     -0.19,     0.07))
 _trfX(_ITRF2020_, _ITRF92_,   epoch=_E(2015),
                               xform=_X(  14.5,     -1.9,   -85.9,    3.27,     _0_0,      _0_0,      0.36),
@@ -1459,7 +1459,7 @@ _trfX(_ITRF2000_, _ITRF93_,   epoch=_E(1988),
                               xform=_X(  12.7,      6.5,   -20.9,    1.95,     -0.39,     0.8,      -1.14),
                               rates=_R(  -2.9,     -0.2,    -0.6,    0.01,     -0.11,    -0.19,      0.07))
 _trfX(_ITRF2000_, _ITRF92_,   epoch=_E(1988),
-                              xform=_X(   1.47,     1.35,   -1.39,   0.75,     _0_0,     _0_0,      -0.18),
+                              xform=_X(   1.47,     1.35,   -1.39,  _0_75,     _0_0,     _0_0,      -0.18),
                               rates=_R(  _0_0,     -0.06,   -0.14,   0.01,     _0_0,     _0_0,       0.02))
 _trfX(_ITRF2000_, _ITRF91_,   epoch=_E(1988),
                               xform=_X(   26.7,    27.5,   -19.9,    2.15,     _0_0,     _0_0,      -0.18),

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

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: pygeodesy
-Version: 25.9.9
+Version: 25.10.10
 Summary: Pure Python geodesy tools
 Home-page: https://GitHub.com/mrJean1/PyGeodesy
 Author: Jean M. Brouwers
@@ -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: Oct 10, 2025.*
 
 License
 =======

{pygeodesy-25.9.9.dist-info → pygeodesy-25.10.10.dist-info}/RECORD

@@ -1,13 +1,13 @@
 pygeodesy/LICENSE,sha256=YfgAiyxOwY6P9Kkb1_5XN81nueTLrpb3Ffkv3EuPgFU,1144
-pygeodesy/__init__.py,sha256=IWSswTDqZR-HX68gyWXtq3a9e7f6DCQks4QpwJzZKTc,42337
+pygeodesy/__init__.py,sha256=-QifNyOMaiddBdni6VVqzLQtjjev0PS_AdfsyVXqz2k,42328
 pygeodesy/__main__.py,sha256=HpOBmuhI7yVP1FIAEugvd9Lb54kzew-nM78be0r4SC4,5597
 pygeodesy/albers.py,sha256=mb3YbVvoBq8a7AytT0HeVxe8DGEvx1KFN2Usl2ksKwk,30908
-pygeodesy/azimuthal.py,sha256=acwIA1B7jmg1oDpJBuFitnOVHUWRUk7Y73yI8zqIGrc,49854
+pygeodesy/azimuthal.py,sha256=GI7nM7RWDwOh5DKYrrZBZDHcP9aSJgCBSdgDm2ovLvY,49849
 pygeodesy/basics.py,sha256=sKIyFcBNz6UzzY2YRr0bt1WuQLnEFUMfLKW7z4vCOII,32915
 pygeodesy/booleans.py,sha256=4ROe0VvrGNW64HUrz8nKG0APs3XA_M3tV8IO631k2M8,73683
 pygeodesy/cartesianBase.py,sha256=KX74yZE-fM7wHDnrV5OkAEyKwyZJy5o9VuzVjb1caZQ,45355
 pygeodesy/clipy.py,sha256=a0azaK7YDmGDU8mVVKe2MvUiPEhqZ9XJ-9u0T1iHKbo,27779
-pygeodesy/constants.py,sha256=iOuKYv9_4-puEiUX6EuuKmBxhwXZLl2s_bLCefiAgAk,20652
+pygeodesy/constants.py,sha256=thlUllMuL_3ZY2hI3xPen1IdW0tAfQnQztMHsv01RBg,20525
 pygeodesy/css.py,sha256=vdQax5m7BMSpr7PUyqedCSvx73ilWuhXQTp39i8FOfA,25793
 pygeodesy/datums.py,sha256=Vp1fP7sSV1EuRbuO2scnZ3qiTZD1CvL4O12KSRlf1SU,34068
 pygeodesy/dms.py,sha256=egtF9gjReovgaSKJ_27KHZjvo4vLDdLEZ7yXSf19EZs,42202
@@ -22,11 +22,11 @@ pygeodesy/ellipsoidalKarney.py,sha256=W_Liw5MxTqJJGsxMi50ui8grRSPKwkG2HH-Mv_mXQY
 pygeodesy/ellipsoidalNvector.py,sha256=HLFjX6o26q2fzVCjLSL0PdEh96xhGtjfzo2V_k54j4A,30146
 pygeodesy/ellipsoidalVincenty.py,sha256=zeTvid09NGuYJ1WJ__a3ApCBGSzX0A-Rhs8-sQXMNQQ,26273
 pygeodesy/ellipsoids.py,sha256=hdBT-Vl5Ukq51s19brzctT0vZfh8qZYlbUjmrrQPE-Q,108817
-pygeodesy/elliptic.py,sha256=R_-SY6W1kzRg-dc4xIM0jBOOUwUCTx46r4ovrPQxYIs,45518
+pygeodesy/elliptic.py,sha256=uWxzmr7woRKqhsBWm6iGTwag5Ga0L184szkPWN3diSg,45515
 pygeodesy/epsg.py,sha256=VIBidzAEuoBF4LdlJ4r3WtpL0Ymhru8kf4fAgfFr5ok,8220
 pygeodesy/errors.py,sha256=qML13IDmvSaVrWp5h9DeHw1L-9J7DTwn86wLsP_sOgw,32498
-pygeodesy/etm.py,sha256=_buzd1LyldAPa-RTFlohwkkFChTO10Yyk5IrP1kZ6kE,46797
-pygeodesy/fmath.py,sha256=IbXQiGoj5MH-fbdBC1BJQMx_1WVFr05Vzpt78LrhbT4,37788
+pygeodesy/etm.py,sha256=swF0QL52uo_DK073hymix8IitMyArOf-n0LL98BGwRU,46796
+pygeodesy/fmath.py,sha256=5gPHRqD5DbPqrM8BFRFVo1IYumKUVWoTcvEtzc8BFiE,37812
 pygeodesy/formy.py,sha256=7-CdPiagdVrbDMjWieX_guN4RvFVCYKOM3yDKSFPRWU,69777
 pygeodesy/frechet.py,sha256=dBnQ5ePuaUdKv62-e9RfovFDxELsZPQEKMLww9NOf7Q,35588
 pygeodesy/fstats.py,sha256=M6UOrwmeGp6kjYuYfc6Iq9ibMJAFNm6D4C5GnvlJeuM,28348
@@ -35,14 +35,14 @@ pygeodesy/gars.py,sha256=AXioT4Lh1cwrd7fF6dHcVJw-9vu2oAnZrzo_ka2Y3Ts,11761
 pygeodesy/geodesici.py,sha256=V2FW_xQUwm6HfiVoZoBXhipsAV0EjXHnAc_wUG78kAQ,74889
 pygeodesy/geodesicw.py,sha256=xDxxxr4HOT-8Mh85zD6QZzyO-mfqUavfwQMStYrsF8c,30029
 pygeodesy/geodsolve.py,sha256=_s6Wk13yTI8eTsaMpnzjqacNXLjPh_HoevmjECtrRjE,28668
-pygeodesy/geohash.py,sha256=1P9-k9sMN-ZVFr2iWbQpFi1lg3uNmChHFCIfETkgFJU,40153
-pygeodesy/geoids.py,sha256=OQvOVutr8R1CalStgxmNu2Mk3oMl4zpCY_OKh0KHpI8,86125
+pygeodesy/geohash.py,sha256=bqJMuA0-TH2vgOd7VaW1q4dw8KL43oCmV3eSjP90Bak,40124
+pygeodesy/geoids.py,sha256=m0EqTAKCuOVxrRpoPdmgymu0i0ufV2xpdt_laJLFqqE,86129
 pygeodesy/hausdorff.py,sha256=KGB0CIeNkBqR0sXIbRJvcgXPrM-4Mpi2aBOqDeZjqZ0,32269
-pygeodesy/heights.py,sha256=_Pn-6w4gGzdXCoA6m72rJm4hb-Fe2vxc8RiRRjq5iTI,41002
-pygeodesy/internals.py,sha256=HWSgcAuhoLHpTpOv0efGKd0Gb0rZqQbeapsr_DXWggk,24717
+pygeodesy/heights.py,sha256=q7ej8RO-iZ5hfNi3OxwO3w_Opk4xHAZIzHbcbFTO3CA,41011
+pygeodesy/internals.py,sha256=8fddKdxThLLo17u3koRo4jW8g1yao1igEJMEnm-NamU,24711
 pygeodesy/interns.py,sha256=QTP_6rStQ283DE-rL7rCAcfEffuiAXustBGkCY_DHJ0,23508
 pygeodesy/iters.py,sha256=nKVp_LUQyTqOz_OLzY25g6-_xlKG7lmAnjKKMXESmiM,20345
-pygeodesy/karney.py,sha256=AlEwFlnj-gLwfGjoHGyyDj111tZE5_7pz19avsuGVK0,41626
+pygeodesy/karney.py,sha256=GJuruLHO7IjhwXYNh507EbHO5VFNXwhRpzIkaGCec7c,41652
 pygeodesy/ktm.py,sha256=9eZ5Wk2HHhRl0IyypJ4guhhxGJLriWZm0Wfb3Z9Gf7Y,27270
 pygeodesy/latlonBase.py,sha256=eF8tr0fP9nc8JMGbJ0JyeFfECZyjlBJbtIchKRdU4B4,75317
 pygeodesy/lazily.py,sha256=bLCBDH13GY1LqVzoO6xmaeGsRDrBLHr-EqxUiSvOESo,47678
@@ -50,7 +50,7 @@ pygeodesy/lcc.py,sha256=nU_aWXgsxwB53c3CKM-dBQoDRoypkrY-tVmKmU9yQQ8,25785
 pygeodesy/ltp.py,sha256=N3TMipTWKUOeZYzn3UeJJqe5J4Yk788390ohQM1ykAE,50893
 pygeodesy/ltpTuples.py,sha256=t-n7p9A3JCVasv8LKCV4PkOTmXJZi91BkBFm8pBZM5o,59101
 pygeodesy/mgrs.py,sha256=ZgCBH_M_YFxuFu5ukVOVetk-uftv-4GJP-QJrtoqtcU,29842
-pygeodesy/named.py,sha256=Ivv_O2N-haps9GwrfMQuhM-zenTkceyNGw5hT5b8coU,52777
+pygeodesy/named.py,sha256=_dlImd3OzHS1vqm673e7CEJVXkbczaAwckw8H8SGoEM,52787
 pygeodesy/namedTuples.py,sha256=JlO8uWmR65RkT7Rdgc9knfI0-bqJeA-tm2ZKsesztVU,30598
 pygeodesy/nvectorBase.py,sha256=r-O-ErwTXqDfLg4zYqhB4IqTJc8P8rbJlMm0XfQN2Gs,29927
 pygeodesy/osgr.py,sha256=4mFUcn1vMOju3PK1rKohPtIlJZrBGPJrcRjFl-CoxqE,30852
@@ -62,8 +62,8 @@ pygeodesy/solveBase.py,sha256=mZdud0LveosWEXbufVzTMM9y6wuObS0O2clDX2LmRgE,18864
 pygeodesy/sphericalBase.py,sha256=hq6dyeeIHBy_T0Q6PIwSvmo1U70Ir5KBeMTwQplzqHQ,32071
 pygeodesy/sphericalNvector.py,sha256=bu0yB4PEq1cs6BS8QFqtljWpVEw7ESwq2yILOzGDuAc,54948
 pygeodesy/sphericalTrigonometry.py,sha256=PqxNMSK8DI7TJuOtkQ5afzJtLfDzOkVe_tRdLTq9x14,64398
-pygeodesy/streprs.py,sha256=-4LErmfPLtGt1tLDBHhJam-zYq8p2Z2C5wDCMELEml4,23599
-pygeodesy/trf.py,sha256=gpsOJVER1KhITOtyiENQ4EnRONinLkzbujsfzzxHVv8,119209
+pygeodesy/streprs.py,sha256=UYlk5v-GEJHMDiHZsqLssY3h6J5LbFw9dsCsKJbiERs,23684
+pygeodesy/trf.py,sha256=M0kzyclDJ5q4wCr-2HDIj6_J42FsrffnFGiyPkgHEmE,119216
 pygeodesy/triaxials.py,sha256=V0MSYJxRgYDhZcF7HXnu7oFzRyWJkhwc6scuNjc6VVg,64046
 pygeodesy/units.py,sha256=nH6uNqtaILtGRXqCcHVYBTfhGrv-WFt3w1aso6PA0kM,35419
 pygeodesy/unitsBase.py,sha256=cVNKG24Fov73gLFPGwGdjvXmM-8fd3MDufVubABAMeE,14141
@@ -104,8 +104,8 @@ pygeodesy/geodesicx/_C4_27.py,sha256=grqtju94Nu9V_eQPqe_7aRbLCYflxRuK127pPOPwuDA
 pygeodesy/geodesicx/_C4_30.py,sha256=I4ZuTy64yOq1IVgRtc_kqOL0HnKzXeTea7kTFSolcJs,201921
 pygeodesy/geodesicx/__init__.py,sha256=1wUYni2aucT66bNkOqp1ZOix6a4rlBRHJimUDUWIgt4,2468
 pygeodesy/geodesicx/__main__.py,sha256=ZWx2vPPaQEXwFBQWOjKYrjM9Deb0b-_778PxwiM1Rww,4330
-pygeodesy/geodesicx/gx.py,sha256=Agny4JdMTBRkOPskWjpNzc3609ezJWRWVvLdFCm3uPE,60352
-pygeodesy/geodesicx/gxarea.py,sha256=UgniZh-WGteC86mMT7tWO_AG05gYJ5n4136cCk8hitc,19571
+pygeodesy/geodesicx/gx.py,sha256=uEXulVciKipCz28VMQqhtV0kLuy71vy-EFF42GWWi-Q,60438
+pygeodesy/geodesicx/gxarea.py,sha256=a_yhNsjjRM0E9-WYBbhHjaefjLuAzhNQxTHiv0LnvO0,19588
 pygeodesy/geodesicx/gxbases.py,sha256=qb1tYsh2nIvyXoVhU-pdmolANHj9tHZxiH2f7Ls3Jkw,6319
 pygeodesy/geodesicx/gxline.py,sha256=9aUbLEok0ewg8j1UpZ1Ud7X0zH8Mu6QoUO2bjdpmEag,28284
 pygeodesy/rhumb/__init__.py,sha256=MByyXXwg0kJRHckEfaxWMW_2i4U15rRRgN7iTzeOK_I,2207
@@ -113,7 +113,7 @@ pygeodesy/rhumb/aux_.py,sha256=2FbR_0Y0ofOROay2Z16oiLILhd7G50LUiatkQhBwEMA,16117
 pygeodesy/rhumb/bases.py,sha256=R-QLV8a0KZUvaI1wEkyAG4yO-IQNNK4I2O7WRIMBJuo,54218
 pygeodesy/rhumb/ekx.py,sha256=LlMjVJqb-Q5rJnexJi-0TPRkUZkJaPQ7fVRypxkxNjo,24058
 pygeodesy/rhumb/solve.py,sha256=z8z_XYObTgz7w1skNLNcLBpe-EO_r0H4sVcZGlBcWnc,24005
-pygeodesy-25.9.9.dist-info/METADATA,sha256=zUFYEcoqpmFzxV3TQ7Pe_bXu__9khC3lO21AD_FVny4,20186
-pygeodesy-25.9.9.dist-info/WHEEL,sha256=Kh9pAotZVRFj97E15yTA4iADqXdQfIVTHcNaZTjxeGM,110
-pygeodesy-25.9.9.dist-info/top_level.txt,sha256=cEQPatCXzKZqrivpULC5V5fuy9_V_bAwaP_gUGid7pQ,10
-pygeodesy-25.9.9.dist-info/RECORD,,
+pygeodesy-25.10.10.dist-info/METADATA,sha256=F_luIoHZU0UkQQxINpuceZKW7FOMvNxy3MOgTxdyocc,20126
+pygeodesy-25.10.10.dist-info/WHEEL,sha256=Kh9pAotZVRFj97E15yTA4iADqXdQfIVTHcNaZTjxeGM,110
+pygeodesy-25.10.10.dist-info/top_level.txt,sha256=cEQPatCXzKZqrivpULC5V5fuy9_V_bAwaP_gUGid7pQ,10
+pygeodesy-25.10.10.dist-info/RECORD,,