403Webshell
Server IP : 35.80.110.71  /  Your IP : 216.73.216.221
Web Server : Apache/2.4.58 (Ubuntu)
System : Linux ip-172-31-21-44 6.17.0-1019-aws #19~24.04.1-Ubuntu SMP Tue Jun 23 18:53:06 UTC 2026 x86_64
User : ubuntu ( 1000)
PHP Version : 8.3.31
Disable Function : NONE
MySQL : OFF  |  cURL : ON  |  WGET : ON  |  Perl : ON  |  Python : OFF  |  Sudo : ON  |  Pkexec : OFF
Directory :  /lib/python3/dist-packages/mpmath/functions/__pycache__/

Upload File :
current_dir [ Writeable ] document_root [ Writeable ]

 

Command :


[ Back ]     

Current File : /lib/python3/dist-packages/mpmath/functions/__pycache__/zeta.cpython-312.pyc
�

�"`%����ddlmZddlmZddlmZmZmZed+d��Zed,d��Z	ed��Z
ed	��Zgd
�Zd�Z
ed-d��Zed
��Zed��Zed��Zed��Zed��Zed��Zed.d��Zd�Zd�Zd�Zd�Zed��Zed.d��Zed.d��Zed��Zed��Zed/d��Zed0d ��Z d!�Z!d"�Z"edgdfd#��Z#edgdfd$��Z$d%�Z%d&�Z&d'�Z'd(�Z(ed1d)��Z)ed*��Z*y)2�)�print_function�)�xrange�)�defun�
defun_wrapped�defun_staticc�"����	��j����j����dkr�jd�St�d�r
�j}n	ix}�_�dk(r0�dk(r
�j�S�|vr|�\}}|�j
k\r|�Sd�	���	�fd�}�j
}	�dkDr'd�_�j
|d�jgd�	��	|d
zt�dz�z�_�j
|d�jgd�	�}�j���zd�zz�j���dzz�dzzz
d|z�z�	zz}|�_�dk(r"�j��r�j
|f|�<|�S#|�_wxYw)
Nrz&Stieltjes constants defined for n >= 0�stieltjes_cacherc���|�z}|�jz
�j��j|zz
��zzd|dzzz�jd�jz|z�dz
z}�j	|��zS�Nrr)�j�ln�exp�pi�_re)�x�xa�v�a�ctx�mag�ns   �����7/usr/lib/python3/dist-packages/mpmath/functions/zeta.py�fzstieltjes.<locals>.fsw���
�q�S��
����X�s�v�v�a����a��i�(�!�+�+�Q�r�1�u�W�5�s�w�w�q����x��z�7J�1�7L�M���w�w�q�z�C����2��)�	maxdegree�
��?r)�convert�
bad_domain�hasattrr�euler�prec�quad�inf�intr�isint)
rrrrr'�sr�origrrs
```      @r�	stieltjesr.s�������A��A����A��A��1�u��~�~�F�G�G��s�%�&��-�-��02�2��#�-��A�v���6��I�I�:�����%�a�(�G�D�!��s�x�x���r�	�
�C� ��8�8�D�
�
�r�6��C�H��(�(�1�q����k�Q�(�7�C��"�9�s�1�c�6�{�*����H�H�Q��3�7�7��r�H�2���F�F�1�I�q�L�!�A�#�������Q�q�S�!1�1�Q�3�!7�7�!�A�#�a�%��)�C������A�v�#�)�)�A�,�!�h�h��]����
�2�I�����s
�+B)F�	Fc�,�t|�}||jk(s||jk(r=|dkr,||jk(r|dk(r|jS|jS|jS|dk(r�|j	|�rY|jdd|zz�}|jdd|zz
�}|j
|j�dz|zd||z
zz
S|j|�r|S|j	|jdd|zz��|j
|j�dz|zz
S|dkDr�d|dz
z|j|dz
dd|zz
�z}d|dz
z|j|dz
dd|zz�z}|j	|�r4|dk(r'd|j|j�zd||zzzSd||zzS|dk(r6|jd|j|j�zd||zzz�S|jd||zz�Sy)Nrr��?y�?���r��)r*r)�ninf�zero�_im�loggammarr�isinf�	polygamma�logr)r�t�
derivative�dr�bs      r�siegelthetar>,s����J��A�	
�c�g�g���c�h�h���q�5��C�H�H�}��a���x�x���7�7�N��8�8�O��A�v��7�7�1�:����T�$�q�&�[�)�A����T�$�q�&�[�)�A��F�F�3�6�6�N�?�1�$�Q�&��q��s��3�3��y�y��|����7�7�3�<�<��T�!�V��4�5����s�v�v��q�8H��8J�J�J��1�u�
�a��c�N�3�=�=��1��d�4��6�k�:�:��
�Q�q�S�M�#�-�-��!��T�$�q�&�[�9�9���7�7�1�:��A�v��C�G�G�C�F�F�O�+�D�!�A�#�J�6�6��Q�q�S�z�!��A�v��w�w�t�C�G�G�C�F�F�O�3�D�!�A�#�J�>�?�?��w�w�t�Q�q�S�z�*�*�
rc	�����d�jz�jd�jd�zdzd�jzz�z�z}�j	��fd�|�S)Nrr�c�F���j|��j�zz
S�N)r>r)r:rrs ��r�<lambda>zgrampoint.<locals>.<lambda>Ss���#�/�/�!�"4�S�V�V�A�X�"=�r)rr�lambertw�e�findroot)rr�gs`` r�	grampointrHNsU���	
�#�&�&������3�<�<��1��Q���3�5�5��(9�:�:�;�;�A��<�<�=�q�A�Arc
�����
����������t|jdd��}�j|�}�j|�}�j	|�}�j
}	t
|�d|zkDr>|dz|kr6�j||�}�j|�r�j|�S|S�xj
dz
c_�j�j|��}�jd�j|zz��|dk(r1|�z}|�_�j|�r�j|�S|�S�jd�j|zzd����j|d���|dk(rD�j|z���zzz}|�_�j|�r�j|�S|�S�jd�j|zzd����j|d����dz�j�zz
�|dk(rM�����fd	�}	�j|	d�}||z}|�_�j|�r�j|�S|�S�xj
d
z
c_�jd�j|zzd����j|d����dzd�jz�z�zz
�z
�|dk(r\�������fd�}	�j|	d�}�j|z|z}|�_�j|�r�j|�S|�S�jd�j|zzd
����j|d
��������fd�}	�j|	d��
|d
k(rQ��
��������f
d�}	�j|	d�}||z}|�_�j|�r�j|�S|�S|d
kDr�fd�}
�j|
||d
z
��Sy#t$rY��GwxYw)Nr;r��r�r"r�r;c� ��d�z�z���zgS�Nr�)�comb1�theta1�z�z1�z2s�����r�termszsiegelz.<locals>.termszs����b�D��K��Q�u�W�-�-rr!rc�2��d�z�zd�z�z���zzgS)NrrO)rP�comb2rQrRrSrT�z3s�������rrUzsiegelz.<locals>.terms�s(����v�X�b�[�!�B�$�u�*�b��5��j�9�9r�c�v���dzd�jz�dzz�zd�dzzd�z�z�j�zgS)NrY���r�������r)rrQ�theta2�theta3�theta4s�����rrUzsiegelz.<locals>.terms�sL�����	�2�c�e�e�8�F�A�I�-�f�4�b����l��v�I�f��c�e�e�F�l�,�	,rc�j�
�d�dzz�zd�jz�z�zd�z�zd�z�z�	��zgS)N�rr[rYr^)
rW�comb3rrQr_rRrSrTrX�z4s
����������rrUzsiegelz.<locals>.terms�sL����v�q�y�[��^�R����X�b�[��%7��6��"���2��e��R��5��*�
*rc�*���j|d��S)NrYrL)�siegelz)rrs �rrCzsiegelz.<locals>.<lambda>�s���c�k�k�!��k�2�r)r)r*�getr#rr5r'�abs�rs_z�
_is_real_type�NotImplementedError�expjr>�zetar�sum_accurately�diff)rr:�kwargsr<�t1�t2r'r�e1rU�hrPrWrdrQr_r`rarRrSrTrXres`          @@@@@@@@@@@@rrgrgVs������F�J�J�|�Q�'�(�A����A��A�	�����B�	�����B��8�8�D�
��r�7�S��X��"�a�%�"�*�����A��A�� � ��#��w�w�q�z�!��H��H�H��N�H�	���#�/�/�!�$�	%�B�����S�U�U�1�W���A��A�v��q�D��������Q���7�7�1�:���r�	�	���#�c�e�e�A�g�+�!��	,�B�
�_�_�Q�1�_�
-�F��A�v�
�U�U�2�X�r�!�F�(�{�
#��������Q���7�7�1�:���r�	�	���#�c�e�e�A�g�+�!��	,�B�
�_�_�Q�1�_�
-�F��A�I�c�e�e�F�l�"�E��A�v�	.�	.����u�a�(���S��U��������Q���7�7�1�:���r�	��H�H��N�H�	���#�c�e�e�A�g�+�!��	,�B�
�_�_�Q�1�_�
-�F��A�I�a����g�f�n�V�+�+�F�2�E��A�v�	:�	:����u�a�(���e�e�V�B�Y�q�[��������Q���7�7�1�:���r�	�	���#�c�e�e�A�g�+�!��	,�B�
�_�_�Q�1�_�
-�F�,�,�
���u�a�(�E��A�v�	*�	*�
���u�a�(��
��T��������Q���7�7�1�:���r�	��1�u�2���x�x��1��!��x�$�$�
��y�
��
�s�)AO3�6O3�3	P�?P)dg�����D,@g|��c�5@g����9@gl���l>@gY�w@@gjD���B@gp`���uD@g,���ݩE@ge��ȨH@g��P��H@g�N�~3|J@g�C�9L@g7��k�M@g�1��wjN@g?a��3GP@g�����P@gM��>�bQ@gY`�OLR@g���R@gX�GEIS@g�Xe���S@gi��C�T@gH�&Q/U@gܩ�7�U@g�K`z�3V@g�wG{W@g��ꞯ�W@gr�x��W@g˻;I2�X@gS�WTY@gts7o�Y@g�~�x�\Z@gEA����Z@gK���[@g|j��[@g�Z�~�\@g�U-�]@gB�/��]@gg�� �W^@g�'�٘�^@g��o_@g�CY�_@gڡ���2`@g#�-X�b`@gg�v�`@g��S5�`@g��ӝ�Ca@g�{�wa@g�:i��a@g��=��a@g�l@b@g'B/K�mb@g�;up��b@gH�ٵ��b@gլJ� c@g����c@gZU�x�c@gҍd3�c@gu��&d@g���`d@g�ĸ�/�d@g�2����d@g�1E#e@g�L�.=e@gOТN+�e@g��BV"�e@g�'�: f@g��x�Lf@g��K�S}f@g��c��f@g���g@goLc<)3g@gstUSgg@g�.�+Q�g@g6�p^�h@g����"h@g%�-!~hh@g@T�#�h@go��j}�h@g
��x(i@g�㸆�Oi@gxۙ��i@g5�\��i@gk���i@g늁�r2j@g�qߋvj@g�'"�j@g֍d��j@g�h�lk@g�ѿ)bk@g�y��k@g�Wȭk@g=I�X9l@g`Lewl@g���x|ml@g+P�̪l@g�����l@g��bn��l@g��]06m@g�J}Ɛm@c��ddl}|j|�}|j�D�cgc]
}t|���}}t	|d�dk(sJ�|t
ddycc}w)Nr�)�urllib�urlopen�	readlines�float�round�_zeta_zeros)�urlrxr<r�Ls     r�_load_zeta_zerosr��sT������s��A��;�;�=�)�a��q��)�A�)���1��;�"�����K��N��	*s�Ac	�l�t|�}|dkr |j|�j�S|dk(rtd��|t	t
�kDr|dkrt
|�|t	t
�kDrtd��|jd|j|jt
|dz
��S)Nrzn must be nonzeroi��zn too large for zetazerosr"r)r*�zetazero�	conjugate�
ValueError�lenr}r�rl�mpcrFrg)rrr~s   r�oldzetazeror��s����A��A��1�u��|�|�Q�B��)�)�+�+��A�v��,�-�-��3�{����V������3�{���!�"=�>�>��7�7�3����S�[�[�+�a��c�2B�C�D�Drc	��|dk(r|jSt|�dkDrZ|j|�}d|j|j|��z}t|�t|�|jzkr|St|�dkr8|xj
t
|jt|�d��z
c_|jx}}|j|�}d}t|�t|�|jzkDrO||z|z}||||j|dz�zzz
}|dz
}t|�t|�|jzkDr�O|S)Nr��r"g{�G�z�?rr)r4ri�li�sqrt�epsr'r*r9�oner�	_zeta_int)rrrr=r,r:�u�ks        r�riemannrr��s!���A�v��x�x��
�1�v��}��F�F�1�I������s�x�x��{�#�#���q�6�C��F�3�7�7�N�"��H�
�1�v��}����C�����Q���*�*�+�+���G�G�O�A�����q�	�A�	�A�

