pygeodesy 24.4.12__py2.py3-none-any.whl → 24.4.24__py2.py3-none-any.whl

{PyGeodesy-24.4.12.dist-info → PyGeodesy-24.4.24.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PyGeodesy
-Version: 24.4.12
+Version: 24.4.24
 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 12, 2024.*
+*Last updated: April 24, 2024.*
 
 License
 =======

{PyGeodesy-24.4.12.dist-info → PyGeodesy-24.4.24.dist-info}/RECORD

@@ -1,13 +1,13 @@
 pygeodesy/LICENSE,sha256=YfgAiyxOwY6P9Kkb1_5XN81nueTLrpb3Ffkv3EuPgFU,1144
-pygeodesy/__init__.py,sha256=VQ4Ax0S1ZenUk2bZtcfv8F-OWwKawIzU85WpPHWDrRE,40824
+pygeodesy/__init__.py,sha256=td1SZErPJBK8oIAHrj4sVvK-n3PjnJbckJQjWhyNHf0,40824
 pygeodesy/__main__.py,sha256=qMFG3caY8ZXWu6uGiemzyT4OqTFZnsFtlxcGCAgkVJw,4637
-pygeodesy/albers.py,sha256=mx1Bu05G52WuDaZdEHSpnN9bmRagem4vsbjYodQeuiM,31102
+pygeodesy/albers.py,sha256=g2AVlpV8JO2AYFCthoIbRC2h1OJqjb9P3hpwF0C3TI8,30994
 pygeodesy/azimuthal.py,sha256=tc4JxgLi-0jzU08m4Bvi-t-kzHXYPeGuzL3j_tyVFUA,50125
 pygeodesy/basics.py,sha256=_BMYLzGKA6OIS3qv25HfHy7MIh0kAbmkzyScb59clUs,28160
 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
@@ -20,16 +20,16 @@ pygeodesy/ellipsoidalGeodSolve.py,sha256=xfNcb70AuPbIiJUqhisIfg6pCQBqXyQffBNuESE
 pygeodesy/ellipsoidalKarney.py,sha256=FoOsfNf96oHuGa9k4AB42y_UUBuOkh4xMW7PFGF63SA,19246
 pygeodesy/ellipsoidalNvector.py,sha256=vHVKSfR86LnIapGlwl5_uCaEuQGdm0e89mwE6Ay-394,30323
 pygeodesy/ellipsoidalVincenty.py,sha256=Uk93B8GFcMnWiwKUOA84TlLwbBlvzVnsAXAx0cU7cn4,26004
-pygeodesy/ellipsoids.py,sha256=cykBbfMMDGWk65rfGlYGm1s_Wuh_NZD1_9bgKelPJ5s,119945
-pygeodesy/elliptic.py,sha256=J9n6YdKeGQbcwgDKXByGRKMdAyDt1V8xKCPE_3RsGWM,42419
+pygeodesy/ellipsoids.py,sha256=SxMsivNnE9ibUjQGd0d24H0SzoSjh6jfhLWjasOWpao,119946
+pygeodesy/elliptic.py,sha256=XRVpmpdm3hbztM4P-qhWKQizbgsXDUxWujgmOrIpgiQ,42428
 pygeodesy/epsg.py,sha256=ldHoLWqJWR4FUiBVnDTtrI_e7TNjjtr9OkxDlri1g5E,8165
-pygeodesy/errors.py,sha256=fRGrdnb81Y9KB4Y1ua1nvB9zuuiv8IL6c6YN_cb3kBo,27280
+pygeodesy/errors.py,sha256=mUaacHJKcSCowlcbzamaBMDqjI2dqsTM5O4ozIITtAk,27233
 pygeodesy/etm.py,sha256=joEhU2lw9ehRv101nRpZvXq8hOYAEzdcYFpif9_obfk,44585
-pygeodesy/fmath.py,sha256=cmObX2axb8BfPyhaSDGFFVRs63kBP4ZWhWzGzPhl6OQ,30286
+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=9ekLTz59nSON-PHtgQp_m6FxnYj3UOUCagDtoT7h8AE,69283
+pygeodesy/fsums.py,sha256=Am8CRdltomItdCpdBUDyc1xP53KmrCXXs9rBTa19hs4,73370
 pygeodesy/gars.py,sha256=fCiWBJ4kOJxPfNOxadX-OzBGDXj7C9g02NuGHiZa_88,11342
 pygeodesy/geodesicw.py,sha256=5yjJ2rLekSsjT7e-Km6v592ZcFlA41flQP_E42jU9Sw,26901
 pygeodesy/geodsolve.py,sha256=3GLI2_4gglzTpZwKchNzAv-XhWFBzWRKjcWsjCC3dvI,21900
@@ -42,18 +42,18 @@ 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=pLyWvdPoIdAQ8mQXoiGVNSkjBDew0T1Nkvv6OWj8Dv0,48624
+pygeodesy/lazily.py,sha256=YeOQQ2Tu0LTD-AMDABc0bemeLRKTteTEJXJDl3_7vZU,48645
 pygeodesy/lcc.py,sha256=yNpmAdOwv3HNh2ZLopF5QtvDokeQYCmnxSc8UFUPHO4,25737
-pygeodesy/ltp.py,sha256=AEzriv6LyqjJwROF_6kLWjsu0_J_cjVhNMNlsjr4ac4,48579
+pygeodesy/ltp.py,sha256=8se8dYG7LG__Ru0FoK939DJBF9GCQGeJik_Un1Im8Xw,48609
 pygeodesy/ltpTuples.py,sha256=iKzvDfv33HfmEevjNISUS4c9KXIkiISEkVlQopyKgfY,55769
 pygeodesy/mgrs.py,sha256=GJ-vG54E2zSLtROyqJ_kj8fHH-8jWGf3kVe7gb8993o,28978
 pygeodesy/named.py,sha256=2Kg4tv883V80ZVUEgFzDeBmIEMkig_VUbJuufC0lhN0,48009
 pygeodesy/namedTuples.py,sha256=msEYvQKj4oCW35Kx80RpeQXY85A1G0O3GEqsysLaV2M,27288
 pygeodesy/nvectorBase.py,sha256=WJY5ylZ_HCeM1wGL7ryEecJphEYibp-EAWL0piKlMG8,28766
-pygeodesy/osgr.py,sha256=5JB6b_nvyZ6rFMlC7ZB45mvYtjc_aiLb4DPdeaui-KY,30995
+pygeodesy/osgr.py,sha256=voC1rbOe4dnS_EjaDaRDRUQ4MoG7puMBEMV6lcAuIBI,31029
 pygeodesy/points.py,sha256=QNAAF0r7iqSXFr_6aQm-wW7f6vU3KusWVhXEMVeVo58,64442
 pygeodesy/props.py,sha256=ZnXRuSOkfkeJY3Of-iT06FNWZkTYgLZQBPBNF6aS0jE,21421
-pygeodesy/resections.py,sha256=WGGsWgN9LCcJ9YCo-lbgtuA9VqEWnYvoe3OFOJ4nYFk,43619
+pygeodesy/resections.py,sha256=nkfBxYrsEdyq2JqAbPrRWKplJBzFkIXLvEZGB4gPu_A,43550
 pygeodesy/simplify.py,sha256=FT3AzXiFjFfvwJ-X2q8iA9cQOUVdBT1vD_BIskM3gK8,25223
 pygeodesy/solveBase.py,sha256=cBe3QQWEaJlPWUlXoiYo1eFMh4z6GDZbQyt8ZiwMfos,16687
 pygeodesy/sphericalBase.py,sha256=sH2WwmY_RL0STyr1OhtCDKvGhO0Ks7o3V-MHLddpiHQ,31944
@@ -61,7 +61,7 @@ pygeodesy/sphericalNvector.py,sha256=TBs_nid_DTcox7bBegh_Djv4zZRGRY4mH6dd1xIHGms
 pygeodesy/sphericalTrigonometry.py,sha256=Q9qsqzIK43xBken2G8uP9snPXPuBoZ20pMHqbT_IjG0,64216
 pygeodesy/streprs.py,sha256=-29aUO7dYe-YWfTjflljesxGUp7IfEAjk5dJLHWMPVs,23413
 pygeodesy/trf.py,sha256=4gVtF780EEZEqWwMBDe4x1AEkIz3hB8h4l0_0NXWM5U,118903
-pygeodesy/triaxials.py,sha256=Yet0J9Fr4FbTEIL9_p3ziollRR3PzUN2-qt7Z9ZSODU,61537
+pygeodesy/triaxials.py,sha256=fkUcWkXkykZNbTq62LcdnEjhNu9q-e2UaSpRotdzuu8,61546
 pygeodesy/units.py,sha256=ofFTYc7mF9wOpfvpYHN83dGszMM18rTdsxfEXVzaIko,38705
 pygeodesy/unitsBase.py,sha256=xRZlNlO7S47Mwy0FGBg_cefvB4fEtWFM5dAE-Nlw-NE,13111
 pygeodesy/ups.py,sha256=07AG-Rb1D4EvmN3g-xk3ZxKzyaNYvOzztcyYbx7t8jc,23229
@@ -69,7 +69,7 @@ pygeodesy/utily.py,sha256=jP82RmCPCVghxTQv4ZRuIS89afmsD-roz2iafHfyHXA,36128
 pygeodesy/utm.py,sha256=0qF5Kjw8Wf9zxGWOaCF-LSweibat-EkHvvmy7dFO_KE,30981
 pygeodesy/utmups.py,sha256=u2hlWSnmZKwR-ApgKw16R9p75sW_yaqQAmbcBUE8Zfg,13007
 pygeodesy/utmupsBase.py,sha256=KetB8R0uRYU6X8nLG3JkBHUr1OfNlChKa0isgvm8Yq8,18807
-pygeodesy/vector2d.py,sha256=7li4n13cpEm6qNVqlwoE7RH6Cjs9HpbAXBrPX8F4QmI,34963
+pygeodesy/vector2d.py,sha256=IAgmn-9-P6zxh4rond2k0tuf26miVwe-eHqm5f7sWTU,35003
 pygeodesy/vector3d.py,sha256=Q_brX_URZg_Eln01y8G8_NAEuPbbOI_IgO2wBJk1N0U,44970
 pygeodesy/vector3dBase.py,sha256=wJpEcgHlJ5Qu0Us7CKuNDepZQzo6a-RvqZbA2-O4-ko,35583
 pygeodesy/webmercator.py,sha256=1IXiAHBXTm5kMCp5gEI67oBJeRXD7BiPVI7bcRRX1U4,14918
@@ -82,11 +82,11 @@ pygeodesy/auxilats/__main__.py,sha256=TVN8evCKEzBKdANU5-Ucps8bZySqBTQmEGU4hy_U1c
 pygeodesy/auxilats/auxAngle.py,sha256=3tTYEuiX2YmEOmoyqGPYrSN3gLXbcKWKFH1TNt3NqbI,16738
 pygeodesy/auxilats/auxDLat.py,sha256=U3BHF6AkW5--nrwH9whL5wZVIZLH_keRLfUz6c6KBx4,10950
 pygeodesy/auxilats/auxDST.py,sha256=bdNEh3TdABWH7qW6ngZx5SnzDM3B5rtPLvYhvCnFsKc,10497
-pygeodesy/auxilats/auxLat.py,sha256=nefYrJK7C_SlzGX2ju2_0l1lJzCWtJVeoY9YV0saxeU,31271
+pygeodesy/auxilats/auxLat.py,sha256=U6thEwMEoUm4OFb2_JviEtEn2-nqyVjNbII7oyMf0hM,31282
 pygeodesy/auxilats/auxily.py,sha256=ZiCFedYcr32O6TWQ_dcnwMrc9TgogyuhNSbR77WMDpg,8240
-pygeodesy/deprecated/__init__.py,sha256=FJz1aCk2T2ta2S0iA6F92mIA1NlHW2nIzm5kY2oVzNA,2815
+pygeodesy/deprecated/__init__.py,sha256=ZRGqgwSdvtKSYZdw09-LPuGa_i0jMHWKIwrPG3r1hHs,2815
 pygeodesy/deprecated/bases.py,sha256=nxUtxscAg5cb7HEzTSmPJFafLc9Ad2SL5ovwkYaU5z8,1655
-pygeodesy/deprecated/classes.py,sha256=NvZsfSzpdNPTBOA0H02-yYZBdnUBSenGFase2ABu-yg,11202
+pygeodesy/deprecated/classes.py,sha256=PXgss-2kEu3bzaxvSeVZJGdcRqsdi_6_aRqsESyUQKQ,11662
 pygeodesy/deprecated/consterns.py,sha256=W-rQiw443_zCFEUzjCGykGND-i5P_ghEIyJUSGg7-Nc,1908
 pygeodesy/deprecated/datum.py,sha256=KfuMA64LYGMj1voW2uai2RGcugtWdpvlOL6UHJvpTeQ,1875
 pygeodesy/deprecated/functions.py,sha256=-Ckor43KBiCDa7rbwY5kHvrs5xwuixTlw1TZpNSBgrU,13647
@@ -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.12.dist-info/METADATA,sha256=f38SpV-U-BTthJlrzADhDwe_IPavzJf_dj3biK_CRvA,19366
-PyGeodesy-24.4.12.dist-info/WHEEL,sha256=DZajD4pwLWue70CAfc7YaxT1wLUciNBvN_TTcvXpltE,110
-PyGeodesy-24.4.12.dist-info/top_level.txt,sha256=cEQPatCXzKZqrivpULC5V5fuy9_V_bAwaP_gUGid7pQ,10
-PyGeodesy-24.4.12.dist-info/RECORD,,
+PyGeodesy-24.4.24.dist-info/METADATA,sha256=IGu3-2Wu7HfzR8MYGeOU_oWysKdkDfMkndw45lFqWvY,19366
+PyGeodesy-24.4.24.dist-info/WHEEL,sha256=DZajD4pwLWue70CAfc7YaxT1wLUciNBvN_TTcvXpltE,110
+PyGeodesy-24.4.24.dist-info/top_level.txt,sha256=cEQPatCXzKZqrivpULC5V5fuy9_V_bAwaP_gUGid7pQ,10
+PyGeodesy-24.4.24.dist-info/RECORD,,

pygeodesy/__init__.py

@@ -586,7 +586,7 @@ else:
     _init__all__ = False
 
 from pygeodesy.interns import _DOT_, _version2  # PYCHOK import
-__version__ = '24.04.12'
+__version__ = '24.04.24'
 # see setup.py for similar logic
 version     = _DOT_(*_version2(__version__, n=3))
 

pygeodesy/albers.py

@@ -22,7 +22,7 @@ from pygeodesy.constants import EPS0, EPS02, _EPSqrt as _TOL, \
 from pygeodesy.datums import _ellipsoidal_datum, _WGS84
 from pygeodesy.errors import _ValueError, _xkwds
 from pygeodesy.fmath import hypot, hypot1, sqrt3
-from pygeodesy.fsums import Fsum, fsum1f_
+from pygeodesy.fsums import Fsum, _Fsum1f_, fsum1f_
 from pygeodesy.interns import NN, _COMMASPACE_, _datum_, _gamma_, _k0_, \
                              _lat_, _lat1_, _lat2_, _lon_, _name_, _not_, \
                              _negative_, _scale_, _SPACE_, _x_, _y_
@@ -38,7 +38,7 @@ from pygeodesy.utily import atan1, atan1d, degrees360, sincos2, sincos2d, \
 from math import atan2, atanh, degrees, fabs, radians, sqrt
 
 __all__ = _ALL_LAZY.albers
-__version__ = '24.04.07'
+__version__ = '24.04.14'
 
 _k1_    = 'k1'
 _NUMIT  =   8  # XXX 4?
@@ -54,14 +54,6 @@ def _ct2(s, c):
     return c, (s / c)
 
 
-def _Fsum1(a, b, c=_0_0):  # floats=True
-    '''(INTERNAL) C{Fsum} 1-primed.
-    '''
-    S = Fsum()
-    S._facc_(_1_0, a, b, c, _N_1_0, up=False)
-    return S
-
-
 def _Ks(**name_k):
     '''(INTERNAL) Scale C{B{k} >= EPS0}.
     '''
@@ -388,12 +380,13 @@ class _AlbersBase(_NamedBase):
 
         den = hypot(nx, y1) + nrho0  # 0 implies origin with polar aspect
         if den:
-            drho = _Fsum1(x * nx, y_ * nrho0 * _N_2_0, y_ * ny).fover(den / k0)
+            drho = _Fsum1f_(x * nx, y_ * nrho0 * _N_2_0, y_ * ny).fover(den / k0)
             # dsxia = scxi0 * dsxi
-            t  +=  drho  * n0
+            t  +=  drho  * n0  # k0 below
             d_  = (nrho0 + t) * drho * self._scxi0_  # / (qZ * E.a2)
-            d_2 = (txi * _2_0 - d_) * d_ + _1_0
-            txi = (txi - d_) / (sqrt(d_2) if d_2 > EPS02 else EPS0)
+            t_  =  txi - d_
+            d_  = (txi + t_) * d_ + _1_0
+            txi =  t_ / (sqrt(d_) if d_ > EPS02 else EPS0)
 
         ta  = self._tanf(txi)
         lat = atan1d(s * ta)
@@ -487,11 +480,11 @@ class _AlbersBase(_NamedBase):
                                        (ca2, sa2, ta2, scb22))
 
         dsxi = ((esa12 / esa1_2) + _DatanheE(sa2, sa1, E)) * dsn_2 / self._qx
