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

pygeodesy/azimuthal.py

@@ -40,7 +40,7 @@ altitude in Earth radii<https://WikiPedia.org/wiki/Azimuthal_equidistant_project
 #/media/File:Comparison_azimuthal_projections.svg>}.
 '''
 # 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, _xinstanceof  # from .ellipsoidalBase
 from pygeodesy.constants import EPS, EPS0, EPS1, NAN, isnon0, _umod_360, \
@@ -71,7 +71,7 @@ from pygeodesy.utily import asin1, atan1, atan2, atan2b, atan2d, \
 from math import acos, degrees, fabs, sin, sqrt
 
 __all__ = _ALL_LAZY.azimuthal
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _EPS_K         = _EPStol * _0_1  # Karney's eps_ or _EPSmin * _0_1?
 _over_horizon_ = 'over horizon'

pygeodesy/booleans.py

@@ -15,7 +15,7 @@ C{reverse-difference}, C{sum} and C{union}.
       <http://www.Inf.USI.CH/hormann/papers/Greiner.1998.ECO.pdf>}.
 '''
 # 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, isodd, issubclassof, map2, \
                              _xscalar,  typename
@@ -45,7 +45,7 @@ from pygeodesy.utily import fabs, _unrollon, _Wrap
 # from math import fabs  # from .utily
 
 __all__ = _ALL_LAZY.booleans
-__version__ = '25.04.30'
+__version__ = '25.05.12'
 
 _0EPS  =  EPS  # near-zero, positive
 _EPS0  = -EPS  # near-zero, negative

pygeodesy/cartesianBase.py

@@ -28,7 +28,7 @@ from pygeodesy.interns import _ellipsoidal_, _spherical_  # PYCHOK used!
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS
 from pygeodesy.named import _name2__, _Pass
 from pygeodesy.namedTuples import LatLon4Tuple, _NamedTupleTo , Vector3Tuple, \
-                                  Vector4Tuple,  Bearing2Tuple  # PYCHOK .sphericalBase
+                                  Vector4Tuple
 # from pygeodesy.nvectorBase import _N_vector  # _MODS
 from pygeodesy.props import deprecated_method, Property, Property_RO, property_doc_, \
                             property_RO, _update_all
@@ -44,7 +44,7 @@ from pygeodesy.vector3d import Vector3d, _xyzhdlln4
 # from math import degrees, fabs, radians, sqrt  # from .fmath, .utily
 
 __all__ = _ALL_LAZY.cartesianBase
-__version__ = '25.04.28'
+__version__ = '25.05.07'
 
 _r_         = 'r'
 _resections = _MODS.into(resections=__name__)

pygeodesy/clipy.py

@@ -10,7 +10,7 @@ I{Greiner-Hormann} and L{clipSH} and L{clipSH3} I{Sutherland-Hodgeman}.
 .
 '''
 # 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,  typename
 from pygeodesy.constants import EPS, _0_0, _1_0
@@ -31,7 +31,7 @@ from pygeodesy.units import Bool, FIx, HeightX, Lat, Lon, Number_
 # from math import fabs  # from .fmath
 
 __all__ = _ALL_LAZY.clipy
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _fj_       = 'fj'
 _original_ = 'original'

pygeodesy/constants.py

@@ -8,7 +8,7 @@ L{pygeodesy.isnear1}, L{pygeodesy.isnear90}, L{pygeodesy.isneg0}, L{pygeodesy.is
 L{pygeodesy.isnon0} and L{pygeodesy.remainder}.
 '''
 # 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, iscomplex, isint
 from pygeodesy.errors import _xError, _xError2, _xkwds_get1, _xkwds_item2
@@ -26,7 +26,7 @@ except ImportError:  # Python 2-
     _inf, _nan = float(_INF_), float(_NAN_)
 
 __all__ = _ALL_LAZY.constants
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 
 def _copysign_0_0(y):

pygeodesy/datums.py

