pygeodesy 25.4.8__py2.py3-none-any.whl → 25.4.25__py2.py3-none-any.whl

pygeodesy/__init__.py

@@ -2,7 +2,7 @@
 # -*- coding: utf-8 -*-
 
 u'''A pure Python implementation of geodesy tools for various ellipsoidal and spherical earth
-models using precision trigonometric, vector-based, exact, elliptic, iterative and approximate
+models using precision exact, elliptic, trigonometric, vector-based, iterative and approximate
 methods for geodetic (lat-/longitude), geocentric (U{ECEF<https://WikiPedia.org/wiki/ECEF>}
 cartesian) and certain U{triaxial ellipsoidal<https://GeographicLib.SourceForge.io/1.44/triaxial.html>}
 coordinates.
@@ -17,14 +17,7 @@ and a geocentric B{C{Cartesian}} class with methods and functions to compute dis
 forward and reverse azimuth, initial and final bearing, intermediate and nearest points, intersections of geodesic,
 great circle and rhumb lines, circle intersections and secants, U{3-point resections
 <https://WikiPedia.org/wiki/Position_resection_and_intersection>}, triangulation, trilateration (by intersection,
-by overlap and in 3-D), conversions and unrolling, among other things.  For more information and further details see
-the U{documentation<https://mrJean1.GitHub.io/PyGeodesy>}, the descriptions of U{Latitude/Longitude
-<https://www.Movable-Type.co.UK/scripts/latlong.html>}, U{Vincenty
-<https://www.Movable-Type.co.UK/scripts/latlong-vincenty.html>} and U{Vector-based
-<https://www.Movable-Type.co.UK/scripts/latlong-vectors.html>} geodesy, the original U{JavaScript source
-<https://GitHub.com/ChrisVeness/geodesy>} or U{docs<https://www.Movable-Type.co.UK/scripts/geodesy/docs>}
-and I{Karney}'s Python U{geographiclib<https://PyPI.org/project/geographiclib>} and C++ U{GeographicLib
-<https://GeographicLib.SourceForge.io/C++/doc/index.html>}.
+by overlap and in 3-D), conversions and unrolling, among other things.
 
 Also included are modules for conversions to and from U{Cassini-Soldner
 <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1CassiniSoldner.html>},
@@ -55,11 +48,20 @@ points (or a U{NumPy array<https://docs.SciPy.org/doc/numpy/reference/generated/
 including implementations of the U{Ramer-Douglas-Peucker<https://WikiPedia.org/wiki/
 Ramer-Douglas-Peucker_algorithm>}, the U{Visvalingam-Whyatt<https://hydra.Hull.ac.UK/
 resources/hull:8338>} and the U{Reumann-Witkam<https://psimpl.SourceForge.net/reumann-witkam.html>}
-algorithms and modified versions of the former.  Other classes provide I{boolean} operations between
-(composite) polygons or U{interpolate <https://docs.SciPy.org/doc/scipy/reference/interpolate.html>}
-the L{height<pygeodesy.heights>} of C{LatLon} points and L{Geoid<pygeodesy.geoids>} models, compute
-various U{Fréchet<https://WikiPedia.org/wiki/Frechet_distance>} or U{Hausdorff<https://WikiPedia.org/
-wiki/Hausdorff_distance>} distances.
+algorithms and modified versions of the former.
+
+Plus modules and classes to U{interpolate<https://docs.SciPy.org/doc/scipy/reference/interpolate.html>} the
+L{height<pygeodesy.heights>} of C{LatLon} points and L{Geoid<pygeodesy.geoids>} models, compute various U{Fréchet
+<https://WikiPedia.org/wiki/Frechet_distance>} or U{Hausdorff<https://WikiPedia.org/wiki/Hausdorff_distance>}
+distances or perform I{boolean} operations between (composite) polygons.
+
+For further details see the U{documentation<https://mrJean1.GitHub.io/PyGeodesy>}, the descriptions of
+U{Latitude/Longitude<https://www.Movable-Type.co.UK/scripts/latlong.html>}, U{Vincenty
+<https://www.Movable-Type.co.UK/scripts/latlong-vincenty.html>} and U{Vector-based
+<https://www.Movable-Type.co.UK/scripts/latlong-vectors.html>} geodesy, the original U{JavaScript source
+<https://GitHub.com/ChrisVeness/geodesy>} or U{docs<https://www.Movable-Type.co.UK/scripts/geodesy/docs>}
+and I{Karney}'s Python U{geographiclib<https://PyPI.org/project/geographiclib>} and C++ U{GeographicLib
+<https://GeographicLib.SourceForge.io/C++/doc/index.html>}.
 
 Installation
 ============
@@ -83,15 +85,16 @@ Dependencies
 ============
 
 Installation of I{Karney}'s Python package U{geographiclib<https://PyPI.org/project/geographiclib>}
-is optional, but required to use modules L{ellipsoidalKarney} and L{css}, L{azimuthal} classes
-L{EquidistantKarney} and L{GnomonicKarney} and the L{HeightIDWkarney} interpolator.
+is optional, but required for module L{ellipsoidalKarney}, L{azimuthal} classes L{EquidistantKarney}
+and L{GnomonicKarney} and the L{HeightIDWkarney} interpolator.
 
 Both U{numpy<https://PyPI.org/project/numpy>} and U{scipy<https://PyPI.org/project/scipy>} must be
 installed for most L{Geoid...<pygeodesy.geoids>} and L{Height...<pygeodesy.heights>} interpolators,
 except L{GeoidKarney} and the L{HeightIDW...<pygeodesy.heights>} ones.
 
-Functions and C{LatLon} methods L{circin6}, L{circum3}, L{circum4_}, L{soddy4}, L{trilaterate3d2} do and
-C{trilaterate5} and modules L{auxilats} and L{rhumb} may require U{numpy<https://PyPI.org/project/numpy>}.
+Functions and C{LatLon} methods L{circin6}, L{circum3}, L{circum4_} and L{soddy4} and functions
+L{triaxum5} and L{trilaterate3d2} require U{numpy<https://PyPI.org/project/numpy>} to be installed,
+modules L{auxilats} and L{rhumb} may need U{numpy<https://PyPI.org/project/numpy>}.
 
 Modules L{ellipsoidalGeodSolve} and L{geodsolve} and L{azimuthal} classes L{EquidistantGeodSolve}
 and L{GnomonicGeodSolve} depend on I{Karney}'s C++ utility U{GeodSolve
@@ -122,7 +125,7 @@ C{epydoc --html --no-private --no-source --name=pygeodesy --url=... -v pygeodesy
 Tests
 =====
 
-The tests ran with Python 3.13.2 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0),
+The tests ran with Python 3.13.3 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0),
 Python 3.12.7 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0,
 U{numpy<https://PyPI.org/project/numpy>} 2.1.0, U{scipy<https://PyPI.org/project/scipy>} 1.14.1,
 U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.3,
@@ -140,20 +143,20 @@ U{GeoConvert<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.3,
 U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.3,
 U{IntersectTool<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.3 and
 U{RhumbSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.3), all in 64-bit on
-macOS 15.4 Sequoia.
+macOS 15.4.1 Sequoia.
 
 All tests ran with and without C{lazy import} for Python 3 and with command line option C{-W default} and
 env variable C{PYGEODESY_WARNINGS=on} for all Python versions.  The results of those tests are included in
 the distribution files.
 
 Test coverage has been measured with U{coverage<https://PyPI.org/project/coverage>} 7.6.1 using Python
-3.13.2, 3.12.7, 3.11.5 and 3.10.8.  The complete coverage report in HTML and a PDF summary are included in
+3.13.3, 3.12.7, 3.11.5 and 3.10.8.  The complete coverage report in HTML and a PDF summary are included in
 the distribution files.
 
-Python 3.13.2, 3.12.7, 3.11.5 and 3.10.8 run on Apple M1 Silicon (C{arm64}), I{natively}.  Python 2.7.18 runs
+Python 3.13.3, 3.12.7, 3.11.5 and 3.10.8 run on Apple M4 Si (C{arm64}), I{natively}.  Python 2.7.18 runs
 on Intel (C{x86_64}) or Intel I{emulation} ("C{arm64_x86_64}", see function L{machine<pygeodesy.machine>}).
 
-The tests also ran with Python 3.13.2 (and U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0) on
+The tests also ran with Python 3.13.3 (and U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0) on
 U{Debian 12<https://Cirrus-CI.com/github/mrJean1/PyGeodesy/master>} in 64-bit only, with Python 3.12.5 (and
 U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0) on U{Windows 2019Server
 <https://CI.AppVeyor.com/project/mrJean1/pygeodesy>} in 64-bit only and with Python 2.7.18 (and U{geographiclib
@@ -163,7 +166,7 @@ in 64- and 32-bit.
 A single-File and single-Directory application with C{pygeodesy} has been bundled using U{PyInstaller
 <https://PyPI.org/project/pyinstaller>} 3.4 and 64-bit Python 3.7.3 on macOS 10.13.6 High Sierra.
 
-Previously, the tests were run with Python 3.13.0-1, 3.12.0-6, 3.11.2-4, 3.10.1-7, 3.9.6, 3.9.1, 3.8.7, 3.7.1, 2.7.15,
+Previously, the tests were run with Python 3.13.0-2, 3.12.0-6, 3.11.2-4, 3.10.1-7, 3.9.6, 3.9.1, 3.8.7, 3.7.1, 2.7.15,
 U{PyPy<https://PyPy.org>} 7.3.12 (Python 3.10.12), 7.3.1 (Python 3.6.9) and U{PyPy<https://PyPy.org>} 7.1.1 (Python
 2.7.13) (and U{geographiclib <https://PyPI.org/project/geographiclib>} 1.52, U{numpy<https://PyPI.org/project/numpy>}
 1.16.3, 1.16.4, 1.16.6, 1.19.0, 1.19.4, 1.19.5 or 1.22.4 and U{scipy<https://PyPI.org/project/scipy>} 1.2.1, 1.4.1,
@@ -187,7 +190,7 @@ All Python source code has been statically U{checked<https://GitHub.com/ActiveSt
 Python/546532_PyChecker_postprocessor>} with U{PyChecker<https://PyPI.org/project/pychecker>}, U{PyFlakes
 <https://PyPI.org/project/pyflakes>}, U{PyCodeStyle<https://PyPI.org/project/pycodestyle>} (formerly Pep8) and
 U{McCabe<https://PyPI.org/project/mccabe>} using Python 2.7.18 and with U{Flake8<https://PyPI.org/project/flake8>}
-using Python 3.11.5, both in 64-bit on macOS 15.4 Sequoia.
+using Python 3.11.5, both in 64-bit on macOS 15.4.1 Sequoia.
 
 For a summary of all I{Karney}-based functionality in C{pygeodesy}, see module U{karney
 <https://mrJean1.GitHub.io/PyGeodesy/docs/pygeodesy.karney-module.html>}.
@@ -583,7 +586,7 @@ if _init__all__ and not _lazy_import2:  # import and set __all__
                     raise ImportError('missing %s%s: %s' % (_attribute_, s, t))
                 ns.extend(attrs)
 # XXX       if ps:  # check that mod is a _pygeodesy_ module
-# XXX           m = globalocals[mod]  # assert(_DUNDER_nameof(m) == mod)
+# XXX           m = globalocals[mod]  # assert(typename(m) == mod)
 # XXX           f = getattr(m, _dunder_file_, _NN)
 # XXX           d = _os_path.dirname(_os_path.abspath(f)) if f else pygeodesy_abspath
 # XXX           p = getattr(m, _dunder_package_, _NN) or _pygeodesy_
@@ -599,7 +602,7 @@ else:
 
 from pygeodesy.internals import _version2,  _DOT_  # PYCHOK import
 # from pygeodesy.interns import _DOT_  # from .internals
-__version__ = '25.04.08'
+__version__ = '25.04.25'
 # see setup.py for similar logic
 version     = _DOT_(*_version2(__version__, n=3))
 

pygeodesy/__main__.py

@@ -5,7 +5,7 @@ u'''Print L{pygeodesy} version, etc. using C{python -m pygeodesy}.
 '''
 
 __all__ = ()
-__version__ = '24.10.14'
+__version__ = '25.04.12'
 
 from os.path import basename, dirname
 
@@ -16,7 +16,7 @@ def _main():  # PYCHOK no cover
         from pygeodesy import constants, _isfrozen, pygeodesy_abspath, version
         from pygeodesy.basics import _xcoverage,_xgeographiclib, _xnumpy, _xscipy
         from pygeodesy.internals import _name_version, printf, _usage, _versions
-        from pygeodesy.interns import NN, _COMMASPACE_, _DEPRECATED_, _DOT_, _DUNDER_all_, \
+        from pygeodesy.interns import NN, _COMMASPACE_, _DALL_, _DEPRECATED_, _DOT_, \
                                      _EQUAL_, _pygeodesy_abspath_, _version_
         from pygeodesy.lazily import _all_deprecates, _all_imports, isLazy
 
@@ -28,7 +28,7 @@ def _main():  # PYCHOK no cover
                                  ('isLazy',                       isLazy),
                                  ('_isfrozen',                 _isfrozen),
                                  ('_floats',      len(constants._floats)),
-                                 (_DUNDER_all_,      len(_all_imports())),
+                                 (_DALL_,            len(_all_imports())),
                                  (_DEPRECATED_,   len(_all_deprecates()))))
 
         def _nv(_xpkg, p):

pygeodesy/albers.py

@@ -15,7 +15,7 @@ and the Albers Conical Equal-Area examples on pp 291-294.
 # make sure int/int division yields float quotient, see .basics
 from __future__ import division as _; del _  # PYCHOK semicolon
 
-from pygeodesy.basics import neg, neg_
+from pygeodesy.basics import _isin, neg, neg_
 from pygeodesy.constants import EPS0, EPS02, _EPSqrt as _TOL, \
                                _0_0, _0_5, _1_0, _N_1_0, _2_0, \
                                _N_2_0, _4_0, _6_0, _90_0, _N_90_0
@@ -38,10 +38,10 @@ from pygeodesy.utily import atan1, atan1d, atan2, degrees360, sincos2, \
 from math import atanh, degrees, fabs, radians, sqrt
 
 __all__ = _ALL_LAZY.albers
-__version__ = '24.11.24'
+__version__ = '25.04.14'
 
 _k1_    = 'k1'
-_NUMIT  =   8  # XXX 4?
+_NUMIT  =   9  # XXX 4?
 _NUMIT0 =  41  # XXX 21?
 _TERMS  =  31  # XXX 16?
 _TOL0   =  sqrt3(_TOL)
@@ -115,7 +115,7 @@ class _AlbersBase(_NamedBase):
         '''(INTERNAL) New C{AlbersEqualArea...} instance.
         '''
         qZ = self._qZ
-        if datum not in (None, self._datum):
+        if not _isin(datum, None, self._datum):
             self._datum = _ellipsoidal_datum(datum, **name)
             qZ = _qZx(self)
         elif qZ is None:
@@ -143,10 +143,8 @@ class _AlbersBase(_NamedBase):
         ca2, ta2 = _ct2(sa2, ca2)
 
         par1 = fabs(ta1 - ta2) < EPS02  # ta1 == ta2
-        if par1 or polar:
-            ta0, C = ta2, _1_0
-        else:
-            ta0, C = self._ta0C2(ca1, sa1, ta1, ca2, sa2, ta2)
+        ta0, C = (ta2, _1_0) if par1 or polar else \
+                  self._ta0C2(ca1, sa1, ta1, ca2, sa2, ta2)
 
         self._lat0 = _Lat(lat0=self._sign * atan1d(ta0))
         self._m02  =  m02 = _1_x21(E.f1 * ta0)
@@ -155,8 +153,8 @@ class _AlbersBase(_NamedBase):
             self._polar = True
 #           self._nrho0 = self._m0 = _0_0
         else:  # m0 = nrho0 / E.a
-            self._m0    = sqrt(m02)
-            self._nrho0 = self._m0 * E.a
+            self._m0    = t = sqrt(m02)
+            self._nrho0 = t * E.a
         t = self._txi0  = self._txif(ta0)
         h = hypot1(t)
         s = self._sxi0  = t / h
@@ -476,12 +474,13 @@ class _AlbersBase(_NamedBase):
         axi, bxi, sxi = self._a_b_sxi3((ca1, sa1, ta1, scb12),
                                        (ca2, sa2, ta2, scb22))
 
-        dsxi = ((esa12 / esa1_2) + _DatanheE(sa2, sa1, E)) * dsn_2 / self._qx
-        C = _Fsum1f_(sxi * dtb12 / dsxi, scb22, scb12).fover(scb22 * scb12 * _2_0)
+        dsxi  = ((esa12 / esa1_2) + _DatanheE(sa2, sa1, E)) * dsn_2 / self._qx
+        C     = _Fsum1f_(sxi * dtb12 / dsxi, scb22, scb12).fover(scb22 * scb12 * _2_0)
 
         S     = _Fsum1f_(sa1, sa2, sa12)
         axi  *= (S * e2 + _1).fover(S + _1, raiser=False)
-        bxi  *= _Fsum1f_(sa1, sa2, esa12).fover(esa1_2) * e2 + _D2atanheE(sa1, sa2, E) * E.e21
+        bxi  *= _Fsum1f_(  sa1, sa2, esa12).fover(esa1_2) * e2 + \
+                _D2atanheE(sa1, sa2, E) * E.e21
         s1_qZ = (axi * self._qZ - bxi) * dsn_2 / dtb12
         ta0   =  self._ta0(s1_qZ, (ta1 + ta2) * _0_5, E)
         return ta0, C
@@ -489,23 +488,19 @@ class _AlbersBase(_NamedBase):
     def _tanf(self, txi):  # in .Ellipsoid.auxAuthalic
         '''(INTERNAL) Function M{tan-phi from tan-xi}.
         '''
-        tol = _tol(_TOL, txi)
-        e2  =  self.ellipsoid.e2
-        qx  =  self._qx
-
-        ta    =  txi
-        _Ta2  =  Fsum(ta).fsum2f_
-        _txif =  self._txif
-        _1    = _1_0
+        e2   =  self.ellipsoid.e2
+        ta   =  txi
+        tol  = _tol(_TOL, ta)
+        _Ta2 =  Fsum(ta).fsum2f_
         for self._iteration in range(1, _NUMIT):  # max 2, mean 1.99
             # dtxi / dta = (scxi / sca)^3 * 2 * (1 - e^2)
             #            / (qZ * (1 - e^2 * sa^2)^2)
-            ta2  =  ta**2
-            sca2 = _1 + ta2
-            txia = _txif(ta)
-            s3qx =  sqrt3(sca2 / (txia**2 + _1)) * qx  # * _1_x21(txia)
-            eta2 = (_1 - e2 * ta2 / sca2)**2
-            ta, d = _Ta2((txi - txia) * s3qx * eta2)
+            ta2   =  ta**2
+            sca2  = _1_0 + ta2
+            txia  =  self._txif(ta)
+            s3    =  sqrt3(sca2 / (txia**2 + _1_0))  # * _1_x21(txia)
+            s3   *= (e2 * ta2 / sca2 - _1_0)**2 * self._qx
+            ta, d = _Ta2((txi - txia) * s3)
             if fabs(d) < tol:
                 return ta
         raise AlbersError(Fmt.no_convergence(d, tol), txt=repr(self))
@@ -546,20 +541,18 @@ class _AlbersBase(_NamedBase):
         '''(INTERNAL) Function M{tan-xi from tan-phi}.
         '''
         E   =  self.ellipsoid
-        _1  = _1_0
-
         ca2 = _1_x21(ta)
         sa  =  sqrt(ca2) * fabs(ta)  # enforce odd parity
-        sa1 = _1 + sa
+        sa1 = _1_0 + sa
 
         es1    =  sa  * E.e2
-        es1m1  =  sa1 * (_1 - es1)
-        es1p1  =  sa1 / (_1 + es1)
-        es2m1  = _1   -  sa * es1
+        es1m1  =  sa1 * (_1_0 - es1)
+        es1p1  =  sa1 / (_1_0 + es1)
+        es2m1  = _1_0 -    sa * es1
         es2m1a =  es2m1 * E.e21  # e2m
-        s =  sqrt((ca2 / (es1p1 * es2m1a) + _atanheE(ca2 / es1m1, E))
-                       * (es1m1 / es2m1a  + _atanheE(es1p1, E)))
-        t = _Fsum1f_(sa / es2m1, _atanheE(sa, E)).fover(s)
+        s  = ca2 / (es1p1 * es2m1a) + _atanheE(ca2 / es1m1, E)
+        s *=        es1m1 / es2m1a  + _atanheE(      es1p1, E)
+        t  = _Fsum1f_( sa / es2m1,    _atanheE(         sa, E)).fover(sqrt(s))
         return neg(t) if ta < 0 else t
 
 

pygeodesy/auxilats/_CX_4.py

@@ -6,13 +6,13 @@ u'''Coefficients for C{_AUXLATITUDE_ORDER} 4 from I{Karney}'s C++ class U{AuxLat
 transcoded to a Python C{_Rdict[auxout][auxin]} of C{_Rtuple}s.
 
 Copyright (C) Charles Karney (2022-2024) Karney@Alum.MIT.edu> and licensed under the
