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

pygeodesy/__init__.py

@@ -8,8 +8,8 @@ geodetic (lat-/longitude), geocentric (U{ECEF<https://WikiPedia.org/wiki/ECEF>}
 <https://GeographicLib.SourceForge.io/1.44/triaxial.html>} coordinates.
 
 Transcoded in part from U{JavaScript originals<https://GitHub.com/ChrisVeness/geodesy>} by I{Chris Veness (C)
-2005-2024} and from several U{C++ classes<https://GeographicLib.SourceForge.io/C++/doc/annotated.html>} by I{Charles
-F. F. Karney (C) 2008-2024} and published under the same U{MIT License<https://OpenSource.org/licenses/MIT>}**.
+2005-2025} and from several U{C++ classes<https://GeographicLib.SourceForge.io/C++/doc/annotated.html>} by I{Charles
+F. F. Karney (C) 2008-2025} and published under the same U{MIT License<https://OpenSource.org/licenses/MIT>}**.
 
 There are four modules for ellipsoidal earth models, C{ellipsoidalExact}, C{-Karney}, C{-Vincenty} and C{-Nvector}
 and two for spherical ones, C{sphericalTrigonometry} and C{-Nvector}.  Each module provides a geodetic B{C{LatLon}}
@@ -101,9 +101,15 @@ and L{GnomonicGeodSolve} depend on I{Karney}'s C++ utility U{GeodSolve
 <https://GeographicLib.SourceForge.io/C++/doc/GeodSolve.1.html>} to be executable and set with
 env variable C{PYGEODESY_GEODSOLVE} or with property L{Ellipsoid.geodsolve}.
 
+Triaxial geodesic classes L{Geodesic3Solve} and L{GeodesicLine3Solve} in module L{geod3solve} need
+I{Karney}'s C++ utility U{Geod3Solve<https://GeographicLib.SourceForge.io/C++/doc/Geod3Solve.1.html>}
+to be executable and set with env variable C{PYGEODESY_GEOD3SOLVE} or with property
+L{Geodesic3Solve.Geod3Solve<geod3solve._Geodesic3SolveBase.Geod3Solve>}.
+
 Class L{Intersectool} in module L{geodesici} needs I{Karney}'s C++ utility U{IntersectTool
 <https://GeographicLib.SourceForge.io/C++/doc/IntersectTool.1.html>} to be executable and set with
-env variable C{PYGEODESY_INTERSECTTOOL}.
+env variable C{PYGEODESY_INTERSECTTOOL} or with property L{Intersectool.IntersectTool
+<geodesici.Intersectool.IntersectTool>}.
 
 To compare C{MGRS} results from modules L{mgrs} and C{testMgrs} with I{Karney}'s C++ utility
 U{GeoConvert<https://GeographicLib.SourceForge.io/C++/doc/GeoConvert.1.html>}, the latter must
@@ -126,38 +132,40 @@ 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.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{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.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.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),
-Python 3.11.5 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0,
+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.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,
-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 26.1 Tahoe.
+U{GeoConvert<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.7,
+U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.7,
+U{Geod3Solve<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), all in 64-bit on
+macOS 26.2 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.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 M4 Si (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.13.9 (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.12.8 (and
+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.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>}
@@ -166,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,
@@ -187,10 +195,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.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.1 Tahoe.
+2.7.18, both in 64-bit on macOS 26.2 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>}.
@@ -210,7 +218,7 @@ The following environment variables are observed by C{pygeodesy}:
  - C{PYGEODESY_FSUM_RESIDUAL} - see module L{fsums<pygeodesy.fsums>} and method L{RESIDUAL<pygeodesy.Fsum.RESIDUAL>}.
  - C{PYGEODESY_GEOCONVERT} - see module L{mgrs<pygeodesy.mgrs>}.
  - C{PYGEODESY_GEODSOLVE} - see module L{geodsolve<pygeodesy.geodsolve>}.
- - C{PYGEODESY_GEOD3SOLVE} - see module L{geodsolve<pygeodesy.geod3solve>}.
+ - C{PYGEODESY_GEOD3SOLVE} - see module L{geod3solve<pygeodesy.geod3solve>}.
  - C{PYGEODESY_INTERSECTTOOL} - see module L{geodesici<pygeodesy.geodesici>}.
  - C{PYGEODESY_LAZY_IMPORT} - see module L{lazily<pygeodesy.lazily>} and variable L{isLazy<pygeodesy.isLazy>}.
  - C{PYGEODESY_NOTIMPLEMENTED} - C{__special__} methods return C{NotImplemented} if set to "std".
@@ -618,7 +626,7 @@ else:
 
 from pygeodesy.internals import _version2,  _DOT_  # noqa: E402
 # from pygeodesy.interns import _DOT_  # from .internals
-__version__ = '25.12.12'
+__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
@@ -332,6 +337,8 @@ R_VM  = _Radius(R_VM=6366707.0194937)  # PYCHOK aViation/naVigation earth radius
 
 _INF_NAN_NINF = {INF, NAN, NINF, _inf, _nan}
 _pos_self     = _1_0.__pos__() is _1_0  # PYCHOK in .fsums, .vector3dBase
+_1_3rd        =  1 / 3  # in .fmath
+_2_3rd        =  2 / 3  # in .fmath, .formy
 
 
 def _0_0s(n):

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

@@ -67,10 +67,10 @@ from __future__ import division as _; del _  # noqa: E702 ;
 # from pygeodesy.albers import AlbersEqualAreaCylindrical  # _MODS
 from pygeodesy.basics import copysign0, isbool, _isin, isint,  typename
 from pygeodesy.constants import EPS, EPS_2, EPS0, EPS02, EPS1, INF, NINF, \
-                               _over, PI_2, PI_3, PI4, R_M, R_MA, R_FM, _EPSqrt, \
+                               _EPSqrt, PI_2, PI_3, PI2, PI4, R_M, R_MA, R_FM, \
                                _EPStol as _TOL, _floatuple as _T, _isfinite, \
-                               _0_0s, _0_0, _0_5, _1_0, _1_EPS, _2_0, _4_0, _90_0, \
-                               _0_25, _3_0  # PYCHOK used!
+                               _over, _0_0s, _0_0, _0_5, _1_0, _1_EPS, _2_0, \
+                               _4_0, _90_0, _0_25, _3_0  # PYCHOK used!
 from pygeodesy.errors import _AssertionError, IntersectionError, _ValueError, _xattr, _xkwds_not
 from pygeodesy.fmath import cbrt, cbrt2, fdot, Fhorner, fpowers, hypot, hypot_, \
                             hypot1, hypot2, sqrt3,  Fsum
@@ -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.11.26'
+__version__ = '26.01.13'
 
 _f_0_0    = Float(f =_0_0)  # zero flattening
 _f__0_0   = Float(f_=_0_0)  # zero inverse flattening
@@ -111,19 +111,19 @@ def _aux(lat, inverse, auxLat, clip=90):
     return Lat(lat, clip=clip, name=_lat_ if inverse else typename(auxLat))
 
 
-def _s2_c2(phi):
+def _sin2cos2(rad):
     '''(INTERNAL) Return 2-tuple C{(sin(B{phi})**2, cos(B{phi})**2)}.
     '''
-    if phi:
-        s2 = sin(phi)**2
+    if rad:
+        s2 = sin(rad)**2
         if s2 > EPS:
             c2 = _1_0 - s2
             if c2 > EPS:
                 if c2 < EPS1:
                     return s2, c2
             else:
-                return _1_0, _0_0  # phi == PI_2
-    return _0_0, _1_0  # phi == 0
+                return _1_0, _0_0  # rad == PI_2
+    return _0_0, _1_0  # rad == 0
 
 
 class a_f2Tuple(_NamedTuple):
@@ -871,6 +871,12 @@ class Ellipsoid(_NamedEnumItem):
 
     equatoradius = a  # Requatorial
 
+    @Property_RO
+    def equatorimeter(self):
+        '''Get the ellipsoid's I{equatorial} perimeter (C{meter}).
+        '''
+        return Meter(equatorimeter=self.a * PI2)
+
     def e2s(self, s):
         '''Compute norm M{sqrt(1 - e2 * s**2)}.
 
@@ -1134,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)
@@ -1188,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
 
@@ -1381,6 +1389,12 @@ class Ellipsoid(_NamedEnumItem):
 
     polaradius = b  # Rpolar
 
+    @property_RO
+    def polarimeter(self):
+        '''Get the ellipsoid's I{polar}, meridional perimeter (C{meter}).
+        '''
+        return Meter(polarimeter=self.L * _4_0)
+
 #   Q = A  # I{meridian arc unit} C{Q}, the mean, meridional length I{per radian}
 
     @deprecated_Property_RO
@@ -1484,7 +1498,7 @@ class Ellipsoid(_NamedEnumItem):
         r, p = self.a, Phid(lat)
         if p and self.f:
             if fabs(p) < PI_2:
-                s2, c2 = _s2_c2(p)
+                s2, c2 = _sin2cos2(p)
                 # R == sqrt((a2**2 * c2 + b2**2 * s2) / (a2 * c2 + b2 * s2))
                 #   == sqrt(a2**2 * (c2 + (b2 / a2)**2 * s2) / (a2 * (c2 + b2 / a2 * s2)))
                 #   == sqrt(a2 * (c2 + (b2 / a2)**2 * s2) / (c2 + (b2 / a2) * s2))
@@ -1709,7 +1723,7 @@ class Ellipsoid(_NamedEnumItem):
         '''(INTERNAL) Helper for C{rocAzimuth} and C{rocBearing}.
         '''
         if self.f:
-            s2, c2 = _s2_c2(radians(deg))
+            s2, c2 = _sin2cos2(radians(deg))
             m, n = self.roc2_(Phid(lat))
             if n < m:  # == n / (c2 * n / m + s2)
                 c2 *= n / m
@@ -1751,7 +1765,7 @@ class Ellipsoid(_NamedEnumItem):
         # ... requires 1 or 2 sqrt
         g = self.b
         if self.f:
-            s2, c2 = _s2_c2(Phid(lat))
+            s2, c2 = _sin2cos2(Phid(lat))
             g = _over(g, c2 + self.b2_a2 * s2)
         return Radius(rocGauss=g)
 
@@ -1898,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
@@ -2421,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,
@@ -2438,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