�a�&�3�q�6�#�'�'�>�
!�
��E�A�I��	�Q�!�c�m�m�A�a�C�(�(�
)�)��	�Q����a�&�3�q�6�#�'�'�>�
!�
�Hrc�X�t|�}|dkryt|j|��S)Nrr)r*r��list_primes)rrs  r�primepir��s)���A��A��1�u���s���q�!�"�"rc�f�t|�}|dkr|jjS|dkr*|jj|j	|��S|j|�}|j
|d��|j|d��zdz|jd��z}|j|jj|�|z
jd��}|j|jj|�|zjd��}|jj||g�S)Nria
r�)�roundingr@r<)
r*�_ivr4�mpfr�r�r�rr�floorr�ceilr=)rr�mid�errrr=s      r�primepi2r��s����A��A��1�u��w�w�|�|���4�x��w�w�{�{�3�;�;�q�>�*�*�

�&�&��)�C�

�(�(�1�c�(�
"�3�6�6�!�S�6�#9�
9�!�
;�C�F�F�C�F�<P�
P�C��	�	�3�7�7�;�;�s�#�C�'�*�*�S�	�9�A����#�'�'�+�+�c�"�3�&�)�)�C��8�A��7�7�;�;��!�u��rc������j��r�S�j��dkrtd���dk(r�jS�dk(r&�j	�j
�j�S�j��}|�jkDrd�zS�jt|�z����fd�}�j|�S)Nrz.prime zeta function defined only for re(s) > 0rr"c3���K��j}d}	|dz
}�j|�}|s���_|�j�j|�z��z|z}|sy|�_|���X�w)Nrr)r'�moebiusrrn)r-r�r�r:rr,�wps    ���rrUzprimezeta.<locals>.termssx������8�8�D��A���Q����K�K��N��������c�f�f�S�X�X�a��c�]�+�+�A�-����������s�A(A+)
�isnan�rer�r)r�r3rr'r*ro)rr,�rrUr�s``  @r�	primezetar��s����
�y�y��|���
�v�v�a�y�A�~��I�J�J��A�v��w�w���C�x��w�w�s�x�x����(�(����q�	�A��3�8�8�|��A�v�
�
�X�X��A��
��	�"���e�$�$rc�:����t����dkrtd���dk(s
�dk(r�dkDr�j��S�dk(r*�jdd�z
�dz
�j��zS�dkr=�dk(r�dzS�dk(r�dz
S�dk(rd�z�dz
zdzdzS�dk(r���dz
zdzzS�j	��r��zS�j��r�St
��dkDr���fd	�}�j|���zzS���fd
�}�j|�S)Nrz-Bernoulli polynomials only defined for n >= 0rr"rrrcg�?c3���K��j}|���j�z}d}|�kr=|�dz|z
z|z|z}|dkDr|dzs|�j|�z��|dz
}|�kr�<yy�wr
)r��	bernoulli)r:r�r�rrrRs   ���rrUzbernpoly.<locals>.terms6sz��������A��G�����	�A��A��q�&��q��s�1�u�I�a�K��M���A��!�a�%��C�M�M�!�,�,�,��Q���	�q�&�s�A"A(�&A(c3���K��j�����j}d}|�krB|�dz|z
z|z�z}�|z
}|dkDr|dzs|�j|�z��|dz
}|�kr�Ayy�wr
)r�r�)r:r��mrrrRs   ���rrUzbernpoly.<locals>.termsBs�������-�-��"�"����A��A��q�&��q��s�1�u�I�a�K�!�O���a�C���A��!�a�%��C�M�M�!�,�,�,��Q����q�&�s�A'A-�+A-)r*r�r��ldexpr7r�riro)rrrRrUs``` r�bernpolyr�!s<���	�A��A��1�u��H�I�I��A�v�!�q�&�Q��U��}�}�Q����C�x��	�	�!�A�a�C� ��"�C�M�M�!�$4�4�4��A�v���6�!�q�&�=���6�!�c�'�>���6�1�Q�3��!��9�Q�;��/�)���6�!�Q��#��Y�s�]�+�+�
�y�y��|��A�v�
�
�y�y��|���
�1�v��z�		��!�!�%�(�1�a�4�/�/�		��!�!�%�(�(rc������t����dkrtd���dkr!�dk(r�dzS�dk(r�dz
S�dk(r��dz
zS�j��r��zS�j��r�S�dz}�dk(r3d�j	d|�dz
z�j|�z|z�dzzS�dk(r3d�j	d|�dz
z�j|�z|z�dzzS�dk(r_�dzr�jS�dks'��jd�z�z�jd	zkr"�j	�j����S���fd
�}�j|�|zS)Nrz)Euler polynomials only defined for n >= 0rrr"����dg.e�T>��?r0c3���K��j}d}�jd�dz�}	�|z
dz}|dkDr|dzsd|z
�j|�z|z��|dz
}|�kDry|�z�|z
dzz|z}|dz}�P�w)Nrrrr")r�r�r�)r:r��wrrrrRs    ���rrUzeulerpoly.<locals>.termsgs�������G�G��
���I�I�a��!������!��A��A���E�a�!�e���s�C�M�M�!�,�,�Q�.�.�
��F�A��1�u���!��Q�q�S��U��A�
�A�
��H�A��s�A5A8)r*r�r7r�r�r�r4rr'�	_eulernumro)rrrRr�rUs```  r�	eulerpolyr�Nsv����A��A��1�u��D�E�E��A�v���6�!�q�&�=���6�!�c�'�>���6�!�Q�q�S�'�>�
�y�y��|��!�t��
�y�y��|���	�!��A��A�v��3�9�9�Q�q�>�!�#�$�S�]�]�1�%5�5�a�7�!�Q�$�>�>��A�v��#�)�)�A�a�.��"�#�C�M�M�!�$4�4�Q�6��A��=�=��C�x��q�5��8�8�O��s�7�a����
�1��-�-�����
�=��9�9�S�]�]�1�-��r�2�2�����e�$�q�(�(rFc��t|�}|rt|j|��S|dkr |j|j|��S|dzr|jS|j	|j|d�|�S)Nr�rr")r*r�r�r4r�r�)rr�exacts   r�eulernumr�vsm���A��A���3�=�=��#�$�$��3�w��w�w�s�}�}�Q�'�(�(��1�u��x�x���9�9�S�]�]�1�S�)�1�-�-rc��|j�}|j}d}|}	|||zz}||z
}t|�|kr	|S||z}|dz
}�))Nr)r�r4ri)rr,rR�tol�lr��zk�terms        r�polylog_seriesr��sg���7�7�(�C����A�	�A�	
�B�
��A�q�D�y��	�T�	���t�9�s�?��
�H�	�a���	�Q���
rc���|dkr|dzSd|jz}||z|j|�z|j||j|�|z�z}|j	|�r|dkr|j|�}|j
|�dks(|j
|�dk(rF|j|�dk\r2|||j|�|dz
zz|j|dz
�zz}|S)Nry@r)r�facr�rrkrr5)rrrR�twopijrs     r�polylog_continuationr��s����1�u���s�
�
�#�&�&�[�F�	���
�3�7�7�1�:�����Q����q�	�&�0@� A�A�A�
������A���G�G�A�J��
�w�w�q�z�A�~�#�'�'�!�*��/�c�g�g�a�j�A�o�	�V�C�F�F�1�I��!��$�
$�S�W�W�Q�q�S�\�
1�1���Hrc��|j�}|dkDr�|j}|j|�}|j}d}	||z
dk7r?|j	||z
�|z|j|�z}|rt
|�|krn||z
}||z}|dz
}�R||j|�|dz
z|j|dz
�z|j|dz
�|j|j|��z
zz
}n�|dkr�|j|�|j|�|dz
zz}|j|�}|j}	d}
	|j|
