pygeodesy 24.7.7__py2.py3-none-any.whl → 24.7.24__py2.py3-none-any.whl

pygeodesy/karney.py

@@ -4,14 +4,12 @@
 u'''Wrapper around several C{geomath.Math} functions from I{Karney}'s Python package U{geographiclib
 <https://PyPI.org/project/geographiclib>}, provided that package is installed.
 
-The I{wrapped} class methods return a L{GDict} instance offering access to the C{dict} items
-either by C{key} or by C{attribute} name.
-All methods of the I{wrapped} L{Geodesic<geodesicw.Geodesic>} and L{GeodesicLine<geodesicw.GeodesicLine>}
+Methods of the I{wrapped} L{Geodesic<geodesicw.Geodesic>} and L{GeodesicLine<geodesicw.GeodesicLine>}
 classes return a L{GDict} instance offering access to the C{dict} items either by C{key} or by
 C{attribute} name.
 
-With env variable C{PYGEODESY_GEOGRAPHICLIB} left undefined or set to C{"2"}, modules L{geodesicx},
-L{geodesicw} and this module will use U{GeographicLib 2.0<https://GeographicLib.SourceForge.io/C++/doc/>}
+With env variable C{PYGEODESY_GEOGRAPHICLIB} left undefined or set to C{"2"}, modules L{geodesicw},
+L{geodesicx} and this module will use U{GeographicLib 2.0+<https://GeographicLib.SourceForge.io/C++/doc/>}
 and newer transcoding, otherwise C{1.52} or older.
 
 Karney-based functionality
@@ -60,8 +58,11 @@ Karney-based functionality
   - L{GnomonicExact}, L{GnomonicGeodSolve}, L{GnomonicKarney} -- U{Gnomonic
     <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1Gnomonic.html>}
 
+  - L{Intersector} -- U{Intersect
+    <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1Intersect.html>} from I{GeographicLib 2.3+}
+
   - L{JacobiConformal} -- U{JacobiConformal
-    <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1JacobiConformal.html#details>}
+    <https://geographiclib.sourceforge.io/C++/doc/classGeographicLib_1_1experimental_1_1JacobiConformal.html>}
 
   - L{KTransverseMercator} - U{TransverseMercator
     <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1TransverseMercator.html>}
@@ -94,13 +95,14 @@ are I{transcoded} from C++ classes in I{Karney}'s U{GeographicLib<https://Geogra
 
 2. These C{pygeodesy} modules and classes
 
-  - L{ellipsoidalGeodSolve}, L{ellipsoidalKarney}, L{geodsolve}, L{karney}, L{rhumb.solve}
+  - L{ellipsoidalGeodSolve}, L{ellipsoidalKarney}, L{geodesici}, L{geodsolve}, L{karney}, L{rhumb.solve}
   - L{EquidistantKarney}, L{FrechetKarney}, L{GeodesicSolve}, L{GeodesicLineSolve}, L{GnomonicGeodSolve},
-    L{GnomonicKarney}, L{HeightIDWkarney}
+    L{GnomonicKarney}, L{HeightIDWkarney}, L{Intersectool}
 
-are or use I{wrappers} around I{Karney}'s Python U{geographiclib<https://PyPI.org/project/geographiclib>}
-C{geodesic}, C++ utility U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/GeodSolve.1.html>} or
-C++ utility U{RhumbSolve<https://GeographicLib.SourceForge.io/C++/doc/RhumbSolve.1.html>}.
+are or use I{wrappers} around I{Karney}'s Python U{geographiclib<https://PyPI.org/project/geographiclib>} or
+C++ utility U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/GeodSolve.1.html>},
+U{IntersectTool<https://GeographicLib.SourceForge.io/C++/doc/IntersectTool.1.html>} or
+U{RhumbSolve<https://GeographicLib.SourceForge.io/C++/doc/RhumbSolve.1.html>}.
 
 3. All C{pygeodesy} functions and methods to compute I{ellipsoidal} intersections, nearest points and trilaterations
 
@@ -139,12 +141,11 @@ in C{pygeodesy} are based on I{Karney}'s post U{Area of a spherical polygon
 # make sure int/int division yields float quotient, see .basics
 from __future__ import division as _; del _  # PYCHOK semicolon
 
-from pygeodesy.basics import _copysign, int1s, isint, itemsorted, neg, unsigned0, \
-                             _xgeographiclib, _zip,  _version_info
+from pygeodesy.basics import _copysign, isint, neg, unsigned0, _xgeographiclib, \
+                             _zip,  _version_info
 from pygeodesy.constants import NAN, _isfinite as _math_isfinite, _0_0, \
                                _1_16th, _1_0, _2_0, _180_0, _N_180_0, _360_0
-from pygeodesy.errors import GeodesicError, _ValueError, _xkwds, _xkwds_get1, \
-                            _xkwds_kwds, _xkwds_not
+from pygeodesy.errors import GeodesicError, _ValueError, _xkwds, _xkwds_get1
 from pygeodesy.fmath import cbrt, fremainder, norm2
 # from pygeodesy.internals import _version_info  # from .basics
 from pygeodesy.interns import NN, _2_, _a12_, _area_, _azi1_, _azi2_, _azi12_, \
@@ -153,7 +154,7 @@ from pygeodesy.interns import NN, _2_, _a12_, _area_, _azi1_, _azi2_, _azi12_, \
                              _UNDER_, _X_,  _BAR_  # PYCHOK used!
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS, _getenv
 from pygeodesy.named import ADict, _NamedBase, _NamedTuple, notImplemented, _Pass
-from pygeodesy.props import deprecated_method, Property_RO
+from pygeodesy.props import deprecated_method, Property_RO, property_ROnce
 from pygeodesy.units import Bearing as _Azi, Degrees as _Deg, Lat, Lon, \
                             Meter as _M, Meter2 as _M2, Number_
 from pygeodesy.utily import atan2d, sincos2d, tand, _unrollon,  fabs
@@ -161,7 +162,7 @@ from pygeodesy.utily import atan2d, sincos2d, tand, _unrollon,  fabs
 # from math import fabs  # from .utily
 
 __all__ = _ALL_LAZY.karney
-__version__ = '24.07.09'
+__version__ = '24.07.16'
 
 _K_2_0      = _getenv('PYGEODESY_GEOGRAPHICLIB', _2_) == _2_
 _perimeter_ = 'perimeter'
@@ -208,39 +209,39 @@ class Area3Tuple(_NamedTuple):  # in .geodesicx.gxarea
     _Units_ = ( Number_, _M,          _M2)
 
 
-class Caps(object):  # PYCHOK
+class Caps(object):
     '''(INTERNAL) Overriden by C{Caps} below.
     '''
-    EMPTY          =  0        # formerly aka NONE
-    _CAP_1         =  1 << 0   # for goedesici/-w
-    _CAP_1p        =  1 << 1   # for goedesici/-w
-    _CAP_2         =  1 << 2
-    _CAP_3         =  1 << 3   # for goedesici/-w
-#   _CAP_4         =  1 << 4
-#   _CAP_ALL       =  0x1F
-#   _CAP_MASK      = _CAP_ALL
-    LATITUDE       =  1 << 7   # compute latitude C{lat2}
-    LONGITUDE      =  1 << 8   # compute longitude C{lon2}  _CAP_3
-    AZIMUTH        =  1 << 9   # azimuths C{azi1} and C{azi2}
-    DISTANCE       =  1 << 10  # compute distance C{s12}  _CAP_1
-    DISTANCE_IN    =  1 << 11  # allow distance C{s12} in Direct  _CAP_1 | _CAP_1p
-    REDUCEDLENGTH  =  1 << 12  # compute reduced length C{m12}    _CAP_1 | _CAP_2
-    GEODESICSCALE  =  1 << 13  # compute geodesic scales C{M12} and C{M21}  _CAP_1 | _CAP_2
-    AREA           =  1 << 14  # compute area C{S12}  _CAP_4
-
-    STANDARD       =  AZIMUTH | DISTANCE | DISTANCE_IN | LATITUDE | LONGITUDE
-    ALL            =  0x7F80   # without LONG_UNROLL, LINE_OFF, REVERSE2 and _DEBUG_*
-
-    _DIRECT3       =  AZIMUTH  | LATITUDE | LONGITUDE | _CAP_3   # for goedesicw only
-    _INVERSE3      =  AZIMUTH  | DISTANCE | _CAP_1   # for goedesicw only
-    _STD           =  STANDARD | _CAP_3   | _CAP_1   # for goedesicw only
-    _STD_LINE      = _STD      | _CAP_2   | _CAP_1p  # for goedesici/-w
-
-    LINE_CAPS      = _STD_LINE | REDUCEDLENGTH | GEODESICSCALE  # .geodesici only
-
-    LINE_OFF       =  1 << 15  # Line without updates from parent geodesic or rhumb
-    LONG_UNROLL    =  1 << 16  # unroll C{lon2} in .Direct and .Position
-    REVERSE2       =  1 << 17  # reverse C{azi2}
+    EMPTY         =  0        # formerly aka NONE
+    _CAP_1        =  1 << 0   # for goedesici/-w
+    _CAP_1p       =  1 << 1   # for goedesici/-w
+    _CAP_2        =  1 << 2
+    _CAP_3        =  1 << 3   # for goedesici/-w
+    _CAP_4        =  1 << 4   # for goedesicw
+    _CAP_ALL      =  0x1F
+#   _CAP_MASK     = _CAP_ALL
+    LATITUDE      =  1 << 7   # compute latitude C{lat2}
+    LONGITUDE     =  1 << 8  | _CAP_3  # compute longitude C{lon2}
+    AZIMUTH       =  1 << 9   # azimuths C{azi1} and C{azi2}
+    DISTANCE      =  1 << 10 | _CAP_1  # compute distance C{s12}
+    DISTANCE_IN   =  1 << 11 | _CAP_1 | _CAP_1p  # allow distance C{s12} in Direct
+    REDUCEDLENGTH =  1 << 12 | _CAP_1 | _CAP_2  # compute reduced length C{m12}
+    GEODESICSCALE =  1 << 13 | _CAP_1 | _CAP_2  # compute geodesic scales C{M12} and C{M21}
+    AREA          =  1 << 14 | _CAP_4  # compute area C{S12}
+
+    STANDARD      =  AZIMUTH | DISTANCE | LATITUDE | LONGITUDE
+
+    _DIRECT3      =  AZIMUTH  | LATITUDE | LONGITUDE  # for goedesicw only
+    _INVERSE3     =  AZIMUTH  | DISTANCE  # for goedesicw only
+    _STD          =  STANDARD  # for goedesicw only
+    _STD_LINE     = _STD      | DISTANCE_IN  # for goedesici/-w
+
+    LINE_CAPS     = _STD_LINE | REDUCEDLENGTH | GEODESICSCALE  # .geodesici only
+    LONG_UNROLL   =  1 << 15  # unroll C{lon2} in .Direct and .Position
+#                 =  1 << 16  # unused
+    LINE_OFF      =  1 << 17  # Line without updates from parent geodesic or rhumb
+    REVERSE2      =  1 << 18  # reverse C{azi2}
+    ALL           =  0x7F80 | _CAP_ALL  # without LONG_UNROLL, LINE_OFF, REVERSE2 and _DEBUG_*
 
     LATITUDE_LONGITUDE         = LATITUDE | LONGITUDE
     LATITUDE_LONGITUDE_AREA    = LATITUDE | LONGITUDE | AREA
@@ -248,35 +249,38 @@ class Caps(object):  # PYCHOK
     AZIMUTH_DISTANCE           = AZIMUTH | DISTANCE
     AZIMUTH_DISTANCE_AREA      = AZIMUTH | DISTANCE | AREA
 
-    _ANGLE_ONLY    =  1 << 18  # angular distance C{a12} only
-    _S_CALPs_      =  1 << 19  # (INTERNAL) GeodesicExact._GenInverse
+    _SALP_CALPs_   =  1 << 19  # (INTERNAL) GeodesicExact._GenInverse
 
     _DEBUG_AREA    =  1 << 20  # (INTERNAL) include Line details
     _DEBUG_DIRECT  =  1 << 21  # (INTERNAL) include Direct details
     _DEBUG_INVERSE =  1 << 22  # (INTERNAL) include Inverse details
     _DEBUG_LINE    =  1 << 23  # (INTERNAL) include Line details
     _DEBUG_ALL     = _DEBUG_AREA | _DEBUG_DIRECT | _DEBUG_INVERSE | \
-                     _DEBUG_LINE | _ANGLE_ONLY | _S_CALPs_
+                     _DEBUG_LINE | _SALP_CALPs_
 
-    _OUT_ALL       =  ALL
+    _OUT_ALL       =  ALL  # see geographiclib.geodesiccapabilities.py
     _OUT_MASK      =  ALL | LONG_UNROLL | REVERSE2 | _DEBUG_ALL
 
     _AZIMUTH_LATITUDE_LONGITUDE           =  AZIMUTH | LATITUDE | LONGITUDE
+    _AZIMUTH_LATITUDE_LONG_UNROLL         =  AZIMUTH | LATITUDE | LONG_UNROLL
     _DEBUG_DIRECT_LINE                    = _DEBUG_DIRECT | _DEBUG_LINE
 #   _DISTANCE_IN_OUT                      =  DISTANCE_IN & _OUT_MASK  # == DISTANCE_IN in .gx, .gxline
-    _LINE                                 =  AZIMUTH | LATITUDE | LONG_UNROLL
     _REDUCEDLENGTH_GEODESICSCALE          =  REDUCEDLENGTH | GEODESICSCALE
 #   _REDUCEDLENGTH_GEODESICSCALE_DISTANCE =  REDUCEDLENGTH | GEODESICSCALE | DISTANCE
 
-    def toStr(self, Csk, sep=_BAR_):
-        '''Return a C{Caps} or C{outmask} as C{str} or tuple of C{str}s.
+    def items(self):
+        '''Yield all I{public} C{Caps} as 2-tuple C{(NAME, mask)}.
+        '''
+        for n, C in Caps.__class__.__dict__.items():
+            if isint(C) and not n.startswith(_UNDER_) \
+                        and n.replace(_UNDER_, NN).isupper():
+                yield n, C
+
+    def toStr(self, Cask, sep=_BAR_):
+        '''Return C{Caps} or an C{outmask} as C{str} or tuple of C{str}s.
         '''
