+
    im                         ^ RI 5 ^ RIHtHtHt R tR tR tR tR t	R t
R	 tR
 tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR t R t!R  t"R!# )"    )*)round_up
from_floatmpf_zeta_intc                  >    \        ^ ^
4      \        R4      8X  g   Q hR# )r   N      )r   r        {/Users/tonyclaw/.openclaw/workspace/skills/math-calculator/venv/lib/python3.14/site-packages/mpmath/tests/test_gammazeta.pytest_zeta_int_bugr      s    2*T"2222r
   c                  l   \        ^ 4      R8X  g   Q h\        ^4      R8X  g   Q h\        ^4      R8X  g   Q h\        ^4      R8X  g   Q h\        ^4      R8X  g   Q h\        ^4      R8X  g   Q h\        ^4      R8X  g   Q h\        ^4      R8X  g   Q h\        ^
4      R8X  g   Q h\        ^4      R8X  g   Q h\        ^4      R8X  g   Q h\        ^4      w  rV R,          R8X  g   Q hVR8X  g   Q h\        R4      w  rV R,          R8X  g   Q hVR8X  g   Q h^\        n        \        ^ 4      ^8X  g   Q h\        ^4      R8X  g   Q h\        ^4      P	                  R4      '       g   Q h\        ^4      ^ 8X  g   Q h\        ^4      P	                  R4      '       g   Q h\        ^4      ^ 8X  g   Q h\        ^4      P	                  R4      '       g   Q h\        \        ^
4      4      R8X  g   Q h\        \        ^4      4      R8X  g   Q h\        \        R4      4      R8X  g   Q h\        \        R4      4      R	8X  g   Q h^2\        n        \        \        ^
4      4      R
8X  g   Q h\        \        ^4      4      R8X  g   Q h\        \        R4      4      R8X  g   Q h\        \        R4      4      R8X  g   Q hR\        n        \        ^
4      P	                  \        ^4      ^B,          4      '       g   Q hR\        n        \        ^
4      P	                  \        ^4      ^B,          4      '       g   Q h^\        n        R# ) r   is	i	   l   kL3 iqz0.0757575757575758z7.62772793964344e+267z-5.82229431461335e+376755z1.19570355039953e+676752584z50.075757575757575757575757575757575757575757575757576z87.6277279396434392486994969020496121553385863373331e+267z<-5.8222943146133508236497045360612887555320691004308e+376755z>1.1957035503995297272263047884604346914602088317782e+676752584iP  N)   r   )   )r      )r   r   )r      )r   *   )   B   )iMi
  )i[  i      d(	 r   gUUUUUU?ggaa?順 i)bernfracmpdps	bernoulliaestrmpf)pqs     r   test_bernoullir"      s   A;%A;&   A;%A;%A;'!!!A;%A;&   A;'!!!B<6!!!B<;&&&B<;&&&C=DAv:""";;D>DAv:###	>>BFQ<1Q<4Q<??4    Q<1Q<??6""""Q<1Q<??5!!!!y}!5555y~"9999y $????y!"&CCCCBFy}!XXXXy~"\\\\y $bbbby!"&ffffBFR=CF2I&&&&BFR=CF2I&&&&BFr
   c                  .   ^\         n        \        ^ R*4      P                  ^4      '       g   Q h\        ^ ^ 4      P                  ^4      '       g   Q h\        ^ R4      P                  ^4      '       g   Q h\        ^ R4      P                  ^4      '       g   Q h\        ^ ^4      P                  ^4      '       g   Q h\        ^ ^4      P                  ^4      '       g   Q h\        ^R*4      P                  R4      '       g   Q h\        ^^ 4      P                  R4      '       g   Q h\        ^R4      P                  ^ 4      '       g   Q h\        ^R4      P                  R4      '       g   Q h\        ^^4      P                  R4      '       g   Q h\        ^^4      P                  R4      '       g   Q h\        ^R*4      P                  R4      '       g   Q h\        ^^ 4      P                  R4      '       g   Q h\        ^R4      P                  R	4      '       g   Q h\        ^R4      P                  R