@@ -65,7 +65,7 @@ datum, q.v. U{"A Guide to Coordinate Systems in Great Britain", Section 6
 @var Transforms.WGS84: Transform(name='WGS84', tx=0.0, ty=0.0, tz=0.0, s1=1.0, rx=0.0, ry=0.0, rz=0.0, s=0.0, sx=0.0, sy=0.0, sz=0.0)
 '''
 # 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, islistuple, map2, neg, _xinstanceof, _zip
 from pygeodesy.constants import R_M, _float as _F, _0_0, _1_0, _2_0, _8_0, _3600_0
@@ -94,7 +94,7 @@ from pygeodesy.units import _isRadius, Radius_,  radians
 # import operator as _operator  # from .fmath
 
 __all__ = _ALL_LAZY.datums
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _a_ellipsoid_ = _UNDER_(_a_, _ellipsoid_)
 _BD72_        = 'BD72'

pygeodesy/deprecated/__init__.py

@@ -10,13 +10,13 @@ bases} or C{from pygeodesy.deprecated import bases}.  Likewise for C{datum} and
 C{nvector}.
 '''
 
-from pygeodesy.deprecated.bases import *  # PYCHOK not pygeodesy.__init__
-from pygeodesy.deprecated.datum import *  # PYCHOK not pygeodesy.__init__
-from pygeodesy.deprecated.nvector import *  # PYCHOK not pygeodesy.__init__
+from pygeodesy.deprecated.bases import *  # noqa: F403 not pygeodesy.__init__
+from pygeodesy.deprecated.datum import *  # noqa: F403 not pygeodesy.__init__
+from pygeodesy.deprecated.nvector import *  # noqa: F403 not pygeodesy.__init__
 
-from pygeodesy.deprecated.classes import *  # PYCHOK expected
-from pygeodesy.deprecated.consterns import *  # PYCHOK expected
-from pygeodesy.deprecated.functions import *  # PYCHOK expected
+from pygeodesy.deprecated.classes import *  # noqa: F403
+from pygeodesy.deprecated.consterns import *  # noqa: F403
+from pygeodesy.deprecated.functions import *  # noqa: F403
 
 from pygeodesy.lazily import _ALL_ATTRS, _ALL_DEPRECATED, _lazy_import_as, _unLazy0  # _lazy_import_star
 
@@ -27,11 +27,11 @@ __all__ = (_ALL_DEPRECATED.deprecated_bases +
            _ALL_DEPRECATED.deprecated_classes +
            _ALL_DEPRECATED.deprecated_consterns +
            _ALL_DEPRECATED.deprecated_functions)
-__version__ = '24.12.31'
+__version__ = '25.05.12'
 
 if _unLazy0:
     from pygeodesy.deprecated import bases, datum, nvector, rhumbBase, \
-                                     rhumbaux, rhumbsolve, rhumbx  # PYCHOK expected
+                                     rhumbaux, rhumbsolve, rhumbx  # noqa: F401 PYCHOK expected
     __all__ += _ALL_ATTRS(_ALL_DEPRECATED.deprecated)  # DEPRECATED modules
 
 else:  # lazily import modules and exported attrs

pygeodesy/deprecated/bases.py

@@ -4,12 +4,12 @@
 u'''DEPRECATED on 2021.02.10, use (INTERNAL) module L{pygeodesy.latlonBase} instead.
 '''
 
-from pygeodesy.iters import points2  # PYCHOK exported
+from pygeodesy.iters import points2  # noqa: F401
 from pygeodesy.latlonBase import LatLonBase as _LatLonBase
 from pygeodesy.lazily import _ALL_DEPRECATED
 
 __all__ = _ALL_DEPRECATED.deprecated_bases
-__version__ = '25.04.09'
+__version__ = '25.05.12'
 
 
 class LatLonHeightBase(_LatLonBase):  # PYCHOK no cover

pygeodesy/deprecated/rhumbBase.py

@@ -4,10 +4,10 @@
 u'''DEPRECATED on 2023.11.26, use module L{pygeodesy.rhumb.bases} instead.
 '''
 
-from pygeodesy.rhumb.bases import *  # PYCHOK expected
+from pygeodesy.rhumb.bases import *  # noqa: F403
 
 __all__ = ()  # _ALL_LAZY.rhumb._bases
-__version__ = '23.11.26'
+__version__ = '25.05.12'
 
 # **) MIT License
 #

pygeodesy/deprecated/rhumbaux.py

@@ -4,11 +4,11 @@
 u'''DEPRECATED on 2023.11.26, use module L{pygeodesy.rhumb.aux_} instead.
 '''
 
-from pygeodesy.rhumb.aux_ import *  # PYCHOK expected
+from pygeodesy.rhumb.aux_ import *  # noqa: F403
 from pygeodesy.lazily import _ALL_LAZY
 
 __all__ = _ALL_LAZY.rhumb_aux_  # attrs NOT DEPRECATED
-__version__ = '23.12.02'
+__version__ = '25.05.12'
 
 # **) MIT License
 #

pygeodesy/deprecated/rhumbsolve.py

@@ -4,11 +4,11 @@
 u'''DEPRECATED on 2023.11.26, use module L{pygeodesy.rhumb.solve} instead.
 '''
 
-from pygeodesy.rhumb.solve import *  # PYCHOK expected
+from pygeodesy.rhumb.solve import *  # noqa: F403
 from pygeodesy.lazily import _ALL_LAZY
 
 __all__ = _ALL_LAZY.rhumb_solve  # attrs NOT DEPRECATED
-__version__ = '23.11.26'
+__version__ = '25.05.12'
 
 # **) MIT License
 #

pygeodesy/deprecated/rhumbx.py

@@ -4,11 +4,11 @@
 u'''DEPRECATED on 2023.11.26, use module L{pygeodesy.rhumb.ekx} instead.
 '''
 
-from pygeodesy.rhumb.ekx import *  # PYCHOK expected
+from pygeodesy.rhumb.ekx import *  # noqa: F403
 from pygeodesy.lazily import _ALL_LAZY
 
 __all__ = _ALL_LAZY.rhumb_ekx  # attrs NOT DEPRECATED
-__version__ = '23.11.30'
+__version__ = '25.05.12'
 
 # **) MIT License
 #

pygeodesy/ecef.py

@@ -60,9 +60,8 @@ plane} as opposed to I{geocentric} (ECEF) ones.
 from pygeodesy.basics import copysign0, _isin, isscalar, issubclassof, neg, map1, \
                             _xinstanceof, _xsubclassof,  typename  # _args_kwds_names
 from pygeodesy.constants import EPS, EPS0, EPS02, EPS1, EPS2, EPS_2, INT0, PI, PI_2, \