-        C = _Fsum1(sxi * dtb12 / dsxi, scb22, scb12).fover(scb22 * scb12 * _2_0)
+        C = _Fsum1f_(sxi * dtb12 / dsxi, scb22, scb12).fover(scb22 * scb12 * _2_0)
 
-        sa12  =  fsum1f_(sa1, sa2, sa12)
-        axi  *= (sa12 * e2 + _1) / (sa12 + _1)
-        bxi  *= _Fsum1(sa1, sa2, esa12).fover(esa1_2) * e2 + _D2atanheE(sa1, sa2, E) * E.e21
+        S     = _Fsum1f_(sa1, sa2, sa12)
+        axi  *= (S * e2 + _1).fover(S + _1, raiser=False)
+        bxi  *= _Fsum1f_(sa1, sa2, esa12).fover(esa1_2) * e2 + _D2atanheE(sa1, sa2, E) * E.e21
         s1_qZ = (axi * self._qZ - bxi) * dsn_2 / dtb12
         ta0   =  self._ta0(s1_qZ, (ta1 + ta2) * _0_5, E)
         return ta0, C
@@ -571,7 +564,7 @@ class _AlbersBase(_NamedBase):
         es2m1a =  es2m1 * E.e21  # e2m
         s =  sqrt((ca2 / (es1p1 * es2m1a) + _atanheE(ca2 / es1m1, E))
                        * (es1m1 / es2m1a  + _atanheE(es1p1, E)))
