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

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.23'
+__version__ = '25.05.12'
 
 _f_0_0    = Float(f =_0_0)  # zero flattening
 _f__0_0   = Float(f_=_0_0)  # zero inverse flattening
@@ -140,12 +140,14 @@ class a_f2Tuple(_NamedTuple):
            @arg f: Flattening (C{scalar} < 1, negative for I{prolate}).
            @kwarg name: Optional C{B{name}=NN} (C{str}).
 
-           @return: An L{a_f2Tuple}C{(a, f)} instance.
+           @return: An L{a_f2Tuple}C{(a, f)}.
 
            @raise UnitError: Invalid B{C{a}} or B{C{f}}.
 
-           @note: C{abs(B{f}) < EPS} is forced to C{B{f}=0}, I{spherical}.
-                  Negative C{B{f}} produces a I{prolate} ellipsoid.
+           @note: C{abs(B{f}) < }L{EPS<pygeodesy.constants.EPS>} is
+                  forced to C{B{f}=0}, I{spherical}.
+
+           @note: Negative C{B{f}} produces a I{prolate} ellipsoid.
         '''
         a = Radius_(a=a)  # low=EPS, high=None
         f = Float_( f=f, low=None, high=EPS1)
@@ -359,11 +361,16 @@ class Ellipsoid(_NamedEnumItem):
         '''
         return a_f2Tuple(self.a, self.f, name=self.name)
 
+    a_f2 = a_f  # synonym
+
     @Property_RO
     def A(self):
         '''Get the UTM I{meridional (or rectifying)} radius (C{meter}).
 
-           @see: I{Meridian arc unit} U{Q<https://StudyLib.net/doc/7443565/>}.
+           @note: C{A * PI / 2} ≈= L{L<Ellipsoid.L>}, the I{quarter meridian}.
+
+           @see: I{Meridian arc unit} U{Q<https://StudyLib.net/doc/7443565/>},
+                 the mean, meridional length I{per radian}.
         '''
         A, n = self.a, self.n
         if n:
@@ -375,6 +382,11 @@ class Ellipsoid(_NamedEnumItem):
                 # A *= fhorner(n**2, 1048576, 262144, 16384, 4096, 1600, 784, 441) / 1048576) / (1 + n)
                 A = Radius(A=Fhorner(n**2, 1048576, 262144, 16384, 4096, 1600, 784, 441).fover(d))
         return A
+#       # Moritz, H. <https://Geodesy.Geology.Ohio-State.EDU/course/refpapers/00740128.pdf>
+#       # q =  4 / self.rocPolar
+#       # Q = (1 - 3 / 4 * e'2 + 45 / 64 * e'4 - 175 / 256 * e'6 + 11025 / 16384 * e'8) * rocPolar
+#       #   = (4 + e'2 * (-3 + e'2 * (45 / 16 + e'2 * (-175 / 64 + e'2 * 11025 / 4096)))) / q
+#       return Radius(Q=Fhorner(self.e22, 4, -3, 45 / 16, -175 / 64, 11025 / 4096).fover(q))
 
     @Property_RO
     def _albersCyl(self):
@@ -386,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
@@ -605,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
@@ -1339,23 +1351,7 @@ class Ellipsoid(_NamedEnumItem):
 
     polaradius = b  # Rpolar
 
-#   @Property_RO
-#   def Q(self):
-#       '''Get the I{meridian arc unit} C{Q}, the mean, meridional length I{per radian} C({float}).
-#
-#          @note: C{Q * PI / 2} ≈ C{L}, the I{quarter meridian}.
-#
-#          @see: Property C{A} and U{Engsager, K., Poder, K.<https://StudyLib.net/doc/7443565/
-#                a-highly-accurate-world-wide-algorithm-for-the-transverse...>}.
-#       '''
-#       n = self.n
-#       d = (n + _1_0) / self.a
-#       return Float(Q=Fhorner(n**2, _1_0, _0_25, _1_16th, _0_25).fover(d) if d else self.b)
-
-#       # Moritz, H. <https://Geodesy.Geology.Ohio-State.EDU/course/refpapers/00740128.pdf>
-#       # Q = (1 - 3/4 * e'2 + 45/64 * e'4 - 175/256 * e'6 + 11025/16384 * e'8) * rocPolar
-#       #   = (4 + e'2 * (-3 + e'2 * (45/16 + e'2 * (-175/64 + e'2 * 11025/4096)))) * rocPolar / 4
-#       return Fhorner(self.e22, 4, -3, 45 / 16, -175 / 64, 11025 / 4096).fover(4 / self.rocPolar)
+#   Q = A  # I{meridian arc unit} C{Q}, the mean, meridional length I{per radian}
 
     @deprecated_Property_RO
     def quarteradius(self):  # PYCHOK no cover