-        s = []
-        for c, C in itemsorted(self.__class__.__dict__):
-            if isint(C) and (Csk & C) and int1s(C) == 1 \
-                        and (C in (Caps.REVERSE2, Caps._S_CALPs_)
-                               or c.replace(_UNDER_, NN).isupper()):
-                s.append(c)
+        s = (n for n, C in sorted(Caps.items())
+                        if C and (Cask & C) == C)  # and int1s(C) <= 3
         return sep.join(s) if sep else tuple(s)
 
 Caps = Caps()  # PYCHOK singleton
@@ -287,7 +291,7 @@ C{AREA} - compute area C{S12},
 
 C{AZIMUTH} - include azimuths C{azi1} and C{azi2},
 
-C{DISTANCE} - compute distance C{s12},
+C{DISTANCE} - compute distance C{s12} and angular distance C{a12},
 
 C{DISTANCE_IN} - allow distance C{s12} in C{.Direct},
 
@@ -311,7 +315,8 @@ and C{ALL} - all of the above.
 
 C{STANDARD} = C{AZIMUTH | DISTANCE | DISTANCE_IN | LATITUDE | LONGITUDE}'''
 
-_KEY2Caps = dict(azi2=Caps.AZIMUTH,  # see GDict._unCaps
+_key2Caps = dict(a12 =Caps.DISTANCE,  # in GDict._unCaps
+                 azi2=Caps.AZIMUTH,
                  lat2=Caps.LATITUDE,
                  lon2=Caps.LONGITUDE,
                  m12 =Caps.REDUCEDLENGTH,
@@ -322,7 +327,7 @@ _KEY2Caps = dict(azi2=Caps.AZIMUTH,  # see GDict._unCaps
 
 
 class _CapsBase(_NamedBase):  # in .auxilats, .geodesicx.gxbases
-    '''(INTERNAL) Base class for C{[_]Geodesic*Exact}.
+    '''(INTERNAL) Base class for C{Geodesic*}, C{Geodesic*Exact}, C{Intersectool} and C{Rhumb*Base}.
     '''
     ALL           = Caps.ALL
     AREA          = Caps.AREA
