pygeodesy 25.9.9__py2.py3-none-any.whl → 25.11.5__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.7 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.1),
+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,27 +137,23 @@ 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.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.
+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.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
+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.7 (and U{geographiclib<https://PyPI.org/project/geographiclib>} 2.1) on
@@ -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-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,
+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>}.
@@ -605,7 +605,7 @@ else:
 
 from pygeodesy.internals import _version2,  _DOT_  # noqa: E402
 # from pygeodesy.interns import _DOT_  # from .internals
-__version__ = '25.09.09'
+__version__ = '25.11.05'
 # see setup.py for similar logic
 version     = _DOT_(*_version2(__version__, n=3))
 

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/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.09.04'
+__version__ = '25.10.30'
 
 _0EPS  =  EPS  # near-zero, positive
 _EPS0  = -EPS  # near-zero, negative
@@ -1434,17 +1434,14 @@ class _CompositeGH(_CompositeBase):
 
     def _clip(self, corners, s_entry, c_entry, Clas=None,
                            **closed_inull_raiser_xtend_eps):
-        # Clip this polygon with another one, C{corners}.
-
-        # Core of Greiner/Hormann's algorithm, enhanced U{Correia's
-        # <https://GitHub.com/helderco/univ-polyclip>} implementation***
-        # and extended to optionally handle so-called "degenerate cases"
+        # Clip this polygon with another one, C{corners}, using
+        # Greiner/Hormann's algorithm extended to optionally
+        # handle some so-called "degenerate cases"
         S = self
         C = self._class(corners, closed_inull_raiser_xtend_eps,
                                  raiser=False, xtend=False)
         bt = C._bottom_top_eps2
         lr = C._left_right_eps2
-
         # 1. find intersections
         for s1, s2, Sc in S._edges3(**closed_inull_raiser_xtend_eps):
             if not (_outside(s1.x, s2.x, *lr) or
@@ -1507,13 +1504,13 @@ class _CompositeGH(_CompositeBase):
 
     @property
     def xtend(self):
-        '''Get the option to handle I{degenerate cases} (C{bool}).
+        '''Get the option to handle some I{degenerate cases} (C{bool}).
         '''
         return self._xtend
 
     @xtend.setter  # PYCHOK setter!
     def xtend(self, xtend):
-        '''Set the option to handle I{degenerate cases} (C{bool}).
+        '''Set the option to handle some I{degenerate cases} (C{bool}).
         '''
         self._xtend = bool(xtend)
 
@@ -1825,9 +1822,7 @@ class BooleanFHP(_CompositeFHP, _BooleanBase):
 class BooleanGH(_CompositeGH, _BooleanBase):
     '''I{Composite} class providing I{boolean} operations between two
        I{composites} using the U{Greiner-Hormann<http://www.Inf.USI.CH/
-       hormann/papers/Greiner.1998.ECO.pdf>} algorithm and U{Correia
-       <https://GitHub.com/helderco/univ-polyclip>}'s implementation,
-       modified and extended.
+       hormann/papers/Greiner.1998.ECO.pdf>} algorithm, extended.
 
        The supported operations between (composite) polygon A and B are:
 
@@ -2003,23 +1998,3 @@ __all__ += _ALL_DOCS(_BooleanBase, _Clip,
 # OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,
 # ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
 # OTHER DEALINGS IN THE SOFTWARE.
-
-# ***) GNU GPL 3
-#
-# Copyright (C) 2011-2012 Helder Correia <Helder.MC@Gmail.com>
-#
-# This program is free software: you can redistribute it and/or
-# modify it under the terms of the GNU General Public License as
-# published by the Free Software Foundation, either version 3 of
-# the License, or any later version.
-#
-# This program is distributed in the hope that it will be useful,
-# but WITHOUT ANY WARRANTY; without even the implied warranty of
-# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
-# GNU General Public License for more details.
-#
-# You should have received a copy of the GNU General Public License
-# along with this program.  If not, see <http://www.GNU.org/licenses/>.
-#
-# You should have received the README file along with this program.
-# If not, see <https://GitHub.com/helderco/univ-polyclip>.

pygeodesy/clipy.py

@@ -31,7 +31,7 @@ from pygeodesy.units import Bool, FIx, HeightX, Lat, Lon, Number_
 # from math import fabs  # from .fmath
 
 __all__ = _ALL_LAZY.clipy
-__version__ = '25.05.12'
+__version__ = '25.10.30'
 
 _fj_       = 'fj'
 _original_ = 'original'
@@ -237,9 +237,8 @@ def clipCS4(points, lowerleft, upperight, closed=False, inull=False):
 
 
 class ClipFHP4Tuple(_NamedTuple):
-    '''4-Tuple C{(lat, lon, height, clipid)} for each point of the
-       L{clipFHP4} result with the C{lat}-, C{lon}gitude, C{height}
-       and C{clipid} of the polygon or clip.
+    '''4-Tuple C{(lat, lon, height, clipid)} for each point of the L{clipFHP4} result
+       with the C{lat}-, C{lon}gitude, C{height} and C{clipid} of the polygon or clip.
 
        @note: The C{height} is a L{HeightX} instance if this point is
               an intersection, otherwise a L{Height} or C{int(0)}.
@@ -275,14 +274,13 @@ def clipFHP4(points, corners, closed=False, inull=False, raiser=False, eps=EPS):
                    as the B{C{points}} and B{C{corners}} coordinates).
 
        @return: Yield a L{ClipFHP4Tuple}C{(lat, lon, height, clipid)} for each
-                clipped point.  The result may consist of several clips, each
-                a (closed) polygon with a unique C{clipid}.
+                clipped point.  The result may consist of several I{clips}, each
+                a (closed) polygon with a unique C{clipid} identifier.
 
        @raise ClipError: Insufficient B{C{points}} or B{C{corners}} or an open clip.
 
        @see: U{Forster, Hormann and Popa<https://www.ScienceDirect.com/science/
-             article/pii/S259014861930007X>}, class L{BooleanFHP} and function
-             L{clipGH4}.
+             article/pii/S259014861930007X>} and class L{BooleanFHP}.
     '''
     P = _MODS.booleans._CompositeFHP(points, kind=_points_, raiser=raiser,
                                                    eps=eps, name__=clipFHP4)
@@ -292,12 +290,11 @@ def clipFHP4(points, corners, closed=False, inull=False, raiser=False, eps=EPS):
 
 
 class ClipGH4Tuple(ClipFHP4Tuple):
-    '''4-Tuple C{(lat, lon, height, clipid)} for each point of the
-       L{clipGH4} result with the C{lat}-, C{lon}gitude, C{height}
-       and C{clipid} of the polygon or clip.
+    '''4-Tuple C{(lat, lon, height, clipid)} for each point of the L{clipGH4} result
+       with the C{lat}-, C{lon}gitude, C{height} and C{clipid} of the polygon or clip.
 
-       @note: The C{height} is a L{HeightX} instance if this is
-              an intersection, otherwise a L{Height} or C{int(0)}.
+       @note: The C{height} is a L{HeightX} instance if this is an intersection,
+              otherwise a L{Height} or C{int(0)}.
     '''
     _Names_ = ClipFHP4Tuple._Names_
     _Units_ = ClipFHP4Tuple._Units_
@@ -305,8 +302,7 @@ class ClipGH4Tuple(ClipFHP4Tuple):
 
 def clipGH4(points, corners, closed=False, inull=False, raiser=True, xtend=False, eps=EPS):
     '''Clip one or more polygons against a clip region or box using the U{Greiner-Hormann
-       <http://www.Inf.USI.CH/hormann/papers/Greiner.1998.ECO.pdf>} algorithm, U{Correia
-       <https://GitHub.com/helderco/univ-polyclip>}'s implementation modified and extended.
+       <http://www.Inf.USI.CH/hormann/papers/Greiner.1998.ECO.pdf>} algorithm, extended.
 
        @arg points: The polygon points and clips (C{LatLon}[]).
        @arg corners: Three or more points defining the clip regions (C{LatLon}[])
@@ -320,19 +316,17 @@ def clipGH4(points, corners, closed=False, inull=False, raiser=True, xtend=False
                    as the B{C{points}} and B{C{corners}} coordinates).
 
        @return: Yield a L{ClipGH4Tuple}C{(lat, lon, height, clipid)} for each
-                clipped point.  The result may consist of several clips, each
-                a (closed) polygon with a unique C{clipid}.
+                clipped point.  The result may consist of several I{clips}, each
+                a (closed) polygon with a unique C{clipid} identifier.
 
        @raise ClipError: Insufficient B{C{points}} or B{C{corners}}, an open clip,
                          a I{degenerate case} or I{unhandled} intersection.
 
        @note: To handle I{degenerate cases} like C{point-edge} and C{point-point}
-              intersections, use function L{clipFHP4}.
+              intersections I{properly}, use function L{clipFHP4}.
 
-       @see: U{Greiner-Hormann<https://WikiPedia.org/wiki/Greiner–Hormann_clipping_algorithm>},
-             U{Ionel Daniel Stroe<https://Davis.WPI.edu/~matt/courses/clipping/>}, I{Correia}'s
-             U{univ-polyclip<https://GitHub.com/helderco/univ-polyclip>}, class L{BooleanGH}
-             and function L{clipFHP4}.
+       @see: U{Greiner-Hormann<https://WikiPedia.org/wiki/Greiner–Hormann_clipping_algorithm>}
+             and class L{BooleanGH}.
     '''
     S = _MODS.booleans._CompositeGH(points, raiser=raiser, xtend=xtend, eps=eps,
                                             name__=clipGH4, kind=_points_)

pygeodesy/constants.py

@@ -11,7 +11,7 @@ L{pygeodesy.isnon0} and L{pygeodesy.remainder}.
 from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import _copysign, isbool, iscomplex, isint, signBit
-from pygeodesy.errors import _xError, _xError2, _xkwds_get1, _xkwds_item2
+from pygeodesy.errors import _ValueError, _xError, _xkwds_get1, _xkwds_item2
 # from pygeodesy.fsums import _isFsum_2Tuple  # _MODS
 from pygeodesy.internals import _0_0, _100_0, typename
 from pygeodesy.interns import _DMAIN_, _INF_, _NAN_
@@ -26,7 +26,7 @@ except ImportError:  # Python 2-
     _inf, _nan = float(_INF_), float(_NAN_)
 
 __all__ = _ALL_LAZY.constants
-__version__ = '25.09.09'
+__version__ = '25.10.15'
 
 
 def _copysign_0_0(y):
@@ -47,7 +47,7 @@ def _copysignINF(y):
     return NINF if signBit(y) else INF
 
 
-def _flipsign(x, y):
+def _flipsign(x, y=-1):
     '''(INTERNAL) Negate C{x} for negative C{y}.
     '''
     return (-x) if signBit(y) else x
@@ -60,65 +60,78 @@ def _Float(**name_arg):
     return Float(_float(arg), name=n)
 
 
-def _Radius(**name_arg):
-    '''(INTERNAL) New named, cached C{Radius}.
+def _float(x):  # in .datums, .ellipsoids, ...
+    '''(INTERNAL) Cache initial C{float}s.
     '''
-    n, arg = _xkwds_item2(name_arg)
-    return Radius(_float(arg), name=n)
+    f = float(x)
+    return _floats.setdefault(f, f)  # PYCHOK del _floats
 
 
-def float_(*fs, **sets):  # sets=False
-    '''Get scalars as C{float} or I{intern}'ed C{float}.
+def float_(x, sets=False):
+    '''Get scalar as C{float} or I{intern}'ed C{float}.
 
-       @arg fs: One more values (C{scalar}), all positional.
-       @kwarg sets: Use C{B{sets}=True} to C{intern} each
-                    B{C{fs}}, otherwise don't C{intern}.
+       @arg x: The scalar (C{scalar}).
+       @kwarg sets: Use C{True} to C{intern} the B{C{f}},
+                    otherwise don't (C{bool}).
 
-       @return: A single C{float} if only one B{C{fs}} is
-                given, otherwise a tuple of C{float}s.
+       @return: A C{float}.
 
-       @raise TypeError: Some B{C{fs}} is not C{scalar}.
+       @raise ValueError: Invalid B{C{x}}.
     '''
-    fl = []
-    _a =  fl.append
-    _f = _floats.setdefault if _xkwds_get1(sets, sets=False) else \
-         _floats.get
     try:
-        for i, f in enumerate(fs):
-            f = float(f)
-            _a(_f(f, f))
-    except Exception as x:
-        E, t = _xError2(x)
-        fs_i =  Fmt.SQUARE(fs=i)
-        raise E(fs_i, f, txt=t)
-    return fl[0] if len(fl) == 1 else tuple(fl)
-
-
-def _float(f):  # in .datums, .ellipsoids, ...
-    '''(INTERNAL) Cache initial C{float}s.
-    '''
-    f = float(f)
-    return _floats.setdefault(f, f)  # PYCHOK del _floats
+        f = float(x)
+        if f:
+            _f = _floats.setdefault if sets else _floats.get
+            f  = _f(f, f)
+        else:
+            f  = _N_0_0 if signBit(f) else _0_0
+    except Exception as X:
+        raise _ValueError(x=x, cause=X)
+    return f
 
 
 def float0_(*xs):
     '''Yield C{B{x}s} as a non-NEG0 C{float}.
     '''
     for x in xs:
-        yield float(x) if x else _0_0
+        yield float(x) or _0_0  # if x else _0_0
 
 
-def _float0(f):  # in .auxilats.auxily, .resections, .vector3dBase, ...
+def _float0(x):  # in .auxilats.auxily, .resections, .vector3dBase, ...
     '''(INTERNAL) Return C{float(B{f})} or C{INT0}.
     '''
-    if f:
-        f =  float(f)
+    if x:
+        f =  float(x)
         f = _floats.get(f, f)
-    elif f is not INT0:
-        f =  float(f) or _0_0  # force None, NN error
+    elif x is INT0:
+        f =  x
+    else:
+        f =  float(x) or _0_0  # force None, NN error
     return f
 
 
+def floats_(*xs, **sets):  # sets=False
+    '''Yield each scalar as C{float} or I{intern}'ed C{float}.
+
+       @arg xs: One more values (C{scalar}), all positional.
+       @kwarg sets: Use C{B{sets}=True} to C{intern} each
+                    B{C{fs}}, otherwise don't C{intern}.
+
+       @raise ValueError: Some B{C{xs}} is not C{scalar}.
+    '''
+    if sets:
+        sets = _xkwds_get1(sets, sets=False)
+    _f = _floats.setdefault if sets else _floats.get
+    try:
+        for i, x in enumerate(xs):
+            f = float(x)
+            yield _f(f, f) if f else \
+                  (_N_0_0 if signBit(f) else _0_0)  # preserve NEG0
+    except Exception as X:
+        xs_i = Fmt.SQUARE(xs=i)
+        raise _ValueError(xs_i, x, cause=X)
+
+
 def _floatuple(*fs):
     '''(INTERNAL) Cache a tuple of C{float}s.
     '''
@@ -172,6 +185,13 @@ def _1_over(x):
         return  NINF if isneg0(x) else INF
 
 
+def _Radius(**name_arg):
+    '''(INTERNAL) New named, cached C{Radius}.
+    '''
+    n, arg = _xkwds_item2(name_arg)
+    return Radius(_float(arg), name=n)
+
+
 _floats  = {}     # PYCHOK floats cache, in .__main__
 # _float = float  # PYCHOK expected
 # del _floats     # XXX zap floats cache never
@@ -186,6 +206,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
@@ -236,40 +257,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
 _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

pygeodesy/deprecated/__init__.py

@@ -27,7 +27,7 @@ __all__ = (_ALL_DEPRECATED.deprecated_bases +
            _ALL_DEPRECATED.deprecated_classes +
            _ALL_DEPRECATED.deprecated_consterns +
            _ALL_DEPRECATED.deprecated_functions)
-__version__ = '25.05.12'
+__version__ = '25.10.25'
 
 if _unLazy0:
     from pygeodesy.deprecated import bases, datum, nvector, rhumbBase, \

pygeodesy/deprecated/classes.py

@@ -18,11 +18,12 @@ from pygeodesy.namedTuples import Forward4Tuple as _Forward4Tuple, \
                                   UtmUps5Tuple  as _UtmUps5Tuple,  _NamedTuple
 from pygeodesy.props import deprecated_class, deprecated_method
 from pygeodesy.resections import TriAngle5Tuple as _TriAngle5Tuple
+from pygeodesy.triaxials import ConformalSphere, ConformalTriaxial, Conformal2Tuple
 from pygeodesy.trf import TRFXform7Tuple as _TRFXform7Tuple
 from pygeodesy.units import Bearing, Int, Lamd, Lat, Lon, Meter, Phid
 
 __all__ = _ALL_DEPRECATED.deprecated_classes
-__version__ = '25.04.11'
+__version__ = '25.10.25'
 
 
 class _Deprecated_NamedTuple(_NamedTuple):
@@ -235,6 +236,28 @@ def HeightIDWcosineForsytheAndoyerLambert(knots, **kwds):  # PYCHOK no cover
     return HeightIDWcosineForsytheAndoyerLambert(knots, **kwds)
 
 
+class JacobiConformal(ConformalTriaxial):
+    '''DEPRECATED on 2025.10.25, use class L{ConformalTriaxial}.'''
+    def __init__(self, *args, **kwds):  # PYCHOK no cover
+        deprecated_class(self.__class__)
+        ConformalTriaxial.__init__(self, *args, **kwds)
+
+
+class JacobiConformalSpherical(ConformalSphere):
+    '''DEPRECATED on 2025.10.25, use class L{ConformalSphere}.'''
+    def __init__(self, *args, **kwds):  # PYCHOK no cover
+        deprecated_class(self.__class__)
+        ConformalSphere.__init__(self, *args, **kwds)
+
+
+class Jacobi2Tuple(Conformal2Tuple):
+    '''DEPRECATED on 25.10.25, use class L{Conformal2Tuple}.
+    '''
+    def __new__(cls, *args, **kwds):
+        deprecated_class(cls)
+        return Conformal2Tuple.__new__(cls, *args, **kwds)
+
+
 class Lam_(Lamd):
     '''DEPRECATED on 2024.06.15, use class L{Lamd}.'''
     def __init__(self, *args, **kwds):  # PYCHOK no cover