@@ -1382,16 +1378,13 @@ class Ellipsoid(_NamedEnumItem):
                  <https://WikiPedia.org/wiki/Earth_radius>}.
         '''
         return Radius(R2=sqrt(self.c2) if self.f else self.a)
-
-    Rauthalic = R2
-
-#   @Property_RO
-#   def R2(self):
 #       # Moritz, H. <https://Geodesy.Geology.Ohio-State.EDU/course/refpapers/00740128.pdf>
 #       # R2 = (1 - 2/3 * e'2 + 26/45 * e'4 - 100/189 * e'6 + 7034/14175 * e'8) * rocPolar
 #       #    = (3 + e'2 * (-2 + e'2 * (26/15 + e'2 * (-100/63 + e'2 * 7034/4725)))) * rocPolar / 3
 #       return Fhorner(self.e22, 3, -2, 26 / 15, -100 / 63, 7034 / 4725).fover(3 / self.rocPolar)
 
+    Rauthalic = R2
+
     @Property_RO
     def R2x(self):
         '''Get the I{authalic} earth radius (C{meter}), M{sqrt(c2x)}.

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.18'
+__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
@@ -66,13 +66,27 @@ class Fdot(Fsum):
         self._facc_dot(n, a, b, **kwds)
 
 
+class Fdot_(Fdot):
+    '''Precision dot product.
+    '''
+    def __init__(self, *xys, **start_name_f2product_nonfinites_RESIDUAL):
+        '''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}),
+                     all positional.
+
+           @see: Class L{Fdot<Fdot.__init__>} for further details.
+        '''
+        Fdot.__init__(self, xys[0::2], *xys[1::2], **start_name_f2product_nonfinites_RESIDUAL)
+
+
 class Fhorner(Fsum):
     '''Precision polynomial evaluation using the Horner form.
     '''
     def __init__(self, x, *cs, **incx_name_f2product_nonfinites_RESIDUAL):
-        '''New L{Fhorner} form evaluation of polynomial M{sum(cs[i] * x**i for
-           i=0..n)} with in- or decreasing exponent M{sum(... i=n..0)}, where C{n
-           = len(cs) - 1}.
+        '''New L{Fhorner} form evaluation of polynomial M{sum(cs[i] * x**i for i=0..n)}
+           with in- or decreasing exponent M{sum(... i=n..0)}, where C{n = len(cs) - 1}.
 
            @arg x: Polynomial argument (C{scalar}, an L{Fsum} or L{Fsum2Tuple}).
            @arg cs: Polynomial coeffients (each C{scalar}, an L{Fsum} or L{Fsum2Tuple}),
@@ -425,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.
@@ -446,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.
 
@@ -454,7 +468,8 @@ def fdot_(*xys, **start_f2product_nonfinites):
 
        @return: Dot product (C{float}).
     '''
-    return fdot(xys[0::2], *xys[1::2], **start_f2product_nonfinites)
+    D = Fdot_(*xys, **_xkwds(start_f2product_nonfinites, nonfinites=True))
+    return float(D)
 
 
 def fdot3(xs, ys, zs, **start_f2product_nonfinites):
@@ -555,7 +570,7 @@ def fidw(xs, ds, beta=2):
 
 
 try:
-    from math import fma as _fma
+    from math import fma as _fma  # in .resections
 except  ImportError:  # PYCHOK DSPACE!
 
     def _fma(x, y, z):  # no need for accuracy
@@ -972,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:
@@ -982,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.21'
 
 _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

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):