@@ -473,8 +478,8 @@ class GDict(ADict):  # XXX _NamedDict
         t  = tuple(_g(self, n, dflt) for n in nTuple._Names_)
         return nTuple(t, iteration=self._iteration)
 
-    def _2X(self, gl, _2X=_X_):  # .Intesectool, .Intersector
-        '''(INTERNAL) Rename attr tail from C{-2} to C{-X} or C{-M}.
+    def _2X(self, gl, _2X=_X_):  # .Intersectool, .Intersector
+        '''(INTERNAL) Rename C{-2} attr to C{-X} or C{-M}.
         '''
         X = GDict(self)
         for n in (_lat2_, _lon2_, _azi2_, _s12_, _a12_):
@@ -488,8 +493,8 @@ class GDict(ADict):  # XXX _NamedDict
     def _unCaps(self, outmask):  # in .geodsolve
         '''(INTERNAL) Remove superfluous items.
         '''
-        for k, m in _KEY2Caps.items():
-            if k in self and not (outmask & m):
+        for k, C in _key2Caps.items():
+            if k in self and (outmask & C) != C:
                 self.pop(k)  # delattr(self, k)
         return self
 
@@ -531,13 +536,13 @@ class _kWrapped(object):  # in .geodesicw
         <https://PyPI.org/project/geographiclib>} classes.
     '''
 