4      '       g   Q h\        ^^4      P                  R4      '       g   Q h\        ^^4      P                  R4      '       g   Q h\        ^R*4      P                  R+4      '       g   Q h\        ^^ 4      P                  ^ 4      '       g   Q h\        ^R4      P                  ^ 4      '       g   Q h\        ^R4      P                  R4      '       g   Q h\        ^^4      P                  ^ 4      '       g   Q h\        ^^4      P                  ^4      '       g   Q h\        ^R*4      P                  R4      '       g   Q h\        ^^ 4      P                  R4      '       g   Q h\        ^R4      P                  R4      '       g   Q h\        ^R4      P                  R4      '       g   Q h\        ^^4      P                  R4      '       g   Q h\        ^^4      P                  R4      '       g   Q h\        ^R*4      P                  R,4      '       g   Q h\        ^^ 4      P                  ^ 4      '       g   Q h\        ^R4      P                  ^ 4      '       g   Q h\        ^R4      P                  R4      '       g   Q h\        ^^4      P                  ^ 4      '       g   Q h\        ^^4      P                  ^4      '       g   Q h\        ^
R*4      P                  R4      '       g   Q h\        ^
^ 4      P                  R4      '       g   Q h\        ^
R4      P                  R4      '       g   Q h\        ^
R4      P                  R4      '       g   Q h\        ^
^4      P                  R4      '       g   Q h\        ^
^4      P                  R4      '       g   Q h\        ^R*4      P                  R-4      '       g   Q h\        ^^ 4      P                  ^ 4      '       g   Q h\        ^R4      P                  ^ 4      '       g   Q h\        ^R4      P                  R4      '       g   Q h\        ^^4      P                  ^ 4      '       g   Q h\        ^^4      P                  ^4      '       g   Q h\	        ^ R*4      P                  ^4      '       g   Q h\	        ^ ^ 4      P                  ^4      '       g   Q h\	        ^ R4      P                  ^4      '       g   Q h\	        ^ R4      P                  ^4      '       g   Q h\	        ^ ^4      P                  ^4      '       g   Q h\	        ^ ^4      P                  ^4      '       g   Q h\	        ^R*4      P                  R4      '       g   Q h\	        ^^ 4      P                  R4      '       g   Q h\	        ^R4      P                  ^ 4      '       g   Q h\	        ^R4      P                  R4      '       g   Q h\	        ^^4      P                  R4      '       g   Q h\	        ^^4      P                  R4      '       g   Q h\	        ^R*4      P                  ^4      '       g   Q h\	        ^^ 4      P                  ^ 4      '       g   Q h\	        ^R4      P                  R4      '       g   Q h\	        ^R4      P                  R4      '       g   Q h\	        ^^4      P                  ^ 4      '       g   Q h\	        ^^4      P                  ^4      '       g   Q h\	        ^R*4      P                  R4      '       g   Q h\	        ^^ 4      P                  R4      '       g   Q h\	        ^R4      P                  ^ 4      '       g   Q h\	        ^R4      P                  R4      '       g   Q h\	        ^^4      P                  R4      '       g   Q h\	        ^^4      P                  R4      '       g   Q h\	        ^R*4      P                  ^4      '       g   Q h\	        ^^ 4      P                  ^ 4      '       g   Q h\	        ^R4      P                  R4      '       g   Q h\	        ^R4      P                  R4      '       g   Q h\	        ^^4      P                  ^ 4      '       g   Q h\	        ^^4      P                  ^4      '       g   Q h\	        ^R*4      P                  R4      '       g   Q h\	        ^^ 4      P                  R4      '       g   Q h\	        ^R4      P                  ^ 4      '       g   Q h\	        ^R4      P                  R4      '       g   Q h\	        ^^4      P                  R4      '       g   Q h\	        ^^4      P                  R4      '       g   Q h\	        ^
R*4      P                  ^4      '       g   Q h\	        ^
^ 4      P                  ^ 4      '       g   Q h\	        ^
R4      P                  R4      '       g   Q h\	        ^
R4      P                  R4      '       g   Q h\	        ^
^4      P                  ^ 4      '       g   Q h\	        ^
^4      P                  ^4      '       g   Q h\	        ^R*4      P                  R 4      '       g   Q h\	        ^^ 4      P                  R!4      '       g   Q h\	        ^R4      P                  ^ 4      '       g   Q h\	        ^R4      P                  R"4      '       g   Q h\	        ^^4      P                  R#4      '       g   Q h\	        ^^4      P                  R$4      '       g   Q h\        R%R%4      P                  R&4      '       g   Q h\        ^R'4      P                  R(4      '       g   Q h\	        ^R'4      P                  R.4      '       g   Q hR)# )/   z1/2z3/4z-3/2z-1/2z1/4z3/2z13/6z1/6z-1/12z-1/48z-3/64z119/30z-1/30z7/240z7/3840z25/1024z665/66z5/66z-2555/33792z-2555/34603008z-555731/4194304z-1/4z-3/16z-9/4z-11/64z9/4z5/16z57/256z361/1024z-50521/1024z-36581523/1048576z-699/4z691/4z-512343611/4194304z-691/4z699/4'  z5.8196915936323387117e+39999g     1@gJrNr   g,Qr%)r   r   bernpolyr   	eulerpolyr	   r
   r   test_bernpoly_eulerpolyr+   4   sl   BFAb>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2222Bq>V$$$$Bq>X&&&&Br?c""""Bq>QBu  ####Bu  !2333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!!-0000R!!"56666R?a    R?a    Rx((((R?g&&&&R!!!$$$$R!!"67777R?h''''R?g&&&&E% ##$BCCCCC  !:;;;;S!!"<====r
   c                     ^\         n        \        R4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h^p \	        ^^4       F  p\        V4      V 8X  g   Q hW,          p K  	  R F=  p^p V\         n        \	        ^^i4       F  p\        V4      V 8X  g   Q hW,          p K  	  K?  	  ^d\         n        \        R4      P                  \        \        4      4      '       g   Q h^\         n        \        ^ 4      \        ^ 4      u;8X  d	   ^8X  g   Q h Q h\        ^4      ^8X  g   Q h\        \        \        4      4      '       g   Q h\        R	4      P                  R
4      '       g   Q h\        ^ 4      ^ 8X  g   Q h\        R4      ^ 8X  g   Q h\        ^4      R8X  g   Q h\        ^4      R8X  g   Q h\        R4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h\        R4      P                  \        RR4      4      '       g   Q hR# )r$         ?g~?@-C6?g	!@i,  z1.0201917073881354535e612      ?gf<9b1?L  z4.8579168073569433667e2866g      ?1e10000z 2.302485092994045684017991e100041e10000jz!-1.570796326794896619231322e10000Nr   gkﴑ[gQy      ?      ?y}?'މ^eÿr   y      {Gz?yT})ۿOf[X@y      4@      >@y<lsfsuWA)   iX  y       @       @yaaJ
!z¶%@)r   r   gammar   rangesqrtpi	factorialfacisnannanrgammaloggammampc)factir   s      r   
test_gammarA      s   BF;>>/0000=12222:==45555;>>01111<??45555 ;>>JKKKK?LMMMM=EFFFF D1b\Qx4	  q#A8t###ID   BF:==b""""BFQ<3q6&Q&&&&&Q<1s;>>67777!9>>":??!9!9D>MNNNNI!!"DEEEEJ""3'JKm#noooor
   c                     ^\         n        \        ^
