pygeodesy 25.11.5__py2.py3-none-any.whl → 25.12.31__py2.py3-none-any.whl

pygeodesy/ellipsoidalBase.py

@@ -1153,7 +1153,7 @@ __all__ += _ALL_DOCS(CartesianEllipsoidalBase, LatLonEllipsoidalBase)
 
 # **) MIT License
 #
-# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2026 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/ellipsoidalBaseDI.py

@@ -888,7 +888,7 @@ del _1_0, _0_01
 
 # **) MIT License
 #
-# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2026 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/ellipsoidalExact.py

@@ -315,7 +315,7 @@ __all__ += _ALL_OTHER(Cartesian, LatLon,  # classes
 
 # **) MIT License
 #
-# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2026 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/ellipsoidalGeodSolve.py

@@ -314,7 +314,7 @@ __all__ += _ALL_OTHER(Cartesian, LatLon,  # classes
 
 # **) MIT License
 #
-# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2026 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/ellipsoidalKarney.py

@@ -354,7 +354,7 @@ __all__ += _ALL_OTHER(Cartesian, LatLon,  # classes
 
 # **) MIT License
 #
-# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2026 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/ellipsoidalNvector.py

@@ -660,7 +660,7 @@ __all__ += _ALL_OTHER(Cartesian, LatLon, Ned, Nvector,  # classes
 
 # **) MIT License
 #
-# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2026 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/ellipsoidalVincenty.py

@@ -573,7 +573,7 @@ __all__ += _ALL_DOCS(Cartesian, LatLon, intersecant2,  # from .ellipsoidalBaseDI
 
 # **) MIT License
 #
-# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2026 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/ellipsoids.py

@@ -67,10 +67,10 @@ from __future__ import division as _; del _  # noqa: E702 ;
 # from pygeodesy.albers import AlbersEqualAreaCylindrical  # _MODS
 from pygeodesy.basics import copysign0, isbool, _isin, isint,  typename
 from pygeodesy.constants import EPS, EPS_2, EPS0, EPS02, EPS1, INF, NINF, \
-                               _over, PI_2, PI_3, PI4, R_M, R_MA, R_FM, _EPSqrt, \
+                               _EPSqrt, PI_2, PI_3, PI2, PI4, R_M, R_MA, R_FM, \
                                _EPStol as _TOL, _floatuple as _T, _isfinite, \
-                               _0_0s, _0_0, _0_5, _1_0, _1_EPS, _2_0, _4_0, _90_0, \
-                               _0_25, _3_0  # PYCHOK used!
+                               _over, _0_0s, _0_0, _0_5, _1_0, _1_EPS, _2_0, \
+                               _4_0, _90_0, _0_25, _3_0  # PYCHOK used!
 from pygeodesy.errors import _AssertionError, IntersectionError, _ValueError, _xattr, _xkwds_not
 from pygeodesy.fmath import cbrt, cbrt2, fdot, Fhorner, fpowers, hypot, hypot_, \
                             hypot1, hypot2, sqrt3,  Fsum
@@ -88,7 +88,7 @@ from pygeodesy.namedTuples import Distance2Tuple, Vector3Tuple, Vector4Tuple
 from pygeodesy.props import deprecated_Property_RO, Property_RO, property_doc_, \
                             deprecated_property_RO, property_RO, property_ROver
 from pygeodesy.streprs import Fmt, fstr, instr, strs, unstr
-# from pygeodesy.triaxials import _hartzell3  # _MODS
+# from pygeodesy.triaxials.triaxial5 import _hartzell3, _plumbTo3  # _MODS
 from pygeodesy.units import Azimuth, Bearing, Distance, Float, Float_, Height, Lamd, Lat, \
                             Meter, Meter2, Meter3, Phi, Phid, Radius, Radius_, Scalar
 from pygeodesy.utily import atan1, atan1d, atan2b, degrees90, m2radians, radians2m, sincos2d
@@ -96,11 +96,12 @@ 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  # as _tan
 
 __all__ = _ALL_LAZY.ellipsoids
-__version__ = '25.08.31'
+__version__ = '25.12.14'
 
 _f_0_0    = Float(f =_0_0)  # zero flattening
 _f__0_0   = Float(f_=_0_0)  # zero inverse flattening
 # see U{WGS84_f<https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1Constants.html>}
+# _f_WGS84  = Float(f =_1_0 / (298257223563 / 1000000000))  # 0.003_352_810_664_747_480_5
 _f__WGS84 = Float(f_=_1_0 / (1000000000 / 298257223563))  # 298.257223562_999_97 vs 298.257223563
 
 
@@ -110,19 +111,19 @@ def _aux(lat, inverse, auxLat, clip=90):
     return Lat(lat, clip=clip, name=_lat_ if inverse else typename(auxLat))
 
 
-def _s2_c2(phi):
+def _sin2cos2(rad):
     '''(INTERNAL) Return 2-tuple C{(sin(B{phi})**2, cos(B{phi})**2)}.
     '''
-    if phi:
-        s2 = sin(phi)**2
+    if rad:
+        s2 = sin(rad)**2
         if s2 > EPS:
             c2 = _1_0 - s2
             if c2 > EPS:
                 if c2 < EPS1:
                     return s2, c2
             else:
-                return _1_0, _0_0  # phi == PI_2
-    return _0_0, _1_0  # phi == 0
+                return _1_0, _0_0  # rad == PI_2
+    return _0_0, _1_0  # rad == 0
 
 
 class a_f2Tuple(_NamedTuple):
@@ -870,6 +871,12 @@ class Ellipsoid(_NamedEnumItem):
 
     equatoradius = a  # Requatorial
 
+    @Property_RO
+    def equatorimeter(self):
+        '''Get the ellipsoid's I{equatorial} perimeter (C{meter}).
+        '''
+        return Meter(equatorimeter=self.a * PI2)
+
     def e2s(self, s):
         '''Compute norm M{sqrt(1 - e2 * s**2)}.
 
@@ -1133,7 +1140,7 @@ class Ellipsoid(_NamedEnumItem):
                  methods L{Ellipsoid.height4} and L{Triaxial.hartzell4}.
         '''
         try:
-            v, d, i = _MODS.triaxials._hartzell3(pov, los, self._triaxial)
+            v, d, i = _MODS.triaxials.triaxial5._hartzell3(pov, los, self._triaxial)
         except Exception as x:
             raise IntersectionError(pov=pov, los=los, cause=x)
         return Vector4Tuple(v.x, v.y, v.z, d, iteration=i, name__=self.hartzell4)
@@ -1187,7 +1194,7 @@ class Ellipsoid(_NamedEnumItem):
                 v = v.times_(t, t, 0)  # force z=0.0
                 h = x - a  # equatorial
             else:  # normal in 1st quadrant
-                x, y, i = _MODS.triaxials._plumbTo3(x, y, self)
+                x, y, i = _MODS.triaxials.triaxial5._plumbTo3(x, y, self)
                 t, v = v, v.times_(x, x, y)
                 h = t.minus(v).length
 
@@ -1380,6 +1387,12 @@ class Ellipsoid(_NamedEnumItem):
 
     polaradius = b  # Rpolar
 
+    @property_RO
+    def polarimeter(self):
+        '''Get the ellipsoid's I{polar}, meridional perimeter (C{meter}).
+        '''
+        return Meter(polarimeter=self.L * _4_0)
+
 #   Q = A  # I{meridian arc unit} C{Q}, the mean, meridional length I{per radian}
 
     @deprecated_Property_RO
@@ -1483,7 +1496,7 @@ class Ellipsoid(_NamedEnumItem):
         r, p = self.a, Phid(lat)
         if p and self.f:
             if fabs(p) < PI_2:
-                s2, c2 = _s2_c2(p)
+                s2, c2 = _sin2cos2(p)
                 # R == sqrt((a2**2 * c2 + b2**2 * s2) / (a2 * c2 + b2 * s2))
                 #   == sqrt(a2**2 * (c2 + (b2 / a2)**2 * s2) / (a2 * (c2 + b2 / a2 * s2)))
                 #   == sqrt(a2 * (c2 + (b2 / a2)**2 * s2) / (c2 + (b2 / a2) * s2))
@@ -1708,7 +1721,7 @@ class Ellipsoid(_NamedEnumItem):
         '''(INTERNAL) Helper for C{rocAzimuth} and C{rocBearing}.
         '''
         if self.f:
-            s2, c2 = _s2_c2(radians(deg))
+            s2, c2 = _sin2cos2(radians(deg))
             m, n = self.roc2_(Phid(lat))
             if n < m:  # == n / (c2 * n / m + s2)
                 c2 *= n / m
@@ -1750,7 +1763,7 @@ class Ellipsoid(_NamedEnumItem):
         # ... requires 1 or 2 sqrt
         g = self.b
         if self.f:
-            s2, c2 = _s2_c2(Phid(lat))
+            s2, c2 = _sin2cos2(Phid(lat))
             g = _over(g, c2 + self.b2_a2 * s2)
         return Radius(rocGauss=g)
 
@@ -2411,7 +2424,7 @@ Ellipsoids._assert(  # <https://WikiPedia.org/wiki/Earth_ellipsoid>
     WGS60          = _lazy('WGS60',          *_T(6378165.0,    6356783.28695944,  298.3)),
     WGS66          = _lazy('WGS66',          *_T(6378145.0,    6356759.76948868,  298.25)),
     WGS72          = _lazy(_WGS72_,          *_T(6378135.0,   _0_0,               298.26)),  # b=6356750.52
-    WGS84          = _lazy(_WGS84_,          *_T(R_MA,        _0_0,           _f__WGS84)),  # GPS b=6356752.3142451793
+    WGS84          = _lazy(_WGS84_,          *_T(R_MA,        _0_0,           _f__WGS84)),   # b=6356752.3142451793
 #   U{NOAA/NOS/NGS/inverse<https://GitHub.com/noaa-ngs/inverse/blob/main/invers3d.f>}
     WGS84_NGS      = _lazy('WGS84_NGS',      *_T(R_MA,        _0_0,               298.257222100882711243162836600094))
 )
@@ -2475,7 +2488,7 @@ if __name__ == _DMAIN_:
 
 # **) MIT License
 #
-# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2026 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/elliptic.py

@@ -1269,7 +1269,7 @@ class _RJfma(object):
 
 # **) MIT License
 #
-# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2026 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/epsg.py

@@ -223,7 +223,7 @@ __all__ += _ALL_DOCS(decode2, encode)
 
 # **) MIT License
 #
-# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2026 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/errors.py

@@ -28,7 +28,7 @@ from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS, _PYTHON_X_DEV
 from copy import copy as _copy
 
 __all__ = _ALL_LAZY.errors  # _ALL_DOCS('_InvalidError', '_IsnotError')  _under
-__version__ = '25.10.30'
+__version__ = '25.12.02'
 
 _argument_   = 'argument'
 _basics      = _MODS.into(basics=__name__)
@@ -214,7 +214,8 @@ class CrossError(_ValueError):
 
 class GeodesicError(_ValueError):
     '''Error raised for convergence or other issues in L{geodesicx<pygeodesy.geodesicx>},
-       L{geodesicw<pygeodesy.geodesicw>} or L{karney<pygeodesy.karney>}.
+       L{geodesicw<pygeodesy.geodesicw>}, L{geodsolve<pygeodesy.geodsolve>} or
+       L{karney<pygeodesy.karney>}.
     '''
     pass
 
@@ -626,7 +627,10 @@ def _xattr(obj, **name_default):
     '''
     if len(name_default) == 1:
         for n, d in name_default.items():
-            return getattr(obj, n, d)
+            try:  # obj is tuple
+                return getattr(obj, n, d)
+            except ValueError:
+                return d
     raise _xAssertionError(_xattr, obj, **name_default)
 
 
@@ -928,7 +932,7 @@ def _xStrError(*Refs, **name_value_Error):  # in .gars, .geohash, .wgrs
 
 # **) MIT License
 #
-# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2026 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/etm.py

@@ -1197,7 +1197,7 @@ if __name__ == _DMAIN_:  # MCCABE 16
 
 # **) MIT License
 #
-# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2026 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/fmath.py

@@ -10,8 +10,8 @@ from __future__ import division as _; del _  # noqa: E702 ;
 from pygeodesy.basics import _copysign, copysign0, isbool, isint, isodd, \
                               isscalar, len2, map1, _xiterable,  typename
 from pygeodesy.constants import EPS0, EPS02, EPS1, NAN, PI, PI_2, PI_4, \
-                               _0_0, _0_125, _0_25, _0_5, _1_0, _1_5, \
-                               _copysign_0_0, isfinite, remainder
+                               _0_0, _0_125, _0_25, _1_3rd, _0_5, _2_3rd, \
+                               _1_0, _1_5, _copysign_0_0, isfinite, remainder
 from pygeodesy.errors import _IsnotError, LenError, _TypeError, _ValueError, \
                              _xError, _xkwds, _xkwds_pop2, _xsError
 from pygeodesy.fsums import _2float, Fsum, fsum, _isFsum_2Tuple,  Fmt, unstr
@@ -25,13 +25,11 @@ from math import fabs, sqrt  # pow
 import operator as _operator  # in .datums, .trf, .utm
 
 __all__ = _ALL_LAZY.fmath
-__version__ = '25.09.15'
+__version__ = '25.12.23'
 
 # sqrt(2) - 1 <https://WikiPedia.org/wiki/Square_root_of_2>
 _0_4142  =  0.41421356237309504880  # ~ 3_730_904_090_310_553 / 9_007_199_254_740_992
-_1_3rd   = _1_0 / 3
-_1_6th   = _1_0 / 6
-_2_3rd   = _1_3rd * 2
+_1_6th   =  1 / 6
 _h_lt_b_ = 'abs(h) < abs(b)'
 
 
@@ -260,7 +258,6 @@ def bqrt(x):
 
 try:
     from math import cbrt as _cbrt  # Python 3.11+
-
 except ImportError:  # Python 3.10-
 
     def _cbrt(x):
@@ -314,8 +311,9 @@ def euclid(x, y):
        @see: Function L{euclid_}.
     '''
     x, y = abs(x), abs(y)  # NOT fabs!
-    return (x + y * _0_4142) if x > y else \
-           (y + x * _0_4142)  # * _0_5 before 20.10.02
+    if x < y:
+        x, y = y, x
+    return x + y * _0_4142  # * _0_5 before 20.10.02
 
 
 def euclid_(*xs):
@@ -332,10 +330,17 @@ def euclid_(*xs):
     e = _0_0
     for x in sorted(map(abs, xs)):  # NOT fabs, reverse=True!
         # e = euclid(x, e)
-        if e < x:
-            e, x = x, e
         if x:
-            e += x * _0_4142
+            if e < x:
+                e, x = x, e
+            x *= _0_4142
+#           s  = e + x
+#           if e < x:  # like .fsums._2sum
+#               x -= s  # e = (x - s) + e
+#           else:
+#               e -= s  # e = (e - s) + x
+            e += x
+#           e += s
     return e
 
 
@@ -1042,7 +1047,7 @@ def sqrt3(x):
 
 
 def sqrt_a(h, b):
-    '''Compute C{I{a}} side of a right-angled triangle from
+    '''Compute the C{I{a}} side of a right-angled triangle from
        C{sqrt(B{h}**2 - B{b}**2)}.
 
        @arg h: Hypotenuse or outer annulus radius (C{scalar}).
@@ -1115,7 +1120,7 @@ def zqrt(x):
 
 # **) MIT License
 #
-# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2026 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),