-                               _0_0, _0_0001, _0_01, _0_5, _1_0, _1_0_1T, _N_1_0, \
-                               _2_0, _N_2_0, _3_0, _4_0, _6_0, _60_0, _90_0, _N_90_0, \
-                               _100_0, _copysign_1_0,  isnon0  # PYCHOK used!
+                               _0_0, _0_5, _1_0, _1_0_1T, _2_0, _N_2_0, _3_0, _4_0, \
+                               _6_0, _90_0, _N_90_0, _copysign_1_0,  isnon0  # PYCHOK used!
 from pygeodesy.datums import _ellipsoidal_datum, _WGS84,  a_f2Tuple, _EWGS84
 from pygeodesy.ecefLocals import _EcefLocal
 # from pygeodesy.ellipsoids import a_f2Tuple, _EWGS84  # from .datums
@@ -90,7 +89,7 @@ from pygeodesy.utily import atan1, atan1d, atan2, atan2d, degrees90, degrees180,
 from math import cos, degrees, fabs, radians, sqrt
 
 __all__ = _ALL_LAZY.ecef
-__version__ = '25.04.28'
+__version__ = '25.05.07'
 
 _Ecef_    = 'Ecef'
 _prolate_ = 'prolate'

pygeodesy/ellipsoidalBase.py