4       U u. uF  p \        V 4      NK  	  up . RO8X  g   Q h\        R4      P	                  R4      '       g   Q h\        R4      P	                  R4      '       g   Q h\        ^24      P	                  R4      '       g   Q h\        R4      P	                  R	4      '       g   Q h\        \
        4      \
        8X  g   Q h\        \        \
        ) 4      4      '       g   Q hR# u up i )
r$   g!a?@FN)
r   r   r         r$   0   i   i  i  r'   gUUUUUU?r(   ggjVQy      ?      ?yHO?j6.Vֿ)r   r   r5   fac2r   infr:   )ns    r   	test_fac2rJ      s    BF"2Y'YDGY'+HHHH8;;t9<<    8;;-....	?LMMMM9sd (s   C>c                  P
   ^\         n        Rp Rp\        p\        R.RR.4      ^ 8X  g   Q h\        R^ .R.4      \        8X  g   Q h\        \        R.R.4      V! RV ,           4      V! RV ,           4      ,          ,
          4      R8  g   Q h\        \        RR.R^ .4      V! RV ,           4      V! RV ,           4      ,          V! RV ,           4      ,          V! ^ V ,           4      ,          ,
          4      R8  g   Q h\        ^^ 4      ^8X  g   Q h\        R^4      R8X  g   Q h\        R^4      ^8X  g   Q h\        \        \        4      P                  \        ^\        ,          4      \        \        4      ,          4      '       g   Q h\        RR4      P                  R4      '       g   Q h\        R^ 4      ^8X  g   Q h\        R^4      R8X  g   Q h\        ^^4      ^8X  g   Q h\        ^^4      ^ 8X  g   Q h\        ^ ^ 4      ^8X  g   Q h\        ^^ 4      ^8X  g   Q h\        ^ R4      ^ 8X  g   Q h\        ^^4      ^8X  g   Q h\        ^^4      ^
8X  g   Q h\        ^^4      ^
8X  g   Q h\        ^^4      ^8X  g   Q h\        R^ 4      ^8X  g   Q h\        RR4      ^8X  g   Q h\        RR4      R	8X  g   Q h\        R
^4      R
8X  g   Q h\        R
^4      R8X  g   Q h\        ^^4      ^8X  g   Q h\        ^ ^ 4      \        8X  g   Q h\        ^^ 4      \        8X  g   Q h\        RR4      \        8X  g   Q h\        R^4      P                  R4      '       g   Q h\        RR4      P                  \        ^,          4      '       g   Q hR\        ^
^d4      ,          P                  R4      '       g   Q h\        \        \        4      ^ 8X  g   Q h\        \        \        ) \        4      4      '       g   Q h\        \        R\        4      4      '       g   Q h\        \        ^ \        4      4      '       g   Q h\        \        R4      \        R\        4      u;8X  d	   ^ 8X  g   Q h Q h\        \        R4      \        8X  g   Q h\        \        R4      \        ) 8X  g   Q h\        RR\        ^,          ,
          4      P                  R4      '       g   Q h\        RR4      ^ 8X  g   Q h\        R^4      P                  R4      '       g   Q h\        RR4      P                  R4      '       g   Q hR# )r$   g:0yE>r.   g      @z-255.5815971722918z-0.5119253100282322z-0.1952720278805729485g      @      ?g     @@r0   i"9	 grP@r/   z18.157330562703710339Nr   r&   r'   gUUUUUU?l     I5       r         ?       @y+5?\K9ʵ?gUUUUUUտ)r   r   r4   	gammaprodrH   absrfjr   ffbinomialbetar7   r:   )hepGs      r   test_gamma_quotientsr[      sh   BFA	BAbT2b'"a'''bVRD!S(((y"rd#a1ga1go56===y"R"Q(1RT71RT7?1RT7+B1QqS6+IIJTQQQa7a<<c!9b8r>>a7::eAaCjq)****"#89<<=UVVVVb8q==b8r>>a7b==a7a<<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""""43<##$9::::C=Asd3    baC=DSM.Q.....D>S   D>cT!!!R!V NOOOOS>Q1:==$%:;>>?VWWWWr
   c                     ^\         n        \        ^4      P                  \        ^,          ^,          4      '       g   Q h\        R4      P                  \        ^,          ^,          4      '       g   Q h\        \        ^4      4      P                  \        ^,          ^,          4      '       g   Q h\        ^d4      P                  ^4      '       g   Q h\        ^ 4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h\        R4      P                  \        ^4      ) ^,          4      '       g   Q h\        R4      ^ 8X  g   Q h\        R4      P                  \        ^4      ^x,          4      '       g   Q h\        R4      ^ 8X  g   Q h\        R4      ^ 8X  g   Q h\        \        \        4      4      '       g   Q h\        R4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h\        \        RR4      4      P                  ^ 4      '       g   Q h\        \        RR4      4      P                  ^ 4      '       g   Q h\        \        RR4      4      P                  ^ 4      '       g   Q h\        \        RR4      4      P                  R4      '       g   Q h\        \        RR4      4      P                  R4      '       g   Q h^2\         n        Rp \        \        RV 4      4      P                  ^ R	4      '       g   Q h^\         n        \        R4      ^
^",          ,          P                  R4      '       g   Q h\        ^ R4      ^ 8X  g   Q h\        ^ ^ 4      R8X  g   Q h\        ^ R^4      P                  R4      '       g   Q h\        RR
4      P                  R4      '       g   Q hR# )r$   g       @r/   KH9g"fD,@gOc5@g9@Ww'&l7z8236.5242296658162058024755079556629786895294952121891237gaU,$Db6y              ?Nr   gs~]r   rN   r&   rM   KHy      N      @y9A!@.@g9B.ֿrO   yN*i6?)	r   r   zetar   r7   r>   r   r:   r;   )ims    r   	test_zetarc      s   BF7::b!eai    9<<A	""""A<??2q519%%%%9<<???7::d9<<,----8;;Awrz""""8q==8;;s1vcz""""8q==:??c;>>$<??4    C/0144Q7777C/0144Q7777C/0144Q7777E% !$$T****F5!"%%d++++BF	CBC  E****BFL2r6!%%&LMMMM#;!!9#a=01111T?OPPPPr
   c                  N   ^\         n        \        R4      ^ 8X  g   Q h\        R4      ^ 8X  g   Q h\        R4      ^ 8X  g   Q h\        ^ 4      R8X  g   Q h\        R4      R8X  g   Q h\        R4      R8X  g   Q h\        R4      R8X  g   Q h\        R4      R8X  g   Q h\        ^4      P                  \	        ^4      4      '       g   Q h\        ^4      P                  \
        ^,          ^,          4      '       g   Q h\        ^
4      P                  ^I\
        ^
,          ,          R,          4      '       g   Q h\        ^24      ^8  g   Q h\        ^<RR7      ^8  g   Q h\        ^<RR7      ^8X  g   Q h\        RRR7      ^8  g   Q h\        RRR7      ^8X  g   Q h\        R4      \        ^4      8X  g   Q hRp \        V 4      P                  ^^^V ,
          ,          ,
          \        V 4      ,          4      '       g   Q hRp \        V 4      P                  ^^^V ,
          ,          ,
          \        V 4      ,          4      '       g   Q h\        R4      P                  R	4      '       g   Q h\        R
4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h\        \        RR4      4      P                  R4      '       g   Q h\        \        RR4      4      P                  R4      '       g   Q hR# )r$   r/   r-   g  z\Ai jh droundingur%   g@o.>Vg?g^s?r]   r^   NrN   rM   r_   r   r&   g      r'   i      @              @      @y            @      Yr`   )r   r   altzetar   logr7   ra   r>   )ss    r   test_altzetaro      s   BF2;!2;!4=A1:2;$2;&   2;$3<=(((1:==Q    1:==Qr""""2;>>"RV)G+,,,,2;??2$q(((2$)))53'!+++53'1,,,4=GAJ&&&A1:==!A!H*d1g-....A1:==!A!H*d1g-....6?677773<??.////5>S!!!!6?c""""3uU#$'',,,,3ve$%((----r
   c                  D   ^\         n        \        \        4      ^8X  g   Q h^2\         n        \        ^d4      P	                  R4      '       g   Q h\        ^(\
        ,          4      P	                  R4      '       g   Q hR\         n        \        R4      p ^\         n        \        V ^,
          4      R8X  g   Q h\        \
        R,          \        R7      ^8  g   Q h\        R\        R7      ^8  g   Q h\        \
        R,          4      ^8X  g   Q h\        R4      ^8X  g   Q hR	# )
