pygeodesy 24.10.24__py2.py3-none-any.whl → 24.12.12__py2.py3-none-any.whl

pygeodesy/fsums.py

@@ -64,7 +64,7 @@ from math import fabs, isinf, isnan, \
                  ceil as _ceil, floor as _floor  # PYCHOK used! .ltp
 
 __all__ = _ALL_LAZY.fsums
-__version__ = '24.10.22'
+__version__ = '24.12.02'
 
 from pygeodesy.interns import (
   _PLUS_     as _add_op_,  # in .auxilats.auxAngle
@@ -100,7 +100,7 @@ def _2finite(x, _isfine=_isfinite):  # in .fstats
 
 
 def _2float(index=None, _isfine=_isfinite, **name_x):  # in .fmath, .fstats
-    '''(INTERNAL) Raise C{TypeError} or C{Overflow-/ValueError} if not finite.
+    '''(INTERNAL) Raise C{TypeError} or C{Overflow-/ValueError} if C{x} not finite.
     '''
     n, x = name_x.popitem()  # _xkwds_item2(name_x)
     try:
@@ -163,9 +163,11 @@ except  ImportError:  # PYCHOK DSPACE! Python 3.12-
 
         _2n_d = None  # redef
 
-    def _fmaX(r, *a_b_c):  # like Python 3.13+ I{Modules/mathmodule.c}:
-        # raise a ValueError for a NAN result from non-NAN C{a_b_c}s or an
-        # OverflowError for a non-NAN non-finite from all finite C{a_b_c}s.
+    def _fmaX(r, *a_b_c):  # PYCHOK no cover
+        # handle non-finite as Python 3.13+ C-function U{math_fma_impl<https://
+        # GitHub.com/python/cpython/blob/main/Modules/mathmodule.c#L2305>}:
+        #  raise a ValueError for a NAN result from non-NAN C{a_b_c}s or an
+        #  OverflowError for a non-NAN non-finite from all finite C{a_b_c}s.
         if isnan(r):
             def _x(x):
                 return not isnan(x)
@@ -212,11 +214,11 @@ except  ImportError:  # PYCHOK DSPACE! Python 3.12-
         return map(_2split3, xs)
 
 
-def f2product(*two):
+def f2product(two=None):
     '''Turn accurate I{TwoProduct} multiplication on or off.
 
-       @arg two: If C{True}, turn I{TwoProduct} on, if C{False} off or
-                 if C{None} or omitted, keep the current setting.
+       @kwarg two: If C{True}, turn I{TwoProduct} on, if C{False} off or
+                   if C{None} or omitted, keep the current setting.
 
        @return: The previous setting (C{bool}).
 
@@ -226,12 +228,12 @@ def f2product(*two):
              equivalent, slower implementation when not available.
     '''
     t = Fsum._f2product
-    if two and two[0] is not None:
-        Fsum._f2product = bool(two[0])
+    if two is not None:
+        Fsum._f2product = bool(two)
     return t
 
 
-def _Fsumf_(*xs):  # in .auxLat, .ltp, ...
+def _Fsumf_(*xs):  # in .auxLat, ...
     '''(INTERNAL) An C{Fsum(xs)}, all C{scalar}, an L{Fsum} or L{Fsum2Tuple}.
     '''
     return Fsum()._facc_scalarf(xs, up=False)
@@ -318,14 +320,14 @@ def _NonfiniteError(x):
            _ValueError    if isnan(x) else _AssertionError)
 
 
-def nonfiniterrors(*raiser):
+def nonfiniterrors(raiser=None):
     '''Throw C{OverflowError} and C{ValueError} exceptions for or
        handle I{non-finite} C{float}s as C{inf}, C{INF}, C{NINF},
        C{nan} and C{NAN} in summations and multiplications.
 
-       @arg raiser: If C{True}, throw exceptions, if C{False} handle
-                    I{non-finites} or if C{None} or omitted, leave
-                    the setting unchanged.
+       @kwarg raiser: If C{True}, throw exceptions, if C{False} handle
+                      I{non-finites} or if C{None} or omitted, leave
+                      the setting unchanged.
 
        @return: Previous setting (C{bool}).
 
@@ -333,8 +335,8 @@ def nonfiniterrors(*raiser):
               C{nan} and C{NAN} a C{ValueError}.
     '''
     d = Fsum._isfine
-    if raiser and raiser[0] is not None:
-        Fsum._isfine = {} if bool(raiser[0]) else Fsum._nonfinites_isfine_kwds[True]
+    if raiser is not None:
+        Fsum._isfine = {} if bool(raiser) else Fsum._nonfinites_isfine_kwds[True]
     return (False if d is Fsum._nonfinites_isfine_kwds[True] else
          _xkwds_get1(d, _isfine=_isfinite) is _isfinite) if d else True
 
@@ -399,7 +401,7 @@ def _residue(other):
     return r
 
 
-def _s_r(s, r):
+def _s_r2(s, r):
     '''(INTERNAL) Return C{(s, r)}, I{ordered}.
     '''
     if _isfinite(s):
@@ -413,24 +415,6 @@ def _s_r(s, r):
     return s, r
 
 
-def _2s_r(other):
-    '''(INTERNAL) Return 2-tuple C{(other, r)} with C{other} as C{int},
-       C{float} or C{as-is} and C{r} the residual of C{as-is} or 0.
-    '''
-    if _isFsum_2Tuple(other):
-        s, r = other._fint2
-        if r:
-            s, r = other._nfprs2
-            if r:  # PYCHOK no cover
-                s = other  # L{Fsum} as-is
-    else:
-        r = 0
-        s = other  # C{type} as-is
-        if isint(s, both=True):
-            s = int(s)
-    return s, r
-
-
 def _strcomplex(s, *args):
     '''(INTERNAL) C{Complex} 2- or 3-arg C{pow} error as C{str}.
     '''
@@ -487,6 +471,24 @@ def _threshold(threshold=_0_0, **kwds):
         raise ResidualError(threshold=threshold, cause=x)
 
 
+def _2tuple2(other):
+    '''(INTERNAL) Return 2-tuple C{(other, r)} with C{other} as C{int},
+       C{float} or C{as-is} and C{r} the residual of C{as-is} or 0.
+    '''
+    if _isFsum_2Tuple(other):
+        s, r = other._fint2
+        if r:
+            s, r = other._nfprs2
+            if r:  # PYCHOK no cover
+                s = other  # L{Fsum} as-is
+    else:
+        r = 0
+        s = other  # C{type} as-is
+        if isint(s, both=True):
+            s = int(s)
+    return s, r
+
+
 class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats, ...
     '''Precision floating point summation, I{running} summation and accurate multiplication.
 
@@ -645,10 +647,6 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
         f = self._copy_2(self.__floordiv__)
         return f._floordiv(other, _floordiv_op_)
 
-    def __format__(self, *other):  # PYCHOK no cover
-        '''Not implemented.'''
-        return _NotImplemented(self, *other)
-
     def __ge__(self, other):
         '''Return C{(B{self} >= B{other})}, see C{__eq__}.
         '''
@@ -1171,6 +1169,14 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
         return self._fadd(xs[0], **up) if len(xs) == 1 else \
                self._facc(xs, **up)  # origin=1?
 
+    def _facc_dot(self, n, xs, ys, **kwds):  # in .fmath
+        '''(INTERNAL) Accumulate C{fdot(B{xs}, *B{ys})}.
+        '''
+        if n > 0:
+            _f = Fsum(**kwds)
+            self._facc(_f(x).fmul(y) for x, y in zip(xs, ys))  # PYCHOK attr?
+        return self
+
     def _facc_neg(self, xs, **up_origin):
         '''(INTERNAL) Accumulate more C{xs}, negated.
         '''
@@ -1462,7 +1468,7 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
         try:
             s, r = self._fprs2
             if r:
-                f  =  self._f2mul(self.fma, other1, **nonfinites)
+                f  =  self._f2mul(self.fma, (other1,), **nonfinites)
                 f +=  other2
             elif _residue(other1) or _residue(other2):
                 fs = _2split3s(_fs(op, other1))
@@ -1501,36 +1507,41 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
     @deprecated_method
     def f2mul(self, *others, **raiser):
         '''DEPRECATED on 2024.09.13, use method L{f2mul_<Fsum.f2mul_>}.'''
-        return self._fset(self.f2mul_(*others, **raiser))
+        return self._fset(self.f2mul_(others, **raiser))
 
-    def f2mul_(self, *others, **nonfinites):  # in .fmath.f2mul
+    def f2mul_(self, *others, **f2product_nonfinites):  # in .fmath.f2mul
         '''Return C{B{self} * B{other} * B{other} ...} for all B{C{others}} using cascaded,
-           accurate multiplication like with L{f2product<Fsum.f2product>} set to C{True}.
+           accurate multiplication like with L{f2product<Fsum.f2product>}C{(B{True})}.
 
            @arg others: Multipliers (each C{scalar}, an L{Fsum} or L{Fsum2Tuple}), all
                         positional.
-           @kwarg nonfinites: Use C{B{nonfinites}=True} or C{False}, to override both
-                              L{nonfinites<Fsum.nonfinites>} and the L{nonfiniterrors}
-                              default (C{bool}).
+           @kwarg f2product_nonfinites: Use C{B{f2product=False}} to override the default
+                            C{True} and C{B{nonfinites}=True} or C{False}, to override
+                            settings L{nonfinites<Fsum.nonfinites>} and L{nonfiniterrors}.
 
            @return: The cascaded I{TwoProduct} (L{Fsum} or C{float}).
 
            @see: U{Equations 2.3<https://www.TUHH.De/ti3/paper/rump/OzOgRuOi06.pdf>}
         '''
-        return self._f2mul(self.f2mul_, *others, **nonfinites)
+        return self._f2mul(self.f2mul_, others, **f2product_nonfinites)
 
-    def _f2mul(self, where, *others, **nonfinites_raiser):
+    def _f2mul(self, where, others, f2product=True, **nonfinites_raiser):
         '''(INTERNAL) See methods C{fma} and C{f2mul_}.
         '''
-        f  = self._copy_2(where)
-        ps = f._ps
-        if ps and others:
-            op = where.__name__
-            try:
+        f = _Psum(self._ps, f2product=f2product, name=where.__name__)
+        if others and f:
+            if f.f2product():
+                def _pfs(f, ps):
+                    return _2products(f, _2split3s(ps))
+            else:
+                def _pfs(f, ps):  # PYCHOK redef
+                    return (f * p for p in ps)
+
+            op, ps = where.__name__, f._ps
+            try:  # as if self.f2product(True)
                 for other in others:  # to pinpoint errors
                     for p in self._ps_other(op, other):
-                        pfs   = _2products(p, _2split3s(ps))
-                        ps[:] =  f._ps_acc([], pfs, up=False)
+                        ps[:] = f._ps_acc([], _pfs(p, ps), update=False)
                 f._update()
             except TypeError as X:
                 raise self._ErrorX(X, op, other)
@@ -1568,7 +1579,7 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
             elif self.is_integer():
                 # return an exact C{int} for C{int}**C{int}
                 i, _ =  self._fint2  # assert _ == 0
-                x, r = _2s_r(other)  # C{int}, C{float} or other
+                x, r = _2tuple2(other)  # C{int}, C{float} or other
                 f = self._Fsum_as(i)._pow_Fsum(other, op, **raiser_RESIDUAL) if r else \
                     self._pow_2_3(i, x,        other, op, **raiser_RESIDUAL)
             else:  # mod[0] is None, power(self, other)
@@ -1634,9 +1645,9 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
 #               Fsum._ps_max = max(Fsum._ps_max, n)
                 if n > 2:
                     r = self._ps_1sum(s)
-                    return Fsum2Tuple(*_s_r(s, r))
+                    return Fsum2Tuple(*_s_r2(s, r))
             if n > 1:  # len(ps) == 2
-                s, r  = _s_r(*_2sum(*ps, **self._isfine))
+                s, r  = _s_r2(*_2sum(*ps, **self._isfine))
                 ps[:] = (r, s) if r else (s,)
             elif ps:  # len(ps) == 1
                 s = ps[0]
@@ -1662,9 +1673,8 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
     def fset_(self, *xs):
         '''Apply C{B{self}.partials = Fsum(*B{xs}).partials}.
 
-           @arg xs: Optional, new values (each C{scalar} or
-                    an L{Fsum} or L{Fsum2Tuple} instance), all
-                    positional.
+           @arg xs: Optional, new values (each C{scalar} or an L{Fsum}
+                    or L{Fsum2Tuple} instance), all positional.
 
            @return: This instance, replaced (C{Fsum}).
 
@@ -1839,7 +1849,7 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
             s, r = _facc(xs, **facc_kwds)._fprs2
             if _isfinite(s):  # _fsum(_1primed((s, -p, r, -q))
                 d, r = _2sum(s - p, r - q, _isfine=_isOK)
-                r, _ = _s_r(d, r)
+                r, _ = _s_r2(d, r)
             return s, (r if _isfinite(r) else _NONFINITEr)
         else:
             return p, _0_0
@@ -2174,7 +2184,7 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
                 return s
 
             b = _s(*(b._fprs2 if m is None else b._fint2))
-            x = _s(*_2s_r(x))
+            x = _s(*_2tuple2(x))
 
         try:
             # 0**INF == 0.0, 1**INF == 1.0, -1**2.3 == -(1**2.3)
@@ -2309,11 +2319,11 @@ class Fsum(_Named):  # sync __methods__ with .vector3dBase.Vector3dBase, .fstats
                     return (p * f for f in fs)
 
             for p in ps:
-                for f in _pfs(p, fs):
-                    yield f if _isfine(f) else _nfError(f)
+                for x in _pfs(p, fs):
+                    yield x if _isfine(x) else _nfError(x)
 
-        fs = _psfs(self._ps, factors, **self._isfine)
-        f  = _Psum(self._ps_acc([], fs, up=False), name=op)
+        xs = _psfs(self._ps, factors, **self._isfine)
+        f  = _Psum(self._ps_acc([], xs, up=False), name=op)
         return f
 
     @property_RO
@@ -2597,7 +2607,7 @@ class Fsum2Tuple(_NamedTuple):  # in .fstats
 
     @Property_RO
     def _Fsum(self):  # this C{Fsum2Tuple} as L{Fsum}, in .fstats
-        s, r = _s_r(*self)
+        s, r = _s_r2(*self)
         ps = (r, s) if r else (s,)
         return _Psum(ps, name=self.name)
 
@@ -2784,7 +2794,7 @@ def _xs(xs, _X=_X_ps, _x=float, _isfine=_isfinite,  # defaults for Fsum._facc
     i, x = 0, xs
     try:
         for i, x in enumerate(_xiterable(xs)):
-            if isinstance(x, _Fsum_2Tuple_types):
+            if _isFsum_2Tuple(x):
                 for p in _X(x):
                     yield p if _isfine(p) else _nfError(p)
             else:
@@ -2806,7 +2816,7 @@ def _xsum(which, xs, nonfinites=None, primed=0, **floats):  # origin=0
         nonfinites = _xkwds_get1(floats, floats=nonfinites)
     elif nonfinites is None:
         nonfinites =  not nonfiniterrors()
-    fs = _xs(xs, **_x_isfine(nonfinites, which=which))
+    fs = _xs(xs, **_x_isfine(nonfinites, which=which))  # PYCHOK yield
     return _fsum(_1primed(fs) if primed else fs)
 
 
@@ -2837,7 +2847,7 @@ if __name__ == '__main__':
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- 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/gars.py

@@ -345,7 +345,7 @@ __all__ += _ALL_DOCS(decode3,  # functions
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- 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/geodesici.py

@@ -51,12 +51,12 @@ from pygeodesy.solveBase import _SolveCapsBase,  pairs
 # from pygeodesy.streprs import Fmt, unstr  # from .ellipsoids
 from pygeodesy.units import Azimuth as Azi, Degrees, Float, Int, \
                            _isDegrees,  Lat, Lon, Meter, Meter_
-from pygeodesy.utily import sincos2,  atan2, fabs, radians
+from pygeodesy.utily import atan2, sincos2,  fabs, radians
 
-# from math import atan2, ceil as _ceil, fabs, radians  # .fsums, .utily
+# from math import ceil as _ceil, fabs, radians  # .fsums, .utily
 
 __all__ = _ALL_LAZY.geodesici
-__version__ = '24.10.22'
+__version__ = '24.11.24'
 
 _0t     =  0,  # int
 _1_1t   = -1, +1
@@ -1774,7 +1774,7 @@ if __name__ == '__main__':  # MCCABE 14
 
 # **) MIT License
 #
-# Copyright (C) 2024-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2024-2025 -- 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/geodesicw.py

@@ -39,7 +39,7 @@ from contextlib import contextmanager
 # from math import fabs  # from .utily
 
 __all__ = _ALL_LAZY.geodesicw
-__version__ = '24.10.14'
+__version__ = '24.11.02'
 
 _plumb_ = 'plumb'
 _TRIPS  =  65
@@ -94,14 +94,14 @@ class _gWrapped(_kWrapped):
                 if name:
                     self._name, _ = _name2__(name, _or_nameof=E)
 
-            def ArcDirect(self, lat1, lon1, azi1, a12, outmask=Caps._STD):
+            def ArcDirect(self, lat1, lon1, azi1, a12, outmask=Caps._STD):  # PYCHOK no cover
                 '''Return the C{_Geodesic.ArcDirect} result as L{GDict}.
                 '''
                 with _wargs(self, lat1, lon1, azi1, a12, outmask) as args:
                     d = _Geodesic.ArcDirect(self, *args)
                 return GDict(d)
 
-            def ArcDirectLine(self, lat1, lon1, azi1, a12, caps=Caps._STD_LINE, **name):
+            def ArcDirectLine(self, lat1, lon1, azi1, a12, caps=Caps._STD_LINE, **name):  # PYCHOK no cover
                 '''Return the C{_Geodesic.ArcDirectLine} as I{wrapped} C{GeodesicLine}.
                 '''
                 return self._GenDirectLine(lat1, lon1, azi1, True, a12, caps, **name)
@@ -445,7 +445,7 @@ _wrapped = _gWrapped()  # PYCHOK singleton, .ellipsoids, .test/base.py
 
 class Geodesic(_gWrapped):  # overwritten by 1st instance
     '''I{Wrapper} around I{Karney}'s class U{geographiclib.geodesic.Geodesic
-       <https://geographiclib.sourceforge.io/Python/doc/code.html>}.
+       <https://GeographicLib.SourceForge.io/Python/doc/code.html>}.
     '''
     def __new__(unused, a_ellipsoid=_EWGS84, f=None, **name):
         '''Return a I{wrapped} C{geodesic.Geodesic} instance from I{Karney}'s
@@ -465,7 +465,7 @@ class Geodesic(_gWrapped):  # overwritten by 1st instance
 
 class GeodesicLine(_gWrapped):  # overwritten by 1st instance
     '''I{Wrapper} around I{Karney}'s class U{geographiclib.geodesicline.GeodesicLine
-       <https://geographiclib.sourceforge.io/Python/doc/code.html>}.
+       <https://GeographicLib.SourceForge.io/Python/doc/code.html>}.
     '''
     def __new__(unused, geodesic, lat1, lon1, azi1, caps=Caps._STD_LINE, **name):
         '''Return a I{wrapped} C{geodesicline.GeodesicLine} instance from I{Karney}'s
@@ -519,13 +519,12 @@ def _Intersecant2(gl, lat0, lon0, radius, tol=_TOL, form=F_D):  # MCCABE in LatL
     n  = _DUNDER_nameof(_Intersecant2)[1:]
     _P =  gl.Position
     _I =  gl.geodesic.Inverse
-    _a =  fabs
 
     def _R3(s):
         # radius, intersection, etc. at distance C{s}
         P = _P(s)
         d = _I(lat0, lon0, P.lat2, P.lon2)
-        return _a(d.s12), P, d
+        return fabs(d.s12), P, d
 
     def _bisect2(s, c, Rc, r, tol, _R3):
         _s = Fsum(c).fsumf_
@@ -536,18 +535,20 @@ def _Intersecant2(gl, lat0, lon0, radius, tol=_TOL, form=F_D):  # MCCABE in LatL
                 break
         else:  # b >>> s and c >>> s
             raise ValueError(Fmt.no_convergence(b, s))
-        __2 = _0_5  # Rb > r > Rc
+        # Rb > r > Rc
         for i in range(_TRIPS):  # 47-48
-            s = (b + c) * __2
+            s = (b + c) * _0_5
             R, P, d = _R3(s)
             if Rb > R > r:
                 b, Rb = s, R
             elif Rc < R < r:
                 c, Rc = s, R
-            t = _a(b - c)
-            if t < tol:  # or _a(R - r) < tol:
+#           else:
+#               break
+            t = fabs(b - c)
+            if t < tol:  # or fabs(R - r) < tol:
                 break
-        else:  # t = min(t, _a(R - r))
+        else:  # t = min(t, fabs(R - r))
             raise ValueError(Fmt.no_convergence(t, tol))
         i += C.iteration  # combine iterations
         P.set_(lat0=lat0, lon0=lon0, azi0=d.azi1, iteration=i,
@@ -558,13 +559,13 @@ def _Intersecant2(gl, lat0, lon0, radius, tol=_TOL, form=F_D):  # MCCABE in LatL
     # one the plumb, pseudo-rhumb line to the other
     C = _PlumbTo(gl, lat0, lon0, tol=tol)
     try:
-        a = _a(C.s02)  # distance between centers
+        a = fabs(C.s02)  # distance between centers
         if a < r:
             c =  C.s12  # distance along pseudo-rhumb line
             h = _copysign(r, c)  # past half chord length
             P, p = _bisect2( h, c, a, r, tol, _R3)
             Q, q = _bisect2(-h, c, a, r, tol, _R3)
-            if _a(p - q) < max(EPS, tol):
+            if fabs(p - q) < max(EPS, tol):
                 Q = P
         elif a > r:
             raise ValueError(_too_(Fmt.distant(a)))
@@ -587,7 +588,7 @@ def _PlumbTo(gl, lat0, lon0, est=None, tol=_TOL):
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- 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/geodesicx/_C4_24.py

@@ -4,7 +4,7 @@
 u'''A Python version of part of I{Karney}'s C++ module U{GeodesicExactC4
 <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1GeodesicExactC4.html>}.
 
-Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2008-2023)
+Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2008-2024)
 and licensed under the MIT/X11 License.  For more information, see the
 U{GeographicLib<https://GeographicLib.SourceForge.io>} documentation.
 '''
@@ -1676,7 +1676,7 @@ del _g, _Gfloats, _f, _f2
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- 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/geodesicx/_C4_27.py

@@ -4,7 +4,7 @@
 u'''A Python version of part of I{Karney}'s C++ module U{GeodesicExactC4
 <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1GeodesicExactC4.html>}.
 
-Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2008-2023)
+Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2008-2024)
 and licensed under the MIT/X11 License.  For more information, see the
 U{GeographicLib<https://GeographicLib.SourceForge.io>} documentation.
 '''
@@ -2372,7 +2372,7 @@ del _g, _Gfloats, _f, _f2
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- 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/geodesicx/_C4_30.py

@@ -4,7 +4,7 @@
 u'''A Python version of part of I{Karney}'s C++ module U{GeodesicExactC4
 <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1GeodesicExactC4.html>}.
 
-Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2008-2023)
+Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2008-2024)
 and licensed under the MIT/X11 License.  For more information, see the
 U{GeographicLib<https://GeographicLib.SourceForge.io>} documentation.
 '''
@@ -3278,7 +3278,7 @@ del _g, _Gfloats, _f, _f2
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- 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/geodesicx/__init__.py

@@ -4,7 +4,7 @@
 u'''A pure Python version of I{Karney}'s C++ classes U{GeodesicExact
 <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1GeodesicExact.html>}
 and U{GeodesicLineExact
-<https://geographiclib.sourceforge.io/C++/doc/classGeographicLib_1_1GeodesicLine.html>}.
+<https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1GeodesicLine.html>}.
 
 For more details, see the C++ U{GeographicLib<https://GeographicLib.SourceForge.io/C++/doc/index.html>}
 documentation, especially the U{Class List<https://GeographicLib.SourceForge.io/C++/doc/annotated.html>}
@@ -23,11 +23,11 @@ from pygeodesy.karney import Caps, GeodesicError
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY
 
 __all__ = _ALL_LAZY.geodesicx + _ALL_DOCS(Caps, GeodesicError)
-__version__ = '24.09.06'
+__version__ = '24.11.02'
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- 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/geodesicx/__main__.py

@@ -65,7 +65,7 @@ _main(C4order=int(argv[1])) if len(argv) == 2 and argv[1].isdigit() else _main()
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- 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/geodesicx/gx.py

@@ -8,7 +8,7 @@ Class L{GeodesicExact} follows the naming, methods and return values
 of class C{Geodesic} from I{Karney}'s Python U{geographiclib
 <https://GitHub.com/geographiclib/geographiclib-python>}.
 
-Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2008-2023)
+Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2008-2024)
 and licensed under the MIT/X11 License.  For more information, see the
 U{GeographicLib<https://GeographicLib.SourceForge.io>} documentation.
 '''
@@ -58,12 +58,13 @@ from pygeodesy.lazily import _ALL_DOCS, _ALL_MODS as _MODS
 from pygeodesy.namedTuples import Destination3Tuple, Distance3Tuple
 from pygeodesy.props import deprecated_Property, Property, Property_RO, property_RO
 # from pygeodesy.streprs import Fmt  # from .fmath
-from pygeodesy.utily import atan2d as _atan2d_reverse, _unrollon, _Wrap, wrap360
+from pygeodesy.utily import atan2, atan2d as _atan2d_reverse, _unrollon, \
+                           _Wrap, wrap360
 
-from math import atan2, copysign, cos, degrees, fabs, radians, sqrt
+from math import copysign, cos, degrees, fabs, radians, sqrt
 
 __all__ = ()
-__version__ = '24.08.13'
+__version__ = '24.11.24'
 
 _MAXIT1 = 20
 _MAXIT2 = 10 + _MAXIT1 + MANT_DIG  # MANT_DIG == C++ digits
@@ -1346,7 +1347,7 @@ __all__ += _ALL_DOCS(GeodesicExact, GeodesicLineExact)
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- 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/geodesicx/gxarea.py

@@ -10,7 +10,7 @@ Class L{GeodesicAreaExact} is intended to work with instances
 of class L{GeodesicExact} and of I{wrapped} class C{Geodesic},
 see module L{pygeodesy.karney}.
 
-Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2008-2023)
+Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2008-2024)
 and licensed under the MIT/X11 License.  For more information, see the
 U{GeographicLib<https://GeographicLib.SourceForge.io>} documentation.
 '''
@@ -518,7 +518,7 @@ __all__ += _ALL_DOCS(GeodesicAreaExact, PolygonArea)
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- 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/geodesicx/gxbases.py

@@ -3,7 +3,7 @@
 
 u'''(INTERNAL) Private L{geodesicx} base class, functions and constants.
 
-Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2008-2023)
+Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2008-2024)
 and licensed under the MIT/X11 License.  For more information, see the
 U{GeographicLib<https://GeographicLib.SourceForge.io>} documentation.
 '''
@@ -158,7 +158,7 @@ def _xnC4(**name_nC4):
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- 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"),