K i mddlddlmZmZmZdZdZdZdZdZ dZ d Z d Z d Z d Zd ZdZdZdZdZdZdZdZdZdZdZdZdZdZdZdZdZdZ dZ!d Z"y!)")*)round_up from_float mpf_zeta_intc8tddtdk(sJy)Nr )rra/mnt/ssd/data/python-lab/Trading/venv/lib/python3.12/site-packages/mpmath/tests/test_gammazeta.pytest_zeta_int_bugr s 2 *T"2 22 2r ctddk(sJtddk(sJtddk(sJtddk(sJtd d k(sJtd dk(sJtd d k(sJtdd k(sJtddk(sJtddk(sJtddk(sJtd\}}|dzdk(sJ|dk(sJtd\}}|dzdk(sJ|dk(sJdt_tddk(sJtddk(sJtdj dsJtddk(sJtd j dsJtd dk(sJtd j d sJt tdd!k(sJt td"d#k(sJt td$d%k(sJt td&d'k(sJd(t_t tdd)k(sJt td"d*k(sJt td$d+k(sJt td&d,k(sJdt_tdj t d d-z sJd.t_tdj t d d-z sJdt_y)/Nr)rr)r)r)rr)rr)r*r)rB )iMi )i[id( is i lkL3iqr gUUUUUU?ggaa?z0.0757575757575758z7.62772793964344e+267順z-5.82229431461335e+376755iz1.19570355039953e+6767525842z50.075757575757575757575757575757575757575757575757576z87.6277279396434392486994969020496121553385863373331e+267z<-5.8222943146133508236497045360612887555320691004308e+376755z>1.1957035503995297272263047884604346914602088317782e+676752584riP)bernfracmpdps bernoulliaestrmpf)pqs r test_bernoullir,s A;%   A;&  A;%   A;%   A;' !! ! A;%   A;&  A;' !! ! B<6 !! ! B<; && & B<; && & C=DAq v: "" " ;; D>DAq v: ## # >> BF Q<1   Q<4   Q  Q   Aa=  A   Ae    "" " Ae    "" " Aa=  A   Aa=  A   Ab>  V $$ $ Aa=  F ## # Ae    "" " Ae    && & Aa=  E "" " Aa=  E "" " Ab>  V $$ $ Aa=  E "" " Ae    (( ( Ae    (( ( Aa=  E "" " Aa=  F ## # Ab>  R  Aa=  A   Ae    "" " Ae    (( ( Aa=  A   Aa=  A   Ab>  X && & Aa=  G $$ $ Ae    (( ( Ae    )) ) Aa=  G $$ $ Aa=  H %% % Ab>  R  Aa=  A   Ae    "" " Ae    ** * Aa=  A   Aa=  A   Br?  h '' ' Bq>  V $$ $ Bu   // / Bu  !1 22 2 Bq>  V $$ $ Bq>  X && & Br?  c "" " Bq>  Q   Bu   ## # Bu  !2 33 3 Bq>  Q   Bq>  R  Qr?  a  Qq>  Q   Qu   ## # Qu   ## # Qq>  Q   Qq>  Q   Qr?  f %% % Qq>  V $$ $ Qu   ## # Qu   '' ' Qq>  U ## # Qq>  U ## # Qr?  a  Qq>  Q   Qu   (( ( Qu   )) ) Qq>  Q   Qq>  Q   Qr?  f %% % Qq>  U ## # Qu   ## # Qu   ** * Qq>  V $$ $ Qq>  U ## # Qr?  a  Qq>  Q   Qu   (( ( Qu   ** * Qq>  Q   Qq>  Q   Qr?  f %% % Qq>  V $$ $ Qu   ## # Qu   ,, , Qq>  U ## # Qq>  U ## # R   q !! ! R?  a  R  ! !- 00 0 R  ! !"5 66 6 R?  a  R?  a  R   x (( ( R?  g && & R  ! !! $$ $ R  ! !"6 77 7 R?  h '' ' R?  g && & E% # #$B CC C C  !: ;; ; S  ! !"< == =r c(dt_tdjdsJtdjdsJtdjdsJtdjd sJtd jd sJtd jd sJtdjdsJtdjdsJd}t ddD]}t||k(sJ||z}dD]5}d}|t_t ddD]}t||k(sJ||z}7dt_tdjt t sJdt_tdtdcxk(rdk(sJJtddk(sJtttsJtdjdsJtddk(sJtddk(sJtddk(sJtddk(sJtd jd!sJtd"jd#sJtd$jtd%d#sJy)&Nr?g~? @-C6?g !@i,z1.0201917073881354535e612r gkﴑ[ gQgf<9b1?y??y}?'މ^eÿy{Gz?yT})ۿOf[X@y4@>@y>/ 00 0 =  1 22 2 :==4 55 5 ;>>0 11 1 >J KK K ?  L MM M =  E FF F D 1b\Qx4  q# A8t# ## AID  BF :==b "" " BF Q<3q6 &Q && && & Q<1   s   ;>>6 77 7 !9>> ":?? !9   !9   D>  M NN N I  ! !"D EE E J  " "3'JKm#n oo or cdt_tdDcgc] }t|c}gdk(sJtdj dsJtdj dsJtdj d sJtd j d sJtt t k(sJt tt sJycc}w) Nrr) rrrrrr0r;iir/gUUUUUU?r1ggjVQr"g!a?@Fy??yHO?j6.Vֿ)r$r%rBfac2r'infrG)ns r test_fac2rTs BF"2Y 'DG '+H HH H 8;;t   9<<  8;;- .. .  ?  L MM M 9   sd   (sCcdt_d}d}t}tdgddgdk(sJtddgdgtk(sJt tdgd g|d|z|d |zz z dksJt tddgd dg|d|z|d|zz|d |zz |d|zz z dksJt d dd k(sJt d d d k(sJt dd dk(sJt ttjtd tzttz sJt ddjdsJtd dd k(sJtd d d k(sJtdd dk(sJtd ddk(sJtddd k(sJtd dd k(sJtdddk(sJtd d d k(sJtdd dk(sJtdd dk(sJtddd k(sJtddd k(sJtd dd k(sJtdddk(sJtdd dk(sJtdd dk(sJtd d d k(sJtddtk(sJtd dtk(sJtddtk(sJtdd jdsJtdd jtdz sJdtddzjdsJtttdk(sJttt tsJttdtsJttdtsJttd td tcxk(rdk(sJJttd!tk(sJttd"t k(sJtd#dtd z z jd$sJtd"d dk(sJtdd jd%sJtddjd&sJy)'Nrg:0yE>r9rr.rr/rrg@rz-255.5815971722918z-0.5119253100282322z-0.1952720278805729485rrrg@?g@@r>i"9 gUUUUUU?lI5r<grP@r=r ?@y+5?\K9ʵ?gUUUUUUտz18.157330562703710339)r$r%rA gammaprodrRabsrfjr'ffbinomialbetarDrG)hepGs r test_gamma_quotientsrhsG BF A B A bT2b' "a '' ' bVRD !S (( ( y"rd#a1ga1go5 6 == = y"R"Q(1RT71RT7?1RT7+B1QqS6+II JT QQ Q a7a<< c!9   b8r>> a7::eAaCjq) ** * "#8 9 < <=U VV V b8q== b8r>> a7b== a7a<< Aa=A   Aa=A   Ab>Q   Aa=A   Aa=B   Aa=B   Aa=A   Bq>Q   Br?a   C  '' ' D t ## # D v %% % !9>> !9   !9   2;#   A;>>$   C=  BrE "" " 43<  # #$9 :: : C=A   sd3  b   a   C=DSM .Q .. .. . D>S  D>cT !! ! R!V    N OO O S>Q   1:==   $%: ; > >?V WW Wr cVdt_tdjtdzdz sJtdjtdzdz sJtt djtdzdz sJtdjdsJtdjdsJtd jd sJtd jt d d z sJtd dk(sJtdjt ddz sJtddk(sJtddk(sJtttsJtdjdsJtdjdsJtt d djdsJtt d djdsJtt d djdsJtt ddjdsJtt ddjdsJdt_d}tt d |jddsJdt_tdddzz jdsJtdd dk(sJtddd k(sJtdd djdsJtd d!jd"sJy)#Nrrrg@r<rrr r=gs~]rrrWr.xrVKH9KHg"fD,@gOc5@g9@Ww'&l7r"z8236.5242296658162058024755079556629786895294952121891237gaU,$Db6yN@r"y9A!@.@g9B.ֿr\y?yN*i6?) r$r%zetar'rDrKr)rGrH)ims r test_zetarrs BF 7::b!eai  9<<A "" " A>$   Vg?g^s?rlrmrn)r$r%altzetar'logrDrprK)ss r test_altzetar sn BF 2;!   2;!   4=A   1:   2;$   2;&  2;$   3<= (( ( 1:==Q  1:==Qr "" " 2;>>"RV)G+ ,, , 2;?? 2 $q (( ( 2 $ )) ) 53 '! ++ + 53 '1 ,, , 4=GAJ && & A 1:==!A!H*d1g- .. . A 1:==!A!H*d1g- .. . 6?  6 77 7 3  S !! ! 6?  c "" " 3uU# $ ' ' ,, , 3ve$ % ( ( -- -r cdt_ttdk(sJdt_tdj dsJtdt zj dsJdt_td }dt_t |dz d k(sJtt d zt dkDsJtd t dkDsJtt d zdk(sJtd dk(sJy)Nrrr"r<z61.0000000000000000000000000000007888609052210118073522(z61.0000000000000000000000000000000000000148407238666182r2iz1.02363019598118e-9934rrwi )r$r%rprRr'rDr(r)vs r test_zeta_huger>s BF 9>> BF 9<<P QQ Q 2;>>R SS S BF U A BF qs8/ // / 4( +a // / x (1 ,, , 4=A   :??r cdt_t ddzz}dt_tt |dk(sJdt_tt |dk(sJdt_y)Nrrrz=2.55880492708712e+1233536161668617575553892558646631323374078r"z`2.5588049270871154960875033337384432038436330847333e+1233536161668617575553892558646631323374078)r$r%rDr(rp)as r test_zeta_negativerMs[ BF BF A BF tAw<Z ZZ Z BF tAw<} }} } BFr c * dt_d}d}|dtddcxk(rtdk(sJJ|tj dsJ|t j dsJ|t dzj d sJ|tt zj d sJ|t t z j d sJt|dj d sJ|d j dsJ|dj dsJ|dj dsJ|dj dsJ|ttk(sJ|tdk(sJtdtdk(sJ|tj dsJ|t j dsJ|tt zj dsJ|t t z j dsJdtdddtzt zzzj dsJ|ddtzt zzj dsJt|tsJt|t sJ|dj d sJ|d!j |dsJ|d"j d#sJttdtdtsJttdtdtsJttdtdt sJttdtdtsJttdtdtsJttdtdt sJttdtttsJttdtttsJttdtt t sJd$t_td%j td&d'sJdt_y)(Nrctd|SNrpsizs r z test_polygamma..X S1Xr ctd|SNrrrs r rz test_polygamma..Yrr rrguTE?g,$7@rgdGDE@y)\Z[?ÕޛD?ybabJr?N<g?Vcb?绽|=g& _rng\)c=y|=|=yF< _ _AyWw'&l7Ww'&l7y\)c-ȥ\)c-Hrg a#?g^jJ@yY/?Iyϐ=ӿ{?i@Brry˜hVKٿy_%F @@K?r|g69u@rzyY@y]vCq@˼@ry?z%0.46317279488182026118963809283042317z$2.4821070143037957102007677817351115) r$r%rdigammarDr'rafindrootrRrGrHrK)psi0psi1s r test_polygammarVs BF D D 7c!Ah ,'!* ,, ,, , 8;;* ++ + 9<<* ++ + A;>>, -- - 1:==C DD D A;>>B CC C D!    2 33 3 ;>>, -- - ;>>0 11 1   !H II I   !P QQ Q 9   9>> q:?? 8;;* ++ + 9<<* ++ + 1:==C DD D A;>>E FF F #a"R%' " " & &'P QQ Q "R%' ?  F GG G c   sd   F FF F s2bz?G GG Gr cdt_d}d}d}|djt dt dzz sJ|djt sJ|djdt dzzsJ|djt dzdz sJ|d jt dzd t zzsJ|djd tzsJd t_d tzdtdzt dzzz}dt_tddj|sJtddjt dzsJy)Nrctd|Srrrs r rz+test_polygamma_identities..rr ctd|Srrrs r rz+test_polygamma_identities..rr ctd|S)Nrrrs r rz+test_polygamma_identities..rr r=rrrr8rrWrXiJrrg?) r$r%r'eulerr~rDcatalanaperyrCr)rrpsi2rys r test_polygamma_identitiesrs- BF D D D 9<<qQx (( ( 7::uf   9<<BE "" " 7::b!eAg   :==Q7* ++ + 7::bh   BF U 1T!W9RU?"A BF q;>>!   q:==Q  r cTd}dt_dtddz}dtddz}|d|z |d|zz |dd|zzz|dd|z zz}dtzt dt dztzz}|j |sJd t_y) Nctd|Srrrs r rz'test_foxtrot_identity..rr r"rrrrr=r r)r$r%fractionrDsechrCr')rrbxys r test_foxtrot_identityrs D BF hqmA hqmA c!e tDF|#d3!9o5S!A#YGA "T#d1g+b. !!A 447N7 BFr cdt_ttdtdk(sJttdtdt zzdk(sJttdtdt zt zzdk(sJdt_ttdtdk(sJttdtdt zzd k(sJttdtdt zt zzd k(sJy) Nr<r"zf-1344100348958402765749252447726432491812.641985273160531055707095989227897753035823152397679626136483zw-0.00000000000000000189793739550804321623512073101895801993019919886375952881053090844591920308111549337295143780341396z(-0.00000000000000005225169411521692489752254721556835657523758895106315132447859377385233700094871256507814151956624433 - 0.00000000000000001813157041407010184702414110218205348527862196327980417757665282244728963891298080199341480881811613j)rz-1.34410034895841e+39z-1.89793739550804e-18z.(-5.2251694115217e-17 - 1.81315704140701e-17j))r$r%r(rrDerar r r test_polygamma_high_orderrs BF s2r{  H HH H s2rBqDy! "'` `` ` s2rBqDF{# $*^ __ _BF s2r{ 6 66 6 s2rBqDy! "&= == = s2rBqDF{# $(X XX Xr cdt_tddk(sJtddk(sJtddk(sJtdjdsJtdjd sJtd d zjd sJtd jddt dzz sJtt t k(sJtddk(sJtdjdsJy)NrrrrrZrgUUUUUU?rg_b7@rrgS@r=y@y?r]yk4q?@w?)r$r%harmonicr'r~rRr r r test_harmonicrs BF A;!   A;!   A;#   A;>>* ++ + F   4 55 5 BH  !2 33 3 C=  AaAhJ '' ' C=C   D>V ## # D>  D EE Er cdt_tddz }dt_tt |dk(sJdt_tt |dk(sJdt_y)Nrrz6.26075321389519e+12458010678r"z@6.2607532138951929201303779291707455874010420783933e+12458010678r$r%r)r(rArs r test_gamma_huge_1rsY BF F aA BF uQx=; ;; ; BF uQx=^ ^^ ^ BFr cdt_tddzdz }dt_tt |dk(sJdt_tt |dk(sJy) Nrrr<rzw1.82341134776679e+51729974693233641689901335581750771368291828240422018860515119656908623426021308685461258226190190661r"z1.82341134776678875374414910350027596939980412984e+51729974693233641689901335581750771368291828240422018860515119656908623426021308685461258226190190661rrs r test_gamma_huge_2rsm BF BG rA BF uQx= 3 44 4BF uQx= T UU Ur crdt_ddzdzddztzdz z}dt_t|}t |j dk(sJt |j d k(sJd t_t|}t |j d k(sJt |j d k(sJy) NrrPrFrrze-6.82925203918106e+2636286142112569524501781477865238132302397236429627932441916056964386399485392600zd8.54647143678418e+2636286142112569524501781477865238132302397236429627932441916056964386399485392600r"z-6.8292520391810548460682736226799637356016538421817e+2636286142112569524501781477865238132302397236429627932441916056964386399485392600z8.5464714367841748507479306948130687511711420234015e+2636286142112569524501781477865238132302397236429627932441916056964386399485392600r$r%rarAr(realimagrrs r test_gamma_huge_3rs BF B! b"fQhl"A BF aA qvv;  !! ! qvv;    BF aA qvv; C DD D qvv; B CC Cr cd}dt_tt|dk(sJdt_tt|dk(sJdt_y)Ny@v@rz2(8.95783268539713e+5164 - 1.94678798329735e+5164j)r"zx(8.9578326853971339570292952697675570822206567327092e+5164 - 1.9467879832973509568895402139429643650329524144794e+5164j))r$r%r(rArs r test_gamma_huge_4rsZA BF uQx=< == = BF uQx=   BFr cvdt_ddztzdz }dt_t|}t |j dk(sJt |j dk(sJdt_t|}t |j d k(sJt |j d k(sJdt_y) NrrrurrzO-3.27753899634941e-227396058973640224580963937571892628368354580620654233316839zN-7.1519888950416e-227396058973640224580963937571892628368354580620654233316841r"zr-3.2775389963494132168950056995974690946983219123935e-227396058973640224580963937571892628368354580620654233316839zr-7.1519888950415979749736749222530209713136588885897e-227396058973640224580963937571892628368354580620654233316841rrs r test_gamma_huge_5rs BF B QA BF aA qvv;k kk k qvv;j jj j BF aA qvv; - .. . qvv; - .. .BFr crdt_dttddzz z}dt_t |}t |j dk(sJt |jdk(sJdt_t |}t |j dk(sJt |jd k(sJy) Nr<rrrrz2.0z1.84556867019693e-1000r"z91.8455686701969342787869758198351951379156813281202e-1000)r$r%rar)rAr(rr)rrs r test_gamma_huge_7rs BF Ac"gtmOA BF aA qvv;%   qvv;2 22 2 BF aA qvv;%   qvv;V VV Vr cHdt_tdjtsJdt_tdjdsJtdjdsJtdjd sJtd jd sJdt_tdjd sJtdjd sJtdjdsJtd jdsJy)Nrrrz -0.07281584548367672486058637587rz!-0.009690363192872318484530386035rz 0.002053834420303345866160046543rz 0.002325370065467300057468170178g*gF|؃gT3`?g; c?)r$r% stieltjesr'rr r r test_stieltjesrs BF Q> > Q ?? ? Q> > Q> > BF Q  3 44 4 1c?  3 44 4 1d   6 77 7 1d   N OO O 1c?  J KK K 1d   6 77 7 1c$q'k " % %&T UU U 2s   t $$ $ 2t    ** * 1c?  3 44 4 2s   u %% % 2r?  g && & 3qs7C!G $ ' '(X YY Y 4?  d4j )) ) 4    ~ .. . 3+ , / /0B CC C 3    M NN N 5"  Q -- - 5$  " "4 (( (5 GsN#cdt_tdddzjdsJtdddzjdsJtdd jd sJtd d jd sJt dddk(sJt dddk(sJt dddk(sJt dt dk(sJt dt dk(sJt t dtsJt dddk(sJt dddk(sJt dt t k(sJt dt t k(sJt t dtsJt dd dk(sJt dddk(sJt ddk(sJt dddk(sJt dddk(sJt dddk(sJt ddjdsJdt_t dddk(sJt dddk(sJt ddk(sJdt_t dddk(sJt d djttd zsJy)!Nrrrrgo?rg8?ry@y )fɳ?@rg @g* 0(@rrriiiii0r"z3.4746902914629720259e74rlv reh@3 VBFv'r.y UR|s&Trl.e0u&Y:j+I<1 |6vJbJlrw_3dtL#m, `3 .v(*A=Vy4@y@j #MB(3VjwBr=) r$r%polyexpr'bellrRrGrHsincrDr r r test_bell_polyexprUsb BF Ae V # ' ' 66 6 Ae V # ' ' 55 5 1R=  I JJ J 2c?  2 33 3 !9>> !9>> !9>> #;!   3$<1   a   !9>> !9>> #;#   3$> 8v   1:   B<8 ## # #;* $$ $ 3>$r#v, '' 'r cdt_tdjdsJtdjdsJtdtk(sJttdk(sJt tt sJy)Nrry2k?-DT! @rgMҵ?rr)r$r% primezetar'rRrGrHr r r test_primezetar{so BF S>  K LL L QQ   3  r cdt_tdjdsJtdjdsJtddjd sJtd djd sJtd d jd sJtdd jdsJd"d}|tjddsJ|tjddsJ|tjddsJ|tjdd dsJ|tjd dsJ|tjd d d sJ|tjdd dsJ|tjdddsJ|tjdddsJ|tjdddsJ|tjdd dsJ|tjdddsJ|tjddd sJ|tjddd!sJy)#Nry?j@y&6#+?–@y?j@y5g?7-m^I?y?.Ar derivativey@V%.aes1Q3x#c!f*$$r y?jy&6#+?–y?jy5g?7-m^Iy’z?ی^p@y. F'@ p 1y$P??y,ς?%@lUS>ryԂ]5MԺ$k@y!&@\_Ury(1Sg@y-6 ڿ+~y "虺m"'= UA@yِ@[@MomykLVy"DЫ@)gư>)r$r%rpr'fp)r's r test_rs_zetarsE BF     N OO O   !N OO O   + . ./] ^^ ^  q ) , ,-[ \\ \  q ) , ,-^ __ _   + . ./] ^^ ^% bggk"$R SS S bggl#%R SS S bggh!Q RR R bggh1g-/] ^^ ^ bggfoO PP P bggfg+-^ __ _ bggkag02_ `` ` bggkag02a bb b bggkag02` aa a bggkag02a bb b bgglqg13c dd d bgglqg13` aa a bgglqg13a bb b bgglqg13a bb br cdt_tdjdsJtddjdsJtddjdsJtdtzd jd sJy) Nrr!g1/[㳄@rrgKKrgFhnqryQ1@a߻w)r$r%siegelzr'rar r r test_siegelzrsz BF 6?  4 55 5 6a ( + +,C DD D 6a ( + +,C DD D 6!8q ) , ,-[ \\ \r cdt_tddd}tddd}tddd}t|j dsJt|j d sJt|}|j j d sJ|j j d sJd t_tdd d}tddd}tddd}t|j dsJt|j dsJt|}|j j dsJ|j j dsJy)Nrrz1e-10T)exactz-1e-10z1e-10jg# _Bga_gox?g _rz1e-50z-1e-50z1e-50j1e50z-1e50z40.57721566490153286060651209008240243104215933593992)r$r%faddrpr'rr)s1s2s3rs r test_zeta_near_1rs+ BF a %B a &B a &B 8;;7 88 8 8;;6 77 7 RA 6699+ ,, , 66993 44 4 BF a %B a &B a &B 8;;v   8;;w   RA 6699K LL L 6699W  r N)#mpmath mpmath.libmprrrr r,r6rNrTrhrrrrrrrrrrrrrrrrrrrrrrrrrr r r rs;;3+Ze>N%pN/Xb!QF.< +ZG  "  Y F UC( " W >**)<$(L!c6]r