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

pygeodesy/fsums.py

@@ -39,8 +39,8 @@ results may differ from Python's C{math.fsum} results.
 # make sure int/int division yields float quotient, see .basics
 from __future__ import division as _; del _  # PYCHOK semicolon
 
-from pygeodesy.basics import isbool, iscomplex, isint, isscalar, \
-                            _signOf, itemsorted, signOf, _xiterable
+from pygeodesy.basics import _gcd, isbool, iscomplex, isint, isscalar, \
+                             _signOf, itemsorted, signOf, _xiterable
 from pygeodesy.constants import INF, INT0, MANT_DIG, NEG0, NINF, _0_0, \
                                _1_0, _N_1_0, _isfinite, _pos_self, \
                                Float, Int
@@ -64,20 +64,21 @@ from math import fabs, isinf, isnan, \
                  ceil as _ceil, floor as _floor  # PYCHOK used! .ltp
 
 __all__ = _ALL_LAZY.fsums
-__version__ = '24.12.02'
+__version__ = '25.01.12'
 
 from pygeodesy.interns import (
   _PLUS_     as _add_op_,  # in .auxilats.auxAngle
+  _DSLASH_   as _floordiv_op_,
   _EQUAL_    as _fset_op_,
   _RANGLE_   as _gt_op_,
   _LANGLE_   as _lt_op_,
   _PERCENT_  as _mod_op_,
   _STAR_     as _mul_op_,
   _NOTEQUAL_ as _ne_op_,
+  _DSTAR_    as _pow_op_,
   _DASH_     as _sub_op_,  # in .auxilats.auxAngle
   _SLASH_    as _truediv_op_
 )
-_floordiv_op_ = _truediv_op_ * 2  # _DSLASH_
 _divmod_op_   = _floordiv_op_ + _mod_op_
 _F2PRODUCT    = _getPYGEODESY('FSUM_F2PRODUCT')
 _iadd_op_     = _add_op_ + _fset_op_  # in .auxilats.auxAngle, .fstats
@@ -86,7 +87,6 @@ _isub_op_     = _sub_op_ + _fset_op_  # in .auxilats.auxAngle
 _NONFINITEr   = _0_0  # NOT INT0!
 _NONFINITES   = _getPYGEODESY('FSUM_NONFINITES')
 _non_zero_    = 'non-zero'
-_pow_op_      = _mul_op_ * 2  # _DSTAR_
 _RESIDUAL_0_0 = _getPYGEODESY('FSUM_RESIDUAL', _0_0)
 _significant_ = 'significant'
 _threshold_   = 'threshold'
@@ -96,7 +96,7 @@ def _2finite(x, _isfine=_isfinite):  # in .fstats
     '''(INTERNAL) return C{float(x)} if finite.
     '''
     return (float(x) if _isfine(x)  # and isscalar(x)
-            else _nfError(x))
+                   else _nfError(x))
 
 
 def _2float(index=None, _isfine=_isfinite, **name_x):  # in .fmath, .fstats
@@ -132,27 +132,23 @@ except  ImportError:  # PYCHOK DSPACE! Python 3.12-
 
     if _F2PRODUCT and _F2PRODUCT != _std_:
         # backward to PyGeodesy 24.09.09, with _fmaX
+        from pygeodesy.basics import _integer_ratio2
 
         def _fma(*a_b_c):  # PYCHOK no cover
             # mimick C{math.fma} from Python 3.13+,
             # the same accuracy, but ~14x slower
-            (na, da), (nb, db), (nc, dc) = map(_2n_d, a_b_c)
+            (na, da), (nb, db), (nc, dc) = map(_integer_ratio2, a_b_c)
             n  = na * nb * dc
             n += da * db * nc
             d  = da * db * dc
             try:
-                n, d = _n_d2(n, d)
+                n, d = _n_d2(n,  d)
                 r    = float(n / d)
             except OverflowError:  # "integer division result too large ..."
                 r = NINF if (_signOf(n, 0) * _signOf(d, 0)) < 0 else INF
             return r if _isfinite(r) else _fmaX(r, *a_b_c)  # "overflow in fma"
-
-        def _2n_d(x):  # PYCHOK no cover
-            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
     else:
+        _integer_ratio2 = None  # redef, in Fsum.is_math_fma
 
         def _fma(a, b, c):  # PYCHOK redef
             # mimick C{math.fma} from Python 3.13+,
