pygeodesy 25.10.10__py2.py3-none-any.whl → 25.12.12__py2.py3-none-any.whl

pygeodesy/units.py

@@ -820,6 +820,11 @@ def _isMeter(obj, iscalar=True):
     return isinstance(obj, _Meters) or (iscalar and _isScalar(obj))
 
 
+def _isRadians(obj, iscalar=True):
+    # Check for valid radian types.
+    return isinstance(obj, _Radians) or (iscalar and _isScalar(obj))
+
+
 def _isRadius(obj, iscalar=True):
     # Check for valid earth radius types.
     return isinstance(obj, _Radii) or (iscalar and _isScalar(obj))
@@ -872,7 +877,7 @@ __all__ += _ALL_DOCS(_NamedUnit)
 
 # **) 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/unitsBase.py

@@ -350,7 +350,7 @@ __all__ += _ALL_DOCS(_NamedUnit)
 
 # **) 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/ups.py

@@ -22,7 +22,7 @@ other than C{"std"}, zones C{'A'} and C{'Y'} are used for negative, west longitu
 '''
 
 # from pygeodesy.basics import neg as _neg  # from .dms
-from pygeodesy.constants import EPS, EPS0, _EPSmin as _Tol90, \
+from pygeodesy.constants import EPS, EPS0, _EPSmin as _Tol90, _K0_UPS, \
                                 isnear90, _0_0, _0_5, _1_0, _2_0
 from pygeodesy.datums import _ellipsoidal_datum, _WGS84
 from pygeodesy.dms import degDMS, _neg, parseDMS2
@@ -53,7 +53,6 @@ __version__ = '25.04.14'
 
 _BZ_UPS  = _envPYGEODESY('UPS_POLES', _std_) == _std_
 _Falsing =  Meter(2000e3)  # false easting and northing (C{meter})
-_K0_UPS  =  Float(_K0_UPS= 0.994)  # scale factor at central meridian
 _K1_UPS  =  Float(_K1_UPS=_1_0)    # rescale point scale factor
 
 
@@ -509,7 +508,7 @@ def upsZoneBand5(lat, lon, strict=True, **name):
 
 # **) 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/utily.py

@@ -16,7 +16,7 @@ from pygeodesy.constants import EPS, EPS0, NAN, PI, PI2, PI_2, PI_4, PI_6, R_M,
                                _M_KM, _M_NM, _M_SM, _0_0, _0_5, _1_0, _N_1_0, \
                                _10_0, _90_0, _180_0, _360_0, _copysign_0_0, \
                                _copysignINF, _float, _isfinite, isnan, isnear0, \
-                               _over_1, _umod_360, _umod_PI2
+                               _over_1, _umod_360, _umod_PI2, OVERFLOW
 from pygeodesy.errors import _ValueError, _xkwds, _xkwds_get
 from pygeodesy.internals import _Enum, _passargs, typename
 from pygeodesy.interns import _edge_, _radians_, _semi_circular_, _SPACE_
@@ -28,7 +28,7 @@ from math import acos, asin, asinh, atan2 as _atan2, cos, degrees, fabs, \
                  radians, sin, sinh, tan as _tan  # pow
 
 __all__ = _ALL_LAZY.utily
-__version__ = '25.09.09'
+__version__ = '25.11.09'
 
 # sqrt(3) <https://WikiPedia.org/wiki/Square_root_of_3>
 _COS_30,  _SIN_30 = 0.86602540378443864676, _0_5  # sqrt(3) / 2
@@ -807,15 +807,16 @@ def _sin0cos2(q, r, sign):
     return s, c
 
 
-def SinCos2(x):
+def SinCos2(x, unit=Radians):
     '''Get C{sin} and C{cos} of I{typed} angle.
 
        @arg x: Angle (L{Degrees}, L{Radians} or scalar C{radians}).