r$   z61.0000000000000000000000000000007888609052210118073522z61.0000000000000000000000000000000000000148407238666182r%   i  z1.02363019598118e-9934r   rf   i  N)r   r   ra   rH   r   r7   r   r   )vs    r   test_zeta_hugerr   >  s    BF9>>BF9<<PQQQQ2;>>RSSSSBFUABFqs8////4(+a///x(1,,,4=A:??r
   c                      ^\         n        \        ) ^
^(,          ,          p ^\         n        \        \	        V 4      4      R8X  g   Q h^2\         n        \        \	        V 4      4      R8X  g   Q h^\         n        R# )   z=2.55880492708712e+1233536161668617575553892558646631323374078z`2.5588049270871154960875033337384432038436330847333e+1233536161668617575553892558646631323374078N)r   r   r7   r   ra   )as    r   test_zeta_negativerv   M  s[    BF	BF
ABFtAw<ZZZZBFtAw<}}}}BFr
   c            	         ^\         n        R p R pV ! ^4      \        ^ ^4      u;8X  d   \        ^4      8X  g   Q h Q hV ! \        4      P                  R4      '       g   Q hV ! \        ) 4      P                  R4      '       g   Q hV ! \        ) ^,           4      P                  R4      '       g   Q hV ! \        \        ,           4      P                  R4      '       g   Q hV ! \        ) \        ,
          4      P                  R4      '       g   Q h\        V ^4      P                  R4      '       g   Q hV ! R4      P                  R4      '       g   Q hV ! R4      P                  R4      '       g   Q hV ! R4      P                  R4      '       g   Q hV ! R4      P                  R4      '       g   Q hV ! \        4      \        8X  g   Q hV! \        4      ^ 8X  g   Q h\        ^\        4      ^ 8X  g   Q hV! \        4      P                  R	4      '       g   Q hV! \        ) 4      P                  R
4      '       g   Q hV! \        \        ,           4      P                  R4      '       g   Q hV! \        ) \        ,
          4      P                  R4      '       g   Q hR\        ^^^