-        t = _Fsum1(sa / es2m1, _atanheE(sa, E)).fover(s)
+        t = _Fsum1f_(sa / es2m1, _atanheE(sa, E)).fover(s)
         return neg(t) if ta < 0 else t
 
 

pygeodesy/auxilats/auxLat.py

@@ -27,7 +27,7 @@ from pygeodesy.datums import _ellipsoidal_datum, _WGS84,  Ellipsoid
 from pygeodesy.elliptic import Elliptic as _Ef
 from pygeodesy.errors import AuxError, _xkwds, _xkwds_get, _Xorder
 # from pygeodesy.fmath import cbrt  # from .karney
-from pygeodesy.fsums import Fsum, _sum
+from pygeodesy.fsums import Fsum, _Fsumf_, _sum
 from pygeodesy.karney import _2cos2x, _polynomial,  _ALL_DOCS, cbrt, _MODS
 from pygeodesy.interns import NN, _DOT_, _UNDER_  # _earth_
 # from pygeodesy.lazily import _ALL_DOCS, _ALL_MODS as _MODS  # from .karney
@@ -44,7 +44,7 @@ except ImportError:  # Python 3.11-
         return pow(_2_0, x)
 
 __all__ = ()
-__version__ = '23.12.10'
+__version__ = '24.04.14'
 
 _TRIPS = 1024  # XXX 2 or 3?
 
