pygeodesy 25.1.5__py2.py3-none-any.whl → 25.4.8__py2.py3-none-any.whl

pygeodesy/__init__.py

@@ -7,7 +7,7 @@ methods for geodetic (lat-/longitude), geocentric (U{ECEF<https://WikiPedia.org/
 cartesian) and certain U{triaxial ellipsoidal<https://GeographicLib.SourceForge.io/1.44/triaxial.html>}
 coordinates.
 
-Transcoded partially from U{JavaScript originals<https://GitHub.com/ChrisVeness/geodesy>} by I{Chris Veness (C)
+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>}**.
 
@@ -64,19 +64,19 @@ wiki/Hausdorff_distance>} distances.
 Installation
 ============
 
-To install PyGeodesy, type C{python[3] -m pip install PyGeodesy} or C{python[3] -m easy_install
-PyGeodesy} in a terminal or command window.
+To install C{pygeodesy}, type C{python[3] -m pip install pygeodesy} or C{python[3] -m easy_install
+pygeodesy} in a terminal or command window.
 
-If the wheel C{PyGeodesy-yy.m.d-py2.py3-none-any.whl} is missing in U{PyPI Download files<https://
-PyPI.org/project/PyGeodesy/#files>}, download the file from U{GitHub/dist<https://GitHub.com/mrJean1/
+If the wheel C{pygeodesy-yy.m.d-py2.py3-none-any.whl} is missing in U{PyPI Download files<https://
+PyPI.org/project/pygeodesy/#files>}, download the file from U{GitHub/dist<https://GitHub.com/mrJean1/
 PyGeodesy/tree/master/dist>}.  Install that with C{python[3] -m pip install <path-to-downloaded-wheel>}
 and verify with C{python[3] -m pygeodesy}.
 
-Alternatively, download C{PyGeodesy-yy.m.d.tar.gz} from U{PyPI<https://PyPI.org/project/PyGeodesy>}
+Alternatively, download C{pygeodesy-yy.m.d.tar.gz} from U{PyPI<https://PyPI.org/project/pygeodesy>}
 or U{GitHub<https://GitHub.com/mrJean1/PyGeodesy>}, C{unzip} the downloaded file, C{cd} to
-directory C{Pygeodesy-yy.m.d} and type C{python[3] setup.py install}.
+directory C{pygeodesy-yy.m.d} and type C{python[3] setup.py install}.
 
-To run all PyGeodesy tests, type C{python[3] test/run.py} or type C{python[3] test/unitTestSuite.py}
+To run all C{pygeodesy} tests, type C{python[3] test/run.py} or type C{python[3] test/unitTestSuite.py}
 before or after installation.
 
 Dependencies
@@ -113,16 +113,16 @@ env variable C{PYGEODESY_RHUMBSOLVE} or with property L{Ellipsoid.rhumbsolve}.
 Documentation
 =============
 
-In addition to the C{pygeodesy} package, the U{PyGeodesy<https://PyPI.org/project/PyGeodesy>}
+In addition to the C{pygeodesy} package, the U{pygeodesy<https://PyPI.org/project/pygeodesy>}
 U{distribution files<https://GitHub.com/mrJean1/PyGeodesy/tree/master/dist>} contain the tests,
 the test results (on macOS only) and the complete U{documentation<https://mrJean1.GitHub.io/PyGeodesy>}
 (generated by U{Epydoc<https://PyPI.org/project/epydoc>} using command line:
-C{epydoc --html --no-private --no-source --name=PyGeodesy --url=... -v pygeodesy}).
+C{epydoc --html --no-private --no-source --name=pygeodesy --url=... -v pygeodesy}).
 
 Tests
 =====
 
-The tests ran with Python 3.13.1 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0),
+The tests ran with Python 3.13.2 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0),
 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.3,
@@ -139,20 +139,21 @@ U{numpy<https://PyPI.org/project/numpy>} 1.16.6, U{scipy<https://PyPI.org/projec
 U{GeoConvert<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.3,
 U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.3,
 U{IntersectTool<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.3 and
-U{RhumbSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.3), all in 64-bit on macOS 14.6.1 Sonoma.
+U{RhumbSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.3), all in 64-bit on
+macOS 15.4 Sequoia.
 
 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.1, 3.12.7, 3.11.5 and 3.10.8.  The complete coverage report in HTML and a PDF summary are included in
+3.13.2, 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.
 
-Python 3.13.1, 3.12.7, 3.11.5 and 3.10.8 run on Apple M1 Silicon (C{arm64}), I{natively}.  Python 2.7.18 runs
+Python 3.13.2, 3.12.7, 3.11.5 and 3.10.8 run on Apple M1 Silicon (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.1 (and U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0) on
+The tests also ran with Python 3.13.2 (and U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0) on
 U{Debian 12<https://Cirrus-CI.com/github/mrJean1/PyGeodesy/master>} in 64-bit only, with Python 3.12.5 (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
@@ -162,7 +163,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, 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-1, 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,
@@ -171,7 +172,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 14.0-5 Sonoma, 13.0-5.2 Ventura, 12.1-6 Monterey, 11.0-5.2-6.1
+<https://PyPI.org/project/numpy>} 1.11.3) on macOS 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
@@ -186,7 +187,7 @@ All Python source code has been statically U{checked<https://GitHub.com/ActiveSt
 Python/546532_PyChecker_postprocessor>} 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 and with U{Flake8<https://PyPI.org/project/flake8>}
-using Python 3.11.5, both in 64-bit on macOS 14.6.1 Sonoma.
+using Python 3.11.5, both in 64-bit on macOS 15.4 Sequoia.
 
 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>}.
@@ -197,7 +198,7 @@ characters and if so, I{not} C{unicode}.
 Env variables
 =============
 
-The following environment variables are observed by C{PyGeodesy}:
+The following environment variables are observed by C{pygeodesy}:
 
  - C{PYGEODESY_EXCEPTION_CHAINING} - see module L{errors<pygeodesy.errors>}.
  - C{PYGEODESY_FMT_FORM} - see module L{dms<pygeodesy.dms>}.
@@ -355,7 +356,7 @@ SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.}
 
 @var pygeodesy_abspath: Fully qualified C{pygeodesy} directory name (C{str}).
 
-@var version: Normalized C{PyGeodesy} version (C{str}).
+@var version: Normalized C{pygeodesy} version (C{str}).
 '''
 
 import os.path as _os_path
@@ -598,7 +599,7 @@ else:
 
 from pygeodesy.internals import _version2,  _DOT_  # PYCHOK import
 # from pygeodesy.interns import _DOT_  # from .internals
-__version__ = '25.01.05'
+__version__ = '25.04.08'
 # see setup.py for similar logic
 version     = _DOT_(*_version2(__version__, n=3))
 

pygeodesy/auxilats/_CX_Rs.py

@@ -11,8 +11,10 @@ MIT/X11 License.  For more information, see <https:#GeographicLib.SourceForge.io
 # make sure int/int division yields float quotient, see .basics
 from __future__ import division as _; del _  # PYCHOK semicolon
 
-from pygeodesy.constants import _floats
+# from pygeodesy.basics import _splituple  # _MODS
+from pygeodesy.constants import _floats as _Cs
 from pygeodesy.errors import _AssertionError,  _MODS
+# from pygeodesy.internals import _sizeof  # _MODS
 from pygeodesy.interns import NN, MISSING, _COMMA_, _duplicate_, _NL_, \
                              _QUOTE3_, _SLASH_,  _ELLIPSIS4_  # PYCHOK used!
 # from pygeodesy.lazily import _ALL_MODS as _MODS  # from .errors
@@ -20,7 +22,7 @@ from pygeodesy.named import ADict,  Property_RO
 # from pygeodesy.props import Property_RO  # from .named
 
 __all__ = ()
-__version__ = '24.09.04'
+__version__ = '25.04.08'
 
 
 class _Rcoeffs(ADict):
@@ -31,17 +33,19 @@ class _Rcoeffs(ADict):
         '''
         try:
             if not isinstance(coeffs, dict):
-                raise _RdictError(coeffs=type(coeffs))
+                raise _RsError(coeffs=type(coeffs))
             n = 0
             for k, d in coeffs.items():
+                if not isinstance(k, _Rkey):
+                    raise _RsError(k=type(k))
                 if not isinstance(d, _Rdict):
-                    raise _RdictError(k=k, d=type(d))
+                    raise _RsError(k=k, d=type(d))
                 n += d.n
 
             ADict.__init__(self, coeffs)
             self.set_(ALorder=ALorder, n=n)  # in .validate
         except Exception as x:
-            raise _RdictError(ALorder=ALorder, cause=x)
+            raise _RsError(ALorder=ALorder, cause=x)
 
     def bnuz4(self):  # in .auxilats.__main__  # PYCHOK no cover
         # get C{(strB, number, unique, floatB)} rationals
@@ -53,8 +57,8 @@ class _Rcoeffs(ADict):
             t  =  R._tuple
             z += _zB(t)  # Float
             # assert R.Rdict is None
-            n += len(t)
-            u += sum(1 for f in t if f in _floats)
+            n +=  len(t)
+            u +=  sum(1 for f in t if f in _Cs)
         return b, n, (n - u), z
 
     def items(self):  # string-ify keys  # PYCHOK no cover
@@ -74,48 +78,48 @@ class _Rcoeffs(ADict):
 #       for R in self._Rtuples():
 #           assert isinstance(R, _Rtuple)
         if aL != a or lenAux != n:
-            raise _RdictError(aL=aL, ALorder=a, lenAux=lenAux, n=n)
+            raise _RsError(aL=aL, ALorder=a, lenAux=lenAux, n=n)
         return self
 
 
 class _Rdict(dict):  # in ._CX_#, .auxLat, .rhumb.aux_
     '''(INTERNAL) Dict of C{_Rtuple}s.
     '''
-    n = 0  # sum(R.k_n_k[1] for r in Rtuples)
+    n = 0  # sum(R.k_n_k[1] for R in Rtuples)
 
     def __init__(self, nt, *Rtuples):
         '''New C{_Rdict}.
         '''
         if not Rtuples:
-            raise _RdictError(Rtuples=MISSING)
+            raise _RsError(Rtuples=MISSING)
 
         for R in Rtuples:
             if not isinstance(R, _Rtuple):
-                raise _RdictError(R, R=type(R))
+                raise _RsError(R, R=type(R))
             k, n, _ = R.k_n_rs
             if k in self:
-                raise _RdictError(_duplicate_, k=k)
+                raise _RsError(_duplicate_, k=k)
             R.Rdict = self
             self[k] = R  # overwritten in self._floatuple
             self.n += n
         if self.n != nt:
-            raise _RdictError(n=n, nt=nt)
+            raise _RsError(n=n, nt=nt)
 
     def _floats(self, rs):
         # yield floats from a string of comma-separated rationals
         def _p_q(p=NN, q=1, *x):
             return (NN if x else p), q
 
-        _get = _floats.get
-        for r in NN(*rs.split()).split(_COMMA_):
+        _get = _Cs.get
+        for r in _MODS.basics._splituple(rs):
             p, q = _p_q(*r.split(_SLASH_))
             if p:
                 f = int(p) / int(q)  # fractions.Fraction?
                 if not isinstance(f, float):
-                    raise _RdictError(rs, f=f, p=p, q=q, r=r)
-                yield _get(f, f)  # from .constants?
+                    raise _RsError(rs, f=f, p=p, q=q, r=r)
+                yield _get(f, f)  # from .constants
             else:
-                raise _RdictError(rs, r=r)
+                raise _RsError(rs, r=r)
 
     def _floatuple(self, Rtuple):
         # return a tuple of floats from an C{_Rtuple}
@@ -125,13 +129,24 @@ class _Rdict(dict):  # in ._CX_#, .auxLat, .rhumb.aux_
         # @see: <https://StackOverflow.com/questions/10632839/>
         #       and <https://realPython.com/python-flatten-list/>
         if len(t) != n:
-            raise _RdictError(*rs, len=len(t), n=n)
-        self[k] = t
+            raise _RsError(*rs, len=len(t), n=n)
+        self[k] = t  # replace _Rtuple with tuple instance
         return t
 
 
-class _RdictError(_AssertionError):
-    '''(INTERNAL) For C{_Rdict} issues.
+class _Rkey(int):
+    '''(INTERNAL) For C{_Rcoeffs}, C{_Rdict} and C{_Rtuple} keys.
+    '''
+    def __new__(cls, k):
+        if not isinstance(k, int):
+            raise _RsError(k=type(k))
+        if not 0 <= k <= 8:  # 0.._MODS.auxilats.auxily.Aux.N + 2
+            raise _RsError(k=k)
+        return int.__new__(cls, k)
+
+
+class _RsError(_AssertionError):
+    '''(INTERNAL) For C{_Rcoeffs}, C{_Rdict or} C{_Rtuple} issues.
     '''
     def __init__(self, *rs, **kwds_cause):  # PYCHOK no cover
         if rs:
@@ -147,7 +162,7 @@ class _RdictError(_AssertionError):
 class _Rtuple(list):  # MUST be list, NOT tuple!
     '''(INTERNAL) I{Pseudo-tuple} of float rationals used in C{_Rdict}s.
     '''
-    Rdict  = None
+    Rdict  = None  # set by _Rdict.__init__
     k_n_rs = None, 0, ()
 
     def __init__(self, k, n, *rs):
@@ -156,20 +171,26 @@ class _Rtuple(list):  # MUST be list, NOT tuple!
            where C{p} and C{q} are C{int} digits only.
         '''
         try:
+            if not isinstance(k, _Rkey):
+                raise _RsError(k=type(k))
+            if not (isinstance(n, int) and n > 0):
+                raise _RsError(n=type(n))
             if not rs:
-                raise _RdictError(rs=MISSING)
+                raise _RsError(rs=MISSING)
             for t in rs:
                 if not isinstance(t, str):
-                    raise _RdictError(rs=type(t))
-            if not (isinstance(n, int) and n > 0):
-                raise _RdictError(n=type(n))
+                    raise _RsError(rs=type(t))
             self.k_n_rs = k, n, rs
         except Exception as x:
-            raise _RdictError(*rs, k=k, n=n, cause=x)
+            raise _RsError(*rs, k=k, n=n, cause=x)
 
     def __getitem__(self, i):
         return self._tuple[i]
 
+    if _MODS.sys_version_info2 < (3, 0):
+        def __getslice__(self, i, j):
+            return self._tuple[slice(i, j)]
+
     def __iter__(self):
         return iter(self._tuple)
 
@@ -178,26 +199,24 @@ class _Rtuple(list):  # MUST be list, NOT tuple!
 
     @Property_RO
     def _tuple(self):
-        # build the C{tuple} once, replace C{_Rdict}
-        # item at C{key} with the C{tuple} and fill
-        # this C{_Rlist} with the C{tuple} values
-        # for the initial __getitem__ retrieval[s]
+        # build the C{tuple} I{once} and replace
+        # C{_Rdict[key]} item with the C{tuple}
         try:
-            k, n, rs = self.k_n_rs
+            k, n, rs = self.k_n_rs  # for except ...
             t = self.Rdict._floatuple(self)
-            self[:] = t  # MUST copy into self!
+#           self[:] = t  # MUST copy into self?
         except Exception as x:
             if len(rs) > 1 and _QUOTE3_ in str(x):
                 rs = rs[0], _ELLIPSIS4_
-            raise _RdictError(k=k, n=n, rs=rs, cause=x)
+            raise _RsError(k=k, n=n, rs=rs, cause=x)
         del self.Rdict, self.k_n_rs  # trash refs
         return t
 
     def append(self, arg):
-        raise _RdictError(append=arg)
+        raise _RsError(append=arg)
 
     def extend(self, arg):
-        raise _RdictError(extend=arg)
+        raise _RsError(extend=arg)
 
 # **) MIT License
 #

pygeodesy/auxilats/__init__.py

@@ -29,7 +29,7 @@ from pygeodesy.lazily import _ALL_OTHER
 # no modules: auxAngle, auxDLat, auxDST, auxily, auxLat
 __all__ = _ALL_OTHER(Aux, AuxAngle, AuxDLat, AuxDST, AuxLat,
                      AuxBeta, AuxChi, AuxMu, AuxPhi, AuxTheta, AuxXi)
-__version__ = '24.09.04'
+__version__ = '25.01.15'
 
 # **) MIT License
 #

pygeodesy/auxilats/__main__.py

@@ -5,7 +5,7 @@ u'''Print L{auxilats} version, etc. using C{python -m pygeodesy.auxilats}.
 '''
 
 __all__ = ()
-__version__ = '24.09.04'
+__version__ = '25.04.04'
 
 
 def _main(**ALorder):  # PYCHOK no cover
@@ -36,7 +36,7 @@ def _main(**ALorder):  # PYCHOK no cover
         print(_usage(__file__))
 
 
-from sys import argv  # .internals._isPyChecker
+from sys import argv  # .internals._isPyChOK
 _main(ALorder=int(argv[1])) if len(argv) == 2 and argv[1].isdigit() else _main()
 
 # % python3.12 -m pygeodesy.auxilats 8

pygeodesy/auxilats/auxLat.py

@@ -18,7 +18,7 @@ from __future__ import division as _; del _  # PYCHOK semicolon
 from pygeodesy.auxilats.auxAngle import AuxAngle, AuxBeta, AuxChi, _AuxClass, \
                                         AuxMu, AuxPhi, AuxTheta, AuxXi
 from pygeodesy.auxilats.auxily import Aux, _sc, _sn
-from pygeodesy.auxilats._CX_Rs import _Rdict, _Rtuple
+from pygeodesy.auxilats._CX_Rs import _Rdict, _Rkey, _Rtuple
 from pygeodesy.basics import _reverange, _xinstanceof,  _passarg
 from pygeodesy.constants import INF, MAX_EXP, MIN_EXP, NAN, PI_2, PI_4, _EPSqrt, \
                                _0_0, _0_0s, _0_1, _0_5, _1_0, _2_0, _3_0, _360_0, \
@@ -48,7 +48,7 @@ except ImportError:  # Python 3.11-
         return pow(_2_0, x)
 
 __all__ = ()
-__version__ = '24.09.03'
+__version__ = '25.01.15'
 
 _TRIPS = 1024  # XXX 2 or 3?
 
@@ -829,16 +829,16 @@ def _Newton(tphi, Zeta, _toZeta, **name):
 
 
 _AR2Coeffs = _Rdict(18,
-    _Rtuple(4, 4, '4/315,  4/105,  4/15,  -1/3'),
-    _Rtuple(6, 6, '4/1287, 4/693,  4/15,   4/105, 4/315, -1/3'),
-    _Rtuple(8, 8, '4/3315, 4/2145, 4/1287, 4/693, 4/315,  4/105, 4/15, -1/3'))
+    _Rtuple(_Rkey(4), 4, '4/315,  4/105,  4/15,  -1/3'),
+    _Rtuple(_Rkey(6), 6, '4/1287, 4/693,  4/15,   4/105, 4/315, -1/3'),
+    _Rtuple(_Rkey(8), 8, '4/3315, 4/2145, 4/1287, 4/693, 4/315,  4/105, 4/15, -1/3'))
 
 _RRCoeffs = _Rdict(9,
-    _Rtuple(4, 2,  '1/64,    1/4'),
-    _Rtuple(6, 3,  '1/256,   1/64,  1/4'),
-    _Rtuple(8, 4, '25/16384, 1/256, 1/64, 1/4'))  # PYCHOK used!
+    _Rtuple(_Rkey(4), 2,  '1/64,    1/4'),
+    _Rtuple(_Rkey(6), 3,  '1/256,   1/64,  1/4'),
+    _Rtuple(_Rkey(8), 4, '25/16384, 1/256, 1/64, 1/4'))  # PYCHOK used!
 
-del _Rdict, _Rtuple
+del _Rdict, _Rkey, _Rtuple
 # assert set(_AR2Coeffs.keys()) == set(_RRCoeffs.keys())
 
 # AuxLat._Lmax = max(_AR2Coeffs.keys())  # == max(ALorder)

pygeodesy/auxilats/auxily.py

@@ -15,6 +15,7 @@ under the MIT/X11 License.  For more information, see the U{GeographicLib
 from __future__ import division as _; del _  # PYCHOK semicolon
 
 # from pygeodesy import auxilats  # _MODS
+from pygeodesy.auxilats._CX_Rs import _Rkey
 from pygeodesy.constants import INF, NAN, isinf, isnan, _0_0, _0_5, _1_0, \
                                _copysign_1_0, _over, _1_over
 from pygeodesy.errors import AuxError
@@ -26,19 +27,19 @@ from pygeodesy.utily import atan1
 from math import asinh, copysign
 
 __all__ = ()
-__version__ = '24.09.03'
+__version__ = '25.01.15'
 
 
 class Aux(object):
     '''Enum-style Aux names.
     '''
     _coeffs    = {}
-    GEOGRAPHIC = PHI   = GEODETIC = 0
-    PARAMETRIC = BETA  = REDUCED  = 1
-    GEOCENTRIC = THETA = 2  # all ...
-    RECTIFYING = MU    = 3  # use n^2
-    CONFORMAL  = CHI   = 4  # use n
-    AUTHALIC   = XI    = 5  # use n
+    GEOGRAPHIC = PHI   =  GEODETIC = _Rkey(0)
+    PARAMETRIC = BETA  =  REDUCED  = _Rkey(1)
+    GEOCENTRIC = THETA = _Rkey(2)  # all ...
+    RECTIFYING = MU    = _Rkey(3)  # use n^2
+    CONFORMAL  = CHI   = _Rkey(4)  # use n
+    AUTHALIC   = XI    = _Rkey(5)  # use n
     N  = 6
     N2 = 36
 
@@ -228,6 +229,7 @@ def _sn(tx):
 
 
 __all__ += _ALL_DOCS(Aux.__class__)
+del _Rkey
 
 # **) MIT License
 #

pygeodesy/basics.py

@@ -37,7 +37,7 @@ from math import copysign as _copysign
 # import inspect as _inspect  # _MODS
 
 __all__ = _ALL_LAZY.basics
-__version__ = '24.12.31'
+__version__ = '25.01.17'
 
 _below_           = 'below'
 _list_tuple_types = (list, tuple)
@@ -207,6 +207,22 @@ def _enumereverse(iterable):
         yield j, iterable[j]
 
 
+try:
+    from math import gcd as _gcd
+except ImportError:  # 3.4-
+
+    def _gcd(a, b):  # PYCHOK redef
+        # <https://WikiPedia.org/wiki/Greatest_common_divisor>
+        a, b = int(a), int(b)
+        if b > a:
+            a, b = b, a
+#       if b <= 0:
+#           return 1
+        while b:
+            a, b = b, (a % b)
+        return a
+
+
 def halfs2(str2):
     '''Split a string in 2 halfs.
 
@@ -222,6 +238,15 @@ def halfs2(str2):
     return str2[:h], str2[h:]
 
 
+def _integer_ratio2(x):  # PYCHOK no cover
+    '''(INTERNAL) Return C{B{x}.as_interger_ratio()}.
+    '''
+    try:  # int.as_integer_ratio in 3.8+
+        return  x.as_integer_ratio()
+    except (AttributeError, OverflowError, TypeError, ValueError):
+        return (x if isint(x) else float(x)), 1
+
+
 def int1s(x):  # PYCHOK no cover
     '''Count the number of 1-bits in an C{int}, I{unsigned}.
 
@@ -371,7 +396,7 @@ def isiterable(obj, strict=False):
        @return: C{True} if C{iterable}, C{False} otherwise.
     '''
     # <https://PyPI.org/project/isiterable/>
-    return isiterabletype(obj) if strict else hasattr(obj, '__iter__')  # map, range, set
+    return bool(isiterabletype(obj)) if strict else hasattr(obj, '__iter__')  # map, range, set
 
 
 def isiterablen(obj, strict=False):
@@ -383,7 +408,7 @@ def isiterablen(obj, strict=False):
        @return: C{True} if C{iterable} with C{len}gth, C{False} otherwise.
     '''
     _has = isiterabletype if strict else hasattr
-    return _has(obj, '__len__') and _has(obj, '__getitem__')
+    return bool(_has(obj, '__len__') and _has(obj, '__getitem__'))
 
 
 def isiterabletype(obj, method='__iter__'):
@@ -763,9 +788,13 @@ def _splituple(strs, *sep_splits):  # in .mgrs, ...
     '''(INTERNAL) Split a C{comma}- or C{whitespace}-separated
        string into a C{tuple} of stripped C{str}ings.
     '''
-    t = (strs.split(*sep_splits) if sep_splits else
-         strs.replace(_COMMA_, _SPACE_).split()) if strs else ()
-    return tuple(s.strip() for s in t if s)
+    if sep_splits:
+        t = (t.strip() for t in strs.split(*sep_splits))
+    else:
+        t = strs.strip()
+        if t:
+            t = t.replace(_COMMA_, _SPACE_).split()
+    return tuple(t) if t else ()
 
 
 def unsigned0(x):

pygeodesy/etm.py

@@ -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__ = '24.11.24'
+__version__ = '25.01.12'
 
 _OVERFLOW = _1_EPS**2  # ~2e+31
 _TAYTOL   =  pow(EPS,  0.6)
@@ -196,7 +196,7 @@ class Etm(Utm):
     def _toLLEB(self, unfalse=True, **unused):  # PYCHOK signature
         '''(INTERNAL) Compute (ellipsoidal) lat- and longitude.
         '''
-        xTM, d = self.exactTM, self.datum
+        d, xTM = self.datum, self.exactTM
         # double check that this and exactTM's ellipsoid match
         if xTM._E != d.ellipsoid:  # PYCHOK no cover
             t = repr(d.ellipsoid)
@@ -1022,7 +1022,7 @@ class ExactTransverseMercator(_NamedBase):
             g_k += atan1d(tau), degrees(lam)
         return g_k  # or (g, k, lat, lon)
 
-_allPropertiesOf_n(22, ExactTransverseMercator, Property_RO)  # PYCHOK assert
+_allPropertiesOf_n(22, ExactTransverseMercator, Property_RO)  # PYCHOK assert _ROs = ...
 del _0_1, _allPropertiesOf_n, EPS, _1_EPS, _EWGS84
 
 

pygeodesy/fmath.py

@@ -23,7 +23,7 @@ from math import fabs, sqrt  # pow
 import operator as _operator  # in .datums, .trf, .utm
 
 __all__ = _ALL_LAZY.fmath
-__version__ = '24.12.31'
+__version__ = '25.01.09'
 
 # sqrt(2) - 1 <https://WikiPedia.org/wiki/Square_root_of_2>
 _0_4142  =  0.41421356237309504880  # ... ~ 3730904090310553 / 9007199254740992
@@ -981,7 +981,7 @@ def _root(x, p, where):
             raise ValueError(_negative_)
     except Exception as X:
         raise _xError(X, unstr(where, x))
-    return _0_0
+    return _0_0 if p else _1_0
 
 
 def sqrt0(x, Error=None):