pygeodesy 25.5.5__py2.py3-none-any.whl → 25.5.25__py2.py3-none-any.whl

pygeodesy/namedTuples.py

@@ -10,7 +10,7 @@ of C{_NamedTuple} defined in C{pygeodesy.named}.
 
 from pygeodesy.basics import isinstanceof, issubclassof, map1, _xinstanceof
 # from pygeodesy.cartesianBase import CartesianBase  # _MODS
-# from pygeodesy.constants import INT0  # from .units
+from pygeodesy.constants import INT0
 # from pygeodesy.dms import toDMS  # _MODS
 from pygeodesy.errors import _TypeError, _xattr, _xkwds, _xkwds_not
 from pygeodesy.interns import NN, _1_, _2_, _a_, _A_, _area_, _angle_, _b_, _B_, \
@@ -26,10 +26,10 @@ from pygeodesy.props import deprecated_property_RO, property_RO
 from pygeodesy.units import Band, Bearing, Degrees, Degrees2, Easting, FIx, \
                             Height, Int, Lam, Lat, Lon, Meter, Meter2, \
                             Northing, Number_, Phi, Precision_, Radians, \
-                            Radius, Scalar, Str,  INT0
+                            Radius, Scalar, Str
 
 __all__ = _ALL_LAZY.namedTuples
-__version__ = '24.11.22'
+__version__ = '25.05.12'
 
 # __DUNDER gets mangled in class
 _closest_     = 'closest'

pygeodesy/nvectorBase.py

@@ -16,7 +16,7 @@ from pygeodesy.constants import EPS, EPS0, EPS1, EPS_2, R_M, \
 # from pygeodesy.datums import _spherical_datum  # from .formy
 from pygeodesy.errors import IntersectionError, _ValueError, VectorError, \
                             _xattrs, _xkwds, _xkwds_pop2
-from pygeodesy.fmath import fdot, fidw, hypot  # PYCHOK fdot shared
+from pygeodesy.fmath import fidw, hypot
 from pygeodesy.fsums import Fsum, fsumf_
 from pygeodesy.formy import _isequalTo,  _spherical_datum
 # from pygeodesy.internals import _under  # from .named
@@ -38,7 +38,7 @@ from pygeodesy.vector3d import Vector3d, _xyzhdlln4
 from math import degrees, fabs, sqrt
 
 __all__ = _ALL_LAZY.nvectorBase
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 
 class NvectorBase(Vector3d):  # XXX kept private

pygeodesy/osgr.py