@@ -286,7 +286,7 @@ class AuxLat(AuxAngle):
                     atphi  = _asinh_2(tphi)  # atanh(sphi)
                     t      = _asinh_2(em1 * (tphi * scphib)) / em1
                     try:
-                        Dg =  Fsum(atphi, atphib, t, e * t)
+                        Dg = _Fsumf_(atphi, atphib, t, e * t)
                     except ValueError:  # Fsum(NAN) exception
                         Dg = _sum((atphi, atphib, t, e * t))
                     e *= atphib

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

pygeodesy/deprecated/classes.py

@@ -22,7 +22,7 @@ from pygeodesy.trf import TRFXform7Tuple as _TRFXform7Tuple
 from pygeodesy.units import Bearing, Int, Lat, Lon, Meter
 
 __all__ = _ALL_DEPRECATED.deprecated_classes
-__version__ = '24.02.02'
+__version__ = '24.04.17'
 
 
 class _Deprecated_NamedTuple(_NamedTuple):
@@ -98,6 +98,19 @@ def EcefCartesian(*args, **kwds):
     return EcefCartesian_(*args, **kwds)
 
 
+def Fn_rt(root, *xs, **name_RESIDUAL_raiser):  # PYCHOK no cover
+    '''DEPRECATED on 2024.04.17, use class L{Froot}.
+    '''
+    Froot = _MODS.fmath.Froot
+
+    class Fn_rt(Froot):  # PYCHOK no cover
+        '''DEPRECATED on 2024.04.17, use class L{Froot}.
+        '''
+        def __init__(self, root, *xs, **name_RESIDUAL_raiser):
+            deprecated_class(self.__class__)
+            Froot.__init__(self, root, *xs, **name_RESIDUAL_raiser)
+
+
 def HeightIDW(knots, **kwds):  # PYCHOK no cover
     '''DEPRECATED, use class L{HeightIDWeuclidean}.'''
     HeightIDWeuclidean = _MODS.heights.HeightIDWeuclidean