@@ -161,17 +157,15 @@ except  ImportError:  # PYCHOK DSPACE! Python 3.12-
             r   = _fsum(_2products(a, b3s, c))
             return r if _isfinite(r) else _fmaX(r, a, b, c)
 
-        _2n_d = None  # redef
-
     def _fmaX(r, *a_b_c):  # PYCHOK no cover
-        # handle non-finite as Python 3.13+ C-function U{math_fma_impl<https://
-        # GitHub.com/python/cpython/blob/main/Modules/mathmodule.c#L2305>}:
-        #  raise a ValueError for a NAN result from non-NAN C{a_b_c}s or an
-        #  OverflowError for a non-NAN non-finite from all finite C{a_b_c}s.
+        # handle non-finite fma result as Python 3.13+ C-function U{math_fma_impl
+        # <https://GitHub.com/python/cpython/blob/main/Modules/mathmodule.c#L2305>}:
+        # raise a ValueError for a NAN result from non-NAN C{a_b_c}s otherwise an
+        # OverflowError for a non-finite, non-NAN result from all finite C{a_b_c}s.
         if isnan(r):
             def _x(x):
                 return not isnan(x)
-        else:  # non-NAN non-finite
+        else:  # non-finite, non-NAN
             _x = _isfinite
         if all(map(_x, a_b_c)):
             raise _nfError(r, unstr(_fma, *a_b_c))
@@ -179,10 +173,8 @@ except  ImportError:  # PYCHOK DSPACE! Python 3.12-
 
     def _2products(x, y3s, *zs):  # PYCHOK in _fma, ...
         # yield(x * y3 for y3 in y3s) + yield(z in zs)
-        # TwoProduct U{Algorithm 3.3
-        # <https://www.TUHH.De/ti3/paper/rump/OgRuOi05.pdf>}
-        # also in Python 3.13+ C{Modules/mathmodule.c} under
-        # #ifndef UNRELIABLE_FMA ... #else ... #endif
+        # TwoProduct U{Algorithm 3.3<https://www.TUHH.De/ti3/paper/rump/OgRuOi05.pdf>}, also
+        # in Python 3.13+ C{Modules/mathmodule.c} under #ifndef UNRELIABLE_FMA ... #else ...
         _, a, b = _2split3(x)
         for y, c, d in y3s:
             y *= x
@@ -282,25 +274,17 @@ def _isOK_or_finite(x, _isfine=_isfinite):
     return _isfine(x)  # C{bool}
 
 