-    @Property_RO
+    @property_ROnce
     def geographiclib(self):
         '''Lazily import C{geographiclib}, provided the U{geographiclib
            <https://PyPI.org/project/geographiclib>} package is installed,
            otherwise raise a C{LazyImportError}.
         '''
-        g = _xgeographiclib(self.__class__, 1, 49)
+        g = _xgeographiclib(self.__class__.__module__, 1, 49)
         from geographiclib.geodesic import Geodesic
         g.Geodesic = Geodesic
         from geographiclib.geodesicline import GeodesicLine
@@ -546,7 +551,7 @@ class _kWrapped(object):  # in .geodesicw
         g.Math = Math
         return g
 
-    @Property_RO
+    @property_ROnce
     def Math(self):
         '''Wrap the C{geomath.Math} class, provided the U{geographiclib
            <https://PyPI.org/project/geographiclib>} package is installed,
@@ -599,7 +604,7 @@ class Rhumb8Tuple(_GTuple):
 
            @kwarg dflt: Default value for missing items (C{any}).
            @kwarg a12_m12_M12_M21_salp1_calp1_salp2_calp2: Optional keyword
-                     arguments to specify or override L{Inverse10Tuple} items.
+                      arguments to specify or override L{Inverse10Tuple} items.
 
            @return: L{Inverse10Tuple}C{(a12, s12, salp1, calp1, salp2, calp2,
                     m12, M12, M21, S12)}.
@@ -724,14 +729,15 @@ def _isfinite(x):  # mimick geomath.Math.isfinite
         return _math_isfinite(x)  # and fabs(x) <= _MAX
 
 
-def _llz2line(line, **llz2):
-    '''(INTERNAL) Set C{line.lat2, .lon2, .azi2} from C{llz2}.
+def _llz2gl(gl, **llz2):  # see .geodesici._llz2G
+    '''(INTERNAL) Set C{gl.lat2, .lon2, .azi2} from C{llz2}.
     '''
     if llz2:
-        llz2 = _xkwds_not(None, **_xkwds_kwds(llz2, lat2=None, lon2=None, azi2=None))
-        if llz2:
-            line.__dict__.update(llz2)
-    return line
+        for n in (_lat2_, _lon2_, _azi2_):  # _lat1_, _lon1_, _azi1_
+            v = llz2.get(n, None)
+            if v is not None:
+                setattr(gl, n, v)
+    return gl
 
 
 def _norm2(x, y):  # mimick geomath.Math.norm

pygeodesy/ktm.py

@@ -67,7 +67,7 @@ from cmath import polar
 from math import atan2, asinh, cos, cosh, degrees, fabs, sin, sinh, sqrt, tanh
 
 __all__ = _ALL_LAZY.ktm
-__version__ = '24.06.11'
+__version__ = '24.07.16'
 
 
 class KTMError(_ValueError):
@@ -433,10 +433,7 @@ def _cma(a, b0, b1, Cn):
        @see: CPython function U{_Py_c_prod<https://GitHub.com/python/
              cpython/blob/main/Objects/complexobject.c>}.
 
-       @note: Python function C{cmath.fsum} is no longer available,
-              but stil mentioned in Note 4 of the comments before
-              CPython function U{math_fsum<https://GitHub.com/python/
-              cpython/blob/main/Modules/mathmodule.c>}
+       @note: Python function C{cmath.fsum} is no longer available.
     '''
     r = fsum1f_(a.real * b0.real, -a.imag * b0.imag, -b1.real, Cn)
     j = fsum1f_(a.real * b0.imag,  a.imag * b0.real, -b1.imag)

pygeodesy/latlonBase.py

@@ -41,7 +41,7 @@ from pygeodesy.namedTuples import Bounds2Tuple, LatLon2Tuple, PhiLam2Tuple, \
                                   Trilaterate5Tuple
 # from pygeodesy.nvectorBase import _N_vector_  # _MODS
 from pygeodesy.props import deprecated_method, Property, Property_RO, \
-                            property_RO, _update_all
+                            property_RO, property_ROnce, _update_all
 # from pygeodesy.streprs import Fmt, hstr  # from .named, _MODS
 from pygeodesy.units import _isDegrees, _isRadius, Distance_, Lat, Lon, \
                              Height, Radius, Radius_, Scalar, Scalar_
@@ -54,7 +54,7 @@ from contextlib import contextmanager
 from math import asin, cos, degrees, fabs, radians
 
 __all__ = _ALL_LAZY.latlonBase
-__version__ = '24.06.11'
+__version__ = '24.07.12'
 
 _formy = _MODS.into(formy=__name__)
 
@@ -475,12 +475,11 @@ class LatLonBase(_NamedBase):
         r = func_(phi2, self.phi, lam21, datum=D)
         return r * (D.ellipsoid.a if radius is None else radius)
 
-    @property_RO
+    @property_ROnce
     def Ecef(self):
         '''Get the ECEF I{class} (L{EcefKarney}), I{once}.
         '''
-        LatLonBase.Ecef = E = _MODS.ecef.EcefKarney  # overwrite property_RO
-        return E
+        return _MODS.ecef.EcefKarney
 
     @Property_RO
     def _Ecef_forward(self):

pygeodesy/lazily.py

@@ -284,8 +284,8 @@ _ALL_LAZY = _NamedEnum_RO(_name='_ALL_LAZY',
                           fsums=_i('Fsum', 'DivMod2Tuple', 'Fsum2Tuple', 'ResidualError',
                                    'fsum', 'fsum_', 'fsumf_', 'fsum1', 'fsum1_', 'fsum1f_'),
                            gars=_i('Garef', 'GARSError'),
-                      geodesici=_i('Intersectool', 'Intersectool5Tuple', 'Middle5Tuple',
-                                   'Intersector',  'Intersector5Tuple',  'XDict'),
+                      geodesici=_i('Intersectool', 'Intersectool5Tuple', 'Intersect7Tuple',
+                                   'Intersector',  'Intersector5Tuple',  'Middle5Tuple', 'XDict'),
                       geodesicw=_i('Geodesic', 'GeodesicLine', 'Geodesic_WGS84'),
                       geodesicx=_i('gx', 'gxarea', 'gxbases', 'gxline',  # modules
                                    'GeodesicAreaExact', 'GeodesicExact', 'GeodesicLineExact', 'PolygonArea'),
@@ -343,7 +343,8 @@ _ALL_LAZY = _NamedEnum_RO(_name='_ALL_LAZY',
                                    'areaOf', 'boundsOf', 'centroidOf', 'fractional',
                                    'isclockwise', 'isconvex', 'isconvex_', 'isenclosedBy', 'ispolar',
                                    'luneOf', 'nearestOn5', 'perimeterOf', 'quadOf'),
-                          props=_i('Property', 'Property_RO', 'property_RO', 'property_doc_',
+                          props=_i('Property', 'Property_RO', 'property_doc_',
+                                   'property_RO', 'property_ROnce', 'property_ROver',
                                    'deprecated_class', 'deprecated_function', 'deprecated_method',
                                    'deprecated_Property_RO', 'deprecated_property_RO', 'DeprecationWarnings'),
                      resections=_i('Collins5Tuple', 'ResectionError', 'Survey3Tuple', 'Tienstra7Tuple',
@@ -521,7 +522,7 @@ class _ALL_MODS(_internals._MODS_Base):
 _internals._MODS = _ALL_MODS = _ALL_MODS()  # PYCHOK singleton
 
 __all__ = _ALL_LAZY.lazily
-__version__ = '24.07.09'
+__version__ = '24.07.18'
 
 
 def _ALL_OTHER(*objs):

pygeodesy/ltp.py

@@ -33,8 +33,8 @@ from pygeodesy.ltpTuples import Attitude4Tuple, ChLVEN2Tuple, ChLV9Tuple, \
                                 ChLVyx2Tuple, _XyzLocals4, _XyzLocals5, Xyz4Tuple
 from pygeodesy.named import _name__, _name2__, _NamedBase, notOverloaded
 from pygeodesy.namedTuples import LatLon3Tuple, LatLon4Tuple, Vector3Tuple
-from pygeodesy.props import Property, Property_RO, property_doc_, property_RO, \
-                           _update_all
+from pygeodesy.props import Property, Property_RO, property_doc_, \
+                            property_ROver, _update_all
 from pygeodesy.streprs import Fmt, strs, unstr
 from pygeodesy.units import Bearing, Degrees, _isHeight, Meter
 from pygeodesy.utily import cotd, _loneg, sincos2d, sincos2d_, tand, tand_, \
@@ -44,7 +44,7 @@ from pygeodesy.vector3d import _ALL_LAZY, Vector3d
 # from math import fabs, floor as _floor  # from .fmath, .fsums
 
 __all__ = _ALL_LAZY.ltp
-__version__ = '24.06.11'
+__version__ = '24.07.12'
 
 _height0_ = _height_ + _0_
 _narrow_  = 'narrow'
@@ -722,15 +722,15 @@ class _ChLV(object):
         t = Y_X_h_lat_lon_h + (self, self._t0, None)  # PYCHOK _t0
         return ChLV9Tuple(t, name=name)
 
-    @property_RO
+    @property_ROver
     def _enh_n_h(self):
         '''(INTERNAL) Get C{ChLV*.reverse} args[1:4] names, I{once}.
         '''
-        _ChLV._enh_n_h = t = _args_kwds_names(_ChLV.reverse)[1:4]  # overwrite property_RO
+        t = _args_kwds_names(_ChLV.reverse)[1:4]
         # assert _args_kwds_names( ChLV.reverse)[1:4] == t
         # assert _args_kwds_names(ChLVa.reverse)[1:4] == t
         # assert _args_kwds_names(ChLVe.reverse)[1:4] == t
-        return t
+        return t  # overwrite propertyROver
 
     def forward(self, latlonh, lon=None, height=0, M=None, **name):  # PYCHOK no cover
         '''Convert WGS84 geodetic to I{Swiss} projection coordinates.  I{Must be overloaded}.