+       @kwarg unit: The C{B{x}} unit (L{Degrees}, L{Radians}).
 
        @return: 2-Tuple (C{sin(B{x})}, C{cos(B{x})}).
     '''
-    return sincos2d(x) if isinstanceof(x, Degrees, Degrees_) else (
-#          sincos2(x)  if isinstanceof(x, Radians, Radians_) else
+    return sincos2d(x) if unit is Degrees or unit is degrees or isinstanceof(x, Degrees, Degrees_) else (
+#          sincos2(x)  if unit is Radians or unit is radians or isinstanceof(x, Radians, Radians_) else
            sincos2(Radians(x)))  # assume C{radians}
 
 
@@ -996,13 +997,14 @@ def tan_(*rads, **raiser_kwds):
         yield _nonfinite(tan_, r, **raiser_kwds)
 
 
-def tand(deg, **raiser_kwds):
+def tand(deg, **raiser_clamp_kwds):
     '''Return the C{tangent} of an angle in C{degrees}.
 
        @arg deg: Angle (C{degrees}).
-       @kwarg raiser_kwds: Use C{B{raiser}=False} to avoid
-                     ValueErrors and optional, additional
-                     ValueError keyword argments.
+       @kwarg raiser_clamp_kwds: Use C{B{raiser}=False} to avoid
+                     ValueErrors, C{B{clamp}=}L{OVERFLOW} and
+                     optional, additional ValueError keyword
+                     argments.
 
        @return: C{tan(B{deg})}.
 
@@ -1011,10 +1013,10 @@ def tand(deg, **raiser_kwds):
     try:
         return _tanu(*sincos2d(deg))
     except ZeroDivisionError:
-        return _nonfinite(tand, deg, **raiser_kwds)
+        return _nonfinite(tand, deg, **raiser_clamp_kwds)
 
 
-def tand_(*degs, **raiser_kwds):
+def tand_(*degs, **raiser_clamp_kwds):
     '''Yield the C{tangent} of angle(s) in C{degrees}.
 
        @arg degs: One or more angles (each in C{degrees}).
@@ -1025,9 +1027,9 @@ def tand_(*degs, **raiser_kwds):
     '''
     try:
         for d in degs:
-            yield _tanu(*sincos2d(d))
+            yield _tanu(*sincos2d(d), **raiser_clamp_kwds)
     except ZeroDivisionError:
-        yield _nonfinite(tand_, d, **raiser_kwds)
+        yield _nonfinite(tand_, d, **raiser_clamp_kwds)
 
 
 def tanPI_2_2(rad):
@@ -1041,7 +1043,7 @@ def tanPI_2_2(rad):
             NAN if isnan(rad) else _copysign(_90_0, rad))
 
 
-def _tanu(s, c):
+def _tanu(s, c, clamp=OVERFLOW, **unused):
     '''(INTERNAL) Helper for functions C{_cotu}, C{sincostan3},
        C{sincostan3d}, C{tan}, C{tan_}, C{tand} and C{tand_}.
     '''
@@ -1051,6 +1053,8 @@ def _tanu(s, c):
         raise ZeroDivisionError()
     else:
         t = _over_1(s, c)
+        if clamp:
+            t = min(clamp, max(t, -clamp))
     return t
 
 
@@ -1383,7 +1387,7 @@ def yard2m(yards):
 
 # **) 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/utm.py

@@ -733,7 +733,7 @@ def utmZoneBand5(lat, lon, cmoff=False, **name):
 
 # **) 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/utmups.py

@@ -294,7 +294,7 @@ def utmupsZoneBand5(lat, lon, cmoff=False, **name):
 
 # **) 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/utmupsBase.py

@@ -586,7 +586,7 @@ __all__ += _ALL_DOCS(UtmUpsBase)
 
 # **) 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/vector2d.py

@@ -641,7 +641,7 @@ def triaxum5(points, useZ=True):
     A = []
     for i, p in enumerate(ps):
         v = _otherV3d(useZ=useZ, i=i, points=p)
-        A.append(v.x2y2z2)
+        A.append(v.x2y2z23)
 
     with _numpy(triaxum5, n=n) as _np:
         A = _np.array(A)
@@ -884,7 +884,7 @@ def _tri5perturbs(eps, r):
 
 # **) 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/vector3d.py

@@ -914,7 +914,7 @@ __all__ += _ALL_DOCS(intersections2, sumOf, Vector3dBase)
 
 # **) 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/vector3dBase.py

