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

pygeodesy/geodesicx/gxarea.py

@@ -19,19 +19,20 @@ from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import _copysign, isodd, unsigned0
 from pygeodesy.constants import NAN, _0_0, _0_5, _720_0
+from pygeodesy.fmath import _fma
 from pygeodesy.internals import printf, typename
 # from pygeodesy.interns import _COMMASPACE_  # from .lazily
 from pygeodesy.karney import Area3Tuple, _diff182, GeodesicError, \
-                            _norm180, _remainder, _sum3
+                            _norm180, _remainder, _sum2
 from pygeodesy.lazily import _ALL_DOCS,  _COMMASPACE_
 from pygeodesy.named import ADict, callername, _NamedBase,  pairs
 from pygeodesy.props import Property, Property_RO, property_RO
 # from pygeodesy.streprs import pairs  # from .named
 
-from math import fmod as _fmod
+from math import fabs, fmod as _fmod
 
 __all__ = ()
-__version__ = '25.09.16'
+__version__ = '25.12.23'
 
 
 class GeodesicAreaExact(_NamedBase):
@@ -427,30 +428,39 @@ class _Accumulator(_NamedBase):
            @kwarg y: Initial value (C{scalar}).
            @kwarg name: Optional C{B{name}=NN} (C{str}).
         '''
-        if isinstance(y, _Accumulator):
-            self._s, self._t, self._n = y._s, y._t, 1
-        elif y:
-            self._s, self._n = float(y), 1
+        self._s_t(*_s_t2(y))
+        self._n = 1
         if name:
             self.name = name
 
-    def Add(self, y):
-        '''Add a value.
+    def Add(self, *ys):
+        '''Add one or more scalars or L{_Accumulator}s.
 
            @return: Current C{sum}.
-        '''
-        self._n += 1
-        self._s, self._t, _ = _sum3(self._s, self._t, y)
-        return self._s  # current .Sum()
+
+           @see: C++ U{Accumulator.Add<https://GeographicLib.sourceforge.io/C++/doc/Accumulator_8hpp_source.html>}
+                 for more details about Karney's and Shewchuk's addition.
+        '''
+        # _Accumulator().Add(1, 1e20, 2, 100, 5000, -1e20) ... 5103.0
+        s, t = self._s, self._t
+        for y in ys:
+            for y in _s_t2(y):
+                if y:
+                    t, u = _sum2(t, y)
+                    s, t = _sum2(s, t)
+                    if s:  # accumlate u in t
+                        t += u
+                    else:  # s == 0 implies t == 0
+                        s  = u
+                    self._n += 1
+        return self._s_t(s, t)
 
     def Negate(self):
         '''Negate sum.
 
            @return: Current C{sum}.
         '''
-        self._s = s = -self._s
-        self._t =     -self._t
-        return s  # current .Sum()
+        return self._s_t(-self._s, -self._t)
 
     @property_RO
     def num(self):
@@ -459,22 +469,27 @@ class _Accumulator(_NamedBase):
         return self._n
 
     def Remainder(self, y):
-        '''Remainder on division by B{C{y}}.
+        '''Remainder of division by B{C{y}}.
 
            @return: Remainder of C{sum} / B{C{y}}.
         '''
-        self._s = _remainder(self._s, y)
-#       self._t = _remainder(self._t, y)
-        self._n = -1
-        return self.Add(_0_0)
+        return self._s_t(_remainder(self._s, y),
+                         _remainder(self._t, y))
 
     def Reset(self, y=0):
-        '''Set value from argument.
+        '''Reset from scalar or L{_Accumulator}.
         '''
-        self._n, self._s, self._t = 0, float(y), _0_0
+        self._s_t(*_s_t2(y))
+        self._n = 0
 
     Set = Reset
 
+    def _s_t(self, s, t=0):
+        if t and fabs(s) < fabs(t):
+            s, t = t, s
+        self._s, self._t = s, t
+        return s
+
     def Sum(self, y=0):
         '''Return C{sum + B{y}}.
 
@@ -488,6 +503,17 @@ class _Accumulator(_NamedBase):
             s = self
         return s._s
 
+    def Times(self, y):
+        '''Multiply by a scalar.
+
+           @return: Current C{sum}.
+        '''
+        s  = d = self._s
+        s *= y
+        d = _fma(y, d, -s)
+        t = _fma(y, self._t, d)
+        return self._s_t(s, t)  # current .Sum()
+
     def toStr(self, prec=6, sep=_COMMASPACE_, **unused):  # PYCHOK signature
         '''Return this C{_Accumulator} as string.
 
@@ -502,11 +528,15 @@ class _Accumulator(_NamedBase):
         return sep.join(pairs(d, prec=prec))
 
 
+def _s_t2(y):
+    return (y._s, y._t) if isinstance(y, _Accumulator) else (float(y),)  # PYCHOK OK
+
+
 __all__ += _ALL_DOCS(GeodesicAreaExact, PolygonArea)
 
 # **) 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/geodesicx/gxbases.py

@@ -165,7 +165,7 @@ def _xnC4(**name_nC4):
 
 # **) 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/geodesicx/gxline.py

@@ -667,7 +667,7 @@ __all__ += _ALL_DOCS(_GeodesicLineExact)
 
 # **) 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/geodsolve.py

@@ -3,33 +3,53 @@
 
 u'''Wrapper to invoke I{Karney}'s U{GeodSolve
 <https://GeographicLib.SourceForge.io/C++/doc/GeodSolve.1.html>} utility
-as an (exact) geodesic, but intended I{for testing purposes only}.
+as an (exact) geodesic, but intended I{mainly for testing purposes}.
 
-Set env variable C{PYGEODESY_GEODSOLVE} to the (fully qualified) path
-of the C{GeodSolve} executable.
+Set env variable C{PYGEODESY_GEODSOLVE} to the (fully qualified) path of
+the C{GeodSolve} executable or use property L{GeodesicSolve.GeodSolve
+<geodsolve._GeodesicSolveBase.GeodSolve>}.
 '''
 
 from pygeodesy.basics import _xinstanceof  # typename
-# from pygeodesy.constants import NAN, _0_0  # from .karney
+from pygeodesy.constants import NAN, _0_0
+# from pygeodesy.errors import GeodesicError  # from .karney
 # from pygeodesy.geodesicx import GeodesicAreaExact  # _MODS
-from pygeodesy.interns import _DMAIN_, NN, _UNDER_
-from pygeodesy.karney import Caps, GeodesicError, GeodSolve12Tuple, \
-                            _sincos2d, _Xables,  _0_0, NAN
+from pygeodesy.interns import _a12_, _azi1_, _azi2_, _DMAIN_, _lat1_, \
+                              _lat2_, _lon1_, _lon2_, _m12_, _M12_, \
+                              _M21_, NN, _s12_, _S12_, _UNDER_
+from pygeodesy.karney import Caps, _GTuple, _sincos2d, _Xables,  _Azi, \
+                            _Deg, GeodesicError, _Lat, _Lon, _M, _M2
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS
-from pygeodesy.namedTuples import Destination3Tuple, Distance3Tuple
+from pygeodesy.namedTuples import Destination3Tuple, Distance3Tuple, \
+                          _Pass
 from pygeodesy.props import Property, Property_RO, property_RO
 from pygeodesy.solveBase import _SolveGDictBase, _SolveGDictLineBase
+# from pygeodesy.units import Azimuth as _Azi, Degrees as _Deg, \
+#                             Meter as _M, Meter2 as _M2  # from .karney
 from pygeodesy.utily import _unrollon, _Wrap, wrap360
 
 __all__ = _ALL_LAZY.geodsolve
-__version__ = '25.09.03'
+__version__ = '25.12.31'
+
+
+class GeodSolve12Tuple(_GTuple):
+    '''12-Tuple C{(lat1, lon1, azi1, lat2, lon2, azi2, s12, a12, m12, M12, M21, S12)} with
+       angles C{lat1}, C{lon1}, C{azi1}, C{lat2}, C{lon2} and C{azi2} and arc C{a12} all in
+       C{degrees}, initial C{azi1} and final C{azi2} forward azimuths in bearings, degrees
+       from North, distance C{s12} and reduced length C{m12} in C{meter}, area C{S12} in
+       C{meter} I{squared} and geodesic scale factors C{M12} and C{M21}, both C{scalar}, see
+       U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/GeodSolve.1.html>}.
+    '''
+    # from GeodSolve --help option -f ... lat1 lon1 azi1 lat2 lon2 azi2 s12 a12 m12 M12 M21 S12
+    _Names_ = (_lat1_, _lon1_, _azi1_, _lat2_, _lon2_, _azi2_, _s12_, _a12_, _m12_, _M12_, _M21_, _S12_)
+    _Units_ = (_Lat,   _Lon,   _Azi,   _Lat,   _Lon,   _Azi,   _M,    _Deg,  _Pass, _Pass, _Pass, _M2)
 
 
 class _GeodesicSolveBase(_SolveGDictBase):
     '''(INTERNAL) Base class for L{GeodesicSolve} and L{GeodesicLineSolve}.
     '''
     _Error         =  GeodesicError
-    _Names_Direct  = \
+    _Names_Direct  = _Names_Distance = \
     _Names_Inverse =  GeodSolve12Tuple._Names_
     _Xable_name    = _Xables.GeodSolve.__name__  # typename
     _Xable_path    = _Xables.GeodSolve()
@@ -65,19 +85,6 @@ class _GeodesicSolveBase(_SolveGDictBase):
         '''
         self._setXable(path)
 
