pygeodesy 25.12.31__py2.py3-none-any.whl → 26.2.2__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.02.02'
 # 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/auxilats/auxDLat.py

@@ -26,7 +26,7 @@ from pygeodesy.utily import atan1, atan2  # from .auxilats.auxily
 from math import cos, sin, sqrt
 
 __all__ = ()
-__version__ = '25.05.12'
+__version__ = '26.01.20'
 
 
 class AuxDLat(AuxLat):
@@ -107,10 +107,9 @@ class AuxDLat(AuxLat):
             t2  = _1_0 + t**2
             Dt *= _2_0 / t2
             sk2 = (d * Dt)**2 * k2
-            d2  = _1_0 - sk2
             c2  = ((_1_0 - t) * (_1_0 + t) / t2)**2 if t else _1_0
             # E(z)/sin(z)
-            Dt *= _Ef._RFRD(c2, d2, _1_0, sk2) - k2 * sx * sy
+            Dt *= _Ef._fRF3RD(c2, _1_0, sk2) - k2 * sx * sy
         return Dt
 
     def DIsometric(self, Phi1, Phi2):

pygeodesy/auxilats/auxLat.py

@@ -48,7 +48,7 @@ except ImportError:  # Python 3.11-
         return pow(_2_0, x)
 
 __all__ = ()
-__version__ = '25.05.12'
+__version__ = '26.01.20'
 
 _TRIPS = 1024  # XXX 2 or 3?
 
@@ -621,19 +621,18 @@ class AuxLat(AuxAngle):
             ka,  kb  = kb,  ka
             ka1, kb1 = kb1, ka1
             sb,  cb  = cb,  sb
-        # now a, b = larger/smaller semiaxis
-        # Beta measured from larger semiaxis
+        # now a, b = larger/smaller semi-axis
+        # Beta measured from larger semi-axis
         # kb, ka = modulus-squared for distance from Beta = 0, pi/2
         # NB kb <= 0; 0 <= ka <= 1
         # sa = b*E(Beta, sqrt(kb))
         # sb = a*E(Beta',sqrt(ka))
         # 1 - ka * (1 - sb2) = 1 - ka + ka*sb2
-        sb2 =  sb**2
-        cb2 =  cb**2
-        da2 =  ka1 + ka * sb2
-        db2 = _1_0 - kb * sb2
+        sb2 = sb**2
+        cb2 = cb**2
+        da2 = ka1 + ka * sb2
         # DLMF Eq. 19.25.9
-        my = b * sb * _Ef._RFRD(cb2, db2, _1_0, kb * sb2)
+        my = b * sb * _Ef._fRF3RD(cb2, _1_0, kb * sb2)
         # DLMF Eq. 19.25.10 with complementary angles
         mx = a * cb * (_Ef.fRF(sb2,  da2, _1_0) * ka1 +
             ka * cb2 * _Ef.fRD(sb2, _1_0,  da2, _3_0) * ka1 +

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