-MIT/X11 License.  For more information, see <https:#GeographicLib.SourceForge.io>.
+MIT/X11 License.  For more information, see <https://GeographicLib.SourceForge.io>.
 '''
 from pygeodesy.auxilats.auxily import Aux
 from pygeodesy.auxilats._CX_Rs import _Rcoeffs, _Rdict, _Rtuple
 
 __all__ = ()
-__version__ = '24.09.03'
+__version__ = '25.04.09'
 
 _coeffs_4 = _Rcoeffs(4, {  # GEOGRAPHICLIB_AUXLATITUDE_ORDER == 4
     Aux.PHI: _Rdict(38,

pygeodesy/auxilats/_CX_6.py

@@ -6,13 +6,13 @@ u'''Coefficients for C{_AUXLATITUDE_ORDER} 6 from I{Karney}'s C++ class U{AuxLat
 transcoded to a Python C{_Rdict[auxout][auxin]} of C{_Rtuple}s.
 
 Copyright (C) Charles Karney (2022-2024) Karney@Alum.MIT.edu> and licensed under the
-MIT/X11 License.  For more information, see <https:#GeographicLib.SourceForge.io>.
+MIT/X11 License.  For more information, see <https://GeographicLib.SourceForge.io>.
 '''
 from pygeodesy.auxilats.auxily import Aux
 from pygeodesy.auxilats._CX_Rs import _Rcoeffs, _Rdict, _Rtuple
 
 __all__ = ()
-__version__ = '24.09.02'
+__version__ = '25.04.09'
 
 _coeffs_6 = _Rcoeffs(6, {  # GEOGRAPHICLIB_AUXLATITUDE_ORDER == 6
     Aux.PHI: _Rdict(78,

pygeodesy/auxilats/_CX_8.py

@@ -6,13 +6,13 @@ u'''Coefficients for C{_AUXLATITUDE_ORDER} 8 from I{Karney}'s C++ class U{AuxLat
 transcoded to a Python C{_Rdict[auxout][auxin]} of C{_Rtuple}s.
 
 Copyright (C) Charles Karney (2022-2024) Karney@Alum.MIT.edu> and licensed under the
