L iuddlZddlZddlmZddlZddlmZm Z m Z m Z m Z m Z mZmZmZmZmZmZmZmZmZmZmZmZmZmZmZmZmZmZm Z m!Z!m"Z"m#Z#m$Z$m%Z%m&Z&m'Z'm(Z(m)Z)m*Z*m+Z+m,Z,m-Z-m.Z.m/Z/m0Z0m1Z1m2Z2m3Z3m4Z4m5Z5m6Z6m7Z7m8Z8m9Z9m:Z:m;Z;mZ>m?Z?m@Z@mAZAmBZBmCZCmDZDmEZEmFZFmGZGmHZHmIZImJZJmKZKmLZLmMZMmNZNmOZOmPZPmQZQmRZRmSZSmTZTmUZUddlVmWZWddlXmYZYejjdZ\e\jdZ^ejje^5Z`eje`Zbddde\jdZcejjec5Z`eje`Zdddde\jd Zeejjee5Z`eje`Zfdddd Zgd Zhd Zid ZjdZkdZldZmdZndZodZpdZqdZrdZsdZtdZudZvdZwdZxdZydZzdZ{dZ|d Z}d!Z~d"Zd#Zd$Zd%Zd&Zd'Zd(Zd)Zd*Zd+Zd,Zd-Zd.Zd/Zd0Zd1Zd2Zd3Zd4Zd5Zd6Zgegeod7d8d9d:;egepdd:;egepd?d=d>d@;egesddAdBCegesd?d=d>dDdBCegeqddAdBCegeqd?d=d>dDdBCege{ddE;ege{d?d=d>dF;ege|ddG;ege|d?d=d>dF;egerddJ;ege)dId=d>dJ;ege)dKd=d>dL;ege*dMd8dNdO;ege*dPd8dNdL;ege*dQd8dNdG;ege*dRd8dNdS;ege+dTd8d9dU;ege?dVdWd9dX;egedVdWdYdX;ege@dVdZdNdX;egedVdZd[dX;eged\d=d>dJ;eged]d=d>dL;ege>d^dWd9dX;eged^dWdYdX;egeFd_dWd9d`;eged_dWdYd`;egeGdadbd>dc;egedadbd9d`;egedddedege(dfddeege(dfdgdeege(dhddedF;ege(dhdgdedJ;ege(diddedj;ege(didgdedj;ege(dkddedj;ege(dkdgdedl;egejdmddeegekdnd=d>dl;egeldoddeegemdpd=d>dl;ege7dqddeege7dqdgdedJ;ege8dqdd>dO;ege7drddeege7drdgdeege8drdd>ds;ege7dtddeege7dtdgdedJ;ege8dtdd>ege9duddeege:dvddeege:dwddedxyege;dzded>dJ;ege;dzd{d>dX;egeege=dd=d>dl;ege"dddeege"dddeege"dddeege"dddeege"ddded;ege"dddeds;ege"ddgded;ege"ddgded;ege"ddgded;ege"ddgded;ege"ddgded;ege"ddgdedF;ege#ddd>d;ege#ddd>d;ege#ddd>d;ege#ddd>d;ege#ddd>d;ege#ddd>ege$dd=dNd;ege$dd=dNdF;ege$dd=dNdF;ege$dd=dNd;egedd=dNdJ;egedd=dNdF;egedd=dNdF;egedd=dNd;ege%dd=d9dJ;ege%dd=d9dF;ege%dd=d9ds;ege%dd=d9d:;egedd=d9dJ;egedd=d9dF;egedd=d9ds;egedd=d9d:;egeddbd>dX;egeddbd9dX;egedd=d>dF;egedd=d>egedd=d>egedd=d9dF;egedd=d9egedd=d9ege}dd=d>dG;ege~dd=d>d;egedd=d>degedd8d9dF;egedd8d9degedddeegeddd>egedd=d>d;egeddd>ddegedd=d>d;egeddd>d;egedd=d>d;egeddd>ddegedd=d>d;egeddd>d;egedd=d>d;egeddd>d;egedd=d>d;egeddd>d;egedd=d>d;egeddd>d;egedd=d>dJ;egeddd>dJ;egedd=d>d:;egeddd>d;egedd=d>d;ege dd=d>d;ege ddd>d;ege dd=d>d;ege ddd>d;ege!dd=d>d;ege!ddd>d;ege!dd=d>d;ege!ddd>d;egedd=d>d;egedd=d>d;egedd=d>d;egeddd>d;egedd=d>d;egeddd>d;egedd=d>d`;egeddd>d`;egendddedyegendddedyegendddedyegendddedyege&dd=d>d:;ege&dd=d>dl;ege&dd=d>d:;ege'dd=d9d;ege'dd=d9dl;ege'dd=d9dl;egedd=d>dDd«egedd=d>dì;egedd=d>egedd=d>dĬegedd=d9dO;egedd=d9egeAddWd9dƬ;egeAddWd9ddɬ«egeddddd̈́d΄dτdЄfѫegedd=d>dJ;egedd=d>dԬ;egeHdd=d>dA;egeHdd=d>d׬;ege2dd=d>d٬;ege3dd8d9d٬;ege3dd8d9d٬;ege3dd8d9d٬;ege3dd8d9d٬;ege3dd8d9d߬;ege3dd8d9d;ege4dd8d9d٬;ege4dd8d9d٬;ege4dd8d9d٬;ege4dd8d9d٬;ege4dd8d9d٬;ege5dd8d9d;ege5dd8d9d;ege5dd8d9d;ege5dd8d9d٬;ege5dd8d9d٬;ege6dddYdddfZej*j,ej*j/deedZeheIdd=d>dJdJeheJdd=d9dJdJeheudd8d9ddJehevdd8dYddJehewdd8d9ddJehexdd8dYddJeheydd8dNddJehezdd8d[ddJgZej*j/deedZeie1dd=d>eie0dddeeie.dd=d>eiedd=ddl;eiedd=ddl;eie$dd=d>d;eie%dd=d>d:;eieRddd[ddJeieQdddNddJeieUdd8d9d:;g Zej*j/deedZej<fdZy#1swY PxYw#1swY xYw#1swY xYw(N)suppress_warnings)Nlpnlpmnlpmvlqnlqmnsph_harm eval_legendre eval_hermite eval_laguerreeval_genlaguerrebinomcbrtexpm1log1pzetajnjvjvpynyvyvpivivpknkvkvpgammagammalngammainc gammaincc gammaincinv gammainccinvdigammabetabetainc betaincinvpochellipe ellipeincellipkellipkm1 ellipkincelliprcelliprdelliprfelliprgelliprjerferfcerfinverfcinvexp1expiexpnbdtrikbtdtriabtdtribchndtrgdtrgdtrcgdtrixgdtribnbdtrikpdtrikowens_t mathieu_a mathieu_b mathieu_cem mathieu_semmathieu_modcem1mathieu_modsem1mathieu_modcem2mathieu_modsem2 ellip_harm ellip_harm_2 spherical_jn spherical_yn wright_bessel)IntegrationWarning)FuncDatazscipy.special.tests.dataz boost.npzzgsl.npzz local.npzcT|jd|t|t|g|i|SNdataname) setdefaultrSDATASETS_BOOSTfuncrVakws c/mnt/ssd/data/python-lab/Trading/venv/lib/python3.12/site-packages/scipy/special/tests/test_data.pydatar^,,MM*h' D.2 =Q =" ==cT|jd|t|t|g|i|SrU)rWrS DATASETS_GSLrYs r]data_gslrc1s,MM*h' D,x0 ;1 ; ;;r`cT|jd|t|t|g|i|SrU)rWrSDATASETS_LOCALrYs r] data_localrf6r_r`ct||zSN)r+ks r]ellipk_rkD !A#;r`c t|||zSrh)r-frjs r] ellipkinc_rpH Q! r`ct||zSrh)r)ris r]ellipe_rsLrlr`c t|||zSrh)r*rns r] ellipeinc_ruPrqr`ct|dSN?)r)xs r]zeta_rzTs 2;r`cBt||jt|Srh)rastypeint)numurys r]assoc_legendre_p_boost_rXs BIIcNA &&r`ctd||SNr)r)r~rys r]legendre_p_via_assoc_r\s 2q>r`ct5}|jtt|j d|ddcdddS#1swYyxYw)Ncategorylr)rfilterDeprecationWarningrr|nrysups r]lpn_r_sJ  , . /188C=!$Q'+,,,s 7A  AcDt|jd|ddS)Nrrr)rr|rrys r]lqn_rds! qxx}a  #B ''r`ct5}|jttd||ddcdddS#1swYyxYw)Nrrrr)rrrrrs r]legendre_p_via_lpmnrgsC  & . /Aq!}Q%&&&s )>Ac(td||ddS)Nrr)rrs r]legendre_q_via_lqmnrls 1a= D !!r`cJt|||dztjz dSNr)rGnppimqrys r]mathieu_ce_radro" q!QsU255[ )! ,,r`cJt|||dztjz dSr)rHrrrs r]mathieu_se_radrsrr`cpt|||dtjtjdz zSNr)rIrsqrtrrs r]mathieu_mc1_scaledrws. 1a #A &q)9 99r`cpt|||dtjtjdz zSr)rJrrrrs r]mathieu_ms1_scaledr}, 1a #A &q)9 99r`cpt|||dtjtjdz zSr)rKrrrrs r]mathieu_mc2_scaledrrr`cpt|||dtjtjdz zSr)rLrrrrs r]mathieu_ms2_scaledrrr`c8t|jd|SNrr r|rs r]eval_legendre_ldr # **r`c8t|jd|SNdrrs r]eval_legendre_ddrrr`c8t|jd|Sr)r r|rs r]eval_hermite_ldr  q ))r`c8t|jd|Srr r|rs r]eval_laguerre_ldrrr`c8t|jd|Srrrs r]eval_laguerre_ddrrr`c:t|jd||Srr r|rr[rys r]eval_genlaguerre_lddr AHHSM1a 00r`c:t|jd||Srrrs r]eval_genlaguerre_dddrrr`c"td|z ||SN)r:yrps r] bdtrik_comprs !A#q! r`c"td|z ||Sr)r;)rbrys r] btdtria_compr 1Q31 r`c"t|d|z |Sr)r<)r[rrys r] btdtrib_comprs 1ac1 r`ctd||Srw)r>rrys r]gdtr_rs Q?r`ctd||Srw)r?rs r]gdtrc_rs a r`ctd||Srwr@rrs r]gdtrix_r #q! r`c"td|d|z SNrxrrrs r] gdtrix_comprs #q!A# r`ctd||SrwrArs r]gdtrib_rrr`c"tdd|z |Srrrs r] gdtrib_comprs #qsA r`c"td|z ||Sr)rBrs r] nbdtrik_comprrr`c td|z |Sr)rC)rrs r] pdtrik_comprs !A#q>r`c dt||z Srwr(zrs r]poch_rs a r`c"dt|| z Srwrrs r] poch_minusrs a! r`c8t|jd|Sr)rOr|rs r] spherical_jn_rrr`c8t|jd|Sr)rPr|rs r] spherical_yn_rrr`ct5}|jtt||||}dddj|j fS#1swY!xYw)Nr)rrrr realimag)rrthetaphirrs r] sph_harm_rsS  ' . / Q5# &' FFAFF ''s %AAcTt|d|zz}|j|jfSN?)rrrryrrs r]cexpm1r% a"Q$hA 66166>r`cTt|d|zz}|j|jfSr)rrrrs r]clog1prrr`z%assoc_legendre_p_ipp-assoc_legendre_p)rrrgdy=)rtolzlegendre_p_ipp-legendre_p)rrrz%legendre_p_large_ipp-legendre_p_largeg [n;=gvIh%,=F)r vectorizedg7QU=gt0=gvIh%L=g+=g-=zbeta_exp_data_ipp-beta_exp_datagvIh%<=zbeta_med_data_ipp-beta_med_datag-a=z%ibeta_small_data_ipp-ibeta_small_datag [nz#ibeta_inva_data_ipp-ibeta_inva_data)rrrg:0y5>)rrrzbinomial_data_ipp-binomial_dataz+binomial_large_data_ipp-binomial_large_dataz,binomial_quantile_ipp-binomial_quantile_dataz>negative_binomial_quantile_ipp-negative_binomial_quantile_datag& .1>z*poisson_quantile_ipp-poisson_quantile_data)rrgA:)>zcbrt_data_ipp-cbrt_datarzdigamma_data_ipp-digamma_datayz%digamma_neg_data_ipp-digamma_neg_dataz'digamma_root_data_ipp-digamma_root_datagV瞯r J QUr`) param_filterzexpint_1_data_ipp-expint_1_datarzexpinti_data_ipp-expinti_dataz+expinti_data_double_ipp-expinti_data_doublez'expinti_data_long_ipp-expinti_data_longz'expint_small_data_ipp-expint_small_datazexpint_data_ipp-expint_dataztest_gamma_data_ipp-near_0ztest_gamma_data_ipp-near_1ztest_gamma_data_ipp-near_2ztest_gamma_data_ipp-near_m10ztest_gamma_data_ipp-near_m55gYLɞ=ztest_gamma_data_ipp-factorialsg& .!>g|=g|=z'igamma_small_data_ipp-igamma_small_datag+z3bessel_i_prime_int_data_ipp-bessel_i_prime_int_datagt@=gYnz'bessel_j_int_data_ipp-bessel_j_int_dataz+bessel_j_large_data_ipp-bessel_j_large_datagp ~=zbessel_j_data_ipp-bessel_j_dataz3bessel_j_prime_int_data_ipp-bessel_j_prime_int_dataz7bessel_j_prime_large_data_ipp-bessel_j_prime_large_datagdy=z'bessel_k_int_data_ipp-bessel_k_int_datazbessel_k_data_ipp-bessel_k_dataz3bessel_k_prime_int_data_ipp-bessel_k_prime_int_datagt =z+bessel_k_prime_data_ipp-bessel_k_prime_datag(3=z#bessel_y01_data_ipp-bessel_y01_dataz!bessel_yn_data_ipp-bessel_yn_dataz!bessel_yv_data_ipp-bessel_yv_datag|=z-bessel_yv_prime_data_ipp-bessel_yv_prime_datazzeta_data_ipp-zeta_datac |dkDSrr r s r]r r rr`zzeta_neg_data_ipp-zeta_neg_datac |dkDSrr r s r]r r rr`z!zeta_1_up_data_ipp-zeta_1_up_datac |dkDSrr r s r]r r rr`z'zeta_1_below_data_ipp-zeta_1_below_datac |dkDSrr r s r]r r rr`z-gamma_inv_small_data_ipp-gamma_inv_small_dataz!gamma_inv_data_ipp-gamma_inv_dataz)gamma_inv_big_data_ipp-gamma_inv_big_datazgdtrix unflow some points)rrg [nznccs_big_ipp-nccs_biggMb@?zchndtr inaccurate some pointsz)spherical_harmonic_ipp-spherical_harmonic)rrrr)rrcBtj|jdSN?ronesshapers r]r r RWWQWWc%:r`cBtj|jdSrrrs r]r r rr`c\tj|dtjzk|dk\S)Nrrr logical_andrrs r]r r s!R^^A"%%Ka%Hr`cVtj|tjk|dk\Srr rs r]r r  sR^^AIqAv%Fr`)rrz#sph_bessel_data_ipp-sph_bessel_dataz%sph_neumann_data_ipp-sph_neumann_datagV瞯=zowens_t_ipp-owens_tz)owens_t_large_data_ipp-owens_t_large_datag-=z!ellint_rc_data_ipp-ellint_rc_datagV瞯 mathieu_mc_msct|yrhr'r)s r]test_gslr/ $r`ellipkinc_neg_mr,ellipeinc_neg_mlog1p_expm1_complex)rrr r!ellip)rrrrrrQct|yrhr'r)s r] test_localr6r0r`ct5}|jtdtjd5|j |ddddddy#1swYxYw#1swYyxYw)z Boost testz,The occurrence of roundoff error is detectedignore)all)dtypeN)rrrRrerrstatecheck)r$r:rs r]r(r(sc  $ %'UV [[X & $ JJUJ # $$$ $ $$$s"-A(A A(A% !A((A1)importlib.resources importlibnumpyr numpy.testingrpytest scipy.specialrrrrrr r r r r rrrrrrrrrrrrrrrrrrr r!r"r#r$r%r&r'r(r)r*r+r,r-r.r/r0r1r2r3r4r5r6r7r8r9r:r;r<r=r>r?r@rArBrCrDrErFrGrHrIrJrKrLrMrNrOrPrQscipy.integraterRscipy.special._testutilsrS resourcesfiles_datadirjoinpath _boost_npzas_fileroloadrX_gsl_npzrb _local_npzrer^rcrfrkrprsrurzrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr BOOST_TESTSmark thread_unsafe parametrizereprr* GSL_TESTSr/ LOCAL_TESTSr6float64r(r r`r]rVs+ /-    $ $%? @   { +   , RWWQZN    Y '  *a2771:L  { +   , RWWQZN > < >', (& "--: :::++*++11**  | $&M ae %| "$? AE # | "$K AG % | "= AEe 5| "I AEe 5| T. AEe 5| T: AEe 5| : AE #!|$ F AE #%|( : AE #)|, F AE #-|2 T. AEe 53|6 T: AEe 57|: "= AEe 5;|> "I AEe 5?|D T4 AE #E|H T4 AE #I|L T4 AE #M|R W= ae %S|V W1 ae %W|Z W9 ae %[|^ W= ae %_|d Z< ad $e|j W; ad $k|n \@ ad $o|t W; ad $u|x \@ ad $y|~ U5 AE #|B UA AE #C|H VC ad $I|L [H ad $M|R WV ad $S|V \ M ad $W|^ VA AD "_|b [F AD "c|h T,a3i|l W5q!|@ Z:E15QA|B W7A>C|D Z_|` W?A* -a|f T4aGg|h T4b!$Gi|j T2AquEk|l T@!QUSm|n TquUE|F V:E15QG|H V:E15QI|J V:E15QK|N W=uadSO|P [A5!RVWQ|T UI AE #U|X UG A Y|\ UI A ]|` ZN AE #a|d ZL A e|h ZN A i|n _3 AE #o|t 8 AE #u|x 8 A$E Gy|| !#< ae %}|@ !#< a&G IA|F U;QBG|H U;QBI|L R2 AE #M|P R2 QU 1Q|T R: AD "U|X R: QU $Y|^ SG AG %_|b SG QW6 3c|h R:E15Qi|j R:FAERk|l R>quUm|n R>PUVo|r R:E15Qs|t R:FAERu|v R2E15Iw|x R2FAEJy|| SG AE #}|@ SG QU $A|D SK AE #E|H SK QU $I|N R:E15QO|R R:E15QS|T R:FAERU|V R2E15IW|X R2FAEJY|\ SG AE #]|` SG QU $a|d S?PUVe|f S?QVWg|j R6quMk|l R4eQUKm|p R4eQUKq|r R4faeLs|t R4eQUKu|v R4faeLw|z SA QT #{|~ SA ad $|D U-q!* -E|H U5q!* -I|L U7A* -M|P U=q!* -Q|V [I AE #W|Z [= AE #[|^ [E AE #_|d \J AE #e|h \> AE #i|l \F AE #m|r WE AE0K Ms|v W9 AE #w|z WA A {|~ [I A$< >|B [= AE #C |F [E A G |L V_ ad $M |P V, ad1P RQ |V YC E::HFH IW |d ]A AE #e |h ]C AE #i |n W+ QU $o |r WA QU $s |Z W961 [ |^ W99a _ |b W7A c |f W7A g |j W7A k || W7A } |@ W7A A |D W99a E |H W7A I |L W7A M |P W7A Q |T W7A U |d W/A e |h W7A i |l W7A m |p W7A q |t W7A u |x W9<e/1 3y | ~ $78 L&!%eLL&!%eL )QTPUV)QTPUV#_ADu 6#_ADu 6 #_ADu 6#_ADu 6  (56 y+VQ7xQ*y+VQ7v,eUGv,eUGxVQU;y+vqu=|W_ae%Pz7?AEN}oy!%H  $78!jj$U    s$u' u4v'u14u>v