pygeodesy 24.8.24__py2.py3-none-any.whl → 24.9.24__py2.py3-none-any.whl

pygeodesy/booleans.py

@@ -24,10 +24,10 @@ from pygeodesy.errors import ClipError, _IsnotError, _TypeError, \
 from pygeodesy.fmath import favg, hypot, hypot2
 # from pygeodesy.fsums import fsum1  # _MODS
 from pygeodesy.interns import NN, _BANG_, _clip_, _clipid_, _COMMASPACE_, \
-                             _composite_, _DOT_, _e_, _ELLIPSIS_, _few_, \
-                             _height_, _lat_, _LatLon_, _lon_, _not_, \
-                             _points_, _SPACE_, _too_, _X_, _x_, \
-                             _B_, _d_, _R_  # PYCHOK used!
+                             _composite_, _DOT_, _duplicate_, _e_, \
+                             _ELLIPSIS_, _few_, _height_, _lat_, _LatLon_, \
+                             _lon_, _not_, _points_, _SPACE_, _too_, _X_, \
+                             _x_,  _B_, _d_, _R_  # PYCHOK used!
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS
 from pygeodesy.latlonBase import LatLonBase, \
                                  LatLon2Tuple, Property_RO, property_RO
@@ -43,7 +43,7 @@ from pygeodesy.utily import fabs, _unrollon, _Wrap
 # from math import fabs  # from .utily
 
 __all__ = _ALL_LAZY.booleans
-__version__ = '24.06.15'
+__version__ = '24.09.23'
 
 _0_EPS =  EPS  # near-zero, positive
 _EPS_0 = -EPS  # near-zero, negative
@@ -51,12 +51,11 @@ _1_EPS = _1_0 + EPS   # near-one, over
 _EPS_1 = _1_0 - EPS   # near-one, under
 _10EPS =  EPS * 10    # see ._2Abs, ._10eps
 
-_alpha_     = 'alpha'
-_boolean_   = 'boolean'
-_case_      = 'case'
-_corners_   = 'corners'
-_duplicate_ = 'duplicate'
-_open_      = 'open'
+_alpha_   = 'alpha'
+_boolean_ = 'boolean'
+_case_    = 'case'
+_corners_ = 'corners'
+_open_    = 'open'
 
 
 def _Enum(txt, enum):  # PYCHOK unused
@@ -1118,7 +1117,7 @@ class _CompositeBase(_Named):
 
     def _sum1(self, _a_p, *args, **kwds):  # in .karney, .points
         # Sum the area or perimeter of all clips
-        return _MODS.fsums.fsum1((_a_p(c, *args, **kwds) for c in self._clips), floats=True)
+        return _MODS.fsums.fsum1((_a_p(c, *args, **kwds) for c in self._clips))
 
     def _sum2(self, LL, _a_p, *args, **kwds):  # in .sphericalNvector, -Trigonometry
         # Sum the area or perimeter of all clips
@@ -1128,7 +1127,7 @@ class _CompositeBase(_Named):
             for v in clip:
                 yield _LL(v.lat, v.lon)  # datum=Sphere
 