-MIT/X11 License.  For more information, see <https:#GeographicLib.SourceForge.io>.
+MIT/X11 License.  For more information, see <https://GeographicLib.SourceForge.io>.
 '''
 from pygeodesy.auxilats.auxily import Aux
 from pygeodesy.auxilats._CX_Rs import _Rcoeffs, _Rdict, _Rtuple
 
 __all__ = ()
-__version__ = '24.09.03'
+__version__ = '25.04.09'
 
 _coeffs_8 = _Rcoeffs(8, {  # GEOGRAPHICLIB_AUXLATITUDE_ORDER == 8
     Aux.PHI: _Rdict(132,

pygeodesy/auxilats/_CX_Rs.py

@@ -6,13 +6,13 @@ Python versions of coefficients from I{Karney}'s C++ class U{AuxLatitude
 <https://GeographicLib.SourceForge.io/C++/doc/classGeographicLib_1_1AuxLatitude.html>}.
 
 Copyright (C) Charles Karney (2022-2024) Karney@Alum.MIT.edu> and licensed under the
-MIT/X11 License.  For more information, see <https:#GeographicLib.SourceForge.io>.
+MIT/X11 License.  For more information, see <https://GeographicLib.SourceForge.io>.
 '''
 # make sure int/int division yields float quotient, see .basics
 from __future__ import division as _; del _  # PYCHOK semicolon
 
 # from pygeodesy.basics import _splituple  # _MODS