|z
dz�}|r4||	z|j|
�|
|z
dzzz}t
|�|krn||z}|	|z}	|
dz
}
�Xt�|j|�r|dkr|j|�}|S)Nrr)r�r4rr�rnr�ri�harmonicr�r�rkr)rrrRr�r��logz�logmzr�r��logkzr�r=s            r�polylog_unitcircler��s����7�7�(�C��1�u��H�H���v�v�a�y������
����!���z��x�x��!��}�u�,�s�w�w�q�z�9���C��I��O���T�	���T�M�E�
��F�A��	
�S�V�V�A�Y��1��
�c�g�g�a��c�l�
*�C�L�L��1��,=�c�f�f�c�f�f�Q�i�Z�>P�,P�
Q�Q��	
�Q���G�G�Q�B�K�#�&�&��)��q��s�+�+���v�v�a�y������
����
�
�a��c�!�e�$�A����w�����
�A�a�C��E� 2�3���t�9�s�?���T�	���T�M�E�
��F�A����
������A���G�G�A�J���Hrc�T�|j}|j|�}t|�dks�|j}d|z
}|j|�d|jz|zz}|j|�||z|j
|d|z�z||z|j
|d|z
�zzzd|jz|zzSd}d}	|j
||z
�|z}	t|	�|jkrn||	z
}|dz
}||z}||z}�E|jd|z
�||dz
zz|zS)N�rrr"r)r4rrirr�gammarnr�)
rr,rRrr�r�yr:r�r�s
          r�polylog_generalr��s=�����A����q�	�A��q�6�A�:��E�E��
�a�C���F�F�A�2�J��#�&�&���
�#���y�y��|�Q��T�#�(�(�1�S��U�"3�3�a�!��e�C�H�H�Q�s�1�u�<M�6M�M�N�PQ�RU�RX�RX�PX�[\�}�\�\�	�A�	�A�
��x�x��!��}�q� ���t�9�s�w�w���	�T�	��	�Q���	�Q���	�Q�����9�9�Q�q�S�>�A�2��1��+�%��)�)rc�v�|j|�}|j|�}|dk(r|j|�S|dk(r|j|�S|dk(r|d|z
zS|dk(r|jd|z
�S|dk(r|d|z
dzzSt	|�dks|j|�st	|�dkr
t
|||�St	|�dk\r8|j|�r'd|dzzt
||d|z�zt|||�zS|j|�rt|t|�|�St|||�S)Nr���rr��?g�������?gffffff�?)r#rn�altzetarrir+r�r�r�r*r�)rr,rRs   r�polylogr��s1�����A��A����A��A��A�v��x�x��{���B�w����A�����A�v��!�A�#�w���A�v����q��s��|���B�w��!�A�#���z��
�1�v��~�c�i�i��l�s�1�v��|��c�1�a�(�(�
�1�v��}����1���a��c�{�>�#�q�!�A�#�6�6�9M�c�ST�VW�9X�X�X�
�y�y��|�!�#�s�1�v�q�1�1��3��1�%�%rc��|j|�r|dkrt|�dzdk(r|dzS|r|j|�}n|j|�}|j	|�r2|j	|�r!|j|j
||��Sd|z}d|j
||�|j
||�z
zS)Nrrrr1)r+r*�expjpirmrk�imr��rr,rRrrr=s      r�clsinr��s���
�y�y��|��A��#�a�&�1�*��/���s�
�	��J�J�q�M���H�H�Q�K��
������ 1� 1�!� 4��v�v�c�k�k�!�A�&�'�'�	�!��A��C�K�K��!�$�s�{�{�1�Q�'7�7�8�8rc��|j|�r|dkrt|�dzdk(r|dzS|r|j|�}n|j|�}|j	|�r2|j	|�r!|j|j
||��Sd|z}d|j
||�|j
||�zzS)Nrrrr")r+r*r�rmrkr�r�r�s      r�clcosr��s���
�y�y��|��A��#�a�&�1�*��/���s�
�	��J�J�q�M���H�H�Q�K��
������ 1� 1�!� 4��v�v�c�k�k�!�A�&�'�'�	�!��A�����A�a� �3�;�;�q��#3�3�4�4rc�j�	|j|fi|��S#t$r|j|�cYSwxYwrB)�_altzetarl�_altzeta_generic)rr,rqs   rr�r�s=��'��s�|�|�A�(��(�(���'��#�#�A�&�&�'�s��2�2c��|dk(r|jd|zzS|jdd|z
�|j|�zS)Nrrr)�ln2�powm1rn)rr,s  rr�r�
s@���A�v��w�w��1��}���I�I�a��1���������+�+rNc�|�t|�}|dk(r|s|s	|j|fi|��S|j|�}|j}|jd�}|jd�}|s(|s&|j
d�|j|�dz
S|dk(r�|dk7r�t|j|��}	t|j|��}
t|	�d|zkDr
d|
z|kr|d	ks|d
k(r*		|rtd�|j||fi|��||_S|dk(r|jSt|�}||jk(r@|j|�|jk(r|dk(r|jS|j S|dzS|j#|�rd|zS|j|�d
|jzkDr'|dk(r"|s |j|j%d
|�zS|j&|||fi|���S#t$rY���wxYw#t$r|rtd�YnwxYw	||_��#||_wxYw)Nr�method�verboser"rzeuler-maclaurinrJr!rYzriemann-siegelz4zeta: Attempting to use the Riemann-Siegel algorithmz0zeta: Could not use the Riemann-Siegel algorithmr)r*�_zetarlr#r'rhr��_convert_paramrir5r�print�rs_zetar)r�r�r4r��power�_hurwitz)rr,rr;r�rqr<r'r�r�r��absss            rrnrns,���J��A��A�v�q�F�	��3�9�9�Q�)�&�)�)�	���A��A��8�8�D�
�Z�Z��
!�F��j�j��#�G�
�
��w�w�s�|�c�0�0��3�A�6�6�6��A�v�&�-�-�
������_��
������_���r�7�S��X��"�R�%�$�,�:��?��&�&�

 ����T�U�&�3�;�;�q�*�?��?� ����A�v��w�w���q�6�D��s�w�w���6�6�!�9������A�v��w�w���8�8�O���s�
�	���4����s�
�
�v�v�a�y�1�S�X�X�:��!�q�&���w�w����1�q�b�)�)�)��C�L�L��A�q�+�F�+�+�+��]#�	��	��6+����P�Q���������4���s5�G<�5 H�<	H	�H	�H%�"H2�$H%�%H2�2	H;c��|j}|jd�}	d}|xj|z
c_|j|�\}}|j|�dkr$|rt	d�	t|||||�||_S|rt	d�	||z|_t|||||dz|�\}	}
|j|	�|j|	|
z�z
}|r%t	d|	�t	d|
�t	d	|d
�||kr|	|
z||_Std|zt|dzd
|z��}||jdd
|z�kDr|jd����#t$rYnwxYw|s��t	d���#||_wxYw)Nr�r!rz#zeta: Attempting reflection formulazzeta: Reflection formula failedz)zeta: Using the Euler-Maclaurin algorithmzTerm 1:zTerm 2:z
Cancellation:�bitsrr�r��maxpreczzeta: too much cancellation)r'rhr�r�r��_hurwitz_reflectionrl�_hurwitz_emr�max�min�
NoConvergence)rr,rr<rqr'r��	extraprec�atype�T1�T2�cancellations            rr�r�Fs����8�8�D��j�j��#�G���	����I����%�%�a�(���5��6�6�!�9�q�=���;�<�
�*�3��1�a��?�,���#��=�>���i�'�C�H� ��a��A�t�B�w��@�F�B���7�7�2�;�����B���7�L���i��$��i��$��o�|�V�<��i�'��B�w����	 ��)��S���1A�3�t�8�-L�M�	��v�z�z�)�S��X�>�>��+�+�,I�J�J���
'�
��
����7�8��$���s>�AE5�,E�BE5�A	E5�	E#� E5�"E#�#E5�)E5�5	E>c�6���
��|dk7rt��j|�}|}�j|�r+t|�}|dkr�j	d|z
|�|dz
zS|dk(s|dk(st�d|z
�d}d}	|��j��dkDr'�dz�|�|zz}|	dz}	�j��dkDr�'�j��dkr'|�|zz
}�dz
�|	dz
}	�j��dkr�'	|j
\}
�
|
|	�
zz
}
d|
cxkr�
ksJ�J��j
���
�fd�td�
dz�D��}|d�j��zd�jz�
z�zzz}||z
}|S#|t|�k(sJ�t|�}
d�
Y��xYw)Nrr�Q�Zc3��K�|]6}�j�dzd|z�zz
��j�|�f�z���8y�w)rN)�cospir�)�.0r�r=r�qr:s  ����r�	<genexpr>z&_hurwitz_reflection.<locals>.<genexpr>�sF������
����1�Q�3�q��s�1�u�9�%�c�l�l�1�a��U�&;�;��s�<?r)
rlr��isnpintr*r��_mpq_�fsum�ranger�r)rr,rr<r��res�negsrr�shift�prGr=r�r:s`           @@@rr�r�ks�����A�v�!�!�

�&�&��)�C�
�2�D�
�{�{�1�~���H����6��<�<��!��Q�'�1�Q�3�/�/��S�L�E�S�L�!�!�	�!��A�	�A�
�E�	�A�

�&�&��)�a�-�	�Q���	�Q��W���
��
���&�&��)�a�-��&�&��)�q�.�	�Q��W���	�Q���
��
���&�&��)�q�.�
��w�w���1�
��q��L�A���;�Q�;��;��;������q��1����	�A���3�9�9�Q�<���3�6�6��!��a��	'�'�A���F�A��H����C��F�{��{���F��
��s�2E7�7Fc
��|j|�}|}d}|dz}|}	d}
|j|�rt|j|��}|dz
}	|j	|||z||z
dz
|g�dd}|
|z
}
||z}
|j|
�}||z}|g}d|z}d|
z}|
|z}|r"|j
|dz||z�||dzzz}nd||
|zzz}|d|z|zz
}dg}|}d}td|	dz�D�]}d|z}|dk(rdg}n
|dz
|dz
g}|D]�}t||dz�}||kr|j|d|z�dg|dzz}t|�D]}d|z
|z
||z||<�td|dz�D]}||xx||dz
z
||dz
zz
cc<� |}||z}��|j||�|z|j|�z|z}||z
}|j|�|kr|
d|z|zfcS||dz|dzzz}��|r td||d|j�d	|�||dz}}|j|�dkr|	|	dzz
}	���)
Nrrrr"rr�z
Sum range:zterm magnitude�	tolerance)r#r+r*r�_zetasumr�gammaincrr��appendr�fdotr�rr�r�) rr,rr<r'r�r��M1�M2�N�lsum�s1r��M2a�logM2a�logM2ad�logs�logr�rM2a�M2as�tailsum�Ur��factr�j2�updsr��D�Un�ir:s                                 rr�r��s������A��A��%�C�	
�B�	
���B�
�A��D�
�y�y��|������
�O��	
�1��B�
��L�L��B�q�D�"�R�%��'�A�3�/��2�1�5��
	
��	����d���������!�)���y����x����u���a�R�y����l�l�1�Q�3��6�	�2�R�!�A�#�Y�>�G��"�s�R�i��(�G��3��=�4�'�'��
�C�������q�!�A�#��	"�A��1��B��A�v��s���1��b��d�|���
����!�A�#�J����6��K�K��R��4��0��S�!�A�#�Y�����8�A�Q�q�S��U�A�a�D�L�B�q�E�8���!�A�#��A�A��1��!�Q�q�S�'�1�Q�q�S�6�1A�(A��A����T�	��
�����D�!�A�%��
�
�b�(9�9�D�5�A�A��q�L�G��w�w�q�z�C���b�1�W�w�.�.�.��R��T�B�q�D�M�!�D�)	"�*��,��B�(8�#�'�'�!�*�k�SV�W��R��T�B���6�6�!�9�q�=�
��A��I�A�arc�8���
��t�j|��d�jzkr	�j|�|||�S�j|d���|dgk7}t
|�dk(}|sg|s*�j��fd�t|dz�D��ggfS|r9|d�
�j���
�fd�t|dz�D��}d�
z|zggfSt|�}	|st|	dz�}|D�
cgc]}
�j��}}
|r|D�
cgc]}
�j��}}
ng}t|dz�D]�}
�|
z}|�z}|r!�j�j||zz�}|rx�j|�}|r)||	z}|dxx||zz
cc<|s�^|dxx|zz
cc<�o�j}|D])�
|�
xx||zz
cc<|r|�
xx|zz
cc<||z}�+��|dxx|z
cc<|s��|dxxz
cc<��||fS#t$rY���wxYwcc}
wcc}
w)	a�
    Returns [xd0,xd1,...,xdr], [yd0,yd1,...ydr] where

    xdk = D^k     ( 1/a^s     +  1/(a+1)^s      +  ...  +  1/(a+n)^s     )
    ydk = D^k conj( 1/a^(1-s) +  1/(a+1)^(1-s)  +  ...  +  1/(a+n)^(1-s) )

    D^k = kth derivative with respect to s, k ranges over the given list of
    derivatives (which should consist of either a single element
    or a range 0,1,...r). If reflect=False, the ydks are not computed.
    r"T�r�rrc3�.�K�|]}�|z�z���y�wrBrO)r�r�rrs  ��rr�z_zetasum.<locals>.<genexpr>�s�����>�Q�a��c�D�[�>�s�c3�^�K�|]$}�j�|z��z�|z�zz���&y�wrB)r)r�r�rrr<rs  ����rr�z_zetasum.<locals>.<genexpr>�s.�����K�!�����!���a��1�Q�3��+�5�K�s�*-r�)rir�r'�
_zetasum_fastrl�fnegr�r�rr�rr4�conjr�r)rr,rr�derivatives�reflect�have_derivatives�have_one_derivativer�maxdr<�xs�ysr�r��xterm�yterm�logwr:rs` `       `        @rrr�sB����3�6�6�!�9�~��c�h�h��&�	��$�$�Q��1�k�7�C�C��8�8�A�T�8�"�D�"�q�c�)���k�*�a�/�����H�H�>�&��1��+�>�>�?��C�C���A��A����K�v�a��c�{�K�K�A��!�G�a�K�=�"�$�$��{��D���D��F�m��'�	(�q�#�(�(�	(�B�	(�� +�
,�1�c�h�h�
,��
,�
��
�A�a�C�[���
��E���T�	����H�H�S�W�W��E�	�2�3�E���F�F�1�I�:�D�"��t�|���1�����%����q�E�U�T�\�)�E��G�G��$��A��q�E�U�Q�Y�&�E���1�����*����I�A�	�
�q�E�U�N�E���1�����-�.�r�6�M��W#�	��	��
)��
,s�H�4H�H�	H�Hc	���|j|�}t|�}t|�}|dkDrtd��|j}	|xjdz
c_|dk(rWd}|D]:}|s�|dk7s�d}|j
�}	|xjd|dzzzc_||	z
}�<|r|j�||_S|j}
td|dz�D]y}|||zs�|dk(rI|
|||z|j|||fd�|j|||f�|j|�zz
zz
}
�Z|
|||z|j|||f�zz
}
�{|
||zz}
||_|
�S#||_wxYw)Nrzarbitrary order derivativesr!rTF)r#r�r*rlr'r�r)r4rrnr9)rr,�chir;r�r<r'�	have_polerrurRrs            r�	dirichletr3s������A��A��C��A��J��A��1�u�!�"?�@�@��8�8�D�����B�����6��I��
����a�� %�I�����A��H�H��1�Q�3��'�H���F�A�
�����x����
�H�H���q��1���	7�A��1�Q�3�x���6���Q�q�S��S�X�X�a�!�A���%:�����Q�q�E�*�3�7�7�1�:�5�&6�7�7�A���Q�q�S��C�H�H�Q��1��$6�6�6�A�
	7�	
