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

pygeodesy/ellipsoidalKarney.py

@@ -42,7 +42,7 @@ from pygeodesy.props import deprecated_method, Property_RO
 # from math import fabs  # from .karney
 
 __all__ = _ALL_LAZY.ellipsoidalKarney
-__version__ = '24.08.13'
+__version__ = '25.05.27'
 
 
 class Cartesian(CartesianEllipsoidalBase):
@@ -110,13 +110,17 @@ class LatLon(LatLonEllipsoidalBaseDI):
         return LatLonEllipsoidalBaseDI.toCartesian(self, **kwds)
 
 
-def areaOf(points, datum=_WGS84, wrap=True):
+def areaOf(points, datum=_WGS84, wrap=True, polar=False):
     '''Compute the area of an (ellipsoidal) polygon or composite using I{Karney}'s
        U{geographiclib<https://PyPI.org/project/geographiclib>} package.
 
        @arg points: The polygon points (L{LatLon}[], L{BooleanFHP} or L{BooleanGH}).
        @kwarg datum: Optional datum (L{Datum}).
        @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the B{C{points}} (C{bool}).
+       @kwarg polar: Use C{B{polar}=True} if the polygon encloses a pole (C{bool}), see
+                     function L{ispolar<pygeodesy.points.ispolar>} and U{area of a polygon
+                     enclosing a pole<https://GeographicLib.SourceForge.io/C++/doc/
+                     classGeographicLib_1_1GeodesicExact.html#a3d7a9155e838a09a48dc14d0c3fac525>}.
 
        @return: Area (C{meter}, same as units of the B{C{datum}}'s ellipsoid axes, I{squared}).
 
@@ -131,12 +135,8 @@ def areaOf(points, datum=_WGS84, wrap=True):
 
        @see: Functions L{pygeodesy.areaOf}, L{ellipsoidalExact.areaOf}, L{ellipsoidalGeodSolve.areaOf},
              L{sphericalNvector.areaOf} and L{sphericalTrigonometry.areaOf}.
-
-       @note: The U{area of a polygon enclosing a pole<https://GeographicLib.SourceForge.io/
-              C++/doc/classGeographicLib_1_1GeodesicExact.html#a3d7a9155e838a09a48dc14d0c3fac525>}
-              can be found by adding half the datum's ellipsoid surface area to the polygon's area.
     '''