-from pygeodesy.constants import _floats as _Cs
+from pygeodesy.constants import _floats as _constants_floats
 from pygeodesy.errors import _AssertionError,  _MODS
 # from pygeodesy.internals import _sizeof  # _MODS
 from pygeodesy.interns import NN, MISSING, _COMMA_, _duplicate_, _NL_, \
@@ -22,7 +22,7 @@ from pygeodesy.named import ADict,  Property_RO
 # from pygeodesy.props import Property_RO  # from .named
 
 __all__ = ()
-__version__ = '25.04.08'
+__version__ = '25.04.11'
 
 
 class _Rcoeffs(ADict):
@@ -58,7 +58,7 @@ class _Rcoeffs(ADict):
             z += _zB(t)  # Float
             # assert R.Rdict is None
             n +=  len(t)
-            u +=  sum(1 for f in t if f in _Cs)
+            u +=  sum(1 for f in t if f in _constants_floats)
         return b, n, (n - u), z
 
     def items(self):  # string-ify keys  # PYCHOK no cover
@@ -74,7 +74,7 @@ class _Rcoeffs(ADict):
 
     def _validate(self, aL, lenAux):
         # in .auxily.Aux._CXcoeffs(al, Aux.len(aL))
-        a, n = self.ALorder, self.n  # PYCHOK Adict!
+        a, n = self.ALorder, self.n  # PYCHOK ADict!
 #       for R in self._Rtuples():
 #           assert isinstance(R, _Rtuple)
         if aL != a or lenAux != n:
@@ -110,7 +110,7 @@ class _Rdict(dict):  # in ._CX_#, .auxLat, .rhumb.aux_
         def _p_q(p=NN, q=1, *x):
             return (NN if x else p), q
 
-        _get = _Cs.get
+        _get = _constants_floats.get
         for r in _MODS.basics._splituple(rs):
             p, q = _p_q(*r.split(_SLASH_))
             if p:
@@ -125,7 +125,7 @@ class _Rdict(dict):  # in ._CX_#, .auxLat, .rhumb.aux_
         # return a tuple of floats from an C{_Rtuple}
         k, n, rs = Rtuple.k_n_rs
         t = tuple(f for m in map(self._floats, rs)
-                        for f in m)  # ... yield f
+                        for f in m)  # == yield f
         # @see: <https://StackOverflow.com/questions/10632839/>
         #       and <https://realPython.com/python-flatten-list/>
         if len(t) != n:
@@ -188,8 +188,8 @@ class _Rtuple(list):  # MUST be list, NOT tuple!
         return self._tuple[i]
 
     if _MODS.sys_version_info2 < (3, 0):
-        def __getslice__(self, i, j):
-            return self._tuple[slice(i, j)]
+        def __getslice__(self, *i_j):  # PYCHOK 3 args
+            return self._tuple[slice(*i_j)]
 
     def __iter__(self):
         return iter(self._tuple)