@@ -10,9 +10,9 @@ A pure Python implementation of vector-based functions by I{(C) Chris Veness
 
 from pygeodesy.basics import _copysign, islistuple, isscalar, \
                               map1, map2, _signOf, _zip
-from pygeodesy.constants import EPS, EPS0, INT0, PI, PI2, \
-                               _1_0, isnear0, isnear1, isneg0, \
-                               _copysignINF, _float0, _pos_self
+from pygeodesy.constants import EPS, EPS0, INT0, PI, PI2, _1_0, \
+                               _copysignINF, isnear0, isnear1, \
+                               _flipsign, _float0, _pos_self
 from pygeodesy.errors import CrossError, VectorError, _xcallable, _xError
 from pygeodesy.fmath import euclid_, fdot, fdot_, hypot_, hypot2_  # _MODS.fmath.fma
 from pygeodesy.interns import _coincident_, _colinear_, _COMMASPACE_, _xyz_
@@ -29,7 +29,7 @@ from pygeodesy.utily import atan2, sincos2,  fabs
 from math import ceil as _ceil, floor as _floor, trunc as _trunc
 
 __all__ = _ALL_LAZY.vector3dBase
-__version__ = '25.08.31'
+__version__ = '25.10.25'
 
 
 class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
@@ -112,9 +112,20 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
 
     cmp = __cmp__
 
-    def __divmod__(self, other):  # PYCHOK no cover
-        '''Not implemented.'''
-        return _NotImplemented(self, other)
+    def __divmod__(self, scalar):  # PYCHOK no cover
+        '''Apply C{scalar} divmod to this vector's components.
+
+           @arg scalar: Divisor (C{scalar}).
+
+           @return: 2-Tuple C{(d, m)} each a L{Vector3d}.
+        '''
+        s = Scalar(divisor=scalar)
+        d, m = [], []
+        for x in self.xyz3:
+            q, r = divmod(x, s)
+            d.append(q)
+            m.append(r)
+        return self.classof(d), self.classof(m)
 
     def __eq__(self, other):
         '''Is this vector equal to an other vector?
@@ -132,15 +143,22 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         return _NotImplemented(self)  # must return C{float}
 
     def __floor__(self):  # PYCHOK no cover
-        '''Return a vector with the C{floor} of these components.
+        '''Return a vector with the C{floor} of each components.
 
            @return: Floor-ed (L{Vector3d}).
         '''
         return self._mapped(_floor)
 
-    def __floordiv__(self, other):  # PYCHOK no cover
-        '''Not implemented.'''
-        return _NotImplemented(self, other)
+    def __floordiv__(self, scalar):  # PYCHOK no cover
+        '''Floor-divide this vector by a scalar, C{this // B{scalar}}.
+
+           @arg scalar: The divisor (C{scalar}).
+
+           @return: Floored quotient (L{Vector3d}).
+
+           @raise TypeError: Non-scalar B{C{scalar}}.
+        '''
+        return self.floorDividedBy(scalar)
 
     def __ge__(self, other):
         '''Is this vector longer than or equal to an other vector?
@@ -183,9 +201,14 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self._xyz(self.plus(other))
 
-    def __ifloordiv__(self, other):  # PYCHOK no cover
-        '''Not implemented.'''
-        return _NotImplemented(self, other)
+    def __ifloordiv__(self, scalar):  # PYCHOK no cover
+        '''Floor-divide this vector by a scalar I{in-place}, C{this //= B{scalar}}.
+
+           @arg scalar: The divisor (C{scalar}).
+
+           @raise TypeError: Non-scalar B{C{scalar}}.
+        '''
+        return self._xyz(self.floorDividedBy(scalar))
 
     def __imatmul__(self, other):  # PYCHOK Python 3.5+
         '''Cross multiply this and an other vector I{in-place}, C{this @= B{other}}.
@@ -215,9 +238,16 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''Not implemented, see method C{ints}.'''
         return _NotImplemented(self)  # must return C{int}
 
-    def __ipow__(self, other, *mod):  # PYCHOK no cover
-        '''Not implemented.'''
-        return _NotImplemented(self, other, *mod)
+    def __ipow__(self, scalar, *mod):  # PYCHOK no cover
+        '''Raise each component I{in-place} as C{pow(C, B{scalar})}.
+
+           @arg scalar: Exponent (C{scalar}).
+
+           @return: Power (L{Vector3d}).
+
+           @raise TypeError: Non-scalar B{C{scalar}}.
+        '''
+        return self._xyz(self.pow(scalar, *mod))
 
     def __isub__(self, other):
         '''Subtract an other vector from this one I{in-place}, C{this -= B{other}}.
@@ -280,9 +310,19 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self.cross(other)
 
-    def __mod__(self, other):  # PYCHOK no cover
-        '''Not implemented.'''
-        return _NotImplemented(self, other)
+    def __mod__(self, modulus):  # PYCHOK no cover
+        '''Apply C{scalar} modulo to this vector's components.
+
+           @arg modulus: Modulus (C{scalar}).
+
+           @return: Modulo (L{Vector3d}).
+        '''
+        m = Scalar(modulus=modulus)
+
+        def _mod(x):
+            return x % m
+
+        return self._mapped(_mod)
 
     def __mul__(self, scalar):
         '''Multiply this vector by a scalar, C{this * B{scalar}}.
@@ -318,13 +358,29 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self if _pos_self else self.copy()
 
-    def __pow__(self, other, *mod):  # PYCHOK no cover
-        '''Not implemented.'''
-        return _NotImplemented(self, other, *mod)
+    def __pow__(self, scalar, *mod):  # PYCHOK no cover
+        '''Return a vector with components as C{pow(C, B{scalar})}.
+
+           @arg scalar: Exponent (C{scalar}).
+
+           @return: Power (L{Vector3d}).
+
+           @raise TypeError: Non-scalar B{C{scalar}}.
+        '''
+        return self.pow(scalar, *mod)
 
-    __radd__ = __add__  # PYCHOK no cover
+    def __radd__(self, other):  # PYCHOK no cover
+        '''Add this vector to an other vector, C{B{other} + this}.
+
+           @arg other: The other vector (L{Vector3d}).
+
+           @return: Sum (L{Vector3d}).
+
+           @raise TypeError: Incompatible B{C{other}} C{type}.
+        '''
+        return self.others(other).plus(self)
 
-    def __rdivmod__ (self, other):  # PYCHOK no cover
+    def __rdivmod__(self, other):
         '''Not implemented.'''
         return _NotImplemented(self, other)
 
@@ -333,7 +389,7 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
 #       '''
 #       return self.toRepr()
 
-    def __rfloordiv__(self, other):  # PYCHOK no cover
+    def __rfloordiv__(self, other):
         '''Not implemented.'''
         return _NotImplemented(self, other)
 
@@ -348,11 +404,11 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self.others(other).cross(self)
 
-    def __rmod__(self, other):  # PYCHOK no cover
+    def __rmod__(self, other):
         '''Not implemented.'''
         return _NotImplemented(self, other)
 
-    __rmul__ = __mul__
+    __rmul__ = __mul__  # scalar other
 
     def __round__(self, *ndigits):  # PYCHOK no cover
         '''Return a vector with these components C{rounded}.
@@ -361,8 +417,11 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
 
            @return: Rounded (L{Vector3d}).
         '''
+        def _rnd(x):
+            return round(x, *ndigits)
+
         # <https://docs.Python.org/3.12/reference/datamodel.html?#object.__round__>
-        return self.classof(*(round(_, *ndigits) for _ in self.xyz3))
+        return self._mapped(_rnd)
 
     def __rpow__(self, other, *mod):  # PYCHOK no cover
         '''Not implemented.'''
@@ -379,9 +438,9 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self.others(other).minus(self)
 
-    def __rtruediv__(self, scalar):  # PYCHOK no cover
+    def __rtruediv__(self, other):
         '''Not implemented.'''
-        return _NotImplemented(self, scalar)
+        return _NotImplemented(self, other)
 
 #   def __str__(self):
 #       '''Return the default C{str(self)}.
@@ -527,22 +586,39 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         self._crosserrors = bool(raiser)
 
-    def dividedBy(self, divisor):
+    def dividedBy(self, scalar):
         '''Divide this vector by a scalar.
 
-           @arg divisor: The divisor (C{scalar}).
+           @arg scalar: The divisor (C{scalar}).
 
            @return: New, scaled vector (L{Vector3d}).
 
-           @raise TypeError: Non-scalar B{C{divisor}}.
+           @raise TypeError: Non-scalar B{C{scalar}}.
 
-           @raise VectorError: Invalid or zero B{C{divisor}}.
+           @raise VectorError: Invalid or zero B{C{scalar}}.
         '''
-        d = Scalar(divisor=divisor)
+        d = Scalar(divisor=scalar)
         try:
             return self._times(_1_0 / d)
         except (ValueError, ZeroDivisionError) as x:
-            raise VectorError(divisor=divisor, cause=x)
+            raise VectorError(divisor=scalar, cause=x)
+
+    def dividedBy_(self, scalar_x, *y_z):
+        '''Divide this vector by separate X, Y and Z divisors.
+
+           @arg scalar_x: X divisor (C{scalar}) or a vector's
+                          X, Y, and Z components as divisors
+                          (C{Cartesian}, L{Ecef9Tuple}, C{Nvector},
+                          L{Vector3d}, L{Vector3Tuple}, L{Vector4Tuple}).
+           @arg y_z: Y and Z divisors (C{scalar}, C{scalar}), ignored
+                     if B{C{scalar_x}} is not C{scalar}.
+
+           @return: New, scaled vector (L{Vector3d}).
+
+           @raise VectorError: Invalid B{C{scalar_x}} or B{C{y_z}}.
+        '''
+        x, y, z = _xyz3(self.dividedBy_, scalar_x, *y_z)
+        return self.classof(self.x / x, self.y / y, self.z / z)
 
     def dot(self, other):
         '''Compute the dot (scalar) product of this and an other vector.
@@ -595,6 +671,43 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self._mapped(_float0)
 
+    def floorDividedBy(self, scalar):
+        '''Floor-divide this vector by a scalar.
+
+           @arg scalar: The divisor (C{scalar}).
+
+           @return: New, scaled vector (L{Vector3d}).
+
+           @raise TypeError: Non-scalar B{C{scalar}}.
+
+           @raise VectorError: Invalid or zero B{C{scalar}}.
+        '''
+        d = Scalar(divisor=scalar)
+        try:
+            def _floor_d(x):
+                return x // d
+
+            return self._mapped(_floor_d)
+        except (ValueError, ZeroDivisionError) as x:
+            raise VectorError(divisor=scalar, cause=x)
+
+    def floorDividedBy_(self, scalar_x, *y_z):
+        '''Floor-divide this vector by separate X, Y and Z divisors.
+
+           @arg scalar_x: X divisor (C{scalar}) or a vector's
+                          X, Y, and Z components as divisors
+                          (C{Cartesian}, L{Ecef9Tuple}, C{Nvector},
+                          L{Vector3d}, L{Vector3Tuple}, L{Vector4Tuple}).
+           @arg y_z: Y and Z divisors (C{scalar}, C{scalar}), ignored
+                     if B{C{scalar_x}} is not C{scalar}.
+
+           @return: New, scaled vector (L{Vector3d}).
+
+           @raise VectorError: Invalid B{C{scalar_x}} or B{C{y_z}}.
+        '''
+        x, y, z = _xyz3(self.floorDividedBy_, scalar_x, *y_z)
+        return self.classof(self.x // x, self.y // y, self.z // z)
+
     @Property
     def _fromll(self):
         '''(INTERNAL) Get the latlon reference (C{LatLon}) or C{None}.
@@ -617,11 +730,8 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
             y =  y / z
 #           z = _1_0
         else:
-            if isneg0(z):
-                x = -x
-                y = -y
-            x = _copysignINF(x)
-            y = _copysignINF(y)
+            x = _copysignINF(_flipsign(x, z))
+            y = _copysignINF(_flipsign(y, z))
 #           z =  NAN
         return self.classof(x, y, _1_0)
 
@@ -705,7 +815,7 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
 
            @see: Properties L{length2} and L{euclid}.
         '''
-        return Float(length=hypot_(self.x, self.y, self.z))
+        return Float(length=hypot_(*self.xyz3))
 
     @Property_RO
     def length2(self):  # __dict__ value overwritten by Property_RO C{_united}
@@ -713,7 +823,7 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
 
            @see: Property L{length} and method C{equirectangular}.
         '''
-        return Float(length2=hypot2_(self.x, self.y, self.z))
+        return Float(length2=hypot2_(*self.xyz3))
 
     def _mapped(self, func):
         '''(INTERNAL) Apply C{func} to all components.
@@ -815,12 +925,28 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
            @arg y_z: Y and Z components (C{scalar}, C{scalar}),
                      ignored if B{C{other_x}} is not C{scalar}.
 
-           @return: New, vectiorial vector (L{Vector3d}).
+           @return: New, vectorial vector (L{Vector3d}).
 
            @raise ValueError: Invalid B{C{other_x}} or B{C{y_z}}.
         '''
         return self._plus(*_xyz3(self.plus_, other_x, *y_z))
 
+    def pow(self, scalar, *mod):
+        '''Raise each X, Y and Z to C{scalar} power.
+
+           @arg scalar: Exponent (C{scalar}).
+
+           @return: Power (L{Vector3d}).
+
+           @raise TypeError: Non-scalar B{C{scalar}}.
+        '''
+        p = Scalar(power=scalar)
+
+        def _pow(x):
+            return pow(x, p, *mod)
+
+        return self._mapped(_pow)
+
     def rotate(self, axis, theta, fma=False):
         '''Rotate this vector around an axis by a specified angle.
 
@@ -863,6 +989,16 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''DEPRECATED, use method C{rotate}.'''
         return self.rotate(axis, theta)  # PYCHOK no cover
 
+    def _roty(self, *pos, **name):
+        '''(INTERNAL) Prolate rotation, for C{+1 if pos else -1}.
+        '''
+        v = self.copy(**name)
+        if pos:
+            x, _, z = v.xyz3
+            _update_all(v, needed=3)
+            v._x, v._z = (-z, x) if pos[0] else (z, -x)
+        return v
+
     def times(self, factor):
         '''Multiply this vector by a scalar.
 
@@ -989,8 +1125,10 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
     def _xyz(self, x_xyz, *y_z):
         '''(INTERNAL) Set the C{_x}, C{_y} and C{_z} attributes.
         '''
-        _update_all(self, needed=3)
-        self._x, self._y, self._z = _xyz3(_xyz_, x_xyz, *y_z)
+        xyz = _xyz3(_xyz_, x_xyz, *y_z)
+        if self.xyz3 != xyz:
+            _update_all(self, needed=3)
+            self._x, self._y, self._z = xyz
         return self
 
     @property_RO
@@ -999,9 +1137,15 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
         '''
         return self.x, self.y, self.z
 
-    @property_RO
+    @Property_RO
     def x2y2z2(self):
-        '''Get the X, Y and Z components I{squared} (3-tuple C{(x**2, y**2, z**2)}).
+        '''Get the X, Y and Z components, I{squared} (L{Vector3Tuple}).
+        '''
+        return _MODS.namedTuples.Vector3Tuple(*self.x2y2z23, name=self.name)
+
+    @property_RO
+    def x2y2z23(self):
+        '''Get the X, Y and Z components, I{squared} (3-tuple C{(x**2, y**2, z**2)}).
         '''
         return self.x**2, self.y**2, self.z**2
 
@@ -1037,7 +1181,7 @@ class Vector3dBase(_NamedBase):  # sync __methods__ with .fsums.Fsum
 
 
 def _xyz3(where, x_xyz, *y_z):  # in .cartesianBase._rtp3
-    '''(INTERNAL) Get , Y and Z as 3-tuple C{(x, y, z)}.
+    '''(INTERNAL) Get X, Y and Z as 3-tuple C{(x, y, z)}.
     '''
     try:
         xyz3 = map1(_float0, x_xyz, *y_z) if y_z else (  # islistuple for Vector*Tuple
@@ -1052,7 +1196,7 @@ __all__ += _ALL_DOCS(Vector3dBase)
 
 # **) 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/webmercator.py

@@ -353,7 +353,7 @@ def toWm(latlon, lon=None, earth=R_MA, Wm=Wm, **name_Wm_kwds_radius):
 
 # **) 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/wgrs.py

@@ -425,7 +425,7 @@ __all__ += _ALL_DOCS(decode3, decode5,  # functions
 
 # **) 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"),