-    return fabs(_polygon(datum.ellipsoid.geodesic, points, True, False, wrap))
+    return fabs(_polygon(datum.ellipsoid.geodesic, points, True, False, wrap, polar))
 
 
 def intersection3(start1, end1, start2, end2, height=None, wrap=False,  # was=True
@@ -240,19 +240,19 @@ def intersections2(center1, radius1, center2, radius2, height=None, wrap=False,
                            equidistant=equidistant, tol=tol, LatLon=LatLon, **LatLon_kwds)
 
 
-def isclockwise(points, datum=_WGS84, wrap=True):
+def isclockwise(points, datum=_WGS84, wrap=True, polar=False):
     '''Determine the direction of a path or polygon using I{Karney}'s
        U{geographiclib<https://PyPI.org/project/geographiclib>} package.
 
        @arg points: The path or polygon points (C{LatLon}[]).
        @kwarg datum: Optional datum (L{Datum}).
-       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the
-                    B{C{points}} (C{bool}).
+       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the B{C{points}} (C{bool}).
+       @kwarg polar: Use C{B{polar}=True} if the C{B{points}} enclose a pole (C{bool}),
+                     see function U{ispolar<pygeodeys.points.ispolar>}.
 
        @return: C{True} if B{C{points}} are clockwise, C{False} otherwise.
 
-       @raise ImportError: Package U{geographiclib
-                           <https://PyPI.org/project/geographiclib>}
+       @raise ImportError: Package U{geographiclib<https://PyPI.org/project/geographiclib>}
                            not installed or not found.
 
        @raise PointsError: Insufficient number of B{C{points}}.
@@ -263,7 +263,7 @@ def isclockwise(points, datum=_WGS84, wrap=True):
 
        @see: L{pygeodesy.isclockwise}.
     '''
-    a = _polygon(datum.ellipsoid.geodesic, points, True, False, wrap)
+    a = _polygon(datum.ellipsoid.geodesic, points, True, False, wrap, polar)
     if a < 0:
         return True
     elif a > 0:
@@ -322,34 +322,29 @@ def perimeterOf(points, closed=False, datum=_WGS84, wrap=True):
     '''Compute the perimeter of an (ellipsoidal) polygon or composite using I{Karney}'s
        U{geographiclib<https://PyPI.org/project/geographiclib>} package.
 
-       @arg points: The polygon points (L{LatLon}[], L{BooleanFHP} or
-                    L{BooleanGH}).
+       @arg points: The polygon points (L{LatLon}[], L{BooleanFHP} or L{BooleanGH}).
        @kwarg closed: Optionally, close the polygon (C{bool}).
        @kwarg datum: Optional datum (L{Datum}).
-       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the
-                    B{C{points}} (C{bool}).
+       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the B{C{points}} (C{bool}).
 
-       @return: Perimeter (C{meter}, same as units of the B{C{datum}}'s
-                ellipsoid axes).
+       @return: Perimeter (C{meter}, same as units of the B{C{datum}}'s ellipsoid axes).
 
-       @raise ImportError: Package U{geographiclib
-                           <https://PyPI.org/project/geographiclib>}
+       @raise ImportError: Package U{geographiclib<https://PyPI.org/project/geographiclib>}
                            not installed or not found.
 
        @raise PointsError: Insufficient number of B{C{points}}.
 
-       @raise TypeError: Some B{C{points}} are not L{LatLon} or C{B{closed}=False}
-                         with B{C{points}} a composite.
+       @raise TypeError: Some B{C{points}} are not L{LatLon} or C{B{closed}=False} with
+                         B{C{points}} a composite.
 
-       @raise ValueError: Invalid C{B{wrap}=False}, unwrapped, unrolled
-                          longitudes not supported or C{B{closed}=False}
-                          with C{B{points}} a composite.
+       @raise ValueError: Invalid C{B{wrap}=False}, unwrapped, unrolled longitudes not
+                          supported or C{B{closed}=False} with C{B{points}} a composite.
 
        @see: Functions L{pygeodesy.perimeterOf}, L{ellipsoidalExact.perimeterOf},
              L{ellipsoidalGeodSolve.perimeterOf}, L{sphericalNvector.perimeterOf}
              and L{sphericalTrigonometry.perimeterOf}.
     '''
-    return _polygon(datum.ellipsoid.geodesic, points, closed, True, wrap)
+    return _polygon(datum.ellipsoid.geodesic, points, closed, True, wrap, False)
 
 
 __all__ += _ALL_OTHER(Cartesian, LatLon,  # classes

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.26'
 
 _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)
@@ -378,6 +378,14 @@ def _Ecef():
     return E
 
 
+def _ellipsoidalOf(_Of, points, **kwds):  # helper for DEPRECATED areaOf and perimeterOf
+    try:
+        r = getattr(_MODS.ellipsoidalKarney, _Of.__name__)(points, **kwds)
+    except ImportError:  # no geographiclib
+        r = getattr(_MODS.ellipsoidalExact, _Of.__name__)(points, **kwds)
+    return r
+
+
 def _sincos22(sa):
     # 2-Tuple C{(sin(a), cos(a)**2)}
     ca2 = _1_0 - sa**2
@@ -396,13 +404,10 @@ def _sincostan3r(a, f):
 
 
 @deprecated_function
-def areaOf(points, **datum_wrap):
+def areaOf(points, **datum_wrap_polar):
     '''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 _ellipsoidalOf(areaOf, points, **datum_wrap_polar)
 
 
 def intersection3(start1, end1, start2, end2, height=None, wrap=False,  # was=True
@@ -558,10 +563,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 _ellipsoidalOf(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.28'
 
 _TolRD  =  zqrt(EPS * 0.002)
 _TolRF  =  zqrt(EPS * 0.030)
@@ -165,9 +165,9 @@ class Elliptic(_Named):
                   that case, we have C{Π(φ, 0, k) = F(φ, k), G(φ, 0, k) =
                   E(φ, k)} and C{H(φ, 0, k) = F(φ, k) - D(φ, k)}.
         '''
-        self.reset(k2=k2, alpha2=alpha2, kp2=kp2, alphap2=alphap2)
         if name:
             self.name = name
+        self.reset(k2=k2, alpha2=alpha2, kp2=kp2, alphap2=alphap2)
 
     @Property_RO
     def alpha2(self):
@@ -760,15 +760,20 @@ class Elliptic(_Named):
                   that these conditions are met to enable accuracy to be
                   maintained, e.g., when C{k} is very close to unity.
         '''
+        def _1p2(kp2, k2):
+            return (_1_0 - k2) if kp2 is None else kp2
+
+        def _S(**kwds):
+            return Scalar_(Error=EllipticError, **kwds)
+
         if self.__dict__:
             _update_all(self, _Named.iteration._uname, Base=Property_RO)
 
-        self._k2  = Scalar_(k2=k2, Error=EllipticError, low=None, high=_1_0)
-        self._kp2 = Scalar_(kp2=((_1_0 - k2) if kp2 is None else kp2), Error=EllipticError)
+        self._k2      = _S(k2 = k2, low=None, high=_1_0)
+        self._kp2     = _S(kp2=_1p2(kp2, k2))  # low=_0_0
 
-        self._alpha2  = Scalar_(alpha2=alpha2, Error=EllipticError, low=None, high=_1_0)
-        self._alphap2 = Scalar_(alphap2=((_1_0 - alpha2) if alphap2 is None else alphap2),
-                                Error=EllipticError)
+        self._alpha2  = _S(alpha2 = alpha2, low=None, high=_1_0)
+        self._alphap2 = _S(alphap2=_1p2(alphap2, alpha2))  # low=_0_0
 
         # Values of complete elliptic integrals for k = 0,1 and alpha = 0,1
         #         K     E     D
@@ -1120,11 +1125,11 @@ def _rC(unused, x, y):
     '''(INTERNAL) Defined only for C{y != 0} and C{x >= 0}.
     '''
     d = x - y
-    if d < 0:  # catch NaN
+    if x < y:  # catch NaN
         # <https://DLMF.NIST.gov/19.2.E18>
-        d = -d
+        d = y - x
         r = atan(sqrt(d / x)) if x > 0 else PI_2
-    elif d == 0:  # XXX d < EPS0? or EPS02 or _EPSmin
+    elif x == y:  # XXX d < EPS0? or EPS02 or _EPSmin
         d, r = y, _1_0
     elif y > 0:  # <https://DLMF.NIST.gov/19.2.E19>
         r = asinh(sqrt(d / y))  # atanh(sqrt((x - y) / x))

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'

pygeodesy/frechet.py

@@ -90,9 +90,9 @@ from pygeodesy.interns import _DMAIN_, _DOT_, _n_, NN, _units_
 # from pygeodesy.iters import points2 as _points2  # from .points
 from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS, _FOR_DOCS
 from pygeodesy.named import _name2__, _Named, _NamedTuple, _Pass
-# from pygeodesy.namedTuples import PhiLam2Tuple  # from .points
-from pygeodesy.points import _distanceTo, _fractional, PhiLam2Tuple, \
-                              points2 as _points2,  radians
+from pygeodesy.namedTuples import PhiLam2Tuple
+from pygeodesy.points import _distanceTo, _fractional, \
+                              points2 as _points2, radians
 from pygeodesy.props import Property, property_doc_, property_RO
 from pygeodesy.units import FIx, Float, Number_
 from pygeodesy import unitsBase as _unitsBase  # _Str_..., _xUnit, _xUnits
@@ -101,7 +101,7 @@ from collections import defaultdict as _defaultdict
 # from math import radians  # from .points
 
 __all__ = _ALL_LAZY.frechet
-__version__ = '25.04.21'
+__version__ = '25.05.26'
 
 _formy = _MODS.into(formy=__name__)
 
@@ -129,7 +129,7 @@ class Frechet(_Named):
        be overloaded.
     '''
     _datum = _WGS84
-    _func  =  None  # formy function/property
+#   _func  =  None  # formy function/property
     _f1    =  1
     _kwds  = {}     # func_ options
     _n1    =  0
@@ -153,7 +153,7 @@ class Frechet(_Named):
            @raise FrechetError: Insufficient number of B{C{point1s}} or an invalid
                                 B{C{point1}}, B{C{fraction}} or B{C{units}}.
         '''
-        name, kwds = _name2__(**name__kwds)  # name__=self.__class__
+        name, kwds = _name2__(**name__kwds)  # name__=type(self)
         if name:
             self.name = name
 

pygeodesy/fstats.py

@@ -6,7 +6,7 @@ L{pygeodesy.Fsum}, precision floating point summation and accurate
 multiplication.
 '''
 # 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 isscalar, isodd, _xinstanceof, \
                             _xiterable, _xsubclassof, _zip,  typename
@@ -24,7 +24,7 @@ from pygeodesy.named import _name__, _Named, _NotImplemented, \
 # from pygeodesy.streprs import Fmt  # from .fmath
 
 __all__ = _ALL_LAZY.fstats
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 
 def _sampled(n, sample):