-    def toStr(self, **prec_sep):  # PYCHOK signature
-        '''Return this C{GeodesicSolve} as string.
-
-           @kwarg prec_sep: Keyword argumens C{B{prec}=6} and C{B{sep}=", "}
-                       for the C{float} C{prec}ision, number of decimal digits
-                       (0..9) and the C{sep}arator string to join.  Trailing
-                       zero decimals are stripped for B{C{prec}} values of 1
-                       and above, but kept for negative B{C{prec}} values.
-
-           @return: GeodesicSolve items (C{str}).
-        '''
-        return _SolveGDictBase._toStr(self, GeodSolve=self.GeodSolve, **prec_sep)
-
     @Property_RO
     def _u_option(self):
         return ('-u',) if self.unroll else ()
@@ -203,6 +210,15 @@ class GeodesicSolve(_GeodesicSolveBase):
         gl._s13 = r.s12  # gl.SetDistance(r.s12)
         return gl
 
+    def toStr(self, **prec_sep_other):  # PYCHOK signature
+        '''Return this C{GeodesicSolve} as string.
+
+           @kwarg prec_sep: See L{toStr<pygeodesy.solveBase._SolveBase.toStr>}.
+
+           @return: GeodesicSolve items (C{str}).
+        '''
+        return _SolveGDictBase.toStr(self, GeodSolve=self.GeodSolve, **prec_sep_other)
+
 
 class GeodesicLineSolve(_GeodesicSolveBase, _SolveGDictLineBase):
     '''Wrapper to invoke I{Karney}'s U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/GeodSolve.1.html>}
@@ -271,7 +287,7 @@ class GeodesicLineSolve(_GeodesicSolveBase, _SolveGDictLineBase):
            @return: A C{GDict} with 12 items C{lat1, lon1, azi1, lat2, lon2,
                     azi2, m12, a12, s12, M12, M21, S12}.
         '''