\        ,          \        ,          ,           4      ,          P                  R4      '       g   Q hV ! ^^
\        ,          \        ,          ,           4      P                  R4      '       g   Q h\        V ! \        4      4      '       g   Q h\        V ! \        ) 4      4      '       g   Q hV ! R4      P                  R4      '       g   Q hV ! R4      P                  V ! ^4      4      '       g   Q hV ! R4      P                  R4      '       g   Q h\        \        ^\        ^ \        4      4      4      '       g   Q h\        \        ^\        ^ \        4      4      4      '       g   Q h\        \        ^\        ^ \        ) 4      4      4      '       g   Q h\        \        ^\        ^\        4      4      4      '       g   Q h\        \        ^\        ^\        4      4      4      '       g   Q h\        \        ^\        ^\        ) 4      4      4      '       g   Q h\        \        ^\        \        \        4      4      4      '       g   Q h\        \        ^\        \        \        4      4      4      '       g   Q h\        \        ^\        \        ) \        ) 4      4      4      '       g   Q h^\         n        \        R 4      P                  \        RR4      4      '       g   Q h^\         n        R# )!r$   c                     \        ^ V 4      # r   psizs   &r   <lambda> test_polygamma.<locals>.<lambda>X  
    S1Xr
   c                     \        ^V 4      # r   rz   r|   s   &r   r~   r   Y  r   r
   guTE?g,$7@gdGDE@g?Vcb?绽|=r^   ga#?g^jJ@g69u@z%0.46317279488182026118963809283042317z$2.4821070143037957102007677817351115Ny)\Z[?ÕޛD?ybabJr?N<g& _g\)c=y|=|=yF<	 _    _AyWw'&l7Ww'&l7y\)c-ȥ\)c-HyY/?Iyϐ=ӿ{?i@B yhVKٿy_%F@@K?rk   ri   y      Y      @y]vCq@˼@y            ?)r   r   r{   digammar7   r   rT   findrootrH   r:   r;   r>   )psi0psi1s     r   test_polygammar   V  s%   BFDD7c!Ah,'!*,,,,, 8;;*++++9<<*++++A;>>,----1:==CDDDDA;>>BCCCCD! 23333;>>,----;>>01111  !HIIII  !PQQQQ99>>q:??8;;*++++9<<*++++1:==CDDDDA;>>EFFFF#a"R%'	""&&'PQQQQ"R%'	?FGGGGcsd<??,----:==a!!!!=HIIIIQs1Sz"####Qs1Sz"####Qs1cT{#$$$$Qs1Sz"####Qs1Sz"####Qs1cT{#$$$$Qs3s|$%%%%Qs3s|$%%%%QsC4~&''''BF8$KMs tuuuuBFr
   c                      ^d\         n        \        \        ^ \        4      4      R8X  g   Q h\        \        ^
\        4      4      R8X  g   Q hR# )d   zf0.9772133079420067332920694864061823436408346099943256380095232865318105924777141317302075654362928734zf-12.98876181434889529310283769414222588307175962213707170773803550518307617769657562747174101900659238N)r   r   r   r{   r7   r	   r
   r   test_polygamma_high_precr     sS    BFs1Ry>  F  F  F  Fs2bz?  G  G  G  Gr
   c                     ^\         n        R p R pR pV ! R4      P                  \        ) ^\	        ^4      ,          ,
          4      '       g   Q hV ! ^4      P                  \        ) 4      '       g   Q hV! R4      P                  R\
        ^,          ,          4      '       g   Q hV! ^4      P                  \
        ^,          ^,          4      '       g   Q hV! R4      P                  \
        ^,          ^\        ,          ,           4      '       g   Q hV! ^4      P                  R\        ,          4      '       g   Q h^\         n        R\        ,          ^\        ^4      ,          \
        ^,          ,          ,           p^\         n        \        ^R	4      P                  V4      '       g   Q h\        ^R4      P                  \
        ^,          4      '       g   Q hR# )
r$   c                     \        ^ V 4      # ry   rz   r|   s   &r   r~   +test_polygamma_identities.<locals>.<lambda>  r   r
   c                     \        ^V 4      # r   rz   r|   s   &r   r~   r     r   r
   c                     \        ^V 4      # )r   rz   r|   s   &r   r~   r     r   r
   r/   r-   NrN   iJg?)
