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

pygeodesy/internals.py

@@ -132,7 +132,7 @@ class _MODS_Base(object):
     def name(self):
         '''Get this name (C{str}).
         '''
-        return typename(self.__class__)
+        return typename(type(self))
 
     @_Property_RO
     def nix2(self):  # PYCHOK no cover
@@ -686,7 +686,7 @@ def _versions(sep=_SPACE_):
 
 
 __all__ = tuple(map(typename, (machine, print_, printf, typename)))
-__version__ = '25.04.14'
+__version__ = '25.05.26'
 
 if __name__ == _DMAIN_:
 

pygeodesy/interns.py

@@ -443,15 +443,15 @@ _LR_PAIRS  = {_LANGLE_:  _RANGLE_,
 
 __all__ = (_NN_,  # NOT MISSING!
             Str_.__name__)  # classes
-__version__ = '25.04.12'
+__version__ = '25.05.26'
 
 if __name__ == _DMAIN_:
 
-    def _main():
+    def _main(globalocals):
         from pygeodesy import itemsorted, printf
 
         t = b = 0
-        for n, v in itemsorted(locals(), asorted=False, reverse=True):
+        for n, v in itemsorted(globalocals, asorted=False, reverse=True):
             if n.endswith(_UNDER_) and n.startswith(_UNDER_) and \
                                    not n.startswith(_DUNDER_):
                 t += 1
@@ -459,10 +459,10 @@ if __name__ == _DMAIN_:
                 m  = n[1:-1]
                 if m != v and m.replace(_UNDER_, _SPACE_) != v:
                     printf('%4d: %s = %r', t, n, v)
-        n = len(locals())
+        n = len(globalocals)
         printf('%4d (%d) names, %s chars total, %.2f chars avg', t, n, b, float(b) / t, nl=1)
 
-    _main()
+    _main(globals())  # or locals()
 
 # **) MIT License
 #

pygeodesy/karney.py

@@ -166,7 +166,7 @@ from pygeodesy.utily import atan2d, sincos2d, tand, _unrollon,  fabs
 # from math import fabs  # from .utily
 
 __all__ = _ALL_LAZY.karney
-__version__ = '25.05.12'
+__version__ = '25.05.28'
 
 _2_4_       = '2.4'
 _K_2_0      = _getenv(_PYGEODESY_ENV(typename(_xgeographiclib)[2:]), _2_)
@@ -240,6 +240,8 @@ class Caps(object):
 
        C{LONG_UNROLL} - unroll C{lon2} in method C{.Direct},
 
+       C{NONFINITONAN} - set all C{GDict} items to C{NAN} iff any argument is C{non-finite}.
+
        C{REDUCEDLENGTH} - compute the reduced length C{m12},
 
        C{REVERSE2} - reverse azimuth C{azi2} by 180 degrees,
@@ -274,7 +276,7 @@ class Caps(object):
 
     LINE_CAPS     =  STANDARD_LINE | REDUCEDLENGTH | GEODESICSCALE  # .geodesici only
     LONG_UNROLL   =  1 << 15  # unroll C{lon2} in .Direct and .Position
-#                 =  1 << 16  # unused
+    NONFINITONAN  =  1 << 16  # see method GDict._toNAN
     LINE_OFF      =  1 << 17  # Line without updates from parent geodesic or rhumb
     REVERSE2      =  1 << 18  # reverse C{azi2}
     ALL           =  0x7F80 | _CAP_ALL  # without LONG_UNROLL, LINE_OFF, REVERSE2 and _DEBUG_*
@@ -327,9 +329,12 @@ if _FOR_DOCS:  # PYCHOK force ...
 else:
     Caps = Caps()  # PYCHOK singleton
 
-_key2Caps = dict(a12 =Caps.DISTANCE,  # in GDict._unCaps
+_key2Caps = dict(a12 =Caps.DISTANCE,  # in GDict._toNAN, -._unCaps
+                 azi1=Caps.AZIMUTH,
                  azi2=Caps.AZIMUTH,
+                 lat1=Caps.LATITUDE,
                  lat2=Caps.LATITUDE,
+                 lon1=Caps.LONGITUDE,
                  lon2=Caps.LONGITUDE,
                  m12 =Caps.REDUCEDLENGTH,
                  M12 =Caps.GEODESICSCALE,
@@ -353,6 +358,7 @@ class _CapsBase(_NamedBase):  # in .auxilats, .geodesicx.gxbases
     LINE_OFF      = Caps.LINE_OFF
     LONGITUDE     = Caps.LONGITUDE
     LONG_UNROLL   = Caps.LONG_UNROLL
+    NONFINITONAN  = Caps.NONFINITONAN
     REDUCEDLENGTH = Caps.REDUCEDLENGTH
     STANDARD      = Caps.STANDARD
     STANDARD_LINE = Caps.STANDARD_LINE  # for geodesici
@@ -449,11 +455,14 @@ class GDict(ADict):  # XXX _NamedDict
         '''
         return self._toTuple(Inverse10Tuple, dflt)
 
-    def _toNAN(self, outmask):  # .GeodesicLineExact._GenPosition
+    def _toNAN(self, outmask):  # .GeodesicExact._GDistInverse, .GeodesicLineExact._GenPosition
         '''(INTERNAL) Convert this C{GDict} to all C{NAN}s.
         '''
-        d = dict((n, NAN) for n in GeodSolve12Tuple._Names_)
-        return self.set_(**d)._unCaps(outmask)
+        if (outmask & Caps.NONFINITONAN):
+            d = dict((k, NAN) for k, C in _key2Caps.items()
+                                       if (outmask & C) == C)
+            self.set_(**d)
+        return self
 
     @deprecated_method
     def toRhumb7Tuple(self, dflt=NAN):  # PYCHOK no cover
@@ -543,7 +552,7 @@ class Inverse10Tuple(_GTuple):
                                    **updates)  # PYCHOK indent
 
 
-class _kWrapped(object):  # in .geodesicw
+class _kWrapped(_CapsBase):  # in .geodesicw
     '''(INTERNAL) Wrapper for some of I{Karney}'s U{geographiclib
         <https://PyPI.org/project/geographiclib>} classes.
     '''
@@ -554,13 +563,15 @@ class _kWrapped(object):  # in .geodesicw
            <https://PyPI.org/project/geographiclib>} package is installed,
            otherwise raise a C{LazyImportError}.
         '''
-        g = _xgeographiclib(type(self).__module__, 1, 49)
+        g = _xgeographiclib(type(self), 1, 49)
         from geographiclib.geodesic import Geodesic
         g.Geodesic = Geodesic
         from geographiclib.geodesicline import GeodesicLine
         g.GeodesicLine = GeodesicLine
         from geographiclib.geomath import Math
         g.Math = Math
+#       from geographiclib.polygonarea import PolygonArea
+#       g.PolygonArea = PolygonArea  # see below
         return g
 
     @property_ROnce
@@ -586,6 +597,11 @@ class _kWrapped(object):  # in .geodesicw
         return (_2_4_ if _K_2_4 else
                (_2_   if _K_2_0 else _1_)) if self.Math else NN
 
+    @property_ROnce
+    def _PolygonArea(self):  # lazy import
+        from geographiclib.polygonarea import PolygonArea
+        return PolygonArea
+
 _wrapped = _kWrapped()  # PYCHOK singleton, .datum, .test/base.py
 
 
@@ -830,7 +846,7 @@ def _norm180(deg):  # mimick geomath.Math.AngNormalize
         return d
 
 
-def _polygon(geodesic, points, closed, line, wrap):
+def _polygon(geodesic, points, closed, line, wrap, polar):
     '''(INTERNAL) Compute the area or perimeter of a polygon,
         using a L{GeodesicExact}, L{GeodesicSolve} or (if the
         C{geographiclib} package is installed) a C{Geodesic}
@@ -848,7 +864,7 @@ def _polygon(geodesic, points, closed, line, wrap):
         if not closed:  # closed only
             raise _ValueError(closed=closed, points=_composite_)
 
-        return points._sum1(_a_p, closed, line, wrap)
+        return points._sum1(_a_p, closed, line, wrap, polar)
 
     gP = geodesic.Polygon(line)
     _A = gP.AddPoint
@@ -867,7 +883,7 @@ def _polygon(geodesic, points, closed, line, wrap):
         _A(p0.lat, p0.lon)
 
     # gP.Compute returns (number_of_points, perimeter, signed area)
-    return gP.Compute(False, True)[1 if line else 2]
+    return gP.Compute(reverse=False, sign=True, polar=polar)[1 if line else 2]
 
 
 try:

pygeodesy/lcc.py

@@ -50,7 +50,7 @@ from pygeodesy.utily import atan1, degrees90, degrees180, sincos2, tanPI_2_2
 from math import atan, fabs, log, radians, sin, sqrt
 
 __all__ = _ALL_LAZY.lcc
-__version__ = '25.05.12'
+__version__ = '25.05.26'
 
 _E0_   = 'E0'
 _N0_   = 'N0'
@@ -457,7 +457,7 @@ class Lcc(_NamedBase):
 #                                  h=self.height, conic=self.conic,
 #                                  name=self._name__(name))
 #       args, kwds = _args_kwds(**kwds)
-#       return self.__class__(*args, **kwds)  # .classof
+#       return type(self)(*args, **kwds)  # .classof
 
     @Property_RO
     def easting(self):

pygeodesy/props.py

@@ -26,7 +26,7 @@ from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS, \
 from functools import wraps as _wraps
 
 __all__ = _ALL_LAZY.props
-__version__ = '25.05.21'
+__version__ = '25.05.26'
 
 _class_       = 'class'
 _DNL_         = _NL_ * 2  # PYCHOK used!
@@ -347,7 +347,7 @@ class property_RO(_PropertyBase):
             if Clas:  # overrides inst.__class__
                 d = Clas[0].__dict__.get(uname, MISSING)
             else:
-                d = getattr(inst.__class__, uname, MISSING)
+                d = getattr(type(inst), uname, MISSING)
 #               if d is MISSING:  # XXX superfluous
 #                   for c in type(inst).__mro__[:-1]:
 #                       if uname in c.__dict__:
@@ -403,7 +403,7 @@ class property_ROver(_property_RO___):
         '''
         v = self.method(inst)
         n = self.name
-        C = inst.__class__
+        C = type(inst)
         for c in C.__mro__:  # [:-1]
             if getattr(c, n, None) is self:
                 setattr(c, n, v)  # overwrite property_ROver

pygeodesy/sphericalBase.py

@@ -40,7 +40,7 @@ from pygeodesy.utily import acos1, asin1, atan2b, atan2d, degrees90, \
 from math import cos, fabs, log, sin, sqrt
 
 __all__ = _ALL_LAZY.sphericalBase
-__version__ = '25.05.12'
+__version__ = '25.05.26'
 
 
 class CartesianSphericalBase(CartesianBase):
@@ -592,7 +592,7 @@ def _intersecant2(c, r, p, b, radius=R_M, exact=False, height=None, wrap=False):
         p = _unrollon(c, p, wrap=wrap)
     nonexact = exact is None
 
-    if not isinstanceof(r, c.__class__, p.__class__):
+    if not isinstanceof(r, type(c), type(p)):
         r = Radius_(circle=r)
     elif nonexact:
         r = c.distanceTo(r, radius=radius, wrap=wrap)
@@ -601,7 +601,7 @@ def _intersecant2(c, r, p, b, radius=R_M, exact=False, height=None, wrap=False):
     else:
         raise _ValueError(exact=exact)
 
-    if not isinstanceof(b, c.__class__, p.__class__):
+    if not isinstanceof(b, type(c), type(p)):
         b = Bearing(b)
     elif nonexact:
         b = p.initialBearingTo(b, wrap=wrap)

pygeodesy/sphericalNvector.py

@@ -62,7 +62,7 @@ from pygeodesy.utily import atan2, degrees360, sincos2, sincos2_, sincos2d, \
 # from math import fabs  # from utily
 
 __all__ = _ALL_LAZY.sphericalNvector
-__version__ = '25.05.12'
+__version__ = '25.05.27'
 
 _lines_ = 'lines'
 
@@ -782,18 +782,19 @@ class Nvector(NvectorBase):
 _Nv00 = LatLon(_0_0, _0_0, name=_Nv00_)  # reference instance (L{LatLon})
 
 
-def areaOf(points, radius=R_M, wrap=False):
+def areaOf(points, radius=R_M, wrap=False, polar=False):
     '''Calculate the area of a (spherical) polygon or composite (with
        great circle arcs joining consecutive points).
 
-       @arg points: The polygon points or clips (C{LatLon}[],
-                    L{BooleanFHP} or L{BooleanGH}).
+       @arg points: The polygon points or clips (C{LatLon}[], L{BooleanFHP}
+                    or L{BooleanGH}).
        @kwarg radius: Mean earth radius (C{meter}) or C{None}.
-       @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 of C{B{polar}=True} not implemented (C{bool}).
 
-       @return: Polygon area (C{meter} I{squared}, same units as
-                B{C{radius}}, or C{radians} if C{B{radius} is None}).
+       @return: Polygon area (C{meter} I{squared}, same units as B{C{radius}},
+                or C{radians} if C{B{radius} is None}).
 
        @raise PointsError: Insufficient number of B{C{points}}.
 
@@ -802,6 +803,9 @@ def areaOf(points, radius=R_M, wrap=False):
        @see: Functions L{pygeodesy.areaOf}, L{sphericalTrigonometry.areaOf}
              and L{ellipsoidalKarney.areaOf}.
     '''
+    if polar:
+        notImplemented(None, polar=polar, up=2)
+
     def _interangles(ps, w):  # like .karney._polygon
         Ps = _Nv00.PointsIter(ps, loop=2, wrap=w)
         # use vector to 1st point as plane normal for sign of α
@@ -824,15 +828,14 @@ def areaOf(points, radius=R_M, wrap=False):
     if _MODS.booleans.isBoolean(points):
         r = points._sum2(LatLon, areaOf, radius=None, wrap=wrap)
     else:
-        # sum interior angles: depending on whether polygon is cw or ccw,
-        # angle between edges is π−α or π+α, where α is angle between
-        # great-circle vectors; so sum α, then take n·π − |Σα| (cannot
-        # use Σ(π−|α|) as concave polygons would fail)
+        # sum interior angles: depending on whether polygon is cw or
+        # ccw, angle between edges is π−α or π+α, where α is angle
+        # between great-circle vectors; so sum α, then take n·π − |Σα|
+        # (cannot use Σ(π−|α|) as concave polygons would fail)
         s = fsum(_interangles(points, wrap))
         # using Girard’s theorem: A = [Σθᵢ − (n−2)·π]·R²
-        # (PI2 - abs(s) == (n*PI - abs(s)) - (n-2)*PI)
-        r = fabs(PI2 - fabs(s))
-    return r if radius is None else (r * Radius(radius)**2)
+        r = fabs(PI2 - fabs(s))  # == n*PI - abs(s) - (n-2)*PI
+    return r if radius is None else (Radius(radius)**2 * r)
 
 
 def intersecant2(center, circle, point, other, **radius_exact_height_wrap):

pygeodesy/sphericalTrigonometry.py

@@ -57,7 +57,7 @@ from pygeodesy.vector3d import sumOf, Vector3d
 from math import asin, cos, degrees, fabs, radians, sin
 
 __all__ = _ALL_LAZY.sphericalTrigonometry
-__version__ = '25.05.12'
+__version__ = '25.05.28'
 
 _PI_EPS4 = PI - EPS4
 if _PI_EPS4 >= PI:
@@ -758,20 +758,21 @@ class LatLon(LatLonSphericalBase):
 _T00 = LatLon(0, 0, name='T00')  # reference instance (L{LatLon})
 
 
-def areaOf(points, radius=R_M, wrap=False):  # was=True
-    '''Calculate the area of a (spherical) polygon or composite (with the
-       points joined by great circle arcs).
+def areaOf(points, radius=R_M, wrap=False, polar=False):  # was wrap=True
+    '''Calculate the area of a (spherical) polygon or composite (with the points joined by
+       great circle arcs).
 
-       @arg points: The polygon points or clips (L{LatLon}[], L{BooleanFHP}
-                    or L{BooleanGH}).
-       @kwarg radius: Mean earth radius, ellipsoid or datum (C{meter},
-                      L{Ellipsoid}, L{Ellipsoid2}, L{Datum} or L{a_f2Tuple})
-                      or C{None}.
-       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the B{C{points}}
-                    (C{bool}).
+       @arg points: The polygon points or clips (L{LatLon}[], L{BooleanFHP} or L{BooleanGH}).
+       @kwarg radius: Mean earth radius, ellipsoid or datum (C{meter}, L{Ellipsoid},
+                      L{Ellipsoid2}, L{Datum} or L{a_f2Tuple}) or C{None}.
+       @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: Polygon area (C{meter} I{quared}, same units as B{C{radius}}
-                or C{radians} if C{B{radius} is None}).
+       @return: Polygon area (C{meter} I{quared}, same units as B{C{radius}} or C{radians} if
+                C{B{radius} is None}).
 
        @raise PointsError: Insufficient number of B{C{points}}.
 
@@ -779,13 +780,11 @@ def areaOf(points, radius=R_M, wrap=False):  # was=True
 
        @raise ValueError: Invalid B{C{radius}} or semi-circular polygon edge.
 
-       @note: The area is based on I{Karney}'s U{'Area of a spherical
-              polygon'<https://MathOverflow.net/questions/97711/
-              the-area-of-spherical-polygons>}, 3rd Answer.
+       @note: The area is based on I{Karney}'s U{'Area of a spherical polygon'
+              <https://MathOverflow.net/questions/97711/ the-area-of-spherical-polygons>}, 3rd Answer.
 
-       @see: Functions L{pygeodesy.areaOf}, L{sphericalNvector.areaOf},
-             L{pygeodesy.excessKarney}, L{ellipsoidalExact.areaOf} and
-             L{ellipsoidalKarney.areaOf}.
+       @see: Functions L{pygeodesy.areaOf}, L{sphericalNvector.areaOf}, L{pygeodesy.excessKarney},
+             L{ellipsoidalExact.areaOf} and L{ellipsoidalKarney.areaOf}.
     '''
     if _MODS.booleans.isBoolean(points):
         return points._sum2(LatLon, areaOf, radius=radius, wrap=wrap)
@@ -822,7 +821,7 @@ def areaOf(points, radius=R_M, wrap=False):  # was=True
             p1, z1 = p2, z2
 
     R = abs(R * _2_0)
-    if  abs(D) < _90_0:  # ispolar(points)
+    if  abs(D) < _90_0 or polar:  # ispolar(points)
         R = abs(R - PI2)
     if radius:
         a  =  degrees(A.fover(len(A)))  # mean lat

pygeodesy/vector2d.py

@@ -36,7 +36,7 @@ from contextlib import contextmanager
 # from math import fabs, sqrt  # from .fmath
 
 __all__ = _ALL_LAZY.vector2d
-__version__ = '25.05.09'
+__version__ = '25.05.26'
 
 _cA_        = 'cA'
 _cB_        = 'cB'
@@ -437,7 +437,7 @@ class _numpy(object):  # see also .formy._idllmn6, .geodesicw._wargs, .latlonBas
     def np(self):
         '''Import numpy 1.10+ once.
         '''
-        return _xnumpy(self.__class__, 1, 10)
+        return _xnumpy(type(self), 1, 10)
 
     def null_space2(self, A, rcond=None):
         '''Return the C{null_space} and C{rank} of matrix B{C{A}}.

{pygeodesy-25.5.25.dist-info → pygeodesy-25.5.28.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: pygeodesy
-Version: 25.5.25
+Version: 25.5.28
 Summary: Pure Python geodesy tools
 Home-page: https://GitHub.com/mrJean1/PyGeodesy
 Author: Jean M. Brouwers
@@ -164,7 +164,7 @@ PyCodeStyle_ (formerly Pep8) and McCabe_ using Python 2.7.18, both in 64-bit on
 
 For a summary of all *Karney*-based functionality in ``pygeodesy``, see module karney_.
 
-*Last updated: May 25, 2025.*
+*Last updated: May 28, 2025.*
 
 License
 =======

{pygeodesy-25.5.25.dist-info → pygeodesy-25.5.28.dist-info}/RECORD

@@ -1,14 +1,14 @@
 pygeodesy/LICENSE,sha256=YfgAiyxOwY6P9Kkb1_5XN81nueTLrpb3Ffkv3EuPgFU,1144
-pygeodesy/__init__.py,sha256=lJoIhQqDE5sqmoU2P-S8CbYP2AR1CEMSqZvNu9ulO4U,42290
-pygeodesy/__main__.py,sha256=B7CImXa4fm8lEqjTTp5AEalqIhldNpghdHtJkyMdyh0,5522
+pygeodesy/__init__.py,sha256=C2E8VxEx5MOANjeIZ12lgXtsJSbJdYh1Ee05IkobA1w,42290
+pygeodesy/__main__.py,sha256=HpOBmuhI7yVP1FIAEugvd9Lb54kzew-nM78be0r4SC4,5597
 pygeodesy/albers.py,sha256=mb3YbVvoBq8a7AytT0HeVxe8DGEvx1KFN2Usl2ksKwk,30908
 pygeodesy/azimuthal.py,sha256=acwIA1B7jmg1oDpJBuFitnOVHUWRUk7Y73yI8zqIGrc,49854
 pygeodesy/basics.py,sha256=ErpbLgHNerGPKiOIlWKdKzq6lO0dhoZjaywTnCn8PyE,32905
-pygeodesy/booleans.py,sha256=ncA4SV-2TX78CUsRKRgFh43g5NhkW_Iy1mnnudA3Nog,73683
+pygeodesy/booleans.py,sha256=njAJlAoVZEcQ5jp8wwyuVo8-jFPE3PFpAF3k3dQihBM,73683
 pygeodesy/cartesianBase.py,sha256=tvb1pGUTvSkIAInNhrnHXjBCwnQrVn6a-H2fipoapKE,45365
 pygeodesy/clipy.py,sha256=a0azaK7YDmGDU8mVVKe2MvUiPEhqZ9XJ-9u0T1iHKbo,27779
-pygeodesy/constants.py,sha256=xU8_niBOG4VFdWEqCi4XARiNJEAmprSTPCNkUrG3Bec,19899
-pygeodesy/css.py,sha256=Vn4MSpCagj55L_SM47AVX6V5sPglG17X1JJiU1rW1yA,25797
+pygeodesy/constants.py,sha256=usdHWtYiy8iaPXxXau2vsViL9zC8I4Jp2fziiEmp2g4,19947
+pygeodesy/css.py,sha256=vdQax5m7BMSpr7PUyqedCSvx73ilWuhXQTp39i8FOfA,25793
 pygeodesy/datums.py,sha256=Vp1fP7sSV1EuRbuO2scnZ3qiTZD1CvL4O12KSRlf1SU,34068
 pygeodesy/dms.py,sha256=CKXt2EOYFBluQV5q6cFYppJl3dCvvEB_0E84J0ki6ws,42202
 pygeodesy/ecef.py,sha256=DcWNhnOZC7-Olns9c1XI9p4MmuC00KN5rZ5W1UQVUxw,59009
@@ -16,37 +16,37 @@ pygeodesy/ecefLocals.py,sha256=zoxGUJeYX2UXzR9r0jaByrvfrrpmrvAhjyQsz433RjA,7794
 pygeodesy/elevations.py,sha256=WrWPZNgkZX9Ys_1ij_gyE4jzdIv85JYAaR2cK3vc4gQ,10879
 pygeodesy/ellipsoidalBase.py,sha256=n6VHIaRZ6UvC1H36qolKxwlHJX5BZysvgE3iir-PunE,55743
 pygeodesy/ellipsoidalBaseDI.py,sha256=PbSJkUUHpRGYgnnciMD-sWqTCReZv_FAc8SUonF3ZbI,38544
-pygeodesy/ellipsoidalExact.py,sha256=YVGVUmBTYu4NuYZdIx2f6zCCXM6Tqbv2nlkXQBqhvus,16973
-pygeodesy/ellipsoidalGeodSolve.py,sha256=wHv9UHYqnoU_PH8zgRJ-kA1-o-7zZLs1-O_zZHSQb1U,16666
-pygeodesy/ellipsoidalKarney.py,sha256=5RfQzNJFq4VklDnr3TRo6kh29rcQ4D7XkXUr_2Zqf38,18721
+pygeodesy/ellipsoidalExact.py,sha256=kKnQV3gz293dPowsRx5aI0CJ05a2TWGpEcWoJCxut30,17111
+pygeodesy/ellipsoidalGeodSolve.py,sha256=w5TubUdqYwvrtxfIwP2Scr1FzISq2szx1HX88TUa42E,17107
+pygeodesy/ellipsoidalKarney.py,sha256=W_Liw5MxTqJJGsxMi50ui8grRSPKwkG2HH-Mv_mXQYk,18839
 pygeodesy/ellipsoidalNvector.py,sha256=HLFjX6o26q2fzVCjLSL0PdEh96xhGtjfzo2V_k54j4A,30146
-pygeodesy/ellipsoidalVincenty.py,sha256=5syxUVS5RbTz9ZATU4ViegJPFBpdZWt6i8ZGVInype4,26122
+pygeodesy/ellipsoidalVincenty.py,sha256=zeTvid09NGuYJ1WJ__a3ApCBGSzX0A-Rhs8-sQXMNQQ,26273
 pygeodesy/ellipsoids.py,sha256=v964yZIL2nCc8MJgt3nzV6pczaQsQubXmX_sAh_nrBw,107765
-pygeodesy/elliptic.py,sha256=JkZzRytShhZsLnS-t031UA49DZnKXjob6ocGQZKI-GA,45020
+pygeodesy/elliptic.py,sha256=nob5Iby3RIj1tnsZJ1AzzxDLwPJm48M9Cf_kIAwU6l8,45035
 pygeodesy/epsg.py,sha256=VIBidzAEuoBF4LdlJ4r3WtpL0Ymhru8kf4fAgfFr5ok,8220
 pygeodesy/errors.py,sha256=OUEJ51M3c-JH1n0BBbRtPTCDX0I00rH_BQPYSRoKrqo,32043
 pygeodesy/etm.py,sha256=NDtpmbQmbGieV64iJ-EDIAvD9ryjrW0hiD3wBBtxULk,46744
 pygeodesy/fmath.py,sha256=sIN7mW2KOiK0yU93mvmHo1OOTkGxKvwRinvyr_Qdf6I,37689
 pygeodesy/formy.py,sha256=7-CdPiagdVrbDMjWieX_guN4RvFVCYKOM3yDKSFPRWU,69777
-pygeodesy/frechet.py,sha256=k-6pWb3WQAXXIUdAzG2AnFxvnmK7k-GHj4CoYdzDNbg,35592
+pygeodesy/frechet.py,sha256=dBnQ5ePuaUdKv62-e9RfovFDxELsZPQEKMLww9NOf7Q,35588
 pygeodesy/fstats.py,sha256=M6UOrwmeGp6kjYuYfc6Iq9ibMJAFNm6D4C5GnvlJeuM,28348
 pygeodesy/fsums.py,sha256=3pXl8rfDPeKzXHEhtf9qNzBtUme4ckskDx3HWfJbfyg,103988
 pygeodesy/gars.py,sha256=AXioT4Lh1cwrd7fF6dHcVJw-9vu2oAnZrzo_ka2Y3Ts,11761
 pygeodesy/geodesici.py,sha256=wPUDlRPOcj3MvvWA0CtnZvC30t1onKyZYvR51iCjgks,74891
-pygeodesy/geodesicw.py,sha256=dafScZ4nVuC_vlHXuYlSan32CsIEiL9uNcJYfnY1NtQ,29004
+pygeodesy/geodesicw.py,sha256=xDxxxr4HOT-8Mh85zD6QZzyO-mfqUavfwQMStYrsF8c,30029
 pygeodesy/geodsolve.py,sha256=fFlyNlkEOpkJaeddEuKwfhtqaZjmVs0egNGM5cwFFoc,29101
 pygeodesy/geohash.py,sha256=1P9-k9sMN-ZVFr2iWbQpFi1lg3uNmChHFCIfETkgFJU,40153
-pygeodesy/geoids.py,sha256=Uke9Lvb0QzMicWyE_H6dwYoCGhby3-1rCAQ-TIoCSgo,86129
-pygeodesy/hausdorff.py,sha256=3iDV4FetRrHcz2D1ddswx56XNi18zzCE2HMTDDoK6sI,32273
-pygeodesy/heights.py,sha256=FfBP1-Zw9XCjp3sBrTqtgZr_KKeejne2cizW0Lhkg_U,41014
-pygeodesy/internals.py,sha256=0jfelLhG50VY0hgy2qzqUvEFDJrcLDfpI_eNs4Go6BI,24028
-pygeodesy/interns.py,sha256=IBM_RuZDmP__4TUx9MLCU3L-Vx5uZz0VBi2eoH0hV44,23467
+pygeodesy/geoids.py,sha256=OQvOVutr8R1CalStgxmNu2Mk3oMl4zpCY_OKh0KHpI8,86125
+pygeodesy/hausdorff.py,sha256=KGB0CIeNkBqR0sXIbRJvcgXPrM-4Mpi2aBOqDeZjqZ0,32269
+pygeodesy/heights.py,sha256=_Pn-6w4gGzdXCoA6m72rJm4hb-Fe2vxc8RiRRjq5iTI,41002
+pygeodesy/internals.py,sha256=y3qT-ZbNDc-7md_T87qtVk1x2rDVtJbP5G0P-HgBlVA,24024
+pygeodesy/interns.py,sha256=QTP_6rStQ283DE-rL7rCAcfEffuiAXustBGkCY_DHJ0,23508
 pygeodesy/iters.py,sha256=nKVp_LUQyTqOz_OLzY25g6-_xlKG7lmAnjKKMXESmiM,20345
-pygeodesy/karney.py,sha256=xlkRlmg66XRdL1ale2sWHhxfostQU86II4qwNyk-ewM,40938
+pygeodesy/karney.py,sha256=IzWyGwX0kRpwrVFBpCzogd2Q8PLvCdCwHheE-zbCQZE,41629
 pygeodesy/ktm.py,sha256=3CLwbRvk-sr4qs8R9fdJ_ejokgZEjKHULV7wzzuonHU,27259
 pygeodesy/latlonBase.py,sha256=B7wNE-IJf-Sm3NKhCXdrgwRbH_4dBGluZlQdxTJqp8s,75320
 pygeodesy/lazily.py,sha256=2FzJkIQK7LftAXXYxMbqTQUAWth7T3syIhqm8Z9bhuo,47609
-pygeodesy/lcc.py,sha256=N31gphU0nSrc-i2CeJYux7CgbHrTfHwAXajta-mV57s,25790
+pygeodesy/lcc.py,sha256=7Di2Y5VIDL5drdSce4ABWGT3MtZSosaLKd68US6788s,25786
 pygeodesy/ltp.py,sha256=N3TMipTWKUOeZYzn3UeJJqe5J4Yk788390ohQM1ykAE,50893
 pygeodesy/ltpTuples.py,sha256=t-n7p9A3JCVasv8LKCV4PkOTmXJZi91BkBFm8pBZM5o,59101
 pygeodesy/mgrs.py,sha256=ZgCBH_M_YFxuFu5ukVOVetk-uftv-4GJP-QJrtoqtcU,29842
@@ -55,13 +55,13 @@ pygeodesy/namedTuples.py,sha256=JlO8uWmR65RkT7Rdgc9knfI0-bqJeA-tm2ZKsesztVU,3059
 pygeodesy/nvectorBase.py,sha256=XnA_DlqbgYPpF7y_BQHQVmntpI3M5fe0q8h3zX-PKbs,29930
 pygeodesy/osgr.py,sha256=4mFUcn1vMOju3PK1rKohPtIlJZrBGPJrcRjFl-CoxqE,30852
 pygeodesy/points.py,sha256=2vQpgqxJiC7Cgn03pqiezTM4I-CTquv9tQZi4uYWjpw,64497
-pygeodesy/props.py,sha256=ALaBNSPLpaDjFTtwdYeTVqMZhGqLLiERTCAmTrGrjd0,24608
+pygeodesy/props.py,sha256=GCBqrkMfpbHuca2Sh68WSGDplIeZLvCDDXUewyvEiqI,24600
 pygeodesy/resections.py,sha256=vCbnEyNghGkblrPA58-MspaxkZcInnvRgS-oSBIe1hM,43796
 pygeodesy/simplify.py,sha256=O04NSb0ezXM8z3lCcajXbetZuum6PURAOPYTsJ27gxI,22068
 pygeodesy/solveBase.py,sha256=I2oW7zaPoVUkD7IB9Ro0paU-T1vie1lomg83f71twDY,18468
-pygeodesy/sphericalBase.py,sha256=T3vZJV7FSb7Y_cZJDZeIuBL_qh-YvRNBNOuQABwXPp4,32087
-pygeodesy/sphericalNvector.py,sha256=uS40EpEQfVaz3-fwFBNfHg4gOkk1CWcm1tHrvfpQbmc,54824
-pygeodesy/sphericalTrigonometry.py,sha256=fATw2WUWvY3smtg0vwIczIsl9GamsjwJ3zEwM-subv0,64110
+pygeodesy/sphericalBase.py,sha256=hq6dyeeIHBy_T0Q6PIwSvmo1U70Ir5KBeMTwQplzqHQ,32071
+pygeodesy/sphericalNvector.py,sha256=bu0yB4PEq1cs6BS8QFqtljWpVEw7ESwq2yILOzGDuAc,54948
+pygeodesy/sphericalTrigonometry.py,sha256=8q3sevqdeJ-mOukkXLsXKAg06hDxaBhiPHhrdQn3UHY,64411
 pygeodesy/streprs.py,sha256=-4LErmfPLtGt1tLDBHhJam-zYq8p2Z2C5wDCMELEml4,23599
 pygeodesy/trf.py,sha256=gpsOJVER1KhITOtyiENQ4EnRONinLkzbujsfzzxHVv8,119209
 pygeodesy/triaxials.py,sha256=V0MSYJxRgYDhZcF7HXnu7oFzRyWJkhwc6scuNjc6VVg,64046
@@ -72,7 +72,7 @@ pygeodesy/utily.py,sha256=ZJKuwoo97K7GE56-cA8bGrvyz9mN9OoOMnSOEOw0mj0,40025
 pygeodesy/utm.py,sha256=kcXliftP9q0nTxviIOixVyvResK2Jydj-G52CXJiSuQ,31123
 pygeodesy/utmups.py,sha256=31RD8l7M-vcRXU9_wnkb_PNUcKjOr1HAmwjOndjkgXE,13176
 pygeodesy/utmupsBase.py,sha256=PptZ4Mm8BaDUCA8v1sK54o8idw_24kq6loAuYv0bejE,22730
-pygeodesy/vector2d.py,sha256=djVevQCHZWfzbwmPIZRpPFeKY6pkrOAsydrYAW4pqK8,40309
+pygeodesy/vector2d.py,sha256=7R5ciPxXEX98LmMWNqwbds8Y6vpsHKEeCzcBVDRiq6g,40305
 pygeodesy/vector3d.py,sha256=sC7u9tgFFWoWHFxqtO-dvKVuiEzraw50TDU_QpdDbYw,43045
 pygeodesy/vector3dBase.py,sha256=BHenp3NGpfhVbJ6JMgdNLx6vAcwDWDJIxl9hAdp-AVw,36301
 pygeodesy/webmercator.py,sha256=z3Ft5TeHc9FoLA8zQhzTrSuShHSz0aHjanbvj53ohmo,15016
@@ -104,16 +104,16 @@ pygeodesy/geodesicx/_C4_27.py,sha256=grqtju94Nu9V_eQPqe_7aRbLCYflxRuK127pPOPwuDA
 pygeodesy/geodesicx/_C4_30.py,sha256=I4ZuTy64yOq1IVgRtc_kqOL0HnKzXeTea7kTFSolcJs,201921
 pygeodesy/geodesicx/__init__.py,sha256=KXnzEGNtOtOke01y72dkIgT6sAwOP3zd6wgXho1F3T4,2468
 pygeodesy/geodesicx/__main__.py,sha256=ugi6c-4y91WOMJIL0N-wgXhoGhV90QJIFji2MaIGLiA,4329
-pygeodesy/geodesicx/gx.py,sha256=nvErvXV6ELc5_B2qBffFMf_n6qbk_g4ewChEuzbgVKI,60791
-pygeodesy/geodesicx/gxarea.py,sha256=zedg2zYlQv4DK0gwr3JeQWdtWcFuLSjzls8ddrk4Av0,19461
-pygeodesy/geodesicx/gxbases.py,sha256=EF9IL_oQBSlF70_ROag-DRBoU-ASWKlX-l-wfrUEQYs,6130
-pygeodesy/geodesicx/gxline.py,sha256=WaptQVCDqdVoyhLOP2e0g5fq-UuDOjz86B1bEjLs1SY,27653
+pygeodesy/geodesicx/gx.py,sha256=MdnQat0Q4b6h9jlsOCsJLpmpiDSU84vgip7Gny9nBlw,60971
+pygeodesy/geodesicx/gxarea.py,sha256=1z-CQzi3vEls9Rs-okn003cdR9o_fWtO3lORbk_jpu8,19506
+pygeodesy/geodesicx/gxbases.py,sha256=R9IQXwCK6KOPPcxI6qvdDbVcyuqfYJiZy5jcmrfkJ5A,6390
+pygeodesy/geodesicx/gxline.py,sha256=JbFzLraITpW4xoOx8qCEzkF_5z8AgDSywZj0tcSexsg,27873
 pygeodesy/rhumb/__init__.py,sha256=DsE5a16CUYWE5YTuo1U8mnDCtykuWw7VeGx23Yny_VU,2207
 pygeodesy/rhumb/aux_.py,sha256=NeHkNLS2ULh3sTb3eBTtmRih9SZNyVJKMtJj7imKkuI,16117
 pygeodesy/rhumb/bases.py,sha256=P5AQyL-MTbJaO4oKH_Y7qcusQxaHoAii0hJVQ9mT7hU,54250
 pygeodesy/rhumb/ekx.py,sha256=GMY6ik4fgiRrjek-42cpZxU3U2bdoAcY6ZcbZRckpII,24026
 pygeodesy/rhumb/solve.py,sha256=z8z_XYObTgz7w1skNLNcLBpe-EO_r0H4sVcZGlBcWnc,24005
-pygeodesy-25.5.25.dist-info/METADATA,sha256=mvciPhcidpH1u2mDSmMdSRQykcl7CDCqanD46bFFzM0,20102
-pygeodesy-25.5.25.dist-info/WHEEL,sha256=Kh9pAotZVRFj97E15yTA4iADqXdQfIVTHcNaZTjxeGM,110
-pygeodesy-25.5.25.dist-info/top_level.txt,sha256=cEQPatCXzKZqrivpULC5V5fuy9_V_bAwaP_gUGid7pQ,10
-pygeodesy-25.5.25.dist-info/RECORD,,
+pygeodesy-25.5.28.dist-info/METADATA,sha256=fowmWp0l0feRG1MkGeJGcJxfcyDdvoKZu6TrePRyaNY,20102
+pygeodesy-25.5.28.dist-info/WHEEL,sha256=Kh9pAotZVRFj97E15yTA4iADqXdQfIVTHcNaZTjxeGM,110
+pygeodesy-25.5.28.dist-info/top_level.txt,sha256=cEQPatCXzKZqrivpULC5V5fuy9_V_bAwaP_gUGid7pQ,10
+pygeodesy-25.5.28.dist-info/RECORD,,