-        return self._GDictInvoke(self._cmdArc, self._Names_Direct, a12)._unCaps(outmask)
+        return self._GDictInvoke(self._cmdArc, self._Names_Distance, a12)._unCaps(outmask)
 
     @Property_RO
     def azi1(self):
@@ -315,14 +331,14 @@ class GeodesicLineSolve(_GeodesicSolveBase, _SolveGDictLineBase):
         self._notImplemented(lat0, lon0, **kwds)
 
     def Position(self, s12, outmask=Caps.STANDARD):
-        '''Find the position on the line given B{C{s12}}.
+        '''Find the position on the line given at distance B{C{s12}}.
 
            @arg s12: Distance from the first point to the second (C{meter}).
 
            @return: A C{GDict} with 12 items C{lat1, lon1, azi1, lat2, lon2,
                     azi2, m12, a12, s12, M12, M21, S12}, possibly C{a12=NAN}.
         '''
-        return self._GDictInvoke(self._cmdDistance, self._Names_Direct, s12)._unCaps(outmask)
+        return self._GDictInvoke(self._cmdDistance, self._Names_Distance, s12)._unCaps(outmask)
 
     @Property_RO
     def s13(self):
@@ -361,19 +377,14 @@ class GeodesicLineSolve(_GeodesicSolveBase, _SolveGDictLineBase):
 #           _update_all(self)
         return self._a13  # NAN for GeodesicLineExact without Cap.DISTANCE_IN
 
-    def toStr(self, **prec_sep):  # PYCHOK signature
+    def toStr(self, **prec_sep_other):  # PYCHOK signature
         '''Return this C{GeodesicLineSolve} as string.
 