@@ -10,7 +10,7 @@ and published under the same MIT Licence**, see for example U{latlon-ellipsoidal
 <https://www.Movable-Type.co.UK/scripts/geodesy/docs/latlon-ellipsoidal.js.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.azimuthal import EquidistantExact, EquidistantKarney  # _MODS
 from pygeodesy.basics import _isin, _xinstanceof
@@ -46,7 +46,7 @@ from pygeodesy.units import Epoch, _isDegrees, Radius_, _1mm as _TOL_M
 # from math import fabs  # from .latlonBase
 
 __all__ = _ALL_LAZY.ellipsoidalBase
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 
 class CartesianEllipsoidalBase(CartesianBase):

pygeodesy/ellipsoidalBaseDI.py

@@ -5,7 +5,7 @@ u'''(INTERNAL) Private, ellipsoidal Direct/Inverse geodesy base
 class C{LatLonEllipsoidalBaseDI} and functions.
 '''
 # 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.azimuthal import _Equidistants  # _MODS
 from pygeodesy.basics import isLatLon, _xsubclassof,  typename
@@ -36,7 +36,7 @@ from pygeodesy.utily import m2km, unroll180, _unrollon, _unrollon3, _Wrap, wrap3
 from math import degrees, radians
 
 __all__ = _ALL_LAZY.ellipsoidalBaseDI
-__version__ = '25.04.21'
+__version__ = '25.05.23'
 
 _polar__ = 'polar?'
 _TRIPS   =  33   # _intersect3, _intersects2, _nearestOn interations, 6..9 sufficient?
@@ -115,14 +115,15 @@ class LatLonEllipsoidalBaseDI(LatLonEllipsoidalBase):
             r = self._Direct2Tuple(LL, height, r)
         return r
 
-    def _Direct2Tuple(self, LL, height, r):
+    def _Direct2Tuple(self, LL, height, r):  # in .ellipsoidalVincenty.py
         '''(INTERNAL) Helper for C{._Direct} result L{Destination2Tuple}.
         '''
         h =  self._heigHt(height)
-        d = _xkwds_not(None, datum=self.datum, name=self.name,
+        n =  self.name
+        d = _xkwds_not(None, datum=self.datum, height=h, name=n,
                              epoch=self.epoch, reframe=self.reframe)
-        d = LL(*_Wrap.latlon(r.lat, r.lon), height=h, **d)
-        return Destination2Tuple(d, wrap360(r.final), name=self.name)
+        d = LL(*_Wrap.latlon(r.lat, r.lon), **d)
+        return Destination2Tuple(d, wrap360(r.final), name=n)
 
     def distanceTo(self, other, wrap=False, **unused):  # radius=R_M
         '''Compute the distance between this and an other point along

pygeodesy/ellipsoidalExact.py

@@ -11,12 +11,12 @@ L{GeodesicExact}, L{GeodesicAreaExact} and L{GeodesicLineExact}.
 
 # from pygeodesy.datums import _WGS84  # from .ellipsoidalBase
 from pygeodesy.ellipsoidalBase import CartesianEllipsoidalBase, \
-                                     _nearestOn, _WGS84
+                                     _nearestOn,  _WGS84
 from pygeodesy.ellipsoidalBaseDI import LatLonEllipsoidalBaseDI, \
                                        _intersection3, _intersections2, \
                                        _TOL_M, intersecant2
 # from pygeodesy.errors import _xkwds  # from .karney
-from pygeodesy.karney import fabs, _polygon,  Property_RO, _xkwds
+from pygeodesy.karney import _polygon,  fabs, Property_RO, _xkwds
 from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS, _ALL_OTHER
 from pygeodesy.points import _areaError, ispolar  # PYCHOK exported
 # from pygeodesy.props import Property_RO  # from .karney
@@ -24,7 +24,7 @@ from pygeodesy.points import _areaError, ispolar  # PYCHOK exported
 # from math import fabs  # from .karney
 
 __all__ = _ALL_LAZY.ellipsoidalExact
-__version__ = '24.11.06'
+__version__ = '25.05.12'
 
 
 class Cartesian(CartesianEllipsoidalBase):

pygeodesy/ellipsoidalNvector.py

@@ -21,7 +21,7 @@ See also I{Kenneth Gade} U{'A Non-singular Horizontal Position Representation'
 The Journal of Navigation (2010), vol 63, nr 3, pp 395-417.
 '''
 # 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, issubclassof, map2, _xinstanceof, \
                              _xsubclassof
@@ -32,7 +32,7 @@ from pygeodesy.ellipsoidalBase import CartesianEllipsoidalBase, \
                                      _nearestOn, LatLonEllipsoidalBase, \
                                      _TOL_M,  _Wrap
 from pygeodesy.errors import _xkwds, _xkwds_pop2
-# from pygeodesy.fmath import fdot  # from .nvectorBase
+from pygeodesy.fmath import fdot,  fabs
 # from pygeodesy.formy import _isequalTo  # _MODS
 from pygeodesy.interns import _Nv00_, _COMMASPACE_,  _pole_  # PYCHOK used!
 from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS, _ALL_OTHER
@@ -41,7 +41,7 @@ from pygeodesy.ltpTuples import Aer as _Aer, Ned as _Ned, Ned4Tuple, \
                                 sincos2d_, _xnamed
 # from pygeodesy.named import _xnamed  # from .ltpTuples
 from pygeodesy.nvectorBase import LatLonNvectorBase, NorthPole, NvectorBase, \
-                                  sumOf as _sumOf,  fabs, fdot
+                                  sumOf as _sumOf
 from pygeodesy.props import deprecated_class, deprecated_function, \
                             deprecated_method, Property_RO, property_RO
 from pygeodesy.streprs import Fmt, fstr, _xzipairs
@@ -51,7 +51,7 @@ from pygeodesy.units import Bearing, Distance, Height, Scalar
 # from math import fabs  # from .nvectorBase
 
 __all__ = _ALL_LAZY.ellipsoidalNvector
-__version__ = '25.04.21'
+__version__ = '25.05.12'
 
 
 class Ned(_Ned):

pygeodesy/ellipsoidalVincenty.py

@@ -48,7 +48,7 @@ or by converting to anothor datum:
     >>> p = p.toDatum(Datums.OSGB36)
 '''
 # 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.constants import EPS, EPS0, _0_0, _1_0, _2_0, _3_0, _4_0, _6_0
 # from pygeodesy.ecef import EcefVeness  # _MODS
@@ -59,7 +59,7 @@ from pygeodesy.ellipsoidalBaseDI import LatLonEllipsoidalBaseDI, \
 # from pygeodesy.ellipsoidalExact import areaOf, perimeterOf  # _MODS
 # from pygeodesy.ellipsoidalKarney import areaOf, perimeterOf  # _MODS
 from pygeodesy.errors import _and, _ValueError, _xkwds
-from pygeodesy.fmath import Fpolynomial, hypot, hypot1
+from pygeodesy.fmath import fdot_, Fpolynomial, hypot, hypot1
 from pygeodesy.interns import _ambiguous_, _antipodal_, _COLONSPACE_, \
                               _to_, _SPACE_,  _limit_  # PYCHOK used!
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS
@@ -76,7 +76,7 @@ from pygeodesy.utily import atan2, atan2b, atan2d, sincos2, sincos2d, \
 from math import cos, degrees, fabs, radians, tan as _tan
 
 __all__ = _ALL_LAZY.ellipsoidalVincenty
-__version__ = '25.04.21'
+__version__ = '25.05.23'
 
 _antipodal_to_ = _SPACE_(_antipodal_, _to_)
 
@@ -230,13 +230,13 @@ class LatLon(LatLonEllipsoidalBaseDI):
             t = self._no_convergence(e)
             raise VincentyError(t, txt=repr(self))  # self.toRepr()
 
-        t = s1 * ss - c1 * cs * cb
+        t = fdot_(s1, ss, -c1, cs * cb)
         # final bearing (reverse azimuth +/- 180)
         d = atan2b(sa, -t)
         if llr:
             b = cb * ss
-            a = atan2d(s1 * cs + c1 * b, hypot(sa, t) * E.b_a)
-            b = atan2d(sb * ss, -s1 * b + c1 * cs) + self.lon \
+            a = atan2d(fdot_(s1, cs, c1, b), hypot(sa, t) * E.b_a)
+            b = atan2d(sb * ss, fdot_(-s1, b, c1, cs)) + self.lon \
               - degrees(_Dl(f, ca2, sa, s, cs, ss, c2sm))
             t = Destination3Tuple(a, wrap180(b), d)
             r = self._Direct2Tuple(self.classof, height, t)
@@ -349,6 +349,13 @@ def _c2sm2(c2sm):
     return c2sm**2 * _2_0 - _1_0
 
 
+def _ellipsoidal():  # helper for areaOf and perimeterOf
+    try:
+        return _MODS.ellipsoidalKarney
+    except ImportError:
+        return _MODS.ellipsoidalExact
+
+
 def _Dl(f, ca2, sa, s, cs, ss, c2sm, dl=_0_0):
     # C{Dl}
     if f and sa:
@@ -399,10 +406,7 @@ def _sincostan3r(a, f):
 def areaOf(points, **datum_wrap):
     '''DEPRECATED, use function L{ellipsoidalExact.areaOf} or L{ellipsoidalKarney.areaOf}.
     '''
-    try:
-        return _MODS.ellipsoidalKarney.areaOf(points, **datum_wrap)
-    except ImportError:
-        return _MODS.ellipsoidalExact.areaOf(points, **datum_wrap)
+    return _ellipsoidal().areaOf(points, **datum_wrap)
 
 
 def intersection3(start1, end1, start2, end2, height=None, wrap=False,  # was=True
@@ -558,10 +562,7 @@ def nearestOn(point, point1, point2, within=True, height=None, wrap=False,
 def perimeterOf(points, **closed_datum_wrap):
     '''DEPRECATED, use function L{ellipsoidalExact.perimeterOf} or L{ellipsoidalKarney.perimeterOf}.
     '''
-    try:
-        return _MODS.ellipsoidalKarney.perimeterOf(points, **closed_datum_wrap)
-    except ImportError:
-        return _MODS.ellipsoidalExact.perimeterOf(points, **closed_datum_wrap)
+    return _ellipsoidal().perimeterOf(points, **closed_datum_wrap)
 
 
 __all__ += _ALL_DOCS(Cartesian, LatLon, intersecant2,  # from .ellipsoidalBaseDI

pygeodesy/ellipsoids.py

@@ -62,7 +62,7 @@ Following is the list of predefined L{Ellipsoid}s, all instantiated lazily.
 @var Ellipsoids.WGS84_NGS: Ellipsoid(name='WGS84_NGS', a=6378137, f=0.00335281, f_=298.2572221, b=6356752.31414035)
 '''
 # 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.albers import AlbersEqualAreaCylindrical  # _MODS
 from pygeodesy.basics import copysign0, isbool, _isin, isint,  typename
@@ -95,7 +95,7 @@ from pygeodesy.utily import atan1, atan1d, atan2b, degrees90, m2radians, radians
 from math import asinh, atan, atanh, cos, degrees, exp, fabs, radians, sin, sinh, sqrt, tan  # as _tan
 
 __all__ = _ALL_LAZY.ellipsoids
-__version__ = '25.04.28'
+__version__ = '25.05.12'
 
 _f_0_0    = Float(f =_0_0)  # zero flattening
 _f__0_0   = Float(f_=_0_0)  # zero inverse flattening
@@ -398,7 +398,7 @@ class Ellipsoid(_NamedEnumItem):
     def AlphaKs(self):
         '''Get the I{Krüger} U{Alpha series coefficients<https://GeographicLib.SourceForge.io/C++/doc/tmseries30.html>} (C{KsOrder}C{-tuple}).
         '''
-        return self._Kseries(  # XXX int/int quotients may require  from __future__ import division as _; del _  # PYCHOK semicolon
+        return self._Kseries(  # XXX int/int quotients may require  from __future__ import division as _; del _  # noqa: E702 ;
             #   n    n**2   n**3      n**4         n**5            n**6                 n**7                     n**8
             _T(1/2, -2/3,   5/16,    41/180,    -127/288,       7891/37800,         72161/387072,        -18975107/50803200),
                  _T(13/48, -3/5,    557/1440,    281/630,   -1983433/1935360,       13769/28800,         148003883/174182400),     # PYCHOK unaligned
@@ -617,7 +617,7 @@ class Ellipsoid(_NamedEnumItem):
     def BetaKs(self):
         '''Get the I{Krüger} U{Beta series coefficients<https://GeographicLib.SourceForge.io/C++/doc/tmseries30.html>} (C{KsOrder}C{-tuple}).
         '''
-        return self._Kseries(  # XXX int/int quotients may require  from __future__ import division as _; del _  # PYCHOK semicolon
+        return self._Kseries(  # XXX int/int quotients may require  from __future__ import division as _; del _  # noqa: E702 ;
             #   n    n**2  n**3     n**4        n**5            n**6                 n**7                   n**8
             _T(1/2, -2/3, 37/96,   -1/360,    -81/512,      96199/604800,     -5406467/38707200,      7944359/67737600),
                   _T(1/48, 1/15, -437/1440,    46/105,   -1118711/3870720,       51841/1209600,      24749483/348364800),      # PYCHOK unaligned

pygeodesy/elliptic.py

@@ -73,7 +73,7 @@ The notation follows U{NIST Digital Library of Mathematical Functions
 U{22<https://DLMF.NIST.gov/22>}.
 '''
 # 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, map2, neg, neg_,  typename
 from pygeodesy.constants import EPS, INF, NAN, PI, PI_2, PI_4, _0_0, \
@@ -99,7 +99,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.04.14'
+__version__ = '25.05.12'
 
 _TolRD  =  zqrt(EPS * 0.002)
 _TolRF  =  zqrt(EPS * 0.030)

pygeodesy/errors.py

@@ -27,7 +27,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.04.21'
+__version__ = '25.05.19'
 
 _argument_   = 'argument'
 _basics      = _MODS.into(basics=__name__)

pygeodesy/etm.py

@@ -61,7 +61,7 @@ control round-off.  For example, C{atanh(sin(phi))} is replaced by C{asinh(tan(p
 which maintains accuracy near C{phi = pi/2}.  Such changes are noted in the code.
 '''
 # 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, map1, neg, neg_, _xinstanceof
 from pygeodesy.constants import EPS, EPS02, PI_2, PI_4, _K0_UTM, \
@@ -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.04.14'
+__version__ = '25.05.12'
 
 _OVERFLOW = _1_EPS**2  # ~2e+31
 _TAYTOL   =  pow(EPS,  0.6)

pygeodesy/fmath.py

@@ -5,7 +5,7 @@ u'''Utilities for precision floating point summation, multiplication,
 C{fused-multiply-add}, polynomials, roots, etc.
 '''
 # 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, copysign0, isbool, isint, isscalar, \
                               len2, map1, _xiterable,  typename
@@ -19,13 +19,13 @@ from pygeodesy.fsums import _2float, Fsum, fsum, _isFsum_2Tuple,  Fmt, unstr
 from pygeodesy.interns import MISSING, _negative_, _not_scalar_
 from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS
 # from pygeodesy.streprs import Fmt, unstr  # from .fsums
-from pygeodesy.units import Int_, _isHeight, _isRadius,  Float_  # PYCHOK for .heights
+from pygeodesy.units import Int_, _isHeight, _isRadius
 
 from math import fabs, sqrt  # pow
 import operator as _operator  # in .datums, .trf, .utm
 
 __all__ = _ALL_LAZY.fmath
-__version__ = '25.04.30'
+__version__ = '25.05.12'
 
 # 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
@@ -37,7 +37,7 @@ class Fdot(Fsum):
     '''Precision dot product.
     '''
     def __init__(self, a, *b, **start_name_f2product_nonfinites_RESIDUAL):
-        '''New L{Fdot} precision dot product M{sum(a[i] * b[i] for i=0..len(a)-1)}.
+        '''New L{Fdot} precision dot product M{start + sum(a[i] * b[i] for i=0..len(a)-1)}.
 
            @arg a: Iterable of values (each C{scalar}, an L{Fsum} or L{Fsum2Tuple}).
            @arg b: Other values (each C{scalar}, an L{Fsum} or L{Fsum2Tuple}), all
@@ -70,7 +70,7 @@ class Fdot_(Fdot):
     '''Precision dot product.
     '''
     def __init__(self, *xys, **start_name_f2product_nonfinites_RESIDUAL):
-        '''New L{Fdot_} precision dot product M{sum(xys[i] * xys[i+1] for i in
+        '''New L{Fdot_} precision dot product M{start + sum(xys[i] * xys[i+1] for i in
            range(0, len(xys), B{2}))}.
 
            @arg xys: Pairwise values (each C{scalar}, an L{Fsum} or L{Fsum2Tuple}),
@@ -439,7 +439,7 @@ def favg(a, b, f=_0_5, nonfinites=True):
 
 
 def fdot(xs, *ys, **start_f2product_nonfinites):
-    '''Return the precision dot product M{sum(xs[i] * ys[i] for i in range(len(xs)))}.
+    '''Return the precision dot product M{start + sum(xs[i] * ys[i] for i in range(len(xs)))}.
 
        @arg xs: Iterable of values (each C{scalar}, an L{Fsum} or L{Fsum2Tuple}).
        @arg ys: Other values (each C{scalar}, an L{Fsum} or L{Fsum2Tuple}), all positional.
@@ -460,7 +460,7 @@ def fdot(xs, *ys, **start_f2product_nonfinites):
 
 
 def fdot_(*xys, **start_f2product_nonfinites):
-    '''Return the (precision) dot product M{sum(xys[i] * xys[i+1] for i in range(0, len(xys), B{2}))}.
+    '''Return the (precision) dot product M{start + sum(xys[i] * xys[i+1] for i in range(0, len(xys), B{2}))}.
 
        @arg xys: Pairwise values (each C{scalar}, an L{Fsum} or L{Fsum2Tuple}), all positional.
 
@@ -987,7 +987,7 @@ def _powers(x, n):
 
 
 def _root(x, p, where):
-    '''(INTERNAL) Raise C{x} to power C{0 < p < 1}.
+    '''(INTERNAL) Raise C{x} to power C{0 <= p < 1}.
     '''
     try:
         if x > 0:
@@ -997,7 +997,7 @@ def _root(x, p, where):
             raise ValueError(_negative_)
     except Exception as X:
         raise _xError(X, unstr(where, x))
-    return _0_0 if p else _1_0
+    return _0_0 if p else _1_0  # x == 0
 
 
 def sqrt0(x, Error=None):

pygeodesy/formy.py

@@ -4,7 +4,7 @@
 u'''Formulary of basic geodesy functions and approximations.
 '''
 # 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  # _args_kwds_count2
 # from pygeodesy.cartesianBase import CartesianBase  # _MODS
@@ -42,7 +42,7 @@ from contextlib import contextmanager
 from math import atan, cos, degrees, fabs, radians, sin, sqrt  # pow
 
 __all__ = _ALL_LAZY.formy
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _RADIANS2 =  radians(_1_0)**2  # degree to radians-squared
 _ratio_   = 'ratio'