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

pygeodesy/latlonBase.py

@@ -25,9 +25,8 @@ from pygeodesy.errors import _AttributeError, IntersectionError, \
 # from pygeodesy.fmath import favg  # _MODS
 # from pygeodesy import formy as _formy  # _MODS.into
 from pygeodesy.internals import _passarg, typename
-from pygeodesy.interns import NN, _COMMASPACE_, _concentric_, _height_, \
-                             _intersection_, _LatLon_, _m_, _negative_, \
-                             _no_, _overlap_, _too_,  _point_  # PYCHOK used!
+from pygeodesy.interns import NN, _COMMASPACE_, _concentric_, _intersection_, \
+                             _LatLon_, _m_, _no_, _overlap_,  _point_  # PYCHOK used!
 # from pygeodesy.iters import PointsIter, points2  # _MODS
 # from pygeodesy.karney import Caps  # _MODS
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS
@@ -49,7 +48,7 @@ from contextlib import contextmanager
 from math import asin, cos, degrees, fabs, radians
 
 __all__ = _ALL_LAZY.latlonBase
-__version__ = '25.04.28'
+__version__ = '25.05.07'
 
 _formy = _MODS.into(formy=__name__)
 

pygeodesy/lazily.py

@@ -510,7 +510,7 @@ class _ALL_MODS(_internals._MODS_Base):
 _internals._MODS = _ALL_MODS = _ALL_MODS()  # PYCHOK singleton
 
 __all__ = _ALL_LAZY.lazily
-__version__ = '25.04.30'
+__version__ = '25.05.19'
 
 
 def _ALL_OTHER(*objs):

pygeodesy/lcc.py