-        return _MODS.fsums.fsum1((_a_p(_lls(c), *args, **kwds) for c in self._clips), floats=True)
+        return _MODS.fsums.fsum1((_a_p(_lls(c), *args, **kwds) for c in self._clips))
 
     def toLatLon(self, LatLon=None, closed=False, **LatLon_kwds):
         '''Yield all (non-duplicate) points and intersections
@@ -1188,7 +1187,7 @@ class _CompositeFHP(_CompositeBase):
             n, b = v, v._label
             if b in L.RIGHT_LEFT_ON:  # next chain
                 while True:
-                    n._label = None  # n.__dict__.pop('_label')
+                    n._label = None  # _xkwds_pop(n.__dict__, _label=None)
                     n = n._next
                     if n is v or n._label is not L.ON_ON:  # n._label and ...
                         break

pygeodesy/clipy.py

@@ -15,22 +15,22 @@ from __future__ import division as _; del _  # PYCHOK semicolon
 # from pygeodesy.basics import len2  # from .fmath
 from pygeodesy.constants import EPS, _0_0, _1_0
 from pygeodesy.errors import _AssertionError, ClipError, PointsError
-from pygeodesy.fmath import fabs, len2
-from pygeodesy.fsums import fsumf_, Property_RO
+from pygeodesy.fmath import fabs, len2,  Fsum
+# from pygeodesy.fsums import Fsum  # from .fmath
 from pygeodesy.interns import NN, _clipid_, _convex_, _DOT_, _end_, _few_, \
                              _fi_, _height_, _i_, _invalid_, _j_, _lat_, \
                              _lon_, _near_, _not_, _points_, _start_, _too_
 from pygeodesy.iters import _imdex2, points2
 from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS
-from pygeodesy.named import _Named, _NamedTuple, _Pass
+from pygeodesy.named import _Named, _NamedTuple, _Pass,  Property_RO
 from pygeodesy.points import areaOf, boundsOf, isconvex_, LatLon_
-# from pygeodesy.props import Property_RO  # from .fsums
+# from pygeodesy.props import Property_RO  # from .named
 from pygeodesy.units import Bool, FIx, HeightX, Lat, Lon, Number_
 
 # from math import fabs  # from .fmath
 
 __all__ = _ALL_LAZY.clipy
-__version__ = '23.05.18'
+__version__ = '23.09.20'
 
 _fj_       = 'fj'
 _original_ = 'original'
@@ -576,10 +576,10 @@ class _SH(_Named):
         x, dx = p1.lon, self._dx
         fy = float(p2.lat - y)
         fx = float(p2.lon - x)
-        d  = fy * dx - fx * dy
+        d  = fy * dx - fx * dy  # fdot((fx, fy), dx, -dy)
         if fabs(d) < EPS:  # PYCHOK no cover
             raise _AssertionError(self._DOT_(self.intersect.__name__))
-        d  = fsumf_(self._xy, -y * dx, x * dy) / d
+        d  = Fsum(self._xy, -y * dx, x * dy).fover(d)
         y += d * fy
         x += d * fx
         return _SHlli(y, x, p1.classof, edge)

pygeodesy/constants.py

@@ -13,9 +13,9 @@ from pygeodesy.basics import _copysign, isbool, iscomplex, isint,  _0_0
 from pygeodesy.errors import _xError, _xError2, _xkwds_get1, _xkwds_item2
 # from pygeodesy.internals import _0_0  # from .basics
 from pygeodesy.interns import _INF_, _NAN_, _UNDER_
-# from pygeodesy.lazily import _ALL_LAZY  # from .unitsBase
+from pygeodesy.lazily import _ALL_MODS as _MODS, _ALL_LAZY
 # from pygeodesy.streprs import Fmt  # from .unitsBase
-from pygeodesy.unitsBase import Float, Int, Radius,  _ALL_LAZY, Fmt
+from pygeodesy.unitsBase import Float, Int, Radius,  Fmt
 
 from math import fabs, isinf, isnan, pi as _pi, sqrt
 try:
@@ -24,7 +24,7 @@ except ImportError:  # Python 2-
     _inf, _nan = float(_INF_), float(_NAN_)
 
 __all__ = _ALL_LAZY.constants
-__version__ = '24.08.18'
+__version__ = '24.09.24'
 
 
 def _copysign_0_0(y):
@@ -100,7 +100,7 @@ def float0_(*xs):
         yield float(x) if x else _0_0
 
 
-def _float0(f):  # in .resections, .vector3dBase, ...
+def _float0(f):  # in .auxilats.auxily, .resections, .vector3dBase, ...
     '''(INTERNAL) Return C{float(B{f})} or C{INT0}.
     '''
     if f:
@@ -181,6 +181,7 @@ _10_0    = _float(  10)       # PYCHOK expected
 _16_0    = _float(  16)       # PYCHOK expected
 _32_0    = _float(  32)       # PYCHOK expected
 _60_0    = _float(  60)       # PYCHOK expected
+_64_0    = _float(  64)       # PYCHOK expected
 _90_0    = _float(  90)       # PYCHOK expected
 _100_0   = _float( 100)       # PYCHOK expected
 _180_0   = _float( 180)       # PYCHOK expected
@@ -224,30 +225,30 @@ except ImportError:  # PYCHOK no cover
     MIN_EXP  =  Int(MIN_EXP   =_log2(MIN))  # -307 base 10
 #   RADIX    =  Int(Radix     =2)           # base
 
-EPS0     = _Float( EPS0   = EPS**2)         # PYCHOK near-/non-zero comparison 4.930381e-32, or EPS or EPS_2
-EPS02    = _Float( EPS02  = EPS**4)         # PYCHOK near-zero-squared comparison 2.430865e-63
-EPS_2    = _Float( EPS_2  = EPS / _2_0)     # PYCHOK ≈ 1.110223024625e-16
-EPS1     = _Float( EPS1   =_1_0 - EPS)      # PYCHOK ≈ 0.9999999999999998
-EPS2     = _Float( EPS2   = EPS * _2_0)     # PYCHOK ≈ 4.440892098501e-16
-EPS4     = _Float( EPS4   = EPS * _4_0)     # PYCHOK ≈ 8.881784197001e-16
-# _1EPS  = _Float(_1EPS   =_1_0 + EPS)      # PYCHOK ≈ 1.0000000000000002
-_1_EPS   = _Float(_1_EPS  =_1_0 / EPS)      # PYCHOK = 4503599627370496.0
-# _2_EPS = _Float(_2_EPS  =_2_0 / EPS)      # PYCHOK = 9007199254740992.0
-_EPS2e4  = _Float(_EPS2e4 = EPS2 * 1.e4)    # PYCHOK ≈ 4.440892098501e-12
-_EPS4e8  = _Float(_EPS4e8 = EPS4 * 1.e8)    # PYCHOK ≈ 8.881784197001e-08
-_EPSmin  = _Float(_EPSmin = sqrt(MIN))      # PYCHOK = 1.49166814624e-154
-_EPSqrt  = _Float(_EPSqrt = sqrt(EPS))      # PYCHOK = 1.49011611938e5-08
-_EPStol  = _Float(_EPStol =_EPSqrt * _0_1)  # PYCHOK = 1.49011611938e5-09 == sqrt(EPS * _0_01)
-
-_89_999_ = _Float(_89_999_= EPS1 * _90_0)   # just below 90.0
+EPS0     = _Float( EPS0  = EPS**2)          # PYCHOK near-/non-zero comparison 4.930381e-32, or EPS or EPS_2
+EPS02    = _Float( EPS02 = EPS**4)          # PYCHOK near-zero-squared comparison 2.430865e-63
+EPS_2    = _Float( EPS_2 = EPS / _2_0)      # PYCHOK ≈ 1.110223024625e-16
+EPS1     = _Float( EPS1  =_1_0 - EPS)       # PYCHOK ≈ 0.9999999999999998
+EPS2     = _Float( EPS2  = EPS * _2_0)      # PYCHOK ≈ 4.440892098501e-16
+EPS4     = _Float( EPS4  = EPS * _4_0)      # PYCHOK ≈ 8.881784197001e-16
+# _1EPS  = _Float(_1EPS  =_1_0 + EPS)       # PYCHOK ≈ 1.0000000000000002
+_1_EPS   = _Float(_1_EPS =_1_0 / EPS)       # PYCHOK = 4503599627370496.0
+# _2_EPS = _Float(_2_EPS =_2_0 / EPS)       # PYCHOK = 9007199254740992.0
+_EPS2e4  = _Float(_EPS2e4= EPS2 * 1.e4)     # PYCHOK ≈ 4.440892098501e-12
+_EPS4e8  = _Float(_EPS4e8= EPS4 * 1.e8)     # PYCHOK ≈ 8.881784197001e-08
+_EPSjam  = _Float(_EPSjam= pow(EPS, 0.75))  # PYCHOK = 1.818989403546e-12
+_EPSmin  = _Float(_EPSmin= sqrt(MIN))       # PYCHOK = 1.49166814624e-154
+_EPSqrt  = _Float(_EPSqrt= sqrt(EPS))       # PYCHOK = 1.49011611938e5-08
+_EPStol  = _Float(_EPStol=_EPSqrt * _0_1)   # PYCHOK = 1.49011611938e5-09 == sqrt(EPS * _0_01)
+
+_89_999  = _Float(_89_999=_90_0 * EPS1)  # just below 90.0
 # <https://Numbers.Computation.Free.FR/Constants/Miscellaneous/digits.html>
-_1__90   = _Float(_1__90  =_1_0 / _90_0)    # PYCHOK = 0.011_111_111_111_111_111_111_111_111_111_111_111_111_111_111_11111
-_2__PI   = _Float(_2__PI  =_2_0 / _pi)      # PYCHOK = 0.636_619_772_367_581_343_075_535_053_490_057_448_137_838_582_96182
+_1__90   = _Float(_1__90 =_1_0 / _90_0)  # PYCHOK = 0.011_111_111_111_111_111_111_111_111_111_111_111_111_111_111_11111
+_2__PI   = _Float(_2__PI =_2_0 / _pi)    # PYCHOK = 0.636_619_772_367_581_343_075_535_053_490_057_448_137_838_582_96182
 
-_1_16th  = _Float(_1_16th =_1_0 / _16_0)  # PYCHOK in .ellipsoids, .karney
-_1_64th  = _Float(_1_64th =_1_0 /  64)    # PYCHOK in .elliptic, pow(2.0, -6)
-_1_3rd   = _Float(_1_3rd  =_1_0 /  _3_0)  # PYCHOK in .fmath
-_1_6th   = _Float(_1_6th  =_1_0 /  _6_0)  # PYCHOK in .fmath
+_1_16th  = _Float(_1_16th=_1_0 / _16_0)  # PYCHOK in .ellipsoids, .karney
+_1_3rd   = _Float(_1_3rd =_1_0 /  _3_0)  # PYCHOK in .fmath
+_1_6th   = _Float(_1_6th =_1_0 /  _6_0)  # PYCHOK in .fmath
 
 _K0_UTM  = _Float(_K0_UTM = 0.9996)  # PYCHOK in .etm, .ktm, .utm, UTM scale at central meridian
 # sqrt(2) <https://WikiPedia.org/wiki/Square_root_of_2>
@@ -319,20 +320,30 @@ def isclose(a, b, rel_tol=1e-12, abs_tol=EPS0):
     return _isclose(a, b, rel_tol=rel_tol, abs_tol=abs_tol)
 
 
+try:
+    from cmath import isfinite as _iscfinite
+except ImportError:  # Python 3.1-
+
+    def _iscfinite(x):  # PYCHOK not self?
+        '''Mimick Python 3.2+ C{cmath.isfinite}.
+        '''
+        return _isfinite(x.real) and _isfinite(x.imag)
+
 try:
     from math import isfinite as _isfinite  # in .ellipsoids, ...
 except ImportError:  # Python 3.1-
 
-    def _isfinite(x):
+    def _isfinite(x):  # PYCHOK not self?
         '''Mimick Python 3.2+ C{math.isfinite}.
         '''
         return not (isinf(x) or isnan(x))
 
 
 def isfinite(obj):
-    '''Check a finite C{scalar} or C{complex} value.
+    '''Check a finite C{scalar}, C{complex}, ... value.
 
-       @arg obj: Value (C{scalar} or C{complex}).
+       @arg obj: Value (C{scalar}, C{complex}, an L{Fsum} or
+                        L{Fsum2Tuple}).
 
        @return: C{False} if B{C{obj}} is C{INF}, C{NINF}
                 or C{NAN}, C{True} otherwise.
@@ -343,7 +354,9 @@ def isfinite(obj):
         return (obj not in _INF_NAN_NINF) and _isfinite(obj)
     except Exception as x:
         if iscomplex(obj):  # _isfinite(complex) thows TypeError
-            return isfinite(obj.real) and isfinite(obj.imag)
+            return _iscfinite(obj)
+        if _MODS.fsums._isFsumTuple(obj):  # OverflowError
+            return obj.is_finite()
         raise _xError(x, Fmt.PAREN(isfinite.__name__, obj))
 
 
@@ -441,7 +454,7 @@ def isnon0(x, eps0=EPS0):
 def _off90(lat):
     '''(INTERNAL) Off 90.0 for .gars and .wgrs.
     '''
-    return max(min(lat, _89_999_), -_89_999_)
+    return max(min(lat, _89_999), -_89_999)
 
 
 try:

pygeodesy/deprecated/__init__.py

@@ -27,7 +27,7 @@ __all__ = (_ALL_DEPRECATED.deprecated_bases +
            _ALL_DEPRECATED.deprecated_classes +
            _ALL_DEPRECATED.deprecated_consterns +
            _ALL_DEPRECATED.deprecated_functions)
-__version__ = '24.07.02'
+__version__ = '24.09.19'
 
 if _unLazy0:
     from pygeodesy.deprecated import bases, datum, nvector, rhumbBase, \

pygeodesy/deprecated/functions.py

@@ -14,7 +14,7 @@ from pygeodesy.props import deprecated_function
 from pygeodesy.units import Number_, Scalar_
 
 __all__ = _ALL_DEPRECATED.deprecated_functions
-__version__ = '24.06.11'
+__version__ = '24.09.19'
 
 _WGS84 = _UTM = object()
 
@@ -193,6 +193,14 @@ def fStrzs(floatstr):  # PYCHOK no cover
     return _MODS.streprs.fstrzs(floatstr)
 
 
+@deprecated_function
+def Fsum2product(*xs, **kwds):
+    '''DEPRECATED on 2024.09.19, use L{Fsum}C{(*B{xs}, B{f2product}=True, ...)}.'''
+    F = _MODS.fsums.Fsum(**kwds)
+    _ = F.f2product(True)
+    return F._facc_args(xs, up=False) if xs else F
+
+
 @deprecated_function
 def hypot3(x, y, z):  # PYCHOK no cover
     '''DEPRECATED, use function L{pygeodesy.hypot_}.'''