pygeodesy 24.4.18__py2.py3-none-any.whl → 24.5.2__py2.py3-none-any.whl

{PyGeodesy-24.4.18.dist-info → PyGeodesy-24.5.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PyGeodesy
-Version: 24.4.18
+Version: 24.5.2
 Summary: Pure Python geodesy tools
 Home-page: https://GitHub.com/mrJean1/PyGeodesy
 Author: Jean M. Brouwers
@@ -158,7 +158,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: April 18, 2024.*
+*Last updated: May 02, 2024.*
 
 License
 =======

{PyGeodesy-24.4.18.dist-info → PyGeodesy-24.5.2.dist-info}/RECORD

@@ -1,13 +1,13 @@
 pygeodesy/LICENSE,sha256=YfgAiyxOwY6P9Kkb1_5XN81nueTLrpb3Ffkv3EuPgFU,1144
-pygeodesy/__init__.py,sha256=Wx1Zo8hUxl6CV1eZZvLbrQSQ6VSB8mHBTmPqxcG0Mpc,40824
+pygeodesy/__init__.py,sha256=LsshXEtTCo0kRJI0vJfZzKu5Lo0Q0j5fiTf1bnCvefs,40733
 pygeodesy/__main__.py,sha256=qMFG3caY8ZXWu6uGiemzyT4OqTFZnsFtlxcGCAgkVJw,4637
 pygeodesy/albers.py,sha256=g2AVlpV8JO2AYFCthoIbRC2h1OJqjb9P3hpwF0C3TI8,30994
 pygeodesy/azimuthal.py,sha256=tc4JxgLi-0jzU08m4Bvi-t-kzHXYPeGuzL3j_tyVFUA,50125
-pygeodesy/basics.py,sha256=_BMYLzGKA6OIS3qv25HfHy7MIh0kAbmkzyScb59clUs,28160
+pygeodesy/basics.py,sha256=Mh70LwAH185gr8aYhdOhfqAa-SkHvAatlJEOFIH7N8Y,28162
 pygeodesy/booleans.py,sha256=HZbwoL-S7Ww9d4C2D1BVqXfmcuqqVpEVSU9_S0uyUyo,74204
 pygeodesy/cartesianBase.py,sha256=I3q29mRdBB3NCDmPoJsJ0QOFfLzkdMWc8X9zG4IwJyA,47264
 pygeodesy/clipy.py,sha256=VU3ynQ1IZ0v5hJlicqD48oW0imRgiL5_ZzRPrIjpfPw,27683
-pygeodesy/constants.py,sha256=ypxYsWB6tNAVgrD9QVpukJAPscaLgHtPLqPuot1gUfU,19174
+pygeodesy/constants.py,sha256=RxO8dMe_3AEfmZyKaRIy_44QnOTdEpzgYwimN0_w9Qs,19113
 pygeodesy/css.py,sha256=sKXsahUiyruDcUk-tGjA6mxq-xzwBoBKxKo9_b2uBmY,25394
 pygeodesy/datums.py,sha256=gJZPgV4bELZvZ8Sj2zE3MBysVtsLxqsN8zm0xjOKvpo,33851
 pygeodesy/dms.py,sha256=op3MU-59CoJQRdybnu21aVM9wtocd_-XFNAZFqmozSo,44439
@@ -25,11 +25,11 @@ pygeodesy/elliptic.py,sha256=XRVpmpdm3hbztM4P-qhWKQizbgsXDUxWujgmOrIpgiQ,42428
 pygeodesy/epsg.py,sha256=ldHoLWqJWR4FUiBVnDTtrI_e7TNjjtr9OkxDlri1g5E,8165
 pygeodesy/errors.py,sha256=mUaacHJKcSCowlcbzamaBMDqjI2dqsTM5O4ozIITtAk,27233
 pygeodesy/etm.py,sha256=joEhU2lw9ehRv101nRpZvXq8hOYAEzdcYFpif9_obfk,44585
-pygeodesy/fmath.py,sha256=9jjm_QxrQKNtOUxgNk02j-UVMFVjC1cDPaEqXznKgzk,30716
+pygeodesy/fmath.py,sha256=cRisVRkuarEwu4pKOCWmBa-98lV2PP9oJVkkODQLc2E,32480
 pygeodesy/formy.py,sha256=Rces4Q5ecED0WVfuI9mQzeyHE54LDowYfI6faBHpyeA,74536
 pygeodesy/frechet.py,sha256=qgee-ISBtT7Ov3rJkcd_t-WeXTbeNoMSQuMa0j3MyQc,33512
-pygeodesy/fstats.py,sha256=r8O2aIknHuoHsW8gZUIY42zpHZy7syHvCB7uCVtPpkc,25559
-pygeodesy/fsums.py,sha256=QjV9B-YGbhFALMPEQMKT8-1j7JFDWnabakAdZDvFsSs,72501
+pygeodesy/fstats.py,sha256=3IUmOJclp3XIzTLKVrwj24N0yCRvVMPSD2-BV7CY3_g,25559
+pygeodesy/fsums.py,sha256=xGA9nNeAHndWhs8pXJ4n46ZlNPIT2sn2f_6HmFOTnPs,76188
 pygeodesy/gars.py,sha256=fCiWBJ4kOJxPfNOxadX-OzBGDXj7C9g02NuGHiZa_88,11342
 pygeodesy/geodesicw.py,sha256=5yjJ2rLekSsjT7e-Km6v592ZcFlA41flQP_E42jU9Sw,26901
 pygeodesy/geodsolve.py,sha256=3GLI2_4gglzTpZwKchNzAv-XhWFBzWRKjcWsjCC3dvI,21900
@@ -37,12 +37,12 @@ pygeodesy/geohash.py,sha256=HwukecfkWvyfz9LHKfoA-WddQHYot_mB04o0lyWrKyU,32382
 pygeodesy/geoids.py,sha256=EwZJ50qRXeZb1591N99hnamM_EC2NP6MIhLb5swVHf4,80605
 pygeodesy/hausdorff.py,sha256=LJULGbDop4xqiWeRH5-Uayvk5WVCZQpxK4PCYu4kqSs,34192
 pygeodesy/heights.py,sha256=tjr-5eb8UXn2MAp6-157sb6qxPREFk5dtDjHtBgIgZ8,42528
-pygeodesy/interns.py,sha256=JqIAp_-TOzijDShKSmKIarz009vIXtqfB-JoD8ILvKo,30177
+pygeodesy/interns.py,sha256=ZeNGk2WIHLunBGSeOkU4Ihc9CNSfVDFZGNkbBQy2_8M,30120
 pygeodesy/iters.py,sha256=CfU365eE9F4oWZbBx0qTEvfodMhpa8rGOepLIv67xF8,20181
 pygeodesy/karney.py,sha256=lYdsSHis8xQJVS1V0F3cphx6Dhgd77xTlMOdVR8MN8I,34994
 pygeodesy/ktm.py,sha256=sCvbLvJItavlruilAyjeZ0sOZx2yJumzGe_UiIVbGi4,27315
 pygeodesy/latlonBase.py,sha256=LimCMj91J3Q4D_rTu7VN7c63cLL_KzxXagFQfedobmc,77730
-pygeodesy/lazily.py,sha256=4ZxBVj1Yu4Db-dw_e6oxXx7EtDSrgdGBV-ypUefMPGY,48633
+pygeodesy/lazily.py,sha256=YeOQQ2Tu0LTD-AMDABc0bemeLRKTteTEJXJDl3_7vZU,48645
 pygeodesy/lcc.py,sha256=yNpmAdOwv3HNh2ZLopF5QtvDokeQYCmnxSc8UFUPHO4,25737
 pygeodesy/ltp.py,sha256=8se8dYG7LG__Ru0FoK939DJBF9GCQGeJik_Un1Im8Xw,48609
 pygeodesy/ltpTuples.py,sha256=iKzvDfv33HfmEevjNISUS4c9KXIkiISEkVlQopyKgfY,55769
@@ -59,9 +59,9 @@ pygeodesy/solveBase.py,sha256=cBe3QQWEaJlPWUlXoiYo1eFMh4z6GDZbQyt8ZiwMfos,16687
 pygeodesy/sphericalBase.py,sha256=sH2WwmY_RL0STyr1OhtCDKvGhO0Ks7o3V-MHLddpiHQ,31944
 pygeodesy/sphericalNvector.py,sha256=TBs_nid_DTcox7bBegh_Djv4zZRGRY4mH6dd1xIHGms,56963
 pygeodesy/sphericalTrigonometry.py,sha256=Q9qsqzIK43xBken2G8uP9snPXPuBoZ20pMHqbT_IjG0,64216
-pygeodesy/streprs.py,sha256=-29aUO7dYe-YWfTjflljesxGUp7IfEAjk5dJLHWMPVs,23413
+pygeodesy/streprs.py,sha256=NIiLVzXLnJR7w5fNTkivpbQN5XdZYJIc189Iv_bjpYw,23513
 pygeodesy/trf.py,sha256=4gVtF780EEZEqWwMBDe4x1AEkIz3hB8h4l0_0NXWM5U,118903
-pygeodesy/triaxials.py,sha256=fkUcWkXkykZNbTq62LcdnEjhNu9q-e2UaSpRotdzuu8,61546
+pygeodesy/triaxials.py,sha256=nYFQcLimivcxpeXMj9a83iWuENDU1gI-2zNgCkZzL0U,61549
 pygeodesy/units.py,sha256=ofFTYc7mF9wOpfvpYHN83dGszMM18rTdsxfEXVzaIko,38705
 pygeodesy/unitsBase.py,sha256=xRZlNlO7S47Mwy0FGBg_cefvB4fEtWFM5dAE-Nlw-NE,13111
 pygeodesy/ups.py,sha256=07AG-Rb1D4EvmN3g-xk3ZxKzyaNYvOzztcyYbx7t8jc,23229
@@ -79,7 +79,7 @@ pygeodesy/auxilats/_CX_6.py,sha256=x4hrHxunjqqkEsZbVW1dyn3R5mR8EWYOuX77RzDEyh4,1
 pygeodesy/auxilats/_CX_8.py,sha256=t-OO48t4cCOsOgMi1FoyPw2DQqybVZFsgQMqAYP-jKY,25304
 pygeodesy/auxilats/__init__.py,sha256=yOa5RDYOOed9Z4dvcKlatH0JBgohAJ8ADXjU9mZGTfQ,2881
 pygeodesy/auxilats/__main__.py,sha256=TVN8evCKEzBKdANU5-Ucps8bZySqBTQmEGU4hy_U1cw,2980
-pygeodesy/auxilats/auxAngle.py,sha256=3tTYEuiX2YmEOmoyqGPYrSN3gLXbcKWKFH1TNt3NqbI,16738
+pygeodesy/auxilats/auxAngle.py,sha256=LEa1GVR1EOYnTxR-mBLUj4mGoLO6jRrpsyG4v7N8W_0,16727
 pygeodesy/auxilats/auxDLat.py,sha256=U3BHF6AkW5--nrwH9whL5wZVIZLH_keRLfUz6c6KBx4,10950
 pygeodesy/auxilats/auxDST.py,sha256=bdNEh3TdABWH7qW6ngZx5SnzDM3B5rtPLvYhvCnFsKc,10497
 pygeodesy/auxilats/auxLat.py,sha256=U6thEwMEoUm4OFb2_JviEtEn2-nqyVjNbII7oyMf0hM,31282
@@ -109,7 +109,7 @@ pygeodesy/rhumb/aux_.py,sha256=W4HkgoHAQz_aWaysfsmOA5010nCGvBfUs2Q-bymnHYU,16660
 pygeodesy/rhumb/bases.py,sha256=kzU_Dgt4FNPMgTg5rqbw-HiNpflDPKpmq9jhPbQmR4U,53851
 pygeodesy/rhumb/ekx.py,sha256=lF3tZ-ZY9fPJV8y1kgHW-7EOZCyb3gJr-kR-jj5Fbf8,23871
 pygeodesy/rhumb/solve.py,sha256=NZfwx7xv5UriQs7A0c7ZhoxxVUeAT15UwXK_jtwEXQw,27802
-PyGeodesy-24.4.18.dist-info/METADATA,sha256=kwVbSazBaM7cCqqiNZ0-SQbbnG7-Y3vLtmAh_yGMK3s,19366
-PyGeodesy-24.4.18.dist-info/WHEEL,sha256=DZajD4pwLWue70CAfc7YaxT1wLUciNBvN_TTcvXpltE,110
-PyGeodesy-24.4.18.dist-info/top_level.txt,sha256=cEQPatCXzKZqrivpULC5V5fuy9_V_bAwaP_gUGid7pQ,10
-PyGeodesy-24.4.18.dist-info/RECORD,,
+PyGeodesy-24.5.2.dist-info/METADATA,sha256=CCvKPsW0UUmlWkrNk0W6dlZb8J--2lIS-wHNpMJGPwk,19363
+PyGeodesy-24.5.2.dist-info/WHEEL,sha256=DZajD4pwLWue70CAfc7YaxT1wLUciNBvN_TTcvXpltE,110
+PyGeodesy-24.5.2.dist-info/top_level.txt,sha256=cEQPatCXzKZqrivpULC5V5fuy9_V_bAwaP_gUGid7pQ,10
+PyGeodesy-24.5.2.dist-info/RECORD,,

pygeodesy/__init__.py

@@ -194,7 +194,6 @@ The following environment variables are observed by C{PyGeodesy}:
 
  - C{PYGEODESY_EXCEPTION_CHAINING} - see module L{pygeodesy.errors}.
  - C{PYGEODESY_FMT_FORM} - see module L{pygeodesy.dms}.
- - C{PYGEODESY_FSUM_PARTIALS} - see module L{pygeodesy.fsums} and class L{pygeodesy.Fsum}.
  - C{PYGEODESY_FSUM_RESIDUAL} - see module L{pygeodesy.fsums} and class L{pygeodesy.Fsum}.
  - C{PYGEODESY_GEOCONVERT} - see module L{pygeodesy.mgrs}.
  - C{PYGEODESY_GEODSOLVE} - see module L{pygeodesy.geodsolve}.
@@ -586,7 +585,7 @@ else:
     _init__all__ = False
 
 from pygeodesy.interns import _DOT_, _version2  # PYCHOK import
-__version__ = '24.04.18'
+__version__ = '24.05.02'
 # see setup.py for similar logic
 version     = _DOT_(*_version2(__version__, n=3))
 

pygeodesy/auxilats/auxAngle.py

@@ -20,8 +20,8 @@ from pygeodesy.constants import EPS, _INF_NAN_NINF, MAX, NAN, _0_0, _0_5, _1_0,
                                _copysign_1_0, _over, _pos_self, isfinite, isnan
 # from pygeodesy.errors import AuxError, _xkwds_get  # from .auxily, .basics
 from pygeodesy.fmath import hypot,  unstr
-from pygeodesy.fsums import _add_op_, _isub_op_, _sub_op_,  _iadd_op_, _Named, NN
-# from pygeodesy.interns import NN, _iadd_op_  # from .fsums
+from pygeodesy.fsums import _add_op_, _iadd_op_, _isub_op_, _sub_op_,  _Named, NN
+# from pygeodesy.interns import NN  # from .fsums
 # from pygeodesy.named import _Named  # from .fsums
 from pygeodesy.lazily import _ALL_DOCS, _ALL_MODS as _MODS
 from pygeodesy.props import Property, Property_RO, property_RO, _update_all
@@ -32,7 +32,7 @@ from pygeodesy.utily import atan2d, sincos2, sincos2d
 from math import asinh, atan2, copysign, degrees, fabs, radians, sinh
 
 __all__ = ()
-__version__ = '23.12.02'
+__version__ = '24.04.26'
 
 _0_INF_NAN_NINF = (0, _0_0) + _INF_NAN_NINF
 _MAX_2          =  MAX * _0_5  # PYCHOK used!

pygeodesy/basics.py

@@ -34,7 +34,7 @@ from math import copysign as _copysign
 import inspect as _inspect
 
 __all__ = _ALL_LAZY.basics
-__version__ = '24.04.04'
+__version__ = '24.04.28'
 
 _0_0                  =  0.0  # in .constants
 _below_               = 'below'
@@ -576,7 +576,7 @@ def signBit(x):
 def _signOf(x, ref):  # in .fsums
     '''(INTERNAL) Return the sign of B{C{x}} versus B{C{ref}}.
     '''
-    return +1 if x > ref else (-1 if x < ref else 0)
+    return (-1) if x < ref else (+1 if x > ref else 0)
 
 
 def signOf(x):

pygeodesy/constants.py

@@ -23,7 +23,7 @@ except ImportError:  # Python 2-
     _inf, _nan = float(_INF_), float(_NAN_)
 
 __all__ = _ALL_LAZY.constants
-__version__ = '24.04.04'
+__version__ = '24.04.23'
 
 
 def _copysign_0_0(y):
@@ -247,7 +247,6 @@ _1_16th  = _Float(_1_16th =_1_0 / _16_0)  # PYCHOK in .ellipsoids, .karney
 _1_64th  = _Float(_1_64th =_1_0 /  64)    # PYCHOK in .elliptic, pow(2.0, -6)
 _1_3rd   = _Float(_1_3rd  =_1_0 /  _3_0)  # PYCHOK in .fmath
 _1_6th   = _Float(_1_6th  =_1_0 /  _6_0)  # PYCHOK in .fmath
-_2_3rd   = _Float(_2_3rd  =_2_0 /  _3_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>

pygeodesy/fmath.py

@@ -6,17 +6,14 @@ u'''Utilities using precision floating point summation.
 # make sure int/int division yields float quotient, see .basics
 from __future__ import division as _; del _  # PYCHOK semicolon
 
-from pygeodesy.basics import _copysign, copysign0, isint, len2
+from pygeodesy.basics import _copysign, copysign0, isbool, isint, isscalar, len2
 from pygeodesy.constants import EPS0, EPS02, EPS1, NAN, PI, PI_2, PI_4, \
-                               _0_0, _0_125, _0_25, _0_5, _1_0, _1_3rd, \
-                               _1_5, _1_6th, _2_0, _2_3rd, _3_0, \
-                               _copysign_0_0, _isfinite, _over, remainder
+                               _0_0, _0_125, _1_6th, _0_25, _1_3rd, _0_5, _1_0, \
+                               _1_5, _copysign_0_0, _isfinite, _over, remainder
 from pygeodesy.errors import _IsnotError, LenError, _TypeError, _ValueError, \
                              _xError, _xkwds_get, _xkwds_pop2
-from pygeodesy.fsums import _2float, Fsum, _fsum, fsum, fsum1_, _1primed, \
-                             Fmt, unstr
-from pygeodesy.interns import MISSING, _few_, _h_, _invokation_, _negative_, \
-                             _not_scalar_, _SPACE_, _too_
+from pygeodesy.fsums import _2float, Fsum, fsum, fsum1_, _1primed,  Fmt, unstr
+from pygeodesy.interns import MISSING, _few_, _negative_, _not_scalar_, _too_
 from pygeodesy.lazily import _ALL_LAZY, _sys_version_info2
 # from pygeodesy.streprs import Fmt, unstr  # from .fsums
 from pygeodesy.units import Int_, _isHeight, _isRadius,  Float_  # PYCHOK for .heights
@@ -25,21 +22,14 @@ from math import fabs, sqrt  # pow
 import operator as _operator  # in .datums, .trf, .utm
 
 __all__ = _ALL_LAZY.fmath
-__version__ = '24.04.17'
+__version__ = '24.04.24'
 
 # sqrt(2) <https://WikiPedia.org/wiki/Square_root_of_2>
 _0_4142  =  0.41421356237309504880  # ... sqrt(2) - 1
+_2_3rd   = _1_3rd * 2
 _h_lt_b_ = 'abs(h) < abs(b)'
 
 
-def _Fsum__init__(inst, raiser=MISSING, **name_RESIDUAL):
-    '''(INTERNAL) Init an C{Fsum} instance.
-    '''
-    Fsum.__init__(inst, **name_RESIDUAL)  # PYCHOK self
-    inst._fset_ps(_0_0)
-    return {} if raiser is MISSING else dict(raiser=raiser)
-
-
 class Fdot(Fsum):
     '''Precision dot product.
     '''
@@ -104,7 +94,7 @@ class Fhorner(Fsum):
 
 
 class Fhypot(Fsum):
-    '''Precision summation and hypotenuse, default C{power=2}.
+    '''Precision summation and hypotenuse, default C{root=2}.
     '''
     def __init__(self, *xs, **root_name_RESIDUAL_raiser):
         '''New L{Fhypot} hypotenuse of (the I{root} of) several components.
@@ -189,7 +179,7 @@ class Froot(Fsum):
         try:
             raiser = _Fsum__init__(self, **name_RESIDUAL_raiser)
             if xs:
-                self. fadd(xs)
+                self.fadd(xs)
             self._fset(self.root(root, **raiser))
         except Exception as X:
             raise self._ErrorXs(X, xs, root=root)
@@ -203,7 +193,7 @@ class Fcbrt(Froot):
 
            @see: Class L{Froot} for further details.
         '''
-        Froot.__init__(self, _3_0, *xs, **name_RESIDUAL_raiser)
+        Froot.__init__(self, 3, *xs, **name_RESIDUAL_raiser)
 
 
 class Fsqrt(Froot):
@@ -214,15 +204,24 @@ class Fsqrt(Froot):
 
            @see: Class L{Froot} for further details.
         '''
-        Froot.__init__(self, _2_0, *xs, **name_RESIDUAL_raiser)
+        Froot.__init__(self, 2, *xs, **name_RESIDUAL_raiser)
+
+
+def _Fsum__init__(inst, raiser=MISSING, **name_RESIDUAL):
+    '''(INTERNAL) Init an C{F...} instance above.
+    '''
+    Fsum.__init__(inst, **name_RESIDUAL)  # PYCHOK self
+    inst._fset_ps(_0_0)
+    return {} if raiser is MISSING else dict(raiser=raiser)
 
 
 def bqrt(x):
-    '''Return the 4-th, I{bi-quadratic} or I{quartic} root, M{x**(1 / 4)}.
+    '''Return the 4-th, I{bi-quadratic} or I{quartic} root, M{x**(1 / 4)},
+       preserving C{type(B{x})}.
 
-       @arg x: Value (C{scalar}).
+       @arg x: Value (C{scalar} or L{Fsum} instance).
 
-       @return: I{Quartic} root (C{float}).
+       @return: I{Quartic} root (C{float} or L{Fsum}).
 
        @raise ValueError: Negative B{C{x}}.
 
@@ -232,70 +231,76 @@ def bqrt(x):
 
 
 try:
-    from math import cbrt  # Python 3.11+
-
-    def cbrt2(x):
-        '''Compute the cube root I{squared} M{x**(2/3)}.
-        '''
-        return cbrt(x)**2  # cbrt(-0.0*2) == -0.0
+    from math import cbrt as _cbrt  # Python 3.11+
 
 except ImportError:  # Python 3.10-
 
-    def cbrt(x):
-        '''Compute the cube root M{x**(1/3)}.
-
-           @arg x: Value (C{scalar}).
-
-           @return: Cubic root (C{float}).
-
-           @see: Functions L{cbrt2} and L{sqrt3}.
+    def _cbrt(x):
+        '''(INTERNAL) Compute the I{signed}, cube root M{x**(1/3)}.
         '''
         # <https://archive.lib.MSU.edu/crcmath/math/math/r/r021.htm>
         # simpler and more accurate than Ken Turkowski's CubeRoot, see
         # <https://People.FreeBSD.org/~lstewart/references/apple_tr_kt32_cuberoot.pdf>
-        return _copysign(pow(fabs(x), _1_3rd), x)  # cbrt(-0.0) == -0.0
+        return _copysign(pow(fabs(x), _1_3rd), x)  # to avoid complex
 
-    def cbrt2(x):  # PYCHOK attr
-        '''Compute the cube root I{squared} M{x**(2/3)}.
 
-           @arg x: Value (C{scalar}).
+def cbrt(x):
+    '''Compute the cube root M{x**(1/3)}, preserving C{type(B{x})}.
 
-           @return: Cube root I{squared} (C{float}).
+       @arg x: Value (C{scalar} or L{Fsum} instance).
 
-           @see: Functions L{cbrt} and L{sqrt3}.
-        '''
-        return pow(fabs(x), _2_3rd)  # XXX pow(fabs(x), _1_3rd)**2
+       @return: Cubic root (C{float} or L{Fsum}).
+
+       @see: Functions L{cbrt2} and L{sqrt3}.
+    '''
+    if isinstance(x, Fsum):
+        r = (-(-x).pow(_1_3rd)) if x < 0 else x.pow(_1_3rd)
+    else:
+        r = _cbrt(x)
+    return r  # cbrt(-0.0) == -0.0
+
+
+def cbrt2(x):  # PYCHOK attr
+    '''Compute the cube root I{squared} M{x**(2/3)}, preserving C{type(B{x})}.
+
+       @arg x: Value (C{scalar} or L{Fsum} instance).
+
+       @return: Cube root I{squared} (C{float} or L{Fsum}).
+
+       @see: Functions L{cbrt} and L{sqrt3}.
+    '''
+    return abs(x).pow(_2_3rd) if isinstance(x, Fsum) else _cbrt(x**2)
 
 
 def euclid(x, y):
     '''I{Appoximate} the norm M{sqrt(x**2 + y**2)} by
        M{max(abs(x), abs(y)) + min(abs(x), abs(y)) * 0.4142...}.
 
-       @arg x: X component (C{scalar}).
-       @arg y: Y component (C{scalar}).
+       @arg x: X component (C{scalar} or L{Fsum} instance).
+       @arg y: Y component (C{scalar} or L{Fsum} instance).
 
-       @return: Appoximate norm (C{float}).
+       @return: Appoximate norm (C{float} or L{Fsum}).
 
        @see: Function L{euclid_}.
     '''
-    x, y = fabs(x), fabs(y)
+    x, y = abs(x), abs(y)  # NOT fabs!
     if x < y:
         x, y = y, x
     return x + y * _0_4142  # XXX * _0_5 before 20.10.02
 
 
 def euclid_(*xs):
-    '''I{Appoximate} the norm M{sqrt(sum(x**2 for x in xs))}
-       by cascaded L{euclid}.
+    '''I{Appoximate} the norm M{sqrt(sum(x**2 for x in xs))} by
+       cascaded L{euclid}.
 
-       @arg xs: X arguments, positional (C{scalar}s).
+       @arg xs: X arguments, positional (C{scalar}s or L{Fsum} instances).
 
-       @return: Appoximate norm (C{float}).
+       @return: Appoximate norm (C{float} or L{Fsum}).
 
        @see: Function L{euclid}.
     '''
     e = _0_0
-    for x in sorted(map(fabs, xs)):  # XXX not reverse=True
+    for x in sorted(map(abs, xs)):  # NOT fabs, reverse=True!
         # e = euclid(x, e)
         if e < x:
             e, x = x, e
@@ -388,8 +393,8 @@ def fatan2(y, x):
 def favg(v1, v2, f=_0_5):
     '''Return the average of two values.
 
-       @arg v1: One value (C{scalar}).
-       @arg v2: Other value (C{scalar}).
+       @arg v1: One value (C{scalar} or L{Fsum} instance).
+       @arg v2: Other value (C{scalar} or L{Fsum} instance).
        @kwarg f: Optional fraction (C{float}).
 
        @return: M{v1 + f * (v2 - v1)} (C{float}).
@@ -606,6 +611,39 @@ except ImportError:
         return freduce(_operator.mul, xs, start)
 
 
+def frandoms(n, seeded=None):
+    '''Generate C{n} (long) lists of random C{floats}.
+
+       @arg n: Number of lists to generate (C{int}, non-negative).
+       @kwarg seeded: If C{scalar}, use C{random.seed(B{seeded})} or
+                      if C{True}, seed using today's C{year-day}.
+
+       @see: U{Hettinger<https://GitHub.com/ActiveState/code/tree/master/recipes/
+             Python/393090_Binary_floating_point_summatiaccurate_full/recipe-393090.py>}.
+    '''
+    from random import gauss, random, seed, shuffle
+
+    if seeded is None:
+        pass
+    elif seeded and isbool(seeded):
+        from time import localtime
+        seed(localtime().tm_yday)
+    elif isscalar(seeded):
+        seed(seeded)
+
+    c = (7, 1e100, -7, -1e100, -9e-20, 8e-20) * 7
+    for _ in range(n):
+        s  = 0
+        t  = list(c)
+        _a = t.append
+        for _ in range(n * 8):
+            v = gauss(0, random())**7 - s
+            _a(v)
+            s += v
+        shuffle(t)
+        yield t
+
+
 def frange(start, number, step=1):
     '''Generate a range of C{float}s.
 
@@ -715,8 +753,8 @@ if _sys_version_info2 < (3, 8):  # PYCHOK no cover
                   computed as M{hypot_(*((c1 - c2) for c1, c2 in zip(p1, p2)))},
                   provided I{p1} and I{p2} have the same, non-zero length I{n}.
         '''
-        h, x2 = _h_x2(xs, hypot_)
-        return (h * sqrt(x2)) if x2 else _0_0
+        _, R = _h_xs2(xs, True, hypot_)
+        return float(R)
 
 elif _sys_version_info2 < (3, 10):
     # In Python 3.8 and 3.9 C{math.hypot} is inaccurate, see
@@ -732,13 +770,8 @@ elif _sys_version_info2 < (3, 10):
 
            @return: C{sqrt(B{x}**2 + B{y}**2)} (C{float}).
         '''
-        if x:
-            h = sqrt(x**2 + y**2) if y else fabs(x)
-        elif y:
-            h = fabs(y)
-        else:
-            h = _0_0
-        return h
+        return (float(Fhypot(x, y, raiser=False)) if y else
+                fabs(x)) if x else fabs(y)
 
     from math import hypot as hypot_  # PYCHOK in Python 3.8 and 3.9
 else:
@@ -746,29 +779,32 @@ else:
     hypot_ = hypot
 
 
-def _h_x2(xs, which):
+def _h_xs2(xs, pot_, which):
     '''(INTERNAL) Helper for L{hypot_} and L{hypot2_}.
     '''
     n, xs = len2(xs)
     if n > 0:
-        h = float(max(map(fabs, xs)))
+        h = float(max(map(abs, xs)))  # NOT fabs!
         if h < EPS0:
-            x2 = _0_0
+            R = _0_0
         elif n > 1:
-            _h = (_1_0 / h) if h != _1_0 else _1_0
-            x2 = _fsum(_1primed((x * _h)**2 for x in xs))
+            if pot_:
+                if h != _1_0:
+                    xs = ((x / h) for x in xs)
+                R = Fhypot(*xs, raiser=False) * h
+            else:
+                R = Fpowers(2, *xs)
         else:
-            x2 = _1_0
-        return h, x2
+            R = h if pot_ else (h**2)
+        return h, R
 
-    t = Fmt.PAREN(which.__name__, xs)
-    raise _ValueError(t, txt=_too_(_few_))
+    raise _ValueError(unstr(which, xs), txt=_too_(_few_))
 
 
 def hypot1(x):
     '''Compute the norm M{sqrt(1 + x**2)}.
 
-       @arg x: Argument (C{scalar}).
+       @arg x: Argument (C{scalar} or L{Fsum} instance).
 
        @return: Norm (C{float}).
     '''
@@ -778,12 +814,12 @@ def hypot1(x):
 def hypot2(x, y):
     '''Compute the I{squared} norm M{x**2 + y**2}.
 
-       @arg x: X argument (C{scalar}).
-       @arg y: Y argument (C{scalar}).
+       @arg x: X argument (C{scalar} or L{Fsum} instance).
+       @arg y: Y argument (C{scalar} or L{Fsum} instance).
 
-       @return: C{B{x}**2 + B{y}**2} (C{float}).
+       @return: C{B{x}**2 + B{y}**2} (C{float} or L{Fsum}).
     '''
-    if fabs(x) < fabs(y):
+    if abs(x) < abs(y):  # NOT fabs!
         x, y = y, x
     if x:
         h2 = x**2
@@ -795,9 +831,10 @@ def hypot2(x, y):
 
 
 def hypot2_(*xs):
-    '''Compute the I{squared} norm C{sum(x**2 for x in B{xs})}.
+    '''Compute the I{squared} norm C{fsum(x**2 for x in B{xs})}.
 
-       @arg xs: X arguments (C{scalar}s), all positional.
+       @arg xs: X arguments (C{scalar}s or L{Fsum} instances),
+                all positional.
 
        @return: Squared norm (C{float}).
 
@@ -807,8 +844,8 @@ def hypot2_(*xs):
 
        @see: Function L{hypot_}.
     '''
-    h, x2 = _h_x2(xs, hypot2_)
-    return (h**2 * x2) if x2 else _0_0
+    _, R = _h_xs2(xs, False, hypot2_)
+    return float(R)
 
 
 def _map_mul(a, b, where):
@@ -860,12 +897,13 @@ def norm_(*xs):
               or zero norm.
     '''
     try:
+        i  = x = h = None
         h  = hypot_(*xs)
         _h = (_1_0 / h) if h else _0_0
         for i, x in enumerate(xs):
             yield x * _h
-    except Exception as e:
-        raise _xError(e, Fmt.SQUARE(xs=i), x, _h_, h)
+    except Exception as X:
+        raise _xError(X, Fmt.SQUARE(xs=i), x, h=h)
 
 
 def _powers(x, n):
@@ -880,34 +918,40 @@ def _powers(x, n):
 def _root(x, p, where):
     '''(INTERNAL) Raise C{x} to power C{0 < p < 1}.
     '''
-    if x < 0:
-        t = _SPACE_(_invokation_, where.__name__)
-        raise _ValueError(unstr(t, x), txt=_negative_)
-    return pow(x, p) if x else _0_0
+    try:
+        if x > 0:
+            return pow(x, p)
+        elif x < 0:
+            raise ValueError(_negative_)
+    except Exception as X:
+        raise _xError(X, unstr(where, x))
+    return _0_0
 
 
 def sqrt0(x, Error=None):
-    '''Return the square root iff C{B{x} >} L{EPS02}.
+    '''Return the square root C{sqrt(B{x})} iff C{B{x} > }L{EPS02},
+       preserving C{type(B{x})}.
 
-       @arg x: Value (C{scalar}).
+       @arg x: Value (C{scalar} or L{Fsum} instance).
        @kwarg Error: Error to raise for negative B{C{x}}.
 
-       @return: Square root (C{float}) or C{0.0}.
+       @return: Square root (C{float} or L{Fsum}) or C{0.0}.
 
        @note: Any C{B{x} < }L{EPS02} I{including} C{B{x} < 0}
               returns C{0.0}.
     '''
     if Error and x < 0:
-        raise Error(Fmt.PAREN(sqrt=x))
-    return sqrt(x) if x > EPS02 else (_0_0 if x < EPS02 else EPS0)
+        raise Error(unstr(sqrt0, x))
+    return _root(x, _0_5, sqrt0) if x > EPS02 else (_0_0 if x < EPS02 else EPS0)
 
 
 def sqrt3(x):
-    '''Return the square root, I{cubed} M{sqrt(x)**3} or M{sqrt(x**3)}.
+    '''Return the square root, I{cubed} M{sqrt(x)**3} or M{sqrt(x**3)},
+       preserving C{type(B{x})}.
 
-       @arg x: Value (C{scalar}).
+       @arg x: Value (C{scalar} or L{Fsum} instance).
 
-       @return: Square root I{cubed} (C{float}).
+       @return: Square root I{cubed} (C{float} or L{Fsum}).
 
        @raise ValueError: Negative B{C{x}}.
 
@@ -956,11 +1000,12 @@ def sqrt_a(h, b):
 
 
 def zcrt(x):
-    '''Return the 6-th, I{zenzi-cubic} root, M{x**(1 / 6)}.
+    '''Return the 6-th, I{zenzi-cubic} root, M{x**(1 / 6)},
+       preserving C{type(B{x})}.
 
-       @arg x: Value (C{scalar}).
+       @arg x: Value (C{scalar} or L{Fsum} instance).
 
-       @return: I{Zenzi-cubic} root (C{float}).
+       @return: I{Zenzi-cubic} root (C{float} or L{Fsum}).
 
        @see: Functions L{bqrt} and L{zqrt}.
 
@@ -970,11 +1015,12 @@ def zcrt(x):
 
 
 def zqrt(x):
-    '''Return the 8-th, I{zenzi-quartic} or I{squared-quartic} root, M{x**(1 / 8)}.
+    '''Return the 8-th, I{zenzi-quartic} or I{squared-quartic} root,
+       M{x**(1 / 8)}, preserving C{type(B{x})}.
 
-       @arg x: Value (C{scalar}).
+       @arg x: Value (C{scalar} or L{Fsum} instance).
 
-       @return: I{Zenzi-quartic} root (C{float}).
+       @return: I{Zenzi-quartic} root (C{float} or L{Fsum}).
 
        @see: Functions L{bqrt} and L{zcrt}.
 

pygeodesy/fstats.py

@@ -12,8 +12,8 @@ from pygeodesy.basics import isodd, islistuple, _xinstanceof, \
 from pygeodesy.constants import _0_0, _2_0, _3_0, _4_0, _6_0,  _xError
 # from pygeodesy.errors import _xError  # from .constants
 from pygeodesy.fmath import hypot2,  sqrt
-from pygeodesy.fsums import _2float, Fsum,  Fmt
-from pygeodesy.interns import NN, _iadd_op_, _invalid_, _other_, _SPACE_
+from pygeodesy.fsums import _2float, Fsum, _iadd_op_,  Fmt
+from pygeodesy.interns import NN, _invalid_, _other_, _SPACE_
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY
 from pygeodesy.named import _Named, _NotImplemented, property_RO
 # from pygeodesy.props import property_RO  # from .named
@@ -22,7 +22,7 @@ from pygeodesy.named import _Named, _NotImplemented, property_RO
 # from math import sqrt  # from .fmath
 
 __all__ = _ALL_LAZY.fstats
-__version__ = '24.04.07'
+__version__ = '24.04.26'
 
 _Floats =  Fsum, float
 _Scalar = _Floats + (int,)  # XXX basics._Ints is ABCMeta in Python 2