pygeodesy/ellipsoids.py

@@ -70,8 +70,8 @@ from pygeodesy.constants import EPS, EPS0, EPS02, EPS1, INF, NINF, PI4, PI_2, PI
                                _0_0s, _0_0, _0_5, _1_0, _1_EPS, _2_0, _4_0, _90_0, \
                                _0_25, _3_0  # PYCHOK used!
 from pygeodesy.errors import _AssertionError, IntersectionError, _ValueError, _xattr, _xkwds_not
-from pygeodesy.fmath import cbrt, cbrt2, fdot, Fhorner, fpowers, Fsum, hypot, hypot_, \
-                            hypot1, hypot2, sqrt3
+from pygeodesy.fmath import cbrt, cbrt2, fdot, Fhorner, fpowers, hypot, hypot_, \
+                            hypot1, hypot2, sqrt3,  Fsum
 # from pygeodesy.fsums import Fsum  # from .fmath
 from pygeodesy.interns import NN, _a_, _Airy1830_, _AiryModified_, _b_, _Bessel1841_, _beta_, \
                              _Clarke1866_, _Clarke1880IGN_, _DOT_, _f_, _GRS80_, _height_, \
@@ -93,7 +93,7 @@ from pygeodesy.utily import atan1, atan1d, atan2b, degrees90, m2radians, radians
 from math import asinh, atan, atanh, cos, degrees, exp, fabs, radians, sin, sinh, sqrt, tan
 
 __all__ = _ALL_LAZY.ellipsoids
