pygeodesy 25.12.31__py2.py3-none-any.whl → 26.1.16__py2.py3-none-any.whl

pygeodesy/__init__.py

@@ -132,17 +132,17 @@ Tests
 =====
 
 The tests ran with Python 3.14.2 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.1),
-Python 3.13.9 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.1),
+Python 3.13.11 (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.7 and
 U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.7),
-Python 3.12.7 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0,
+Python 3.12.10 (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.7,
 U{GeodS3olve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.7,
 U{IntersectTool<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.7 and
 U{RhumbSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.7),
-Python 3.11.5 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0,
+Python 3.11.9 (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),
 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,
@@ -158,14 +158,14 @@ env variable C{PYGEODESY_WARNINGS=on} for all Python versions.  The results of t
 the distribution files.
 
 Test coverage has been measured with U{coverage<https://PyPI.org/project/coverage>} 7.10.7 using Python
-3.14.2, 3.13.9 and 3.12.7.  The complete coverage report in HTML and a PDF summary are included in the
+3.14.2, 3.13.11 and 3.12.10.  The complete coverage report in HTML and a PDF summary are included in the
 distribution files.
 
-Python 3.14.2, 3.13.9, 3.12.7 and 3.11.5 run on Apple Si M4 (C{arm64}), I{natively}.  Python 2.7.18 runs
+Python 3.14.2, 3.13.11, 3.12.10 and 3.11.9 run on Apple Si M4 (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.14.2 (and U{geographiclib<https://PyPI.org/project/geographiclib>} 2.1) on
-U{Debian 12<https://Cirrus-CI.com/github/mrJean1/PyGeodesy/master>} in 64-bit only, with Python 3.13.9 (and
+U{Debian 12<https://Cirrus-CI.com/github/mrJean1/PyGeodesy/master>} in 64-bit only, with Python 3.13.11 (and
 U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0) on U{Windows 2019Server
 <https://CI.AppVeyor.com/project/mrJean1/pygeodesy>} in 64-bit only and with Python 2.7.18 (and U{geographiclib
 <https://PyPI.org/project/geographiclib>} 1.52) on U{Windows 10<https://CI.AppVeyor.com/project/mrJean1/pygeodesy>}
@@ -174,7 +174,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-7, 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-9, 3.12.0-7, 3.11.2-5, 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,
@@ -195,7 +195,7 @@ 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.9 and with
+546532_PyChecker_postprocessor>} with U{Ruff<https://GitHub.com/astral-sh/ruff>} using Python 3.13.11 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 26.2 Tahoe.
@@ -626,7 +626,7 @@ else:
 
 from pygeodesy.internals import _version2,  _DOT_  # noqa: E402
 # from pygeodesy.interns import _DOT_  # from .internals
-__version__ = '25.12.31'
+__version__ = '26.01.16'
 # see setup.py for similar logic
 version     = _DOT_(*_version2(__version__, n=3))
 

pygeodesy/albers.py

@@ -38,11 +38,11 @@ from pygeodesy.utily import atan1, atan1d, atan2, degrees360, sincos2, \
 from math import atanh, degrees, fabs, radians, sqrt
 
 __all__ = _ALL_LAZY.albers
-__version__ = '25.05.12'
+__version__ = '26.01.06'
 
 _k1_    = 'k1'
-_NUMIT  =   9  # XXX 4?
-_NUMIT0 =  41  # XXX 21?
+_MAXIT  =   9  # XXX 4?
+_MAXIT0 =  41  # XXX 21?
 _TERMS  =  31  # XXX 16?
 _TOL0   =  sqrt3(_TOL)
 
@@ -425,7 +425,7 @@ class _AlbersBase(_NamedBase):
         _Ta02  =  Fsum(ta0).fsum2f_
         _1, _2 = _1_0, _2_0
         _4, _6 = _4_0, _6_0
-        for self._iteration in range(1, _NUMIT0):  # 4 trips
+        for self._iteration in range(1, _MAXIT0):  # 4 trips
             ta02  =  ta0**2
             sca02 =  ta02 + _1
             sca0  =  sqrt(sca02)
@@ -492,7 +492,7 @@ class _AlbersBase(_NamedBase):
         ta   =  txi
         tol  = _tol(_TOL, ta)
         _Ta2 =  Fsum(ta).fsum2f_
-        for self._iteration in range(1, _NUMIT):  # max 2, mean 1.99
+        for self._iteration in range(1, _MAXIT):  # max 2, mean 1.99
             # dtxi / dta = (scxi / sca)^3 * 2 * (1 - e^2)
             #            / (qZ * (1 - e^2 * sa^2)^2)
             ta2   =  ta**2

pygeodesy/cartesianBase.py

@@ -43,7 +43,7 @@ from pygeodesy.vector3d import Vector3d, _xyzhdlln4
 # from math import degrees, fabs, radians, sqrt  # from .fmath, .utily
 
 __all__ = _ALL_LAZY.cartesianBase
-__version__ = '25.11.07'
+__version__ = '26.01.06'
 
 _r_         = 'r'
 _resections = _MODS.into(resections=__name__)
@@ -361,7 +361,7 @@ class CartesianBase(Vector3d, _EcefLocal):
         d = self.datum if earth is None else earth
         if normal and d is self.datum:
             r = self._height4
-        elif isinstance(d, _MODS.triaxials.triaxial5.Triaxial_):
+        elif isinstance(d, _MODS.triaxials.Triaxial_):
             r = d.height4(self, normal=normal)
             try:
                 d = d.toEllipsoid(name=n)

pygeodesy/constants.py

@@ -26,7 +26,7 @@ except ImportError:  # Python 2-
     _inf, _nan = float(_INF_), float(_NAN_)
 
 __all__ = _ALL_LAZY.constants
-__version__ = '25.12.06'
+__version__ = '26.01.14'
 
 
 def _copysign_0_0(y):
@@ -295,6 +295,11 @@ OVERFLOW = _Float(OVERFLOW=_1_0 / EPS0)   # PYCHOK = 2.028240960365e+31
 # _1SQRT2= _Float(_1SQRT2 =sqrt(_2_0) + 1)
 # 0.707106781186547_524_400_844_362_104_849_039_284_835_937_688_474_036_588_339_868_99
 _SQRT2_2 = _Float(_SQRT2_2=sqrt(_0_5))  # PYCHOK = 0.707106781186547_5 == sqrt(2) / 2
+# sqrt(3) <https://WikiPedia.org/wiki/Square_root_of_3>
+# 1.732050807568877_293_527_446_341_505_872_366_942_805_253_810_380_628_055_806
+_SQRT3   = _Float(_SQRT3  =sqrt(_3_0))  # PYCHOK = 1.732050807568877_2 == sqrt(3)
+# 0.866025403784438_646_763_723_170_752_936_183_471_402_626_905_190_314_027_903
+_SQRT3_2 = _Float(_SQRT3_2=sqrt(_0_75))  # PYCHOK = 0.866025403784438_6 == sqrt(3) / 2
 
 INF   =  Float(INF =_inf)    # PYCHOK INFinity, see function L{isinf}, L{isfinite}, NOT _Float!
 INT0  =  Int(  INT0= 0)      # PYCHOK unique int(0) instance, see .fsums, useZ=False

pygeodesy/datums.py

@@ -94,7 +94,7 @@ from pygeodesy.units import _isRadius, Radius_,  radians
 # import operator as _operator  # from .fmath
 
 __all__ = _ALL_LAZY.datums
-__version__ = '25.05.12'
+__version__ = '26.01.13'
 
 _a_ellipsoid_ = _UNDER_(_a_, _ellipsoid_)
 _BD72_        = 'BD72'
@@ -723,13 +723,10 @@ assert _WGS84.ellipsoid is _EWGS84
 
 if __name__ == _DMAIN_:
 
-    from pygeodesy.interns import _COMMA_, _NL_, _NLATvar_
-    from pygeodesy import printf
-
-    # __doc__ of this file, force all into registery
-    for r in (Datums, Transforms):
-        t = [NN] + r.toRepr(all=True, asorted=True).split(_NL_)
-        printf(_NLATvar_.join(i.strip(_COMMA_) for i in t))
+    from pygeodesy.internals import _pregistry
+    # __doc__ of this file, force all into registry
+    _pregistry(Datums)
+    _pregistry(Transforms)
 
 # **) MIT License
 #

pygeodesy/ellipsoids.py

@@ -96,7 +96,7 @@ from pygeodesy.utily import atan1, atan1d, atan2b, degrees90, m2radians, radians
 from math import asinh, atan, atanh, cos, degrees, exp, fabs, radians, sin, sinh, sqrt, tan  # as _tan
 
 __all__ = _ALL_LAZY.ellipsoids
-__version__ = '25.12.14'
+__version__ = '26.01.13'
 
 _f_0_0    = Float(f =_0_0)  # zero flattening
 _f__0_0   = Float(f_=_0_0)  # zero inverse flattening
@@ -1140,7 +1140,8 @@ class Ellipsoid(_NamedEnumItem):
                  methods L{Ellipsoid.height4} and L{Triaxial.hartzell4}.
         '''
         try:
-            v, d, i = _MODS.triaxials.triaxial5._hartzell3(pov, los, self._triaxial)
+            m = _MODS._triaxials_triaxial5
+            v, d, i = m._hartzell3(pov, los, self._triaxial)
         except Exception as x:
             raise IntersectionError(pov=pov, los=los, cause=x)
         return Vector4Tuple(v.x, v.y, v.z, d, iteration=i, name__=self.hartzell4)
@@ -1194,7 +1195,8 @@ class Ellipsoid(_NamedEnumItem):
                 v = v.times_(t, t, 0)  # force z=0.0
                 h = x - a  # equatorial
             else:  # normal in 1st quadrant
-                x, y, i = _MODS.triaxials.triaxial5._plumbTo3(x, y, self)
+                m = _MODS._triaxials_triaxial5
+                x, y, i = m._plumbTo3(x, y, self)
                 t, v = v, v.times_(x, x, y)
                 h = t.minus(v).length
 
@@ -1910,8 +1912,8 @@ class Ellipsoid(_NamedEnumItem):
     def _triaxial(self):
         '''(INTERNAL) Get this ellipsoid's un-/ordered C{Triaxial/_}.
         '''
-        a, b, m = self.a, self.b, _MODS.triaxials
-        T = m.Triaxial if a > b else m.Triaxial_
+        a, b, t = self.a, self.b, _MODS.triaxials
+        T = t.Triaxial if a > b else t.Triaxial_
         return T(a, a, b, name=self.name)
 
     @Property_RO
@@ -2433,7 +2435,6 @@ _EWGS84 = Ellipsoids.WGS84  # (INTERNAL) shared
 
 if __name__ == _DMAIN_:
 
-    from pygeodesy.interns import _COMMA_, _NL_, _NLATvar_
     from pygeodesy import nameof, printf
 
     for E in (_EWGS84, Ellipsoids.GRS80,  # NAD83,
@@ -2450,9 +2451,9 @@ if __name__ == _DMAIN_:
         printf('# %s  %s', Ellipsoid.BetaKs.name, fstr(E.BetaKs,  prec=20))
         printf('# %s %s', nameof(Ellipsoid.KsOrder), E.KsOrder)  # property
 
+    from pygeodesy.internals import _pregistry
     # __doc__ of this file, force all into registry
-    t = [NN] + Ellipsoids.toRepr(all=True, asorted=True).split(_NL_)
-    printf(_NLATvar_.join(i.strip(_COMMA_) for i in t))
+    _pregistry(Ellipsoids)
 
 # % python3.13 -m pygeodesy.ellipsoids