pygeodesy 24.3.24__py2.py3-none-any.whl → 24.4.4__py2.py3-none-any.whl

pygeodesy/lazily.py

@@ -284,7 +284,7 @@ _ALL_LAZY = _NamedEnum_RO(_name='_ALL_LAZY',
                                    'FrechetVincentys', 'Frechet6Tuple',
                                    'frechet_'),
                          fstats=_i('Fcook', 'Flinear', 'Fwelford'),
-                          fsums=_i('Fsum', 'Fsum2Tuple', 'ResidualError',
+                          fsums=_i('Fsum', 'DivMod2Tuple', 'Fsum2Tuple', 'ResidualError',
                                    'fsum', 'fsum_', 'fsumf_', 'fsum1', 'fsum1_', 'fsum1f_'),
                            gars=_i('Garef', 'GARSError'),
                       geodesicw=_i('Geodesic', 'GeodesicLine', 'Geodesic_WGS84'),
@@ -504,7 +504,7 @@ class _ALL_MODS(object):
 _ALL_MODS = _ALL_MODS()  # PYCHOK singleton
 
 __all__ = _ALL_LAZY.lazily
-__version__ = '24.02.22'
+__version__ = '24.03.29'
 
 
 def _ALL_OTHER(*objs):

pygeodesy/resections.py

@@ -13,9 +13,10 @@ L{triAngle}, L{triAngle5}, L{triSide}, L{triSide2} and L{triSide4}.
 from __future__ import division as _; del _  # PYCHOK semicolon
 
 from pygeodesy.basics import map1, map2, _zip,  _ALL_LAZY
-from pygeodesy.constants import EPS, EPS0, EPS02, INT0, NEG0, PI, PI2, PI_2, PI_4, \
-                               _0_0, _0_5, _1_0, _N_1_0, _2_0, _N_2_0, _4_0, _16_0, \
-                               _180_0, _360_0, isnear0, _over, _umod_360
+from pygeodesy.constants import EPS, EPS0, EPS02, INT0, PI, PI2, PI_2, PI_4, \
+                               _0_0, _0_5, _1_0, _N_1_0, _2_0, _N_2_0, _4_0, \
+                               _16_0, _180_0, _360_0, _copysign_0_0, isnear0, \
+                               _over, _umod_360
 from pygeodesy.errors import _and, _or, TriangleError, _ValueError, _xcallable, \
                              _xkwds, _xkwds_pop2
 from pygeodesy.fmath import favg, Fdot, fidw, fmean, hypot, hypot2_
@@ -33,7 +34,7 @@ from pygeodesy.vector3d import _otherV3d, Vector3d
 from math import cos, atan2, degrees, fabs, radians, sin, sqrt
 
 __all__ = _ALL_LAZY.resections
-__version__ = '24.03.24'
+__version__ = '24.04.04'
 
 _concyclic_ = 'concyclic'
 _PA_        = 'PA'
@@ -457,7 +458,7 @@ def pierlotx(point1, point2, point3, alpha1, alpha2, alpha3, useZ=False,
 
     def _cot(s, c):  # I{exact} cotangent
         try:
-            return (c / s) if c else (NEG0 if s < 0 else _0_0)
+            return (c / s) if c else _copysign_0_0(s)
         except ZeroDivisionError:
             raise ValueError(_or(_coincident_, _colinear_))
 
@@ -970,44 +971,43 @@ def wildberger3(a, b, c, alpha, beta, R3=min):
     def _s(x):
         return sin(x)**2
 
-    def _vpa(r1, r3, q2, q3, s3):
-        r = r1 * r3 * _4_0
-        n = (r - _F1(r1, r3, -q2).fpow(2)).fover(s3)
-        if n < 0 or isnear0(r):
+    def _vpa(r3, q2, q3, s2, s3):
+        r1 = s2 * q3 / s3
+        r  = r1 * r3 * _4_0
+        n  = (r - _F1(r1, r3, -q2)**2).fover(s3)
+        if n < 0 or r < EPS0:
             raise ValueError(_coincident_)
         return sqrt((n / r) * q3) if n else _0_0
 
     try:
-        a, b, c, da, db = t = map1(float, a, b, c, alpha, beta)
-        if min(t) < 0:
+        a, b, c, da, db = q = map1(float, a, b, c, alpha, beta)
+        if min(q) < 0:
             raise ValueError(_negative_)
 
-        ra, rb = radians(da), radians(db)
-        s1, s2, s3 = s = map1(_s, rb, ra, ra + rb)  # rb, ra!
-        if min(s) < EPS02:
-            raise ValueError(_or(_coincident_, _colinear_))
-
         q1, q2, q3 = q = a**2, b**2, c**2
         if min(q) < EPS02:
             raise ValueError(_coincident_)
 
-        r1 =  s2 * q3 / s3  # s2!
-        r2 =  s1 * q3 / s3  # s1!
-        Qs = _F1(*q)  # == hypot2_(a, b, c)
-        Ss = _F1(*s)  # == fsum1(s)
-        s += (Qs * _0_5),  # tuple!
-        C0 =  Fdot(s, q1, q2, q3, -Ss)
-        r3 =  C0.fover(-s3)
-        d0 =  Qs.fpow(2).fsub_(hypot2_(*q) * _2_0).fmul(s1 * s2).fover(s3)
+        ra, rb = map1(radians, da, db)
+        s1, s2, s3 = s = map1(_s, rb, ra, ra + rb)  # rb, ra!
+        if min(s) < EPS02:
+            raise ValueError(_or(_coincident_, _colinear_))
+
+        q4 =  hypot2_(*q) * _2_0  # a**4 + ...
+        Qs = _F1(*q)   # == hypot2_(a, b, c)
+        d0 = (Qs**2 - q4).fmul(s1 * s2).fover(s3)
+        if d0 < 0:
+            raise ValueError(_negative_)
+        s += _F1(*s),  # == fsum1(s),
+        C0 =  Fdot(s, q1, q2, q3, -Qs * _0_5)
+        r3 =  C0.fover(-s3)  # C0 /= -s3
         if d0 > EPS02:  # > c0
             _xcallable(R3=R3)
             d0 = sqrt(d0)
             r3 = R3(float(C0 + d0), float(C0 - d0))  # XXX min or max
-        elif d0 < 0:
-            raise ValueError(_negative_)
 
-        pa = _vpa(r1, r3, q2, q3, s3)
-        pb = _vpa(r2, r3, q1, q3, s3)
+        pa = _vpa(r3, q2, q3, s2, s3)
+        pb = _vpa(r3, q1, q3, s1, s3)
         pc =  favg(_triSide2(b, pa, ra).a,
                    _triSide2(a, pb, rb).a)
         return Survey3Tuple(pa, pb, pc, name=wildberger3.__name__)

pygeodesy/utily.py

@@ -10,12 +10,14 @@ U{Vector-based geodesy<https://www.Movable-Type.co.UK/scripts/latlong-vectors.ht
 # make sure int/int division yields float quotient, see .basics
 from __future__ import division as _; del _  # PYCHOK semicolon
 
-from pygeodesy.basics import _copysign, isinstanceof, isint, isstr, neg, _passargs
-from pygeodesy.constants import EPS, EPS0, INF, NAN, NEG0, PI, PI2, PI_2, R_M, \
-                               _float as _F, _isfinite, isnan, isnear0, _over, \
-                               _umod_360, _umod_PI2, _M_KM, _M_NM, _M_SM, _0_0, \
-                               _1__90, _0_5, _1_0, _N_1_0, _2__PI, _10_0, _90_0, \
-                               _N_90_0, _180_0, _N_180_0, _360_0, _400_0
+from pygeodesy.basics import _copysign, isinstanceof, isint, isstr, neg, \
+                             _passargs
+from pygeodesy.constants import EPS, EPS0, INF, NAN, PI, PI2, PI_2, R_M, \
+                               _M_KM, _M_NM, _M_SM, _0_0, _1__90, _0_5, _1_0, \
+                               _N_1_0, _2__PI, _10_0, _90_0, _N_90_0, _180_0, \
+                               _N_180_0, _360_0, _400_0, _copysign_0_0, \
+                               _float as _F, _isfinite, isnan, isnear0, \
+                               _over, _umod_360, _umod_PI2
 from pygeodesy.errors import _ValueError, _xkwds, _xkwds_get,  _ALL_LAZY, _MODS
 from pygeodesy.interns import _edge_, _radians_, _semi_circular_, _SPACE_
 # from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS  # from .errors
@@ -25,7 +27,7 @@ from pygeodesy.units import Degrees, Degrees_, Feet, Float, Lam, Lam_, \
 from math import acos, asin, atan2, cos, degrees, fabs, radians, sin, tan  # pow
 
 __all__ = _ALL_LAZY.utily
-__version__ = '24.01.12'
+__version__ = '24.04.04'
 
 # read constant name "_M_Unit" as "meter per Unit"
 _M_CHAIN     = _F(  20.1168)     # yard2m(1) * 22
@@ -683,7 +685,7 @@ def _sin0cos2(q, r, sign):
 #   else:  # r == 0, testUtility failures
 #       t = _0_0, _1_0, _0_0, _N_1_0, _0_0
 #   q &= 3
-    s = t[q]     or (NEG0 if sign < 0 else _0_0)
+    s = t[q]     or _copysign_0_0(sign)
     c = t[q + 1] or _0_0
     return s, c
 
@@ -691,8 +693,7 @@ def _sin0cos2(q, r, sign):
 def SinCos2(x):
     '''Get C{sin} and C{cos} of I{typed} angle.
 
-       @arg x: Angle (L{Degrees}, L{Degrees_}, L{Radians}, L{Radians_}
-               or scalar C{radians}).
+       @arg x: Angle (L{Degrees}, L{Radians} or scalar C{radians}).
 
        @return: 2-Tuple (C{sin(B{x})}, C{cos(B{x})}).
     '''