-__version__ = '24.04.07'
+__version__ = '24.04.14'
 
 _f_0_0    = Float(f =_0_0)  # zero flattening
 _f__0_0   = Float(f_=_0_0)  # zero inverse flattening

pygeodesy/elliptic.py

@@ -97,7 +97,7 @@ from math import asinh, atan, atan2, ceil, cosh, fabs, floor, \
                  radians, sin, sqrt, tanh
 
 __all__ = _ALL_LAZY.elliptic
-__version__ = '24.04.09'
+__version__ = '24.04.14'
 
 _TolRD  =  zqrt(EPS * 0.002)
 _TolRF  =  zqrt(EPS * 0.030)
@@ -1023,13 +1023,12 @@ def _Horner(S, e1, E2, E3, E4, E5, *over):
     #    + 3*E5/26 - E2**3/16 + 3*E3**2/40 + 3*E2*E4/20
     #    + 45*E2**2*E3/272 - 9*(E3*E4+E2*E5)/68)
     # converted to Horner-like form ...
-    F  = Fsum
     e  = e1 * 4084080
     S *= e
-    S += F(E2 * -540540,                               471240).fmul(E5)
-    S += F(E2 *  612612,                E3 * -540540, -556920).fmul(E4)
-    S += F(E2 * -706860, E22 *  675675, E3 *  306306,  680680).fmul(E3)
-    S += F(E2 *  417690, E22 * -255255,               -875160).fmul(E2)
+    S += Fsum(E2 * -540540,                               471240).fmul(E5)
+    S += Fsum(E2 *  612612,                E3 * -540540, -556920).fmul(E4)
+    S += Fsum(E2 * -706860, E22 *  675675, E3 *  306306,  680680).fmul(E3)
+    S += Fsum(E2 *  417690, E22 * -255255,               -875160).fmul(E2)
     S += 4084080
     if over:
         e *= over[0]
