pygeodesy 26.1.16__py2.py3-none-any.whl → 26.2.2__py2.py3-none-any.whl

pygeodesy/__init__.py

@@ -626,7 +626,7 @@ else:
 
 from pygeodesy.internals import _version2,  _DOT_  # noqa: E402
 # from pygeodesy.interns import _DOT_  # from .internals
-__version__ = '26.01.16'
+__version__ = '26.02.02'
 # see setup.py for similar logic
 version     = _DOT_(*_version2(__version__, n=3))
 

pygeodesy/auxilats/auxDLat.py

@@ -26,7 +26,7 @@ from pygeodesy.utily import atan1, atan2  # from .auxilats.auxily
 from math import cos, sin, sqrt
 
 __all__ = ()
-__version__ = '25.05.12'
+__version__ = '26.01.20'
 
 
 class AuxDLat(AuxLat):
@@ -107,10 +107,9 @@ class AuxDLat(AuxLat):
             t2  = _1_0 + t**2
             Dt *= _2_0 / t2
             sk2 = (d * Dt)**2 * k2
-            d2  = _1_0 - sk2
             c2  = ((_1_0 - t) * (_1_0 + t) / t2)**2 if t else _1_0
             # E(z)/sin(z)
-            Dt *= _Ef._RFRD(c2, d2, _1_0, sk2) - k2 * sx * sy
+            Dt *= _Ef._fRF3RD(c2, _1_0, sk2) - k2 * sx * sy
         return Dt
 
     def DIsometric(self, Phi1, Phi2):

pygeodesy/auxilats/auxLat.py

@@ -48,7 +48,7 @@ except ImportError:  # Python 3.11-
         return pow(_2_0, x)
 
 __all__ = ()
-__version__ = '25.05.12'
+__version__ = '26.01.20'
 
 _TRIPS = 1024  # XXX 2 or 3?
 
@@ -621,19 +621,18 @@ class AuxLat(AuxAngle):
             ka,  kb  = kb,  ka
             ka1, kb1 = kb1, ka1
             sb,  cb  = cb,  sb
-        # now a, b = larger/smaller semiaxis
-        # Beta measured from larger semiaxis
+        # now a, b = larger/smaller semi-axis
+        # Beta measured from larger semi-axis
         # kb, ka = modulus-squared for distance from Beta = 0, pi/2
         # NB kb <= 0; 0 <= ka <= 1
         # sa = b*E(Beta, sqrt(kb))
         # sb = a*E(Beta',sqrt(ka))
         # 1 - ka * (1 - sb2) = 1 - ka + ka*sb2
-        sb2 =  sb**2
-        cb2 =  cb**2
-        da2 =  ka1 + ka * sb2
-        db2 = _1_0 - kb * sb2
+        sb2 = sb**2
+        cb2 = cb**2
+        da2 = ka1 + ka * sb2
         # DLMF Eq. 19.25.9
