ó
    ˆ*£hh´  ã                   ó²   • S r SSKr " S S\5      rSSKJr  S rS rS	 rS
 r	SS jr
S/4S jrSS jrSS jrSS jrSS jr\SS j5       r\SS j5       rg)aI  
---------------------------------------------------------------------
.. sectionauthor:: Juan Arias de Reyna <arias@us.es>

This module implements zeta-related functions using the Riemann-Siegel
expansion: zeta_offline(s,k=0)

* coef(J, eps): Need in the computation of Rzeta(s,k)

* Rzeta_simul(s, der=0) computes Rzeta^(k)(s) and Rzeta^(k)(1-s) simultaneously
  for  0 <= k <= der. Used by zeta_offline and z_offline

* Rzeta_set(s, derivatives) computes Rzeta^(k)(s) for given derivatives, used by
  z_half(t,k) and zeta_half

* z_offline(w,k): Z(w) and its derivatives of order k <= 4
* z_half(t,k): Z(t) (Riemann Siegel function) and its derivatives of order k <= 4
* zeta_offline(s): zeta(s) and its derivatives of order k<= 4
* zeta_half(1/2+it,k):  zeta(s)  and its derivatives of order k<= 4

* rs_zeta(s,k=0) Computes zeta^(k)(s)   Unifies zeta_half and zeta_offline
* rs_z(w,k=0)    Computes Z^(k)(w)      Unifies z_offline and z_half
----------------------------------------------------------------------

This program uses Riemann-Siegel expansion even to compute
zeta(s) on points s = sigma + i t  with sigma arbitrary not
necessarily equal to 1/2.

It is founded on a new deduction of the formula, with rigorous
and sharp bounds for the  terms and rest of this expansion.

More information on the papers:

 J. Arias de Reyna, High Precision Computation of Riemann's
 Zeta Function by the Riemann-Siegel Formula I, II

 We refer to them as I, II.

 In them we shall find detailed explanation of all the
 procedure.

The program uses Riemann-Siegel expansion.
This  is useful when t is big, ( say  t > 10000 ).
The precision is limited, roughly it can compute zeta(sigma+it)
with an error less than exp(-c t) for some constant c depending
on sigma.  The program gives an error when the Riemann-Siegel
formula can not compute to the wanted precision.

é    Nc                   ó   • \ rS rSrS rSrg)ÚRSCacheé5   c                 ó   • SS0 0 /U l         g )Nr   é
   )Ú	_rs_cache)Úctxs    ÚT/home/mande/repo/quber/.venv/lib/python3.13/site-packages/mpmath/functions/rszeta.pyÚ__init__ÚRSCache.__init__6   s   € Ø˜B  B˜ˆ�ó    © N)Ú__name__Ú
