pygeodesy 24.7.7__py2.py3-none-any.whl → 24.7.24__py2.py3-none-any.whl

pygeodesy/ellipsoids.py

@@ -83,7 +83,7 @@ from pygeodesy.named import _lazyNamedEnumItem as _lazy, _name__, _NamedEnum, \
                                 _NamedEnumItem, _NamedTuple, _Pass,  _ALL_LAZY, _MODS
 from pygeodesy.namedTuples import Distance2Tuple, Vector3Tuple, Vector4Tuple
 from pygeodesy.props import deprecated_Property_RO, Property_RO, property_doc_, \
-                            deprecated_property_RO, property_RO
+                            deprecated_property_RO, property_RO, property_ROver
 from pygeodesy.streprs import Fmt, fstr, instr, strs, unstr
 # from pygeodesy.triaxials import _hartzell3  # _MODS
 from pygeodesy.units import Bearing_, Distance, Float, Float_, Height, Lamd, Lat, Meter, \
@@ -93,7 +93,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
 
 __all__ = _ALL_LAZY.ellipsoids
-__version__ = '24.06.24'
+__version__ = '24.07.16'
 
 _f_0_0    = Float(f =_0_0)  # zero flattening
 _f__0_0   = Float(f_=_0_0)  # zero inverse flattening
@@ -1029,17 +1029,18 @@ class Ellipsoid(_NamedEnumItem):
                 raise _ValueError(exact=g, ellipsoid=E, txt_not_=self.name)
         return g
 
-    @property_RO
+    @property_ROver
     def _Geodesics(self):
         '''(INTERNAL) Get all C{Geodesic...} classes, I{once}.
         '''
+        t = (_MODS.geodesicx.GeodesicExact,
+             _MODS.geodsolve.GeodesicSolve)
         try:
-            t = _MODS.geodesicw._wrapped.Geodesic,
+            t += (_MODS.geodesicw.Geodesic,
+                  _MODS.geodesicw._wrapped.Geodesic)
         except ImportError:
-            t = ()
-        Ellipsoid._Geodesics = t = (_MODS.geodesicx.GeodesicExact,  # overwrite property_RO
-                                    _MODS.geodsolve.GeodesicSolve) + t
-        return t
+            pass
+        return t  # overwrite propertyROver
 
     @property_RO
     def geodesicw(self):
@@ -1517,14 +1518,13 @@ class Ellipsoid(_NamedEnumItem):
         #     raise _IsnotError(_ellipsoidal_, ellipsoid=self)
         return _MODS.rhumb.ekx.Rhumb(self, name=self.name)
 
-    @property_RO
+    @property_ROver
     def _Rhumbs(self):
         '''(INTERNAL) Get all C{Rhumb...} classes, I{once}.
         '''
         p = _MODS.rhumb
-        Ellipsoid._Rhumbs = t = (p.aux_.RhumbAux,  # overwrite property_RO
-                                 p.ekx.Rhumb, p.solve.RhumbSolve)
-        return t
+        return (p.aux_.RhumbAux,  # overwrite propertyROver
+                p.ekx.Rhumb, p.solve.RhumbSolve)
 
     @property
     def rhumbsolve(self):

pygeodesy/fsums.py

@@ -273,11 +273,17 @@ def _stresidual(prefix, residual, R=0, **mod_ratio):
 def _2sum(a, b):  # by .testFmath
     '''(INTERNAL) Return C{a + b} as 2-tuple (sum, residual).
     '''
+    # Neumaier, A. U{Rundungsfehleranalyse einiger Verfahren zur Summation endlicher
+    # Summen<https://OnlineLibrary.Wiley.com/doi/epdf/10.1002/zamm.19740540106>},
+    # 1974, Zeitschrift für Angewandte Mathmatik und Mechanik, vol 51, nr 1, p 39-51
+    # <https://StackOverflow.com/questions/78633770/can-neumaier-summation-be-sped-up>
     s = a + b
     if _isfinite(s):
         if fabs(a) < fabs(b):
-            b, a = a, b
-        return s, (b - (s - a))
+            r = (b - s) + a
+        else:
+            r = (a - s) + b
+        return s, r
     u = unstr(_2sum, a, b)
     t = Fmt.PARENSPACED(_not_finite_, s)
     raise _OverflowError(u, txt=t)
@@ -1909,10 +1915,10 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase
 
            @raise ResidualError: Invalid B{C{threshold}}.
 
-           @note: L{ResidualError}s may be thrown if the non-zero I{ratio}
-                  C{residual / fsum} exceeds the given B{C{threshold}} and
-                  if the C{residual} is non-zero and I{significant} vs the
-                  C{fsum}, i.e. C{(fsum + residual) != fsum} and if optional
+           @note: L{ResidualError}s may be thrown if (1) the non-zero I{ratio}
+                  C{residual / fsum} exceeds the given B{C{threshold}} and (2)
+                  the C{residual} is non-zero and (3) I{significant} vs the
+                  C{fsum}, i.e. C{(fsum + residual) != fsum} and (4) optional
                   keyword argument C{raiser=False} is missing.  Specify a
                   negative B{C{threshold}} for only non-zero C{residual}
                   testing without I{significant}.