-           @kwarg prec_sep: Keyword argumens C{B{prec}=6} and C{B{sep}=", "}
-                       for the C{float} C{prec}ision, number of decimal digits
-                       (0..9) and the C{sep}arator string to join.  Trailing
-                       zero decimals are stripped for B{C{prec}} values of 1
-                       and above, but kept for negative B{C{prec}} values.
+           @kwarg prec_sep: See L{toStr<pygeodesy.solveBase._SolveBase.toStr>}.
 
            @return: GeodesicLineSolve items (C{str}).
         '''
-        return _SolveGDictLineBase._toStr(self, azi1=self.azi1, geodesic=self._solve,
-                                                GeodSolve=self.GeodSolve, **prec_sep)
+        return _SolveGDictLineBase.toStr(self, geodesic=self._solve, **prec_sep_other)
 
 
 __all__ += _ALL_DOCS(_GeodesicSolveBase)
@@ -385,120 +396,78 @@ if __name__ == _DMAIN_:
         from sys import argv
 
         gS = GeodesicSolve(name='Test')
-        gS.verbose = '--verbose' in argv  # or '-v' in argv
+        gS.verbose = v = '--verbose' in argv  # or '-v' in argv
 
         if not _Xables.X_OK(gS.GeodSolve):  # not set
             gS.GeodSolve = _Xables.GeodSolve(_Xables.bin_)
-        printf('version: %s', gS.version)
+        printf('version: %s', gS.version, nt=v)
 
         r = gS.Direct(40.6, -73.8, 51, 5.5e6)
-        printf('Direct: %r', r, nl=1)
-        printf('Direct3: %r', gS.Direct3(40.6, -73.8, 51, 5.5e6))
+        printf('Direct: %r', r)
+        printf('Direct3: %r', gS.Direct3(40.6, -73.8, 51, 5.5e6), nt=v)
 
-        printf('Inverse: %r',  gS.Inverse( 40.6, -73.8, 51.6, -0.5), nl=1)
+        printf('Inverse: %r',  gS.Inverse( 40.6, -73.8, 51.6, -0.5))
         printf('Inverse1: %r', gS.Inverse1(40.6, -73.8, 51.6, -0.5))
-        printf('Inverse3: %r', gS.Inverse3(40.6, -73.8, 51.6, -0.5))
+        printf('Inverse3: %r', gS.Inverse3(40.6, -73.8, 51.6, -0.5), nt=v)
 
         glS = GeodesicLineSolve(gS, 40.6, -73.8, 51, name='LineTest')
+        printf('Line: %r', glS)
         p = glS.Position(5.5e6)
-        printf('Position:    %5s %r', p == r, p, nl=1)
+        printf('Position:    %r %s', p, p == r)
         p = glS.ArcPosition(49.475527)
-        printf('ArcPosition: %5s %r', p == r, p)
+        printf('ArcPosition: %r %s', p, p == r)
 
     _main()
 
 # % python3 -m pygeodesy.geodsolve
 
-# version: /opt/local/bin/GeodSolve: GeographicLib version 2.2
-
-# Direct: GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141173, m12=4844148.703101, M12=0.650911, M21=0.651229, s12=5500000.0, S12=39735075134877.09375)
-# Direct3: Destination3Tuple(lat=51.884565, lon=-1.141173, final=107.189397)
-
-# Inverse: GDict(a12=49.94131, azi1=51.198883, azi2=107.821777, lat1=40.6, lat2=51.6, lon1=-73.8, lon2=-0.5, m12=4877684.602706, M12=0.64473, M21=0.645046, s12=5551759.400319, S12=40041368848742.53125)
-# Inverse1: 49.94131021789904
-# Inverse3: Distance3Tuple(distance=5551759.400319, initial=51.198883, final=107.821777)
-
-# Position:    True  GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141173, m12=4844148.703101, M12=0.650911, M21=0.651229, s12=5500000.0, S12=39735075134877.09375)
-# ArcPosition: False GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141174, m12=4844148.669561, M12=0.650911, M21=0.651229, s12=5499999.948497, S12=39735074737272.734375)
-
-
-# % python3 -m pygeodesy.geodsolve
-
-# version: /opt/local/bin/GeodSolve: GeographicLib version 2.3
-
+# version: /opt/local/bin/GeodSolve: GeographicLib version 2.7
 # Direct: GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141173, m12=4844148.703101, M12=0.650911, M21=0.651229, s12=5500000.0, S12=39735075134877.078125)
 # Direct3: Destination3Tuple(lat=51.884565, lon=-1.141173, final=107.189397)
-
 # Inverse: GDict(a12=49.94131, azi1=51.198883, azi2=107.821777, lat1=40.6, lat2=51.6, lon1=-73.8, lon2=-0.5, m12=4877684.602706, M12=0.64473, M21=0.645046, s12=5551759.400319, S12=40041368848742.53125)
 # Inverse1: 49.94131021789904
 # Inverse3: Distance3Tuple(distance=5551759.400319, initial=51.198883, final=107.821777)
-
-# Position:    False GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141173, s12=5500000.0)
-# ArcPosition: False GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141174, s12=5499999.948497)
-
-
-# % python3 -m pygeodesy.geodsolve --verbose
-
-# GeodesicSolve 'Test' 1: /opt/local/bin/GeodSolve --version (invoke)
-# GeodesicSolve 'Test' 1: /opt/local/bin/GeodSolve: GeographicLib version 2.2 (0)
-# version: /opt/local/bin/GeodSolve: GeographicLib version 2.2
-# GeodesicSolve 'Test' 2: /opt/local/bin/GeodSolve -f -E -p 10 \ 40.600000000000001 -73.799999999999997 51.0 5500000.0 (Direct)
-# GeodesicSolve 'Test' 2: lat1=40.600000000000001, lon1=-73.799999999999997, azi1=51.0, lat2=51.884564505606761, lon2=-1.141172861200829, azi2=107.189397162605886, s12=5500000.0, a12=49.475527463251467, m12=4844148.703101486, M12=0.65091056699808603, M21=0.65122865892196558, S12=39735075134877.094 (0)
-
-# Direct: GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141173, m12=4844148.703101, M12=0.650911, M21=0.651229, s12=5500000.0, S12=39735075134877.09375)
-# GeodesicSolve 'Test' 3: /opt/local/bin/GeodSolve -f -E -p 10 \ 40.600000000000001 -73.799999999999997 51.0 5500000.0 (Direct3)
-# GeodesicSolve 'Test' 3: lat1=40.600000000000001, lon1=-73.799999999999997, azi1=51.0, lat2=51.884564505606761, lon2=-1.141172861200829, azi2=107.189397162605886, s12=5500000.0, a12=49.475527463251467, m12=4844148.703101486, M12=0.65091056699808603, M21=0.65122865892196558, S12=39735075134877.094 (0)
-# Direct3: Destination3Tuple(lat=51.884565, lon=-1.141173, final=107.189397)
-# GeodesicSolve 'Test' 4: /opt/local/bin/GeodSolve -f -E -p 10 -i \ 40.600000000000001 -73.799999999999997 51.600000000000001 -0.5 (Inverse)
-# GeodesicSolve 'Test' 4: lat1=40.600000000000001, lon1=-73.799999999999997, azi1=51.198882845579824, lat2=51.600000000000001, lon2=-0.5, azi2=107.821776735514248, s12=5551759.4003186841, a12=49.941310217899037, m12=4877684.6027061976, M12=0.64472969205948238, M21=0.64504567852134398, S12=40041368848742.531 (0)
-
-# Inverse: GDict(a12=49.94131, azi1=51.198883, azi2=107.821777, lat1=40.6, lat2=51.6, lon1=-73.8, lon2=-0.5, m12=4877684.602706, M12=0.64473, M21=0.645046, s12=5551759.400319, S12=40041368848742.53125)
-# GeodesicSolve 'Test' 5: /opt/local/bin/GeodSolve -f -E -p 10 -i \ 40.600000000000001 -73.799999999999997 51.600000000000001 -0.5 (Inverse1)
-# GeodesicSolve 'Test' 5: lat1=40.600000000000001, lon1=-73.799999999999997, azi1=51.198882845579824, lat2=51.600000000000001, lon2=-0.5, azi2=107.821776735514248, s12=5551759.4003186841, a12=49.941310217899037, m12=4877684.6027061976, M12=0.64472969205948238, M21=0.64504567852134398, S12=40041368848742.531 (0)
-# Inverse1: 49.94131021789904
-# GeodesicSolve 'Test' 6: /opt/local/bin/GeodSolve -f -E -p 10 -i \ 40.600000000000001 -73.799999999999997 51.600000000000001 -0.5 (Inverse3)
-# GeodesicSolve 'Test' 6: lat1=40.600000000000001, lon1=-73.799999999999997, azi1=51.198882845579824, lat2=51.600000000000001, lon2=-0.5, azi2=107.821776735514248, s12=5551759.4003186841, a12=49.941310217899037, m12=4877684.6027061976, M12=0.64472969205948238, M21=0.64504567852134398, S12=40041368848742.531 (0)
-# Inverse3: Distance3Tuple(distance=5551759.400319, initial=51.198883, final=107.821777)
-
-# Position:    True  GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141173, m12=4844148.703101, M12=0.650911, M21=0.651229, s12=5500000.0, S12=39735075134877.09375)
-# ArcPosition: False GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141174, m12=4844148.669561, M12=0.650911, M21=0.651229, s12=5499999.948497, S12=39735074737272.734375)
+# Line: GeodesicLineSolve(geodesic=GeodesicSolve(ellipsoid=Ellipsoid(name='WGS84', a=6378137, f=0.00335281, f_=298.25722356, b=6356752.31424518), GeodSolve='/opt/local/bin/GeodSolve', invokation=6, status=0), invokation=1, status=0)
+# Position:    GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141173, s12=5500000.0) False
+# ArcPosition: GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141174, s12=5499999.948497) False
 
 
 # % python3 -m pygeodesy.geodsolve --verbose
 
 # GeodesicSolve 'Test'@1: /opt/local/bin/GeodSolve --version (invoke)
-# GeodesicSolve 'Test'@1: '/opt/local/bin/GeodSolve: GeographicLib version 2.3' (0, stdout/-err)
-# GeodesicSolve 'Test'@1: /opt/local/bin/GeodSolve: GeographicLib version 2.3 (0)
-# version: /opt/local/bin/GeodSolve: GeographicLib version 2.3
-# GeodesicSolve 'Test'@2: /opt/local/bin/GeodSolve -f -E -p 10 \ 40.600000000000001 -73.799999999999997 51.0 5500000.0 (Direct)
+# GeodesicSolve 'Test'@1: '/opt/local/bin/GeodSolve: GeographicLib version 2.7' (0, stdout/-err)
+# GeodesicSolve 'Test'@1: /opt/local/bin/GeodSolve: GeographicLib version 2.7 (0)
+# version: /opt/local/bin/GeodSolve: GeographicLib version 2.7
+#
+# GeodesicSolve 'Test'@2: /opt/local/bin/GeodSolve -f -E -p 10 \ 40.6 -73.8 51.0 5500000.0 (Direct)
 # GeodesicSolve 'Test'@2: '40.600000000000001 -73.799999999999997 51.000000000000000 51.884564505606761 -1.141172861200843 107.189397162605871 5500000.0000000000 49.475527463251460 4844148.7031014860 0.65091056699808614 0.65122865892196569 39735075134877.078' (0, stdout/-err)
 # GeodesicSolve 'Test'@2: lat1=40.600000000000001, lon1=-73.799999999999997, azi1=51.0, lat2=51.884564505606761, lon2=-1.141172861200843, azi2=107.189397162605871, s12=5500000.0, a12=49.47552746325146, m12=4844148.703101486, M12=0.65091056699808614, M21=0.65122865892196569, S12=39735075134877.078 (0)
-
 # Direct: GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141173, m12=4844148.703101, M12=0.650911, M21=0.651229, s12=5500000.0, S12=39735075134877.078125)
-# GeodesicSolve 'Test'@3: /opt/local/bin/GeodSolve -f -E -p 10 \ 40.600000000000001 -73.799999999999997 51.0 5500000.0 (Direct3)
+# GeodesicSolve 'Test'@3: /opt/local/bin/GeodSolve -f -E -p 10 \ 40.6 -73.8 51.0 5500000.0 (Direct3)
 # GeodesicSolve 'Test'@3: '40.600000000000001 -73.799999999999997 51.000000000000000 51.884564505606761 -1.141172861200843 107.189397162605871 5500000.0000000000 49.475527463251460 4844148.7031014860 0.65091056699808614 0.65122865892196569 39735075134877.078' (0, stdout/-err)
 # GeodesicSolve 'Test'@3: lat1=40.600000000000001, lon1=-73.799999999999997, azi1=51.0, lat2=51.884564505606761, lon2=-1.141172861200843, azi2=107.189397162605871, s12=5500000.0, a12=49.47552746325146, m12=4844148.703101486, M12=0.65091056699808614, M21=0.65122865892196569, S12=39735075134877.078 (0)
 # Direct3: Destination3Tuple(lat=51.884565, lon=-1.141173, final=107.189397)
-# GeodesicSolve 'Test'@4: /opt/local/bin/GeodSolve -f -E -p 10 -i \ 40.600000000000001 -73.799999999999997 51.600000000000001 -0.5 (Inverse)
+#
+# GeodesicSolve 'Test'@4: /opt/local/bin/GeodSolve -f -E -p 10 -i \ 40.6 -73.8 51.6 -0.5 (Inverse)
 # GeodesicSolve 'Test'@4: '40.600000000000001 -73.799999999999997 51.198882845579824 51.600000000000001 -0.500000000000000 107.821776735514248 5551759.4003186813 49.941310217899037 4877684.6027061967 0.64472969205948238 0.64504567852134398 40041368848742.531' (0, stdout/-err)
 # GeodesicSolve 'Test'@4: lat1=40.600000000000001, lon1=-73.799999999999997, azi1=51.198882845579824, lat2=51.600000000000001, lon2=-0.5, azi2=107.821776735514248, s12=5551759.4003186813, a12=49.941310217899037, m12=4877684.6027061967, M12=0.64472969205948238, M21=0.64504567852134398, S12=40041368848742.531 (0)
-
 # Inverse: GDict(a12=49.94131, azi1=51.198883, azi2=107.821777, lat1=40.6, lat2=51.6, lon1=-73.8, lon2=-0.5, m12=4877684.602706, M12=0.64473, M21=0.645046, s12=5551759.400319, S12=40041368848742.53125)
-# GeodesicSolve 'Test'@5: /opt/local/bin/GeodSolve -f -E -p 10 -i \ 40.600000000000001 -73.799999999999997 51.600000000000001 -0.5 (Inverse1)
+# GeodesicSolve 'Test'@5: /opt/local/bin/GeodSolve -f -E -p 10 -i \ 40.6 -73.8 51.6 -0.5 (Inverse1)
 # GeodesicSolve 'Test'@5: '40.600000000000001 -73.799999999999997 51.198882845579824 51.600000000000001 -0.500000000000000 107.821776735514248 5551759.4003186813 49.941310217899037 4877684.6027061967 0.64472969205948238 0.64504567852134398 40041368848742.531' (0, stdout/-err)
 # GeodesicSolve 'Test'@5: lat1=40.600000000000001, lon1=-73.799999999999997, azi1=51.198882845579824, lat2=51.600000000000001, lon2=-0.5, azi2=107.821776735514248, s12=5551759.4003186813, a12=49.941310217899037, m12=4877684.6027061967, M12=0.64472969205948238, M21=0.64504567852134398, S12=40041368848742.531 (0)
 # Inverse1: 49.94131021789904
-# GeodesicSolve 'Test'@6: /opt/local/bin/GeodSolve -f -E -p 10 -i \ 40.600000000000001 -73.799999999999997 51.600000000000001 -0.5 (Inverse3)
+# GeodesicSolve 'Test'@6: /opt/local/bin/GeodSolve -f -E -p 10 -i \ 40.6 -73.8 51.6 -0.5 (Inverse3)
 # GeodesicSolve 'Test'@6: '40.600000000000001 -73.799999999999997 51.198882845579824 51.600000000000001 -0.500000000000000 107.821776735514248 5551759.4003186813 49.941310217899037 4877684.6027061967 0.64472969205948238 0.64504567852134398 40041368848742.531' (0, stdout/-err)
 # GeodesicSolve 'Test'@6: lat1=40.600000000000001, lon1=-73.799999999999997, azi1=51.198882845579824, lat2=51.600000000000001, lon2=-0.5, azi2=107.821776735514248, s12=5551759.4003186813, a12=49.941310217899037, m12=4877684.6027061967, M12=0.64472969205948238, M21=0.64504567852134398, S12=40041368848742.531 (0)
 # Inverse3: Distance3Tuple(distance=5551759.400319, initial=51.198883, final=107.821777)
-
-# Position:    False GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141173, s12=5500000.0)
-# ArcPosition: False GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141174, s12=5499999.948497)
+#
+# Line: GeodesicLineSolve(geodesic=GeodesicSolve(ellipsoid=Ellipsoid(name='WGS84', a=6378137, f=0.00335281, f_=298.25722356, b=6356752.31424518), GeodSolve='/opt/local/bin/GeodSolve', invokation=6, status=0), invokation=1, status=0)
+# Position:    GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141173, s12=5500000.0) False
+# ArcPosition: GDict(a12=49.475527, azi1=51.0, azi2=107.189397, lat1=40.6, lat2=51.884565, lon1=-73.8, lon2=-1.141174, s12=5499999.948497) False
 
 # **) 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/geohash.py

@@ -1098,7 +1098,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/geoids.py

@@ -1912,7 +1912,7 @@ del _intCs, _T, _T0s12  # trash ints cache and map
 
 # **) 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/hausdorff.py

@@ -828,7 +828,7 @@ def randomrangenerator(seed):
 
 # **) 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/heights.py

@@ -1073,7 +1073,7 @@ __all__ += _ALL_DOCS(_HeightBase, _HeightIDW, _HeightNamed)
 
 # **) 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/internals.py

@@ -726,7 +726,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/interns.py

@@ -10,7 +10,7 @@ except NameError:  # Python 3+
     _intern = _sys.intern
 
 _COMMASPACE_  = ', '  # overriden below
-_SUB_PACKAGES = 'auxilats', 'deprecated', 'geodesicx', 'rhumb'  # PYCHOK in ...
+_SUB_PACKAGES = 'auxilats', 'deprecated', 'geodesicx', 'rhumb', 'triaxials'  # PYCHOK in ...
 # ... .lazily, make._dist, MANIFEST, setup.setup, test/bases, test/testModules
 
 
@@ -443,7 +443,7 @@ _LR_PAIRS  = {_LANGLE_:  _RANGLE_,
 
 __all__ = (_NN_,  # NOT MISSING!
             Str_.__name__)  # classes
-__version__ = '25.05.26'
+__version__ = '25.11.07'
 
 if __name__ == _DMAIN_:
 
@@ -466,7 +466,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/iters.py

@@ -527,7 +527,7 @@ __all__ += _ALL_DOCS(_BaseIter)
 
 # **) 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"),