@@ -24,7 +24,7 @@ U{A Guide<https://www.OrdnanceSurvey.co.UK/documents/resources/guide-coordinate-
 and U{Ordnance Survey National Grid<https://WikiPedia.org/wiki/Ordnance_Survey_National_Grid>}.
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import halfs2, isbool, isfloat, map1, \
                             _splituple, _xsubclassof,  typename
@@ -54,7 +54,7 @@ from pygeodesy.utily import degrees90, degrees180, sincostan3, truncate
 from math import cos, fabs, radians, sin, sqrt
 
 __all__ = _ALL_LAZY.osgr
-__version__ = '25.04.12'
+__version__ = '25.05.12'
 
 _equivalent_ = 'equivalent'
 _OSGR_       = 'OSGR'

pygeodesy/points.py

@@ -52,8 +52,7 @@ from pygeodesy.named import classname, _NamedTuple, nameof, \
                             notImplemented, notOverloaded
 from pygeodesy.namedTuples import Bounds2Tuple, Bounds4Tuple, LatLon2Tuple, \
                                   NearestOn3Tuple, NearestOn5Tuple, \
-                                  Point3Tuple, Vector3Tuple, \
-                                  PhiLam2Tuple  # PYCHOK shared
+                                  Point3Tuple, Vector3Tuple
 from pygeodesy.props import Property_RO, property_doc_, property_RO
 from pygeodesy.streprs import Fmt, instr
 from pygeodesy.units import Number_, Radius, Scalar, Scalar_
@@ -63,7 +62,7 @@ from pygeodesy.utily import atan2b, degrees90, degrees180, degrees2m, \
 from math import cos, fabs, fmod as _fmod, radians, sin
 
 __all__ = _ALL_LAZY.points
-__version__ = '25.04.18'
+__version__ = '25.05.12'
 
 _ilat_  = 'ilat'
 _ilon_  = 'ilon'

pygeodesy/props.py

@@ -26,7 +26,7 @@ from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS, \
 from functools import wraps as _wraps
 
 __all__ = _ALL_LAZY.props
-__version__ = '25.04.14'
+__version__ = '25.05.21'
 
 _class_       = 'class'
 _DNL_         = _NL_ * 2  # PYCHOK used!
@@ -460,7 +460,7 @@ def _deprecated(call, kind, qual_d):
        @see: Brett Slatkin, "Effective Python", 2019 page 105, 2nd
              ed, Addison-Wesley.
     '''
-    doc = _DOCof(call)
+    doc = _DEPRECATEDof(call)
 
     @_wraps(call)  # PYCHOK self?
     def _deprecated_call(*args, **kwds):
@@ -517,6 +517,17 @@ def _deprecated_module(name):  # PYCHOK no cover
         _throwarning(_module_, name, _dont_use_)
 
 
+def _DEPRECATEDof(obj):
+    '''(INTERNAL) Get uniform DEPRECATED __doc__ string.
+    '''
+    try:
+        d = obj.__doc__.strip()
+        i = d.find(_DEPRECATED_)
+    except AttributeError:
+        i = -1
+    return _DOT_(_DEPRECATED_, NN) if i < 0 else d[i:]
+
+
 if _WARNINGS_X_DEV:
     class deprecated_property(_PropertyBase):
         '''Decorator for a DEPRECATED C{property} or C{Property}.
@@ -524,7 +535,7 @@ if _WARNINGS_X_DEV:
         def __init__(self, method):
             '''Decorator for a DEPRECATED C{property} or C{Property} getter.
             '''
-            doc = _DOCof(method)
+            doc = _DEPRECATEDof(method)
 
             def _fget(inst):  # PYCHOK no cover
                 '''Get the C{property} or C{Property} value.
@@ -554,7 +565,7 @@ if _WARNINGS_X_DEV:
                     '''Set the C{property} or C{Property} value.
                     '''
                     q = _qualified(inst, self.name)
-                    _throwarning(typename(property), q, _DOCof(method))
+                    _throwarning(typename(property), q, _DEPRECATEDof(method))
                     method(inst, val)
                     # self._update(inst)  # un-cache this item
 
@@ -594,7 +605,7 @@ def deprecated_property_RO(method):
 def _deprecated_RO(method, _RO):
     '''(INTERNAL) Create a DEPRECATED C{property_RO} or C{Property_RO}.
     '''
-    doc = _DOCof(method)
+    doc = _DEPRECATEDof(method)
 
     if _WARNINGS_X_DEV:
 
@@ -614,17 +625,6 @@ def _deprecated_RO(method, _RO):
         return _RO(method, doc=doc)
 
 
-def _DOCof(obj):
-    '''(INTERNAL) Get uniform DEPRECATED __doc__ string.
-    '''
-    try:
-        d = obj.__doc__.strip()
-        i = d.find(_DEPRECATED_)
-    except AttributeError:
-        i = -1
-    return _DOT_(_DEPRECATED_, NN) if i < 0 else d[i:]
-
-
 def _qualified(inst, name):
     '''(INTERNAL) Fully qualify a name.
     '''

pygeodesy/resections.py

@@ -10,7 +10,7 @@ L{triAngle}, L{triAngle5}, L{triSide}, L{triSide2} and L{triSide4}.
        U{Pierlot<http://www.Telecom.ULg.ac.BE/publi/publications/pierlot/Pierlot2014ANewThree>}.
 '''
 # make sure int/int division yields float quotient
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import map1, map2, _zip,  _ALL_LAZY, typename
 from pygeodesy.constants import EPS, EPS0, EPS02, INT0, PI, PI2, PI_2, PI_4, \
@@ -34,7 +34,7 @@ from pygeodesy.vector3d import _otherV3d, Vector3d
 from math import cos, degrees, fabs, radians, sin, sqrt
 
 __all__ = _ALL_LAZY.resections
-__version__ = '25.05.04'
+__version__ = '25.05.12'
 
 _concyclic_ = 'concyclic'
 _PA_        = 'PA'
@@ -798,12 +798,12 @@ def triArea(a, b, c):
     try:
         r, y, x = sorted(map1(float, a, b, c))
         if r > 0:  # r = min(a, b, c)
-            z = r
-            d = x - y
-            r = (z + d) * (z - d)
+            s =  r
+            d =  x - y
+            r = (s + d) * (s - d)
             if r:
-                x += y
-                r *= (x + z) * (x - z)
+                y +=  x
+                r *= (y + s) * (y - s)
         if r < 0:
             raise ValueError(_negative_)
         return sqrt(r / _16_0) if r else _0_0
@@ -927,9 +927,11 @@ def triSide4(radA, radB, c):
         sc = fsum1f_(sa * cb, sb * ca)
         if sc < EPS0 or min(sa, sb) < 0:
             raise ValueError(_invalid_)
-        sc = c / sc
-        return TriSide4Tuple((sa * sc), (sb * sc), rC, (sa * sb * sc),
-                             name=typename(triSide4))
+        sc  = c / sc
+        sa *= sc
+        sd  = sa * sb
+        sb *= sc
+        return TriSide4Tuple(sa, sb, rC, sd, name=typename(triSide4))
 
     except (TypeError, ValueError) as x:
         raise TriangleError(radA=radA, radB=radB, c=c, cause=x)
@@ -948,7 +950,8 @@ def wildberger3(a, b, c, alpha, beta, R3=min):
        @arg alpha: Angle subtended by triangle side B{C{b}} (C{degrees}, non-negative).
        @arg beta: Angle subtended by triangle side B{C{a}} (C{degrees}, non-negative).
        @kwarg R3: Callable to determine C{R3} from C{(R3 - C)**2 = D}, typically standard
-                  Python function C{min} or C{max}, invoked with 2 arguments.
+                  Python function C{min} or C{max}, invoked with 2 arguments, each
+                  C{float} or an L{Fsum<pygeodesy.fsums.Fsum>}.
 
        @return: L{Survey3Tuple}C{(PA, PB, PC)} with distance from survey point C{P} to
                 each of the triangle corners C{A}, C{B} and C{C}, same units as B{C{a}},
@@ -962,19 +965,20 @@ def wildberger3(a, b, c, alpha, beta, R3=min):
              U{Devine Proportions, page 252<http://www.MS.LT/derlius/WildbergerDivineProportions.pdf>}
              and function L{snellius3}.
     '''
-    def _s(x):
+    def _sin2(x):
         return sin(x)**2
 
     def _vpa(r3, q2, q3, s2, s3):
+        r4 = r3 * _4_0
+        s4 = r4 * s2
+        if isnear0(s4):
+            raise ValueError(_coincident_)
         r1 =  s2 * q3 / s3
-        r  =  r1 * r3 * _4_0
-        R  = _Fsumf_(r1, r3, -q2)
-        R *=  R  # -(R**2 ...
-        R -=  r  # ... - r) / s3
-        n  = -R.fover(s3)
-        if n < 0 or r < EPS0:
+        N2 = _Fsumf_(r1, r3, -q2)**2
+        n  = (r4 * r1 - N2).fover(s4)
+        if n < 0:
             raise ValueError(_coincident_)
-        return sqrt((n / r) * q3) if n else _0_0
+        return sqrt(n) if n else _0_0
 
     try:
         a, b, c, da, db = _noneg(a, b, c, alpha, beta)
@@ -983,21 +987,23 @@ def wildberger3(a, b, c, alpha, beta, R3=min):
             raise ValueError(_coincident_)
 
         ra, rb = map1(radians, da, db)
-        s1, s2, s3 = s = map1(_s, rb, ra, ra + rb)  # rb, ra!
+        s1, s2, s3 = s = map1(_sin2, rb, ra, ra + rb)  # rb, ra!
         if min(s) < EPS02:
             raise ValueError(_or(_coincident_, _colinear_))
+        s += _Fsumf_(*s),  # == fsum1(s),
 
         Q  = _Fsumf_(*q)   # == a**2 + b**2 + ...
-        s += _Fsumf_(*s),  # == fsum1(s),
         C0 =  Fdot(s, q1, q2, q3, -Q * _0_5)
         r3 =  C0.fover(-s3)  # C0 /= -s3
+
         Q *=  Q  # Q**2 - 2 * (a**4 + b**4 ...
-        Q -=  hypot2_(*q) *_2_0  # ... + c**4)
-        d0 =  Q.fmul(s1 * s2).fover(s3)
+        Q -=  hypot2_(*q) * _2_0  # ... + c**4)
+        Q *=  s1 * s2  # Q * s1 * s2 / s3
+        d0 =  Q.fover(s3)
         if d0 > EPS02:  # > c0
             _xcallable(R3=R3)
             d0 = sqrt(d0)
-            r3 = R3(float(C0 + d0), float(C0 - d0))  # XXX min or max
+            r3 = R3(C0 + d0, C0 - d0).as_iscalar  # XXX min or max
         elif d0 < (-EPS02):
             raise ValueError(_negative_)
 

pygeodesy/rhumb/aux_.py

@@ -26,7 +26,7 @@ License.  For more information, see the U{GeographicLib<https://GeographicLib.So
        windows/win32/fileio/naming-a-file#naming-conventions>} with and without extension.
 '''
 # make sure int/int division yields float quotient
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.auxilats.auxAngle import AuxMu, AuxPhi,  hypot
 from pygeodesy.auxilats.auxDLat import AuxDLat, _DClenshaw
@@ -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__ = '25.04.12'
+__version__ = '25.05.12'
 
 # DIGITS = (sizeof(real) * 8) bits
 #        = (ctypes.sizeof(ctypes.c_double(1.0)) * 8) bits

pygeodesy/rhumb/bases.py

@@ -21,7 +21,7 @@ Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2014-2024) and lice
 License.  For more information, see the U{GeographicLib<https://GeographicLib.SourceForge.io>} documentation.
 '''
 # make sure int/int division yields float quotient
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import _copysign, itemsorted, unsigned0, _xinstanceof
 from pygeodesy.constants import EPS, EPS0, EPS1, INT0, NAN, _over, \
@@ -44,7 +44,7 @@ from pygeodesy.namedTuples import Distance2Tuple, LatLon2Tuple
 from pygeodesy.props import deprecated_method, Property, Property_RO, \
                             property_RO, _update_all
 from pygeodesy.streprs import Fmt, pairs
-from pygeodesy.units import Float_, Lat, Lon, Meter, Radius_,  Int  # PYCHOK shared
+from pygeodesy.units import Float_, Lat, Lon, Meter, Radius_
 from pygeodesy.utily import acos1, _azireversed, _loneg, sincos2d, sincos2d_, \
                            _unrollon, _Wrap
 from pygeodesy.vector3d import _intersect3d3, Vector3d  # in .Intersection below
@@ -52,7 +52,7 @@ from pygeodesy.vector3d import _intersect3d3, Vector3d  # in .Intersection below
 from math import cos, fabs
 
 __all__ = ()
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _anti_ = _Dash('anti')
 _rls   = []  # instances of C{RbumbLine...} to be updated

pygeodesy/rhumb/ekx.py

@@ -20,7 +20,7 @@ Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2014-2024) and lice
 License.  For more information, see the U{GeographicLib<https://GeographicLib.SourceForge.io>} documentation.
 '''
 # make sure int/int division yields float quotient
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import copysign0, neg
 from pygeodesy.constants import PI_2, _0_0s, _0_0, _0_5, _1_0, \
@@ -31,8 +31,7 @@ from pygeodesy.errors import RhumbError, _xkwds_pop2, _Xorder
 from pygeodesy.fmath import hypot, hypot1
 # from pygeodesy.fsums import fsum1f_  # _MODS
 # from pygeodesy.karney import Caps  # from .rhumb.bases
-from pygeodesy.ktm import KTransverseMercator, _Xs, \
-                         _AlpCoeffs, _BetCoeffs  # PYCHOK used!
+from pygeodesy.ktm import _Xs, _AlpCoeffs, _BetCoeffs  # PYCHOK used!
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS
 from pygeodesy.props import deprecated_method, Property, Property_RO, \
                             property_RO
@@ -43,7 +42,7 @@ from pygeodesy.utily import atan1, sincos2_
 from math import asinh, atan, cos, cosh, radians, sin, sinh, sqrt, tan  # as _tan
 
 __all__ = _ALL_LAZY.rhumb_ekx
-__version__ = '25.04.12'
+__version__ = '25.05.12'
 
 
 class Rhumb(RhumbBase):

pygeodesy/rhumb/solve.py

@@ -9,7 +9,7 @@ as an (exact) rhumb or rhumb line from I{either GeographicLib 2.0 or 2.2+}.
        path of the C{RhumbSolve} executable.
 '''
 from pygeodesy.basics import _xinstanceof
-from pygeodesy.constants import _0_0, _180_0, _N_180_0, _over,  _90_0  # PYCHOK used!
+from pygeodesy.constants import _180_0, _N_180_0, _over,  _0_0  # PYCHOK used!
 from pygeodesy.errors import RhumbError   # PYCHOK used!
 from pygeodesy.interns import NN, _a12_, _azi12_, _DMAIN_, _lat2_, _lon2_, _s12_, \
                              _S12_, _UNDER_
@@ -21,7 +21,7 @@ from pygeodesy.solveBase import _SolveGDictBase, _SolveGDictLineBase
 from pygeodesy.utily import _unrollon, _Wrap, wrap360
 
 __all__ = _ALL_LAZY.rhumb_solve
-__version__ = '25.04.12'
+__version__ = '25.05.07'
 
 
 class _RhumbSolveBase(_SolveGDictBase):

pygeodesy/simplify.py

@@ -65,7 +65,7 @@ See:
  - U{https://PyPI.org/project/simplification}
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 # from pygeodesy.basics import len2  # from .fmath
 from pygeodesy.constants import EPS, R_M, _1_0
@@ -80,7 +80,7 @@ from pygeodesy.units import _ALL_LAZY, _1mm, Radius_
 from math import degrees, fabs, radians
 
 __all__ = _ALL_LAZY.simplify
-__version__ = '24.12.02'
+__version__ = '25.05.12'
 
 
 # try:

pygeodesy/solveBase.py

@@ -20,10 +20,10 @@ from pygeodesy.named import callername, _name2__, notOverloaded
 from pygeodesy.props import Property, Property_RO, property_RO, _update_all
 from pygeodesy.streprs import Fmt, fstr, fstrzs, pairs, strs
 from pygeodesy.units import Precision_
-from pygeodesy.utily import unroll180,  wrap360  # PYCHOK shared
+from pygeodesy.utily import unroll180
 
 __all__ = _ALL_LAZY.solveBase
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _ERROR_ = 'ERROR'
 

pygeodesy/sphericalBase.py

@@ -10,10 +10,10 @@ and published under the same MIT Licence**, see
 U{Latitude/Longitude<https://www.Movable-Type.co.UK/scripts/latlong.html>}.
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import _copysign, isbool, _isin, isinstanceof, map1
-from pygeodesy.cartesianBase import CartesianBase,  Bearing2Tuple
+from pygeodesy.cartesianBase import CartesianBase
 from pygeodesy.constants import EPS, EPS0, PI, PI2, PI_2, R_M, \
                                _0_0, _0_5, _1_0, _180_0, _360_0, \
                                _over, isnear0, isnon0
@@ -24,9 +24,9 @@ from pygeodesy.fmath import favg, fdot, hypot, sqrt_a
 from pygeodesy.interns import _COMMA_, _concentric_, _datum_, _distant_, \
                               _exceed_PI_radians_, _name_, _near_, \
                               _radius_, _too_
-from pygeodesy.latlonBase import LatLonBase,  _trilaterate5  # PYCHOK passed
+from pygeodesy.latlonBase import LatLonBase
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS
-# from pygeodesy.namedTuples import Bearing2Tuple  # from .cartesianBase
+# from pygeodesy.namedTuples import Bearing2Tuple  # _MODS
 from pygeodesy.nvectorBase import NvectorBase,  Fmt
 from pygeodesy.props import deprecated_method, property_doc_, property_RO, \
                            _update_all
@@ -40,7 +40,7 @@ from pygeodesy.utily import acos1, asin1, atan2b, atan2d, degrees90, \
 from math import cos, fabs, log, sin, sqrt
 
 __all__ = _ALL_LAZY.sphericalBase
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 
 class CartesianSphericalBase(CartesianBase):
@@ -157,7 +157,7 @@ class LatLonSphericalBase(LatLonBase):
         # .initialBearingTo is inside .-Nvector and .-Trigonometry
         i = self.initialBearingTo(other, wrap=wrap, raiser=raiser)  # PYCHOK .initialBearingTo
         f = self.finalBearingTo(  other, wrap=wrap, raiser=raiser)
-        return Bearing2Tuple(i, f, name=self.name)
+        return _MODS.namedTuples.Bearing2Tuple(i, f, name=self.name)
 
     @property_doc_(''' this point's datum (L{Datum}).''')
     def datum(self):

pygeodesy/sphericalNvector.py

@@ -31,7 +31,7 @@ Also note that on a spherical earth model, an n-vector is equivalent
 to a normalised version of an (ECEF) cartesian coordinate.
 '''
 # make sure int/int division yields float quosient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import _isin, _xinstanceof,  typename
 from pygeodesy.constants import EPS, EPS0, PI, PI2, PI_2, R_M, \
@@ -62,7 +62,7 @@ from pygeodesy.utily import atan2, degrees360, sincos2, sincos2_, sincos2d, \
 # from math import fabs  # from utily
 
 __all__ = _ALL_LAZY.sphericalNvector
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _lines_ = 'lines'
 

pygeodesy/sphericalTrigonometry.py

@@ -14,7 +14,7 @@ I{(C) Chris Veness 2011-2024} published under the same MIT Licence**, see
 U{Latitude/Longitude<https://www.Movable-Type.co.UK/scripts/latlong.html>}.
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import copysign0, _isin, map1, signOf,  typename
 from pygeodesy.constants import EPS, EPS1, EPS4, PI, PI2, PI_2, PI_4, R_M, \
@@ -34,8 +34,9 @@ from pygeodesy.interns import _1_, _2_, _coincident_, _composite_, _colinear_, \
                               _concentric_, _convex_, _end_, _infinite_, \
                               _invalid_, _line_, _near_, _null_, _parallel_, \
                               _point_, _SPACE_, _too_
+from pygeodesy.latlonBase import _trilaterate5
 from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS, _ALL_OTHER
-# from pygeodesy.nvectorBase import NvectorBase, sumOf  # _MODE
+# from pygeodesy.nvectorBase import NvectorBase, sumOf  # _MODS
 from pygeodesy.namedTuples import LatLon2Tuple, LatLon3Tuple, NearestOn3Tuple, \
                                   Triangle7Tuple, Triangle8Tuple
 from pygeodesy.points import ispolar, nearestOn5 as _nearestOn5, \
@@ -43,10 +44,10 @@ from pygeodesy.points import ispolar, nearestOn5 as _nearestOn5, \
 from pygeodesy.props import deprecated_function, deprecated_method
 from pygeodesy.sphericalBase import _m2radians, CartesianSphericalBase, \
                                     _intersecant2, LatLonSphericalBase, \
-                                    _rads3, _radians2m, _trilaterate5
+                                    _rads3, _radians2m
 # from pygeodesy.streprs import Fmt as _Fmt  # from .points XXX shadowed
 from pygeodesy.units import Bearing_, Height, _isDegrees, _isRadius, Lamd, \
-                              Phid, Radius_, Scalar
+                            Phid, Radius_, Scalar
 from pygeodesy.utily import acos1, asin1, atan1d, atan2, atan2d, degrees90, \
                             degrees180, degrees2m, m2radians, radiansPI2, \
                             sincos2_, tan_2, unrollPI, _unrollon, _unrollon3, \
@@ -56,7 +57,7 @@ from pygeodesy.vector3d import sumOf, Vector3d
 from math import asin, cos, degrees, fabs, radians, sin
 
 __all__ = _ALL_LAZY.sphericalTrigonometry
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _PI_EPS4 = PI - EPS4
 if _PI_EPS4 >= PI:

pygeodesy/trf.py

@@ -93,7 +93,7 @@ from math import ceil as _ceil, fabs
 # import operator as _operator  # from .datums
 
 __all__ = _ALL_LAZY.trf
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _EP0CH    =  Epoch(0, low=0)
 _Es       = {_EP0CH: _EP0CH}  # L{Epoch}s, deleted below
@@ -148,7 +148,7 @@ def _E(epoch):  # deleted below
     '''(INTERNAL) De-dup L{Epochs}s.
     '''
     e = Epoch(_F(epoch))
-    return _Es.setdefault(e, e)  # PYCHOK del
+    return _Es.setdefault(e, e)  # noqa: F821 del
 
 
 def _P(ps, name, _Ps):  # deleted below
@@ -1299,11 +1299,11 @@ def _trfX(n1, n2, raiser=True, **epoch_xform_rates):
 
 
 def _X(*ps):  # deleted below
-    return _P(ps, _xform_, _Xs)  # PYCHOK del
+    return _P(ps, _xform_, _Xs)  # noqa: F821 del
 
 
 def _R(*ps):  # deleted below
-    return _P(ps, _rates_, _Rs)  # PYCHOK del
+    return _P(ps, _rates_, _Rs)  # noqa: F821 del
 
 
 _P_0_0s = TRFXform7Tuple(_0_0s(len(_Names7)), name='unity')

pygeodesy/triaxials.py

@@ -28,7 +28,7 @@ see the U{GeographicLib<https://GeographicLib.SourceForge.io>} documentation.
 @var Triaxials.WGS84_35: Triaxial(name='WGS84_35', a=6378172, b=6378102, c=6356752.314245179, e2ab=0.00002195, e2bc=0.006683478, e2ac=0.006705281, volume=1083207319768789942272, area=510065621722018.125, area_p=510065626587483.3125)
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import _isin, isLatLon, isscalar
 from pygeodesy.constants import EPS, EPS0, EPS02, EPS4, INT0, PI2, PI_3, PI4, \
@@ -59,7 +59,7 @@ from pygeodesy.vector3d import _otherV3d, Vector3d,  _ALL_LAZY, _MODS
 from math import fabs, sqrt
 
 __all__ = _ALL_LAZY.triaxials
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _not_ordered_ = _not_('ordered')
 _omega_       = 'omega'

pygeodesy/units.py

@@ -7,18 +7,17 @@ L{Feet}, L{Meter}, L{Radians}, etc.
 '''
 
 from pygeodesy.basics import isscalar, issubclassof, signOf,  typename
-from pygeodesy.constants import EPS, EPS1, PI, PI2, PI_2, _umod_360, _0_0, \
-                               _0_001,  _0_5, INT0  # PYCHOK for .mgrs, .namedTuples
+from pygeodesy.constants import EPS, EPS1, PI, PI2, PI_2, _umod_360, _0_0, _0_001
 from pygeodesy.dms import F__F, F__F_, S_NUL, S_SEP, parseDMS, parseRad, _toDMS
 from pygeodesy.errors import _AssertionError, TRFError, UnitError, _xattr, _xcallable
 # from pygeodesy.internals import typename  # from .basics
 from pygeodesy.interns import NN, _azimuth_, _band_, _bearing_, _COMMASPACE_, \
                              _degrees_, _degrees2_, _distance_, _E_, _easting_, \
-                             _epoch_, _EW_, _feet_, _height_, _lam_, _lat_, _LatLon_, \
-                             _lon_, _meter_, _meter2_, _N_, _negative_, _northing_, \
+                             _epoch_, _EW_, _feet_, _height_, _lam_, _lat_, _lon_, \
+                             _meter_, _meter2_, _N_, _negative_, _northing_, _radians_, \
                              _NS_, _NSEW_, _number_, _PERCENT_, _phi_, _precision_, \
-                             _radians_, _radians2_, _radius_, _S_, _scalar_, \
-                             _units_, _W_, _zone_,  _std_  # PYCHOK used!
+                             _radians2_, _radius_, _S_, _scalar_, \
+                             _W_, _zone_,  _std_  # PYCHOK used!
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS
 # from pygeodesy.named import _name__  # _MODS
 from pygeodesy.props import Property_RO
@@ -28,7 +27,7 @@ from pygeodesy.unitsBase import Float, Int, _NamedUnit, Radius, Str,  Fmt, fstr
 from math import degrees, isnan, radians
 
 __all__ = _ALL_LAZY.units
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 
 class Float_(Float):
@@ -861,7 +860,7 @@ def _std_repr(*Classes):
     '''
     from pygeodesy.internals import _getenv
     for C in Classes:
-        if hasattr(C, typename(_std_repr)):  # PYCHOK del _std_repr
+        if hasattr(C, typename(_std_repr)):  # noqa: F821 del
             env = 'PYGEODESY_%s_STD_REPR' % (typename(C).upper(),)
             if _getenv(env, _std_).lower() != _std_:
                 C._std_repr = False

pygeodesy/utily.py

@@ -9,7 +9,7 @@ and U{Vector-based geodesy<https://www.Movable-Type.co.UK/scripts/latlong-vector
 and published under the same MIT Licence**.
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import _copysign, isinstanceof, isint, isstr
 from pygeodesy.constants import EPS, EPS0, INF, NAN, PI, PI2, PI_2, R_M, \
@@ -28,7 +28,7 @@ from math import acos, asin, atan2 as _atan2, cos, degrees, fabs, radians, \
                  sin, tan as _tan  # pow
 
 __all__ = _ALL_LAZY.utily
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _G_DEG     = _float(  400.0 / _360_0)  # grades per degree
 _G_RAD     = _float(  400.0 /  PI2)    # grades per radian

pygeodesy/utmupsBase.py

@@ -5,7 +5,7 @@ u'''(INTERNAL) Private class C{UtmUpsBase}, functions and constants
 for modules L{epsg}, L{etm}, L{mgrs}, L{ups} and L{utm}.
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import _copysign, _isin, isint, isscalar, isstr, \
                               neg_, _xinstanceof, _xsubclassof
@@ -32,7 +32,7 @@ from pygeodesy.utily import atan1, _Wrap, wrap360
 from math import cos, degrees, fabs, sin, tan  # copysign as _copysign
 
 __all__ = _ALL_LAZY.utmupsBase
-__version__ = '25.04.26'
+__version__ = '25.05.12'
 
 _UPS_BANDS = _A_, _B_, _Y_, _Z_  # UPS polar bands SE, SW, NE, NW
 # _UTM_BANDS = _MODS.utm._Bands