@@ -1122,14 +1121,14 @@ def _rG2(inst, x, y, PI_=PI_4):  # 2-args
     '''(INTERNAL) Carlson, eqs 2.36 - 2.39.
     '''
     m = -1  # neg!
-    S = _Dsum()
-    # assert not S
+    S = None
     for a, b in _ab2(inst, x, y):  # PYCHOK yield
-        if S:
+        if S is None:  # initial
+            S  = _Dsum()
+            S += (a + b)**2 * _0_5
+        else:
             S += (a - b)**2 *  m
             m *=  2
-        else:  # initial
-            S += (a + b)**2 * _0_5
     return S(PI_).fover(a + b)
 
 

pygeodesy/errors.py

@@ -11,9 +11,9 @@ chaining}, use command line option C{python -X dev} I{OR} set env variable
 C{PYTHONDEVMODE=1} or to any non-empty string I{OR} set env variable
 C{PYGEODESY_EXCEPTION_CHAINING=std} or to any non-empty string.
 '''
-# from pygeodesy.basics import isint, isodd, itemsorted, _xinstanceof, _zip  # _MODS
+# from pygeodesy.basics import isint, isodd, issubclassof, itemsorted, _xinstanceof, _zip  # _MODS
 # from pygeodesy.ellipsoidalBase import CartesianEllipsoidalBase, LatLonEllipsoidalBase  # _MODS
-# from pygeodesy import errors  # _MODS.getattr
+# from pygeodesy import errors  # _MODS, _MODS.getattr
 from pygeodesy.interns import MISSING, NN, _a_, _an_, _and_, _clip_, _COLON_, \
                              _COLONSPACE_, _COMMASPACE_, _datum_, _ellipsoidal_, \
                              _incompatible_, _invalid_, _len_, _not_, _or_, _SPACE_, \
@@ -25,7 +25,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.03.24'
+__version__ = '24.04.18'
 
 _box_        = 'box'
 _limiterrors =  True  # in .formy
@@ -425,7 +425,7 @@ def _error_init(Error, inst, args, fmt_name_value='%s (%r)', txt=NN,
 def _error_under(inst):
     '''(INTERNAL) Remove leading underscore from instance' class name.
     '''
-    n = inst.__class__.__name__
+    n = inst.__class__.__name__  # _tailof?
     if n.startswith(_UNDER_):
         inst.__class__.__name__ = n.lstrip(_UNDER_)
     return inst
@@ -479,8 +479,15 @@ def isError(exc):
                 C{False} if B{C{exc}} is a standard Python error
                 of C{None} if neither.
     '''
