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

pygeodesy/geodsolve.py

@@ -3,33 +3,52 @@
 
 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.
 '''
 
 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.06'
+
+
+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 +84,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 +209,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 +286,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 +330,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 +376,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 +395,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

@@ -498,7 +498,7 @@ def printf(fmt, *args, **nl_nt_prec_prefix__end_file_flush_sep__kwds):
             t =  fmt
     except Exception as x:
         _Error, s = _MODS.errors._xError2(x)
-        _unstr    = _MODS.strepr.unstr
+        _unstr    = _MODS.streprs.unstr
         t = _unstr(printf, fmt, *args, **nl_nt_prec_prefix__end_file_flush_sep__kwds)
         raise _Error(s, txt=t, cause=x)
     try:
@@ -709,7 +709,7 @@ def _versions(sep=_SPACE_):
 
 
 __all__ = tuple(map(typename, (machine, print_, printf, typename)))
-__version__ = '25.10.06'
+__version__ = '25.10.26'
 
 if __name__ == _DMAIN_:
 
@@ -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"),