@@ -22,7 +22,7 @@ and John P. Snyder U{'Map Projections - A Working Manual'<https://Pubs.USGS.gov/
 @var Conics.WRF_Lb: Conic(name='WRF_Lb', lat0=40, lon0=-97, par1=33, par2=45, E0=0, N0=0, k0=1, SP=2, datum=Datum(name='WGS84', ellipsoid=Ellipsoids.WGS84, transform=Transforms.WGS84)
 '''
 # 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, _xinstanceof, _xsubclassof, \
                              typename
@@ -50,7 +50,7 @@ from pygeodesy.utily import atan1, degrees90, degrees180, sincos2, tanPI_2_2
 from math import atan, fabs, log, radians, sin, sqrt
 
 __all__ = _ALL_LAZY.lcc
-__version__ = '25.04.14'
+__version__ = '25.05.26'
 
 _E0_   = 'E0'
 _N0_   = 'N0'
@@ -457,7 +457,7 @@ class Lcc(_NamedBase):
 #                                  h=self.height, conic=self.conic,
 #                                  name=self._name__(name))
 #       args, kwds = _args_kwds(**kwds)
-#       return self.__class__(*args, **kwds)  # .classof
+#       return type(self)(*args, **kwds)  # .classof
 
     @Property_RO
     def easting(self):

pygeodesy/ltp.py

@@ -11,7 +11,7 @@ and L{ChLVe} and L{Ltp}, L{ChLV}, L{LocalError}, L{Attitude} and L{Frustum}.
       <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1LocalCartesian.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 _args_kwds_names, _isin, map1, map2, _xinstanceof, \
                              _xsubclassof,  typename   # .datums
@@ -45,7 +45,7 @@ from pygeodesy.vector3d import _ALL_LAZY, Vector3d
 # from math import fabs, floor as _floor  # from .fmath, .fsums
 
 __all__ = _ALL_LAZY.ltp
-__version__ = '25.05.01'
+__version__ = '25.05.12'
 
 _height0_ = _height_ + _0_
 _narrow_  = 'narrow'
@@ -306,11 +306,11 @@ class Frustum(_NamedBase):
             # left and right corners, or swapped
             if r < EPS:  # no roll
                 r =  a * tan_h_2
-                l = -r  # PYCHOK l is ell
+                l = -r  # noqa: E741 l is ell
             else:  # roll
-                r, l = tand_(r - h_2, r + h_2, roll_hfov=r)  # PYCHOK l is ell
+                r, l = tand_(r - h_2, r + h_2, roll_hfov=r)  # noqa: E741 l is ell
                 r *= -a  # negate right positive
-                l *= -a  # PYCHOK l is ell
+                l *= -a  # noqa: E741 l is ell
             y = a * cotd(e, tilt_vfov=e)
             return (l, y), (r, y)
 

pygeodesy/mgrs.py

@@ -36,7 +36,7 @@ and compare the MGRS results with those from I{Karney}'s utility U{GeoConvert
 '''
 
 from pygeodesy.basics import halfs2, _isin, _splituple, _xinstanceof
-# from pygeodesy.constants import _0_5  # from .units
+from pygeodesy.constants import _0_5
 from pygeodesy.datums import _ellipsoidal_datum, _WGS84
 from pygeodesy.errors import _AssertionError, MGRSError, _parseX, \
                              _ValueError, _xkwds
@@ -49,14 +49,14 @@ from pygeodesy.named import _name2__, _NamedBase, _NamedTuple, _Pass
 from pygeodesy.namedTuples import EasNor2Tuple, UtmUps5Tuple
 from pygeodesy.props import deprecated_property_RO, property_RO, Property_RO
 from pygeodesy.streprs import enstr2, _enstr2m3, Fmt, _resolution10, _xzipairs
-from pygeodesy.units import Easting, Northing, Str, _100km,  _0_5
+from pygeodesy.units import Easting, Northing, Str, _100km
 from pygeodesy.units import _1um, _2000km  # PYCHOK used!
 from pygeodesy.ups import _hemi, toUps8, Ups, _UPS_ZONE
 from pygeodesy.utm import toUtm8, _to3zBlat, Utm, _UTM_ZONE_MAX, _UTM_ZONE_MIN
 # from pygeodesy.utmupsBase import _UTM_ZONE_MAX, _UTM_ZONE_MIN  # from .utm
 
 __all__ = _ALL_LAZY.mgrs
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _AN_    = 'AN'  # default south pole grid tile and band B
 _AtoPx_ = _AtoZnoIO_.tillP

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.26'
 
 _class_       = 'class'
 _DNL_         = _NL_ * 2  # PYCHOK used!
@@ -347,7 +347,7 @@ class property_RO(_PropertyBase):
             if Clas:  # overrides inst.__class__
                 d = Clas[0].__dict__.get(uname, MISSING)
             else:
-                d = getattr(inst.__class__, uname, MISSING)
+                d = getattr(type(inst), uname, MISSING)
 #               if d is MISSING:  # XXX superfluous
 #                   for c in type(inst).__mro__[:-1]:
 #                       if uname in c.__dict__:
@@ -403,7 +403,7 @@ class property_ROver(_property_RO___):
         '''
         v = self.method(inst)
         n = self.name
-        C = inst.__class__
+        C = type(inst)
         for c in C.__mro__:  # [:-1]
             if getattr(c, n, None) is self:
                 setattr(c, n, v)  # overwrite property_ROver
@@ -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.26'
 
 
 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):
@@ -592,7 +592,7 @@ def _intersecant2(c, r, p, b, radius=R_M, exact=False, height=None, wrap=False):
         p = _unrollon(c, p, wrap=wrap)
     nonexact = exact is None
 
-    if not isinstanceof(r, c.__class__, p.__class__):
+    if not isinstanceof(r, type(c), type(p)):
         r = Radius_(circle=r)
     elif nonexact:
         r = c.distanceTo(r, radius=radius, wrap=wrap)
@@ -601,7 +601,7 @@ def _intersecant2(c, r, p, b, radius=R_M, exact=False, height=None, wrap=False):
     else:
         raise _ValueError(exact=exact)
 
-    if not isinstanceof(b, c.__class__, p.__class__):
+    if not isinstanceof(b, type(c), type(p)):
         b = Bearing(b)
     elif nonexact:
         b = p.initialBearingTo(b, wrap=wrap)