pygeodesy/named.py

@@ -34,7 +34,7 @@ from pygeodesy.streprs import attrs, Fmt, lrstrip, pairs, reprs, unstr
 # from pygeodesy.units import _toUnit  # _MODS
 
 __all__ = _ALL_LAZY.named
-__version__ = '24.06.24'
+__version__ = '24.07.12'
 
 _COMMANL_           = _COMMA_ + _NL_
 _COMMASPACEDOT_     = _COMMASPACE_ + _DOT_
@@ -1231,9 +1231,9 @@ def modulename(clas, prefixed=None):  # in .basics._xversion
        @return: The B{C{class}}'s C{[module.]class} name (C{str}).
     '''
     try:
-        n =  clas.__name__
+        n = clas.__name__
     except AttributeError:
-        n = _dunder_name_
+        n = clas if isstr(clas) else _dunder_name_
     if prefixed or (classnaming() if prefixed is None else False):
         try:
             m =  clas.__module__.rsplit(_DOT_, 1)

pygeodesy/nvectorBase.py

@@ -29,7 +29,7 @@ from pygeodesy.named import _xother3,  _under
 from pygeodesy.namedTuples import Trilaterate5Tuple, Vector3Tuple, \
                                   Vector4Tuple,  map1
 from pygeodesy.props import deprecated_method, Property_RO, property_doc_, \
-                            property_RO, _update_all
+                            property_RO, property_ROnce, _update_all
 from pygeodesy.streprs import Fmt, hstr, unstr
 from pygeodesy.units import Bearing, Height, Radius_, Scalar
 from pygeodesy.utily import sincos2d, _unrollon, _unrollon3
@@ -38,7 +38,7 @@ from pygeodesy.vector3d import Vector3d, _xyzhdlln4
 from math import fabs, sqrt
 
 __all__ = _ALL_LAZY.nvectorBase
-__version__ = '24.06.12'
+__version__ = '24.07.12'
 
 
 class NvectorBase(Vector3d):  # XXX kept private
@@ -78,12 +78,11 @@ class NvectorBase(Vector3d):  # XXX kept private
         '''
         return self._datum
 