r   r   r   eulerrm   r7   catalanaperyr6   r{   )r   r   psi2rh   s       r   test_polygamma_identitiesr     sE   BFDDD9<<qQx((((7::uf9<<BE	""""7::b!eAg:==Q7*++++7::bhBFU
1T!W9RU?"ABFq;>>!q:==Qr
   c                     R  p ^2\         n        R\        ^^4      ,          pR\        ^^4      ,          pV ! RV,          4      ) V ! RV,          4      ,
          V ! R^V,           ,          4      ,           V ! R^V,
          ,          4      ,           p^\        ,          \	        R\        ^4      ,          \        ,          4      ,          pVP                  V4      '       g   Q h^\         n        R# )c                     \        ^ V 4      # ry   rz   r|   s   &r   r~   'test_foxtrot_identity.<locals>.<lambda>  r   r
   r/   Nr   r   )r   r   fractionr7   sechr6   r   )r   ru   bxys        r   test_foxtrot_identityr     s     DBF	hqmA	hqmA	c!etDF|#d3!9o5S!A#YGA	"T#d1g+b.!!A4477N7BFr
   c                  n   ^d\         n        \        \        ^2\        4      4      R8X  g   Q h\        \        ^2\        ^\
        ,          ,           4      4      R8X  g   Q h\        \        ^2\        ^\
        ,          \        ,          ,           4      4      R8X  g   Q h^\         n        \        \        ^2\        4      4      R8X  g   Q h\        \        ^2\        ^\
        ,          ,           4      4      R8X  g   Q h\        \        ^2\        ^\
        ,          \        ,          ,           4      4      R8X  g   Q hR# )r   zf-1344100348958402765749252447726432491812.641985273160531055707095989227897753035823152397679626136483zw-0.00000000000000000189793739550804321623512073101895801993019919886375952881053090844591920308111549337295143780341396z(-0.00000000000000005225169411521692489752254721556835657523758895106315132447859377385233700094871256507814151956624433 - 0.00000000000000001813157041407010184702414110218205348527862196327980417757665282244728963891298080199341480881811613j)z-1.34410034895841e+39z-1.89793739550804e-18z.(-5.2251694115217e-17 - 1.81315704140701e-17j)N)r   r   r   r{   r7   erT   r	   r
   r   test_polygamma_high_orderr     s    BFs2r{   H  H  H  Hs2rBqDy!"  '`  `  `  `s2rBqDF{#$ *^ _ _ _ BFs2r{6666s2rBqDy!"&====s2rBqDF{#$(XXXXr
   c                  p   ^\         n        \        ^ 4      ^ 8X  g   Q h\        ^4      ^8X  g   Q h\        ^4      R8X  g   Q h\        ^4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h\        ^
R,          4      P                  R4      '       g   Q h\        R4      P                  ^^\	        ^4      ,          ,
          4      '       g   Q h\        \
        4      \
        8X  g   Q h\        R	4      R
8X  g   Q h\        R4      P                  R4      '       g   Q hR# )r$   rL   g_b7@r   gS@r/   NgUUUUUU?r   y       @        y      ?        rP   yk4q?@w?)r   r   harmonicr   rm   rH   r	   r
   r   test_harmonicr     s    BFA;!A;!A;#A;>>*++++F45555BH  !23333C=AaAhJ''''C=CD>V###D>DEEEEr
   c                      R \         n        \        R4      ^,          p ^\         n        \        \	        V 4      4      R8X  g   Q h^2\         n        \        \	        V 4      4      R8X  g   Q h^\         n        R# )  z6.26075321389519e+12458010678z@6.2607532138951929201303779291707455874010420783933e+12458010678Nr   r   r   r   r   r4   r   s    r   test_gamma_huge_1r     sY    BFFaABFuQx=;;;;BFuQx=^^^^BFr
   c                      R \         n        \        ^
^d,          4      ^,          p ^\         n        \        \	        V 4      4      R8X  g   Q h^2\         n        \        \	        V 4      4      R8X  g   Q hR# )r   zw1.82341134776679e+51729974693233641689901335581750771368291828240422018860515119656908623426021308685461258226190190661z1.82341134776678875374414910350027596939980412984e+51729974693233641689901335581750771368291828240422018860515119656908623426021308685461258226190190661Nr   r   s    r   test_gamma_huge_2r     sm    BFBGrABFuQx=	34 4 4 BFuQx=	TU U Ur
   c                     R \         n        ^
^P,          ^,          ^
^F,          \        ,          ^,          ,           p ^\         n        \        V 4      p\	        VP
                  4      R8X  g   Q h\	        VP                  4      R8X  g   Q h^2\         n        \        V 4      p\	        VP
                  4      R8X  g   Q h\	        VP                  4      R8X  g   Q hR# )r   ze-6.82925203918106e+2636286142112569524501781477865238132302397236429627932441916056964386399485392600zd8.54647143678418e+2636286142112569524501781477865238132302397236429627932441916056964386399485392600z-6.8292520391810548460682736226799637356016538421817e+2636286142112569524501781477865238132302397236429627932441916056964386399485392600z8.5464714367841748507479306948130687511711420234015e+2636286142112569524501781477865238132302397236429627932441916056964386399485392600Nr   r   rT   r4   r   realimagr   r   s     r   test_gamma_huge_3r     s    BF
B!b"fQhl"ABFaAqvv;	 ! ! ! qvv;	      BFaAqvv;	CD D D qvv;	BC C Cr
   c                      Rp ^\         n        \        \        V 4      4      R8X  g   Q h^2\         n        \        \        V 4      4      R8X  g   Q h^\         n        R# )i  z2(8.95783268539713e+5164 - 1.94678798329735e+5164j)zx(8.9578326853971339570292952697675570822206567327092e+5164 - 1.9467879832973509568895402139429643650329524144794e+5164j)Ny      @     v@)r   r   r   r4   r   s    r   test_gamma_huge_4r     sZ    ABFuQx=<= = =BFuQx=	   BFr
   c                     R \         n        ^
^<,          \        ,          ^,          p ^\         n        \        V 4      p\	        VP
                  4      R8X  g   Q h\	        VP                  4      R8X  g   Q h^2\         n        \        V 4      p\	        VP
                  4      R8X  g   Q h\	        VP                  4      R8X  g   Q h^\         n        R# )r   zO-3.27753899634941e-227396058973640224580963937571892628368354580620654233316839zN-7.1519888950416e-227396058973640224580963937571892628368354580620654233316841zr-3.2775389963494132168950056995974690946983219123935e-227396058973640224580963937571892628368354580620654233316839zr-7.1519888950415979749736749222530209713136588885897e-227396058973640224580963937571892628368354580620654233316841Nr   r   s     r   test_gamma_huge_5r     s    BF
B
QABFaAqvv;kkkkqvv;jjjjBFaAqvv;	-. . . qvv;	-. . . BFr
   c                     ^d\         n        ^\        \        ^
4      R,          ,          ,           p ^\         n        \	        V 4      p\        VP                  4      R8X  g   Q h\        VP                  4      R8X  g   Q h^2\         n        \	        V 4      p\        VP                  4      R8X  g   Q h\        VP                  4      R8X  g   Q hR# )r   r   z2.0z1.84556867019693e-1000z91.8455686701969342787869758198351951379156813281202e-1000N)r   r   rT   r   r4   r   r   r   )ru   r   s     r   test_gamma_huge_7r     s    BF	Ac"gtmOABFaAqvv;% qvv;2222BFaAqvv;%qvv;VVVVr
   c                     ^\         n        \        ^ 4      P                  \        54      '       g   Q h^\         n        \        ^4      P                  R4      '       g   Q h\        ^4      P                  R4      '       g   Q h\        ^4      P                  R4      '       g   Q h\        ^4      P                  R4      '       g   Q h^\         n        \        ^4      P                  R4      '       g   Q h\        ^4      P                  R	4      '       g   Q h\        ^4      P                  R4      '       g   Q h\        ^4      P                  R4      '       g   Q hR# )
r$   z -0.07281584548367672486058637587z!-0.009690363192872318484530386035z 0.002053834420303345866160046543z 0.002325370065467300057468170178gT3`?g;c?Ng*gF|؃)r   r   	stieltjesr   r   r	   r
   r   test_stieltjesr     s   BFQ<??E6""""BFQ<??=>>>>Q<??>????Q<??=>>>>Q<??=>>>>BFQ<??;<<<<Q<??<====Q<??;<<<<Q<??<====r
   c                     ^\         n        \        ^ 4      \        R4      u;8X  d	   ^ 8X  g   Q h Q h\        ^4       U u. uF  p \	        V 4      NK  	  up . RO8X  g   Q h\        \	        R4      4      R8X  g   Q h\        \        \        4      4      '       g   Q h\        \	        \        4      4      '       g   Q h\        \        \        4      4      '       g   Q h\        \        4      \        8X  g   Q h\	        \        4      \        8X  g   Q h\        \        4      \        8X  g   Q h\        \	        \        ) 4      4      '       g   Q h\        R4      P                  R4      '       g   Q h\        R4      P                  R	4      '       g   Q h\        ^4       Uu. uF  p\        V4      NK  	  up. R
O8X  g   Q h\        ^ RR4       Uu. uF  p\        V4      NK  	  up. RO8X  g   Q h\        R4      pV^ 8X  d   \        V\        4      '       g   Q h\        R4      pVR8X  d   \        V\        4      '       g   Q hR# u up i u upi u upi )r$   r   z3.24570818422368e+1177245gffffff?gd?Nr   )r   r   r      i   i   i {l     2ft y       @      @ygƿ1?)r   r      l    l  i \&l     rM   )r   r   r   rM   r   r   i y              )r   r   barnesgr5   superfacr   r:   r;   hyperfacrH   r   
isinstancer>   )r@   rI   ru   s      r   test_barnesgr   "  s   BF1:))))))!&q*AHQK*.____x~"====##3<3C=CC=C3$    3<??-....4=EFFFF!&q*AHQK*.\\\\!&qB0AHQK04SSSSA6jC((((A7z!S)))))! + +0s   H?I	I	c                  ^	   ^\         n        R U u. uF  p \        V 4      NK  	  pp V F7  p \        ^V 4      P	                  \        ^V ,
          4      ) 4      '       d   K7  Q h	  V F4  p \        ^ V 4      P	                  V ^V ,
          ,          4      '       d   K4  Q h	  V F;  p \        RV 4      P	                  V ^V ,
          ^,          ,          4      '       d   K;  Q h	  V FI  p \        RV 4      P	                  V ^V ,           ,          ^V ,
          ^,          ,          4      '       d   KI  Q h	  V F^  p \        RV 4      P	                  V ^^V ,          ,           V ^,          ,           ,          ^V ,
          ^,          ,          4      '       d   K^  Q h	  \        ^^4      P	                  R4      '       g   Q h\        ^R4      P	                  R4      '       g   Q h\        ^R4      P	                  R4      '       g   Q h\        ^R4      P	                  R4      '       g   Q h\        ^R4      P	                  R4      '       g   Q h\        ^R4      P	                  R4      '       g   Q h\        ^R4      P	                  R4      '       g   Q h\        ^R\        \        4      ,          4      P	                  R4      '       g   Q h\        RR4      P	                  R4      '       g   Q h\        RR4      P	                  R4      '       g   Q h\        ^R4      P	                  R4      '       g   Q h\        RR4      P	                  R4      '       g   Q h\        RR4      P	                  R4      '       g   Q h\        R\        ^,          ,           R\        ^,          ,           4      P	                  R4      '       g   Q h\        R^4      P	                  \        R4      4      '       g   Q h\        RR4      P	                  \        R4      ) 4      '       g   Q h\        R	R 4      P	                  R!4      '       g   Q h\        R	^4      P	                  R"4      '       g   Q h\        R#R4      P	                  \        ^4      R
,          4      '       g   Q h\        R#R4      P	                  R$4      '       g   Q hR# u up i )%r$   r/   ?gV?y        ?g?i  g=.g?i2  rL   g      +@g      ?NrM   rj   )r   r/   gGz?r   r   rM   y              ?rj   r   rN   r&   y[-F@U޴r   g!fGggy#]ſ`j?y8n+Yd?^n)ӿggᾁyp/M?[]?g jߊg-C6zyB8$P?g ?gtHgWcFyQ K1Q y               iL)r   r   	mpmathifypolylogr   rm   r6   rT   ra   rl   r   )r}   zss     r   test_polylogr   7  s   BF E	F E1)A, EB	F1((#ac(3333R1((AaC1111RA))!QqS1H*5555RA))!QqS'1Q3(*:;;;;RA))!QqsU1a4Z.!A#*ABBBBR1a=HIIII1b>344441c?344441d677771dNOOOO1c?JKKKK1d677771c$q'k"%%&TUUUU2st$$$$2t	****1c?344442su%%%%2r?g&&&&3qs7C!G$''(XYYYY4?d4j))))4~....3+,//0BCCCC3 MNNNN5"  Q----5$""4((((5 
Gs   R*c                     ^\         n        \        ^ R4      R,          P                  R4      '       g   Q h\        ^R4      R,          P                  R4      '       g   Q h\        ^R4      P                  R4      '       g   Q h\        RR4      P                  R4      '       g   Q h\	        ^ ^ 4      ^8X  g   Q h\	        ^ ^4      ^8X  g   Q h\	        ^ ^4      ^8X  g   Q h\	        ^ \
        4      ^8X  g   Q h\	        ^ \
        ) 4      ^8X  g   Q h\        \	        ^ \        4      4      '       g   Q h\	        ^^4      ^8X  g   Q h\	        ^^ 4      ^ 8X  g   Q h\	        ^\
        4      \
        8X  g   Q h\	        ^\
        ) 4      \
        ) 8X  g   Q h\        \	        ^\        4      4      '       g   Q h\	        ^R4      ^ 8X  g   Q h\	        ^^ 4      ^ 8X  g   Q h\	        ^
