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

pygeodesy/__init__.py

@@ -17,7 +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.
+by overlap and in 3-D), 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>},
@@ -51,9 +51,9 @@ resources/hull:8338>} and the U{Reumann-Witkam<https://psimpl.SourceForge.net/re
 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
+L{height<pygeodesy.heights>} of C{LatLon} points and L{Geoid<pygeodesy.geoids>} models, to 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.
+distances or to perform I{boolean} operations between (composite) polygons of C{LatLon} points.
 
 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
@@ -125,7 +125,11 @@ C{epydoc --html --no-private --no-source --name=pygeodesy --url=... -v pygeodesy
 Tests
 =====
 
-The tests ran with Python 3.13.5 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0),
+The tests ran with Python 3.14 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.1),
+Python 3.13.7 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.1),
+U{numpy<https://PyPI.org/project/numpy>} 2.3.3, U{scipy<https://PyPI.org/project/scipy>} 1.16.2,
+U{GeoConvert<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5 and
+U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5),
 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.5,
@@ -133,30 +137,26 @@ U{IntersectTool<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.
 U{RhumbSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5),
 Python 3.11.5 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0,
 U{numpy<https://PyPI.org/project/numpy>} 1.24.2 and U{scipy<https://PyPI.org/project/scipy>} 1.10.1),
-Python 3.10.8 (with U{geographiclib <https://PyPI.org/project/geographiclib>} 2.0,
-U{numpy<https://PyPI.org/project/numpy>} 1.23.3, U{scipy<https://PyPI.org/project/scipy>} 1.9.1,
-U{GeoConvert<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5,
-U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5) and
-Python 2.7.18 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 1.50,
+and with Python 2.7.18 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 1.50,
 U{numpy<https://PyPI.org/project/numpy>} 1.16.6, U{scipy<https://PyPI.org/project/scipy>} 1.2.2,
 U{GeoConvert<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5,
 U{GeodSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5,
 U{IntersectTool<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5 and
 U{RhumbSolve<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.5), all in 64-bit on
-macOS 15.6.1 Sequoia.
+macOS 26.0.1 Tahoe.
 
 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.5, 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.
+Test coverage has been measured with U{coverage<https://PyPI.org/project/coverage>} 7.10.7 using Python
+3.14, 3.13.7 and 3.12.7.  The complete coverage report in HTML and a PDF summary are included in the
+distribution files.
 
-Python 3.13.5, 3.12.7, 3.11.5 and 3.10.8 run on Apple M4 Si (C{arm64}), I{natively}.  Python 2.7.18 runs
+Python 3.14, 3.13.7, 3.12.7 and 3.11.5 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.5 (and U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0) on
+The tests also ran with Python 3.13.7 (and U{geographiclib<https://PyPI.org/project/geographiclib>} 2.1) on
 U{Debian 12<https://Cirrus-CI.com/github/mrJean1/PyGeodesy/master>} in 64-bit only, with Python 3.12.8 (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
@@ -166,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-4, 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-6, 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,
@@ -175,7 +175,7 @@ U{PyPy<https://PyPy.org>} 7.3.12 (Python 3.10.12), 7.3.1 (Python 3.6.9) and U{Py
 1.16.5, 1.16.2, 1.15.2, 1.14.0, 1.13.1, 1.8.0rc1 or 1.6.2 and U{scipy<https://PyPI.org/project/scipy>} 1.5.0), U{PyPy
 <https://PyPy.org>} 7.3.0 (Python 2.7.13 and 3.6.9), U{PyPy<https://PyPy.org>} 6.0.0 (Python 2.7.13 and 3.5.3)
 and U{Intel-Python<https://software.Intel.com/en-us/distribution-for-python>} 3.5.3 (and U{numpy
-<https://PyPI.org/project/numpy>} 1.11.3) on macOS 15.0-5 Sequoia, 14.0-6.1 Sonoma, 13.0-5.2 Ventura, 12.1-6 Monterey,
+<https://PyPI.org/project/numpy>} 1.11.3) on macOS 15.0-6 Sequoia, 14.0-6.1 Sonoma, 13.0-5.2 Ventura, 12.1-6 Monterey,
 11.0-5.2-6.1 Big Sur (aka 10.16), 10.15.3, 10.15.5-7 Catalina, 10.14 Mojave, 10.13.6 High Sierra and 10.12 Sierra,
 MacOS X 10.11 El Capitan and/or MacOS X 10.10 Yosemite, with U{Pythonista<https://OMZ-Software.com/pythonista>}3.2
 (with geographiclib 1.50 or 1.49 and numpy 1.8.0) on iOS 14.4.2, 11.4.1, 12.0-3 on iPad4, iPhone6, iPhone10 and/or
@@ -187,10 +187,10 @@ Notes
 =====
 
 All Python source code has been statically U{checked<https://GitHub.com/ActiveState/code/tree/master/recipes/Python/
-546532_PyChecker_postprocessor>} with U{Ruff<https://GitHub.com/astral-sh/ruff>} using Python 3.13.5 and with
+546532_PyChecker_postprocessor>} with U{Ruff<https://GitHub.com/astral-sh/ruff>} using Python 3.13.7 and 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, both in 64-bit on macOS 15.6.1 Sequoia.
+2.7.18, both in 64-bit on macOS 26.0.1 Tahoe.
 
 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>}.
@@ -329,6 +329,7 @@ SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.}
 @var PI_3:  One third PI, M{PI / 3} (C{float}).
 @var PI4:   Four PI, M{PI * 4} (C{float}).
 @var PI_4:  Quarter PI, M{PI / 4} (C{float}).
+@var PI_6:  One sixth PI, M{PI / 6} (C{float}).
 
 @var R_MA: Equatorial earth radius (C{meter}), WGS84, EPSG:3785.
 @var R_MB: Polar earth radius (C{meter}), WGS84, EPSG:3785.
@@ -604,7 +605,7 @@ else:
 
 from pygeodesy.internals import _version2,  _DOT_  # noqa: E402
 # from pygeodesy.interns import _DOT_  # from .internals
-__version__ = '25.08.25'
+__version__ = '25.10.10'
 # see setup.py for similar logic
 version     = _DOT_(*_version2(__version__, n=3))
 

pygeodesy/auxilats/__init__.py

@@ -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.04.14'
+__version__ = '25.08.28'
 
 # **) MIT License
 #

pygeodesy/auxilats/auxAngle.py

@@ -17,7 +17,7 @@ from __future__ import division as _; del _  # noqa: E702 ;
 from pygeodesy.auxilats.auxily import Aux, _Aux2Greek,  AuxError
 from pygeodesy.basics import map1, map2, _xinstanceof
 from pygeodesy.constants import EPS, _INF_NAN_NINF, MAX, NAN, _0_0, _0_5, _1_0, \
-                               _copysign_1_0, _over, _pos_self, isfinite, isnan
+                               _copysign_1_0, isfinite, isnan, _over, _pos_self
 # from pygeodesy.errors import AuxError  # from .auxilats.auxily
 from pygeodesy.fmath import hypot,  unstr
 from pygeodesy.fsums import _add_op_, _iadd_op_, _isub_op_, _sub_op_
@@ -32,7 +32,7 @@ from pygeodesy.utily import atan2, atan2d, sincos2, sincos2d
 from math import asinh, copysign, degrees, fabs, radians, sinh
 
 __all__ = ()
-__version__ = '25.05.12'
+__version__ = '25.08.31'
 
 _0_INF_NAN_NINF = (0, _0_0) + _INF_NAN_NINF
 _MAX_2          =  MAX * _0_5  # PYCHOK used!
@@ -505,6 +505,7 @@ _AUXClass = {Aux.BETA:  AuxBeta,
              Aux.CHI:   AuxChi,
              Aux.MU:    AuxMu,
              Aux.PHI:   AuxPhi,
+#            Aux.PSI:   AuxPSI,  # Isometric
              Aux.THETA: AuxTheta,
              Aux.XI:    AuxXi}
 
@@ -525,7 +526,7 @@ def _yx2(yx):
     return y, x
 
 
-__all__ += _ALL_DOCS(AuxAngle, AuxBeta, AuxChi, AuxMu, AuxPhi, AuxTheta, AuxXi)
+__all__ += _ALL_DOCS(AuxAngle, *_AUXClass.values())
 
 # **) MIT License
 #

pygeodesy/auxilats/auxily.py

@@ -29,7 +29,7 @@ from pygeodesy.utily import atan1
 from math import asinh, copysign
 
 __all__ = ()
-__version__ = '25.05.12'
+__version__ = '25.08.31'
 
 
 class Aux(object):

pygeodesy/azimuthal.py

@@ -49,7 +49,7 @@ from pygeodesy.ellipsoidalBase import LatLonEllipsoidalBase as _LLEB, \
                                      _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.fmath import euclid, hypot as _hypot,  Fsum
 # from pygeodesy.fsums import Fsum  # from .fmath
 # from pygeodesy.formy import antipode  # _MODS
 # from pygeodesy.internals import typename  # from .karney
@@ -71,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__ = '25.05.12'
+__version__ = '25.09.29'
 
 _EPS_K         = _EPStol * _0_1  # Karney's eps_ or _EPSmin * _0_1?
 _over_horizon_ = 'over horizon'
@@ -147,7 +147,7 @@ class _AzimuthalBase(_NamedBase):
         '''I{Must be overloaded}.'''
         self._notOverloaded(lat, lon, **name)
 
-    def _forward(self, lat, lon, name, _k_t_2):
+    def _forward(self, lat, lon, name, _k_t):
         '''(INTERNAL) Azimuthal (spherical) forward C{lat, lon} to C{x, y}.
         '''
         lat, lon = Lat_(lat), Lon_(lon)
@@ -155,16 +155,16 @@ class _AzimuthalBase(_NamedBase):
         s0, c0 = self._sc0
 
         cb  *= ca
-        k, t = _k_t_2(fdot_(s0, sa, c0, cb))
+        k, t = _k_t(c0 * cb + s0 * sa)
         if t:
             r = k * self.radius
-            y = r * fdot_(c0, sa, -s0, cb)
+            y = r * (c0 * sa - s0 * cb)
             e, n, z, _ = _enzh4(r * sb * ca, y, None)
         else:  # 0 or 180
             e = n = z = _0_0
 
         t = Azimuthal7Tuple(e, n, lat, lon, z, k, self.datum,
-                                       name=self._name__(name))
+                                                  name=self._name__(name))
         return t
 
     def _forwards(self, *lls):
@@ -258,7 +258,7 @@ class _AzimuthalBase(_NamedBase):
             if lea or fabs(c0) > EPS:
                 d = atan2d(e * sc, r * c0 * cc - n * s0 * sc)
             else:
-                d = atan2d(e, (n if s0 < 0 else -n))
+                d = atan2d(e,     (n if s0 < 0 else -n))
             lon = _norm180(self.lon0 + d)
 
         if LatLon is None:
@@ -745,8 +745,8 @@ class _GnomonicBase(_AzimuthalGeodesic):
 
            @arg lat: Latitude of the location (C{degrees90}).
            @arg lon: Longitude of the location (C{degrees180}).
-           @kwarg raiser: Do or don't throw an error (C{bool}) if
-                          the location lies over the horizon.
+           @kwarg raiser: If C{False}, do not throw an error if the location lies
+                          over the horizon (C{bool}).
            @kwarg name: Optional C{B{name}=NN} for the location (C{str}).
 
            @return: An L{Azimuthal7Tuple}C{(x, y, lat, lon, azimuth, scale, datum)}
@@ -760,8 +760,6 @@ class _GnomonicBase(_AzimuthalGeodesic):
            @raise AzimuthalError: Invalid B{C{lat}}, B{C{lon}} or the location lies
                                   over the horizon and C{B{raiser}=True}.
         '''
-        self._iteration = 0
-
         r = self.geodesic.Inverse(self.lat0, self.lon0,
                                   Lat_(lat), Lon_(lon), outmask=self._mask)
         M = r.M21
@@ -776,7 +774,7 @@ class _GnomonicBase(_AzimuthalGeodesic):
             e = n = NAN
 
         t = self._7Tuple(e, n, r, _name__(name), M=M)
-        t._iteraton = 0
+        t._iteration = self._iteration = 0
         return t
 
     def reverse(self, x, y, LatLon=None, **name_LatLon_kwds):  # PYCHOK signature

pygeodesy/basics.py

@@ -37,7 +37,7 @@ from math import copysign as _copysign
 # import inspect as _inspect  # _MODS
 
 __all__ = _ALL_LAZY.basics
-__version__ = '25.04.14'
+__version__ = '25.09.04'
 
 _below_           = 'below'
 _list_tuple_types = (list, tuple)
@@ -201,7 +201,7 @@ def copytype(x, y):
 
 
 def _enumereverse(iterable):
-    '''(INTERNAL) Reversed C{enumberate}.
+    '''(INTERNAL) Reversed C{enumerate}.
     '''
     for j in _reverange(len(iterable)):
         yield j, iterable[j]
@@ -714,11 +714,11 @@ def _reverange(n, stop=-1, step=-1):
 
 
 def signBit(x):
-    '''Return C{signbit(B{x})}, like C++.
+    '''Return C{signbit(B{x})}, like C++, see also L{isneg}.
 
        @return: C{True} if C{B{x} < 0} or C{NEG0} (C{bool}).
     '''
-    return x < 0 or _MODS.constants.isneg0(x)
+    return (x or _copysign(1, x)) < 0
 
 
 def _signOf(x, ref):  # in .fsums

pygeodesy/booleans.py

@@ -45,7 +45,7 @@ from pygeodesy.utily import fabs, _unrollon, _Wrap
 # from math import fabs  # from .utily
 
 __all__ = _ALL_LAZY.booleans
-__version__ = '25.05.26'
+__version__ = '25.09.04'
 
 _0EPS  =  EPS  # near-zero, positive
 _EPS0  = -EPS  # near-zero, negative
@@ -60,26 +60,26 @@ _corners_ = 'corners'
 _open_    = 'open'
 
 
-def _Enum(txt, enum):  # PYCHOK unused
+def _Cnum(txt, enum):  # PYCHOK unused
     return txt  # NN(txt, _TILDE_, enum)
 
 
 class _L(object):  # Intersection labels
-    CROSSING      = _Enum(_X_, 1)  # C++ enum
-    CROSSING_D    = _Enum(_X_ + _d_, 8)
+    CROSSING      = _Cnum(_X_, 1)  # C++ enum
+    CROSSING_D    = _Cnum(_X_ + _d_, 8)
     CROSSINGs     = (CROSSING, CROSSING_D)
-    BOUNCING      = _Enum(_B_, 2)
-    BOUNCING_D    = _Enum(_B_ + _d_, 9)
+    BOUNCING      = _Cnum(_B_, 2)
+    BOUNCING_D    = _Cnum(_B_ + _d_, 9)
     BOUNCINGs     = (BOUNCING, BOUNCING_D) + CROSSINGs
-    LEFT_ON       = _Enum('Lo', 3)
-    ON_ON         = _Enum('oo', 5)
-    ON_LEFT       = _Enum('oL', 6)
-    ON_RIGHT      = _Enum('oR', 7)
-    RIGHT_ON      = _Enum('Ro', 4)
+    LEFT_ON       = _Cnum('Lo', 3)
+    ON_ON         = _Cnum('oo', 5)
+    ON_LEFT       = _Cnum('oL', 6)
+    ON_RIGHT      = _Cnum('oR', 7)
+    RIGHT_ON      = _Cnum('Ro', 4)
     RIGHT_LEFT_ON = (RIGHT_ON, LEFT_ON)
     # Entry/Exit flags
-    ENTRY         = _Enum(_e_, 1)
-    EXIT          = _Enum(_x_, 0)
+    ENTRY         = _Cnum(_e_, 1)
+    EXIT          = _Cnum(_x_, 0)
     Toggle        = {ENTRY: EXIT,
                      EXIT: ENTRY,
                      None: None}
@@ -88,10 +88,10 @@ _L = _L()  # PYCHOK singleton
 
 
 class _RP(object):  # RelativePositions
-    IS_Pm = _Enum('Pm', 2)  # C++ enum
-    IS_Pp = _Enum('Pp', 3)
-    LEFT  = _Enum('L',  0)
-    RIGHT = _Enum(_R_,  1)
+    IS_Pm = _Cnum('Pm', 2)  # C++ enum
+    IS_Pp = _Cnum('Pp', 3)
+    LEFT  = _Cnum('L',  0)
+    RIGHT = _Cnum(_R_,  1)
 
 _RP = _RP()  # PYCHOK singleton
 
@@ -1521,16 +1521,16 @@ class _CompositeGH(_CompositeBase):
 class _EdgeFHP(object):
     # An edge between two L{LatLonFHP} points.
 
-    X_INTERSECT = _Enum('Xi', 1)  # C++ enum
-    X_OVERLAP   = _Enum('Xo', 5)
-    P_INTERSECT = _Enum('Pi', 3)
-    P_OVERLAP   = _Enum('Po', 7)
+    X_INTERSECT = _Cnum('Xi', 1)  # C++ enum
+    X_OVERLAP   = _Cnum('Xo', 5)
+    P_INTERSECT = _Cnum('Pi', 3)
+    P_OVERLAP   = _Cnum('Po', 7)
     Ps          = (P_INTERSECT, P_OVERLAP, X_OVERLAP)
-    Q_INTERSECT = _Enum('Qi', 2)
-    Q_OVERLAP   = _Enum('Qo', 6)
+    Q_INTERSECT = _Cnum('Qi', 2)
+    Q_OVERLAP   = _Cnum('Qo', 6)
     Qs          = (Q_INTERSECT, Q_OVERLAP, X_OVERLAP)
-    V_INTERSECT = _Enum('Vi', 4)
-    V_OVERLAP   = _Enum('Vo', 8)
+    V_INTERSECT = _Cnum('Vi', 4)
+    V_OVERLAP   = _Cnum('Vo', 8)
     Vs          = (V_INTERSECT, V_OVERLAP)
 
     def __init__(self, p1, p2, **unused):

pygeodesy/constants.py

@@ -10,7 +10,7 @@ L{pygeodesy.isnon0} and L{pygeodesy.remainder}.
 # make sure int/int division yields float quotient, see .basics
 from __future__ import division as _; del _  # noqa: E702 ;
 
-from pygeodesy.basics import _copysign, isbool, iscomplex, isint
+from pygeodesy.basics import _copysign, isbool, iscomplex, isint, signBit
 from pygeodesy.errors import _xError, _xError2, _xkwds_get1, _xkwds_item2
 # from pygeodesy.fsums import _isFsum_2Tuple  # _MODS
 from pygeodesy.internals import _0_0, _100_0, typename
@@ -26,25 +26,31 @@ except ImportError:  # Python 2-
     _inf, _nan = float(_INF_), float(_NAN_)
 
 __all__ = _ALL_LAZY.constants
-__version__ = '25.05.26'
+__version__ = '25.09.11'
 
 
 def _copysign_0_0(y):
     '''(INTERNAL) copysign(0.0, y), only C{float}.
     '''
-    return _N_0_0 if y < 0 else _0_0
+    return _N_0_0 if signBit(y) else _0_0
 
 
 def _copysign_1_0(y):
     '''(INTERNAL) copysign(1.0, y), only C{float}.
     '''
-    return _N_1_0 if y < 0 else _1_0
+    return _N_1_0 if signBit(y) else _1_0
 
 
 def _copysignINF(y):
     '''(INTERNAL) copysign(INF, y), only C{float}.
     '''
-    return NINF if y < 0 else INF
+    return NINF if signBit(y) else INF
+
+
+def _flipsign(x, y):
+    '''(INTERNAL) Negate C{x} for negative C{y}.
+    '''
+    return (-x) if signBit(y) else x
 
 
 def _Float(**name_arg):
@@ -137,7 +143,7 @@ def _naninf(p, *xs):
 
 
 def _over(p, q):
-    '''(INTERNAL) Return C{B{p} / B{q}} avoiding C{ZeroDivisionError} exceptions.
+    '''(INTERNAL) Return C{B{p} / B{q}} without C{ZeroDivisionError} exceptions.
     '''
     try:
         return (p / q)  # if p else _copysign_0_0(q)
@@ -145,6 +151,18 @@ def _over(p, q):
         return (_copysignINF(p) if isfinite(p) else NAN) if p else p
 
 
+def _over_1(p, q):
+    '''(INTERNAL) Return C{B{p} / B{q}} with exact C{1.0} and without C{ZeroDivisionError} exceptions.
+    '''
+    if fabs(p) != fabs(q):
+        r = _over(p, q)
+    else:
+        r = _flipsign(p, q)
+        if p:
+            r = _copysign_1_0(r)
+    return r
+
+
 def _1_over(x):
     '''(INTERNAL) Return reciprocal C{1 / B{x}} avoiding C{ZeroDivisionError} exceptions.
     '''
@@ -168,6 +186,7 @@ _0_1     = _float(   0.1)     # PYCHOK expected
 _0_125   = _float(   0.125)   # PYCHOK expected
 _0_25    = _float(   0.25)    # PYCHOK expected
 _0_5     = _float(   0.5)     # PYCHOK expected
+_0_75    = _float(   0.75)     # PYCHOK expected
 _1_0     = _float(   1)       # PYCHOK expected
 _1_0_1T  = _1_0,              # PYCHOK 1-tuple
 _1_5     = _float(   1.5)     # PYCHOK expected
@@ -218,40 +237,35 @@ try:
 #   RADIX    =  Int(  RADIX   =_f_i.radix)     # PYTHON system's float base
     del _f_i
 except ImportError:  # PYCHOK no cover
-    DIG      =  Int(  DIG     =15)          # PYCHOK system's 64-bit float decimal digits
+    DIG      =  Int(  DIG     =15)           # PYCHOK system's 64-bit float decimal digits
     EPS      = _Float(EPS     =2.220446049250313e-16)  # PYCHOK EPSilon 2**-52, M{EPS +/- 1 != 1}
-    MANT_DIG =  Int(  MANT_DIG=53)          # PYCHOK float mantissa bits ≈ 53 (C{int})
+    MANT_DIG =  Int(  MANT_DIG=53)           # PYCHOK float mantissa bits ≈ 53 (C{int})
     MAX      = _Float(MAX     =pow(_2_0,  1023) * (_2_0 - EPS))  # PYCHOK ≈ 10**308
-    MAX_EXP  =  Int(  MAX_ESP =_log2(MAX))  # 308 base 10
+    MAX_EXP  =  Int(  MAX_ESP =_log2(MAX))   # 308 base 10
     MIN      = _Float(MIN     =pow(_2_0, -1021))  # PYCHOK ≈ 10**-308
-    MIN_EXP  =  Int(MIN_EXP   =_log2(MIN))  # -307 base 10
-#   RADIX    =  Int(Radix     =2)           # base
-
-EPS0     = _Float( EPS0  = EPS**2)          # PYCHOK near-/non-zero comparison 4.930381e-32, or EPS or EPS_2
-EPS02    = _Float( EPS02 = EPS**4)          # PYCHOK near-zero-squared comparison 2.430865e-63
-EPS_2    = _Float( EPS_2 = EPS / _2_0)      # PYCHOK ≈ 1.110223024625e-16
-EPS1     = _Float( EPS1  =_1_0 - EPS)       # PYCHOK ≈ 0.9999999999999998
-EPS2     = _Float( EPS2  = EPS * _2_0)      # PYCHOK ≈ 4.440892098501e-16
-EPS4     = _Float( EPS4  = EPS * _4_0)      # PYCHOK ≈ 8.881784197001e-16
-# _1EPS  = _Float(_1EPS  =_1_0 + EPS)       # PYCHOK ≈ 1.0000000000000002
-_1_EPS   = _Float(_1_EPS =_1_0 / EPS)       # PYCHOK = 4503599627370496.0
-# _2_EPS = _Float(_2_EPS =_2_0 / EPS)       # PYCHOK = 9007199254740992.0
-_EPS2e4  = _Float(_EPS2e4= EPS2 * 1.e4)     # PYCHOK ≈ 4.440892098501e-12
-_EPS4e8  = _Float(_EPS4e8= EPS4 * 1.e8)     # PYCHOK ≈ 8.881784197001e-08
-_EPSjam  = _Float(_EPSjam= pow(EPS, 0.75))  # PYCHOK = 1.818989403546e-12
-_EPSmin  = _Float(_EPSmin= sqrt(MIN))       # PYCHOK = 1.49166814624e-154
-_EPSqrt  = _Float(_EPSqrt= sqrt(EPS))       # PYCHOK = 1.49011611938e5-08
-_EPStol  = _Float(_EPStol=_EPSqrt * _0_1)   # PYCHOK = 1.49011611938e5-09 == sqrt(EPS * _0_01)
+    MIN_EXP  =  Int(MIN_EXP   =_log2(MIN))   # -307 base 10
+#   RADIX    =  Int(Radix     =2)            # base
+
+EPS0     = _Float( EPS0  = EPS**2)           # PYCHOK near-/non-zero comparison 4.930381e-32, or EPS or EPS_2
+EPS02    = _Float( EPS02 = EPS**4)           # PYCHOK near-zero-squared comparison 2.430865e-63
+EPS_2    = _Float( EPS_2 = EPS / _2_0)       # PYCHOK ≈ 1.110223024625e-16
+EPS1     = _Float( EPS1  =_1_0 - EPS)        # PYCHOK ≈ 0.9999999999999998
+EPS2     = _Float( EPS2  = EPS * _2_0)       # PYCHOK ≈ 4.440892098501e-16
+EPS4     = _Float( EPS4  = EPS * _4_0)       # PYCHOK ≈ 8.881784197001e-16
+# _1EPS  = _Float(_1EPS  =_1_0 + EPS)        # PYCHOK ≈ 1.0000000000000002
+_1_EPS   = _Float(_1_EPS =_1_0 / EPS)        # PYCHOK = 4503599627370496.0
+# _2_EPS = _Float(_2_EPS =_2_0 / EPS)        # PYCHOK = 9007199254740992.0
+_EPS2e4  = _Float(_EPS2e4= EPS2 * 1.e4)      # PYCHOK ≈ 4.440892098501e-12
+_EPS4e8  = _Float(_EPS4e8= EPS4 * 1.e8)      # PYCHOK ≈ 8.881784197001e-08
+_EPSjam  = _Float(_EPSjam= pow(EPS, _0_75))  # PYCHOK = 1.818989403546e-12
+_EPSmin  = _Float(_EPSmin= sqrt(MIN))        # PYCHOK = 1.49166814624e-154
+_EPSqrt  = _Float(_EPSqrt= sqrt(EPS))        # PYCHOK = 1.49011611938e5-08
+_EPStol  = _Float(_EPStol=_EPSqrt * _0_1)    # PYCHOK = 1.49011611938e5-09 == sqrt(EPS * _0_01)
 
 _89_999  = _Float(_89_999=_90_0 * EPS1)  # just below 90.0
 # <https://Numbers.Computation.Free.FR/Constants/Miscellaneous/digits.html>
-_1__90   = _Float(_1__90 =_1_0 / _90_0)  # PYCHOK = 0.011_111_111_111_111_111_111_111_111_111_111_111_111_111_111_11111
+# _1__90 = _Float(_1__90 =_1_0 / _90_0)  # PYCHOK = 0.011_111_111_111_111_111_111_111_111_111_111_111_111_111_111_11111
 _2__PI   = _Float(_2__PI =_2_0 / _pi)    # PYCHOK = 0.636_619_772_367_581_343_075_535_053_490_057_448_137_838_582_96182
-
-_1_16th  = _Float(_1_16th=_1_0 / _16_0)  # PYCHOK in .ellipsoids, .karney
-_1_3rd   = _Float(_1_3rd =_1_0 /  _3_0)  # PYCHOK in .fmath
-_1_6th   = _Float(_1_6th =_1_0 /  _6_0)  # PYCHOK in .fmath
-
 _K0_UTM  = _Float(_K0_UTM = 0.9996)  # PYCHOK in .etm, .ktm, .utm, UTM scale at central meridian
 # sqrt(2) <https://WikiPedia.org/wiki/Square_root_of_2>
 # 1.414213562373095_048_801_688_724_209_698_078_569_671_875_376_948_073_176_679_737_99
@@ -273,6 +287,7 @@ PI3_2 = _Float(PI3_2=_pi * _1_5)  # PYCHOK PI and a half, M{PI * 3 / 2}
 PI_3  = _Float(PI_3 =_pi / _3_0)  # PYCHOK One third PI, M{PI / 3}
 PI4   = _Float(PI4  =_pi * _4_0)  # PYCHOK Four PI, M{PI * 4}
 PI_4  = _Float(PI_4 =_pi / _4_0)  # PYCHOK Quarter PI, M{PI / 4}
+PI_6  = _Float(PI_6 =_pi / _6_0)  # PYCHOK One sixth PI, M{PI / 6}
 
 R_MA  = _Radius(R_MA=6378137.0)       # PYCHOK equatorial earth radius (C{meter}), WGS84, EPSG:3785
 R_MB  = _Radius(R_MB=6356752.3)       # PYCHOK polar earth radius (C{meter}), WGS84, EPSG:3785
@@ -417,6 +432,17 @@ def isnear90(x, eps90=EPS0):
     return bool(eps90 > (x - _90_0) > -eps90)
 
 
+def isneg(x):
+    '''Check for negative C{x}, including L{NEG0}.
+
+       @arg x: Value (C{scalar}).
+
+       @return: C{True} if C{B{x} < 0 or NEG0},
+                C{False} otherwise.
+    '''
+    return signBit(x)
+
+
 def isneg0(x):
     '''Check for L{NEG0}, negative C{0.0}.
 

pygeodesy/deprecated/functions.py

@@ -4,6 +4,7 @@
 u'''DEPRECATED functions kept for backward compatibility.
 '''
 
+# from pygeodesy.basics import copysign0  # _MODS_
 from pygeodesy.constants import EPS, R_M, float0_
 from pygeodesy.deprecated.classes import ClipCS3Tuple, TriAngle4Tuple, _TriAngle5Tuple
 from pygeodesy.interns import NN, _area_, _COMMASPACE_, _negative_, \

pygeodesy/dms.py

@@ -63,7 +63,7 @@ U{Vector-based geodesy<https://www.Movable-Type.co.UK/scripts/latlong-vectors.ht
 
 from pygeodesy.basics import copysign0, isLatLon, isodd, issequence, isstr, \
                              neg as _neg,  typename  # in .ups
-from pygeodesy.constants import _umod_360, _0_0, _0_5, _60_0, _360_0, _3600_0
+from pygeodesy.constants import _0_0, _0_5, _60_0, _360_0, _3600_0, _umod_360
 from pygeodesy.errors import ParseError, RangeError, _TypeError, _ValueError, \
                             _parseX, rangerrors, _xError, _xkwds,  _envPYGEODESY
 # from pygeodesy.internals import _envPYGEODESY, typename  # from .errors
@@ -86,7 +86,7 @@ except ImportError:  # Python 3+
     from string import ascii_letters as _LETTERS
 
 __all__ = _ALL_LAZY.dms
-__version__ = '25.04.14'
+__version__ = '25.08.31'
 
 _beyond_      = 'beyond'
 _deg_min_     = 'deg+min'

pygeodesy/ecef.py

@@ -62,8 +62,8 @@ plane} as opposed to I{geocentric} (ECEF) ones.
 from pygeodesy.basics import copysign0, _isin, isscalar, issubclassof, neg, map1, \
                             _xinstanceof, _xsubclassof,  typename  # _args_kwds_names
 from pygeodesy.constants import EPS, EPS0, EPS02, EPS1, INT0, PI, PI_2, _0_0, \
-                               _0_5, _1_0, _1_0_1T, _2_0, _3_0, _4_0, _6_0, \
-                               _90_0, _copysign_1_0,  isnon0  # PYCHOK used!
+                               _0_5, _1_0, _1_0_1T, _2_0, _3_0, _4_0, _6_0, _90_0, \
+                               _copysign_1_0,  isnon0  # PYCHOK used!
 from pygeodesy.datums import _ellipsoidal_datum, _WGS84,  a_f2Tuple, _EWGS84
 from pygeodesy.ecefLocals import _EcefLocal
 # from pygeodesy.ellipsoids import a_f2Tuple, _EWGS84  # from .datums
@@ -91,7 +91,7 @@ from pygeodesy.utily import atan1, atan1d, atan2, atan2d, degrees90, degrees180,
 from math import cos, degrees, fabs, radians, sqrt
 
 __all__ = _ALL_LAZY.ecef
-__version__ = '25.08.25'
+__version__ = '25.08.31'
 
 _Ecef_    = 'Ecef'
 _prolate_ = 'prolate'

pygeodesy/ellipsoidalExact.py

@@ -11,20 +11,20 @@ L{GeodesicExact}, L{GeodesicAreaExact} and L{GeodesicLineExact}.
 
 # from pygeodesy.datums import _WGS84  # from .ellipsoidalBase
 from pygeodesy.ellipsoidalBase import CartesianEllipsoidalBase, \
-                                     _nearestOn,  _WGS84
+                                     _nearestOn,  Property_RO, _WGS84
 from pygeodesy.ellipsoidalBaseDI import LatLonEllipsoidalBaseDI, \
                                        _intersection3, _intersections2, \
                                        _TOL_M, intersecant2
 # from pygeodesy.errors import _xkwds  # from .karney
-from pygeodesy.karney import _polygon,  fabs, Property_RO, _xkwds
+from pygeodesy.karney import _polygon,  fabs, _xkwds
 from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS, _ALL_OTHER
 from pygeodesy.points import _areaError, ispolar  # PYCHOK exported
-# from pygeodesy.props import Property_RO  # from .karney
+# from pygeodesy.props import Property_RO  # from .ellipsoidalBase
 
 # from math import fabs  # from .karney
 
 __all__ = _ALL_LAZY.ellipsoidalExact
-__version__ = '25.05.28'
+__version__ = '25.08.28'
 
 
 class Cartesian(CartesianEllipsoidalBase):

pygeodesy/ellipsoidalGeodSolve.py

@@ -12,11 +12,11 @@ L{geodsolve}, a wrapper invoking I{Karney}'s U{GeodSolve
 
 # from pygeodesy.datums import _WGS84  # from .ellipsoidalBase
 from pygeodesy.ellipsoidalBase import CartesianEllipsoidalBase, \
-                                     _nearestOn, _WGS84
+                                     _nearestOn,  Property_RO, _WGS84
 from pygeodesy.ellipsoidalBaseDI import LatLonEllipsoidalBaseDI, _TOL_M, \
                                        _intersection3, _intersections2
 # from pygeodesy.errors import _xkwds  # from .karney
-from pygeodesy.karney import fabs, _polygon,  Property_RO, _xkwds
+from pygeodesy.karney import fabs, _polygon,  _xkwds
 from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS, _ALL_OTHER
 from pygeodesy.points import _areaError, ispolar  # PYCHOK exported
 # from pygeodesy.props import Property_RO  # from .karney
@@ -24,7 +24,7 @@ from pygeodesy.points import _areaError, ispolar  # PYCHOK exported
 # from math import fabs  # from .karney
 
 __all__ = _ALL_LAZY.ellipsoidalGeodSolve
-__version__ = '25.05.28'
+__version__ = '25.08.28'
 
 
 class Cartesian(CartesianEllipsoidalBase):