pygeodesy/auxilats/__init__.py

@@ -14,8 +14,8 @@ License.  For more information, see the U{GeographicLib<https://GeographicLib.So
 @note: Class L{AuxDST} requires U{numpy<https://PyPI.org/project/numpy>} to be installed, version 1.16
        or newer.
 
-@see: U{Auxiliary latitudes<https:#GeographicLib.SourceForge.io/C++/doc/auxlat.html>} and
-      U{On auxiliary latitudes<https:#ArXiv.org/abs/2212.05818>}.
+@see: U{Auxiliary latitudes<https://GeographicLib.SourceForge.io/C++/doc/auxlat.html>} and
+      U{On auxiliary latitudes<https://ArXiv.org/abs/2212.05818>}.
 '''
 
 from pygeodesy.auxilats.auxAngle import AuxAngle, AuxBeta, AuxChi, AuxMu, \
@@ -29,7 +29,7 @@ from pygeodesy.lazily import _ALL_OTHER
 # no modules: auxAngle, auxDLat, auxDST, auxily, auxLat
 __all__ = _ALL_OTHER(Aux, AuxAngle, AuxDLat, AuxDST, AuxLat,
                      AuxBeta, AuxChi, AuxMu, AuxPhi, AuxTheta, AuxXi)
-__version__ = '25.01.15'
+__version__ = '25.04.14'
 
 # **) MIT License
 #

pygeodesy/auxilats/__main__.py

@@ -5,23 +5,22 @@ u'''Print L{auxilats} version, etc. using C{python -m pygeodesy.auxilats}.
 '''
 
 __all__ = ()
-__version__ = '25.04.04'
+__version__ = '25.04.14'
 
 
 def _main(**ALorder):  # PYCHOK no cover
 
     try:
-        from pygeodesy import auxilats
+        from pygeodesy import ADict, auxilats
         from pygeodesy.internals import _fper, _name_version, \
                                          printf, _versions
-        from pygeodesy.interns import _COMMASPACE_, _EQUAL_
+        from pygeodesy.interns import _COMMASPACE_
 
         A  = auxilats.AuxLat(**ALorder)
         Cx = A._CXcoeffs  # PropertyRO: Adict of _Rdicts
         b, n, u, z = Cx.bnuz4()
-        p = dict(ALorder=A.ALorder,CXb=b, CXb_z=_fper(b, z),
-                            CXn=n, CXu=u, CXu_n=_fper(u, n))
-        p = list(_EQUAL_(*t) for t in p.items())
+        p  = ADict(ALorder=A.ALorder, CXb=b, CXb_z=_fper(b, z),
+                               CXn=n, CXu=u, CXu_n=_fper(u, n))._toL()
         try:
             import geographiclib
             p.append(_name_version(geographiclib))
@@ -39,6 +38,9 @@ def _main(**ALorder):  # PYCHOK no cover
 from sys import argv  # .internals._isPyChOK
 _main(ALorder=int(argv[1])) if len(argv) == 2 and argv[1].isdigit() else _main()
 
+# % python3.13 -m pygeodesy.auxilats
+# pygeodesy.auxilats 25.04.14: ALorder=6, CXb=11099, CXb_z=64.1%, CXn=522, CXu=448, CXu_n=85.8%, geographiclib 2.0 (pygeodesy 25.4.24 Python 3.13.3 64bit arm64 macOS 15.4)
+
 # % python3.12 -m pygeodesy.auxilats 8
 # pygeodesy.auxilats 24.09.04: ALorder=8, CXb=20310, CXb_z=71.5%, CXn=888, CXu=780, CXu_n=87.8%, geographiclib 2.0 (pygeodesy 24.9.9 Python 3.12.5 64bit arm64 macOS 14.6.1)
 

pygeodesy/auxilats/auxAngle.py

@@ -32,7 +32,7 @@ from pygeodesy.utily import atan2, atan2d, sincos2, sincos2d
 from math import asinh, copysign, degrees, fabs, radians, sinh
 
 __all__ = ()
-__version__ = '24.11.24'
+__version__ = '25.04.14'
 
 _0_INF_NAN_NINF = (0, _0_0) + _INF_NAN_NINF
 _MAX_2          =  MAX * _0_5  # PYCHOK used!
@@ -207,7 +207,7 @@ class AuxAngle(_Named):
     def _copy_2(self, which):
         '''(INTERNAL) Copy for I{dyadic} operators.
         '''
-        return _Named.copy(self, deep=False, name=which.__name__)
+        return _Named.copy(self, deep=False, name__=which)
 
     def _copy_r2(self, other, which):
         '''(INTERNAL) Copy for I{reverse-dyadic} operators.

pygeodesy/auxilats/auxLat.py

@@ -9,8 +9,8 @@ Copyright (C) U{Charles Karney<mailto:Karney@Alum.MIT.edu>} (2022-2024) and lice
 under the MIT/X11 License.  For more information, see the U{GeographicLib
 <https://GeographicLib.SourceForge.io>} documentation.
 
-@see: U{Auxiliary latitudes<https:#GeographicLib.SourceForge.io/C++/doc/auxlat.html>}
-      U{On auxiliary latitudes<https:#ArXiv.org/abs/2212.05818>}.
+@see: U{Auxiliary latitudes<https://GeographicLib.SourceForge.io/C++/doc/auxlat.html>}
+      U{On auxiliary latitudes<https://ArXiv.org/abs/2212.05818>}.
 '''
 # make sure int/int division yields float quotient, see .basics
 from __future__ import division as _; del _  # PYCHOK semicolon
@@ -19,7 +19,7 @@ from pygeodesy.auxilats.auxAngle import AuxAngle, AuxBeta, AuxChi, _AuxClass, \
                                         AuxMu, AuxPhi, AuxTheta, AuxXi
 from pygeodesy.auxilats.auxily import Aux, _sc, _sn
 from pygeodesy.auxilats._CX_Rs import _Rdict, _Rkey, _Rtuple
-from pygeodesy.basics import _reverange, _xinstanceof,  _passarg
+from pygeodesy.basics import _isin, _reverange, _xinstanceof,  _passarg
 from pygeodesy.constants import INF, MAX_EXP, MIN_EXP, NAN, PI_2, PI_4, _EPSqrt, \
                                _0_0, _0_0s, _0_1, _0_5, _1_0, _2_0, _3_0, _360_0, \
                                _log2, _over, isfinite, isinf, isnan
@@ -48,7 +48,7 @@ except ImportError:  # Python 3.11-
         return pow(_2_0, x)
 
 __all__ = ()
-__version__ = '25.01.15'
+__version__ = '25.04.14'
 
 _TRIPS = 1024  # XXX 2 or 3?
 
@@ -90,7 +90,7 @@ class AuxLat(AuxAngle):
         else:
             name = NN
         try:
-            if a_earth not in (_EWGS84, _WGS84):
+            if not _isin(a_earth, _EWGS84, _WGS84):
                 n = _name__(name, name__=AuxLat)
                 if b is f is None:
                     E = _ellipsoidal_datum(a_earth, name=n).ellipsoid  # XXX raiser=_earth_

pygeodesy/auxilats/auxily.py

@@ -16,10 +16,12 @@ from __future__ import division as _; del _  # PYCHOK semicolon
 
 # from pygeodesy import auxilats  # _MODS
 from pygeodesy.auxilats._CX_Rs import _Rkey
+from pygeodesy.basics import _isin,  typename
 from pygeodesy.constants import INF, NAN, isinf, isnan, _0_0, _0_5, _1_0, \
                                _copysign_1_0, _over, _1_over
 from pygeodesy.errors import AuxError
 from pygeodesy.fmath import hypot1 as _sc, hypot2_
+# from pygeodesy.internals import typename  # from .basics
 from pygeodesy.interns import NN,  _DOT_, _UNDER_  # PYCHOK used!
 from pygeodesy.lazily import _ALL_DOCS,  _ALL_MODS as _MODS  # PYCHOK used!
 from pygeodesy.utily import atan1
@@ -27,7 +29,7 @@ from pygeodesy.utily import atan1
 from math import asinh, copysign
 
 __all__ = ()
-__version__ = '25.01.15'
+__version__ = '25.04.14'
 
 
 class Aux(object):
@@ -58,7 +60,7 @@ class Aux(object):
             _coeffs = Aux._coeffs[aL]
         except KeyError:
             try:  # from pygeodesy.auxilats._CX_x import _coeffs_x as _coeffs
-                _CX_x   = _DOT_(_MODS.auxilats.__name__, _UNDER_('_CX', aL))
+                _CX_x   = _DOT_(typename(_MODS.auxilats), _UNDER_('_CX', aL))
                 _coeffs = _MODS.getattr(_CX_x, _UNDER_('_coeffs', aL))
             except (AttributeError, ImportError, KeyError, TypeError) as x:
                 raise AuxError(ALorder=aL, cause=x)
@@ -92,7 +94,7 @@ class Aux(object):
         return (auxout - auxin) if max(auxin, auxout) < Aux.MU else None
 
     def use_n2(self, aux):
-        return aux not in (Aux.CHI, Aux.XI)
+        return not _isin(aux, Aux.CHI, Aux.XI)
 
 Aux = Aux()  # PYCHOK singleton
 

pygeodesy/azimuthal.py

@@ -42,19 +42,20 @@ altitude in Earth radii<https://WikiPedia.org/wiki/Azimuthal_equidistant_project
 # make sure int/int division yields float quotient, see .basics
 from __future__ import division as _; del _  # PYCHOK semicolon
 
-# from pygeodesy.basics import _xinstanceof  # from .ellipsoidalBase
+# from pygeodesy.basics import _isin, _xinstanceof  # from .ellipsoidalBase
 from pygeodesy.constants import EPS, EPS0, EPS1, NAN, isnon0, _umod_360, \
                                _EPStol, _0_0, _0_1, _0_5, _1_0, _N_1_0, _2_0
 from pygeodesy.ellipsoidalBase import LatLonEllipsoidalBase as _LLEB, \
-                              _xinstanceof
+                                     _isin, _xinstanceof
 from pygeodesy.datums import _spherical_datum, _WGS84
 from pygeodesy.errors import _ValueError, _xdatum, _xkwds
 from pygeodesy.fmath import euclid, fdot_, hypot as _hypot,  Fsum
 # from pygeodesy.fsums import Fsum  # from .fmath
 # from pygeodesy.formy import antipode  # _MODS
+# from pygeodesy.internals import typename  # from .karney
 from pygeodesy.interns import _azimuth_, _datum_, _lat_, _lon_, _scale_, \
                               _SPACE_, _x_, _y_
-from pygeodesy.karney import _norm180
+from pygeodesy.karney import _norm180,  typename
 from pygeodesy.latlonBase import _MODS, LatLonBase as _LLB
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _FOR_DOCS  # ALL_MODS
 from pygeodesy.named import _name__, _name2__, _NamedBase, _NamedTuple, _Pass
@@ -70,7 +71,7 @@ from pygeodesy.utily import asin1, atan1, atan2, atan2b, atan2d, \
 from math import acos, degrees, fabs, sin, sqrt
 
 __all__ = _ALL_LAZY.azimuthal
-__version__ = '24.11.24'
+__version__ = '25.04.14'
 
 _EPS_K         = _EPStol * _0_1  # Karney's eps_ or _EPSmin * _0_1?
 _over_horizon_ = 'over horizon'
@@ -112,7 +113,7 @@ class _AzimuthalBase(_NamedBase):
 
            @raise TypeError: Invalid B{C{datum}}.
        '''
-        if datum not in (None, self._datum):
+        if not _isin(datum, None, self._datum):
             self._datum = _spherical_datum(datum, **name)
         if name:
             self.name = name
@@ -1088,7 +1089,7 @@ class Stereographic(_AzimuthalBase):
     def k0(self, factor):
         '''Set the central scale factor (C{scalar}).
         '''
-        n = Stereographic.k0.fget.__name__  # 'k0', name__=Stereographic.k0.fget
+        n = typename(Stereographic.k0.fget)  # 'k0', name__=Stereographic.k0.fget
         self._k0  = Scalar_(factor, name=n, low=EPS, high=2)  # XXX high=1, 2, other?
         self._k02 = self._k0 * _2_0