pygeodesy 24.9.9__py2.py3-none-any.whl → 24.9.24__py2.py3-none-any.whl

{PyGeodesy-24.9.9.dist-info → PyGeodesy-24.9.24.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PyGeodesy
-Version: 24.9.9
+Version: 24.9.24
 Summary: Pure Python geodesy tools
 Home-page: https://GitHub.com/mrJean1/PyGeodesy
 Author: Jean M. Brouwers
@@ -115,7 +115,7 @@ line: ``epydoc --html --no-private --no-source --name=PyGeodesy --url=... -v pyg
 Tests
 =====
 
-The tests ran with Python 3.13.0rc1, 3.12.5 (with geographiclib_ 2.0, numpy_ 2.1.0, scipy_ 1.14.1,
+The tests ran with Python 3.13.0rc2, 3.12.6 (with geographiclib_ 2.0, numpy_ 2.1.0, scipy_ 1.14.1,
 GeodSolve_ 2.3, IntersectTool_ 2.3 and RhumbSolve_ 2.3), 3.11.5 (with geographiclib_ 2.0, numpy_
 1.24.2 and scipy_ 1.10.1), Python 3.10.8 (with geographiclib_ 2.0, numpy_ 1.23.3, scipy_ 1.9.1,
 GeoConvert_ 2.3, GeodSolve_ 2.3), Python 3.9.6 and Python 2.7.18 (with geographiclib_ 1.50, numpy_
@@ -126,13 +126,13 @@ All tests ran with and without ``lazy import`` for Python 3 and with command lin
 and env variable ``PYGEODESY_WARNINGS=on`` for all Python versions.  The results of those tests are
 included in the distribution files.
 
-Python 3.13.0rc1, 3.12.5, 3.11.5 and 3.10.8 run on Apple M1 Silicon (``arm64``), *natively*.  Python 2.7.18
+Python 3.13.0rc2, 3.12.6, 3.11.5 and 3.10.8 run on Apple M1 Silicon (``arm64``), *natively*.  Python 2.7.18
 runs on Intel (``x86_64``) or Intel *emulation* (\"``arm64_x86_64``\", see function `pygeodesy.machine`_).
 
-Test coverage has been measured with coverage_ 7.6.0 using Python 3.12.5, 3.11.5 and 3.10.8.  The complete
+Test coverage has been measured with coverage_ 7.6.0 using Python 3.12.6, 3.11.5 and 3.10.8.  The complete
 coverage report in HTML and a PDF summary are included in the distribution files.
 
-The tests also ran with Python 3.12.5 (and geographiclib_ 2.0) on `Debian 12`_ in 64-bit only and with
+The tests also ran with Python 3.12.6 (and geographiclib_ 2.0) on `Debian 12`_ in 64-bit only and with
 Python 3.12.3 (and geographiclib_ 2.0) on `Windows 2019Server`_ in 64-bit only and with Python 2.7.18
 (and with geographiclib_ 1.52) on `Windows 10`_ in 64- and 32-bit.
 
@@ -160,7 +160,7 @@ and McCabe_ using Python 2.7.18 and with Flake8_ using Python 3.11.5, both in 64
 
 For a summary of all *Karney*-based functionality in ``pygeodesy``, see module karney_.
 
-*Last updated: Sep 09, 2024.*
+*Last updated: Sep 24, 2024.*
 
 License
 =======

{PyGeodesy-24.9.9.dist-info → PyGeodesy-24.9.24.dist-info}/RECORD

@@ -1,13 +1,13 @@
 pygeodesy/LICENSE,sha256=YfgAiyxOwY6P9Kkb1_5XN81nueTLrpb3Ffkv3EuPgFU,1144
-pygeodesy/__init__.py,sha256=ohpvlHHdgFq_uc-u1gq7deCyTlzB_nq-Lc8tZGWN-Go,42225
+pygeodesy/__init__.py,sha256=gu-FjebeihoF1kw47l9X8UjHwQsI4lWn6n4_BPVu6hE,42475
 pygeodesy/__main__.py,sha256=OupAiDFT_RXFw1Lh6C8bVG9NaonV_Wy6jibN_hA4aJg,5442
 pygeodesy/albers.py,sha256=JXqeWAteV2N601e5R-D7bCsa8qAYv3oWm7M3R4uszXM,31078
 pygeodesy/azimuthal.py,sha256=erI5DTM0-Nr5YC2LnBTVNbUKlrnnXJUwKe0SCQ7-Soc,50223
-pygeodesy/basics.py,sha256=5fuWL_azUd6DxSGUwB4RJERH_ZrBB0vVqDc9wvxpIHo,30116
-pygeodesy/booleans.py,sha256=mWNdmk5cYfcH3v0sEpU3Aw5sFLp_dq5mXJodzenOEm0,73506
+pygeodesy/basics.py,sha256=pY0SAIRoD_4okhVwNGAAhxEa9FIZI07XbyPu_sn6sNM,30121
+pygeodesy/booleans.py,sha256=H1QbPYpelIGUzfOdwWd2hWXxi6vvbrQiEV_rK5cg1LA,73491
 pygeodesy/cartesianBase.py,sha256=jzPJR4_3Yp5jmK7A5inGOXskDebs7YMdjQzINHjzkP8,46935
-pygeodesy/clipy.py,sha256=3ysIEWQ89HjvLnSlZm0nZAt9EN5AQxIs1sLLM8yF9zM,27675
-pygeodesy/constants.py,sha256=vW6CWRa62hHJqmKhlPxI4Eu-cbLhoUYWYTd9ctfFo-8,19240
+pygeodesy/clipy.py,sha256=hzI6IutTRkbssjPPgAjXv9OcgCIIZchjnMWaHin4YAU,27728
+pygeodesy/constants.py,sha256=l5MDpwvrl_BhavQHSPdbqT6kt7YP1O78b6rIl43sVAo,19616
 pygeodesy/css.py,sha256=rPr5EcGm_blVt72h-851ytpVxdPhx1MEqqGs3Ppx3eI,25651
 pygeodesy/datums.py,sha256=-m9-cvUVOaJ92eeFgc3O54lel3WqnVYcS7o6WCXkNyw,34043
 pygeodesy/dms.py,sha256=QxiRtyc8QSiW-egm6Mrcrgmi_deWz3YnYwIm82JS4TU,44443
@@ -23,13 +23,13 @@ pygeodesy/ellipsoidalVincenty.py,sha256=hSQ7s6QGDJSsrZvP7KwM1Q4wJz0YjyFjODF6XPGZ
 pygeodesy/ellipsoids.py,sha256=i7piclEb17hYt5wqL7Q0p5qLuRVvz3yonUSUcvmiozc,108005
 pygeodesy/elliptic.py,sha256=M30EjB1Il_5Y8E94rN5hz4zPhYBQNoCo5Lj0TvCkiRE,45010
 pygeodesy/epsg.py,sha256=Vq_M1g7Z0x5Noqsmt6OCQkPdYJj4jpWObtixL5pAjDA,8220
-pygeodesy/errors.py,sha256=-UnQ1GhkiTR6uBd0z0vd_pFEmjWXlOESRDwMdxY_JD0,30605
+pygeodesy/errors.py,sha256=6aofErGAi6l1xwXqd4PU7YSoPQiutc-gFiFiTamGoU4,31597
 pygeodesy/etm.py,sha256=G6dtN1LQAoft5mcUoxVgxHPnl9s8NBsXl-IlS8WP7jo,45328
-pygeodesy/fmath.py,sha256=7Bokavxcu4oFGMKHCn4JxbETDisbBkWJqJx-KiCtHPM,34374
+pygeodesy/fmath.py,sha256=JEjiBJKjNt0Lm3EAIn_IiRoxtuTZbXNz2nscjJPEzUA,35376
 pygeodesy/formy.py,sha256=PrAlrSASKIXeFvawueyo7kRC-UJJSEQR3APhirHqhtE,75032
 pygeodesy/frechet.py,sha256=F6TYbOEYSOFcll6nvp4GhqafUWXh1umAKb1QaMYQtuM,34372
-pygeodesy/fstats.py,sha256=4eIhzSCMnn-2nVP3Q-5cCWCI-rvH0CE-zkk9hp3Pvn8,28642
-pygeodesy/fsums.py,sha256=Zk65c445gAz7Ld_iXL1zZ3ZtkA_gM2edDOf9nJSEnXE,87433
+pygeodesy/fstats.py,sha256=nV5Kj5DrKO6zjU18rxlbTRwM-m7Rct9GLhOuH3_3hl8,28541
+pygeodesy/fsums.py,sha256=cVNCB90BhvBfTCpZyvvPFVa9aEllzKd-CPJLPPeTFx4,99816
 pygeodesy/gars.py,sha256=KRlyZUv1euEl7ZkDYSZG913GAfFMRG2m2hp37ObYY_Y,11829
 pygeodesy/geodesici.py,sha256=I9e7saiO5lRodOXJk5DSQxladVcGfT0FR6KJSJVCw7Y,75015
 pygeodesy/geodesicw.py,sha256=oEs96Epw839BH3rVpxVdpXuQRBudbOftUsSBUid8fcc,28710
@@ -41,10 +41,10 @@ pygeodesy/heights.py,sha256=Gv3b1xC7I0tBvTO-DdmWVW8A3ZV3RKxY3jveUUZso38,42909
 pygeodesy/internals.py,sha256=TSJueckcn6dQddZUhLR0Vf9XN-Whm85IYmSvhVZbB_s,22550
 pygeodesy/interns.py,sha256=XvNyWODsUd8_vBLX9XLmcU9Bq97L0nahhF-gvLivzKA,23278
 pygeodesy/iters.py,sha256=xPH0ytomvJeKHu717Fjzi3QbRrEikO6n7vqkY_IDz1k,20311
-pygeodesy/karney.py,sha256=2G5Y9FvZlb8cTrDHtHF71n3gxCvGLxlWUgVwXGmAnjA,39033
+pygeodesy/karney.py,sha256=ZkZhbsieGnr2cAIXjzmK-9OJdh72Hjp6Hya1wIoVP3c,38913
 pygeodesy/ktm.py,sha256=eUWYKhoPou1AiWiNz-xI8eDr-q-RvKDY6newAmYP26Q,27239
 pygeodesy/latlonBase.py,sha256=kCK6IcVg4oeGPYkIU0rgKbZLCRo36wC9NtY5OKBlpeQ,79106
-pygeodesy/lazily.py,sha256=jIfuqXPh-Gvg4H2XeGgbKoOGuc49SR_sMVMoBEVq4LQ,47370
+pygeodesy/lazily.py,sha256=NpWnTIAOkIxMorgUMvyR7d3D65-55DqaU98eBnwxcrQ,47411
 pygeodesy/lcc.py,sha256=oDooTzFHnqsMtSKmmmfEbypHpdlPltfo7Mt0pAO4-_A,25669
 pygeodesy/ltp.py,sha256=ZMdfJHLKeDHSoymmjWfoAZTtrdoADf0O6WccRf53L2w,50745
 pygeodesy/ltpTuples.py,sha256=oNTxtjxmHveMDQvPdQCHl-JzSJJGb1jFRrF8SEG_Ojw,58480
@@ -53,14 +53,14 @@ pygeodesy/named.py,sha256=-uqZk54BC4rFtuVh4LwfWty704DVPLEIwI1EzpNImIg,52350
 pygeodesy/namedTuples.py,sha256=e0rV1Zb2tbx8fztXcSkUFsRTmZtye8B5i_Hi4vBUndU,28758
 pygeodesy/nvectorBase.py,sha256=rClYu26Bof9DsM1F60i0eccktNo8xYA_vkThwoY7N4o,28644
 pygeodesy/osgr.py,sha256=3Qon_oWuVrIb1vYnd3ae3L5MDkjp_tngyVdiYHuL4to,30821
-pygeodesy/points.py,sha256=e8ri5Pu5snpN4P6pPbGxmW_yx4-ApD9mOBuR0REoflA,64409
+pygeodesy/points.py,sha256=vtPQdL2kLOenIrpOdGjOYrUWsav3I1JXQEJugMv46eY,64383
 pygeodesy/props.py,sha256=3BIkcaEDwyhSMJle6ufBaji0R0FrD64G3nuR4tPsnrU,23819
-pygeodesy/resections.py,sha256=fFthuI9EoO4WdY25Y_rCDN56gz2mPFsUypR3oGlh8Ks,43556
+pygeodesy/resections.py,sha256=IrPHAHbF9JdBgAVd9LMgXoPgPFRkY8j1G7p0P8y8pdc,43543
 pygeodesy/simplify.py,sha256=vkOtO478sOjZHZBuSFvCaRh72qSDQuzXLQpMK1LKyYI,25232
 pygeodesy/solveBase.py,sha256=oI_-FFjSPI1_yQ598hF6HFKUa4tLsAfbnWtKoQ-TJgE,19206
 pygeodesy/sphericalBase.py,sha256=bDVuArfePzqcothiTTkyisUGrMDTallXKROjihotcn8,32239
-pygeodesy/sphericalNvector.py,sha256=MiWcCy5CdvXsCa3963GeUMEnSYGhBiXzkoAEX_tGX4Q,56952
-pygeodesy/sphericalTrigonometry.py,sha256=F84lPOSLj4Z2-f8xMjpqvTISuBpsjFwdF8fkuaEi5Hs,64077
+pygeodesy/sphericalNvector.py,sha256=Ip0jPv-MXBoHf3f-uXitGqiMz2xynDtE15OSJEl268Q,56913
+pygeodesy/sphericalTrigonometry.py,sha256=cfXuJaJJ6jrJZxQ5anFHhz3ZbgkAB7aiMfuyMtJXqsE,64064
 pygeodesy/streprs.py,sha256=1szfCzFlVgsyJ0H_QFGtzyMfOJpyIZNfj34qvWNACM4,23041
 pygeodesy/trf.py,sha256=CRFHmu0lFNtc41fddXavEkHmgnTFWEqR4RTFodcFKsw,119022
 pygeodesy/triaxials.py,sha256=a9vdXmSqUXOneL0WpLriEFwNMW9EB6MpVbKYXuw_wi8,62562
@@ -87,12 +87,12 @@ pygeodesy/auxilats/auxDLat.py,sha256=3sq7Wc-ykpVbldNh10IZ3tIlJ7DcQFDlmZ4nOhJEgo0
 pygeodesy/auxilats/auxDST.py,sha256=BSHR0M3hQBR1KjF_BUJHrNzGXr2PblYXUQsMD5L1Kcc,10468
 pygeodesy/auxilats/auxLat.py,sha256=7uxe7mtkY-9uSI3z2wmUNBqzCKyFujIpJU1v-8yDgds,31934
 pygeodesy/auxilats/auxily.py,sha256=X1y4_qC40Th7CgmZILqoQCLAUQoyORcR5W-bN4X0W4I,7882
-pygeodesy/deprecated/__init__.py,sha256=cBYqbh5Njkd5SnaJDaMUP6PL66HfLWii_JoLjnKfnP4,2815
+pygeodesy/deprecated/__init__.py,sha256=jSwmipBzyIOYIQhiZxZwsIQlqMXPfNaPfVpppBu7zkw,2815
 pygeodesy/deprecated/bases.py,sha256=nxUtxscAg5cb7HEzTSmPJFafLc9Ad2SL5ovwkYaU5z8,1655
 pygeodesy/deprecated/classes.py,sha256=gktUr5qYXCm8ztH6lfzFgaHKqMiw1tHc1Vw6i6DEvNA,12377
 pygeodesy/deprecated/consterns.py,sha256=W-rQiw443_zCFEUzjCGykGND-i5P_ghEIyJUSGg7-Nc,1908
 pygeodesy/deprecated/datum.py,sha256=s8Hke053RucP5ACvg3gdX8weQfGM0x59pQyHpY11eOs,1875
-pygeodesy/deprecated/functions.py,sha256=EZMZe9WGkcbwmLDoWO2xKJyFlUkWAAMxmlbtCLa28lM,13866
+pygeodesy/deprecated/functions.py,sha256=-spa_5k7yRWGd7s4UBObzf9QpAM0o-tgP4S0wX876JY,14122
 pygeodesy/deprecated/nvector.py,sha256=_dZf84RdmPr7e7i8Yi683R-8Mqn0ipTqIsDhNl29hGo,2117
 pygeodesy/deprecated/rhumbBase.py,sha256=d7YKW--dnRpobB2CSvhFsntWF0B5-HYJqID2pv5lD5k,1375
 pygeodesy/deprecated/rhumbaux.py,sha256=n29xazYMbVG4MJkLjAGu8B88bbcFF5xnvuj9xSxjfIk,1428
@@ -112,7 +112,7 @@ pygeodesy/rhumb/aux_.py,sha256=XN0zlNwhU525NdcHko-G_Igvc7-mzXbT7ypcVIt5ADw,16048
 pygeodesy/rhumb/bases.py,sha256=awpP60K6EOzS753oXM-h9S6paS6FwoCnco8Zt8n81jY,54162
 pygeodesy/rhumb/ekx.py,sha256=JsmdCatWKqsffF8jCh4HDu5IKt_kHq1qUdmP1LMuofY,24044
 pygeodesy/rhumb/solve.py,sha256=d4RujlW8B7BFAHlrwfj2Bup-uXZZTpgKWgI2JJBqUGc,24059
-PyGeodesy-24.9.9.dist-info/METADATA,sha256=NZ7vpJ3lZq-5fmen4MR9YO3yRCcEuJN8zhpO8QtbKe0,19876
-PyGeodesy-24.9.9.dist-info/WHEEL,sha256=fS9sRbCBHs7VFcwJLnLXN1MZRR0_TVTxvXKzOnaSFs8,110
-PyGeodesy-24.9.9.dist-info/top_level.txt,sha256=cEQPatCXzKZqrivpULC5V5fuy9_V_bAwaP_gUGid7pQ,10
-PyGeodesy-24.9.9.dist-info/RECORD,,
+PyGeodesy-24.9.24.dist-info/METADATA,sha256=tIX957Sljzlxz05L-pGuqjt4ZCcfJnLuQ9E0mTwWruw,19877
+PyGeodesy-24.9.24.dist-info/WHEEL,sha256=fS9sRbCBHs7VFcwJLnLXN1MZRR0_TVTxvXKzOnaSFs8,110
+PyGeodesy-24.9.24.dist-info/top_level.txt,sha256=cEQPatCXzKZqrivpULC5V5fuy9_V_bAwaP_gUGid7pQ,10
+PyGeodesy-24.9.24.dist-info/RECORD,,

pygeodesy/__init__.py

@@ -122,7 +122,7 @@ C{epydoc --html --no-private --no-source --name=PyGeodesy --url=... -v pygeodesy
 Tests
 =====
 
-The tests ran with Python 3.13.0rc1, Python 3.12.5 (with U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0,
+The tests ran with Python 3.13.0rc2, Python 3.12.6 (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,
 U{IntersectTool<https://GeographicLib.SourceForge.io/C++/doc/utilities.html>} 2.3 and
@@ -144,13 +144,13 @@ All tests ran with and without C{lazy import} for Python 3 and with command line
 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.0 using Python 3.12.5,
+Test coverage has been measured with U{coverage<https://PyPI.org/project/coverage>} 7.6.0 using Python 3.12.6,
 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.0rc1, 3.12.5, 3.11.5 and 3.10.8 run on Apple M1 Silicon (C{arm64}), I{natively}.  Python 2.7.18 runs on
+Python 3.13.0rc2, 3.12.6, 3.11.5 and 3.10.8 run on Apple M1 Silicon (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.12.5 (and U{geographiclib<https://PyPI.org/project/geographiclib>} 2.0) on U{Debian
+The tests also ran with Python 3.12.6 (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.3 (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
@@ -199,6 +199,8 @@ The following environment variables are observed by C{PyGeodesy}:
 
  - C{PYGEODESY_EXCEPTION_CHAINING} - see module L{errors<pygeodesy.errors>}.
  - C{PYGEODESY_FMT_FORM} - see module L{dms<pygeodesy.dms>}.
+ - C{PYGEODESY_FSUM_F2PRODUCT} - see module L{fsums<pygeodesy.fsums>} and method L{f2product<pygeodesy.Fsum.f2product>}.
+ - C{PYGEODESY_FSUM_NONFINITES} - see module L{fsums<pygeodesy.fsums>} and function L{nonfiniterrors<pygeodesy.nonfiniterrors>}.
  - C{PYGEODESY_FSUM_RESIDUAL} - see module L{fsums<pygeodesy.fsums>} and class L{Fsum<pygeodesy.Fsum>}.
  - C{PYGEODESY_GEOCONVERT} - see module L{mgrs<pygeodesy.mgrs>}.
  - C{PYGEODESY_GEODSOLVE} - see module L{geodsolve<pygeodesy.geodsolve>}.
@@ -595,7 +597,7 @@ else:
 
 from pygeodesy.internals import _version2,  _DOT_  # PYCHOK import
 # from pygeodesy.interns import _DOT_  # from .internals
-__version__ = '24.09.09'
+__version__ = '24.09.24'
 # see setup.py for similar logic
 version     = _DOT_(*_version2(__version__, n=3))
 

pygeodesy/basics.py

@@ -19,6 +19,7 @@ del division
 from pygeodesy.errors import _AttributeError, _ImportError, _NotImplementedError, \
                              _TypeError, _TypesError, _ValueError, _xAssertionError, \
                              _xkwds_get1
+# from pygeodesy.fsums import _isFsumTuple  # _MODS
 from pygeodesy.internals import _0_0, _enquote, _passarg, _version_info
 from pygeodesy.interns import MISSING, NN, _1_, _by_, _COMMA_, _DOT_, _DEPRECATED_, \
                              _ELLIPSIS4_, _EQUAL_, _in_, _invalid_, _N_A_, _not_, \
@@ -37,7 +38,7 @@ from math import copysign as _copysign
 import inspect as _inspect
 
 __all__ = _ALL_LAZY.basics
-__version__ = '24.09.02'
+__version__ = '24.09.12'
 
 _below_               = 'below'
 _list_tuple_types     = (list, tuple)
@@ -460,18 +461,16 @@ def isscalar(obj, both=False):
     '''Is B{C{obj}}ect an C{int} or integer C{float} value?
 
        @arg obj: The object (any C{type}).
-       @kwarg both: If C{True}, check L{Fsum<Fsum.residual>}.
+       @kwarg both: If C{True}, check L{Fsum} and L{Fsum2Tuple}
+                    residuals.
 
-       @return: C{True} if C{int}, C{float} or L{Fsum} with
-                zero residual, C{False} otherwise.
+       @return: C{True} if C{int}, C{float} or C{Fsum/-2Tuple}
+                with zero residual, C{False} otherwise.
     '''
     if isinstance(obj, _Scalars):
-        return not isbool(obj)
-    elif both:  # and isinstance(obj, Fsum)
-        try:
-            return bool(obj.residual == 0)
-        except (AttributeError, TypeError):
-            pass  # XXX float(int(obj)) == obj?
+        return not isbool(obj)  # exclude bool
+    elif both and _MODS.fsums._isFsumTuple(obj):
+        return bool(obj.residual == 0)
     return False
 
 

pygeodesy/booleans.py

@@ -43,7 +43,7 @@ from pygeodesy.utily import fabs, _unrollon, _Wrap
 # from math import fabs  # from .utily
 
 __all__ = _ALL_LAZY.booleans
-__version__ = '24.08.30'
+__version__ = '24.09.23'
 
 _0_EPS =  EPS  # near-zero, positive
 _EPS_0 = -EPS  # near-zero, negative
@@ -1117,7 +1117,7 @@ class _CompositeBase(_Named):
 
     def _sum1(self, _a_p, *args, **kwds):  # in .karney, .points
         # Sum the area or perimeter of all clips
-        return _MODS.fsums.fsum1((_a_p(c, *args, **kwds) for c in self._clips), floats=True)
+        return _MODS.fsums.fsum1((_a_p(c, *args, **kwds) for c in self._clips))
 
     def _sum2(self, LL, _a_p, *args, **kwds):  # in .sphericalNvector, -Trigonometry
         # Sum the area or perimeter of all clips
@@ -1127,7 +1127,7 @@ class _CompositeBase(_Named):
             for v in clip:
                 yield _LL(v.lat, v.lon)  # datum=Sphere
 
-        return _MODS.fsums.fsum1((_a_p(_lls(c), *args, **kwds) for c in self._clips), floats=True)
+        return _MODS.fsums.fsum1((_a_p(_lls(c), *args, **kwds) for c in self._clips))
 
     def toLatLon(self, LatLon=None, closed=False, **LatLon_kwds):
         '''Yield all (non-duplicate) points and intersections
@@ -1187,7 +1187,7 @@ class _CompositeFHP(_CompositeBase):
             n, b = v, v._label
             if b in L.RIGHT_LEFT_ON:  # next chain
                 while True:
-                    n._label = None  # n.__dict__.pop('_label')
+                    n._label = None  # _xkwds_pop(n.__dict__, _label=None)
                     n = n._next
                     if n is v or n._label is not L.ON_ON:  # n._label and ...
                         break

pygeodesy/clipy.py

@@ -15,22 +15,22 @@ from __future__ import division as _; del _  # PYCHOK semicolon
 # from pygeodesy.basics import len2  # from .fmath
 from pygeodesy.constants import EPS, _0_0, _1_0
 from pygeodesy.errors import _AssertionError, ClipError, PointsError
-from pygeodesy.fmath import fabs, len2
-from pygeodesy.fsums import fsumf_, Property_RO
+from pygeodesy.fmath import fabs, len2,  Fsum
+# from pygeodesy.fsums import Fsum  # from .fmath
 from pygeodesy.interns import NN, _clipid_, _convex_, _DOT_, _end_, _few_, \
                              _fi_, _height_, _i_, _invalid_, _j_, _lat_, \
                              _lon_, _near_, _not_, _points_, _start_, _too_
 from pygeodesy.iters import _imdex2, points2
 from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS
-from pygeodesy.named import _Named, _NamedTuple, _Pass
+from pygeodesy.named import _Named, _NamedTuple, _Pass,  Property_RO
 from pygeodesy.points import areaOf, boundsOf, isconvex_, LatLon_
-# from pygeodesy.props import Property_RO  # from .fsums
+# from pygeodesy.props import Property_RO  # from .named
 from pygeodesy.units import Bool, FIx, HeightX, Lat, Lon, Number_
 
 # from math import fabs  # from .fmath
 
 __all__ = _ALL_LAZY.clipy
-__version__ = '23.05.18'
+__version__ = '23.09.20'
 
 _fj_       = 'fj'
 _original_ = 'original'
@@ -576,10 +576,10 @@ class _SH(_Named):
         x, dx = p1.lon, self._dx
         fy = float(p2.lat - y)
         fx = float(p2.lon - x)
-        d  = fy * dx - fx * dy
+        d  = fy * dx - fx * dy  # fdot((fx, fy), dx, -dy)
         if fabs(d) < EPS:  # PYCHOK no cover
             raise _AssertionError(self._DOT_(self.intersect.__name__))
-        d  = fsumf_(self._xy, -y * dx, x * dy) / d
+        d  = Fsum(self._xy, -y * dx, x * dy).fover(d)
         y += d * fy
         x += d * fx
         return _SHlli(y, x, p1.classof, edge)

pygeodesy/constants.py

@@ -13,9 +13,9 @@ from pygeodesy.basics import _copysign, isbool, iscomplex, isint,  _0_0
 from pygeodesy.errors import _xError, _xError2, _xkwds_get1, _xkwds_item2
 # from pygeodesy.internals import _0_0  # from .basics
 from pygeodesy.interns import _INF_, _NAN_, _UNDER_
-# from pygeodesy.lazily import _ALL_LAZY  # from .unitsBase
+from pygeodesy.lazily import _ALL_MODS as _MODS, _ALL_LAZY
 # from pygeodesy.streprs import Fmt  # from .unitsBase
-from pygeodesy.unitsBase import Float, Int, Radius,  _ALL_LAZY, Fmt
+from pygeodesy.unitsBase import Float, Int, Radius,  Fmt
 
 from math import fabs, isinf, isnan, pi as _pi, sqrt
 try:
@@ -24,7 +24,7 @@ except ImportError:  # Python 2-
     _inf, _nan = float(_INF_), float(_NAN_)
 
 __all__ = _ALL_LAZY.constants
-__version__ = '24.08.26'
+__version__ = '24.09.24'
 
 
 def _copysign_0_0(y):
@@ -320,20 +320,30 @@ def isclose(a, b, rel_tol=1e-12, abs_tol=EPS0):
     return _isclose(a, b, rel_tol=rel_tol, abs_tol=abs_tol)
 
 
+try:
+    from cmath import isfinite as _iscfinite
+except ImportError:  # Python 3.1-
+
+    def _iscfinite(x):  # PYCHOK not self?
+        '''Mimick Python 3.2+ C{cmath.isfinite}.
+        '''
+        return _isfinite(x.real) and _isfinite(x.imag)
+
 try:
     from math import isfinite as _isfinite  # in .ellipsoids, ...
 except ImportError:  # Python 3.1-
 
-    def _isfinite(x):
+    def _isfinite(x):  # PYCHOK not self?
         '''Mimick Python 3.2+ C{math.isfinite}.
         '''
         return not (isinf(x) or isnan(x))
 
 
 def isfinite(obj):
-    '''Check a finite C{scalar} or C{complex} value.
+    '''Check a finite C{scalar}, C{complex}, ... value.
 
-       @arg obj: Value (C{scalar} or C{complex}).
+       @arg obj: Value (C{scalar}, C{complex}, an L{Fsum} or
+                        L{Fsum2Tuple}).
 
        @return: C{False} if B{C{obj}} is C{INF}, C{NINF}
                 or C{NAN}, C{True} otherwise.
@@ -344,7 +354,9 @@ def isfinite(obj):
         return (obj not in _INF_NAN_NINF) and _isfinite(obj)
     except Exception as x:
         if iscomplex(obj):  # _isfinite(complex) thows TypeError
-            return isfinite(obj.real) and isfinite(obj.imag)
+            return _iscfinite(obj)
+        if _MODS.fsums._isFsumTuple(obj):  # OverflowError
+            return obj.is_finite()
         raise _xError(x, Fmt.PAREN(isfinite.__name__, obj))
 
 

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__ = '24.07.02'
+__version__ = '24.09.19'
 
 if _unLazy0:
     from pygeodesy.deprecated import bases, datum, nvector, rhumbBase, \

pygeodesy/deprecated/functions.py

@@ -14,7 +14,7 @@ from pygeodesy.props import deprecated_function
 from pygeodesy.units import Number_, Scalar_
 
 __all__ = _ALL_DEPRECATED.deprecated_functions
-__version__ = '24.06.11'
+__version__ = '24.09.19'
 
 _WGS84 = _UTM = object()
 
@@ -193,6 +193,14 @@ def fStrzs(floatstr):  # PYCHOK no cover
     return _MODS.streprs.fstrzs(floatstr)
 
 
+@deprecated_function
+def Fsum2product(*xs, **kwds):
+    '''DEPRECATED on 2024.09.19, use L{Fsum}C{(*B{xs}, B{f2product}=True, ...)}.'''
+    F = _MODS.fsums.Fsum(**kwds)
+    _ = F.f2product(True)
+    return F._facc_args(xs, up=False) if xs else F
+
+
 @deprecated_function
 def hypot3(x, y, z):  # PYCHOK no cover
     '''DEPRECATED, use function L{pygeodesy.hypot_}.'''

pygeodesy/errors.py

@@ -27,7 +27,7 @@ from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS, _getenv, _PYTHON_X_D
 from copy import copy as _copy
 
 __all__ = _ALL_LAZY.errors  # _ALL_DOCS('_InvalidError', '_IsnotError')  _under
-__version__ = '24.07.07'
+__version__ = '24.09.19'
 
 _argument_   = 'argument'
 _box_        = 'box'
@@ -759,14 +759,16 @@ try:
     _ = {}.__or__  # {} | {}  # Python 3.9+
 
     def _xkwds(kwds, **dflts):
-        '''(INTERNAL) Update C{dflts} with specified C{kwds}.
+        '''(INTERNAL) Update C{dflts} with specified C{kwds},
+           i.e. C{copy(kwds).update(dflts)}.
         '''
-        return (dflts | kwds) if kwds else dflts
+        return ((dflts | kwds) if kwds else dflts) if dflts else kwds
 
 except AttributeError:
 
     def _xkwds(kwds, **dflts):  # PYCHOK expected
-        '''(INTERNAL) Update C{dflts} with specified C{kwds}.
+        '''(INTERNAL) Update C{dflts} with specified C{kwds},
+           i.e. C{copy(kwds).update(dflts)}.
         '''
         d = dflts
         if kwds:
@@ -787,6 +789,20 @@ except AttributeError:
 #             setattr(inst, NN(_UNDER_, n), v)
 
 
+def _xkwds_from(orig, *args, **kwds):
+    '''(INTERNAL) Return the items from C{orig} with the keys
+       from C{kwds} and a value not in C{args} and C{kwds}.
+    '''
+    def _items(orig, args, items):
+        for n, m in items:
+            if n in orig:  # n in (orig.keys() & kwds.keys())
+                t = orig[n]
+                if t is not m and t not in args:
+                    yield n, t
+
+    return _items(orig, args, kwds.items())
+
+
 def _xkwds_get(kwds, **name_default):
     '''(INTERNAL) Get a C{kwds} value by C{name} or the
        C{default} if not present.
@@ -809,7 +825,8 @@ def _xkwds_get_(kwds, **names_defaults):
 
 def _xkwds_get1(kwds, **name_default):
     '''(INTERNAL) Get one C{kwds} value by C{name} or the
-       C{default} if not present.
+       C{default} if not present.  Raise an C{_UnexpectedError}
+       with any remaining keyword arguments.
     '''
     v, kwds = _xkwds_pop2(kwds, **name_default)
     if kwds:
@@ -842,6 +859,16 @@ def _xkwds_not(*args, **kwds):
     return dict((n, v) for n, v in kwds.items() if v not in args)
 
 
+def _xkwds_pop(kwds, **name_default):
+    '''(INTERNAL) Pop an item by C{name} from C{kwds} and
+       return its value, otherwise return the C{default}.
+    '''
+    if isinstance(kwds, dict) and len(name_default) == 1:
+        for n, v in name_default.items():
+            return kwds.pop(n, v)
+    raise _xAssertionError(_xkwds_pop, kwds, **name_default)
+
+
 def _xkwds_pop2(kwds, **name_default):
     '''(INTERNAL) Pop a C{kwds} item by C{name} and return the value and
        reduced C{kwds} copy, otherwise the C{default} and original C{kwds}.

pygeodesy/fmath.py

@@ -10,11 +10,11 @@ from pygeodesy.basics import _copysign, copysign0, isbool, isint, isscalar, \
                               len2, map1, _xiterable
 from pygeodesy.constants import EPS0, EPS02, EPS1, NAN, PI, PI_2, PI_4, \
                                _0_0, _0_125, _1_6th, _0_25, _1_3rd, _0_5, _1_0, \
-                               _N_1_0, _1_5, _copysign_0_0, _isfinite, remainder
+                               _N_1_0, _1_5, _copysign_0_0, isfinite, remainder
 from pygeodesy.errors import _IsnotError, LenError, _TypeError, _ValueError, \
                              _xError, _xkwds_get1, _xkwds_pop2
-from pygeodesy.fsums import _2float, Fsum, fsum, fsum1_, _isFsumTuple, _1primed, \
-                            Fmt, unstr
+from pygeodesy.fsums import _2float, Fsum, fsum, fsum1_, _isFsumTuple, \
+                            _1primed, _Psum_,  Fmt, unstr
 from pygeodesy.interns import MISSING, _negative_, _not_scalar_
 from pygeodesy.lazily import _ALL_LAZY, _sys_version_info2
 # from pygeodesy.streprs import Fmt, unstr  # from .fsums
@@ -24,7 +24,7 @@ from math import fabs, sqrt  # pow
 import operator as _operator  # in .datums, .trf, .utm
 
 __all__ = _ALL_LAZY.fmath
-__version__ = '24.09.09'
+__version__ = '24.09.19'
 
 # sqrt(2) - 1 <https://WikiPedia.org/wiki/Square_root_of_2>
 _0_4142  =  0.41421356237309504880  # ... ~ 3730904090310553 / 9007199254740992
@@ -431,8 +431,9 @@ def fdot(a, *b):
 
        @raise LenError: Unequal C{len(B{a})} and C{len(B{b})}.
 
-       @see: Class L{Fdot} and U{Algorithm 5.10 B{DotK}
-             <https://www.TUHH.De/ti3/paper/rump/OgRuOi05.pdf>}.
+       @see: Class L{Fdot}, U{Algorithm 5.10 B{DotK}
+             <https://www.TUHH.De/ti3/paper/rump/OgRuOi05.pdf>} and function
+             C{math.sumprod} in Python 3.12 and later.
     '''
     return fsum(_map_mul(a, b, fdot))
 
@@ -544,16 +545,20 @@ def fidw(xs, ds, beta=2):
     return x
 
 
-def fma(x, y, z):
-    '''Fused-multiply-add, as C{math.fma(x, y, z)} from Python 3.13+.
+def fma(x, y, z, **raiser):
+    '''Fused-multiply-add, using C{math.fma(x, y, z)} from Python 3.13+
+       or an equivalent implementation.
 
-       @arg x: A C{scalar}, an L{Fsum} or L{Fsum2Tuple} instance.
-       @arg y: A C{scalar}, an L{Fsum} or L{Fsum2Tuple} instance.
-       @arg z: A C{scalar}, an L{Fsum} or L{Fsum2Tuple} instance.
+       @arg x: Multiplicand (C{scalar}, an L{Fsum} or L{Fsum2Tuple}).
+       @arg y: Multiplier (C{scalar}, an L{Fsum} or L{Fsum2Tuple}).
+       @arg z: Addend (C{scalar}, an L{Fsum} or L{Fsum2Tuple}).
+       @kwarg raiser: Keyword argument C{B{raiser}=False}, if C{True},
+                      throw an exception, otherwise pass any non-finite
+                      result (C{bool}).
 
-       @return: C{(x * y) + z} (C{float} or an L{Fsum}).
+       @return: C{(x * y) + z} (L{Fsum} or C{float}).
     '''
-    return Fsum(x).fma(y, z).as_iscalar
+    return _Psum_(x).fma(y, z, raiser=raiser).as_iscalar
 
 
 def fmean(xs):
@@ -581,6 +586,22 @@ def fmean_(*xs):
     return fmean(xs)
 
 
+def f2mul_(x, *ys, **raiser):
+    '''Cascaded, accurate multiplication C{B{x} * B{y} * B{y} ...} for all B{C{ys}}.
+
+       @arg x: Multiplicand (C{scalar}, an L{Fsum} or L{Fsum2Tuple}).
+       @arg ys: Multipliers (each C{scalar}, an L{Fsum} or L{Fsum2Tuple}), all
+                positional.
+       @kwarg raiser: Keyword argument C{B{raiser}=False}, if C{True}, throw an
+                      exception, otherwise pass any non-finite result (C{bool}).
+
+       @return: The cascaded I{TwoProduct} (L{Fsum}, C{float} or C{int}).
+
+       @see: U{Equations 2.3<https://www.TUHH.De/ti3/paper/rump/OzOgRuOi06.pdf>}
+    '''
+    return _Psum_(x).f2mul_(*ys, **raiser).as_iscalar
+
+
 def fpolynomial(x, *cs, **over):
     '''Evaluate the polynomial M{sum(cs[i] * x**i for i=0..len(cs))
        [/ over]}.
@@ -615,7 +636,7 @@ def fpowers(x, n, alts=0):
     elif n < 1:
         raise _ValueError(n=n)
 
-    p  = x if isint(x) or _isFsumTuple(x) else _2float(x=x)
+    p  = x if isscalar(x) or _isFsumTuple(x) else _2float(x=x)
     ps = tuple(_powers(p, n))
 
     if alts > 0:  # x**2, x**4, ...
@@ -752,7 +773,7 @@ def fremainder(x, y):
     # On Windows 32-bit with python 2.7, math.fmod(-0.0, 360)
     # == +0.0.  This fixes this bug.  See also Math::AngNormalize
     # in the C++ library, Math.sincosd has a similar fix.
-    if _isfinite(x):
+    if isfinite(x):
         try:
             r = remainder(x, y) if x else x
         except Exception as e:

pygeodesy/fstats.py

@@ -10,9 +10,10 @@ from __future__ import division as _; del _  # PYCHOK semicolon
 from pygeodesy.basics import isscalar, isodd, _xinstanceof, \
                             _xiterable, _xsubclassof, _zip
 from pygeodesy.constants import _0_0, _1_0, _2_0, _3_0, _4_0, _6_0
-from pygeodesy.errors import _AssertionError, _ValueError, _xError
+from pygeodesy.errors import _ValueError, _xError, _xkwds_item2
 from pygeodesy.fmath import Fsqrt,  Fmt
-from pygeodesy.fsums import _2finite, Fsum, _iadd_op_, _isFsumTuple
+from pygeodesy.fsums import _2finite, Fsum, _iadd_op_, \
+                            _isFsumTuple, _xsError
 from pygeodesy.interns import _odd_, _SPACE_
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY
 from pygeodesy.named import _name__, _Named, _NotImplemented, \
@@ -21,24 +22,19 @@ from pygeodesy.named import _name__, _Named, _NotImplemented, \
 # from pygeodesy.streprs import Fmt  # from .fmath
 
 __all__ = _ALL_LAZY.fstats
-__version__ = '24.05.21'
+__version__ = '24.09.23'
 
 
 def _2Floats(**xs):
-    '''(INTERNAL) Yield each value as C{float} or L{Fsum}.
+    '''(INTERNAL) Yield all C{xs} as C{float} or L{Fsum}.
     '''
+    name, xs = _xkwds_item2(xs)
     try:
-        name, xs = xs.popitem()
-    except Exception as X:
-        raise _AssertionError(xs=xs, cause=X)
-
-    try:
-        i = None
+        i, x = None, xs
         for i, x in enumerate(_xiterable(xs)):  # don't unravel Fsums
             yield x._Fsum if _isFsumTuple(x) else _2finite(x)
     except Exception as X:
-        raise _xError(X, name, xs) if i is None else \
-              _xError(X, Fmt.INDEX(name, i), x)
+        raise _xsError(X, xs, i, x, name)
 
 
 def _sampled(n, sample):