-try:
-    from math import gcd as _gcd
-
-    def _n_d2(n, d):
-        '''(INTERNAL) Reduce C{n} and C{d} by C{gcd}.
-        '''
-        if n and d:
-            try:
-                c = _gcd(n, d)
-                if c > 1:
-                    n, d = (n // c), (d // c)
-            except TypeError:  # non-int float
-                pass
-        return n, d
-
-except ImportError:  # 3.4-
-
-    def _n_d2(*n_d):  # PYCHOK redef
-        return n_d
+def _n_d2(n, d):
+    '''(INTERNAL) Reduce C{n} and C{d} by C{gcd}.
+    '''
+    if n and d:
+        try:
+            c = _gcd(n, d)
+            if c > 1:
+                return (n // c), (d // c)
+        except TypeError:  # non-int float
+            pass
+    return n, d
 
 
 def _nfError(x, *args):
@@ -1332,16 +1316,16 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
 
     def _fhorner(self, x, cs, where, incx=True):  # in .fmath
         '''(INTERNAL) Add an L{Fhorner} evaluation of polynomial
-           C{sum(cs[i] * B{x}**i for i=0..len(cs)-1) if B{incx}
-           else sum(... i=len(cs)-1..0)}.
+           C{sum(c * B{x}**i for i, c in _e(cs))} where C{_e =
+           enumerate if B{incx} else _enumereverse}.
         '''
         # assert _xiterablen(cs)
         try:
-            n  = len(cs)
-            H  = self._Fsum_as(name__=self._fhorner)
-            _m = H._mul_Fsum if _isFsum_2Tuple(x) else \
-                 H._mul_scalar
-            if _2finite(x, **self._isfine) and n > 1:
+            n = len(cs)
+            if n > 1 and _2finite(x, **self._isfine):
+                H  = self._Fsum_as(name__=self._fhorner)
+                _m = H._mul_Fsum if _isFsum_2Tuple(x) else \
+                     H._mul_scalar
                 for c in (reversed(cs) if incx else cs):
                     H._fset(_m(x, _mul_op_), up=False)
                     H._fadd(c, up=False)
@@ -1804,8 +1788,8 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
             kwds.update(name_f2product_nonfinites_RESIDUAL)
         f = Fsum(**kwds)
         # assert all(v == self.__dict__[n] for n, v in f.__dict__.items())
-        return f._fset(xs[0], op=_fset_op_) if len(xs) == 1 else (
-               f._facc(xs, up=False) if xs else f)
+        return (f._facc(xs, up=False) if len(xs) > 1 else
+                f._fset(xs[0], op=_fset_op_)) if xs else f
 
     def fsum2(self, xs=(), **name):
         '''Add an iterable's items, summate and return the
@@ -1960,7 +1944,7 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
                     an C{fma} implementation as C{math.fma} or C{None}, a previous
                     C{PyGeodesy} implementation.
         '''
-        return (_2split3s is _passarg) or (False if _2n_d is None else None)
+        return (_2split3s is _passarg) or (False if _integer_ratio2 is None else None)
 
     def is_math_fsum(self):
         '''Are the summation functions L{fsum}, L{fsum_}, L{fsumf_}, L{fsum1},

pygeodesy/gars.py

@@ -14,14 +14,13 @@ Transcoded from I{Charles Karney}'s C++ class U{GARS
       (GARS)<https://Earth-Info.NGA.mil/GandG/coordsys/grids/gars.html>}.
 '''
 
-# from pygeodesy.basics import isstr  # from .named
+from pygeodesy.basics import isstr, _splituple
 from pygeodesy.constants import _off90, _1_over, _0_5, \
                                 _1_0  # PYCHOK used!
 from pygeodesy.errors import _ValueError, _xkwds, _xStrError
-from pygeodesy.interns import NN, _0to9_, _AtoZnoIO_, _COMMA_, \
-                             _INV_, _SPACE_
+from pygeodesy.interns import NN, _0to9_, _AtoZnoIO_, _INV_
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY
-from pygeodesy.named import _name__,  Fmt, isstr, Property_RO
+from pygeodesy.named import _name__,  Fmt, Property_RO
 from pygeodesy.namedTuples import LatLon2Tuple, LatLonPrec3Tuple
 # from pygeodesy.props import Property_RO  # from .named
 # from pygeodesy.streprs import Fmt  # from .named
@@ -30,7 +29,7 @@ from pygeodesy.units import Int_, Lat, Lon, Precision_, Scalar_, Str
 from math import floor
 
 __all__ = _ALL_LAZY.gars
-__version__ = '24.08.13'
+__version__ = '25.01.15'
 
 _Digits  = _0to9_
 _LatLen  =    2
@@ -136,7 +135,7 @@ class Garef(Str):
             if isinstance(lat_gll, Garef):
                 g, ll, p = str(lat_gll), lat_gll.latlon, lat_gll.precision
             elif isstr(lat_gll):
-                ll = lat_gll.replace(_COMMA_, _SPACE_).split()
+                ll = _splituple(lat_gll)
                 if len(ll) > 1:
                     g, ll, p = _encode3(ll[0], ll[1], precision)
                 else:

pygeodesy/geodesicx/__main__.py

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

pygeodesy/geohash.py

@@ -17,14 +17,13 @@ U{Geohashes<https://www.Movable-Type.co.UK/scripts/geohash.html>}.
       U{geohash-js<https://GitHub.com/DaveTroy/geohash-js>}.
 '''
 
-from pygeodesy.basics import isstr, map2
+from pygeodesy.basics import isstr, map2, _splituple
 from pygeodesy.constants import EPS, R_M, _0_0, _0_5, _180_0, _360_0, \
                                _90_0, _N_90_0, _N_180_0  # PYCHOK used!
 from pygeodesy.errors import _ValueError, _xkwds, _xStrError
 # from pygeodesy import formy as _formy  # _MODS
-from pygeodesy.interns import NN, _COMMA_, _DOT_, _E_, _height_, _N_, _NE_, \
-                             _NW_, _radius_, _S_, _SE_, _SPACE_, _SW_, _W_, \
-                             _width_  # _INV_
+from pygeodesy.interns import NN, _DOT_, _E_, _height_, _N_, _NE_, _NW_, \
+                             _radius_, _S_, _SE_, _SW_, _W_, _width_  # _INV_
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS
 from pygeodesy.named import _name__, _NamedDict, _NamedTuple, nameof, _xnamed
 from pygeodesy.namedTuples import Bounds2Tuple, Bounds4Tuple, LatLon2Tuple, \
@@ -38,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__ = '24.10.12'
+__version__ = '25.01.15'
 
 _formy   = _MODS.into(formy=__name__)
 _MASK5   =  16, 8, 4, 2, 1  # PYCHOK used!
@@ -290,7 +289,7 @@ class Geohash(Str):
             if isinstance(lat_ghll, Geohash):
                 gh, ll = str(lat_ghll), lat_ghll.latlon
             elif isstr(lat_ghll):  # "lat, lon" or "geohash"
-                ll = lat_ghll.replace(_COMMA_, _SPACE_).split()
+                ll = _splituple(lat_ghll)
                 if len(ll) > 1:
                     gh, ll = _GH.encode2(ll[0], ll[1], precision, eps)
                 else:

pygeodesy/geoids.py

@@ -89,7 +89,7 @@ courtesy of SBFRF.
 # make sure int/int division yields float quotient, see .basics
 from __future__ import division as _; del _  # PYCHOK semicolon
 
-from pygeodesy.basics import len2, min2, isodd, ub2str as _ub2str
+from pygeodesy.basics import len2, min2, isodd, _splituple, ub2str as _ub2str
 from pygeodesy.constants import EPS, _float as _F, _1_0, _N_90_0, _180_0, \
                                      _N_180_0, _360_0
 from pygeodesy.datums import Datums, _ellipsoidal_datum, _WGS84
@@ -105,7 +105,6 @@ from pygeodesy.interns import NN, _COLONSPACE_, _COMMASPACE_, _E_, _height_, \
                              _in_, _kind_, _lat_, _lon_, _mean_, _N_, _n_a_, \
                              _numpy_, _on_, _outside_, _S_, _s_, _scipy_, \
                              _SPACE_, _stdev_, _tbd_, _W_, _width_, _4_
-from pygeodesy.interns import _COMMA_  # PYCHOK used!
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS, _FOR_DOCS
 from pygeodesy.named import _name__, _Named, _NamedTuple
 # from pygeodesy.namedTuples import LatLon3Tuple  # _MODS
@@ -126,7 +125,7 @@ except ImportError:  # Python 3+
     from io import BytesIO as _BytesIO  # PYCHOK expected
 
 __all__ = _ALL_LAZY.geoids
-__version__ = '24.12.31'
+__version__ = '25.01.15'
 
 _assert_    = 'assert'
 _bHASH_     = b'#'
@@ -912,8 +911,7 @@ def _I(i):
 def _T(cs):
     '''(INTERNAL) Cache a tuple of single C{int} constants.
     '''
-    cs = cs.replace(_COMMA_, _SPACE_).strip()
-    return tuple(map(_I, cs.split()))
+    return tuple(map(_I, _splituple(cs)))
 
 _T0s12 = (_I(0),) * 12  # PYCHOK _T('0, 0, ..., 0')
 

pygeodesy/internals.py

@@ -345,11 +345,12 @@ def _isNix():  # in test/bases
     return _MODS.nix2[0]
 
 
-def _isPyChecker():  # PYCHOK no cover
-    '''(INTERNAL) Is C{PyChecker} running? (C{bool}).
+def _isPyChOK():  # PYCHOK no cover
+    '''(INTERNAL) Is C{PyChecker} running? (C{bool})
     '''
     # .../pychecker/checker.py --limit 0 --stdlib pygeodesy/<mod>/<name>.py
-    return _sys.argv[0].endswith('/pychecker/checker.py')
+    return _sys.argv[0].endswith('/pychecker/checker.py') or \
+            bool(_getPYGEODESY('PYCHOK'))
 
 
 def _isPyPy():  # in test/bases
@@ -701,7 +702,7 @@ def _versions(sep=_SPACE_):
 
 
 __all__ = tuple(map(_DUNDER_nameof, (machine, print_, printf)))
-__version__ = '24.11.06'
+__version__ = '25.04.04'
 
 if _is_DUNDER_main(__name__):  # PYCHOK no cover
 
@@ -714,7 +715,7 @@ if _is_DUNDER_main(__name__):  # PYCHOK no cover
     _main()
 
 # % python3 -m pygeodesy.internals
-# pygeodesy 24.11.11 Python 3.13.0 64bit arm64 macOS 14.6.1 _isfrozen False isLazy 1
+# pygeodesy 25.4.4 Python 3.13.2 64bit arm64 macOS 15.4 _isfrozen False isLazy 1
 
 # **) MIT License
 #

pygeodesy/interns.py

@@ -429,8 +429,8 @@ _SW_       = _S_  + _W_   # PYCHOK negative ones
 _DDOT_     = Str_(_DOT_   * 2)  # PYCHOK OK
 # _DEQUAL_ = Str_(_EQUAL_ * 2)  # PYCHOK OK
 # _DNL_    = Str_(_NL_    * 2)  # PYCHOK OK
-# _DSLASH_ = Str_(_SLASH_ * 2)  # PYCHOK OK
-# _DSTAR_  = Str_(_STAR_  * 2)  # PYCHOK OK
+_DSLASH_   = Str_(_SLASH_ * 2)  # PYCHOK OK
+_DSTAR_    = Str_(_STAR_  * 2)  # PYCHOK OK
 _DUNDER_   = Str_(_UNDER_ * 2)  # PYCHOK OK
 
 _LR_PAIRS  = {_LANGLE_:  _RANGLE_,
@@ -440,7 +440,7 @@ _LR_PAIRS  = {_LANGLE_:  _RANGLE_,
 
 __all__ = (_NN_,  # NOT MISSING!
             Str_.__name__)  # classes
-__version__ = '24.11.27'
+__version__ = '25.01.12'
 
 if __name__ == '__main__':
 

pygeodesy/karney.py

@@ -165,7 +165,7 @@ from pygeodesy.utily import atan2d, sincos2d, tand, _unrollon,  fabs
 # from math import fabs  # from .utily
 
 __all__ = _ALL_LAZY.karney
-__version__ = '24.11.26'
+__version__ = '25.04.08'
 
 _2_4_       = '2.4'
 _K_2_0      = _getenv(_PYGEODESY(_xgeographiclib, 1), _2_)
@@ -880,8 +880,9 @@ except ImportError:  # Python 3.12-
             s, t, _ = _sum3(s * x, t * x, c)
         return s + t
 
-#   def _poly_fma(x, *cs):
-#       S = Fhorner(x, *cs, incx=False)
+#   def _poly_fma(x, s, *cs):
+#       S  = Fhorner(x, *cs, incx=False)
+#       S += s
 #       return float(S)
 
 def _polynomial(x, cs, i, j):  # PYCHOK shared
@@ -893,8 +894,8 @@ def _polynomial(x, cs, i, j):  # PYCHOK shared
     if (i + 1) < j <= len(cs):  # load _Rtuple._tuple
         try:
             r = _wrapped.Math.polyval(j - i - 1, cs, i, x)
-        except AttributeError:
-            r = _poly_fma(x, *cs[i:j])
+        except AttributeError:  # no .Math
+            r = _poly_fma(x, _0_0, *cs[i:j])
     else:
         r = cs[i]
     return float(r)

pygeodesy/lazily.py

@@ -32,7 +32,7 @@ from pygeodesy import internals as _internals, interns as _interns, \
 # from pygeodesy.errors import _error_init, _xkwds_item2  # _ALL_MODS
 from pygeodesy.internals import _caller3, _DUNDER_nameof, _getPYGEODESY, _headof, \
                                 _is_DUNDER_main, printf, _tailof, _versions
-from pygeodesy.interns import NN, _attribute_, _by_, _COLONSPACE_, _COMMASPACE_, \
+from pygeodesy.interns import NN, _1_, _attribute_, _by_, _COLONSPACE_, _COMMASPACE_, \
                              _doesn_t_exist_, _DOT_, _DUNDER_all_, _EQUALSPACED_, \
                              _from_, _HASH_, _immutable_, _line_, _module_, _no_, \
                              _not_, _or_, _pygeodesy_abspath_, _pygeodesy_,  _sys, \
@@ -64,7 +64,7 @@ class LazyAttributeError(AttributeError):
 
 
 class LazyImportError(ImportError):
-    '''Raised if C{lazy import} is not supported, disabled or failed some other way.
+    '''Raised if C{lazy import} is not supported, disabled or failed.
     '''
     def __init__(self, *args, **kwds):
         _ALL_MODS.errors._error_init(ImportError, self, args, **kwds)
@@ -485,7 +485,7 @@ class _ALL_MODS(_internals._MODS_Base):
                 mod, dun = self.errors._xkwds_item2(mod_DUNDER_name)
                 _mod = _UNDER_(NN, mod)
                 d =  self.getmodule(dun)  # '__main__' OK
-                i = _getattribute(d, _mod, dun)
+                i = _getmodattr(d, _mod, dun)
                 assert isinstance(i, _Into)
                 m =  self.getmodule(mod)
                 setattr(d, _mod, m)  # overwrite C{d._mod}
@@ -509,7 +509,7 @@ class _ALL_MODS(_internals._MODS_Base):
 _internals._MODS = _ALL_MODS = _ALL_MODS()  # PYCHOK singleton
 
 __all__ = _ALL_LAZY.lazily
-__version__ = '25.01.05'
+__version__ = '25.01.25'
 
 
 def _ALL_OTHER(*objs):
@@ -592,7 +592,7 @@ def _getattras(attr_as):  # test/testDeprecated
     return as_ or a_.rstrip(_DOT_)
 
 
-def _getattribute(m, name, mod=_pygeodesy_):
+def _getmodattr(m, name, mod=_pygeodesy_):
     '''(INTERNAL) Get attr C{m.name}.
     '''
     try:
@@ -691,7 +691,7 @@ def _lazy_import2(pack):  # MCCABE 14
             if t not in sub_packages:  # invalid module package
                 raise LazyImportError(_DOT_(mod, _DUNDER_package_), t)
             if attr:  # get mod.attr
-                v = _getattribute(v, attr, mod)
+                v = _getmodattr(v, attr, mod)
 
         elif name in (_DUNDER_all_,):  # XXX _DUNDER_dir_, _DUNDER_members_?
             v = _ALL_INIT + tuple(imports.keys())
@@ -785,11 +785,11 @@ def _lazy_init2(pack):
     '''(INTERNAL) Initialize lazy import and set globals C{isLazy} and C{_unLazy0}.
 
        @arg pack: The name of the package (C{str}) performing the imports,
-                  to help resolving relative imports, usually C{__package__}.
+                  to resolve relative imports, usually C{__package__}.
 
        @return: 2-Tuple C{(package, parent)} with the importing C{package}
                 for easy reference within itself and its name aka the
-                C{parent}, same as B{C{pack}}.
+                C(package)'s C{parent}, same as B{C{pack}}.
 
        @raise LazyImportError: Lazy import not supported or not enabled,
                                an import failed or the package name is
@@ -799,12 +799,9 @@ def _lazy_init2(pack):
     '''
     global isLazy, _unLazy0
 
-    z = _getPYGEODESY('LAZY_IMPORT', None)
-    if z is None:  # not set, but ...
-        isLazy = 1  # ... default on 3.7+
-    else:
-        z = z.strip()  # like PYTHONVERBOSE et.al.
-        isLazy = int(z) if z.isdigit() else (1 if z else 0)
+    z = _getPYGEODESY('LAZY_IMPORT', _1_)  # 1 default on 3.7+
+    z =  z.strip()  # like PYTHONVERBOSE et.al.
+    isLazy = int(z) if z.isdigit() else (1 if z else 0)
 
     _unLazy0 = _unlazy or not isLazy  # pre-3.7 or w/o lazy import
 

pygeodesy/rhumb/__init__.py

@@ -9,7 +9,7 @@ u'''Package of lazily imported C{rhumb} modules L{rhumb.aux_}, L{rhumb.ekx} and
 from pygeodesy.lazily import _ALL_LAZY, _ALL_OTHER, _lazy_import_as, _unLazy0
 
 __all__ = _ALL_LAZY.rhumb
-__version__ = '24.11.07'
+__version__ = '25.01.15'
 
 if _unLazy0:  # or _isfrozen
     from pygeodesy.rhumb.aux_ import RhumbAux, RhumbLineAux

pygeodesy/rhumb/aux_.py

@@ -32,7 +32,7 @@ from pygeodesy.auxilats.auxAngle import AuxMu, AuxPhi,  hypot
 from pygeodesy.auxilats.auxDLat import AuxDLat, _DClenshaw
 # from pygeodesy.auxilats.auxDST import AuxDST  # _MODS
 from pygeodesy.auxilats.auxily import _Dlam, _Dp0Dpsi
-from pygeodesy.auxilats._CX_Rs import _Rdict, _Rtuple
+from pygeodesy.auxilats._CX_Rs import _Rdict, _Rkey, _Rtuple
 from pygeodesy.basics import copysign0, _reverange,  _xkwds_get1
 from pygeodesy.constants import EPS_2, MANT_DIG, PI4, isinf, \
                                _0_0, _4_0, _720_0, _log2, _over
@@ -48,7 +48,7 @@ from pygeodesy.rhumb.bases import RhumbBase, RhumbLineBase, \
 from math import ceil as _ceil, fabs, radians
 
 __all__ = _ALL_LAZY.rhumb_aux_
-__version__ = '24.09.02'
+__version__ = '25.01.15'
 
 # DIGITS = (sizeof(real) * 8) bits
 #        = (ctypes.sizeof(ctypes.c_double(1.0)) * 8) bits
@@ -314,25 +314,25 @@ def _RAseries(auxD):
 
 
 _RACoeffs = _Rdict(110,  # Rhumb Area Coefficients in matrix Q
- _Rtuple(4, 10,  # GEOGRAPHICLIB_RHUMBAREA_ORDER == 4
+ _Rtuple(_Rkey(4), 10,  # GEOGRAPHICLIB_RHUMBAREA_ORDER == 4
     '596/2025, -398/945, 22/45, -1/3',
     '1543/4725, -118/315, 1/5',
     '152/945, -17/315',
     '5/252'),
- _Rtuple(5, 15,  # GEOGRAPHICLIB_RHUMBAREA_ORDER == 5
+ _Rtuple(_Rkey(5), 15,  # GEOGRAPHICLIB_RHUMBAREA_ORDER == 5
     '-102614/467775, 596/2025, -398/945, 22/45, -1/3',
     '-24562/155925, 1543/4725, -118/315, 1/5',
     '-38068/155925, 152/945, -17/315',
     '-752/10395, 5/252',
     '-101/17325'),
- _Rtuple(6, 21,  # GEOGRAPHICLIB_RHUMBAREA_ORDER == 6
+ _Rtuple(_Rkey(6), 21,  # GEOGRAPHICLIB_RHUMBAREA_ORDER == 6
     '138734126/638512875, -102614/467775, 596/2025, -398/945, 22/45, -1/3',
     '17749373/425675250, -24562/155925, 1543/4725, -118/315, 1/5',
     '1882432/8513505, -38068/155925, 152/945, -17/315',
     '268864/2027025, -752/10395, 5/252',
     '62464/2027025, -101/17325',
     '11537/4054050'),
- _Rtuple(7, 28,  # GEOGRAPHICLIB_RHUMBAREA_ORDER == 7
+ _Rtuple(_Rkey(7), 28,  # GEOGRAPHICLIB_RHUMBAREA_ORDER == 7
     '-565017322/1915538625, 138734126/638512875, -102614/467775, 596/2025, -398/945, 22/45, -1/3',
     '-1969276/58046625, 17749373/425675250, -24562/155925, 1543/4725, -118/315, 1/5',
     '-58573784/638512875, 1882432/8513505, -38068/155925, 152/945, -17/315',
@@ -340,7 +340,7 @@ _RACoeffs = _Rdict(110,  # Rhumb Area Coefficients in matrix Q
     '-112832/1447875, 62464/2027025, -101/17325',
     '-4096/289575, 11537/4054050',
     '-311/525525'),
- _Rtuple(8, 36,  # GEOGRAPHICLIB_RHUMBAREA_ORDER == 8
+ _Rtuple(_Rkey(8), 36,  # GEOGRAPHICLIB_RHUMBAREA_ORDER == 8
     '188270561816/488462349375, -565017322/1915538625, 138734126/638512875, -102614/467775, 596/2025, -398/945, 22/45, -1/3',
     '2332829602/23260111875, -1969276/58046625, 17749373/425675250, -24562/155925, 1543/4725, -118/315, 1/5',
     '-41570288/930404475, -58573784/638512875, 1882432/8513505, -38068/155925, 152/945, -17/315',
@@ -350,7 +350,7 @@ _RACoeffs = _Rdict(110,  # Rhumb Area Coefficients in matrix Q
     '4193792/723647925, -311/525525',
     '1097653/1929727800')
 )
-del _Rdict, _Rtuple
+del _Rdict, _Rkey, _Rtuple
 
 __all__ += _ALL_DOCS(Caps)