pygeodesy 25.9.9__py2.py3-none-any.whl → 25.11.5__py2.py3-none-any.whl

pygeodesy/elliptic.py

@@ -100,7 +100,7 @@ from math import asin, asinh, atan, ceil, cosh, fabs, floor, radians, \
                  sin, sinh, sqrt, tan, tanh  # tan as _tan
 
 __all__ = _ALL_LAZY.elliptic
-__version__ = '25.09.04'
+__version__ = '25.10.10'
 
 _TolRD  =  zqrt(EPS * 0.002)
 _TolRF  =  zqrt(EPS * 0.030)
@@ -932,7 +932,7 @@ class Elliptic(_Named):
             raise _ellipticError(Elliptic._RFRD, x, y, z, m, cause=X)
         return float(R)
 
-_allPropertiesOf_n(16, Elliptic)  #  # PYCHOK assert, see Elliptic.reset
+_allPropertiesOf_n(16, Elliptic)  # PYCHOK assert, see Elliptic.reset
 
 
 class EllipticError(_ValueError):

pygeodesy/errors.py

@@ -28,7 +28,7 @@ from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS, _PYTHON_X_DEV
 from copy import copy as _copy
 
 __all__ = _ALL_LAZY.errors  # _ALL_DOCS('_InvalidError', '_IsnotError')  _under
-__version__ = '25.09.04'
+__version__ = '25.10.30'
 
 _argument_   = 'argument'
 _basics      = _MODS.into(basics=__name__)

pygeodesy/etm.py

@@ -67,7 +67,7 @@ from pygeodesy.basics import _isin, map1, neg, neg_, _xinstanceof
 from pygeodesy.constants import EPS, EPS02, PI_2, PI_4, _K0_UTM, \
                              _1_EPS, _0_0, _0_1, _0_5, _1_0, _2_0, \
                              _3_0, _90_0, isnear0, isnear90
-from pygeodesy.constants import _4_0  # PYCHOK used!
+from pygeodesy.constants import _0_75, _4_0  # PYCHOK used!
 from pygeodesy.datums import _ellipsoidal_datum, _WGS84,  _EWGS84
 # from pygeodesy.ellipsoids import _EWGS84  # from .datums
 # from pygeodesy.elliptic import Elliptic  # _MODS
@@ -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__ = '25.08.31'
+__version__ = '25.09.15'
 
 _OVERFLOW = _1_EPS**2  # ~2e+31
 _TAYTOL   =  pow(EPS,  0.6)
@@ -410,7 +410,7 @@ class ExactTransverseMercator(_NamedBase):
     def _Ev_3cKE_4(self):
         '''(INTERNAL) Get and cache C{_Ev.cKE * 3 / 4}.
         '''
-        return self._Ev_cKE * 0.75  # _0_75
+        return self._Ev_cKE * _0_75
 
     @Property_RO
     def _Ev_5cKE_4(self):

pygeodesy/fmath.py

@@ -10,8 +10,8 @@ from __future__ import division as _; del _  # noqa: E702 ;
 from pygeodesy.basics import _copysign, copysign0, isbool, isint, isodd, \
                               isscalar, len2, map1, _xiterable,  typename
 from pygeodesy.constants import EPS0, EPS02, EPS1, NAN, PI, PI_2, PI_4, \
-                               _0_0, _0_125, _1_6th, _0_25, _1_3rd, _0_5, _1_0, \
-                               _1_5, _copysign_0_0, isfinite, remainder
+                               _0_0, _0_125, _0_25, _0_5, _1_0, _1_5, \
+                               _copysign_0_0, isfinite, remainder
 from pygeodesy.errors import _IsnotError, LenError, _TypeError, _ValueError, \
                              _xError, _xkwds, _xkwds_pop2, _xsError
 from pygeodesy.fsums import _2float, Fsum, fsum, _isFsum_2Tuple,  Fmt, unstr
