pygeodesy 25.7.25__py2.py3-none-any.whl → 25.9.9__py2.py3-none-any.whl

pygeodesy/__init__.py

@@ -125,7 +125,7 @@ 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.13.7 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.1),
 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,
@@ -143,20 +143,20 @@ 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.5 Sequoia.
+macOS 15.6.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.4, 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.7, 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.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.13.7, 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.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-5, 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-4 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-5 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
@@ -190,7 +190,7 @@ All Python source code has been statically U{checked<https://GitHub.com/ActiveSt
 546532_PyChecker_postprocessor>} with U{Ruff<https://GitHub.com/astral-sh/ruff>} using Python 3.13.5 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.5 Sequoia.
+2.7.18, both in 64-bit on macOS 15.6.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>}.
@@ -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.07.25'
+__version__ = '25.09.09'
 # 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/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/cartesianBase.py

@@ -27,11 +27,11 @@ from pygeodesy.interns import _COMMASPACE_, _datum_, _no_, _phi_
 from pygeodesy.interns import _ellipsoidal_, _spherical_  # PYCHOK used!
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS
 from pygeodesy.named import _name2__, _Pass
-from pygeodesy.namedTuples import LatLon4Tuple, _NamedTupleTo , Vector3Tuple, \
+from pygeodesy.namedTuples import LatLon4Tuple, _NamedTupleTo, Vector3Tuple, \
                                   Vector4Tuple
-# from pygeodesy.nvectorBase import _N_vector  # _MODS
-from pygeodesy.props import deprecated_method, Property, Property_RO, property_doc_, \
-                            property_RO, _update_all
+# from pygeodesy.nvectorBase import _N_Vector  # _MODS
+from pygeodesy.props import deprecated_method, Property, Property_RO, \
+                            property_doc_, property_RO, _update_all
 # from pygeodesy import resections as _resections  # _MODS.into
 # from pygeodesy.streprs import Fmt  # from .fsums
 # from pygeodesy.triaxials import Triaxial_  # _MODS
@@ -39,12 +39,11 @@ from pygeodesy.units import Degrees, Height, _heigHt, _isMeter, Meter, Radians
 from pygeodesy.utily import acos1, atan2, sincos2d, sincos2_,  degrees, radians
 from pygeodesy.vector3d import Vector3d, _xyzhdlln4
 # from pygeodesy.vector3dBase import _xyz3  # _MODS
-# from pygeodesy import ltp  # _MODS
 
 # from math import degrees, fabs, radians, sqrt  # from .fmath, .utily
 
 __all__ = _ALL_LAZY.cartesianBase
-__version__ = '25.05.07'
+__version__ = '25.08.24'
 
 _r_         = 'r'
 _resections = _MODS.into(resections=__name__)
@@ -172,12 +171,13 @@ class CartesianBase(Vector3d, _EcefLocal):
            @raise TypeError: The B{C{datum}} is not a L{Datum}.
         '''
         d = _spherical_datum(datum, name=self.name)
-        if self._datum:  # is not None
-            if d.isEllipsoidal and not self._datum.isEllipsoidal:
+        D =  self._datum
+        if D:  # is not None
+            if d.isEllipsoidal and not D.isEllipsoidal:
                 raise _IsnotError(_ellipsoidal_, datum=datum)
-            elif d.isSpherical and not self._datum.isSpherical:
+            elif d.isSpherical and not D.isSpherical:
                 raise _IsnotError(_spherical_, datum=datum)
-        if self._datum != d:
+        if D != d:  # or (D and D.name != d.name)
             _update_all(self)
             self._datum = d
 
@@ -302,8 +302,6 @@ class CartesianBase(Vector3d, _EcefLocal):
            @note: Include keyword argument C{B{datum}=None} if class B{C{Cartesian}}
                   does not accept a B{C{datum}} keyword agument.
 
-           @raise TriaxialError: No convergence in triaxial root finding.
-
            @raise TypeError: Invalid or undefined B{C{earth}} or C{datum}.
         '''
         n = typename(self.height3)
@@ -410,13 +408,12 @@ class CartesianBase(Vector3d, _EcefLocal):
 
     @Property_RO
     def _N_vector(self):
-        '''(INTERNAL) Get the (C{nvectorBase._N_vector_}).
+        '''(INTERNAL) Get the (C{nvectorBase._N_Vector}).
         '''
-        _N = _MODS.nvectorBase._N_vector_
-        x, y, z, h = self._n_xyzh4(self.datum)
-        return _N(x, y, z, h=h, name=self.name)
+        _N = _MODS.nvectorBase._N_Vector
+        return self._n_xyzh4(self.datum, Vector=_N)
 
-    def _n_xyzh4(self, datum):
+    def _n_xyzh4(self, datum, Vector=Vector4Tuple):
         '''(INTERNAL) Get the n-vector components as L{Vector4Tuple}.
         '''
         def _ErrorEPS0(x):
@@ -434,7 +431,11 @@ class CartesianBase(Vector3d, _EcefLocal):
         p = hypot2(x, y) * E.a2_
         q = z**2 * E.e21 * E.a2_
         r = fsumf_(p, q, -E.e4) / _6_0
+        if isnear0(r):
+            raise _ErrorEPS0(r)
         s = (p * q * E.e4) / (_4_0 * r**3)
+        if s < 0:
+            raise _ErrorEPS0(s)
         t = cbrt(fsumf_(_1_0, s, sqrt(s * (_2_0 + s))))
         if isnear0(t):
             raise _ErrorEPS0(t)
@@ -458,7 +459,7 @@ class CartesianBase(Vector3d, _EcefLocal):
             raise _ErrorEPS0(t)
         h = fsumf_(k, E.e2, _N_1_0) / k * t
         s = e / t  # == e * tmp
-        return Vector4Tuple(x * s, y * s, z / t, h, name=self.name)
+        return Vector(x * s, y * s, z / t, h, name=self.name)
 
     @Property_RO
     def philam(self):
@@ -980,8 +981,8 @@ def xyz2rtp_(x_xyz, y=0, z=0, **name):
        @kwarg name: Optional C{B{name}=NN} (C{str}).
 
        @return: L{RadiusThetaPhi3Tuple}C{(r, theta, phi)} with radial distance C{r} (C{meter},
-                same units as C{x}, C{y} and C{z}), inclination C{theta} (with respect to the
-                positive z-axis) and azimuthal angle C{phi}, both in L{Radians}.
+                same units as C{x}, C{y} and C{z}), inclination C{theta} (polar angle with
+                respect to the positive z-axis) and azimuthal angle C{phi}, both in L{Radians}.
 
        @see: U{Physics convention<https://WikiPedia.org/wiki/Spherical_coordinate_system>}
              (ISO 80000-2:2019), class L{RadiusThetaPhi3Tuple} and function L{xyz2rtp}.

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.09'
 
 
 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.
     '''
@@ -245,7 +263,7 @@ _EPStol  = _Float(_EPStol=_EPSqrt * _0_1)   # PYCHOK = 1.49011611938e5-09 == sqr
 
 _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
@@ -273,6 +291,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 +436,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'