a Z"d?S)@)*)round_up from_float mpf_zeta_intcCstddtdksJdS)Nr )rrrrk/Users/vegardjervell/Documents/master/model/venv/lib/python3.9/site-packages/mpmath/tests/test_gammazeta.pytest_zeta_int_bugsr cCstddksJtddks Jtddks0Jtddks@Jtd d ksPJtd dks`Jtd d kspJtdd ksJtddksJtddksJtddksJtd\}}|ddksJ|dksJtd\}}|ddksJ|dksJdt_tddksJtddks,Jtdds@JtddksRJtd dsfJtd dksxJtd d sJttdd!ksJttd"d#ksJttd$d%ksJttd&d'ksJd(t_ttdd)ksJttd"d*ksJttd$d+ks,Jttd&d,ksBJdt_tdtd d-sdJd.t_tdtd d-sJdt_dS)/Nr)r r )r )r )rr )r r)r *r)rB )iMi )i[id( is i lkL3iqrgUUUUUU?ggaa?z0.0757575757575758z7.62772793964344e+267順z-5.82229431461335e+376755iz1.19570355039953e+6767525842z50.075757575757575757575757575757575757575757575757576z87.6277279396434392486994969020496121553385863373331e+267z<-5.8222943146133508236497045360612887555320691004308e+376755z>1.1957035503995297272263047884604346914602088317782e+676752584riP)ZbernfracmpdpsZ bernoulliaestrmpf)pqrrr test_bernoullisN   r&cCstdt_tdddsJtddds.JtdddsBJtdddsVJtdddsjJtddds~JtdddsJtddd sJtdddsJtddd sJtdddsJtddd sJtddd s Jtddd s"Jtddds8JtdddsNJtddd sdJtddd szJtdddsJtdddsJtdddsJtdddsJtdddsJtdddsJtdddsJtddds*Jtddds@JtdddsVJtdddslJtdddsJtdddsJtdddsJtdddsJtdddsJtdddsJtdddsJtdddsJtddds2JtdddsHJtddds^JtdddstJtdddsJtd dd!sJtd ddsJtd ddsJtd dd"sJtd ddsJtd dd sJtddds$Jtddds:JtdddsPJtdddsfJtddds|JtdddsJtdddsJtddd sJtdddsJtddd sJtdddsJtddd sJtddds,JtdddsBJtddd#sXJtddd$snJtdddsJtdddsJtddd%sJtddd sJtdddsJtddd&sJtddd#sJtddd'sJtddds4JtdddsJJtddd(s`Jtddd)svJtdddsJtdddsJtdddsJtddd sJtdddsJtddd*sJtdddsJtddd s&Jtddds@ycCsdt_ddtdDgdks$Jtdds6Jtdd sHJtd d sZJtd d slJtttks|Jttt sJdS)NrcSsg|] }t|qSr)fac2.0nrrr ztest_fac2..r) r r r rrr0r1iir(gUUUUUU?r*ggjVQrg!a?@Fy??yHO?j6.Vֿ)rr r6r?r!infr:rrrr test_fac2srGcCsdt_d}d}t}tdgddgdks*JtddgdgtksBJttdgd g|d||d |dkstJttddgd dg|d||d||d ||d|dksJtd dd ksJtd d d ksJtdd dksJttttd tttsJtddds4Jt d dd ksHJt d d d ks\Jt dd dkspJt d ddksJt ddd ksJt d dd ksJt dddksJt d d d ksJt dd dksJt dd dksJt ddd ksJt ddd ks$Jt d dd ks8Jt dddksLJt dd dks`Jt dd dkstJt d d d ksJt ddtksJt d dtksJt ddtksJt dd dsJt dd t dsJdt dddsJt ttdks"Jt t t ts8Jt t dtsLJt t dts`Jt td t d tkrdksnJt td!tksJt td"t ksJt d#dtd d$sJt d"d dksJt dd d%sJt ddd&sJdS)'Nrg:0yE>r/r r'rr(rr g@r z-255.5815971722918z-0.5119253100282322z-0.1952720278805729485rrrg@?g@@r4i"9 gUUUUUU?lI5r2grP@r3r?@y+5?\K9ʵ?gUUUUUUտz18.157330562703710339)rr r5Z gammaprodrFabsrfjr!ffZbinomialbetar8r:)hepGrrr test_gamma_quotientss^2N&*rXcCsRdt_tdtdds Jtdtdds:JttdtddsXJtddsjJtdds|Jtd d sJtd td d sJtd dksJtdtddsJtddksJtddksJtttsJtddsJtdds.Jttd ddsHJttd ddsbJttd dds|JttdddsJttdddsJdt_d}ttd |ddsJdt_tddddsJtdd dks Jtddd ks Jtdd dds8Jtd d!d"sNJdS)#Nrr rg@r2r rrr3gs~]r rrIr'xrHKH9KHg"fD,@gOc5@g9@Ww'&l7rz8236.5242296658162058024755079556629786895294952121891237gaU,$Db6yN@r"y9A!@.@g9B.ֿrNy?yN*i6?) rr zetar!r8r<r#r:r;)Zimrrr test_zetas:r`cCs0dt_tddksJtddks&Jtddks6JtddksFJtddksVJtd d ksfJtd dksvJtd d ksJtdtdsJtdtddsJtddtddsJtddksJtddddksJtddddksJtddddks$Jtddddks:JtdtdksPJd}t|ddd|t|s|Jd}t|ddd|t|sJtddsJtd d!sJtd"dsJtd#dsJttd"d$dsJttd#d$ds,JdS)%NrrIrrHrZr3r r.r'gr(igz\Ar r rrIijhr<droundingur+@r@@y@ Yg@o.>Vg?g^s?r[r\r])rr altzetar!logr8r_r<)srrr test_altzeta s8((rmcCsdt_ttdksJdt_tdds.JtdtdsDJdt_td }dt_t|dd kslJttd td dksJtd td dksJttd dksJtd dksJdS)Nrr rr2z61.0000000000000000000000000000007888609052210118073522(z61.0000000000000000000000000000000000000148407238666182r+iz1.02363019598118e-9934rrdi )rr r_rFr!r8r"r)vrrr test_zeta_huge>srpcCsRdt_t dd}dt_tt|dks.Jdt_tt|dksHJdt_dS)Nrrnrz=2.55880492708712e+1233536161668617575553892558646631323374078rz`2.5588049270871154960875033337384432038436330847333e+1233536161668617575553892558646631323374078)rr r8r"r_arrr test_zeta_negativeMsrtcCs>dt_dd}dd}|dtddkr8tdks>nJ|tdsPJ|t dsdJ|t d d s|J|ttd sJ|t td sJt|d d sJ|ddsJ|ddsJ|ddsJ|ddsJ|ttksJ|tdks,Jtdtdks@J|tdsTJ|t dsjJ|ttdsJ|t tdsJdtdd dttdsJ|d dttdsJt |t sJt |t sJ|d d!sJ|d"|ds4J|d#d$sHJt tdt dtsbJt tdt dt s|Jt tdt dt sJt tdt d tsJt tdt d t sJt tdt d t sJt tdt ttsJt tdt t t sJt tdt t t s:JdS)%NrcSs td|SNrpsizrrr XrDz test_polygamma..cSs td|SNr rvrxrrr rzYrDrrguTE?g,$7@r gdGDE@y)\Z[?ÕޛD?ybabJr?N<g?Vcb?绽|=g& _r]g\)c=y|=|=yF< _ _AyWw'&l7Ww'&l7y\)c-ȥ\)c-Hr g a#?g^jJ@yY/?Iyϐ=ӿ{?i@Brry˜hVKٿy_%F @@K?rig69u@rgyY@y]vCq@˼@) rr rwZdigammar8r!rRZfindrootrFr:r;r<)psi0psi1rrr test_polygammaVsJ(& rcCs6dt_ttdtdksJttdtdks2JdS)Nr2rzf0.9772133079420067332920694864061823436408346099943256380095232865318105924777141317302075654362928734rzf-12.98876181434889529310283769414222588307175962213707170773803550518307617769657562747174101900659238)rr r"rwr8rrrr test_polygamma_high_precsrcCsdt_dd}dd}dd}|dt dtds>J|dt sRJ|ddtdslJ|dtdd sJ|d tdd tsJ|dd tsJd t_dtdtdtd}dt_t dd|sJt ddtdsJdS)NrcSs td|Srurvrxrrr rzrDz+test_polygamma_identities..cSs td|Sr{rvrxrrr rzrDcSs td|S)Nr rvrxrrr rzrDr3r r rr.rrIrJiJrrg?) rr r!eulerrkr8catalanZaperyr7rw)r}r~Zpsi2rfrrr test_polygamma_identitiess rcCsdd}dt_dtdd}dtdd}|d| |d ||dd||dd|}dttdtdt}||sJd t_dS) NcSs td|Srurvrxrrr rzrDz'test_foxtrot_identity..rr r rr r3rr)rr fractionr8Zsechr7r!)r}rsbxyrrr test_foxtrot_identitys:rcCsdt_ttdtdksJttdtdtdks:Jttdtdttdks\Jdt_ttdtdksxJttdtdtd ksJttdtdttd ksJdS) Nr2rzf-1344100348958402765749252447726432491812.641985273160531055707095989227897753035823152397679626136483zw-0.00000000000000000189793739550804321623512073101895801993019919886375952881053090844591920308111549337295143780341396z(-0.00000000000000005225169411521692489752254721556835657523758895106315132447859377385233700094871256507814151956624433 - 0.00000000000000001813157041407010184702414110218205348527862196327980417757665282244728963891298080199341480881811613j)rz-1.34410034895841e+39z-1.89793739550804e-18z.(-5.2251694115217e-17 - 1.81315704140701e-17j))rr r"rwr8erRrrrr test_polygamma_high_orders"rcCsdt_tddksJtddks&Jtddks6JtddsHJtdd sZJtd d d spJtd ddtdsJtttksJtddksJtddsJdS)Nrrr r rLrgUUUUUU?rg_b7@rrgS@r3y@y?rOyk4q?@w?)rr Zharmonicr!rkrFrrrr test_harmonicsrcCsPdt_tdd}dt_tt|dks,Jdt_tt|dksFJdt_dS)Nrrz6.26075321389519e+12458010678rz@6.2607532138951929201303779291707455874010420783933e+12458010678rr r#r"r5rrrr test_gamma_huge_1s rcCsNdt_tddd}dt_tt|dks0Jdt_tt|dksJJdS) Nrrr2rzw1.82341134776679e+51729974693233641689901335581750771368291828240422018860515119656908623426021308685461258226190190661rz1.82341134776678875374414910350027596939980412984e+51729974693233641689901335581750771368291828240422018860515119656908623426021308685461258226190190661rrrrr test_gamma_huge_2s  rcCsdt_dddddtd}dt_t|}t|jdksBJt|jd ksTJd t_t|}t|jd kstJt|jd ksJdS) NrrPrFrrze-6.82925203918106e+2636286142112569524501781477865238132302397236429627932441916056964386399485392600zd8.54647143678418e+2636286142112569524501781477865238132302397236429627932441916056964386399485392600rz-6.8292520391810548460682736226799637356016538421817e+2636286142112569524501781477865238132302397236429627932441916056964386399485392600z8.5464714367841748507479306948130687511711420234015e+2636286142112569524501781477865238132302397236429627932441916056964386399485392600rr rRr5r"realimagrrrrr test_gamma_huge_3s$rcCsBd}dt_tt|dksJdt_tt|dks8Jdt_dS)Ny@v@rz2(8.95783268539713e+5164 - 1.94678798329735e+5164j)rzx(8.9578326853971339570292952697675570822206567327092e+5164 - 1.9467879832973509568895402139429643650329524144794e+5164j))rr r"r5rrrr test_gamma_huge_4s  rcCsdt_ddtd}dt_t|}t|jdks6Jt|jdksHJdt_t|}t|jd kshJt|jd kszJdt_dS) NrrrbrrzO-3.27753899634941e-227396058973640224580963937571892628368354580620654233316839zN-7.1519888950416e-227396058973640224580963937571892628368354580620654233316841rzr-3.2775389963494132168950056995974690946983219123935e-227396058973640224580963937571892628368354580620654233316839zr-7.1519888950415979749736749222530209713136588885897e-227396058973640224580963937571892628368354580620654233316841rrrrr test_gamma_huge_5srcCsdt_dttdd}dt_t|}t|jdks:Jt|jdksLJdt_t|}t|jdkslJt|jd ks~JdS) Nr2rrrrz2.0z1.84556867019693e-1000rz91.8455686701969342787869758198351951379156813281202e-1000)rr rRr#r5r"rr)rsrrrr test_gamma_huge_7srcCsdt_tdt sJdt_tdds2JtddsDJtdd sVJtd d shJdt_tdd sJtdd sJtddsJtd dsJdS)Nrrr z -0.07281584548367672486058637587r z!-0.009690363192872318484530386035rz 0.002053834420303345866160046543rz 0.002325370065467300057468170178g*gF|؃gT3`?g; c?)rr Z stieltjesr!rrrrr test_stieltjessrcCs|dt_tdtdkr"dks(nJddtdDgdksFJttdd ksZJtttsjJtttszJtttsJtt t ksJtt t ksJtt t ksJttt sJtd  d sJtd  d sJddtdDgdksJddtdddDgdks4Jtd}|dkrRt |t sVJtd}|dkrtt |t sxJdS)Nrrr cSsg|] }t|qSr)superfac)rAr=rrr rC!rDz test_barnesg..r)r r r ri ii{l2ftrz3.24570818422368e+1177245gffffff?gd?y@@yg ƿ1?cSsg|] }t|qSrhyperfacr@rrr rC,rDr)r r rlli\&lcSsg|] }t|qSrrr@rrr rC-rD)r r r rHrriyrH) rr Zbarnesgr6r"rr:r;rrFr! isinstancer<rrrrr test_barnesgs&" $rcCsdt_dddD}|D]"}td|td| sJq|D] }td||d|s@Jq@|D]$}td||d|dsfJqf|D],}td ||d|d|d sJq|D]8}td ||dd ||dd|d sJqtd d dsJtd dds(Jtddds>JtdddsTJtdddsjJtdddsJtdddsJtddttdsJtd ddsJtd ddsJtd ddsJtd dds Jtd d d!s"Jtd"td d"tdd#sHJtd$dtd$sbJtd$dtd$ s~Jtd%d&d'sJtd%d(d)sJdS)*NrcSsg|] }t|qSr)Z mpmathify)rAryrrr rC5rDz test_polylog..)rr3gGz?rrrHy?rhr rr r rIrr'rry[-F@U޴rg!fG?gV?ggy?y#]ſ`j?g?y8n+Yd?^n)ӿggᾁyp/M?[]?ig jߊg=.g?i2rHg-C6zr3yB8$P?g?rhrLgtHgWcFr,yQ K1Q ) rr Zpolylogr!rkr7rRr_rj)Zzsryrrr test_polylog3s2(&*2>&rcCs8dt_tddddsJtdddds6Jtdd d sJJtd d d s^JtdddkspJtdddksJtdddksJtdtdksJtdt dksJttdtsJtdddksJtdddksJtdttksJtdt t ksJttdts0Jtdd dksDJtdddksXJtddksjJtdddks~JtdddksJtdddksJtdddsJdt_tdddksJtdddksJtddksJdt_tdddksJtd dtt d s4JdS)!Nrrr|rgo?r g8?ry@y )fɳ?@r g @g* 0(@r rriiiii0rz3.4746902914629720259e74rlv reh@3 VBFv'r.y UR|s&Trnl.e0u&Y:j+I<1 |6vJbJlrw_3dtL#m, `3 .v(*A=Vy4@y@j #MB(3VjwBr3) rr Zpolyexpr!ZbellrFr:r;Zsincr8rrrr test_bell_polyexpOs<rcCs^dt_tddsJtdds*Jtdtks:JttdksJJtttsZJdS)Nrry2k?-DT! @rgMҵ?r r)rr Z primezetar!rFr:r;rrrr test_primezetaus rcCsdt_tddsJtdds*Jtdddd s@Jtd ddd sVJtd d dd slJtdd ddsJd$dd}|tddsJ|tddsJ|tddsJ|tjdd ddsJ|td dsJ|tjd d dd sJ|tjdd dds(J|tjddddsBJ|tjdddds\J|tjddddsvJ|tjdd dd sJ|tjdddd!sJ|tjdddd"sJ|tjdddd#sJdS)%Nry?j@y&6#+?–@y?j@y5g?7-m^I?y?.ArZ derivativey@V%cSst|||t|kS)N)rP)rrZtolrrr r!sztest_rs_zeta..aey?jy&6#+?–y?jy5g?7-m^Iy’z?ی^p@y. F'@ p 1y$P??y,ς?%@lUS>r yԂ]5MԺ$k@y!&@\_Ury(1Sg@y-6 ڿ+~y "虺m"'= UA@yِ@[@MomykLVy"DЫ@)r)rr r_r!fp)r!rrr test_rs_zeta}s, rcCsbdt_tddsJtdddds.JtddddsDJtdtd dd s^JdS) Nrrg1/[㳄@r rgKKrgFhnqr yQ1@a߻w)rr Zsiegelzr!rRrrrr test_siegelzs rcCsdt_tdddd}tdddd}tdddd}t|dsBJt|d sTJt|}|jd slJ|jd s|Jd t_tdd dd}tdddd}tdddd}t|dsJt|dsJt|}|jdsJ|jdsJdS)Nrr z1e-10T)exactz-1e-10z1e-10jg# _Bga_gox?g _rz1e-50z-1e-50z1e-50jZ1e50z-1e50z40.57721566490153286060651209008240243104215933593992)rr Zfaddr_r!rr)s1s2Zs3ryrrr test_zeta_near_1s$rN)#ZmpmathZ mpmath.libmprrrr r&r-r>rGrXr`rmrprtrrrrrrrrrrrrrrrrrrrrrrrr s>-g' 1# )      &