-    return True  if isinstance(exc, _XErrors)   else (
-           False if isinstance(exc,  Exception) else None)
+    def _X(exc):
+        X = type(exc)
+        m = X.__module__
+        return _MODS.basics.issubclassof(X, *_XErrors) or \
+               ((m is __name__ or m == __name__) and
+               _tailof(X.__name__).startswith(_UNDER_))
+
+    return True   if isinstance(exc, _XErrors)   else (
+          _X(exc) if isinstance(exc,  Exception) else None)
 
 
 def _IsnotError(*nouns, **name_value_Error_cause):  # name=value [, Error=TypeError, cause=None]
@@ -665,31 +672,23 @@ def _xError2(exc):  # in .constants, .fsums, .lazily, .vector2d
 
        @arg exc: The exception instance (usually, C{Exception}).
     '''
-    m =  __name__  # 'pygeodesy.errors'
-    X =  type(exc)
-    n =  NN(_UNDER_, _tailof(X.__name__))
-    E = _MODS.getattr(m, n, X)  # == _X2Error.get(X, X)
-    if E is X and not isError(exc):
-        E = _NotImplementedError
-        t =  repr(exc)
-    else:
-        t =  str(exc)
-    return E, t
+    x = isError(exc)
+    if x:
+        E =  type(exc)
+    elif x is None:
+        E = _AssertionError
+    else:  # get _Error from Error
+        n =  NN(_UNDER_, _tailof(type(exc).__name__))
+        E = _MODS.getattr(__name__, n, _NotImplementedError)
+        x =  E is not _NotImplementedError
+    return E, (str(exc) if x else repr(exc))
 
 
-_XErrors = _TypeError, _ValueError
+_XErrors = (_AssertionError, _AttributeError,  # some isError's
+            _TypeError, _ValueError, _ZeroDivisionError)
 # map certain C{Exception} classes to the C{_Error}
-# _X2Error = {AssertionError:      _AssertionError,
-#             AttributeError:      _AttributeError,
-#             ImportError:         _ImportError,
-#             IndexError:          _IndexError,
-#             KeyError:            _KeyError,
-#             NameError:           _NameError,
-#             NotImplementedError: _NotImplementedError,
-#             OverflowError:       _OverflowError,
-#             TypeError:           _TypeError,
-#             ValueError:          _ValueError,
-#             ZeroDivisionError:   _ZeroDivisionError}
+# _X2Error = {AssertionError:    _AssertionError, ...
+#             ZeroDivisionError: _ZeroDivisionError}
 
 try:
     _ = {}.__or__  # {} | {}  # Python 3.9+