-        my = b * sb * _Ef._RFRD(cb2, db2, _1_0, kb * sb2)
+        my = b * sb * _Ef._fRF3RD(cb2, _1_0, kb * sb2)
         # DLMF Eq. 19.25.10 with complementary angles
         mx = a * cb * (_Ef.fRF(sb2,  da2, _1_0) * ka1 +
             ka * cb2 * _Ef.fRD(sb2, _1_0,  da2, _3_0) * ka1 +

pygeodesy/elliptic.py

@@ -102,7 +102,7 @@ from math import asin, asinh, atan, ceil, cosh, fabs, floor, radians, \
 # import operator as _operator  # from .fmath
 
 __all__ = _ALL_LAZY.elliptic
-__version__ = '26.01.16'
+__version__ = '26.01.20'
 
 _TolRD  =  zqrt(EPS * 0.002)
 _TolRF  =  zqrt(EPS * 0.030)
@@ -1212,6 +1212,17 @@ class Elliptic(_Named):
         except Exception as X:
             raise _ellipticError(Elliptic.fRF, x, y, z, cause=X)
 
+    @staticmethod
+    def _fRF3RD(x, z, m):  # in .auxilats.AuxDLat.DE, -.AuxLat.Rectifying
+        y = _1_0 - m
+        try:  # float(RF(x, y, z) - RD(x, y, z, 3 / m))
+            R = _RF3(x, y, z)
+            if m:
+                R -= _RD(x, y, z, _3_0 / m)
+        except Exception as X:
+            raise _ellipticError(Elliptic._fRF3RD, x, y, z, m, cause=X)
+        return float(R)
+
     @staticmethod
     def fRG(x, y, z=0):
         '''Symmetric or complete symmetric integral of the second kind
@@ -1246,16 +1257,6 @@ class Elliptic(_Named):
         except Exception as X:
             raise _ellipticError(Elliptic.fRJ, x, y, z, p, cause=X)
 
-    @staticmethod
-    def _RFRD(x, y, z, m):  # in .auxilats.AuxDLat.DE, -.AuxLat.Rectifying
-        try:  # float(RF(x, y, z) - RD(x, y, z, 3 / m))
-            R = _RF3(x, y, z)
-            if m:
-                R -= _RD(x, y, z, _3_0 / m)
-        except Exception as X:
-            raise _ellipticError(Elliptic._RFRD, x, y, z, m, cause=X)
-        return float(R)
-
 _allPropertiesOf_n(16, Elliptic)  # PYCHOK assert, see Elliptic.reset
 
 

pygeodesy/fsums.py

@@ -52,9 +52,8 @@ from pygeodesy.interns import NN, _arg_, _COMMASPACE_, _DMAIN_, _DOT_, _from_, \
 # from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS  # from .named
 from pygeodesy.named import _name__, _name2__, _Named, _NamedTuple, \
                             _NotImplemented,  _ALL_LAZY, _MODS
-from pygeodesy.props import _allPropertiesOf_n, deprecated_method, \
-                             deprecated_property_RO, Property, \
-                             Property_RO, property_RO
+from pygeodesy.props import _allPropertiesOf_n, deprecated_method, Property, \
+                             deprecated_property_RO, Property_RO, property_RO
 from pygeodesy.streprs import Fmt, fstr, unstr
 # from pygeodesy.units import Float, Int  # from .constants
 
@@ -62,7 +61,7 @@ from math import fabs, isinf, isnan, \
                  ceil as _ceil, floor as _floor  # PYCHOK used! .ltp
 
 __all__ = _ALL_LAZY.fsums
-__version__ = '26.01.16'
+__version__ = '26.02.02'
 
 from pygeodesy.interns import (
   _PLUS_     as _add_op_,  # in .auxilats.auxAngle
@@ -197,8 +196,7 @@ except  ImportError:  # PYCHOK DSPACE! Python 3.12-
     _2FACTOR = pow(2, (MANT_DIG + 1) // 2) + _1_0  # 134217729 if MANT_DIG == 53
 
     def _2split3(x):
-        # Split U{Algorithm 3.2
-        # <https://www.TUHH.De/ti3/paper/rump/OgRuOi05.pdf>}
+        # Split U{Algorithm 3.2<https://www.TUHH.De/ti3/paper/rump/OgRuOi05.pdf>}
         a  = c = x * _2FACTOR
         a -= c - x
         b  = x - a
@@ -326,13 +324,31 @@ def nonfiniterrors(raiser=None):
          _xkwds_get1(d, _isfine=_isfinite) is _isfinite) if d else True
 
 
-def _1primed(xs):  # in .fmath
-    '''(INTERNAL) 1-Primed summation of iterable C{xs}
+# def _nsum(xs):
+#     '''(INTERNAL) U{Neumaier summation
+#        <https://StackOverflow.com/questions/78633770/can-neumaier-summation-be-sped-up>},
+#        see IV. Verbessertes Kahan-Babuška-Verfahren.
+#     '''
+#     s = r = _0_0
+#     for x in map(float, xs):
+#         t = s + x
+#         if fabs(x) <= fabs(s):
+#             r += (s - t) + x
+#         else:
+#             r += (x - t) + s
+#         s = t
+#     return s + r
+
+
+def _1primed(xs, *ys):  # in .fmath
+    '''(INTERNAL) 1-Primed summation of iterable C{xs} less any C{ys}
        items, all I{known} to be C{scalar}.
     '''
     yield _1_0
     for x in xs:
         yield x
+    for y in ys:
+        yield -y
     yield _N_1_0
 
 
@@ -485,10 +501,10 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
               i.e. any C{type} having method C{__float__}.
 
        @note: Handling of I{non-finites} as C{inf}, C{INF}, C{NINF}, C{nan} and C{NAN} is
-              determined by function L{nonfiniterrors<fsums.nonfiniterrors>} for the default
-              or by method L{nonfinites<Fsum.nonfinites>} for individual C{Fsum} instances,
-              overruling the default.  For backward compatibility, I{non-finites} raise
-              exceptions by default.
+              determined globally by function L{nonfiniterrors<fsums.nonfiniterrors>} or
+              by method L{nonfinites<Fsum.nonfinites>} for individual C{Fsum} instances,
+              overruling the global setting.  For backward compatibility, I{non-finites}
+              raise exceptions by default.
 
        @see: U{Hettinger<https://GitHub.com/ActiveState/code/tree/master/recipes/Python/
              393090_Binary_floating_point_summatiaccurate_full/recipe-393090.py>},
@@ -1333,11 +1349,10 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
                     H._fadd(c, up=False)
             else:  # x == 0
                 H = cs[0] if n else 0
-            self._fadd(H)
+            return self._fadd(H)
         except Exception as X:
             t = unstr(where, x, *cs, _ELLIPSIS=4, incx=incx)
             raise self._ErrorX(X, _add_op_, t)
-        return self
 
     def _finite(self, other, op=None):
         '''(INTERNAL) Return B{C{other}} if C{finite}.
@@ -2286,7 +2301,7 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
         # assert isscalar(s) and isscalar(x)
         return self._pow_2_3(s, x, other, op, **raiser_RESIDUAL)
 
-    def _ps_acc(self, ps, xs, up=True, **unused):
+    def _ps_acc(self, ps, xs, up=True, **unused):  # in .geoids._Dotf and ._Hornerf
         '''(INTERNAL) Accumulate C{xs} known scalars into list C{ps}.
         '''
         n   =  0
@@ -2356,15 +2371,7 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
     def _ps_1sum(self, *less):
         '''(INTERNAL) Return the partials sum, 1-primed C{less} some scalars.
         '''
-        def _1psls(ps, ls):
-            yield _1_0
-            for p in ps:
-                yield  p
-            for p in ls:
-                yield -p
-            yield _N_1_0
-
-        return _fsum(_1psls(self._ps, less))
+        return _fsum(_1primed(self._ps, *less))
 
     def _raiser(self, r, s, raiser=True, **RESIDUAL):
         '''(INTERNAL) Does ratio C{r / s} exceed the RESIDUAL threshold
@@ -2726,7 +2733,7 @@ try:
 except ImportError:
     _sum =  sum
 
-    def _fsum(xs):  # in .elliptic
+    def _fsum(xs):  # in .elliptic, .geoids
         '''(INTERNAL) Precision summation, Python 2.5-.
         '''
         F = Fsum(name=_fsum.name, f2product=False, nonfinites=True)
@@ -2867,8 +2874,11 @@ if __name__ == _DMAIN_:
         # copied from Hettinger, see L{Fsum} reference
         from pygeodesy import frandoms, printf
 
+#       printf(typename(_sum),  end=_COMMASPACE_)
         printf(typename(_fsum), end=_COMMASPACE_)
         printf(typename(_psum), end=_COMMASPACE_)
+        printf(len(Fsum.__dict__), end=_COMMASPACE_)
+#       printf(len(globals()), end=_COMMASPACE_)
 
         F = Fsum()
         if F.is_math_fsum():

pygeodesy/geoids.py

@@ -89,18 +89,17 @@ courtesy of SBFRF.
 # make sure int/int division yields float quotient, see .basics
 from __future__ import division as _; del _  # noqa: E702 ;
 
-from pygeodesy.basics import _isin, len2, min2, isodd, _splituple, \
-                              ub2str as _ub2str
+from pygeodesy.basics import _isin, len2, min2, isodd, _splituple
 from pygeodesy.constants import EPS, _float as _F, _1_0, _N_90_0, _180_0, \
                                      _N_180_0, _360_0
 from pygeodesy.datums import Datums, _ellipsoidal_datum, _WGS84
 # from pygeodesy.dms import parseDMS2  # _MODS
 from pygeodesy.errors import _incompatible, LenError, RangeError, _SciPyIssue, \
                              _xkwds_pop2
-from pygeodesy.fmath import favg, Fdot, fdot, Fhorner, frange
+from pygeodesy.fmath import favg, frange,  Fsum
 # from pygoedesy.formy import heightOrthometric  # _MODS
-from pygeodesy.heights import _as_llis2, _ascalar, _HeightBase, HeightError, \
-                              _Wrap
+# from pygeodesy.fsums import Fsum  # from .fmath
+from pygeodesy.heights import _as_llis2, _ascalar, _HeightBase, HeightError,  _Wrap
 # from pygeodesy.internals import typename, _version2  # _MODS
 from pygeodesy.interns import NN, _COLONSPACE_, _COMMASPACE_, _DMAIN_, _E_, \
                              _height_, _in_, _kind_, _lat_, _lon_, _mean_, _N_, \
@@ -115,8 +114,8 @@ from pygeodesy.units import Height, Int_, Lat, Lon
 # from pygeodesy.utily import _Wrap  # from .heights
 
 from math import floor as _floor
-# from os import SEEK_CUR, SEEK_SET  # _MODS
-# import os.path  # _MODS
+# import os as _os  # _MODS
+# import os.path  # _os.path
 from struct import calcsize as _calcsize, unpack as _unpack
 try:
     from StringIO import StringIO as _BytesIO  # reads bytes
@@ -124,9 +123,10 @@ try:
 
 except ImportError:  # Python 3+
     from io import BytesIO as _BytesIO  # PYCHOK expected
+    from pygeodesy.basics import ub2str as _ub2str
 
 __all__ = _ALL_LAZY.geoids
-__version__ = '25.09.26'
+__version__ = '26.02.02'
 
 _assert_    = 'assert'
 _bHASH_     = b'#'
@@ -135,9 +135,50 @@ _format_    = '%s %r'
 _header_    = 'header'
 _intCs      = {}  # cache int value, del below
 _lli_       = 'lli'
+_os         = _MODS.os
 _rb_        = 'rb'
 _supported_ = 'supported'
 
+if __debug__:
+    from pygeodesy.fmath import Fdot as _Dotf, Fhorner as _Hornerf
+
+else:  # -OO ... runs GeoidKarney 8+X faster (w/o Fwelford)
+    from pygeodesy.fsums import _fsum
+    from operator import mul as _mul
+
+    class _GKsum(Fsum):  # for GeoidKarney only
+
+        def __iadd__(self, other):
+            xs = other._ps if isinstance(other, _GKsum) else (other,)
+            self._ps_acc(self._ps, xs, up=False)
+            return self
+
+        def __imul__(self, x):
+            self._ps[:] = (x * p for p in self._ps)
+            return self
+
+        fmul = __imul__  # like .fsums.Fsum
+
+        def fsum(self):  # PYCHOK signature
+            return _fsum(self._ps)
+
+    class _Dotf(_GKsum):  # PYCHOK redef
+        '''Fast precision dot product M{sum(a[i] * b[i] for i=0..len(a))}
+           but for C{float} C{a} and C{b} only.
+        '''
+        def __init__(self, a, *b):
+            self._ps = self._ps_acc([], map(_mul, a, b), up=False)
+
+    class _Hornerf(_GKsum):  # PYCHOK redef
+        '''Fast polynomial evaluation M{sum(cs[i] * x**i for i=0..len(cs))}
+           but for C{float} C{x} and C{cs} and C{incx=True} only.
+        '''
+        def __init__(self, x, *cs):  # incx=True
+            self._ps = []
+            for c in reversed(cs):  # multiply-accumulate
+                self *= x  # ps[:] = (x * p for p in ps)
+                self += c  # self._ps_acc(ps, (c,), up=False)
+
 
 class GeoidError(HeightError):
     '''Geoid interpolator C{Geoid...} or interpolation issue.
@@ -484,7 +525,7 @@ class _GeoidBase(_HeightBase):
 
     def _load(self, g, dtype=float, n=-1, offset=0, **sep):  # sep=NN
         # numpy.fromfile, like .frombuffer
-        g.seek(offset, _MODS.os.SEEK_SET)
+        g.seek(offset, _os.SEEK_SET)
         return self.numpy.fromfile(g, dtype, count=n, **sep)
 
     @Property_RO
@@ -566,8 +607,8 @@ class _GeoidBase(_HeightBase):
     def _open(self, geoid, datum, kind, name, smooth):
         # open the geoid file
         try:
-            self._geoid = _MODS.os.path.basename(geoid)
-            self._sizeB = _MODS.os.path.getsize(geoid)
+            self._geoid = _os.path.basename(geoid)
+            self._sizeB = _os.path.getsize(geoid)
             g = open(geoid, _rb_)
         except (IOError, OSError) as x:
             raise GeoidError(geoid=geoid, cause=x)
@@ -774,15 +815,11 @@ class GeoidEGM96(_GeoidBase):
 
     def _g2ll2(self, lat, lon):
         # convert grid (lat, lon) to earth (lat, lon)
-        while lon > _180_0:  # Eastern
-            lon -= _360_0
-        return -lat, lon  # invert lat
+        return -lat, _lonE2lon(lon)  # invert lat
 
     def _ll2g2(self, lat, lon):
         # convert earth (lat, lon) to grid (lat, lon)
-        while lon < 0:  # Eastern
-            lon += _360_0
-        return -lat, lon  # invert lat
+        return -lat, _lon2lonE(lon)  # invert lat
 
     if _FOR_DOCS:
         __call__ = _GeoidBase.__call__
@@ -791,29 +828,27 @@ class GeoidEGM96(_GeoidBase):
 
 class GeoidG2012B(_GeoidBase):
     '''Geoid height interpolator for U{GEOID12B Model
-       <https://www.NGS.NOAA.gov/GEOID/GEOID12B/>} grids U{CONUS
-       <https://www.NGS.NOAA.gov/GEOID/GEOID12B/GEOID12B_CONUS.shtml>},
-       U{Alaska<https://www.NGS.NOAA.gov/GEOID/GEOID12B/GEOID12B_AK.shtml>},
-       U{Hawaii<https://www.NGS.NOAA.gov/GEOID/GEOID12B/GEOID12B_HI.shtml>},
+       <https://Geodesy.NOAA.gov/GEOID/GEOID12B/>} grids U{CONUS
+       <https://Geodesy.NOAA.gov/GEOID/GEOID12B/GEOID12B_CONUS.shtml>},
+       U{Alaska<https://Geodesy.NOAA.gov/GEOID/GEOID12B/GEOID12B_AK.shtml>},
+       U{Hawaii<https://Geodesy.NOAA.gov/GEOID/GEOID12B/GEOID12B_HI.shtml>},
        U{Guam and Northern Mariana Islands
-       <https://www.NGS.NOAA.gov/GEOID/GEOID12B/GEOID12B_GMNI.shtml>},
+       <https://Geodesy.NOAA.gov/GEOID/GEOID12B/GEOID12B_GMNI.shtml>},
        U{Puerto Rico and U.S. Virgin Islands
-       <https://www.NGS.NOAA.gov/GEOID/GEOID12B/GEOID12B_PRVI.shtml>} and
-       U{American Samoa<https://www.NGS.NOAA.gov/GEOID/GEOID12B/GEOID12B_AS.shtml>}
+       <https://Geodesy.NOAA.gov/GEOID/GEOID12B/GEOID12B_PRVI.shtml>} and
+       U{American Samoa<https://Geodesy.NOAA.gov/GEOID/GEOID12B/GEOID12B_AS.shtml>}
        based on C{SciPy} interpolation U{RectBivariateSpline<https://docs.SciPy.org/doc/
        scipy/reference/generated/scipy.interpolate.RectBivariateSpline.html>}, U{interp2d
        <https://docs.SciPy.org/doc/scipy/reference/generated/scipy.interpolate.interp2d.html>}
        or U{bisplrep/-ev<https://docs.scipy.org/doc/scipy/reference/generated/scipy.interpolate.
        bisplrep.html>}.
-
-       Use any of the C{le} (little endian) or C{be} (big endian) C{g2012b*.bin} binary grid files.
     '''
     _datum = Datums.NAD83
 
     def __init__(self, g2012b_bin, datum=Datums.NAD83, kind=3, smooth=0, **name_crop):
         '''New L{GeoidG2012B} interpolator.
 
-           @arg g2012b_bin: A C{GEOID12B} grid file name (C{.bin}).
+           @arg g2012b_bin: A C{GEOID12B} grid file name (C{.bin}, see B{note}).
            @kwarg datum: Optional grid datum (L{Datum}, L{Ellipsoid}, L{Ellipsoid2} or
                          L{a_f2Tuple}), overriding C{NAD83}.
            @kwarg kind: C{scipy.interpolate} order (C{int}), use 1..5 for U{RectBivariateSpline
@@ -827,8 +862,8 @@ class GeoidG2012B(_GeoidBase):
            @kwarg name_crop: Optional geoid C{B{name}=NN} (C{str}) and UNSUPPORTED keyword argument
                        C{B{crop}=None}.
 
-           @raise GeoidError: Invalid B{C{crop}}, B{C{kind}} or B{C{smooth}} or a G2012B grid file
-                              B{C{g2012b_bin}} issue.
+           @raise GeoidError: Invalid B{C{crop}}, B{C{kind}} or B{C{smooth}} or a B{C{g2012b_bin}}
+                              grid file issue.
 
            @raise ImportError: Package C{numpy} or C{scipy} not found or not installed.
 
@@ -840,8 +875,11 @@ class GeoidG2012B(_GeoidBase):
 
            @raise TypeError: Invalid B{C{datum}}.
 
-           @note: Specify C{B{kind}=-1, -3 or -5} to use C{scipy.interpolate.interp2d}
-                  before or C{scipy.interpolate.bisplrep/-ev} since C{Scipy} version 1.14.
+           @note: Only use any of the C{le} (little-endian) or C{be} (big-endian) C{g2012b*.bin}
+                  I{binary} grid files.
+
+           @note: Specify C{B{kind}=-1, -3 or -5} to use C{scipy.interpolate.interp2d} I{before}
+                  or C{scipy.interpolate.bisplrep/-ev} I{since} C{Scipy} version 1.14.
         '''
         crop, name = _xkwds_pop2(name_crop, crop=None)
         if crop is not None:
@@ -871,7 +909,7 @@ class GeoidG2012B(_GeoidBase):
             p.slat, p.wlon, p.dlat, p.dlon = map(float, self._load(g, en_+'f8', 4))
             # read all f4 heights, ignoring the first 4xf8 and 3xi4
             hs = self._load(g, self._endian, n, 44).reshape(p.nlat, p.nlon)
-            p.wlon -= _360_0  # western-most East longitude to earth (..., lon)
+            p.wlon -= _360_0  # western-most lon XXX _lonE2lon
             _GeoidBase.__init__(self, hs, p)
 
         except Exception as x:
@@ -975,12 +1013,12 @@ class GeoidKarney(_GeoidBase):
     _CM = (_T(' 0, -1'),  # 10x12 cubic matrix [i, j] indices
            _T(' 1, -1'),
            _T('-1,  0'),
-           _T(' 0,  0'),
-           _T(' 1,  0'),
+           _T(' 0,  0'),  # _BT[0]
+           _T(' 1,  0'),  # _BT[1]
            _T(' 2,  0'),
            _T('-1,  1'),
-           _T(' 0,  1'),
-           _T(' 1,  1'),
+           _T(' 0,  1'),  # _BT[2]
+           _T(' 1,  1'),  # _BT[3]
            _T(' 2,  1'),
            _T(' 0,  2'),
            _T(' 1,  2'))
@@ -1063,7 +1101,7 @@ class GeoidKarney(_GeoidBase):
         p = self._pgm
         if 0 < x < (p.nlon - 2) and 0 < y < (p.nlat - 2):
             # read 4x4 ushorts, drop the 4 corners
-            S = _MODS.os.SEEK_SET
+            S = _os.SEEK_SET
             e =  self._4endian
             g =  self._egm
             n =  self._4u2B
@@ -1099,7 +1137,7 @@ class GeoidKarney(_GeoidBase):
         y, x = int(_floor(fy)), int(_floor(fx))
         fy -= y
         fx -= x
-        H  = self._ev2d(fy, fx, y, x)  # PYCHOK ._ev2k or ._ev3k
+        H  = self._ev2d(fy, fx, y, x)  # PYCHOK ._ev3k or ._ev2k
         H *= self._pgm.Scale   # H.fmul(self._pgm.Scale)
         H += self._pgm.Offset  # H.fadd(self._pgm.Offset)
         return H.fsum()  # float(H)
@@ -1112,10 +1150,10 @@ class GeoidKarney(_GeoidBase):
             y, x = self._yx_i = yx
             self._yx_t = self._raws(y, x, GeoidKarney._BT)
         t  =  self._yx_t
-        v  = _1_0, -fx, fx
-        H  = Fdot(v, t[0], t[0], t[1]).fmul(_1_0 - fy)  # c = a * (1 - fy)
-        H += Fdot(v, t[2], t[2], t[3]).fmul(fy)  # c += b * fy
-        return H
+        v  = _1_0, (-fx), fx
+        H  = _Dotf(v, t[0], t[0], t[1]).fmul(_1_0 - fy)  # c = a * (1 - fy)
+        H += _Dotf(v, t[2], t[2], t[3]).fmul(fy)  # c += b * fy
+        return H  # Fsum
 
     def _ev3k(self, fy, fx, *yx):
         # compute the cubic 10-tuple and interpolate raw H
@@ -1124,7 +1162,7 @@ class GeoidKarney(_GeoidBase):
         else:
             c0, c3, v  = self._c0c3v(*yx)
             # assert len(c3) == self._nterms
-            self._yx_t = tuple(fdot(v, *r3) / c0 for r3 in c3)
+            self._yx_t = tuple(_Dotf(v, *r3).fover(c0) for r3 in c3)
             self._yx_i = yx
         # GeographicLib/Geoid.cpp Geoid::height(lat, lon) ...
         # real h = t[0] + fx * (t[1] + fx * (t[3] + fx * t[6])) +
@@ -1132,20 +1170,18 @@ class GeoidKarney(_GeoidBase):
         #                 fy * (t[5] + fx *  t[8] + fy * t[9]));
         t  =  self._yx_t
         v  = _1_0, fx, fy
-        H  = Fdot(v, t[5], t[8], t[9])
-        H *= fy
-        H += Fhorner(fx, t[2], t[4], t[7])
-        H *= fy
-        H += Fhorner(fx, t[0], t[1], t[3], t[6])
-        return H
+        H  = _Dotf(v, t[5], t[8], t[9])
+        H *=  fy
+        H += _Hornerf(fx, t[2], t[4], t[7])
+        H *=  fy
+        H += _Hornerf(fx, t[0], t[1], t[3], t[6])
+        return H  # Fsum
 
     _ev2d = _ev3k  # overriden for kind=2, see ._ev_name
 
     def _g2ll2(self, lat, lon):
         # convert grid (lat, lon) to earth (lat, lon), uncropped
-        while lon > _180_0:
-            lon -= _360_0
-        return lat, lon
+        return lat, _lonE2lon(lon)
 
     def _g2yx2(self, lat, lon):
         # convert grid (lat, lon) to grid (y, x) indices
@@ -1188,9 +1224,7 @@ class GeoidKarney(_GeoidBase):
 
     def _ll2g2(self, lat, lon):
         # convert earth (lat, lon) to grid (lat, lon), uncropped
-        while lon < 0:
-            lon += _360_0
-        return lat, lon
+        return lat, _lon2lonE(lon)
 
     def _llh3minmax(self, highest, *lat2):
         # find highest or lowest, takes 10+ secs for egm2008-1.pgm geoid
@@ -1277,7 +1311,7 @@ class GeoidKarney(_GeoidBase):
         p, g = self._pgm, self._egm
         if g:
             b = p.skip + (y * p.nlon + x) * self._u2B
-            g.seek(b, _MODS.os.SEEK_SET)
+            g.seek(b, _os.SEEK_SET)
             return b  # position
         raise GeoidError('closed file', txt=repr(p.egm))  # IOError
 
@@ -1347,7 +1381,7 @@ class GeoidPGM(_GeoidBase):
                   the available memory, especially with 32-bit Python.  To reduce memory
                   usage, set keyword argument B{C{crop}} to the region of interest.  For
                   example C{B{crop}=(20, -125, 50, -65)} covers the U{conterminous US
-                  <https://www.NGS.NOAA.gov/GEOID/GEOID12B/maps/GEOID12B_CONUS_grids.png>}
+                  <https://Geodesy.NOAA.gov/GEOID/GEOID12B/maps/GEOID12B_CONUS_grids.png>}
                   (CONUS), less than 3% of the entire C{egm2008-1.pgm} dataset.
 
            @see: Class L{GeoidKarney} and function L{egmGeoidHeights}.
@@ -1382,8 +1416,7 @@ class GeoidPGM(_GeoidBase):
         if self._cropped:
             lon -= self._lon_of
         else:
-            while lon > _180_0:
-                lon -= _360_0
+            lon = _lonE2lon(lon)
         return lat, lon
 
     def _ll2g2(self, lat, lon):
@@ -1391,8 +1424,7 @@ class GeoidPGM(_GeoidBase):
         if self._cropped:
             lon += self._lon_of
         else:
-            while lon < 0:
-                lon += _360_0
+            lon = _lon2lonE(lon)
         return lat, lon
 
     if _FOR_DOCS:
@@ -1541,8 +1573,8 @@ class _PGM(_Gpars):
         self.slat, self.wlon = self.Origin
         # note, negative .dlat and .rlat since rows
         # are from .slat 90N down in decreasing lat
-        self.dlat, self.dlon = _180_0 / (1 - nlat), _360_0 / nlon
-        self.rlat, self.rlon = (1 - nlat) / _180_0, nlon / _360_0
+        self.dlat, self.dlon = (_180_0 / (1 - nlat)), (_360_0 / nlon)
+        self.rlat, self.rlon = ((1 - nlat) / _180_0), (nlon / _360_0)
 
         # grid corners in earth (lat, lon), .slat = 90, .dlat < 0
         n = float(self.slat)
@@ -1604,11 +1636,11 @@ class _PGM(_Gpars):
         t, c = 0, self._tmpfile()
         # reading (s - n) rows, forward
         for y in range(n, s):  # PYCHOK y unused
-            g.seek(z, _MODS.os.SEEK_SET)
+            g.seek(z, _os.SEEK_SET)
             # Python 2 tmpfile.write returns None
             t += c.write(g.read(r)) or r
             if p:  # wrap around to start of row
-                g.seek(-q, _MODS.os.SEEK_CUR)
+                g.seek(-q, _os.SEEK_CUR)
                 # assert(g.tell() == (z - w * self.u2B))
                 # Python 2 tmpfile.write returns None
                 t += c.write(g.read(p)) or p
@@ -1635,7 +1667,7 @@ class _PGM(_Gpars):
         self.knots = k
         self.skip  = 0  # no header lines in c
 
-        c.seek(0, _MODS.os.SEEK_SET)
+        c.seek(0, _os.SEEK_SET)
         # c = open(c.name, _rb_)  # reopen for numpy 1.8.0-
         return c
 
@@ -1659,11 +1691,11 @@ class _PGM(_Gpars):
         try:
             from tempfile import NamedTemporaryFile as tmpfile
         except ImportError:  # Python 2.7.16-
-            from _MODS.os import tmpfile
-        t = _MODS.os.path.basename(self.pgm)
-        t = _MODS.os.path.splitext(t)[0]
+            from _os import tmpfile  # PYCHOK from
+        t = _os.path.basename(self.pgm)
+        t = _os.path.splitext(t)[0]
         f =  tmpfile(mode='w+b', prefix=t or 'egm')
-        f.seek(0, _MODS.os.SEEK_SET)  # force overwrite
+        f.seek(0, _os.SEEK_SET)  # force overwrite
         return f
 
     @Property_RO
@@ -1703,7 +1735,7 @@ def egmGeoidHeights(GeoidHeights_dat):
         dat = _BytesIO(dat)
 
     try:
-        dat.seek(0, _MODS.os.SEEK_SET)  # reset
+        dat.seek(0, _os.SEEK_SET)  # reset
     except AttributeError as x:
         raise GeoidError(GeoidHeights_dat=type(dat), cause=x)
 
@@ -1711,11 +1743,26 @@ def egmGeoidHeights(GeoidHeights_dat):
         t = t.strip()
         if t and not t.startswith(_bHASH_):
             lat, lon, egm84, egm96, egm2008 = map(float, t.split())
-            while lon > _180_0:  # EasternLon to earth lon
-                lon -= _360_0
+            lon = _lonE2lon(lon)  # Eastern to earth lon
             yield GeoidHeight5Tuple(lat, lon, egm84, egm96, egm2008)
 
 
+def _lonE2lon(lon):
+    '''(INTERNAL) East to earth longitude.
+    '''
+    while lon > _180_0:
+        lon -= _360_0
+    return lon
+
+
+def _lon2lonE(lon):
+    '''(INTERNAL) Earth to East longitude.
+    '''
+    while lon < 0:
+        lon += _360_0
+    return lon
+
+
 __all__ += _ALL_DOCS(_GeoidBase)
 
 if __name__ == _DMAIN_:  # MCCABE 14

pygeodesy/heights.py

@@ -59,7 +59,7 @@ C{>>> h1, h2, ... = hinterpolator.height_(lat1, lon1, lat2, lon2, ...)}
        L{HeightSmoothBiSpline} require package U{scipy<https://SciPy.org>}.
        Classes L{HeightIDWkarney} and L{HeightIDWdistanceTo} -if used with
        L{ellipsoidalKarney.LatLon} points- require I{Karney}'s U{geographiclib
-       <https://PyPI.org/project/geographiclib>} to be installed.
+       <https://PyPI.org/project/geographiclib>} package to be installed.
 
 @note: Errors from C{scipy} are raised as L{SciPyError}s.  Warnings issued
        by C{scipy} can be thrown as L{SciPyWarning} exceptions, provided
@@ -92,7 +92,7 @@ from pygeodesy.units import _isDegrees, Float_, Int_
 # from math import radians  # from .points
 
 __all__ = _ALL_LAZY.heights
-__version__ = '25.09.29'
+__version__ = '26.02.02'
 
 _error_  = 'error'
 _formy   = _MODS.into(formy=__name__)
@@ -160,7 +160,7 @@ def _orderedup(ts, lo=EPS, hi=PI2-EPS):
     return sorted(set(max(lo, min(hi, t)) for t in ts))  # list
 
 
-def _xyhs(wrap=False, _lat=_90_0, _lon=_180_0, **name_lls):
+def _xyhs(wrap=False, _lat=_90_0, _lon=_180_0, height=True, **name_lls):
     # map (lat, lon, h) to (x, y, h) in radians, offset
     # x as 0 <= lon <= PI2 and y as 0 <= lat <= PI
     name, lls = _xkwds_item2(name_lls)
@@ -168,8 +168,9 @@ def _xyhs(wrap=False, _lat=_90_0, _lon=_180_0, **name_lls):
     try:
         for i, ll in enumerate(lls):
             y, x = _w(ll.lat, ll.lon)
-            yield max(_0_0, _r(x + _lon)), \
-                  max(_0_0, _r(y + _lat)), ll.height
+            h = ll.height if height else 0
+            yield (max(_0_0, _r(x + _lon)),
+                   max(_0_0, _r(y + _lat)), h)
     except Exception as x:
         i = Fmt.INDEX(name, i)
         raise HeightError(i, ll, cause=x)
@@ -266,9 +267,9 @@ class _HeightBase(_HeightNamed):  # in .geoids
         # convert lli C{LatLon}s to tuples or C{NumPy} arrays of
         # C{SciPy} sphericals and determine the return type
         atype = self.numpy.array
-        wrap = _xkwds(wrap, wrap=self._wrap)
+        kwds = _xkwds(wrap, wrap=self._wrap, height=False)
         _as, llis   = _as_llis2(llis)
-        xis, yis, _ =  zip(*_xyhs(llis=llis, **wrap))  # PYCHOK yield
+        xis, yis, _ =  zip(*_xyhs(llis=llis, **kwds))  # PYCHOK yield
         return _as, atype(xis), atype(yis), llis
 
     def _ev(self, *args):  # PYCHOK no cover

{pygeodesy-26.1.16.dist-info → pygeodesy-26.2.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: pygeodesy
-Version: 26.1.16
+Version: 26.2.2
 Summary: Pure Python geodesy tools
 Home-page: https://GitHub.com/mrJean1/PyGeodesy
 Author: Jean M. Brouwers
@@ -8,7 +8,7 @@ Author-email: mrJean1@Gmail.com
 Maintainer: Jean M. Brouwers
 Maintainer-email: mrJean1@Gmail.com
 License: MIT
-Keywords: AER AGM Albers altitude Andoyer annulus antipode area Arithmetic-Geometric-Mean attitude Authalic auxiliary azimuth azimuthal azimuth-elevation-range bearing bank Barsky Barth beta bi-quadratic boolean cached Cagnoli cartesian Cassini Cassini-Soldner chord circle-intersections circumcenter circumcircle circumradius clip Cohen Cohen-Sutherland Collins composite conformal conformal-sphere conic constants contact-triangle Cook Correia cosines-law coverage curvature cylindrical datum deprecation deficit development discrete distance Douglas earth east-north-up eccentricity ECEF elevation ellipsoid ellipsoidal-latitude-beta ellipsoidal-longitude-omega elliptic ENU EPSG equal-area equidistant equirectangular ETM ETRF Euclidean even-odd-rule ExactTM excess Farrell Farrell-Barth Ferrari-solution Field-Of-View flattening fma fmath footpoint footprint Forster Forster-Hormann-Popa Forsythe FOV fractional Frechet Fréchet frustum Fsum fused-multiply-add GARS Gauss-Kummer geocentric GeoConvert GeodesicExact geodesy geodetic GeodSolve Geod3Solve GeodTest geographiclib Geohash geoid geoidHeight GeoidHeights georef Girard gnomonic gons grades gradians Greiner Greiner-Hormann Hartzell Hausdorff Haversine heading hectare height Heikkinen Heron Hodgman horizon Hormann Hubeny IDW incenter incirle infix_@_operator inradius intermediate interpolate intersect intersection intersection3d intersections IntersectTool Inverse-Distance-Weighting Isometric ITRF Jacobi Jacobi-Conformal Jarque-Bera Jekel Karney Krueger Krüger kurtosis Lambert latitude law-of-cosines least-squares Lesh L_Huilier LHuilier Liang Liang-Barsky Linderholm-Segal linearize Line-Of-Sight LocalCartesian local-tangent-plane local-x-y-z longitude LOS loxodrome lstsq LTP lune LV03 LV95 mean memoize memoized Mercator Meeus MGRS nearest NED Niemeyer non-finite normalize Norrdine north-east-down numpy n-vector Nvector oblate omega orthographic orthometric-height OSGB OSGR overlap parallel parallel-of-latitude Parametric path-intersection perimeter Peucker Pierlot pitch plumb Point-Of-View polar Popa POV precision-cubic-root precision-hypotenuse precision-powers precision-running-summation precision-square-root precision-summation prolate Pseudo-Mercator pygeodesy PyInstaller PyPy quartic radical radii radius Ramanujan Ramer Ramer-Douglas-Peucker Rectifying Reduced resect resection Rey-Jer Reumann Reumann-Witkam rho-theta-phi rhumb RhumbSolve running-linear-regression running-statistics running-stats running-summation scipy secant semi-perimeter sexagecimal simplify skewness Snellius Snellius-Pothenot Snyder Soddy Soddy-circles Soldner sphere sphere-intersections spherical-deficit spherical-excess spherical-polar spherical-triangle squared-quartic standard-deviation stereographic Sudano surface-area Sutherland Sutherland-Hodgman tangent-circles Terrestrial-Reference-Frame Thomas Tienstra tilt TMcoords TMExact toise transverse TransverseMercatorExact TRF triangle triangulate triaxial triaxial-ellipsoid trigonometry trilaterate trilaterate-2d trilaterate-3d TwoProduct TwoSum umbilic-point unit unroll UPS UTM UTM/UPS variance velocities Veness Vermeille viewing-frustum Vincenty Visvalingam Visvalingam-Whyatt volume  volumetric Web-Mercator Welford WGRS WGS Whyatt Wildberger Witkam winding-number XYZ yaw You zenzi-cubic zenzi-quartic
+Keywords: AER AGM Albers altitude Andoyer annulus antipode area Arithmetic-Geometric-Mean attitude Authalic auxiliary azimuth azimuthal azimuth-elevation-range bearing bank Barsky Barth beta bi-quadratic boolean cached Cagnoli cartesian Cassini Cassini-Soldner chord circle-intersections circumcenter circumcircle circumradius clip Cohen Cohen-Sutherland Collins composite conformal conformal-sphere conic constants contact-triangle Cook Correia cosines-law coverage curvature cylindrical datum deprecation deficit development discrete distance Douglas earth east-north-up eccentricity ECEF elevation ellipse ellipsoid ellipsoidal-latitude-beta ellipsoidal-longitude-omega elliptic ENU EPSG equal-area equidistant equirectangular ETM ETRF Euclidean even-odd-rule ExactTM excess Farrell Farrell-Barth Ferrari-solution Field-Of-View flattening fma fmath footpoint footprint Forster Forster-Hormann-Popa Forsythe FOV fractional Frechet Fréchet frustum Fsum fused-multiply-add GARS Gauss-Kummer geocentric GeoConvert GeodesicExact geodesy geodetic GeodSolve Geod3Solve GeodTest geographiclib Geohash geoid geoidHeight GeoidHeights georef Girard gnomonic gons grades gradians Greiner Greiner-Hormann Hartzell Hausdorff Haversine heading hectare height Heikkinen Heron Hodgman horizon Hormann Hubeny IDW incenter incirle infix_@_operator inradius intermediate interpolate intersect intersection intersection3d intersections IntersectTool Inverse-Distance-Weighting Isometric ITRF Jacobi Jacobi-Conformal Jarque-Bera Jekel Karney Krueger Krüger kurtosis Lambert latitude law-of-cosines least-squares Lesh L_Huilier LHuilier Liang Liang-Barsky Linderholm-Segal linearize Line-Of-Sight LocalCartesian local-tangent-plane local-x-y-z longitude LOS loxodrome lstsq LTP lune LV03 LV95 mean memoize memoized Mercator Meeus MGRS nearest NED Niemeyer non-finite normalize Norrdine north-east-down numpy n-vector Nvector oblate omega orthographic orthometric-height OSGB OSGR overlap parallel parallel-of-latitude Parametric path-intersection perimeter Peucker Pierlot pitch plumb Point-Of-View polar Popa POV precision-cubic-root precision-hypotenuse precision-powers precision-running-summation precision-square-root precision-summation prolate Pseudo-Mercator pygeodesy PyInstaller PyPy quartic radical radii radius Ramanujan Ramer Ramer-Douglas-Peucker Rectifying Reduced resect resection Rey-Jer Reumann Reumann-Witkam rho-theta-phi rhumb RhumbSolve running-linear-regression running-statistics running-stats running-summation scipy secant semi-perimeter sexagecimal simplify skewness Snellius Snellius-Pothenot Snyder Soddy Soddy-circles Soldner sphere sphere-intersections spherical-deficit spherical-excess spherical-polar spherical-triangle squared-quartic standard-deviation stereographic Sudano surface-area Sutherland Sutherland-Hodgman tangent-circles Terrestrial-Reference-Frame Thomas Tienstra tilt TMcoords TMExact toise transverse TransverseMercatorExact TRF triangle triangulate triaxial triaxial-ellipsoid trigonometry trilaterate trilaterate-2d trilaterate-3d TwoProduct TwoSum umbilic-point unit unroll UPS UTM UTM/UPS variance velocities Veness Vermeille viewing-frustum Vincenty Visvalingam Visvalingam-Whyatt volume  volumetric Web-Mercator Welford WGRS WGS Whyatt Wildberger Witkam winding-number XYZ yaw You zenzi-cubic zenzi-quartic
 Classifier: Development Status :: 5 - Production/Stable
 Classifier: Environment :: Console
 Classifier: Intended Audience :: Developers
@@ -171,7 +171,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: Jan 16, 2026.*
+*Last updated: Feb 02, 2026.*
 
 License
 =======

{pygeodesy-26.1.16.dist-info → pygeodesy-26.2.2.dist-info}/RECORD

@@ -1,5 +1,5 @@
 pygeodesy/LICENSE,sha256=YfgAiyxOwY6P9Kkb1_5XN81nueTLrpb3Ffkv3EuPgFU,1144
-pygeodesy/__init__.py,sha256=DmJmEyYy4jSMO4rkA65hK5ZTWBan1UwlXKGAehPbWoc,43790
+pygeodesy/__init__.py,sha256=22CZUTpgF_0-L5ZLo-wpXgBAtgi8sQFXTN6E8P_Csbw,43790
 pygeodesy/__main__.py,sha256=06sdpwiuTP1hPvKOdzjqhuuDqqpicErwnAYwioeiTNw,5597
 pygeodesy/albers.py,sha256=6S0hHJA6MspeKLsQHW0cFUpJcWRL1yGshJxtEgJ9PMQ,30908
 pygeodesy/angles.py,sha256=68sNVXU3QKGZnIYyjWFcKhWlBKSmaduO--rjtoM2sqU,29996
@@ -23,7 +23,7 @@ pygeodesy/ellipsoidalKarney.py,sha256=CxPB91CKA0q5E97TOcErTTQdufdB2dT4KF84QZtfnG
 pygeodesy/ellipsoidalNvector.py,sha256=48N-KtVUn0wi9kXZe12hg0jABWNeaiNXyY__oS-x5I4,30146
 pygeodesy/ellipsoidalVincenty.py,sha256=X4OpMPtoretveNjHlEXKL4ifmAARTRY0yxbIOlR26J8,26273
 pygeodesy/ellipsoids.py,sha256=5UwnIOpDCFrxqQIF9mkQUEEPW2H1YRog_MpYkvcYx7c,109248
-pygeodesy/elliptic.py,sha256=xCllssJSm2O4W2Mo2hb1f65pT-GcaVwQuIpzR84eD7Y,59483
+pygeodesy/elliptic.py,sha256=nUBHDM9v-dQt9WHHKSBRB74Hd5zysG1TyKom8nlSN-o,59505
 pygeodesy/epsg.py,sha256=wS8Mzzk70Fzry6c0fBA16VVez3ghKG2Dhgpo2eN_Ik8,8220
 pygeodesy/errors.py,sha256=76wVrhxev-n-kBOHv3MicHell3ifA2whUp4SRCrqUx4,32633
 pygeodesy/etm.py,sha256=3pQmpjsrnI0mgLwhkLOCOFq_44bkhHwXExf-Keyky2w,46796
@@ -31,16 +31,16 @@ pygeodesy/fmath.py,sha256=OUitWAuSM7OzlbJJ_JahfUjsVrIv5Bnt4U_0pUgXTQ4,37995
 pygeodesy/formy.py,sha256=TSDEyzOInkhmeZYVRjoOl4oO7wCFrE2YLWGTs9Ewa_4,69842
 pygeodesy/frechet.py,sha256=FySLdZgLb_NXd81AILx4G2C9J6dxwI5uYW1MEHuEocg,35588
 pygeodesy/fstats.py,sha256=EuyC9UrJ5K0IlZqM7h3C8C7tYDrTujqoBp9-V1ICPx0,28348
-pygeodesy/fsums.py,sha256=VmxnlFld-3z-_xuksKsQinWmzZyeo6CEMM9c9ZabMCU,105500
+pygeodesy/fsums.py,sha256=dh9tXnSE1BNH7oEwq4Zt2G6BGNy5Bz9wJLCde6hYvW0,105954
 pygeodesy/gars.py,sha256=z7yvE6dRAjtXUAB9yakH8duwF2gBKbcjXUIznNnr8_w,11761
 pygeodesy/geod3solve.py,sha256=ATQmfkg2KU1YwAgQU32t78Fq5mLvpYykLpqH6TH6mXw,22052
 pygeodesy/geodesici.py,sha256=eHkXKDQzeFfJ6lVkGANqnEk5hQR-MkYXQlgXCLzbG3Y,74933
 pygeodesy/geodesicw.py,sha256=bhwSBtHAaNAl487uAgifS42cKvcYTe-8O8M5S44u9_k,30029
 pygeodesy/geodsolve.py,sha256=GqMn4PjXAT7jK0VjxeMv7FgNx-_T_LeAfeMda7wUSqM,25032
 pygeodesy/geohash.py,sha256=LmjvzgrzD0Cj-xPCZNJp3LbVOQ7llaE2sAO0V0SxNOo,40124
-pygeodesy/geoids.py,sha256=Fn3pOCokzQJuEY2UE11Ah5XuiT24WRmmoc3FeNhfpsg,86129
+pygeodesy/geoids.py,sha256=7r6iEefNKNW6LwacSAnaluiBw3QTl32UbtUEME-kEgU,87606
 pygeodesy/hausdorff.py,sha256=_zfkZ4pzUE_6LwFD77RU9NC58aN8imglI7TfapM82Wc,32269
-pygeodesy/heights.py,sha256=T_5ZhtFAskQf-1BHL2Nm5bTYQIMJj-d1UKhW11YPWvs,41011
+pygeodesy/heights.py,sha256=S7l_Nq4fcfSRJEoFb3_sGK1Gm35AlN6zEjK-Xqafrfk,41082
 pygeodesy/internals.py,sha256=WPrXmN1GMnVB3JYJbXpC1cZlIYMyU9v1I1rDeHI8yGQ,24934
 pygeodesy/interns.py,sha256=h7LupB9eEhwgW8mq3R4VMFI3i-49ObgXAEvV2JJmmpM,23521
 pygeodesy/iters.py,sha256=xH719oFB-0z4jCAHeC6uLqyfWPtpQoSBF6Pgnoz5weQ,20345
@@ -85,9 +85,9 @@ pygeodesy/auxilats/_CX_Rs.py,sha256=EV6Py2yxqtf3FAykmdAHa3Rj9V_yfP21BYwayXNk_zU,
 pygeodesy/auxilats/__init__.py,sha256=u3hwLEvkaiD4KU-kfHfi26IbKWZHWLIBVukCFYNpQjo,2883
 pygeodesy/auxilats/__main__.py,sha256=JPDbuG6UnghFtWBbNzb7Uuv9bXDW55ZXXd2PV-6Tc9o,3340
 pygeodesy/auxilats/auxAngle.py,sha256=QrkiKKHpH7gYoLA5WSdePX0198oDB0yiwGSx0erWLMc,16671
-pygeodesy/auxilats/auxDLat.py,sha256=6bamjB_s0RbINNU3V8W9lzsjXl3tWcjkkwLLgusclVM,10956
+pygeodesy/auxilats/auxDLat.py,sha256=djSAVT5CF2D49F6MhCM9utZcYHoXYu4in0iGXDSAeI8,10925
 pygeodesy/auxilats/auxDST.py,sha256=4VkQ6RfpUA0NvTZawP-bMatJMbrFz3lbzyXKwfQmdB4,10464
-pygeodesy/auxilats/auxLat.py,sha256=qWvIdfpvIQUF7jsgJyxeNgkbjlbIPsNs8-uvPMxyJbM,31972
+pygeodesy/auxilats/auxLat.py,sha256=-q5Ks4e5hzn5ZYp1ThiVTm3Qq-E9uvz-o8izNivhp_c,31938
 pygeodesy/auxilats/auxily.py,sha256=_D-7n-73pVdZR459NF0WGmO3eOVEtiyaongKCGPvHkA,8085
 pygeodesy/deprecated/__init__.py,sha256=TsSpb0A6l5BelXOuhHvljy6m-UoR0Sf86CsCXlgfuls,2823
 pygeodesy/deprecated/bases.py,sha256=5cwF5nZ6SmAwib3iEmOJ66yzLfimfl1Ru93U_ohbtkA,1649
@@ -119,7 +119,7 @@ pygeodesy/triaxials/bases.py,sha256=WfF3Up7JFAWrLnvsRpWsncXICjxqbjVZZUbB6EvGFHU,
 pygeodesy/triaxials/conformal3.py,sha256=FrJLQ0m58YvP_heh9jFH5TkhR46kSTzCIgtmtRPMErQ,23883
 pygeodesy/triaxials/triaxial3.py,sha256=mTG8Atw7VbGyXWPoD78M-6XM6rarXCwKlatqWl41Hoo,42865
 pygeodesy/triaxials/triaxial5.py,sha256=nJhVDbNfXSbldL7QCrGAXNsFlIWY-rUAaVLMPcX7x-g,50977
-pygeodesy-26.1.16.dist-info/METADATA,sha256=U4um_8MNCzFuw0CYVx46wj3uI-aK9MCcIjYbykv5M1s,20716
-pygeodesy-26.1.16.dist-info/WHEEL,sha256=Kh9pAotZVRFj97E15yTA4iADqXdQfIVTHcNaZTjxeGM,110
-pygeodesy-26.1.16.dist-info/top_level.txt,sha256=cEQPatCXzKZqrivpULC5V5fuy9_V_bAwaP_gUGid7pQ,10
-pygeodesy-26.1.16.dist-info/RECORD,,
+pygeodesy-26.2.2.dist-info/METADATA,sha256=hBT8G-DrD-R9WaU18rjEVX7wnMaLl4ZNb78BG6zYBdg,20723
+pygeodesy-26.2.2.dist-info/WHEEL,sha256=Kh9pAotZVRFj97E15yTA4iADqXdQfIVTHcNaZTjxeGM,110
+pygeodesy-26.2.2.dist-info/top_level.txt,sha256=cEQPatCXzKZqrivpULC5V5fuy9_V_bAwaP_gUGid7pQ,10
+pygeodesy-26.2.2.dist-info/RECORD,,