4      R8X  g   Q h\	        ^
^4      R8X  g   Q h\	        ^
R4      R8X  g   Q h\	        ^R4      R8X  g   Q h\	        ^2R4      P                  R	4      '       g   Q h^P\         n        \	        ^2R4      R
8X  g   Q h\	        ^(^24      R8X  g   Q h\	        ^J4      R8X  g   Q h^\         n        \	        ^
R4      R8X  g   Q h\	        R^ 4      P                  \        \        R,          4      4      '       g   Q hR# )r$   r   go   ?g8   ?y              @g      @g*
0(@i i z3.4746902914629720259e74l   vreh@3 VBFv'r.y	UR|s&T l   .e0u&Y:j+I<1 |6vJb l   rw_3dtL#m,`3	.v(*A=V y              4@r/   Nr   y)fɳ?@r   iii0y @j#MB (3VjwB)
r   r   polyexpr   bellrH   r:   r;   sincr7   r	   r
   r   test_bell_polyexpr   U  sz   BFAeV#''6666AeV#''55551R=IJJJJ2c?23333!9>>!9>>!9>>#;!3$<1a!9>>!9>>#;#3$<C4a":??!9>>8v1:B<8####;*$$$3<??56666BF3<ffff2;cccc8ffffBF3<8888A;>>$r#v,''''r
   c                  <   ^\         n        \        R4      P                  R4      '       g   Q h\        ^4      P                  R4      '       g   Q h\        ^4      \        8X  g   Q h\        \        4      ^ 8X  g   Q h\        \        \        4      4      '       g   Q hR# )r$   r   gMҵ?Ny2k?-DT!	@)r   r   	primezetar   rH   r:   r;   r	   r
   r   test_primezetar   {  sx    BFS>KLLLLQ<??23333Q<3S>Q3    r
   c                     ^\         n        \        R4      P                  R4      '       g   Q h\        R4      P                  R4      '       g   Q h\        R^R7      P                  R	4      '       g   Q h\        R
^R7      P                  R4      '       g   Q h\        R
^R7      P                  R4      '       g   Q h\        R^R7      P                  R4      '       g   Q hRR lp V ! \        P                  R4      R4      '       g   Q hV ! \        P                  R4      R4      '       g   Q hV ! \        P                  R4      R4      '       g   Q hV ! \        P                  R^R7      R4      '       g   Q hV ! \        P                  R
4      R4      '       g   Q hV ! \        P                  R
^R7      R4      '       g   Q hV ! \        P                  R^R7      R4      '       g   Q hV ! \        P                  R^R7      R4      '       g   Q hV ! \        P                  R^R7      R4      '       g   Q hV ! \        P                  R^R7      R4      '       g   Q hV ! \        P                  R^R7      R4      '       g   Q hV ! \        P                  R^R7      R4      '       g   Q hV ! \        P                  R^R7      R4      '       g   Q hV ! \        P                  R^R7      R4      '       g   Q hR# )r$   
derivativec                 J    \        W,
          4      V\        V4      ,          8  # )N)rR   )r   r   tols   &&&r   r   test_rs_zeta.<locals>.ae  s    13x#c!f*$$r
   Ny      ?     j@y& 6#+?@y      ?     j@y5g?7-m^I?y      ?    .Ay@V%<y      ?    .AyˇD@-2yV!V1y      ?    .y.F'@	p 1@)gư>y      ?     jy& 6#+?y      ?     jy5g?7-m^Iyz?ی^p@y.F'@	p 1y$P??y,ς?%@lUS>yԂ]5MԺ$k@y!&@\_Uy(1Sg@y-6ڿ+~y"虺m"'=UA@yِ@[@MomykLVy"DЫ@)r   r   ra   r   fp)r   s    r   test_rs_zetar     s   BF NOOOO  !NOOOO+../]^^^^
q),,-[\\\\
q),,-^____+../]^^^^%bggk"$RSSSSbggl#%RSSSSbggh!QRRRRbggh1g-/]^^^^bggfoOPPPPbggfg+-^____bggkag02_````bggkag02abbbbbggkag02`aaaabggkag02abbbbbgglqg13cddddbgglqg13`aaaabgglqg13abbbbbgglqg13abbbbr
   c                  N   ^\         n        \        R4      P                  R4      '       g   Q h\        R^R7      P                  R4      '       g   Q h\        R^R7      P                  R4      '       g   Q h\        R\        ,           ^R7      P                  R4      '       g   Q hR# )r$   r   g1/[㳄@r   NgKKgFhnqyQ1@a߻w)r   r   siegelzr   rT   r	   r
   r   test_siegelzr     s    BF6?455556a(++,CDDDD6a(++,CDDDD6!8q),,-[\\\\r
   c                  .   ^\         n        \        ^RRR7      p \        ^RRR7      p\        ^RRR7      p\        V 4      P	                  R4      '       g   Q h\        V4      P	                  R4      '       g   Q h\        V4      pVP
                  P	                  R4      '       g   Q hVP                  P	                  R4      '       g   Q h^\         n        \        ^RRR7      p \        ^R	RR7      p\        ^R
RR7      p\        V 4      P	                  R4      '       g   Q h\        V4      P	                  R4      '       g   Q h\        V4      pVP
                  P	                  R4      '       g   Q hVP                  P	                  R4      '       g   Q hR# )r$   z1e-10T)exactz-1e-10z1e-10jg# _Bgox?z1e-50z-1e-50z1e-50j1e50z-1e50z40.57721566490153286060651209008240243104215933593992Nga_g    _)r   r   faddra   r   r   r   )s1s2s3r}   s       r   test_zeta_near_1r     sC   BF	a	%B	a	&B	a	&B8;;788888;;67777RA6699+,,,,669934444BF	a	%B	a	&B	a	&B8;;v8;;wRA6699KLLLL6699Wr
   N)#mpmathmpmath.libmpr   r   r   r   r"   r+   rA   rJ   r[   rc   ro   rr   rv   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r	   r
   r   <module>r      s     ; ;3+Ze>N%pN/Xb!QF.<+ZG
 "
YF
UC(
"W>**)<$(L!c6]r
   