-    @property_RO
+    @property_ROnce
     def Ecef(self):
         '''Get the ECEF I{class} (L{EcefKarney}), I{once}.
         '''
-        NvectorBase.Ecef = E = _MODS.ecef.EcefKarney  # overwrite property_RO
-        return E
+        return _MODS.ecef.EcefKarney
 
     @property_RO
     def ellipsoidalNvector(self):

pygeodesy/props.py

@@ -25,7 +25,7 @@ from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS, \
 from functools import wraps as _wraps
 
 __all__ = _ALL_LAZY.props
-__version__ = '24.05.26'
+__version__ = '24.07.23'
 
 _class_       = 'class'
 _dont_use_    = _DEPRECATED_ + ", don't use."
@@ -48,11 +48,12 @@ def _allPropertiesOf(Clas_or_inst, *Bases, **excls):
         except AttributeError:
             raise
             S = ()  # not an inst
-    B    = Bases or _PropertyBase
-    _isa = isinstance
+    B =  Bases or _PropertyBase
+    P = _property_RO___
     for C in S:
         for n, p in C.__dict__.items():
-            if _isa(p, B) and p.name == n and n not in excls:
+            if isinstance(p, B) and p.name == n and not (
+               isinstance(p, P) or n in excls):
                 yield p
 
 
@@ -251,9 +252,9 @@ class Property_RO(_PropertyBase):
                  <https://Docs.Python.org/3/library/functools.html#functools.cache>}
                  to I{cache} or I{memoize} the property value.
 