__module__Ú__qualname__Ú__firstlineno__r   Ú__static_attributes__r   r   r
   r   r   5   s   † õ(r   r   é   )Údefunc                 óf  • US-   nUS-  n[        U R                  SUS-   -  5      SU-  S-   U R                  U5      -
  5      nU R                  SU-  5      n[        U R                  SU-  5      U R                  U5      S-   U-   5      nXpl        0 nU R                  US'   U R
                  US	'   [        SSU-  S	-   5       H  n	X‰S	-
     U R
                  -  X‰'   M     US-   U l        0 n
0 n[        SUS	-   5       HO  n	S
U	-  U R                  SU	-  5      -  nU R                  U5      U R                  SU	-  5      -  nXÈSU	-     -  X©'   MQ     [        SSU-  S	-   5       H2  n	U R                  U R                  U	5      -  nUSU	-  -  nXØU	   -  X¹'   M4     SU l        SU R                  U5      -
  n0 n[        SUS	-   5       H‡  n	SU l        [        U R                  SU	S-   -  5      SU	-  U-   5      nUU l        Sn[        SU	S	-   5       H   nUS
U-  U
U   -  USU	-  SU-  -
     -  -  nM"     S
U	S	-   -  U R                  -  U-  Xù'   M‰     0 n[        SUS	-   5       Hx  n	SU l        [        U R                  SU	S-   -  5      SU	-  U-   5      nUU l        Sn[        SU	S	-   5       H&  nUU R                  U	U-
  -  U
U   -  X¹U-
     -  -  nM(     UUU	'   Mz     SU l        SU R                  U5      -
  n[        SSU-  U-   5      nUU l        U R                  U R                  S5      5      S-  nU R                  S5      S-  n0 n[        SU5       H9  n	SU l        [        SSU	-  U-   5      nUU l        UXù   -  UUU	   -  -   USU	-  '   M;     [        S	SU-  S5       H  n	SUU	'   M
     X4UU/$ )a   
Computes the coefficients  `c_n`  for `0\le n\le 2J` with error less than eps

**Definition**

The coefficients c_n are defined by

.. math ::

    \begin{equation}
    F(z)=\frac{e^{\pi i
    \bigl(\frac{z^2}{2}+\frac38\bigr)}-i\sqrt{2}\cos\frac{\pi}{2}z}{2\cos\pi
    z}=\sum_{n=0}^\infty c_{2n} z^{2n}
    \end{equation}

they are computed applying the relation

.. math ::

    \begin{multline}
    c_{2n}=-\frac{i}{\sqrt{2}}\Bigl(\frac{\pi}{2}\Bigr)^{2n}
    \sum_{k=0}^n\frac{(-1)^k}{(2k)!}
    2^{2n-2k}\frac{(-1)^{n-k}E_{2n-2k}}{(2n-2k)!}+\\
    +e^{3\pi i/8}\sum_{j=0}^n(-1)^j\frac{
    E_{2j}}{(2j)!}\frac{i^{n-j}\pi^{n+j}}{(n-j)!2^{n-j+1}}.
    \end{multline}
é   ç       @r   é   é   é   é(   r   r   éÿÿÿÿé   é	   é   Ú2g      Ø?)ÚmaxÚmagÚ	_eulernumÚprecÚoneÚpiÚrangeÚmpfÚfacÚjÚsqrtÚexpjpi)r	   ÚJÚepsÚnewJÚneweps6ÚwpvwÚEÚwppiÚpipowerÚnÚvÚwÚvaÚwaÚwpp1aÚP1Úwpp1ÚsumpÚkÚP2Úwpp2Úwpc0ÚwpcÚmuÚnuÚcs                             r
   Ú_coefrG   N   sì  € ð: ˆQ‰3€DØ�"‰f€Gô ˆs�w‰w�r˜4 ™6‘{Ó# Q t¡V¨A¡X¨c¯g©g°gÓ.>Ñ%>Ó?€Dð 	�‰�a˜‘fÓ€Aô" ˆs�w‰w�r˜$‘wÓ §¡¨£¨q¡°$Ñ!6Ó7€DØ„HØ€GØ—‘€GˆA�JØ—‘€GˆA�JÜ�1�Q�t‘V˜A‘XÖˆØ˜q™S‘\ #§&¡&Ñ(ˆ‹
ñ ð
 �A‰v€C„HØ€AØ€AÜ�1�T˜!‘VŽ_ˆØ�1‰W�s—}‘} Q q¡SÓ)Ñ)ˆØ�W‰W�R‹[˜Ÿ™  1¡›Ñ%ˆØ˜˜!™‘‰_ˆ‹ñ ô �1�Q�t‘V˜A‘XÖˆØ�W‰W�S—W‘W˜Q“ZÑˆØˆq�!‰t‰9ˆØ˜‘
‰]ˆ‹ñ ð €C„HØ�—‘˜Ó Ñ €EØ	€BÜ�1�T˜!‘VŽ_ˆØˆŒÜ�3—7‘7˜2˜q ™s™8Ó$ Q q¡S¨¡YÓ/ˆØˆŒØˆÜ�q˜˜1™–ˆAØ�b˜1‘W  !¡Ñ$ Q q¨¡s¨1¨Q©3¡w¡ZÑ/Ñ/ŠDñ à�a˜‘c‰{˜CŸE™EÑ! $Ñ&ˆ‹ñ ð 
€BÜ�1�T˜!‘VŽ_ˆØˆŒÜ�3—7‘7˜2˜q ™s™8Ó$ Q q¡S¨¡YÓ/ˆØˆŒØˆÜ�q˜˜1™–ˆAØ�S—U‘U˜Q˜q™S‘\ Q q¡TÑ)¨!¨a©C©&Ñ0Ñ0ŠDñ àˆˆ1‹ñ ð €C„HØˆs�w‰w�wÓÑ€DÜ
ˆa��$‘�t‘Ó
€CØ€C„HØ	�‰�#—'‘'˜#“,Ó	 Ñ	!€BØ	�‰�DÓ	˜!Ñ	€BØ€AÜ�1�TŽ]ˆØˆŒÜ�!�A�a‘C˜‘H‹oˆØˆŒØ�B‘E‘˜"˜R ™U™(Ñ"ˆˆ!ˆA‰#‹ñ	 ô
 �1�Q�t‘V˜AÖˆØˆˆ!‹ñ à˜1˜gÐ&Ð&r   c                 ó  ^ • T R                   nXS   ::  a  X#S   :¼  a
  US   US   4$ T R                  R                  n [        T R                  X5      nUT R                  l        T T R                  LaR  [	        U 4S jUS   R                  5        5       5      US'   [	        U 4S jUS   R                  5        5       5      US'   UT R                   S S & T R                   S   T R                   S   4$ ! UT R                  l        f = f)Nr   r   r   r   c              3   óN   >#   • U  H  u  pUTR                  U5      4v •  M     g 7f©N©Úconvert©Ú.0r?   r7   r	   s      €r
   Ú	<genexpr>Úcoef.<locals>.<genexpr>Î   ó    øé € ÐG²©e¨q˜˜#Ÿ+™+ a›.Õ)²ùó   ƒ"%c              3   óN   >#   • U  H  u  pUTR                  U5      4v •  M     g 7frJ   rK   rM   s      €r
   rO   rP   Ï   rQ   rR   )r   Ú_mpr%   rG   ÚdictÚitems)r	   r.   r/   Ú_cacheÚorigÚdatas   `     r
   ÚcoefrZ   Ä   sà   ø€ Ø�]‰]€FØ�1‰Iƒ~˜#¨¡Ó*Ø�a‰y˜& ™)Ð#Ð#Ø�7‰7�<‰<€DðÜ�S—W‘W˜aÓ%ˆàˆ�‰ŒØ
�#—'‘'ÒÜÔG°t¸A±w·}±}´ÓGÓGˆˆQ‰ÜÔG°t¸A±w·}±}´ÓGÓGˆˆQ‰Ø€C‡M�M‘!ÐØ�=‰=˜Ñ˜SŸ]™]¨1Ñ-Ð-Ð-øð ˆ�‰�ús   ¿C4 Ã4Dc                 óˆ  • SU-  U-  nU R                  U5      nSU R                  U R                  5      -  U R                  U5      -   U R                  U5      -   S-  U R                  SU R                  -  5      S-  -
  nSU-  S-
  nU R                  US-   5      U R                  US-  S-   5      -
  S-  U R                  US-   S-  5      -
  n	XhU-  U	-   :  ae  US:”  a_  US-
  nU R                  US-   5      U R                  US-  S-   5      -
  S-  U R                  US-   S-  5      -
  n	XhU-  U	-   :  a  US:”  a  M_  Un
U
$ )Ngb—òk5„_@r   r   r   ç      @r   )Úlnr'   Úloggamma)r	   ÚxAÚxeps4ÚaÚxB1ÚxLÚaux1Úaux2ÚmÚaux3ÚxMs              r
   Úaux_M_Fpri   ß   sA  € ð �r‰>˜%Ñ€DØ�6‰6�$‹<€DØˆc�f‰f�S—V‘V‹nÑ˜SŸV™V C›[Ñ(¨¯©°«Ñ2°AÑ5°s·v±v¸aÀÇÁ¹hÓ7GÈÑ7IÑI€DØ	ˆ"‰ˆQ‰€AØ�<‰<˜˜!™Ó˜SŸ\™\¨!¨C©%°©'Ó2Ñ2°AÑ
5°s·|±|ÀQÀqÁSÈ"ÁHÓ7MÑ
M€DØ�D‘&˜$‘,Ó A a£CØ�‰EˆØ—‘˜Q˜q™SÓ! #§,¡,¨q°©u°Q©wÓ"7Ñ7¸Ñ:¸C¿L¹LÈ!ÈAÉ#ÈrÉÓ<RÑRˆð �D‘&˜$‘,Ó A a¥Cð 
€BØ€Ir   c                 ó~   • USU-  -  nXC-  S-  nX`R                  XuS-  -  US-
  S-  5      -  U-  n[        Xg5      nU$ )Néx  gJ¤R”_@r   r   r   )ÚpowerÚmin)	r	   r_   r`   ra   rb   rh   Úh1Úh2Úh3s	            r
   Úaux_J_neededrq   í   sR   € ð
 
��B‘‰€BØ	‰Ð&Ñ	&€BØ	�i‰i˜ ™E™ B q¡D¨!¡8Ó,Ñ	,¨rÑ	1€BÜ	ˆR‹€BØ€Ir   c                 óÜ7  ^ • T R                   nT R                  U5      nT R                  U5      nSU-
  nST l         T R                  UST R                  -  -  5      nXu-  nXv-  n	T R                  ST R                  U5      S-
  5      n
T R                  ST R                  U	5      S-
  5      nT R                  SU* 5      nUS-  nXÊ-  S-  nXË-  S-  nST l         US:”  a6  Sn[        R                  " SU5      S	-  n[        R                  " SU5      nSnO7S
n[        R                  " SU* 5      S	-  n[        R                  " SU* 5      nSnUS:”  a6  Sn[        R                  " SU5      S	-  n[        R                  " SU5      nSnO7S
n[        R                  " SU* 5      S	-  n[        R                  " SU* 5      nSnST l         SnSU-  T R                  US-  5      -  T R                  UU-  U* 5      -  U:¼  a<  US-   nSU-  T R                  US-  5      -  T R                  UU-  U* 5      -  U:¼  a  M<  [        SU5      nSnSU-  T R                  US-  5      -  T R                  UU-  U* 5      -  U:¼  a<  US-   nSU-  T R                  US-  5      -  T R                  UU-  U* 5      -  U:¼  a  M<  [        SU5      nSU-  SU-  U-  S-  :¼  dT  SU-  S-   U-   S:  dE  [        U5      US-  :”  d3  SU-  SU-  U-  S-  :¼  d!  SU-  S-   U-   S:  d  [        U5      US-  :”  a  UT l         [        S5      e[        UU5      nUSU-  -  nUSU-  -  nUSU-  -  nUSU-  -  n[        T UUUUU5      n[        T UUUUU5      n [        UU 5      n![        T UUUUU5      n"[        T UUUUU 5      n#[        U"U#5      n"Sn$ST R                  -  U$-  T R                  U$S-   5      -  n%U%U":”  a"  U$S-   n$ST R                  -  U%-  U$-  n%U%U":”  a  M"  0 n&[        R                  [        R                  -  U-  U-  n'[        R                  [        R                  -  U-  U-  n([!        SS5       H¦  n)[        R                  " U'U)S-  5      SU-  -  n*[        R                  " U(U)S-  5      SU-  -  n+[        U*U+5      n,T R                  U)S-   5      T R                  U)S-  S-   5      -  n-[        R                  " U-5      n-U,U--  [        UU5      -  U&U)'   M¨     [        SU-  S-
  S5      S-   n.SU$-  n/T R                  SU$-  5      n0[!        SU.5       H9  n)[        U0T R                  U/T R                  U)S-   5      -  U&U)   -  5      5      n0M;     U0T R                  U5      -   S-   T l         T R                  UST R                  -  -  5      nT R#                  U5      n1SSUU1-
  -  -
  n2T R#                  U2T R%                  S5      U0-  -  5      n3U2T R%                  S5      U0-  -  U3-
  n4U4S:  a  U3n3OU3S-   n3T R'                  U3T R%                  S5      U0* -  -  5      n2T R                  T R'                  ST R                  -  5      U$5      T R                  U$S-   5      S-  U$-  -  n50 n60 n7[)        T U$U55      u  n7n8U7R+                  5       n60 n9[!        U!SU-  S-
  5       H  n:SU9U:'   M
     0 n9U0T l         [!        SU!S-   5       H_  n)Sn;[!        SU$-  U)-
  S-
  SS5       H  n<U;U2-  U6U<   -   n;M     U;U9U)'   [!        SSU$-  U)-
  S-
  5       H  n<U<S-   U6U<S-      -  U6U<'   M     Ma     0 n=[        ST R                  SU-  U-  5      5      n>ST R                  S[        U5      -   U-  5      -   T R                  U5      -
  S-
  n?T R-                  ST R                  T R                  -  U-  U-  U-  U-  -  5      S-  n@[!        SU5       H\  n:T R                  T R                  T R                  U:S-
  5      5      5      T R#                  U:W@-  5      -
  U?-   nA[        UAU>5      U=U:'   M^     U=S   S-   T l         SSU-  -
  nB0 nCSUCS'   SUCS '   SUCS!'   SUCS"'   SUCS#'   SUCS$'   SUCS%'   SUCS&'   [!        SU5       GH9  n:U=U:   S-   T l         [!        SSU:-  S-  S-   5       Hî  n<SU:-  SU<-  -
  n)U)S:w  a\  T R.                  U)-  nDUDS-  nEWBUD-  S-  nFU)S-   * nGUGWCSU:S-
  U<S-
  4   -  UEUCSU:S-
  U<4   -  -   UFUCSU:S-
  U<S-
  4   -  -   UCSU:U<4'   Mp  SWCSU:U<4'   [!        SU<5       Hf  nHWCSU:UH4   T R%                  S'5      T R1                  SU<-  SUH-  -
  5      -  T R1                  U<UH-
  5      -  -  nIUCSU:U<4==   SU<UH-
  -  UI-  -  ss'   Mh     Mð     SWCSU:S(4'   SUCSU:S4'   SUCSU:SU:-  S-  S-   4'   GM<     [!        S(US-   5       Hd  nJ[!        S(U5       HQ  n:[!        S)[        SSU:-  S-  S-   5      5       H+  n<WJS:  d  U:S:  d  U<S:  d  U<SU:-  S-  :”  d  M#  SWCWJU:U<4'   M-     MS     Mf     [!        SUS-   5       H‘  nJ[!        SU5       H~  n:U=U:   S-   T l         [!        SSU:-  S-  S-   5       HU  n<SWJ-  S-
  WCUJS-
  U:S-
  U<S-
  4   -  SUU:-   S-
  -  UCUJS-
  U:S-
  U<S-
  4   -  -   n/U/UCUJS-
  U:S-
  U<S-
  4   -
  UCUJU:U<4'   MW     M€     M“     0 nK[        ST R                  SU-  U-  5      5      n>ST R                  S[        U5      -   U-  5      -   T R                  U5      -
  S-
  nLT R-                  ST R                  T R                  -  U-  U-  U-  U-  -  5      S-  nM[!        SU5       H\  n:T R                  T R                  T R                  U:S-
  5      5      5      T R#                  U:WM-  5      -
  WL-   nN[        UNU>5      WKU:'   M^     WKS   S-   T l         SSU-  -
  nO0 nPSUPS'   SUPS '   SUPS!'   SUPS"'   SUPS#'   SUPS$'   SUPS%'   SUPS&'   [!        SU5       GH9  n:WKU:   S-   T l         [!        SSU:-  S-  S-   5       Hî  n<SU:-  SU<-  -
  n)U)S:w  a\  T R.                  U)-  nDUDS-  nEWOUD-  S-  nFU)S-   * nGUGWPSU:S-
  U<S-
  4   -  UEUPSU:S-
  U<4   -  -   UFUPSU:S-
  U<S-
  4   -  -   UPSU:U<4'   Mp  SWPSU:U<4'   [!        SU<5       Hf  nHWPSU:UH4   T R%                  S'5      T R1                  SU<-  SUH-  -
  5      -  T R1                  U<UH-
  5      -  -  nIUPSU:U<4==   SU<UH-
  -  UI-  -  ss'   Mh     Mð     SWPSU:S(4'   SUPSU:S4'   SUPSU:SU:-  S-  S-   4'   GM<     [!        S(US-   5       Hd  nJ[!        S(U5       HQ  n:[!        S)[        SSU:-  S-  S-   5      5       H+  n<WJS:  d  U:S:  d  U<S:  d  U<SU:-  S-  :”  d  M#  SWPWJU:U<4'   M-     MS     Mf     [!        SUS-   5       H‘  nJ[!        SU5       H~  n:WKU:   S-   T l         [!        SSU:-  S-  S-   5       HU  n<SWJ-  S-
  WPUJS-
  U:S-
  U<S-
  4   -  SUU:-   S-
  -  UPUJS-
  U:S-
  U<S-
  4   -  -   n/U/UPUJS-
  U:S-
  U<S-
  4   -
  UPUJU:U<4'   MW     M€     M“     0 nQ0 nRST l         T R                  SUS-   -  5      nET R                  S*US-   -  U-  5      nST R                  UU-  [        R                  " T R                  5      -  5      S-
  nT[!        SU5       Hž  n<WSU<WT-  -
  T R                  T R1                  U<S-   5      5      S-  -   nU[        WEUUT R                  U5      -
  S-   5      S-   S-   S+-  WQU<'   [        UET R                  US-   5      UU-   T R                  U5      -
  S-   5      S-   S-   WRU<'   M      0 nV0 nWST l         T R                  SUS-   -  5      nET R                  S*US-   -  U-  5      nXT R                  UU-  [        R                  " T R                  5      -  5      S-
  nY[!        SU5       H›  n<WXU<WY-  -
  T R                  T R1                  U<S-   5      5      S-  -   nZ[        WEUZT R                  U5      -
  S-   5      S-   S+-  WVU<'   [        UET R                  US-   5      UZ-   T R                  U5      -
  S-   5      S-   S-   WWU<'   M�     0 n[[!        SUS-   5       H:  nJ[!        SU5       H'  n<[!        S(SU<-  S-  S-   5       H  n\SW[WJU<U\4'   M     M)     M<     [!        SUS-   5       H†  nJ[!        SU5       Hs  n<WQU<   T l         [!        SSU<-  S-  S-   5       HM  n\WCWJU<U\4   U9SU<-  SU\-  -
     -  U8SU<-  U\-
     -  W[UJU<U\4'   U[UJU<U\4   ST R2                  -  U\-  -  U[UJU<U\4'   MO     Mu     Mˆ     U 4S, jn]0 n^[!        SUS-   5       H:  nJ[!        SU5       H'  n<[!        S(SU<-  S-  S-   5       H  n\SW^WJU<U\4'   M     M)     M<     [        WQS   WVS   5      S-   T l         W]" U5      n_T R-                  U_5      * n`[!        SUS-   5       H›  na[!        SU5       Hˆ  n<WQU<   T l         [!        SSU<-  S-  S-   5       Hb  n\SW^WaU<U\4'   [!        SUaS-   5       HD  nJT R5                  WaUJ5      T R                  W`UJ5      -  W[UaUJ-
  U<W\4   -  nbW^UaU<U\4==   Ub-  ss'   MF     Md     MŠ     M�     0 nc[!        SUS-   5       H:  nJ[!        SU5       H'  n<[!        S(SU<-  S-  S-   5       H  n\SWcWJU<U\4'   M     M)     M<     [!        SUS-   5       H†  nJ[!        SU5       Hs  n<WVU<   T l         [!        SSU<-  S-  S-   5       HM  n\WPWJU<U\4   U9SU<-  SU\-  -
     -  U8SU<-  U\-
     -  WcUJU<U\4'   UcUJU<U\4   ST R2                  -  U\-  -  UcUJU<U\4'   MO     Mu     Mˆ     0 nd[!        SUS-   5       H:  na[!        SU5       H'  n<[!        S(SU<-  S-  S-   5       H  n\SWdWaU<U\4'   M     M)     M<     [!        SUS-   5       H›  na[!        SU5       Hˆ  n<WVU<   T l         [!        SSU<-  S-  S-   5       Hb  n\SWdWaU<U\4'   [!        SUaS-   5       HD  nJT R5                  WaUJ5      T R                  W`UJ5      -  WcUaUJ-
  U<W\4   -  nbWdUaU<U\4==   Ub-  ss'   MF     Md     MŠ     M�     [        WQS   WVS   5      S-   T l         0 neSUeS'   UeS   U-  UeS'   [        UQS   UVS   5      T l         [!        SU5       H  n<WeU<S-
     UeS   -  UeU<'   M     0 nf[!        SUS-   5       HP  na[!        SU5       H=  n<WQU<   T l         [!        SSU<-  S-  S-   5       H  n\W^WaU<U\4   WeU<   -  WfUaU<U\4'   M     M?     MR     0 ng[!        SUS-   5       HP  na[!        SU5       H=  n<WVU<   T l         [!        SSU<-  S-  S-   5       H  n\WdWaU<U\4   WeU<   -  WgUaU<U\4'   M     M?     MR     0 nh[!        SUS-   5       HP  na[!        SU5       H=  n:WRU:   T l         Sni[!        SSU:-  S-  S-   5       H  n<WiWfWaU:U<4   -  niM     WiWhWaU:4'   M?     MR     0 nj[!        SUS-   5       HP  na[!        SU5       H=  n:WWU:   T l         Sni[!        SSU:-  S-  S-   5       H  n<WiWgWaU:U<4   -  niM     WiWjWaU:4'   M?     MR     0 nkST l         [        R                  " T R                  5      U-  UU-  -  nlST l         T R                  S-US-   S-  -  Ul-  U-  5      nm[        UmT R                  SUS-   -  5      5      nmUmT l         [!        SUS-   5       H3  naSWkUa'   [!        SUS-   5       H  n<WkWa==   WhUaUU<-
  4   -  ss'   M     M5     0 nnST l         [        R                  " T R                  5      U-  UU-  -  noST l         T R                  S-US-   S-  -  Uo-  U-  5      np[        UpT R                  SUS-   -  5      5      npUpT l         [!        SUS-   5       H3  naSWnUa'   [!        SUS-   5       H  n<WnWa==   WjUaUU<-
  4   -  ss'   M     M5     ST l         S[        T R                  [        WkS   5      5      T R                  [        WnS   5      5      5      -  nqUSUq-  -  nrUT R-                  UST R                  -  -  5      -  nsS.T R                  SSUr-  T R                  Xu* 5      -  -   Us-  5      -   ntS.T R                  SSUr-  T R                  Xv* 5      -  -   Us-  5      -   nu[        UtUu5      T l         UST R                  -  -  nvUS-  T R-                  Uv5      -  US-  -
  T R                  S-  -
  nwT R7                  Uw* 5      nxW]" U5      nT R                  Xu* 5      nT R                  Xv* 5      n	SU1S-
  -  U-  Ux-  nySU1S-
  -  U	-  Ux-  nzST l         ST R                  U1T R                  U1SU-
  5      -   T R-                  U15      -  U-  5      -   n{ST R                  U1T R                  U1SU-
  5      -   T R-                  U15      -  U-  5      -   n|[        U{U|5      n}U}S-   T l          T R9                  US[;        U15      S-
  [!        SUS-   5      S/5      u  n~nST l         [        U~U   5      n€[        UkU   Uy-  5      n�[        SUT R                  SSU€-  S0U�-  -   -  5      -   5      n‚U‚T l         0 nƒ[!        SUS-   5       H  naW~Ua   WkUa   Wy-  -   WƒUa'   M     ST l         [        WU   5      n„[        WnU   Wz-  5      n…[        SUT R                  SSU„-  S0U…-  -   -  5      -   5      n†U†T l         0 n‡[!        SUS-   5       H+  naWUa   WnUa   Wz-  -   W‡Ua'   T R=                  U‡Ua   5      U‡Ua'   M-     UT l         WƒW‡4$ )1Nr   r   r   ç      @r\   r   r   r   ç½5°U‚Å@çŸ`<@çYÝê¹ñ?r   ç      à?ç      9@ú2Riemann-Siegel can not compute with such precisionr   é   é   ç     Às@é   éœ  é,   é   r!   r   r    r   é   é   r   ©r   r   éþÿÿÿ©r   r   r   ©r   r   r   ©r   r   r   ©r   r   r„   ©r   r   r   ©r   r   r   ©r   r   r   z1.0r„   éýÿÿÿéD   ç      ø?c                 ó„   >• TR                   nUS-   Tl         TR                  U STR                  -  -  5      nUTl         U$ ©Nr   ©r%   r,   r'   ©ÚtÚwpÚaar	   s      €r
   Útrunc_aÚRzeta_simul.<locals>.trunc_a0  ó=   ø€ Ø�X‰XˆØ˜‘6ˆŒØ�X‰X�a˜˜3Ÿ6™6™‘lÓ#ˆØˆŒØˆ	r   çš™™™™™@r   Té   )r%   Ú_imÚ_rer,   r'   rl   r#   ÚmathÚpowÚgammar"   ÚabsÚNotImplementedErrorri   rq   rm   r(   Úfloorr)   rL   rZ   Úcopyr]   r&   r*   r+   ÚbinomialÚexpjÚ_zetasumÚintÚconj)ˆr	   ÚsÚderÚ	wpinitialr“   ÚxsigmaÚysigmara   ÚxasigmaÚyasigmaÚxA1ÚyA1r/   Úeps1Úxeps2Úyeps2ÚxbÚxcr_   rb   ÚybÚycÚyAÚyB1rc   ÚyLÚLÚxeps3Úyeps3r`   Úyeps4rh   ÚyMÚMrp   Úh4r.   ÚjvalueÚeps5Úxforeps5Úyforeps5rf   Úxaux1Úyaux1rd   re   ÚtwentyÚauxÚwpfpÚNÚpÚnumÚ
differenceÚeps6ÚccÚcontÚpipowersÚFpr6   ÚsumPr?   ÚxwpdÚd1Úxd2ÚxconstÚxd3ÚxpsigmaÚxdÚm1Úc1Úc2Úc3ÚrÚaddrD   ÚywpdÚyd2ÚyconstÚyd3ÚypsigmaÚydÚxwptcoefÚxwptermÚxc2Úxc4Úxc3ÚywptcoefÚywptermÚyc2Úyc4Úyc3Ú	xfortcoefÚellr–   Úxtcoefr•   ÚlaÚchiÚtcoefterÚ	yfortcoefÚytcoefÚavÚxtvÚytvÚxtermÚteÚytermÚxrssumÚxrsboundÚxwprssumÚyrssumÚyrsboundÚywprssumÚA2Úeps8ÚTÚxwps3Úywps3ÚtpiÚargÚUÚxS3ÚyS3ÚxwpsumÚywpsumÚwpsumÚxS1ÚyS1ÚxabsS1ÚxabsS2ÚxwpendÚxrzÚyabsS1ÚyabsS2ÚywpendÚyrzsˆ   `                                                                                                                                       r
   ÚRzeta_simulr  ø   s  ø€ à—‘€Ið 	�‰�‹
€AØ�W‰W�Q‹Z€FØ�‰Z€Fð €C„HØ�‰��A�c—f‘f‘H‘Ó€AØ‰k€GØ‰k€Gð 	�	‰	�!�S—W‘W˜WÓ% aÑ'Ó(€CØ�	‰	�!�S—W‘W˜WÓ% aÑ'Ó(€Cð �)‰)�A˜	�zÓ
"€CØˆr‰6€DØ‰I�b‰L€EØ‰I�b‰L€Eð €C„HØ�ƒzØˆÜ�XŠX�a˜Ó Ñ'ˆä�XŠX�a˜ÓˆØ‰àˆÜ�XŠX�a˜˜Ó  Ñ(ˆÜ�XŠX�a˜˜Ó ˆØˆà�ƒzØˆÜ�XŠX�a˜Ó Ñ'ˆä�XŠX�a˜ÓˆØ‰àˆÜ�XŠX�a˜˜Ó  Ñ(ˆÜ�XŠX�a˜˜Ó ˆØˆð
 €C„HØ	
€BØ
ˆB‰$ˆs�y‰y˜˜C™Ó Ñ
  3§9¡9¨R°©T°2°#Ó#6Ñ
6¸%Ó
?Ø�‰Tˆð ˆB‰$ˆs�y‰y˜˜C™Ó Ñ
  3§9¡9¨R°©T°2°#Ó#6Ñ
6¸%Õ
?ä	ˆQˆr‹€BØ	
€BØ
ˆB‰$ˆs�y‰y˜˜C™Ó Ñ
  3§9¡9¨R°©T°2°#Ó#6Ñ
6¸%Ó
?Ø�‰Tˆð ˆB‰$ˆs�y‰y˜˜C™Ó Ñ
  3§9¡9¨R°©T°2°#Ó#6Ñ
6¸%Õ
?ä	ˆQˆr‹€Bð 
ˆ2‰��1‘�Q‘�s‘Ó	  "¡ Q¡ v¡¨a£´S¸³[À1ÀRÁ4Ó5GØ	
ˆ2‰��1‘�Q‘�s‘Ó	  "¡ Q¡ v¡¨a£´S¸³[À1ÀRÁ4Ó5GØˆŒÜ!Ð"VÓWÐWô 	ˆB�‹€Að �A�b‘D‰\€EØ�A�b‘D‰\€Eð
 �1�R‘4‰L€EØ�1�R‘4‰L€Eô 
�#�r˜5 ! S¨"Ó	-€BÜ	�#�r˜5 ! S¨"Ó	-€BÜˆB�‹€Aô 
�c˜2˜u a¨¨bÓ	1€BÜ	�c˜2˜u a¨¨bÓ	1€BÜ	ˆR�‹€BØ
€AØ�—‘‰h˜‰]˜SŸY™Y q¨¡s›^Ñ+€FØ
�2‹+Øˆa‰CˆØ�C—F‘F‘(˜FÑ" 1Ñ$ˆð �2�+ð 
€DÜ�w‰w”t—w‘w‰˜sÑ" 1Ñ$€HÜ�w‰w”t—w‘w‰˜sÑ" 1Ñ$€HÜ�1�RŽ[ˆÜ—’˜ 1 Q¡3Ó'¨¨b©Ñ1ˆÜ—’˜ 1 Q¡3Ó'¨¨b©Ñ1ˆÜ�5˜%Ó ˆØ�y‰y˜˜1™‹~˜cŸi™i¨¨#©¨c©	Ó2Ñ2ˆÜ�yŠy˜‹ˆØ˜‘9œS  uÓ-Ñ-ˆˆQ‹ñ ô ��1‘�Q‘˜‹^˜AÑ€FØ
ˆq‰&€CØ�7‰7�2�a‘4‹=€DÜ�1�VŽ_ˆÜ�4˜Ÿ™  S§Y¡Y¨q°©s£^Ñ!3°D¸±GÑ!;Ó<Ó=Šñ ð
 �c—g‘g˜a“jÑ  Ñ#€C„HØ�‰��A�c—f‘f‘H‘Ó€AØ�	‰	�!‹€AØ	ˆ!ˆQˆq‰S‰'‰	€Að
 	�	‰	�!�S—W‘W˜S“\ 4Ñ'Ñ(Ó)€CØ�c—g‘g˜c“l DÑ(Ñ)¨#Ñ-€JØ�SÓØ‰à�!‰eˆØ�‰�C˜3Ÿ7™7 3›<¨4¨%Ñ0Ñ1Ó2€Að �9‰9�S—[‘[  3§6¡6¡Ó*¨AÓ.°·	±	¸!¸A¹#³¸qÑ0@ÀÑ0BÑC€Dð 
€BØ€DÜ˜#˜q $Ó'�N€Dˆ(Ø‡y�yƒ{€BØ	€BÜ�1�a˜‘c˜!‘eŽ_ˆØˆˆ1‹ñ à	€BØ€C„HÜ�1�Q�q‘SŽ\ˆØˆÜ�q˜‘s˜1‘u˜Q‘w˜r "Ö%ˆAØ˜1‘H˜r !™uÑ$ŠDñ &àˆˆ1‰ä�q˜˜1™˜Q™˜q™Ö!ˆAØ�q‘S˜2˜a ™c™7‘NˆBˆq‹Eó "ñ ð 
€DÜ	ˆQˆs�w‰w�r˜!‘t˜A‘v‹Ó	€BØ
ˆS�W‰W�aœ˜F›‘m RÑ'Ó(Ñ
(¨¯©°«Ñ
7¸Ñ
9€CØ�V‰V�A�s—v‘v˜cŸf™f‘} Q‘ qÑ(¨Ñ,¨SÑ0Ñ1Ó2°AÑ5€FÜ�1�QŽZˆØ�g‰g�c—h‘h˜sŸy™y¨¨3©Ó/Ó0Ó1°#·)±)¸A¸f¹HÓ2EÑEÀcÑIˆÜ�C˜“ˆˆQ‹ñ ð
 �A‰w�r‰z€C„HØ��6‘‰l€GØ	€BØ€B€v�J˜Q�"�V‘*¨!  5¡	°q¨R°©YØ€B€w�K˜q�2�g‘;¨Q " V¡*¸1°°6±
Ü�1�Q�ZˆØ˜‘7˜2‘:ˆŒÜ�q˜˜1™˜a™ ™Ö"ˆAØ�!‘�A�a‘C‘ˆAØ�!‹tØ—W‘W˜Q‘Y�Ø�q‘D�Ø˜B‘J ‘>�Ø�q‘S�6�Ø˜R  ! A¡# a¨¡c 	™]Ñ*¨2¨b°°1°Q±3°q°©k©>Ñ9¸"¸RÀÀ!ÀAÁ#ÀaÈÁcÀ	¹]Ñ:JÑJ��1�Q�q�5“	à��1�Q�q�5‘	Ü˜q ž�AØ˜1˜Q˜q˜5™	 3§7¡7¨5£>°#·'±'¸!¸A¹#¸aÀ¹c¹'Ó2BÑ#BÀ3Ç7Á7È1ÈQÉ3Ã<Ñ#OÑP�CØ�q˜˜1�u“I 2¨¨1©¡+¨sÑ!2Ñ2•Ió $ñ #ð ˆˆ1ˆQˆrˆ6‰
 �b˜˜1˜R˜‘j°Q B q¨¨1¨Q©3°©6°!©8 |Ô$4ñ ô  �B�s˜1‘uŽoˆÜ�r˜!–ˆAÜ˜2œc ! A a¡C¨¡F¨1¡H›oÖ.�Ø˜“T˜a ›c a¨£c¨a°°!±°Q±­hØ!"�B�r˜!˜A�v“Jó /ó ñ ô
 �A�c˜!‘eŽnˆÜ�q˜–ˆAØ˜A‘w˜r‘zˆCŒHÜ˜1˜Q˜q™S !™V A™XÖ&�Ø�r‘T˜!‘V˜R  1¡ Q q¡S¨¨1© Ñ-Ñ-¨a°¸±¸!±©n¸RÀÀ1ÁÀQÀqÁSÈÈ1ÉÀÑ=MÑ.MÑM�Ø  2 b¨¡d¨1¨Q©3¨q°©s lÑ#3Ñ3��2�a˜�6“
ó 'ó ñ ð 
€DÜ	ˆQˆs�w‰w�r˜!‘t˜A‘v‹Ó	€BØ
ˆS�W‰W�aœ˜F›‘m RÑ'Ó(Ñ
(¨¯©°«Ñ
7¸Ñ
9€CØ�V‰V�A�s—v‘v˜cŸf™f‘} Q‘ qÑ(¨Ñ,¨SÑ0Ñ1Ó2°AÑ5€FÜ�1�QŽZˆØ�g‰g�c—h‘h˜sŸy™y¨¨3©Ó/Ó0Ó1°#·)±)¸A¸f¹HÓ2EÑEÀcÑIˆÜ�C˜“ˆˆQ‹ñ ð
 �A‰w�r‰z€C„HØ��6‘‰l€GØ	€BØ€B€v�J˜Q�"�V‘*¨!  5¡	°q¨R°©YØ€B€w�K˜q�2�g‘;¨Q " V¡*¸1°°6±
Ü�1�Q�ZˆØ˜‘7˜2‘:ˆŒÜ�q˜˜1™˜a™ ™Ö"ˆAØ�!‘�A�a‘C‘ˆAØ�!‹tØ—W‘W˜Q‘Y�Ø�q‘D�Ø˜B‘J ‘>�Ø�q‘S�6�Ø˜R  ! A¡# a¨¡c 	™]Ñ*¨2¨b°°1°Q±3°q°©k©>Ñ9¸"¸RÀÀ!ÀAÁ#ÀaÈÁcÀ	¹]Ñ:JÑJ��1�Q�q�5“	à��1�Q�q�5‘	Ü˜q ž�AØ˜1˜Q˜q˜5™	 3§7¡7¨5£>°#·'±'¸!¸A¹#¸aÀ¹c¹'Ó2BÑ#BÀ3Ç7Á7È1ÈQÉ3Ã<Ñ#OÑP�CØ�q˜˜1�u“I 2¨¨1©¡+¨sÑ!2Ñ2•Ió $ñ #ð ˆˆ1ˆQˆrˆ6‰
 �b˜˜1˜R˜‘j°Q B q¨¨1¨Q©3°©6°!©8 |Ô$4ñ ô" �B�s˜1‘uŽoˆÜ�r˜!–ˆAÜ˜2œc ! A a¡C¨¡F¨1¡H›oÖ.�Ø˜“T˜a ›c a¨£c¨a°°!±°Q±­hØ!"�B�r˜!˜A�v“Jó /ó ñ ô
 �A�c˜!‘eŽnˆÜ�q˜–ˆAØ˜A‘w˜r‘zˆCŒHÜ˜1˜Q˜q™S !™V A™XÖ&�Ø�r‘T˜!‘V˜R  1¡ Q q¡S¨¨1© Ñ-Ñ-¨a°¸±¸!±©n¸RÀÀ1ÁÀQÀqÁSÈÈ1ÉÀÑ=MÑ.MÑM�Ø  2 b¨¡d¨1¨Q©3¨q°©s lÑ#3Ñ3��2�a˜�6“
ó 'ó ñ ð €HØ€GØ€C„HØ	�‰��Q�q‘S‘Ó	€BØ
�'‰'�"�a˜‘c‘(˜2‘+Ó
€CØ
�'‰'�#�a‘%œŸ	š	 #§&¡&Ó)Ñ)Ó
*¨1Ñ
,€CÜ�1�QŽZˆØ�A�c‘E‰k˜#Ÿ'™' #§'¡'¨!¨C©%£.Ó1°"Ñ4Ñ4ˆÜ˜2˜c #§'¡'¨%£.Ñ0°Ñ2Ó3°AÑ5°rÑ9¸3Ñ>ˆ�‰Ü˜"˜SŸW™W Q q¡S›\¨#Ñ-¨c¯g©g°e«nÑ<¸QÑ>Ó?ÀÑAÀ2ÑEˆ�‹
ñ ð €HØ€GØ€C„HØ	�‰��Q�q‘S‘Ó	€BØ
�'‰'�"�a˜‘c‘(˜2‘+Ó
€CØ
�'‰'�#�a‘%œŸ	š	 #§&¡&Ó)Ñ)Ó
*¨1Ñ
,€CÜ�1�QŽZˆØ�A�c‘E‰k˜#Ÿ'™' #§'¡'¨!¨C©%£.Ó1°"Ñ4Ñ4ˆÜ˜B˜s 3§7¡7¨5£>Ñ1°!Ñ3Ó4°bÑ8¸#Ñ=ˆ�‰Ü˜"˜SŸW™W Q q¡S›\¨#Ñ-¨c¯g©g°e«nÑ<¸QÑ>Ó?ÀÑAÀ2ÑEˆ�‹
ñ ð €IÜ�A�c˜!‘eŽnˆÜ�q˜–ˆAÜ˜R  !¡ Q¡ q¡Ö)�Ø$%�	˜"˜Q˜s˜(Ó#ó *ó ñ ô �A�c˜!‘eŽnˆÜ�q˜–ˆAØ ‘{ˆCŒHÜ˜Q˜q ™s A™v a™xÖ(�Ø$& r¨!¨C x¡L°°A°a±C¸¸#¹±I±Ñ$>¸xÈÈ!ÉÈCÉÑ?PÑ$P�	˜"˜Q˜s˜(Ñ#Ø$-¨b°°3¨hÑ$7¸!¸C¿E¹E¹'ÀC¹Ñ$H�	˜"˜Q˜s˜(Ó#ó )ó ñ õð €FÜ�A�c˜!‘eŽnˆÜ�q˜–ˆAÜ˜R  !¡ Q¡ q¡Ö)�Ø!"��r˜!˜C�xÓ ó *ó ñ ô �8˜A‘;˜x¨™{Ó+¨AÑ-€C„HÙ�‹
€BØ
�&‰&�‹*ˆ€Bä�Q�s˜1‘uŽ~ˆÜ�q˜–ˆAØ ‘{ˆCŒHÜ˜Q˜q ™s A™v a™xÖ(�Ø#$��s˜1˜S�yÑ!Ü  3 q¡5ž/�BØ Ÿ\™\¨#¨bÓ1°#·)±)¸B¸rÓ2BÑBÀ9ÈSÐQSÉVÐTUÐVYÈ\ÑCZÑZ�HØ˜3˜q ˜9Ó%¨Ñ1Õ%ó *ó )ó ñ ð €IÜ�A�c˜!‘eŽnˆÜ�q˜–ˆAÜ˜R  !¡ Q¡ q¡Ö)�Ø$%�	˜"˜Q˜s˜(Ó#ó *ó ñ ô �A�c˜!‘eŽnˆÜ�q˜–ˆAØ ‘{ˆCŒHÜ˜Q˜q ™s A™v a™xÖ(�Ø$& r¨!¨C x¡L°°A°a±C¸¸#¹±I±Ñ$>¸xÈÈ!ÉÈCÉÑ?PÑ$P�	˜"˜Q˜s˜(Ñ#Ø$-¨b°°3¨hÑ$7¸!¸C¿E¹E¹'ÀC¹Ñ$H�	˜"˜Q˜s˜(Ó#ó )ó ñ ð €FÜ�Q�s˜1‘uŽ~ˆÜ�q˜–ˆAÜ˜R  !¡ Q¡ q¡Ö)�Ø"#��s˜1˜S�yÓ!ó *ó ñ ô �Q�s˜1‘uŽ~ˆÜ�q˜–ˆAØ ‘{ˆCŒHÜ˜Q˜q ™s A™v a™xÖ(�Ø#$��s˜1˜S�yÑ!Ü  3 q¡5ž/�BØ Ÿ\™\¨#¨bÓ1°#·)±)¸B¸rÓ2BÑBÀ9ÈSÐQSÉVÐTUÐVYÈ\ÑCZÑZ�HØ˜3˜q ˜9Ó%¨Ñ1Õ%ó *ó )ó ñ ô �8˜A‘; ¨¡Ó,¨QÑ.€C„HØ	€BØ€B€q�EØˆq‰E�!‰G€B€q�Eä�8˜A‘;˜x¨™{Ó+€C„HÜ�1�QŽZˆØ�1�Q‘3‘˜"˜Q™%‘ˆˆ1‹ñ ð €CÜ�Q�s˜1‘uŽ~ˆÜ�q˜–ˆAØ ‘{ˆCŒHÜ˜Q˜q ™s A™v a™xÖ(�Ø!'¨¨A¨c¨	Ñ!2°B°q±EÑ!9��C˜˜#�I“ó )ó ñ ð €CÜ�Q�s˜1‘uŽ~ˆÜ�q˜–ˆAØ ‘{ˆCŒHÜ˜Q˜q ™s A™v a™xÖ(�Ø!'¨¨A¨c¨	Ñ!2°B°q±EÑ!9��C˜˜#�I“ó )ó ñ ð €EÜ�Q�s˜1‘uŽ~ˆÜ�q˜–ˆAØ˜q‘zˆCŒHØˆBÜ˜1˜a ™c 1™f Q™hÖ'�Ø�c˜#˜a ˜'‘lÑ"’ñ (àˆE�#�a�%‹Ló ñ ð €EÜ�Q�s˜1‘uŽ~ˆÜ�q˜–ˆAØ˜q‘zˆCŒHØˆBÜ˜1˜a ™c 1™f Q™hÖ'�Ø�c˜#˜a ˜'‘lÑ"’ñ (àˆE�#�a�%‹Ló ñ ð €FØ€C„HÜ�yŠy˜Ÿ™Ó  2Ñ%¨¨1©Ñ-€HØ€C„HØ�w‰w�s˜Q˜q™S 1™H‘~ hÑ.°Ñ6Ó7€HÜ�8˜SŸW™W R¨¨1©¡XÓ.Ó/€HØ€C„HÜ�Q�s˜1‘uŽ~ˆØˆˆs‰Ü�q˜˜1™–ˆAØ�3‹K˜5  Q q¡S ™>Ñ)�Kó ñ ð €FØ€C„HÜ�yŠy˜Ÿ™Ó  2Ñ%¨¨1©Ñ-€HØ€C„HØ�w‰w�s˜Q˜q™S 1™H‘~ hÑ.°Ñ6Ó7€HÜ�8˜SŸW™W R¨¨1©¡XÓ.Ó/€HØ€C„HÜ�Q�s˜1‘uŽ~ˆØˆˆs‰Ü�q˜˜1™–ˆAØ�3‹K˜5  Q q¡S ™>Ñ)�Kó ñ ð €C„HØ	
ŒS�—‘œ˜V A™Y›Ó(¨#¯'©'´#°f¸Q±i³.Ó*AÓBÑ	C€BØ��"‘‰:€DØ	ˆ3�6‰6�!�Q�s—v‘v‘X‘,ÓÑ€AØ�—‘˜!˜Q˜t™V S§Y¡Y¨q°Ó%9Ñ9Ñ9¸1Ñ<Ó=Ñ=€EØ�—‘˜!˜Q˜t™V S§Y¡Y¨q°Ó%9Ñ9Ñ9¸1Ñ<Ó=Ñ=€Eä�5˜%Ó €C„Hà
ˆQˆs�v‰v‰X‰,€CØˆQ‰3�—‘�s“Ñ
˜Q˜q™SÑ
! #§&¡&¨¡(Ñ
*€CØ�‰�#�‹€AÙ�‹
€AØ�i‰i˜˜7Ó#€GØ�i‰i˜˜7Ó#€GØ�!�A‘#‰;˜'Ñ
! AÑ
%€CØ�!�A‘#‰;˜'Ñ
! AÑ
%€Cð €C„HØ�—‘˜!˜CŸI™I a¨¨&©Ó1Ñ1°3·6±6¸!³9Ñ<¸dÑBÓCÑC€FØ�—‘˜!˜CŸI™I a¨¨&©Ó1Ñ1°3·6±6¸!³9Ñ<¸dÑBÓCÑC€FÜ�˜Ó€Eà�b‰y€C„Hðð$ �|‰|˜A˜q¤# a£&¨¡(¬E°!°C¸±E«N¸DÓA�H€Cˆð €C„HÜ��S‘‹]€FÜ�˜‘˜sÑ"Ó#€FÜ��I˜cŸg™g a¨¨6©°!°F±(Ñ):Ñ&;Ó=Ñ=Ó?€Fà€C„HØ
€CÜ�Q�s˜1‘uŽ~ˆØ�s‘8˜F 3™K¨™OÑ+ˆˆC‹ñ ð €C„HÜ��S‘‹]€FÜ�˜‘˜sÑ"Ó#€FÜ��I˜cŸg™g a¨¨6©°!°F±(Ñ):Ñ&;Ó=Ñ=Ó?€Fà€C„HØ
€CÜ�Q�s˜1‘uŽ~ˆØ�s‘8˜F 3™K¨™OÑ+ˆˆC‰Ø—8‘8˜C ™HÓ%ˆˆC‹ñ ð €C„HØ�ˆ8€Or   c                 óV   ^ • [        U5      nT R                  nT R                  U5      nT R                  U5      nST l        T R	                  UST R
                  -  -  5      nT R                  Xv5      nT R                  ST R                  U5      S-
  5      n	T R                  SU* 5      n
U
S-  nX©-  S-  nST l        US:”  a6  Sn[        R                  " SU5      S	-  n[        R                  " SU5      nSnO7S
n[        R                  " SU* 5      S	-  n[        R                  " SU* 5      nSnST l        SnSU-  T R                  US-  5      -  T R                  X×-  U* 5      -  U:¼  a;  US-   nSU-  T R                  US-  5      -  T R                  X×-  U* 5      -  U:¼  a  M;  [        SU5      nSU-  SU-  U-  S-  :¼  d!  SU-  S-   U-   S:  d  [        U5      US-  :”  a  UT l        [        S5      eUSU-  -  nUSU-  -  n[        T UUUUU5      n0 n[        USU-  S-
  5       H  nSUU'   M
     USU-  -  nT R
                  T R
                  -  U-  U-  T R	                  S5      -  [        R                  -  [        R                  -  nUT R                  UUS-  -  US-
  S-  5      -  U-  n[!        UU5      nSnST R
                  -  U-  T R                  US-   5      -  nUU:”  a"  US-   nST R
                  -  U-  U-  nUU:”  a  M"  0 n[        R
                  [        R
                  -  U-  U-  n[        SS5       Hp  n[        R                  " UUS-  5      SU-  -  nT R                  US-   5      T R                  US-  S-   5      -  n [        R                  " U 5      n UU -  U-  UU'   Mr     [!        SU-  S-
  S5      S-   n!SU-  n"T R                  SU-  5      n#[        SU!5       H9  n[        U#T R                  U"T R                  US-   5      -  UU   -  5      5      n#M;     U#T R                  U5      -   S-   T l        T R	                  UST R
                  -  -  5      nT R#                  U5      n$SSUU$-
  -  -
  n%T R#                  U%T R%                  S5      U#-  -  5      n&U%T R%                  S5      U#-  -  U&-
  n'U'S:  a  U&n&OU&S-   n&T R'                  U&T R%                  S5      U#* -  -  5      n%T R                  ST R
                  -  U5      T R                  US-   5      S-  U-  -  n(0 n)0 n*[)        T UU(5      u  n*n+U*R+                  5       n)0 n[        USU-  S-
  5       H  nSUU'   M
     U#T l        [        SUS-   5       H_  nSn,[        SU-  U-
  S-
  SS5       H  n-U,U%-  U)U-   -   n,M     U,UU'   [        SSU-  U-
  S-
  5       H  n-U-S-   U)U-S-      -  U)U-'   M     Ma     0 n.[        ST R                  SU-  U-  5      5      n/ST R                  S[        U5      -   U-  5      -   T R                  U5      -
  S-
  n0T R-                  ST R
                  T R
                  -  U-  U-  U-  U-  -  5      S-  n1[        SU5       H\  nT R                  T R	                  T R                  US-
  5      5      5      T R#                  UU1-  5      -
  U0-   n2[        U2U/5      U.U'   M^     U.S   S-   T l        SSU-  -
  n30 n4SU4S'   SU4S '   SU4S!'   SU4S"'   SU4S#'   SU4S$'   SU4S%'   SU4S&'   [        SU5       GH4  nU.U   S-   T l        [        SSU-  S-  S-   5       Hé  n-SU-  SU--  -
  nUS:w  a\  T R.                  U-  n5U5S-  n6U3U5-  S-  n7US-   * n8U8U4SUS-
  U-S-
  4   -  U6U4SUS-
  U-4   -  -   U7U4SUS-
  U-S-
  4   -  -   U4SUU-4'   Mp  SU4SUU-4'   [        SU-5       Ha  n9U4SUU94   T R.                  T R1                  SU--  SU9-  -
  5      -  T R1                  U-U9-
  5      -  -  n:U4SUU-4==   SU-U9-
  -  U:-  -  ss'   Mc     Më     SU4SUS'4'   SU4SUS4'   SU4SUSU-  S-  S-   4'   GM7     [        S'US-   5       Hd  n;[        S'U5       HQ  n[        S([        SSU-  S-  S-   5      5       H+  n-U;S:  d  US:  d  U-S:  d  U-SU-  S-  :”  d  M#  SU4U;UU-4'   M-     MS     Mf     [        SUS-   5       H‘  n;[        SU5       H~  nU.U   S-   T l        [        SSU-  S-  S-   5       HU  n-SU;-  S-
  U4U;S-
  US-
  U-S-
  4   -  SUU-   S-
  -  U4U;S-
  US-
  U-S-
  4   -  -   n"U"U4U;S-
  US-
  U-S-
  4   -
  U4U;UU-4'   MW     M€     M“     0 n<0 n=ST l        T R                  SUS-   -  5      n6T R                  S)US-   -  U-  5      n7T R                  UU-  [        R                  " T R
                  5      -  5      S-
  n>[        SU5       H›  n-U7U-U>-  -
  T R                  T R1                  U-S-   5      5      S-  -   n8[        U6U8T R                  U5      -
  S-   5      S-   S-   U<U-'   [        U6T R                  US-   5      U8-   T R                  U5      -
  S-   5      S-   S-   U=U-'   M�     0 n?U H:  n;[        SU5       H'  n-[        S'SU--  S-  S-   5       H  n@SU?U;U-U@4'   M     M)     M<     U H†  n;[        SU5       Hs  n-U<U-   T l        [        SSU--  S-  S-   5       HM  n@U4U;U-U@4   USU--  SU@-  -
     -  U+SU--  U@-
     -  U?U;U-U@4'   U?U;U-U@4   ST R2                  -  U@-  -  U?U;U-U@4'   MO     Mu     Mˆ     U 4S* jnA0 nBU H:  nC[        SU5       H'  n-[        S'SU--  S-  S-   5       H  n@SWBWCU-U@4'   M     M)     M<     U<S   S-   T l        WA" U5      nDT R-                  UD5      * nEU HŽ  nC[        SU5       H{  n-U<U-   T l        [        SSU--  S-  S-   5       HU  n@SWBWCU-U@4'   [        SUCS-   5       H7  n;T R5                  WCU;5      WEU;-  -  U?UCU;-
  U-W@4   -  nFWBUCU-U@4==   UF-  ss'   M9     MW     M}     M�     U<S   S-   T l        0 nGSUGS'   UGS   U-  UGS'   U<S   T l        [        SU5       H  n-WGU-S-
     UGS   -  UGU-'   M     0 nHU HP  nC[        SU5       H=  n-U<U-   T l        [        SSU--  S-  S-   5       H  n@WBWCU-U@4   WGU-   -  WHUCU-U@4'   M     M?     MR     0 nIU HP  nC[        SU5       H=  nU=U   T l        SnJ[        SSU-  S-  S-   5       H  n-WJWHWCUU-4   -  nJM     WJWIWCU4'   M?     MR     0 nKST l        [        R                  " T R
                  5      U-  X×-  -  nLST l        T R                  S+US-   S-  -  UL-  U-  5      nM[        UMT R                  SUS-   -  5      5      nMUMT l        U H3  nCSWKUC'   [        SUS-   5       H  n-WKWC==   WIUCUU--
  4   -  ss'   M     M5     ST l        ST R                  WKS   5      -  nNU
SUN-  -  nOUT R-                  UST R
                  -  -  5      -  nPS,T R                  SSUO-  T R                  Xv* 5      -  -   UP-  5      -   nQUQT l        UST R
                  -  -  nRUS-  T R-                  UR5      -  US-  -
  T R
                  S-  -
  nST R7                  US* 5      nTWA" U5      nT R                  Xv* 5      nSU$S-
  -  U-  UT-  nUST l        ST R                  U$T R                  U$SU-
  5      -   T R-                  U$5      -  U-  5      -   nVUVS-   T l         T R9                  US[;        U$5      S-
  U5      S   nWST l        [        UWU   5      nX[        UKU   UU-  5      nY[        SUT R                  SSUX-  S-UY-  -   -  5      -   5      nZUZT l        0 n[U H  nCWWUC   WKUC   WU-  -   W[UC'   M     UT l        W[$ ).aa  
Computes several derivatives of the auxiliary function of Riemann `R(s)`.

**Definition**

The function is defined by

.. math ::

    \begin{equation}
    {\mathop{\mathcal R }\nolimits}(s)=
    \int_{0\swarrow1}\frac{x^{-s} e^{\pi i x^2}}{e^{\pi i x}-
    e^{-\pi i x}}\,dx
    \end{equation}

To this function we apply the Riemann-Siegel expansion.
r   r   r   rs   r\   r   r   r   rt   ru   rv   r   rw   rx   ry   r   rk   rz   r{   r|   r}   r~   r   r€   r   r    r   r�   r‚   r   rƒ   r…   r†   r‡   rˆ   r‰   rŠ   r‹   r„   rŒ   r�   c                 ó„   >• TR                   nUS-   Tl         TR                  U STR                  -  -  5      nUTl         U$ r�   r‘   r’   s      €r
   r–   ÚRzeta_set.<locals>.trunc_aô  r˜   r   r™   r   rš   )r"   r%   r›   rœ   r,   r'   rl   r#   r�   rž   rŸ   r    r¡   ri   r(   Úerm   r¢   r)   rL   rZ   r£   r]   r&   r*   r+   r¤   r¥   r¦   r§   )\r	   r©   Úderivativesrª   r«   r“   Úsigmara   ÚasigmaÚA1r/   r²   Úeps2ÚbrF   ÚAÚB1r¼   Úeps3Úeps4rÁ   rÔ   r6   rn   ro   rp   r.   rÃ   rÄ   Úforeps5rf   rd   re   rÉ   rÊ   rË   rÌ   rÍ   rÎ   rÏ   rÐ   rÑ   rÒ   rÓ   rÕ   r?   Úwpdr×   Úd2ÚconstÚd3ÚpsigmaÚdrÝ   rÞ   rß   rà   rá   râ   rD   ÚwptcoefÚwptermÚc4Úfortcoefrô   r–   Útcoefr÷   r•   rö   rø   rû   ÚtvÚtermrÿ   ÚrssumÚrsboundÚwprssumr  r  r	  Úwps3r  r  r  ÚS3r  ÚS1ÚabsS1ÚabsS2ÚwpendÚrzs\   `                                                                                           r
   Ú	Rzeta_setrE  þ  sc  ø€ ô$ ˆkÓ
€Cð —‘€Ià�‰�‹
€AØ�G‰G�A‹J€Eà€C„HØ�‰��A�c—f‘f‘H‘Ó€AØ�Y‰Y�qÓ €Fà	�‰�1�c—g‘g˜f“o aÑ'Ó	(€Bà
�)‰)�A˜	�zÓ
"€CØˆr‰6€DØ‰8�B‰;€Dð €C„HØˆqƒyØˆÜ�HŠH�Q�uÓ˜gÑ%ˆä�HŠH�Q�uÓˆØ‰àˆÜ�HŠH�Q˜�vÓ˜wÑ&ˆÜ�HŠH�Q˜�vÓˆØˆð €C„HØ	€AØ
ˆA‰#ˆc�i‰i˜˜#™ÓÑ
 §¡¨1©3°¨rÓ!2Ñ
2°dÓ
:Øˆa‰Cˆð ˆA‰#ˆc�i‰i˜˜#™ÓÑ
 §¡¨1©3°¨rÓ!2Ñ
2°dÕ
:äˆAˆa‹€Að 
ˆ1‰��!‘�A‘�c‘	Ó	˜q ™s 1™u U™{¨1›}´#°e³*¸aÀ¹dÓ2BàˆŒÜ!Ð"VÓWÐWð �!�A‘#‰J€Dð
 ��1‘‰:€Dô 	��a˜˜q " aÓ(€AØ	€BÜ�1�a˜‘c˜!‘eŽ_ˆØˆˆ1‹ñ ð 
ˆs�1‰u‰€BØ	�‰�—‘‰�rÑ	˜!Ñ	˜SŸX™X a›[Ñ	(¬¯©Ñ	/´·±Ñ	6€BØ	ˆc�i‰i˜˜A˜q™D™ 1 Q¡3¨¡'Ó*Ñ	*¨QÑ	.€BÜ	ˆR�‹€BØ€AØ�—‘‰h˜‰]˜SŸY™Y q¨¡s›^Ñ+€FØ
�2‹+Øˆa‰CˆØ�C—F‘F‘(˜FÑ" 1Ñ$ˆð �2�+ð 
€DÜ�g‰g”d—g‘g‰o˜bÑ  Ñ"€GÜ�1�RŽ[ˆÜ�xŠx˜  1¡Ó% t¨A¡vÑ.ˆØ�y‰y˜˜1™‹~˜cŸi™i¨¨#©¨c©	Ó2Ñ2ˆÜ�yŠy˜‹ˆØ�t‘)˜D‘.ˆˆQ‹ñ	 ô ��1‘�Q‘˜‹^˜AÑ€FØ
ˆq‰&€CØ�7‰7�2�a‘4‹=€DÜ�1�fÖˆÜ�4˜Ÿ™  S§Y¡Y¨q°©s£^Ñ!3°D¸±GÑ!;Ó<Ó=Šñ ð �c—g‘g˜a“jÑ  2Ñ%€C„HØ�‰��A�c—f‘f‘H‘Ó€AØ�	‰	�!‹€AØ	ˆ!ˆQˆq‰S‰'‰	€Að �)‰)�A�s—w‘w˜q“z 4Ñ'Ñ(Ó
)€CØ�c—g‘g˜a“j $Ñ&Ñ'¨Ñ+€JØ�CÓØ‰à�!‰eˆØ�‰�C˜3Ÿ7™7 1›:¨¨Ñ.Ñ/Ó0€Að �9‰9�Q�s—v‘v‘X˜qÓ! 3§9¡9¨Q¨q©S£>°!Ñ#3°AÑ#5Ñ6€Dð 
€BØ	€DÜ˜#˜q $Ó'�N€Dˆ(Ø	�‰‹€BØ	€BÜ�1�a˜‘c˜!‘eŽ_ˆØˆˆ1‹ñ à€C„HÜ�1�Q�q‘SŽ\ˆØˆÜ�q˜‘s˜1‘u˜Q‘w˜r "Ö%ˆAØ˜1‘H  1¡Ñ%ŠDñ &àˆˆ1‰ä�q˜!˜A™#˜a™% ™'Ö"ˆAØ�q‘S˜B˜q ™s™G‘OˆBˆq‹Eó #ñ ð €CÜ	ˆQ�—‘˜˜1™˜Q™“Ó	 €BØ	ˆC�G‰G�Q”s˜5“z‘\ 1Ñ$Ó%Ñ	% c§g¡g¨d£mÑ	3°AÑ	5€BØ�F‰F�1�c—f‘f˜SŸV™V‘m A‘o aÑ'¨Ñ*¨2Ñ-Ñ.Ó/°Ñ2€EÜ�1�QŽZˆØ�W‰W�S—X‘X˜cŸi™i¨¨#©Ó.Ó/Ó0°·±¸1¸U¹7Ó1CÑCÀBÑFˆÜ�R˜“ˆˆA‹ñ ð
 �1‰v�b‰y€C„HØ��%‘‰[€FØ
€AØ€A€f�I˜1��6‘ q˜a ™h°1¨!¨E©(Ø€A€g�J˜Q�!�G‘*¨!  &¡	°q¨Q¨v©YÜ�1�Q�ZˆØ�q‘6˜"‘9ˆŒÜ�q˜˜1™˜a™ ™Ö"ˆAØ�!‘�A�a‘C‘ˆAØ�1“Ø—W‘W˜Q‘Y�Ø˜‘T�Ø˜R‘i ‘]�Ø˜‘s�V�Ø˜a  ! A¡# a¨¡c 	™l™?¨2¨a°°!°A±#°a°©j©=Ñ8¸¸A¸aÀÀ!ÁÀAÀaÁC¸i¹L¹ÑH��!�A�a�%“à��!�A�a�%‘Ü˜q ž�AØ˜A˜a ˜E™( C§G¡G¨C¯G©G°A°a±C¸¸!¹±GÓ,<Ñ$<¸S¿W¹WÀQÀqÁS»\Ñ$IÑJ�CØ�a˜˜!�e“H "¨¨!©¡¨cÑ 1Ñ1•Hó $ñ #ð ˆˆ!ˆAˆbˆ&‰	˜q�Q�q˜˜2�v‘Y°! ! A a¨¨!©¨Q©¨q© L¤/ñ ô" �B�s˜1‘uŽoˆÜ�r˜!–ˆAÜ˜2œc ! A a¡C¨¡F¨1¡H›oÖ.�Ø˜“T˜a ›c a¨£c¨a°°!±°Q±­hØ !�A�b˜˜1�f“Ió /ó ñ ô �A�c˜!‘eŽnˆÜ�q˜–ˆAØ˜1‘v˜b‘yˆCŒHÜ˜1˜Q˜q™S !™V A™XÖ&�Ø�r‘T˜!‘V˜Q˜r !™t A a¡C¨¨!©˜|™_Ñ,¨Q°°a±¸±	©]¸1¸RÀ¹TÀ!ÀAÁ#ÀaÈÁc¸\¹?Ñ-JÑJ�Ø ! B q¡D¨¨1©¨Q¨q©S L¡/Ñ1��"�Q�q�&“	ó 'ó ñ ð €GØ€FØ€C„HØ	�‰��Q�q‘S‘Ó	€BØ	�‰��Q�q‘S‘˜!‘Ó	€BØ	�‰��A‘”d—i’i §¡Ó'Ñ'Ó	(¨Ñ	*€BÜ�1�QŽZˆØ�!�B‘$‰Y�s—w‘w˜sŸw™w q¨¡u›~Ó.¨rÑ1Ñ1ˆÜ˜˜B˜sŸw™w t›}Ñ,¨QÑ.Ó/°Ñ1°2Ñ5ˆ�‰
Ü˜˜3Ÿ7™7 1 Q¡3›<¨™?¨3¯7©7°4«=Ñ8¸Ñ:Ó;¸AÑ=¸rÑAˆˆq‹	ñ ð €HÛˆÜ�q˜–ˆAÜ˜R  !¡ Q¡ q¡Ö)�Ø#$�˜˜A˜c˜Ó"ó *ó ñ ó
 ˆÜ�q˜–ˆAØ˜q‘zˆCŒHÜ˜Q˜q ™s A™v a™xÖ(�Ø#$ R¨¨# X¡;¨r°!°A±#°a¸±e±)©}Ñ#<¸XÀaÈÁcÈ#ÁgÑ=NÑ#N�˜˜A˜c˜Ñ"Ø#+¨B¨q°¨HÑ#5¸¸#¿%¹%¹À#±~Ñ#F�˜˜A˜c˜Ó"ó )ó ñ õð €EÛˆÜ�q˜–ˆAÜ˜R  !¡ Q¡ q¡Ö)�Ø!"��c˜!˜C�iÓ ó *ó ñ ð �q‰z˜!‰|€C„HÙ	�‹€BØ
�&‰&�‹*ˆ€BãˆÜ�q˜–ˆAØ˜q‘zˆCŒHÜ˜Q˜q ™s A™v a™xÖ(�Ø#$��c˜!˜C�iÑ Ü  3 q¡5ž/�BØ"Ÿ|™|¨C°Ó3°b¸"±fÑ<Ø   R¡¨¨# Ñ.ñ /�Hà˜#˜a ˜)Ó$¨Ñ0Õ$ó *ó )ó ñ ð �q‰z˜A‰~€C„Hð 
€BØ€B€q�EØˆq‰E�!‰G€B€q�Eà�q‰z€C„HÜ�1�QŽZˆØ�1�Q‘3‘˜"˜Q™%‘ˆˆ1‹ñ ð 
€BÛˆÜ�q˜–ˆAØ˜q‘zˆCŒHÜ˜Q˜q ™s A™v a™xÖ(�Ø % c¨!¨C iÑ 0°"°Q±%Ñ 7��3�q˜�9“ó )ó ñ ð €DÛˆÜ�q˜–ˆAØ˜a‘yˆCŒHØˆBÜ˜1˜a ™c 1™f Q™hÖ'�Ø�b˜˜Q˜q˜‘kÑ!’ñ (àˆD��Q�‹Kó ñ ð €EØ€C„HÜ�iŠi˜Ÿ™Ó !Ñ# a¡cÑ*€GØ€C„HØ�g‰g�c˜A˜a™C !™8‘n WÑ,¨tÑ3Ó4€GÜ�'˜3Ÿ7™7 2 q¨¡s¡8Ó,Ó-€GØ€C„HÛˆØˆˆc‰
Ü�q˜˜1™–ˆAØ�#‹J˜$˜s 1 Q¡3˜w™-Ñ'�Jó ñ ð €C„HØ	
ˆS�W‰W�U˜1‘XÓÑ	€BØ��2‘‰;€DØ	ˆC�F‰F�1�a˜Ÿ™‘h‘<Ó Ñ €AØˆs�w‰w˜˜1˜T™6 3§9¡9¨Q¨vÓ#6Ñ6Ñ6¸Ñ9Ó:Ñ:€Dà€C„HØ
ˆQˆs�v‰v‰X‰,€CØˆQ‰3�—‘�s“Ñ
˜Q˜q™SÑ
! #§&¡&¨¡(Ñ
*€CØ�‰�#�‹€AÙ�‹
€AØ�Y‰Y�q˜&Ó!€FØ��1‘‰+˜Ñ	 !Ñ	#€Bð €C„HØ�—‘˜˜3Ÿ9™9 Q q¨¡wÓ/Ñ/°·±¸³Ñ:¸4Ñ?Ó@Ñ@€Eà�r‰z€C„Hðð 
�‰�a˜œC ›F 1™H kÓ	2°1Ñ	5€Bð €C„HÜ��3‘‹L€EÜ��c‘
˜R‘Ó €EÜ��9˜sŸw™w q¨!¨E©'°!°E±'©/Ñ':Ó;Ñ;Ó<€EØ€C„HØ	€BÛˆØ�S‘'˜% ™* R™-Ñ'ˆˆ3‹ñ à€C„HØ€Ir   c                 óò  • U R                  S5      U R                  U-  -   nU R                  nSU l        USU R                  -  -  nUS-   U R	                  SUS-  -  U R                  U5      -  5      -   nUS-   U R	                  SU-  U R                  U5      -  5      -   nX`l        U R                  U5      nXpl        [        X[        US-   5      5      n	US:”  aF  U R                  U R                  SUS-  5      S-  U R                  U R                  5      S-  -
  5      n
US:”  a4  U R                  U R                  U R                  SUS-  5      -  S	-  5      nUS:”  a(  U R                  U R                  SUS-  5      * S
-  5      nUS:”  a5  U R                  U R                  * U R                  SUS-  5      -  S-  5      nU R                  U5      nUS:X  a  SU-  U	S   -  nUS:X  a  SU-  nUW
U	S   -  U	S   -   -  nUS:X  a<  SU-  nU* SU	S   -  W
-  U	S   U
S-  -  -   U	S   -   U R                  U	S   -  W-  -
  -  nUS:X  aZ  SU-  nSU	S   -  W
S-  -  U	S   U
S-  -  -   SU
-  U	S   -  -   nUUSU	S   -  W-  -
  SU	S   -  U
-  U-  -
  U	S   -   U	S   W-  -
  -  nUS	:X  a¶  SU-  nS	U	S   -  W
S-  -  U	S   U
S	-  -  -   SU
S-  -  U	S   -  -   nUSU	S   -  U
-  W-  -
  SU	S   -  U
S-  -  U-  -
  SU	S   -  U-  -
  SU	S   -  U-  U-  -
  nUS	U
-  U	S   -  -   S	U	S   -  W-  -
  S	U	S   -  U
-  U-  -
  U	S	   -   U R                  U	S   -  W-  -   nUU-  nX@l        U R                  W5      $ )z%
z_half(t,der=0) Computes Z^(der)(t)
ú0.5r   r   r   r   rŽ   rz   r   r   r‚   é   ù               @y       €       Àù              @r    ù              (@ù              @)r)   r+   r%   r'   r#   r]   ÚsiegelthetarE  r(   rœ   Úpsir¥   )r	   r“   rª   r©   r«   ÚttÚwpthetaÚwpzÚthetarD  Úps1Úps2Úps3Úps4ÚexpthetaÚzÚzfs                    r
   Úz_halfrZ  r  sÖ  € ð 
‡g�gˆeƒn�S—U‘U˜1‘WÑ€AØ—‘€IØ€C„HØ	
ˆAˆc�f‰f‰H‰€BØ˜‰l˜SŸW™W Q¨¨C©¡[°·±¸³Ñ%;Ó<Ñ<€GØ
�a‰-˜#Ÿ'™' " R¡%¨¯©¨r«
Ñ"2Ó3Ñ
3€CØ„HØ�O‰O˜AÓ€EØ„HÜ	�3œ%  A¡›,Ó	'€BØ
ˆQƒw�c—g‘g˜cŸg™g a¨¨!©›n¨QÑ.°·±¸¿¹³ÀÑ1AÑAÓB�Ø
ˆQƒw�c—g‘g˜cŸe™e C§G¡G¨A¨a°©c£NÑ2°1Ñ4Ó5�Ø
ˆQƒw�c—g‘g˜sŸw™w q¨¨1©›~˜o¨aÑ/Ó0�Ø
ˆQƒw�c—g‘g˜sŸu™u˜f S§W¡W¨Q¨q°©s£^Ñ3°BÑ6Ó7�Ø�x‰x˜‹€HØ
ˆaƒxØˆh‰J�r˜!‘uÑˆØ
ˆaƒxØ�‰[ˆØ��B�q‘E‘	˜"˜Q™%‘Ñ ˆØ
ˆaƒxØ�‰\ˆØˆC��2�a‘5‘˜‘˜R ™U 3¨¡6™\Ñ)¨"¨Q©%Ñ/°·±°b¸±e±¸C±Ñ?Ñ@ˆØ
ˆaƒxØ�‰\ˆØˆb�‰e‰G�C˜‘F‰N˜2˜a™5  a¡™<Ñ'¨¨#©¨b°©e©Ñ3ˆØ��"�R˜‘U‘(˜3‘,‘˜r " Q¡%™x¨™|¨CÑ/Ñ/°°1±Ñ5°b¸±e¸C±iÑ?Ñ@ˆØ
ˆaƒxØˆx‰ZˆØˆb�‰e‰G�C˜‘F‰N˜2˜a™5  a¡™<Ñ'¨¨#¨q©&©°°A±©Ñ6ˆØˆc�"�Q‘%‰i˜‰m˜CÑÑ  2 a¡5¡¨¨a©¡°Ñ 3Ñ3°B°r¸!±u±H¸S±LÑ@ÀÀ2ÀaÁ5ÁÈÁÈSÁÑPˆØ��#‘�b˜‘e‘‰O˜A˜b ™e™G C™KÑ'¨¨"¨Q©%©°©°C©Ñ7¸¸1¹Ñ=¸c¿e¹eÀBÀqÁE¹kÈ#¹oÑMˆØˆq‰DˆØ„HØ�7‰7�1‹:Ðr   c                 ó&	  • U R                   nU R                  U5      nU R                  U5      nSU l         US:”  a  U R                  [	        U5      5      nO*SU R
                  -  US-
  -  [	        SU-
  5      SU-
  -  -  nUS:”  a%  SU R                  USU R
                  -  -  5      -  nOSU-  U-  n[	        SU-
  5      nSU-  [        R                  " U5      -  n	[        SSU R                  U5      -   5      n
SU R                  S	U-  S
U-  U-  -   SU-  U-  U	-  -   5      -   U-   S-   n[        X«5      n
[        SSU R                  SU-  U-  5      -   U-   S-   5      nSU R                  SSU-  -   5      -   U-   S-   nXÀl         U R                  XPR                  X@R                  S5      -
  -  -
  5      nUS:”  aF  U R                  U R                  SUS-  5      5      S-  U R                  U R
                  5      S-  -
  nUS:”  a(  U R                  U R                  SUS-  5      5      * S-  nUS:”  a(  U R                  U R                  SUS-  5      5      * S-  nUS:”  a'  U R                  U R                  SUS-  5      5      S-  nXÐl         [        X[!        US-   5      5      n0 n[!        SUS-   5       H  nU R#                  UU   5      UU'   M     X l         U R%                  SU-  5      nUS:X  a  US   UUS   -  -   nUS:X  a  US   * SUS   -  W-  -
  nUS   UU-  -   nUS:X  a7  SUS   -  W-  SUS   -  US-  -  -   US   -   SUS   -  W-  -   nUS   UU-  -   nUS:X  ac  SUS   -  WS-  -  SUS   -  US-  -  -
  SUS   -  U-  -
  SUS   -  W-  -
  nUSUS   -  U-  U-  -
  US   -
  SUS   -  W-  -   nUS   UU-  -   nUS:X  a²  SUS   -  WS-  -  SUS   -  US-  -  -   SUS   -  US-  -  -   nUSUS   -  U-  W-  -   SUS   -  US-  -  U-  -   nUSUS   -  U-  -   SUS   -  US-  -  -
  SUS   -  U-  -   SUS   -  W-  -
  nUSUS   -  U-  U-  -
  US   -   SUS   -  W-  -
  nUS   UU-  -   nX0l         W$ )z7
zeta_half(s,k=0) Computes zeta^(k)(s) when Re s = 0.5
r   r   r   r   rw   r   r    r   çö(\�Âõ @ç3333335@çÍÌÌÌÌÌô?çš™™™™™@çš™™™™™ñ?rG  r‚   rH  r„   rI  éôÿÿÿrL  rK  é    é   ù              H@rz   )r%   rœ   r›   r,   r    r'   r�   Úlogr"   r#   rM  r+   r)   rN  r]   rE  r(   r¨   r¥   )r	   r©   r?   r«   r$  r“   ÚXÚM1Úabstr	  ÚwpbasicÚwpbasic2rP  ÚwpRrR  rS  rT  rU  rV  r  r  r÷   rW  ÚzvÚzv1s                            r
   Ú	zeta_halfrn  š  s!  € ð —‘€IØ�G‰G�A‹J€EØ�‰�‹
€Aà€C„HàˆqƒyØ�H‰H”S˜“VÓ‰àˆs�v‰v‰X˜˜q™Ñ!¤C¨¨!©£H¨s°5©yÑ$9Ñ9ˆàˆqƒyØˆs�x‰x˜˜1˜SŸV™V™8™Ó%Ñ%‰à�‰U�Q‰Yˆäˆs�1‰u‹:€DØ	ˆ4‰”—’˜“Ñ€Aä�!�A�c—g‘g˜a“j‘LÓ!€GØ�—‘˜˜b™  b¡¨¡Ñ*¨3¨r©6°!©8°A©:Ñ5Ó6Ñ6°yÑ@ÀÑB€HÜ�'Ó$€GÜ�!�Q�s—w‘w˜s 2™v a™xÓ(Ñ(¨Ñ2°1Ñ4Ó5€GØ
ˆC�G‰G�C˜˜!™‘GÓÑ
˜YÑ
& qÑ
(€CØ„HØ�O‰O˜AŸe™e U¯7©7°5«>Ñ%9Ñ:Ñ:Ó;€EØˆ1ƒu�S—W‘W˜SŸW™W Q q¨¡s›^Ó,¨aÑ/°#·&±&¸¿¹³.ÀÑ2BÑBˆcØˆ1ƒu�c—g‘g˜cŸg™g a¨¨!©›nÓ-Ð.¨qÑ0ˆcØˆ1ƒu�c—g‘g˜cŸg™g a¨¨!©›nÓ-Ð.¨qÑ0ˆcØˆ1ƒu�S—W‘W˜SŸW™W Q q¨¡s›^Ó,¨bÑ0ˆcØ„HÜ
�Cœ%  !¡›*Ó
%€CØ
€CÜ�Q�q˜‘sŽ|ˆØ—8‘8˜C ™HÓ%ˆˆC‹ñ à„HØ�x‰x˜˜5™Ó!€HØˆ!ƒtØ�‰V�H˜S ™V‘OÑ#ˆØˆ!ƒtØ�1‰vˆg˜˜#˜a™&™ ™Ñ$ˆØ�‰V�h˜s‘lÑ"ˆØˆ!ƒtØ��A‘‰h�s‰l˜1˜S ™V™8 S¨!¡VÑ,Ñ,¨S°©VÑ3°B°s¸1±v±I¸c±MÑAˆØ�‰V�H˜S‘LÑ ˆØˆ!ƒtØ�#�a‘&‰j˜˜a™Ñ  # a¡&¡¨¨a©¡Ñ/°°#°a±&±¸±Ñ<¸RÀÀAÁ¹YÀs¹]ÑJˆØ�C˜˜A™‘J˜s‘N 3Ñ&Ñ& s¨1¡vÑ-¨a°°A±©h°s©lÑ:ˆØ�‰V�H˜S‘LÑ ˆØˆAƒvØ��Q‘‰i˜˜Q™Ñ  3 q¡6¡	¨#¨q©&Ñ 0Ñ0°°C¸±F±¸3À¹6Ñ1AÑAˆØ�3�s˜1‘v‘:˜c‘> #Ñ%Ñ% c¨#¨a©&¡j°#°q±&Ñ&9¸#Ñ&=Ñ=ˆØ�#�c˜!‘f‘*˜S‘.Ñ   C¨¡F¡¨3°©6Ñ!1Ñ1°!°C¸±F±(¸3±,Ñ>¸qÀÀQÁ¹xÈ¹|ÑKˆØ�"�S˜‘V‘)˜C‘- Ñ#Ñ# C¨¡FÑ*¨2¨c°!©f©9°S©=Ñ8ˆØ�‰V�H˜S‘LÑ ˆØ„HØ€Ir   c                 ó‚
  • U R                   nU R                  U5      nU R                  U5      nSU l         US:”  a  U R                  [	        U5      S5      nODU R                  SU R
                  -  US-
  5      U R                  [	        SU-
  5      SU-
  5      -  nUS:”  a%  SU R                  USU R
                  -  -  5      -  nOSU-  U-  nSU-
  S:”  a%  SU R                  USU R
                  -  -  5      -  nOESU-  U R                  SU R
                  -  U* 5      -  U R                  [	        U5      US-
  5      -  n[	        SU-
  5      n	SU	-  [        R                  " U	5      -  n
[        SSU R                  U5      -   5      nSU R                  S	U-  S
U-  U-  -   SU-  U-  U
-  -   5      -   U-   S-   n[        X¼5      n[        SSU R                  SU-  U-  5      -   U-   S-   5      nSU R                  SSU-  -   5      -   U-   S-   nXÐl         U R                  XPR                  X@R                  S5      -
  -  -
  5      nUnU R                  SU-
  5      nXàl         [        XU5      u  nnUS:”  aP  U R!                  SUS-  5      U R!                  SSU-
  S-  5      -   S-  U R#                  U R
                  5      S-  -
  nUS:”  a>  U R                  U R!                  SUS-  5      U R!                  SSU-
  S-  5      -
  -  S-  nUS:”  a2  U R!                  SUS-  5      U R!                  SSU-
  S-  5      -   * S-  nUS:”  a?  U R                  * U R!                  SUS-  5      U R!                  SSU-
  S-  5      -
  -  S-  nX°l         U R%                  SU-  5      nUS:X  a  US   UUS   -  -   nUS:X  a  US   * SUS   -  W-  -
  nUS   UU-  -   nUS:X  a7  SUS   -  W-  SUS   -  US-  -  -   US   -   SUS   -  W-  -   nUS   UU-  -   nUS:X  ac  SUS   -  WS-  -  SUS   -  US-  -  -
  SUS   -  U-  -
  SUS   -  W-  -
  nUSUS   -  U-  U-  -
  US   -
  SUS   -  W-  -   nUS   UU-  -   nUS:X  a²  SUS   -  WS-  -  SUS   -  US-  -  -   SUS   -  US-  -  -   nUSUS   -  U-  W-  -   SUS   -  US-  -  U-  -   nUSUS   -  U-  -   SUS   -  US-  -  -
  SUS   -  U-  -   SUS   -  W-  -
  nUSUS   -  U-  U-  -
  US   -   SUS   -  W-  -
  nUS   UU-  -   nX0l         W$ )z#
Computes zeta^(k)(s) off the line
r   r   rw   r   r   r   r    r   r\  r]  r^  r_  r`  rG  r‚   rH  rb  r„   rI  ra  rL  rK  rc  rd  rz   )r%   rœ   r›   rl   r    r'   r,   r�   re  r"   r#   rM  r+   r)   r¨   r  rN  r]   r¥   )r	   r©   r?   r«   r$  r“   rf  rg  ÚM2rh  r	  ri  rj  rP  rk  rR  Ús1Ús2r  r  rS  rT  rU  rV  rW  rl  rm  s                              r
   Úzeta_offliners  Ø  sÄ  € ð —‘€IØ�G‰G�A‹J€EØ�‰�‹
€Aà€C„HàˆqƒyØ�I‰I”c˜!“f˜cÓ"‰à�I‰I�a˜Ÿ™‘h  a¡Ó(¨¯©´3°q¸±s³8¸CÀ¹IÓ)FÑFˆà�‹	Øˆs�x‰x˜˜1˜SŸV™V™8™Ó%Ñ%‰à�‰U�Q‰Yˆà	ˆ%‰�!‹Øˆs�x‰x˜˜1˜SŸV™V™8™Ó%Ñ%‰àˆq‰S�—‘˜1˜SŸV™V™8 e VÓ,Ñ,¨S¯Y©Y´s¸1³v¸eÀC¹iÓ-HÑHˆäˆs�1‰u‹:€DØ	ˆ4‰”—’˜“Ñ€Aä�!�A�c—g‘g˜a“j‘LÓ!€GØ�—‘˜˜b™  b¡¨¡Ñ*¨3¨r©6°!©8°A©:Ñ5Ó6Ñ6°yÑ@ÀÑB€HÜ�'Ó$€GÜ�!�Q�s—w‘w˜s 2™v a™xÓ(Ñ(¨Ñ2°1Ñ4Ó5€GØ
ˆC�G‰G�C˜˜!™‘GÓÑ
˜YÑ
& qÑ
(€CØ„HØ�O‰O˜AŸe™e U¯7©7°5«>Ñ%9Ñ:Ñ:Ó;€EØ	
€BØ	�‰�!�B‘$‹€BØ„HÜ˜3 1Ó%�H€CˆØˆ1ƒu�S—W‘W˜Q˜r !™t“_ S§W¡W¨Q°°"±°a©xÓ%8Ñ8¸!Ñ;¸c¿f¹fÀSÇVÁV»nÈQÑ>NÑNˆcØˆ1ƒu�C—E‘E˜3Ÿ7™7 1 R¨¡T›?¨3¯7©7°1°a¸±d¸A±XÓ+>Ñ>Ñ?ÀÑAˆcØˆ1ƒu�c—g‘g˜a  1¡“o c§g¡g¨a°°2±°q±Ó&9Ñ9Ð:¸2Ñ=ˆcØˆ1ƒu�S—U‘U�F˜CŸG™G A b¨¡d›O¨C¯G©G°A°q¸±t¸Q±hÓ,?Ñ?Ñ@ÀÑCˆcØ„HØ�x‰x˜˜5™Ó!€HØˆAƒvØ�‰V�H˜S ™V‘OÑ#ˆØˆAƒvØ�1‰vˆg�a˜˜A™‘h˜s‘lÑ"ˆØ�‰V�H˜S‘LÑ ˆØˆAƒvØ��A‘‰h�s‰l˜1˜S ™V™8 S¨!¡VÑ,Ñ,¨c°!©fÑ4°R¸¸A¹±Y¸s±]ÑBˆØ�‰V�H˜S‘LÑ ˆØˆAƒvØ�#�a‘&‰j˜˜a™Ñ  3 q¡6¡¨#¨q©&¡Ñ0°°3°q±6±¸#±Ñ=¸bÀÀQÁ¹iÈ¹mÑKˆØ�C˜˜A™‘J˜s‘N 3Ñ&Ñ& s¨1¡vÑ-¨a°°A±©h°s©lÑ:ˆØ�‰V�H˜S‘LÑ ˆØˆAƒvØ��Q‘‰i˜˜Q™Ñ  3 q¡6¡	¨#¨q©&Ñ 0Ñ0°°C¸±F±¸3À¹6Ñ1AÑAˆØ�3�s˜1‘v‘:˜c‘> #Ñ%Ñ% c¨#¨a©&¡j°#°q±&Ñ&9¸#Ñ&=Ñ=ˆØ�#�c˜!‘f‘*˜S‘.Ñ   C¨¡F¡¨3°©6Ñ!1Ñ1°!°C¸±F±(¸3±,Ñ>¸qÀÀQÁ¹xÈ¹|ÑKˆØ�"�S˜‘V‘)˜C‘- Ñ#Ñ# C¨¡FÑ*¨2¨c°!©f©9°S©=Ñ8ˆØ�‰V�H˜S‘LÑ ˆØ„HØ€Ir   c           
      óÀ  • U R                  S5      U R                  U-  -   nUnU R                  SU-
  5      nU R                  nSU l        U R	                  U5      S:¼  aN  SU R                  U R                  U5      SU R                  -  -  5      -  nU R                  [        U5      5      nO_SU R                  -  U R	                  U5      S-
  -  [        SU-
  5      SU R	                  U5      -
  -  -  nSU R                  U5      -  U-  nU R	                  U5      S:¼  a4  SU R                  U R                  U5      SU R                  -  -  5      -  n	O]SU R                  U5      -  SU R                  -  U R	                  U5      S-
  -  -  [        SU-
  5      SU R	                  U5      -
  -  -  n	S[        U R                  U5      5      -  n
U R                  U5      nX·U	-   -  nSU-   n[        S	SU R                  U
5      -   U R                  US
SU
-  -   -  5      U-   5      n[        SU R                  SU-  5      U-   5      nU R                  SU-  5      U-   nXðl        U R                  U5      nUU l        [        XU5      u  nnSSU-  -   nSSU-  -
  nUS:”  aG  SU R                  SU5      U R                  SU5      -   -  U R                  U R                  5      S-  -
  nUS:”  a(  SU R                  SU5      U R                  SU5      -
  -  nUS:”  a(  SU R                  SU5      U R                  SU5      -   -  nUS:”  a(  SU R                  SU5      U R                  SU5      -
  -  nXàl        U R                  U5      nUS:X  a  UUS   -  US   U-  -   nU R                  nUS:X  a)  UU-  US   US   W-  -   -  UUS   US   U-  -   -  U-  -
  nUS:X  a[  USUS   -  W-  US   US-  -  -
  US   -
  UUS   -  W-  -   -  nUSUS   -  U-  US   US-  -  -
  US   -
  UUS   -  U-  -
  U-  -   nUS:X  a»  SUS   -  WS-  -  US   US-  -  -
  SUS   -  U-  -
  US-  US   -  W-  -   nUSUS   -  U-  U-  -   US   -
  US   W-  -   U-  U-  nSUS   -  US-  -  US   US-  -  -   SUS   -  U-  -   US-  US   -  U-  -   nUUSUS   -  U-  U-  -   US   -   US   U-  -
  -  U-  nUU-   nUS:X  GaS  SUS   -  WS-  -  US   US-  -  -   S	US   -  US-  -  -   nUSUS   -  U-  W-  -
  SUS   -  US-  -  U-  -
  SUS   -  U-  -
  nUSUS   -  U-  U-  -
  SUS   -  U-  -   SUS   -  W-  -
  SUS   -  U-  U-  -
  nUUS   -   UUS   -  W-  -   nSUS   -  US-  -  US   US-  -  -   S	US   -  US-  -  -   nUSUS   -  U-  U-  -   SUS   -  US-  -  U-  -   SUS   -  U-  -   nUSUS   -  U-  U-  -
  SUS   -  U-  -   SUS   -  U-  -
  SUS   -  U-  U-  -
  nUUS   -   UUS   -  U-  -
  nUU-  UU-  -   nX`l        W$ )z0
Computes Z(w) and its derivatives off the line
rG  r   é#   r   r   rw   r   r   r    é   gR¸…ëQ @g      Ð?y              à?y              À?g      °¿y       €       ¿r„   rŒ   rJ  rK  rL  )r)   r+   r¨   r%   rœ   r,   r›   r'   r    rM  r"   r#   r  rN  r]   r¥   )r	   r8   r?   r©   rq  rr  r«   rg  rf  rp  r	  rd   re   rg   ri  rP  rk  rR  r  r  ÚptaÚptbrS  rT  rU  rV  rW  rl  r+   rm  Úzv2s                                  r
   Ú	z_offlinerz    sA  € ð 	�‰�‹�s—u‘u˜Q‘wÑ€AØ	
€BØ	�‰�!�B‘$‹€BØ—‘€IØ€C„Hð 	�‰�‹�qÓØˆs�x‰x˜Ÿ™ › Q¨¯©¡ZÑ0Ó1Ñ1ˆØ�H‰H”S˜“WÓ‰àˆs�v‰v‰X˜Ÿ™ › Q™Ñ'¬#¨a°©d«)°c¸#¿'¹'À"»+±oÑ*FÑFˆØ�—‘˜“‰_˜QÑˆà�‰�‹�qÓØˆs�x‰x˜Ÿ™ › Q¨¯©¡ZÑ0Ó1Ñ1‰à�—‘˜“‰_˜a §¡™h¨#¯'©'°"«+°a©-Ñ8Ñ8¼¸Q¸r¹T»ÀSÈÏÉÐQSËÁ_Ñ9UÑUˆà	Œ#ˆc�o‰o˜aÓ Ó
!Ñ!€Að �8‰8�A‹;€DØ�B‘‰<€DØˆi‰<€DÜ�!�Q�s—w‘w˜q“z‘\ 3§7¡7¨4°°A°a±C±©=Ó#9¸$Ñ#>Ó?€GÜ�!�C—G‘G˜D ™IÓ& tÑ+Ó,€GØ
�'‰'�!�D‘&‹/˜$Ñ
€Cà„HØ�O‰O˜AÓ€EØ€C„HÜ˜3 Ó#�H€CˆØ
��a‘‰-€CØ
��a‘‰-€CØˆ1ƒu�D˜#Ÿ'™' ! C›.¨¯©°°3«Ñ7Ñ8¸3¿6¹6À#Ç&Á&»>È!Ñ;KÑKˆcØˆ1ƒu�T˜CŸG™G A c›N¨3¯7©7°1°S«>Ñ9Ñ:ˆcØˆ1ƒu�V˜cŸg™g a¨›n¨S¯W©W°Q°s«^Ñ;Ñ<ˆcØˆ1ƒu�V˜cŸg™g a¨›n¨S¯W©W°Q°s«^Ñ;Ñ<ˆcØ„HØ�x‰x˜‹€HØˆAƒvØ�c˜!‘f‰_˜S ™V H™_Ñ,ˆØ�‰€AØˆAƒvØˆx‰Z˜˜Q™  A¡ s¡
Ñ*Ñ+¨A¨s°1©v°c¸!±f¸S±jÑ/@Ñ,AÀ(Ñ,JÑJˆØˆAƒvØ�r˜#˜a™&‘y ‘} S¨¡V¨C°©F¡]Ñ2°3°q±6Ñ9¸!¸CÀ¹F¹(À3¹,ÑFÑGˆØ�"�S˜‘V‘)˜C‘-  A¡ s¨A¡v¡Ñ-¨c°!©fÑ4°Q°s¸1±v±X¸c±\ÑAÀ8ÑKÑKˆØˆAƒvØ��Q‘‰i˜˜Q™Ñ˜s 1™v c¨1¡f™}Ñ,¨Q¨s°1©v©X°c©\Ñ9¸!¸A¹#¸cÀ!¹f¹*ÀS¹.ÑHˆØ�B�s˜1‘v‘I˜c‘M #Ñ%Ñ% c¨!¡fÑ,¨S°©V°C©ZÑ7¸Ñ:¸8ÑCˆØ��A‘‰h�s˜A‘v‰o˜c !™f S¨!¡V™mÑ+¨A¨c°!©f©H°S©LÑ8¸¸1¹¸SÀ¹V¹ÀC¹ÑGˆØ��r˜#˜a™&‘y ‘} SÑ(Ñ(¨#¨a©&Ñ0°°Q±¸±Ñ;Ñ<¸XÑEˆØ�‰WˆØˆA„vØ��A‘‰h�s˜A‘v‰o˜c !™f S¨!¡V™mÑ+¨a°°A±©h°s¸A±v©oÑ=ˆØ�#�c˜!‘f‘*˜S‘. Ñ$Ñ$ R¨¨A©¡Y¨s°A©vÑ%5°cÑ%9Ñ9¸"¸SÀ¹V¹)ÀC¹-ÑGˆØ�!�C˜‘F‘(˜3‘,˜sÑ"Ñ" 1 S¨¡V¡8¨C¡<Ñ/°°#°a±&±¸±Ñ<¸Q¸sÀ1¹v¹XÀc¹\È#Ñ=MÑMˆØ�#�a‘&‰j˜˜3˜q™6™ #™Ñ%ˆØ��A‘‰h�s˜A‘v‰o˜c !™f S¨!¡V™mÑ+¨a°°A±©h°s¸A±v©oÑ=ˆØ�#�c˜!‘f‘*˜S‘. Ñ$Ñ$ R¨¨A©¡Y¨s°A©vÑ%5°cÑ%9Ñ9¸"¸SÀ¹V¹)ÀC¹-ÑGˆØ�!�C˜‘F‘(˜3‘,˜sÑ"Ñ" 1 S¨¡V¡8¨C¡<Ñ/°°#°a±&±¸±Ñ<¸Q¸sÀ1¹v¹XÀc¹\È#Ñ=MÑMˆØ�#�a‘&‰j˜˜3˜q™6™ #™Ñ%ˆØ�c‰\˜#˜h™,Ñ&ˆØ„HØ€Ir   c                 ó6  • US:”  a  [         eU R                  U5      nU R                  U5      o@R                  U5      nUS:  a2  U R	                  U R                  U R	                  U5      U5      5      nU$ US:H  nU(       a  [        XU5      $ [        XU5      $ )Nr   r   rw   )r¡   rL   rœ   r›   r¨   Úrs_zetarn  rs  )r	   r©   Ú
derivativeÚkwargsÚreÚimrX  Úcritical_lines           r
   r|  r|  b  s†   € à�Aƒ~Ü!Ð!Ø�‰�A‹€AØ	�‰�‹€BŸ'™' !›*�RØ	ˆAƒvØ�H‰H�S—[‘[ §¡¨!£¨jÓ9Ó:ˆØˆØ˜3‘Y€MÞÜ˜ Ó,Ð,ä˜C JÓ/Ð/r   c                 óÔ   • U R                  U5      nU R                  U5      o0R                  U5      nUS:  a  [        X* U5      $ US:H  nU(       a  [	        XU5      $ [        XU5      $ )Nr   )rL   rœ   r›   Úrs_zrZ  rz  )r	   r8   r}  r  r€  r�  s         r
   rƒ  rƒ  q  sb   € à�‰�A‹€AØ	�‰�‹€BŸ'™' !›*�RØ	ˆAƒvÜ�C˜˜ZÓ(Ð(Ø˜1‘W€MÞÜ�c˜jÓ)Ð)ä˜ Ó,Ð,r   )r   )Ú__doc__r�   Úobjectr   Ú	functionsr   rG   rZ   ri   rq   r  rE  rZ  rn  rs  rz  r|  rƒ  r   r   r
   Ú<module>r‡     sŽ   ðñ0ód ô(ˆfô (õ ò*t'òl.ò6ò	ôDðL $% #ô qôh&ôP<ô|@ôDFðP ó0ó ð0ð ó	-ó ñ	-r   