@@ -25,10 +25,12 @@ from math import fabs, sqrt  # pow
 import operator as _operator  # in .datums, .trf, .utm
 
 __all__ = _ALL_LAZY.fmath
-__version__ = '25.08.31'
+__version__ = '25.09.15'
 
 # sqrt(2) - 1 <https://WikiPedia.org/wiki/Square_root_of_2>
 _0_4142  =  0.41421356237309504880  # ~ 3_730_904_090_310_553 / 9_007_199_254_740_992
+_1_3rd   = _1_0 / 3
+_1_6th   = _1_0 / 6
 _2_3rd   = _1_3rd * 2
 _h_lt_b_ = 'abs(h) < abs(b)'
 

pygeodesy/geodesicx/gx.py

@@ -65,7 +65,7 @@ from pygeodesy.utily import atan2, atan2d as _atan2d_reverse, _unrollon, \
 from math import copysign, cos, degrees, fabs, radians, sqrt
 
 __all__ = ()
-__version__ = '25.09.02'
+__version__ = '25.09.16'
 
 _MAXIT1 = 20
 _MAXIT2 = 10 + _MAXIT1 + MANT_DIG  # MANT_DIG == C++ digits
@@ -188,7 +188,7 @@ class GeodesicExact(_GeodesicBase):
                  <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1GeodesicExact.html>}
                  and Python U{Geodesic.ArcDirect<https://GeographicLib.SourceForge.io/Python/doc/code.html>}.
         '''
-        return self._GDictDirect(lat1, lon1, azi1, True, a12, outmask)
+        return self._GDictDirect(lat1, lon1, azi1, True, a12, outmask=outmask)
 
     def ArcDirectLine(self, lat1, lon1, azi1, a12, caps=Caps.ALL, **name):
         '''Define a L{GeodesicLineExact} in terms of the I{direct} geodesic problem and as arc length.
@@ -348,7 +348,7 @@ class GeodesicExact(_GeodesicBase):
                  <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1GeodesicExact.html>}
                  and Python U{Geodesic.Direct<https://GeographicLib.SourceForge.io/Python/doc/code.html>}.
         '''
-        return self._GDictDirect(lat1, lon1, azi1, False, s12, outmask)
+        return self._GDictDirect(lat1, lon1, azi1, False, s12, outmask=outmask)
 
     def Direct3(self, lat1, lon1, azi1, s12):  # PYCHOK outmask
         '''Return the destination lat, lon and reverse azimuth
@@ -356,7 +356,7 @@ class GeodesicExact(_GeodesicBase):
 
            @return: L{Destination3Tuple}C{(lat, lon, final)}.
         '''
-        r = self._GDictDirect(lat1, lon1, azi1, False, s12, Caps._AZIMUTH_LATITUDE_LONGITUDE)
+        r = self._GDictDirect(lat1, lon1, azi1, False, s12, outmask=Caps._AZIMUTH_LATITUDE_LONGITUDE)
         return Destination3Tuple(r.lat2, r.lon2, r.azi2)  # no iteration
 
     def _DirectLine(self, ll1, azi12, s12=0, **caps_name):
@@ -710,7 +710,7 @@ class GeodesicExact(_GeodesicBase):
            @return: L{Direct9Tuple}C{(a12, lat2, lon2, azi2,
                                       s12, m12,  M12,  M21, S12)}.
         '''
-        r = self._GDictDirect(lat1, lon1, azi1, arcmode, s12_a12, outmask)
+        r = self._GDictDirect(lat1, lon1, azi1, arcmode, s12_a12, outmask=outmask)
         return r.toDirect9Tuple()
 
     def _GenInverse(self, lat1, lon1, lat2, lon2, outmask=Caps.STANDARD):
@@ -719,7 +719,7 @@ class GeodesicExact(_GeodesicBase):
            @return: L{Inverse10Tuple}C{(a12, s12, salp1, calp1, salp2, calp2,
                                              m12,   M12,   M21,   S12)}.
         '''
-        r = self._GDictInverse(lat1, lon1, lat2, lon2, outmask | Caps._SALP_CALPs_)
+        r = self._GDictInverse(lat1, lon1, lat2, lon2, outmask=outmask | Caps._SALP_CALPs_)
         return r.toInverse10Tuple()
 
     def _Inverse(self, ll1, ll2, wrap, **outmask):
@@ -752,7 +752,7 @@ class GeodesicExact(_GeodesicBase):
                  <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1GeodesicExact.html>} and
                  Python U{Geodesic.InverseLine<https://GeographicLib.SourceForge.io/Python/doc/code.html>}.
         '''
-        return self._GDictInverse(lat1, lon1, lat2, lon2, outmask)
+        return self._GDictInverse(lat1, lon1, lat2, lon2, outmask=outmask)
 
     def Inverse1(self, lat1, lon1, lat2, lon2, wrap=False):
         '''Return the non-negative, I{angular} distance in C{degrees}.
@@ -764,7 +764,8 @@ class GeodesicExact(_GeodesicBase):
         # and .HeightIDWkarney._distances
         if wrap:
             _, lat2, lon2 = _Wrap.latlon3(lat1, lat2, lon2, True)  # _Geodesic.LONG_UNROLL
-        return fabs(self._GDictInverse(lat1, lon1, lat2, lon2, Caps.EMPTY).a12)  # a12 always
+        r = self._GDictInverse(lat1, lon1, lat2, lon2, outmask=Caps.EMPTY)
+        return fabs(r.a12)  # a12 always
 
     def Inverse3(self, lat1, lon1, lat2, lon2):  # PYCHOK outmask
         '''Return the distance in C{meter} and the forward and
@@ -772,7 +773,7 @@ class GeodesicExact(_GeodesicBase):
 
            @return: L{Distance3Tuple}C{(distance, initial, final)}.
         '''
-        r = self._GDictInverse(lat1, lon1, lat2, lon2, Caps.AZIMUTH_DISTANCE)
+        r = self._GDictInverse(lat1, lon1, lat2, lon2, outmask=Caps.AZIMUTH_DISTANCE)
         return Distance3Tuple(r.s12, wrap360(r.azi1), wrap360(r.azi2),
                               iteration=r.iteration)
 
@@ -804,7 +805,7 @@ class GeodesicExact(_GeodesicBase):
                  <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1GeodesicExact.html>} and
                  Python U{Geodesic.InverseLine<https://GeographicLib.SourceForge.io/Python/doc/code.html>}.
         '''
-        r = self._GDictInverse(lat1, lon1, lat2, lon2, caps | Caps._SALP_CALPs_)
+        r = self._GDictInverse(lat1, lon1, lat2, lon2, outmask=caps | Caps._SALP_CALPs_)
         return GeodesicLineExact(self, lat1, lon1, None, caps=caps, _s_calp1=(r.salp1, r.calp1),
                                                             **name)._GenSet(self._debug, **r)
 

pygeodesy/geodesicx/gxarea.py

@@ -31,7 +31,7 @@ from pygeodesy.props import Property, Property_RO, property_RO
 from math import fmod as _fmod
 
 __all__ = ()
-__version__ = '25.06.04'
+__version__ = '25.09.16'
 
 
 class GeodesicAreaExact(_NamedBase):
@@ -74,7 +74,7 @@ class GeodesicAreaExact(_NamedBase):
         self._g_gX = g = geodesic
         # use the class-level Caps since the values
         # differ between GeodesicExact and Geodesic
-        self._mask = g.DISTANCE | g.LATITUDE | g.LONGITUDE
+        self._mask =  g.DISTANCE | g.LATITUDE | g.LONGITUDE
         self._Peri = _Accumulator(name='_Peri')
         if not polyline:  # perimeter and area
             self._mask |=  g.AREA | g.LONG_UNROLL
@@ -188,7 +188,7 @@ class GeodesicAreaExact(_NamedBase):
         '''(INTERNAL) Edge helper.
         '''
         lon1 = self.lon1
-        r = self._g_gX._GDictDirect(self.lat1, lon1, azi, False, s, self._mask)
+        r = self._g_gX._GDictDirect(self.lat1, lon1, azi, False, s, outmask=self._mask)
         if self._Area:  # aka transitDirect
             # Count crossings of prime meridian exactly as
             # int(ceil(lon2 / 360)) - int(ceil(lon1 / 360))
@@ -219,7 +219,7 @@ class GeodesicAreaExact(_NamedBase):
     def _Inverse(self, lat1, lon1, lat2, lon2):
         '''(INTERNAL) Point helper.
         '''
-        r = self._g_gX._GDictInverse(lat1, lon1, lat2, lon2, self._mask)
+        r = self._g_gX._GDictInverse(lat1, lon1, lat2, lon2, outmask=self._mask)
         if self._Area:  # aka transit
             # count crossings of prime meridian as +1 or -1
             # if in east or west direction, otherwise 0

pygeodesy/geohash.py

@@ -37,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__ = '25.04.21'
+__version__ = '25.09.16'
 
 _formy   = _MODS.into(formy=__name__)
 _MASK5   =  16, 8, 4, 2, 1  # PYCHOK used!
@@ -150,10 +150,10 @@ class _GH(object):
         '''
         def _encodes(lat, lon, prec, eps=0):
             s, w, n, e = self.SWNE4
-            E, d, _mid = self.EncodeB32, True, _2mid
+            d, _mid = True, _2mid
             for _ in range(prec):
                 i = 0
-                for _ in range(5):  # len(_MASK5)
+                for _ in _MASK5:
                     i += i
                     if d:  # bisect longitude
                         a = _mid(e, w)
@@ -170,7 +170,7 @@ class _GH(object):
                             s  = a
                             i += 1
                     d = not d
-                yield E[i]
+                yield self.EncodeB32[i]
                 if eps > 0:  # infer prec
                     if _2dab(lon, e, w) < eps and \
                        _2dab(lat, n, s) < eps:
@@ -233,10 +233,10 @@ class _GH(object):
         if not (0 < nc <= _MaxPrec):  # or geohash.startswith(_INV_)
             raise GeohashError(geohash=geohash, len=nc)
         s, w, n, e = self.SWNE4
-        D, d, _mid = self.DecodeB32, True, _2mid
+        d, _mid = True, _2mid
         try:
             for j, c in enumerate(geohash.lower()):
-                i = D[c]
+                i = self.DecodeB32[c]
                 for m in mask5:
                     if d:  # longitude
                         a = _mid(e, w)

pygeodesy/geoids.py

@@ -28,7 +28,7 @@ or L{GeoidPGM} and download a C{geoid} model file, containing locations with
 known heights also referred to as the C{grid knots}.  See the documentation of
 the interpolator class for references to available C{grid} models.
 
-C{>>> from pygeodesy import GeoidEGM96  # or -G2012B, -Karney or -PGM as GeoidXyz}
+C{>>> from pygeodesy import GeoidEGM96 as GeoidXyz  # or GeoidG2012B, -Karney or -PGM}
 
 2. Instantiate an interpolator with the C{geoid} model file and use keyword
 arguments to select different interpolation options
@@ -126,7 +126,7 @@ except ImportError:  # Python 3+
     from io import BytesIO as _BytesIO  # PYCHOK expected
 
 __all__ = _ALL_LAZY.geoids
-__version__ = '25.05.21'
+__version__ = '25.09.26'
 
 _assert_    = 'assert'
 _bHASH_     = b'#'

pygeodesy/heights.py

@@ -21,7 +21,7 @@ C{>>> ...}
 
 2. Select one of the C{Height} classes for height interpolation
 
-C{>>> from pygeodesy import HeightCubic  # or other Height... as HeightXyz}
+C{>>> from pygeodesy import HeightCubic as HeightXyz  # or an other Height... class}
 
 3. Instantiate a height interpolator with the C{knots} and use keyword
 arguments to select different interpolation options
@@ -92,7 +92,7 @@ from pygeodesy.units import _isDegrees, Float_, Int_
 # from math import radians  # from .points
 
 __all__ = _ALL_LAZY.heights
-__version__ = '25.05.26'
+__version__ = '25.09.29'
 
 _error_  = 'error'
 _formy   = _MODS.into(formy=__name__)

pygeodesy/internals.py

@@ -202,7 +202,7 @@ class _MODS_Base(object):
                      (_macOS_,   p.mac_ver),
                      (_Windows_, p.win32_ver),
                      (_Nix,  _MODS.nix_ver),
-                     ('Java',    p.java_ver),
+# removed Py 3.15    ('Java',    p.java_ver),
                      ('uname',   p.uname)):
             v = v()[0]
             if v and n:
@@ -498,7 +498,7 @@ def printf(fmt, *args, **nl_nt_prec_prefix__end_file_flush_sep__kwds):
             t =  fmt
     except Exception as x:
         _Error, s = _MODS.errors._xError2(x)
-        _unstr    = _MODS.strepr.unstr
+        _unstr    = _MODS.streprs.unstr
         t = _unstr(printf, fmt, *args, **nl_nt_prec_prefix__end_file_flush_sep__kwds)
         raise _Error(s, txt=t, cause=x)
     try:
@@ -659,11 +659,10 @@ def _usage_argv(argv0, *args):
     '''(INTERNAL) Return 3-tuple C{(python, '-m', module, *args)}.
     '''
     o = _MODS.os
-    m =  o.path.dirname(argv0)
-    m =  m.replace(o.getcwd(), _ELLIPSIS_) \
-          .replace(o.sep,      _DOT_).strip()
-    b =  o.path.basename(argv0)
-    b, x = o.path.splitext(b)
+    p =  o.path
+    m =  p.dirname(argv0).replace(o.getcwd(), _ELLIPSIS_) \
+                         .replace(o.sep,      _DOT_).strip()
+    b, x = p.splitext(p.basename(argv0))
     if x == '.py' and b != _DMAIN_:
         m = _DOT_(m or _pygeodesy_, b)
     p = NN(_python_, _MODS.sys_version_info2[0])
@@ -710,7 +709,7 @@ def _versions(sep=_SPACE_):
 
 
 __all__ = tuple(map(typename, (machine, print_, printf, typename)))
-__version__ = '25.09.05'
+__version__ = '25.10.26'
 
 if __name__ == _DMAIN_:
 

pygeodesy/karney.py

@@ -145,8 +145,8 @@ from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import _copysign, _isin, isint, neg, unsigned0, \
                              _xgeographiclib, _zip
-from pygeodesy.constants import NAN, _isfinite as _math_isfinite, _0_0, \
-                               _1_16th, _1_0, _2_0, _180_0, _N_180_0, _360_0
+from pygeodesy.constants import NAN, _isfinite as _math_isfinite, \
+                               _0_0, _1_0, _2_0, _180_0, _N_180_0, _360_0
 from pygeodesy.errors import GeodesicError, _ValueError, _xkwds
 from pygeodesy.fmath import cbrt, fremainder, norm2  # Fhorner, Fsum
 from pygeodesy.internals import _getenv, _popen2, _PYGEODESY_ENV, typename, \
@@ -165,8 +165,9 @@ from pygeodesy.utily import atan2d, sincos2d, tand, _unrollon,  fabs
 # from math import fabs  # from .utily
 
 __all__ = _ALL_LAZY.karney
-__version__ = '25.08.31'
+__version__ = '25.09.13'
 
+_1_16th     = _1_0 / 16
 _2_4_       = '2.4'
 _K_2_0      = _getenv(_PYGEODESY_ENV(typename(_xgeographiclib)[2:]), _2_)
 _K_2_4      = _K_2_0 ==  _2_4_
@@ -723,11 +724,10 @@ def _around(x):  # in .utily.sincos2d
     try:
         return _wrapped.Math.AngRound(x)
     except AttributeError:
-        b, a = _1_16th, fabs(x)
-        if a < b:
-            a -= b
-            a += b
-            x = _copysign(a, x)
+        z = _1_16th
+        w =  z - fabs(x)
+        if w > 0:  # don't "simplify" z - (z - x) to x
+            x = _copysign(z - w, x)
         return x
 
 
@@ -772,13 +772,13 @@ def _diff182(deg0, deg, K_2_0=False):
     try:
         return _wrapped.Math.AngDiff(deg0, deg)
     except AttributeError:
-        if K_2_0 or _K_2_0:  # geographiclib 2.0
-            _r, _360 = fremainder, _360_0
-            d, t = _sum2(_r(-deg0, _360),
-                         _r( deg,  _360))
-            d, t = _sum2(_r( d,    _360), t)
+        if K_2_0 or _K_2_0:  # geographiclib 2.0+
+            _r = fremainder
+            d, t = _sum2(_r(-deg0, _360_0),
+                         _r( deg,  _360_0))
+            d, t = _sum2(_r( d,    _360_0), t)
             if _isin(d, _0_0, _180_0, _N_180_0):
-                d = _copysign(d, -t if t else (deg - deg0))
+                d = _copysign(d, (-t) if t else (deg - deg0))
         else:
             _n = _norm180
             d, t = _sum2(_n(-deg0), _n(deg))

pygeodesy/lazily.py

@@ -203,7 +203,7 @@ _ALL_LAZY = _NamedEnum_RO(_name='_ALL_LAZY',
                                    'INF', 'INT0', 'MANT_DIG', 'MAX', 'MAX_EXP', 'MIN', 'MIN_EXP', 'NAN', 'NEG0', 'NINF',
                                    'PI', 'PI2', 'PI_2', 'PI3', 'PI_3', 'PI3_2', 'PI4', 'PI_4', 'PI_6',
                                    'R_FM', 'R_GM', 'R_KM', 'R_M', 'R_MA', 'R_MB', 'R_NM', 'R_QM', 'R_SM', 'R_VM',
-                                   'float_', 'float0_', 'isclose', 'isfinite', 'isinf', 'isint0',
+                                   'float_', 'float0_', 'floats_', 'isclose', 'isfinite', 'isinf', 'isint0',
                                    'isnan', 'isnear0', 'isnear1', 'isnear90', 'isneg', 'isneg0', 'isninf', 'isnon0',
                                    'remainder'),
                          datums=_a('Datum', 'Datums', 'Transform', 'Transforms'),
@@ -346,9 +346,9 @@ _ALL_LAZY = _NamedEnum_RO(_name='_ALL_LAZY',
                                    'lrstrip', 'pairs', 'reprs', 'strs', 'unstr'),
                             trf=_a('RefFrame', 'RefFrames', 'TransformXform', 'TRFXform', 'TRFXform7Tuple',
                                    'date2epoch', 'epoch2date', 'trfTransform0', 'trfTransforms', 'trfXform'),
-                      triaxials=_a('BetaOmega2Tuple', 'BetaOmega3Tuple', 'Jacobi2Tuple',
-                                   'JacobiConformal', 'JacobiConformalSpherical',
-                                   'Triaxial', 'Triaxial_', 'TriaxialError', 'Triaxials', 'hartzell4'),
+                      triaxials=_a('BetaOmega2Tuple', 'BetaOmega3Tuple',
+                                   'ConformalSphere', 'ConformalTriaxial', 'Conformal2Tuple',
+                                   'Triaxial', 'Triaxial_', 'TriaxialError', 'Triaxials', 'hartzell4'),  # 'height4'
                           units=_a('Azimuth', 'Band', 'Bearing', 'Bearing_', 'Bool',
                                    'Degrees', 'Degrees_', 'Degrees2', 'Distance', 'Distance_', 'Easting', 'Epoch',
                                    'Feet', 'FIx', 'Float_', 'Height', 'Height_', 'HeightX', 'Int_',
@@ -394,6 +394,7 @@ _ALL_DEPRECATED = _NamedEnum_RO(_name='_ALL_DEPRECATED',
                                          'HausdorffCosineAndoyerLambert', 'HausdorffCosineForsytheAndoyerLambert',
                                          'HeightIDW', 'HeightIDW2', 'HeightIDW3', 'HeightIDWcosineAndoyerLambert',
                                          'HeightIDWcosineForsytheAndoyerLambert', 'Helmert7Tuple',
+                                         'JacobiConformal', 'JacobiConformalSpherical', 'Jacobi2Tuple',
                                          'Lam_', 'LatLonExact4Tuple', 'NearestOn4Tuple', 'Ned3Tuple',
                                          'Phi_', 'RefFrameError', 'Rhumb7Tuple', 'RhumbOrder2Tuple',
                                          'Transform7Tuple', 'TriAngle4Tuple', 'UtmUps4Tuple', 'XDist'),
@@ -511,7 +512,7 @@ class _ALL_MODS(_internals._MODS_Base):
 _internals._MODS = _ALL_MODS = _ALL_MODS()  # PYCHOK singleton
 
 __all__ = _ALL_LAZY.lazily
-__version__ = '25.09.09'
+__version__ = '25.10.30'
 
 
 def _ALL_OTHER(*objs):

pygeodesy/named.py

@@ -1247,7 +1247,7 @@ def modulename(clas, prefixed=None):  # in .basics._xversion
     '''Return the class name optionally prefixed with the
        module name.
 
-       @arg clas: The class (any C{class}).
+       @arg clas: The class (any C{class} or C{str}).
        @kwarg prefixed: Include the module name (C{bool}), see
                         function C{classnaming}.
 

pygeodesy/streprs.py

@@ -22,7 +22,7 @@ from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS
 from math import fabs, log10 as _log10
 
 __all__ = _ALL_LAZY.streprs
-__version__ = '25.04.14'
+__version__ = '25.10.10'
 
 _at_        = 'at'         # PYCHOK used!
 _EN_PREC    =  6           # max MGRS/OSGR precision, 1 micrometer
@@ -545,22 +545,22 @@ def unstr(where, *args, **kwds_):
        @arg where: Class, function, method (C{type}) or name (C{str}).
        @arg args: Optional positional arguments.
        @kwarg kwds_: Optional keyword arguments, except C{B{_Cdot}=None},
-                     C{B{_ELLIPSIS}=False} and C{B{_fmt}=Fmt.g}.
+                     C{B{_ELLIPSIS}=False}, C{B{_fmt}=Fmt.g} and C{B{_prec}=6}.
 
        @return: Representation (C{str}).
     '''
-    def _C_e_g_kwds3(_Cdot=None, _ELLIPSIS=0, _fmt=Fmt.g, **kwds):
-        return _Cdot, _ELLIPSIS, _fmt, kwds
+    def _C_e_fmt_prec_kwds(_Cdot=None, _ELLIPSIS=0, _fmt=Fmt.g, _prec=6, **kwds):
+        return _Cdot, _ELLIPSIS, dict(fmt=_fmt, prec=_prec), kwds
 
-    C, e, g, kwds = _C_e_g_kwds3(**kwds_)
+    C, e, fmt_prec, kwds = _C_e_fmt_prec_kwds(**kwds_)
     if e and len(args) > (e + 1):
-        t  =  reprs(args[:e], fmt=g)
+        t  =  reprs(args[:e], **fmt_prec)
         t += _ELLIPSIS_,
-        t +=  reprs(args[-1:], fmt=g)
+        t +=  reprs(args[-1:], **fmt_prec)
     else:
-        t  =  reprs(args, fmt=g) if args else ()
+        t  =  reprs(args, **fmt_prec) if args else ()
     if kwds:
-        t += pairs(itemsorted(kwds), fmt=g)
+        t += pairs(itemsorted(kwds), **fmt_prec)
     n = where if isstr(where) else typename(where)  # _NN_
     if C and hasattr(C, n):
         try:  # bound method of class C?

pygeodesy/trf.py

@@ -71,7 +71,7 @@ en/how-to-deal-with-etrs89-datum-and-time-dependent-transformation-parameters-45
 
 from pygeodesy.basics import _isin, map1, neg, isidentifier, isstr, _xinstanceof, \
                              _xscalar,  typename
-from pygeodesy.constants import _float as _F, _0_0s, _0_0, _0_001, _0_5, _1_0
+from pygeodesy.constants import _float as _F, _0_0s, _0_0, _0_001, _0_5, _0_75, _1_0
 from pygeodesy.datums import Datums, _earth_datum, _equall, _GDA2020_, _Names7, \
                             _negastr, Transform, _WGS84,  _EWGS84, _operator
 # from pygeodesy.ellipsoids import _EWGS84  # from .datums
@@ -93,7 +93,7 @@ from math import ceil as _ceil, fabs
 # import operator as _operator  # from .datums
 
 __all__ = _ALL_LAZY.trf
-__version__ = '25.05.12'
+__version__ = '25.09.11'
 
 _EP0CH    =  Epoch(0, low=0)
 _Es       = {_EP0CH: _EP0CH}  # L{Epoch}s, deleted below
@@ -1336,7 +1336,7 @@ _trfX(_ITRF2020_, _ITRF94_,   epoch=_E(2015),
                               xform=_X(   6.5,     -3.9,   -77.9,    3.98,     _0_0,      _0_0,      0.36),
                               rates=_R(   0.1,     -0.6,    -3.1,    0.12,     _0_0,      _0_0,      0.02))
 _trfX(_ITRF2020_, _ITRF93_,   epoch=_E(2015),
-                              xform=_X( -65.8,      1.9,   -71.3,    4.47,     -3.36,     -4.33,     0.75),
+                              xform=_X( -65.8,      1.9,   -71.3,    4.47,     -3.36,     -4.33,    _0_75),
                               rates=_R(  -2.8,     -0.2,    -2.3,    0.12,     -0.11,     -0.19,     0.07))
 _trfX(_ITRF2020_, _ITRF92_,   epoch=_E(2015),
                               xform=_X(  14.5,     -1.9,   -85.9,    3.27,     _0_0,      _0_0,      0.36),
@@ -1459,7 +1459,7 @@ _trfX(_ITRF2000_, _ITRF93_,   epoch=_E(1988),
                               xform=_X(  12.7,      6.5,   -20.9,    1.95,     -0.39,     0.8,      -1.14),
                               rates=_R(  -2.9,     -0.2,    -0.6,    0.01,     -0.11,    -0.19,      0.07))
 _trfX(_ITRF2000_, _ITRF92_,   epoch=_E(1988),
-                              xform=_X(   1.47,     1.35,   -1.39,   0.75,     _0_0,     _0_0,      -0.18),
+                              xform=_X(   1.47,     1.35,   -1.39,  _0_75,     _0_0,     _0_0,      -0.18),
                               rates=_R(  _0_0,     -0.06,   -0.14,   0.01,     _0_0,     _0_0,       0.02))
 _trfX(_ITRF2000_, _ITRF91_,   epoch=_E(1988),
                               xform=_X(   26.7,    27.5,   -19.9,    2.15,     _0_0,     _0_0,      -0.18),