-           @see: Luciano Ramalho, "Fluent Python", O'Reilly, Example 19-24, 2016
-                 p. 636 or Example 22-28, 2022 p. 870 and U{class Property
-                 <https://docs.Python.org/3/howto/descriptor.html>}.
+           @see: Luciano Ramalho, "Fluent Python", O'Reilly, 2016 p. 636
+                 Example 19-24 or 2022 p. 870 Example 22-28 and U{class
+                 Property<https://docs.Python.org/3/howto/descriptor.html>}.
         '''
         _fget = method if _FOR_DOCS else self._fget  # XXX force method.__doc__ to epydoc
         _PropertyBase.__init__(self, method, _fget, self._fset_error)
@@ -355,6 +356,61 @@ class property_RO(_PropertyBase):
                 inst.__dict__.pop(uname)
 
 
+class _property_RO___(_PropertyBase):
+    # No __doc__ on purpose
+
+    def __init__(self, method, doc=NN):  # PYCHOK expected
+        '''New C{property_ROnce} or C{property_ROver}, holding a singleton value as
+           class attribute for all instances of that class and overwriting C{self},
+           the C{property_ROver} instance in the 1st invokation.
+
+           @see: L{property_RO} for further details.
+        '''
+        _PropertyBase.__init__(self, method, self._fget, self._fset_error, doc=doc)
+
+    def _fdel(self, unused):  # PYCHOK no cover
+        '''Silently ignored, always.
+        '''
+        pass
+
+    def _update(self, *unused):  # PYCHOK signature
+        '''(INTERNAL) No-op, ignore updates.
+        '''
+        pass
+
+
+class property_ROnce(_property_RO___):
+    # No __doc__ on purpose
+
+    def _fget(self, inst):
+        '''Get the C{property} value, only I{once} and memoize/cache it.
+        '''
+        try:
+            v = self._val
+        except AttributeError:
+            v = self._val = self.method(inst)
+        return v
+
+
+class property_ROver(_property_RO___):
+    # No __doc__ on purpose
+
+    def _fget(self, inst):
+        '''Get the C{property} value I{once} and overwrite C{self}, this C{property} instance.
+        '''
+        v = self.method(inst)
+        n = self.name
+        C = inst.__class__
+        for c in C.__mro__:  # [:-1]
+            if getattr(c, n, None) is self:
+                setattr(c, n, v)  # overwrite property_ROver
+                break
+        else:
+            n = _DOT_(C.__name__, n)
+            raise _AssertionError(_EQUALSPACED_(n, v))
+        return v
+
+
 class _NamedProperty(property):
     '''Class C{property} with retrievable name.
     '''
@@ -372,16 +428,18 @@ def property_doc_(doc):
 
        @example:
 
-        >>> @property_doc_("documentation text.")
-        >>> def name(self):
-        >>>     ...
+        >>>class Clas(object):
+        >>>
+        >>>     @property_doc_("documentation text.")
+        >>>     def name(self):
+        >>>         ...
         >>>
-        >>> @name.setter
-        >>> def name(self, value):
-        >>>     ...
+        >>>     @name.setter
+        >>>     def name(self, value):
+        >>>         ...
     '''
-    # See Luciano Ramalho, "Fluent Python", O'Reilly, Example 7-23,
-    # 2016 p. 212+, 2022 p. 331+, Example 9-22 and <https://
+    # See Luciano Ramalho, "Fluent Python", O'Reilly, 2016 p. 212+
+    # Example 7-23 or 2022 p. 331+ Example 9-22 and <https://
     # Python-3-Patterns-Idioms-Test.ReadTheDocs.io/en/latest/PythonDecorators.html>
 
     def _documented_property(method):
@@ -396,8 +454,8 @@ def property_doc_(doc):
 def _deprecated(call, kind, qual_d):
     '''(INTERNAL) Decorator for DEPRECATED functions, methods, etc.
 
-       @see: Brett Slatkin, "Effective Python", page 105, 2nd ed,
-             Addison-Wesley, 2019.
+       @see: Brett Slatkin, "Effective Python", 2019 page 105, 2nd
+             ed, Addison-Wesley.
     '''
     doc = _docof(call)
 

pygeodesy/rhumb/solve.py

@@ -21,7 +21,7 @@ from pygeodesy.solveBase import _SolveGDictBase, _SolveGDictLineBase
 from pygeodesy.utily import _unrollon, _Wrap, wrap360
 
 __all__ = _ALL_LAZY.rhumb_solve
-__version__ = '24.06.28'
+__version__ = '24.07.11'
 
 
 class _RhumbSolveBase(_SolveGDictBase):
@@ -173,7 +173,7 @@ class RhumbSolve(_RhumbSolveBase):
         return r
 
     def _GDictInverse(self, lat1, lon1, lat2, lon2, *unused, **floats):  # PYCHOK signature
-        '''(INTERNAL) Get C{_GenInverse}-like result as an 8-item C{GDict}, but I{without} C{_S_CALPs_}.
+        '''(INTERNAL) Get C{_GenInverse}-like result as an 8-item C{GDict}, but I{without} C{_SALP_CALPs_}.
         '''
         i = _RhumbSolveBase._GDictInverse(self, lat1, lon1, lat2, lon2, **floats)
         a = _over(float(i.s12), self._mpd)  # for .Inverse1