�Q��T�	�����
�2�I�����s%�#E�)E�/?E�6&E�A8E�	E(c�f������j}����fd�}�jx}}�j}d}	||kDrC||z
}|	dz
}	�j�j	|	���||	�}t|�}||kDr�Cd}
|j
d�r��j��}d�jdzztd|�z�|dz
zz�j�d�jzz�z�jd��dzz�zt�j�dz��z}
t|
�}
|�|
|	fS)	Nc�L���jd�z��dzzd����zzS)Nr"rT)�regularized)r)rrr�gammr,s ����rrCz&secondzeta_main_term.<locals>.<lambda>7s-���#�,�,�s�1�u�a��a��i�T�,�B�4�1�"�:�M�rrr�errorr"r�rr2)
r�r4r)r��zetazero_memoizedrirhr�rr�r9rr�)
rr,rrqr�r�totsumr��mgrr��sgr7s
```         @r�secondzeta_main_termr=5s)���

�'�'�C�M�A��H�H��F�T�	���B�	�A�
�s�(��$���	�Q����v�v�c�+�+�A�.�/����t��
��Y��
�s�(��C�
�z�z�'��
�V�V�A�Y���#�&�&�2�,��s�1�R�y�(��R��V��4�S�W�W�T�1�S�V�V�8�_�5M�M�
�\�\�$��$��'�	�
*�+�+.�s�y�y��1��~�+>�?���#�h���7�C��?�rc�"�����j}���fd�}�jx}}�j}d}	||kDs|	dkr;||z
}|	dz
}	||	�}|dk(r
�j}nt|�}||kDr�5|	dkr�;|j	d�r|}
|�
|	fS)Nc�P���jdd�z
zd�j|�dzz�dzz�d�j|�z�dz
zz�j|�z�j|�zd�j	d�z�z�j�j
�zzS)Nr"rr0rr�)rr9�mangoldtr�r�r�rrrr,s ���rrCz'secondzeta_prime_term.<locals>.<lambda>Ks����#�,�,�s�A�a�C�y��c�g�g�a�j�!�m�);�a�"�g�)E�F�

�c�g�g�a�j�.�A�a�C�	 �"�"%�,�,�q�/�2�25�(�(�1�+�>�	
�3�9�9�S��U��	�C�H�H�S�V�V�,�	,�.�rr�	rr8)r�r4r)rirh)rr,rrqr�rr:r�r;rr�s```        r�secondzeta_prime_termrCIs����

�'�'�C�	.�A��H�H��F�T�	���B�	�A�
�s�(�a�!�e��$���	�Q�����t���1�9����B��T��B�
�s�(�a�!�e��z�z�'�����7�C��?�rc�������j��rT�j��dkr@tt�j����}|dzs�j	d�|dzzS�j
}���fd�}�j}|d�}�j}d}	||kDr#||z
}|	dz
}	||	�}t|�}||kDr�#�d�zz|z�jd�z�z}
|
S)Nrrz-0.25rc�J��d�z|z|d�zz�j|�zzS)Nr0r")r�rAs ���rrCz%secondzeta_exp_term.<locals>.<lambda>cs)���4��6�A�+��#�a�%�������3�4�rr")
r+r�r*r|r�r�r4r)rir�)rr,rr�r�rr:r�r;rrs```        r�secondzeta_exp_termrF]s����
�y�y��|����q�	�Q����c�f�f�Q�i� �!���1�u��7�7�7�#�q�b�!�e�,�,�

�'�'�C�4�A�
�X�X�F��Q�4�D�	���B�	�A�
�s�(��$���	�Q�����t��
��Y��	
�s�(�
	
�C��E�
�6��#�)�)�C��E�*�*�A��Hrc	������d�dz
zzd�j�j�z�jd�z�zz}�j|�}�xj|z
c_�d�dz
zzd�j�j�z�jd�z�zz}�j
}���fd�}�j}�j}	d}
||
�}t|�}||kDrA||	kr<||z
}|
dz
}
||
�}||z
}|
dz
}
||
�}|}	t|�}||kDr||	kr�<||z
}d�dz
dzz�j�jd�jdzz�z�z�dz
dzzz}
||
|zz}d	}|jd
�r�||kDr||	ks|||kr9�jd�t�jt||z�d��z}||	kDr9�jd�t�jt||	z�d��z}t|�j
dz�}�xj|zc_|�|fS)
Nr"rrYc����j|d�d�j��z|zz�jd|z�z�|zdz
�j|�zzS)Nr�rYr"r)r�r�r�r�rAs ���rrCz*secondzeta_singular_term.<locals>.<lambda>vsZ���#�,�,�q��&��#�(�(�1�+�
��'9�9�
�y�y��Q�����A�#�a�%������+�-�rr��rr�rr8r!��?)r�rr�rr'r�r4r)rir&r9rhr�r*r�)rr,rrq�factorr�r�rr:�mg1rr��mg2�pole�str�s```             r�secondzeta_singular_termrPps6���
��a��c��^�Q�s�x�x����/�/��	�	�#�a�%�0@�@�
A�F������I��H�H�	��H�
��a��c��^�Q�s�x�x����/�/��	�	�#�a�%�0@�@�
A�F�

�'�'�C�	-�A�
�X�X�F�

�'�'�C�	�A��Q�4�D�

�d�)�C�

��)��s�
��$���	�Q�����t���$���	�A�����t�����$�i����)��s�
��d�N�F�
�q��s�b�k�>�3�9�9�S�W�W�R�����	�\�!�^�%<�<�q��s�b�k�I�I�D�	��f��	�B�
�C�
�z�z�'���s������c�z��g�g�b�k�3�s�w�w�s�6�#�:��r�'B�#C�C���S�y��g�g�b�k�3�s�w�w�s�6�#�:��r�'B�#C�C���#�s�w�w�r�z�"���H�H�	��H��3��8�Orc��|j|�}|j|�}|j}|j|�r�|j|�dkr�t	|dz
�|dzkr|j
St
t|j|���}|dzr|j
Sd|dzz|jd|j|d��z
d|dzz�zS|j}	t|||�}t|j|�d	�}|xj|dzz
c_
t|||d
d
��\}	}
}t|||d
d
��\}}
}t!|||d
��\}}t|||�}|
|
z|z}|	|z
|z|z
}|j#d
�rJt%d|	�t%d|d�t%d|�t%d|d�t%d|�t%d|�||_
|j#d�rNt|jt	|��d	�}t|d|zz|jdzd|zz�}|�|fS|�S#||_
wxYw)a�

    Evaluates the secondary zeta function `Z(s)`, defined for
    `\mathrm{Re}(s)>1` by

    .. math ::

        Z(s) = \sum_{n=1}^{\infty} \frac{1}{\tau_n^s}

    where `\frac12+i\tau_n` runs through the zeros of `\zeta(s)` with
    imaginary part positive.

    `Z(s)` extends to a meromorphic function on `\mathbb{C}`  with a
    double pole at `s=1` and  simple poles at the points `-2n` for
    `n=0`,  1, 2, ...

    **Examples**

        >>> from mpmath import *
        >>> mp.pretty = True; mp.dps = 15
        >>> secondzeta(2)
        0.023104993115419
        >>> xi = lambda s: 0.5*s*(s-1)*pi**(-0.5*s)*gamma(0.5*s)*zeta(s)
        >>> Xi = lambda t: xi(0.5+t*j)
        >>> chop(-0.5*diff(Xi,0,n=2)/Xi(0))
        0.023104993115419

    We may ask for an approximate error value::

        >>> secondzeta(0.5+100j, error=True)
        ((-0.216272011276718 - 0.844952708937228j), 2.22044604925031e-16)

    The function has poles at the negative odd integers,
    and dyadic rational values at the negative even integers::

        >>> mp.dps = 30
        >>> secondzeta(-8)
        -0.67236328125
        >>> secondzeta(-7)
        +inf

    **Implementation notes**

    The function is computed as sum of four terms `Z(s)=A(s)-P(s)+E(s)-S(s)`
    respectively main, prime, exponential and singular terms.
    The main term `A(s)` is computed from the zeros of zeta.
    The prime term depends on the von Mangoldt function.
    The singular term is responsible for the poles of the function.

    The four terms depends on a small parameter `a`. We may change the
    value of `a`. Theoretically this has no effect on the sum of the four
    terms, but in practice may be important.

    A smaller value of the parameter `a` makes `A(s)` depend on
    a smaller number of zeros of zeta, but `P(s)`  uses more values of
    von Mangoldt function.

    We may also add a verbose option to obtain data about the
    values of the four terms.

        >>> mp.dps = 10
        >>> secondzeta(0.5 + 40j, error=True, verbose=True)
        main term = (-30190318549.138656312556 - 13964804384.624622876523j)
            computed using 19 zeros of zeta
        prime term = (132717176.89212754625045 + 188980555.17563978290601j)
            computed using 9 values of the von Mangoldt function
        exponential term = (542447428666.07179812536 + 362434922978.80192435203j)
        singular term = (512124392939.98154322355 + 348281138038.65531023921j)
        ((0.059471043 + 0.3463514534j), 1.455191523e-11)

        >>> secondzeta(0.5 + 40j, a=0.04, error=True, verbose=True)
        main term = (-151962888.19606243907725 - 217930683.90210294051982j)
            computed using 9 zeros of zeta
        prime term = (2476659342.3038722372461 + 28711581821.921627163136j)
            computed using 37 values of the von Mangoldt function
        exponential term = (178506047114.7838188264 + 819674143244.45677330576j)
        singular term = (175877424884.22441310708 + 790744630738.28669174871j)
        ((0.059471043 + 0.3463514534j), 1.455191523e-11)

    Notice the great cancellation between the four terms. Changing `a`, the
    four terms are very different numbers but the cancellation gives
    the good value of Z(s).

    **References**

    A. Voros, Zeta functions for the Riemann zeros, Ann. Institute Fourier,
    53, (2003) 665--699.

    A. Voros, Zeta functions over Zeros of Zeta Functions, Lecture Notes
    of the Unione Matematica Italiana, Springer, 2009.
    rr�r�rr@Tr rr�True)r8r�)r8r�zmain term =z    computed usingz
zeros of zetazprime term =z#values of the von Mangoldt functionzexponential term =zsingular term =r8rJ)r#r�r+r�rir)r*r|�fractionr�r'rFr�rr=rCrPrhr�)rr,rrqr�r�r'�t3r�rr�r1�gtrs�r2�pt�t4�r4r�r:r�s                    r�
secondzetar[�s9��x	���A��A����A��A�

�'�'�C�
�y�y��|����q�	�Q���q��s�8�c�$�h���7�7�N���c�f�f�Q�i� �!���q�5��7�7�N��A�2�q�5�M��<�<��#�,�,��r��,�"=� =�a�1�"�Q�$�i�H�I�
J��8�8�D��
 ��a��
+��������A�&�	����I��M�!��)�#�a����O�
��B��*�3�q����P�
��B��)�#�a���?���B�
 ��a��
+����e�B�h���r�E�"�H�R�K���:�:�i� ��-��$��&��O�<��.�"�%��&��,Q�R��&��+��#�R�(����
�z�z�'�������A����"���#�a��d�(�C�G�G�B�J�q�!�t�O�,���r�3�w��
�2�I�����s
�-C5I�	Ic�����	�
��|dk(r��zS|dk(r�j���S�dk(r�j�|�|zS�j��dkr��j��rt	d��t�j
d�j��z
��}�j}�j}t|�D]}||�|z�zzz
}||z}�|�j|��|z�z|zS�j|��	dd��zzz�jd�z
��	z��	�dz
zz|�zzz}�dz�
d�jz����	�
��fd�}|d�j|d�jg�zz
}�j!|�sG�j!��s6�j!��s%�j|�dkr�j#|�}|S)a�

    Gives the Lerch transcendent, defined for `|z| < 1` and
    `\Re{a} > 0` by

    .. math ::

        \Phi(z,s,a) = \sum_{k=0}^{\infty} \frac{z^k}{(a+k)^s}

    and generally by the recurrence `\Phi(z,s,a) = z \Phi(z,s,a+1) + a^{-s}`
    along with the integral representation valid for `\Re{a} > 0`

    .. math ::

        \Phi(z,s,a) = \frac{1}{2 a^s} +
                \int_0^{\infty} \frac{z^t}{(a+t)^s} dt -
                2 \int_0^{\infty} \frac{\sin(t \log z - s
                    \operatorname{arctan}(t/a)}{(a^2 + t^2)^{s/2}
                    (e^{2 \pi t}-1)} dt.

    The Lerch transcendent generalizes the Hurwitz zeta function :func:`zeta`
    (`z = 1`) and the polylogarithm :func:`polylog` (`a = 1`).

    **Examples**

    Several evaluations in terms of simpler functions::

        >>> from mpmath import *
        >>> mp.dps = 25; mp.pretty = True
        >>> lerchphi(-1,2,0.5); 4*catalan
        3.663862376708876060218414
        3.663862376708876060218414
        >>> diff(lerchphi, (-1,-2,1), (0,1,0)); 7*zeta(3)/(4*pi**2)
        0.2131391994087528954617607
        0.2131391994087528954617607
        >>> lerchphi(-4,1,1); log(5)/4
        0.4023594781085250936501898
        0.4023594781085250936501898
        >>> lerchphi(-3+2j,1,0.5); 2*atanh(sqrt(-3+2j))/sqrt(-3+2j)
        (1.142423447120257137774002 + 0.2118232380980201350495795j)
        (1.142423447120257137774002 + 0.2118232380980201350495795j)

    Evaluation works for complex arguments and `|z| \ge 1`::

        >>> lerchphi(1+2j, 3-j, 4+2j)
        (0.002025009957009908600539469 + 0.003327897536813558807438089j)
        >>> lerchphi(-2,2,-2.5)
        -12.28676272353094275265944
        >>> lerchphi(10,10,10)
        (-4.462130727102185701817349e-11 - 1.575172198981096218823481e-12j)
        >>> lerchphi(10,10,-10.5)
        (112658784011940.5605789002 - 498113185.5756221777743631j)

    Some degenerate cases::

        >>> lerchphi(0,1,2)
        0.5
        >>> lerchphi(0,1,-2)
        -0.5

    Reduction to simpler functions::

        >>> lerchphi(1, 4.25+1j, 1)
        (1.044674457556746668033975 - 0.04674508654012658932271226j)
        >>> zeta(4.25+1j)
        (1.044674457556746668033975 - 0.04674508654012658932271226j)
        >>> lerchphi(1 - 0.5**10, 4.25+1j, 1)
        (1.044629338021507546737197 - 0.04667768813963388181708101j)
        >>> lerchphi(3, 4, 1)
        (1.249503297023366545192592 - 0.2314252413375664776474462j)
        >>> polylog(4, 3) / 3
        (1.249503297023366545192592 - 0.2314252413375664776474462j)
        >>> lerchphi(3, 4, 1 - 0.5**10)
        (1.253978063946663945672674 - 0.2316736622836535468765376j)

    **References**

    1. [DLMF]_ section 25.14

    rrz#Lerch transcendent complex infinityrc����j��j|�z�z|�zz
��dz|dzz�z�j�|z�zzSrN)�sin�atan�expm1)r:rrrGrur�r,s ������rrCzlerchphi.<locals>.<lambda>}sR���#�'�'�!�C�H�H�Q�q�S�M�/�!�A�#�-�.�
�Q�$�q�!�t�)�a��#�)�)�A�a�C�.�	(�*�r)rnr�r�r�r�r*r�r4r�r�lerchphirrrr(r)r��chop)rrRr,rr�r�zpowrrrGrur�s` ``     @@@rraras����b	�A�v��a�R�y���A�v��x�x��1�~���A�v��{�{�1�a� �1�$�$�
�v�v�a�y�1�}��;�;�q�>��B�C�C������3�6�6�!�9��%�&���H�H���w�w�����	�A�
���1��q��� �A��A�I�D�	��c�l�l�1�Q��!��,�,�q�0�0����q�	�A�	�1�Q��T�6�
�S�\�\�!�A�#��r�!�t�,���a��c�{�:�Q��T�A�A�A�	�A��A�	�#�&�&��A�	*�	*�A���3�8�8�A��3�7�7�|�$�	$�$�A��6�6�!�9�S�V�V�A�Y�s�v�v�a�y�S�V�V�A�Y��]��H�H�Q�K���Hr)r)r)z2http://www.dtc.umn.edu/~odlyzko/zeta_tables/zeros1)F)rrN)rr)g���Q��?)+�
__future__r�
libmp.backendr�	functionsrrr	r.r>rHrgr}r�r�r�r�r�r�r�r�r�r�r�r�r�r�r�r�r�r�rnr�r�r�rr3r=rCrFrPr[rarOrr�<module>rgs@��%�"�9�9��#��#�J�+��+�B�B��B��J%��J%�Z��.��
E��
E��
��
�,�#��#������%��%�F�*)��*)�X�%)��%)�N�.��.�
�	
�"
�H*�(�&��&�*�
9��
9��
5��
5��'��'��,��,�
�3,��3,�j�"��"�H'
�R=�B�()�s�E�;��;�z��#�!����B�(�(
�&"�H����D�j
��j
r

Youez - 2016 - github.com/yon3zu
LinuXploit