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__/hypergeometric.cpython-312.pyc
�

�"`r�����ddlmZddlmZmZd�ZdZegdfd��Zed��Zed	��Z	ed
��Z
ed��Zed��Zed
��Z
ed��Zed��Zed��Zed��Zed��Zed��Zd�Zed��Zed��Zed��Zed��Zed��Zed��Zed��Zed$d��Zed��Zed��Zed ��Zed!��Z ed"��Z!	ed#��Z"y)%�)�xrange�)�defun�
defun_wrappedc��dx}}d}g}t|�D�],\}}	|	\}
}}}
}}}d}t|
�D],\}}|r�	|j||�dks�!||s�'dx}}d}�.gd�}t||
|g�D]�\}}t|�D]\}}|j|�\}}|dkDr� ||jk(rCd}|dk(r+|D]&}|j	|�s�|t|�k\s�$d}n|r�d||xxdz
cc<�r|dks�x||z
}d}����|r&|d|d|dzkDr|s|j
|���t|�s��)dx}}��/||||fS)NF�T)rrrrr���)�	enumerate�re�
nint_distance�ninf�isnpint�int�append�sum)�ctx�terms�prec�discard_known_zeros�perturb�	recompute�	extraprec�discard�
term_index�term�w_s�c_s�alpha_s�beta_s�a_s�b_s�z�have_singular_nongamma_weight�k�w�
pole_count�
data_index�data�i�x�n�d�ok�us                             �A/usr/lib/python3/dist-packages/mpmath/functions/hypergeometric.py�_check_need_perturbr0s�����G�i��I��G�%�e�,�('��
�D�15�.��S�'�6�3��Q�(-�%��c�N�	9�D�A�q���6�6�#�a�&�>�Q�&�3�q�6�*.�.�G�i�48�1�		9�
�
� )�7�F�C�*@� A�	%��J��!�$��
%���1��(�(��+���1��q�5������=��B�!�Q��!$�&�A�"�{�{�1�~�!�s�1�v�+�%)�� %�&�� ��z�*�a�/�*���V��!��O�I� $�I�-
%�	%�0�:�a�=�:�a�=�:�a�=�3P�#P�1��N�N�:�&�
��_�"&�&�G�i�Q('�R�I�y�'�1�1�a
hypercomb() failed to converge to the requested %i bits of accuracy
using a working precision of %i bits. The function value may be zero or
infinite; try passing zeroprec=N or infprec=M to bound finite values between
2^(-N) and 2^M. Otherwise try a higher maxprec or maxterms.
Tc�d�|j}|j}|j}|j}|dd}	|j	dd�}
|j	d|j|��}||d<|j	d�}|j	d�}
d}d}		|xjdz
c_|j|kDrt
t||jfz��|j}|	dd}||�}|
r7t�td	�td
|j�td|�t||||�\}}}}|xj|z
c_|r�d|vr|d}n-|jrt|jd
z�}n|dz|z}|j|j|�}|dz|zdz|_tt|��D]}||xx|z
cc<|||dzzz
}�|r||�}|r"t!|�D��cgc]\}}||vs�|��}}}|s|j||_Sg}t!|�D�]�\}}|\}}} }!}"}#}$|
r�t�td|dzt|�t|"�t|#�fz�td|j#|�|j#|��td|j#| �|j#|!��td|j#|"�|j#|#��td|j#|$��|j%|j&|"|#|$fi|��g| D�%cgc]}%|j)|%���c}%z|!D�&cgc]}&|j+|&���c}&zt-||�D�'�(cgc]\}'}(|j/|'|(���c}(}'z�})|
rtd|)�|j1|)����t|�dk(r	|s|d}�n�|jr|j3|�}�nx|j3|�}|D�*cgc]}*|j5|*���}+}*t7|+�},|j5|�}-|,|-z
}.|
r1t�td|.d�td|j|z
d�|.|j|z
k}/|�d}0n|,|jz
|k}0|
�d}1n|,|
kD}1|/r|r|j9|.�rn�|/r`|�	|dz
}|}��N|j5||z
�|j5|�|z
krnt|0r
|j}ne|1r
|j:}nVd|vrnQ|dz}|}nGt=t7|.|dz�t7||��}2|xj|2z
c_|
rtd������||_|�Scc}}wcc}%wcc}&wcc}(}'wcc}*w#||_wxYw)N�verboseF�maxprec�zeroprec�infprecrr�
zENTERING hypercomb main loopzprec =�hextra�hmagg333333�?z  Evaluating term %i/%i : %iF%iz
    powersz	    gammaz	    hyperz    zz
    Value:z  Cancellation:�bitsz  Increased precision:�rz*  Must start over with increased precision)r�zerorr
�get�_default_hyper_maxprec�
ValueError�_hypercomb_msg�printr0�_fixed_precisionr�ldexp�one�range�lenr
�nstr�fprod�hyper�gamma�rgamma�zip�powerr�fsum�mag�max�isnan�inf�min)3r�function�paramsr�kwargs�orig�sumvalue�distr
�orig_paramsr3r4r5r6�perturbed_reference_valuer8�orig2rrrrrr9�hr$r)r�evaluated_termsr�	term_datarrrrr r!r"�a�br%�c�vr*�term_magnitudes�
max_magnitude�
sum_magnitude�cancellation�precision_ok�zero_ok�inf_ok�	increments3                                                   r/�	hypercombrl:s����8�8�D��x�x�H����D��8�8�D���)�K��j�j��E�*�G��j�j��C�$>�$>�t�$D�E�G��F�9���z�z�*�%�H��j�j��#�G� $��
�F�w���H�H��N�H��x�x�'�!� ��4����2B�!B�C�C��H�H�E� ��^�F��f�%�E�����4�5��h����)��h��'�#�C���6I�J�
3�G�Y�	�7��H�H�	�!�H���V�#�!�&�>�D��)�)��s�x�x��|�,�D��"�9�v�-�D��I�I�c�g�g��u�-�� �2�:��,�r�1����s�6�{�+�!�A��1�I��N�I�
��A�a�C��L�A�
!�� �&�)��"�/8��/?�T�)�1�d�1�G�CS��T��T���x�x�d���c!�O�)2�5�)9�
*�%�
�I�9B�6��S�'�6�3��Q���G��;�#�A��s�5�z�3�s�8�S��X�F�G�H��,�����
�s�x�x��}�E��+�s�x�x��'8�#�(�(�6�:J�K��+�s�x�x��}�c�h�h�s�m�D��'�3�8�8�A�;�/�
�I�I�y�s�y�y��c�1�?��?�@�+2�3�a�S�Y�Y�q�\�3�4�,2�3�q�S�Z�Z��]�3�4�25�S���>���1�S�Y�Y�q��^�>�?�@����,��*��&�&�q�)�)
*�,�5�z�Q���*�1�-����#�#��8�8�O�4����x�x��0�H�3B�C�a�s�w�w�q�z�C�O�C���0�M��G�G�H�-�M�(�=�8�L�����'��v�>��.����4���H�'�#�(�(�T�/�9�L�����'�#�(�(�2�h�Y�>������&��0���W����<�1H���,�4��b�L�F�08�-���W�W�X�(A�A�B�����)�D�0�1���"�x�x�H���"�w�w�H���v�%���a�K�F�08�-� ��L�$��'� :�C�	�$�<O�P�	����I�%����F�G��i�l����9���kU��&4��3��>��D��b���s^�E(V&�=V�
V�V&�&DV&�V�V&� V�8V&�V�'A7V&�V!�6EV&�V&�&	V/c��|j|�}t|�}t|�}|D�cgc]}|j|���}}|D�cgc]}|j|���}}|jdd�rt|jdd�}	d}
|
|kr[|rY||
}||vrC|	s|j	|d�s-|j|�|j|�|dz}|dz}n|
dz
}
|
|kr|r�Y|dk(r0|dk(r|j||fi|��S|dk(�r#|j|�S|dk(rQ|dk(r|j|||fi|��S|dk(r|j|||fi|��S|dk(r�|j|dd|�S|dk(rh|dk(r|j|||fi|��S|dk(r|j|||fi|��S|dk(r|j|||fi|��S|dk(rd|j|||fi|��S||dzk(r|j|||||fi|��S||dzkDr(|jd	�s|j |||||fi|��St#||z�\}}|j$|||||fi|��Scc}wcc}w)
z0
    Hypergeometric function, general case.
    �	eliminateT�
eliminate_allFrrr��force_series)�convertrF�_convert_paramr=r�remove�_hyp0f1�exp�_hyp1f1�_hyp1f2�_hyp1f0�_hyp2f1�_hyp2f2�_hyp2f3�_hyp2f0�_hypq1fq�
_hyp_borelrL�hypsum)
rr r!r"rV�p�qr`ra�elim_nonpositiver)�coeffs�typess
             r/rIrI�s���
	���A��A��C��A��C��A�*-�
.�Q�3���a� �
.�C�
.�*-�
.�Q�3���a� �
.�C�
.�
�z�z�+�t�$�!�:�:�o�u�=��
���!�e���A��A��C�x�-�S�[�[��1��5F��
�
�1�
��
�
�1�
��Q����Q����Q����!�e��	�A�v�
�!�V�K�C�K�K��Q�9�&�9�9�
�!�V�C�G�G�A�J�&�	
�a��
�!�V�K�C�K�K��S�!�>�v�>�>�
�!�V�K�C�K�K��S�!�>�v�>�>�
�!�V�C�K�K��A��q�	�1�5�5�	
�a��
�!�V�K�C�K�K��S�!�>�v�>�>�
�!�V�K�C�K�K��S�!�>�v�>�>�
�!�V�K�C�K�K��S�!�>�v�>�>�
�!�V�K�C�K�K��S�!�>�v�>�>�	
�a��c���s�|�|�A�q�#�s�A�8��8�8�	
�Q�q�S�����N�3��s�~�~�a��C��a�:�6�:�:��#�c�'�O�M�F�E��3�:�:�a��E�6�1�7��7�7��C/��
.s�I�
Ic�.�|jg|g|fi|��S�N�rI)rrar"rVs    r/�hyp0f1r��s���3�9�9�R���A�'��'�'r1c�0�|j|g|g|fi|��Sr�r��rr`rar"rVs     r/�hyp1f1r��s ���3�9�9�a�S�!��Q�(��(�(r1c�2�|j|g||g|fi|��Sr�r�)r�a1�b1�b2r"rVs      r/�hyp1f2r��s"���3�9�9�b�T�2�b�'�!�-�f�-�-r1c�2�|j||g|g|fi|��Sr�r�)rr`rarbr"rVs      r/�hyp2f1r��s"���3�9�9�a��U�A�3�q�*�6�*�*r1c�4�|j||g||g|fi|��Sr�r�)rr��a2r�r�r"rVs       r/�hyp2f2r��s$���3�9�9�b��W�b��W�Q�0��0�0r1c�6�|j||g|||g|fi|��Sr�r�)rr�r�r�r��b3r"rVs        r/�hyp2f3r�s&���3�9�9�b��W�b��B�Z��3�F�3�3r1c�0�|j||gg|fi|��Sr�r�r�s     r/�hyp2f0r�s ���3�9�9�a��U�2�a�)�&�)�)r1c�6�|j|||g||g|fi|��Sr�r�)rr�r��a3r�r�r"rVs        r/�hyp3f2r�	s&���3�9�9�b��B�Z��B���3�F�3�3r1c��d|z
|zS�Nr�)rr`r"s   r/ryry
s��
�a�C�a�R�=�r1c�J���	�|\\�	}�r�j��}nd}|dk\r�|jd�s�	�j}	�xjd|dzzz
c_�	��fd�}�j|gd��}�j	�	�d�j�j�zz|z}|�_�j�	�r"�j��r�j|�}|�S�jdd	|f�	g�fi|��S#|�_wxYw#�j$rY�6wxYw)
Nr�rq�rc�L���j�	�}�j|z}dd|zz}�j�z
}�jd|z�}||g|dggg��jz
�j�z
gg|f}||g|dggg��jz
�j�z
gg|f}||fS)Nr�r���)�sqrt�j�mpq_1_2rv�mpq_3_2)
r%�jwr.rb�E�T1�T2rarr"s
       ���r/r]z_hyp0f1.<locals>.h!s�������!���A����q��B��1�R�4��A����a��A�����"��
�A��3�q�'�A�b�6�2�r�A�c�k�k�M�3�;�;�q�=�3Q�SU�XY�WY�Z�B��a�&�1�Q�%��R�!�C�K�K�-����Q��1O�QS�UV�W�B��r�6�Mr1T�rqr)rOr=rrlrJr��pi�
_is_real_type�_re�
NoConvergencer�)
rr!r"rV�btype�magzrWr]rcras
` `      @r/rurus"����K�J�Q����w�w�q�z�����q�y����N�3�	�
�8�8�D�
 ����B��q��L�(��"��M�M�!�R�d�M�;���I�I�a�L�!�C�H�H�S�V�V�$4�"4�5�a�7������ � ��#��(9�(9�!�(<��G�G�A�J���2�I��3�:�:�a��U�H�q�c�1�7��7�7��
 ����� � �	��	�s*�D�A)D�.<D�	D
�
D�D"�!D"c�����|\\}}|\\}}�s�j�zS�j��}	|	dk\�r5�j|�r�j|�dk�s�j	��rZ�j|��j|�cxk(r%�j��cxk(rdk(r
�jS�j�zS		�xj|	z
c_�j��dk����fd�}
�j|
||gd��}�j|�r3�j|�r"�j��r�j|�}|��xj|	zc_S�jdd||f||g�fi|��}|S#�j$rYnwxYw	�xj|	zc_�I#�xj|	zc_wxYw)N�rrc����r#�j�j|d���}n�j|�}d�z}|�gd|g|g||z
g|d|z|z
gg|f}�j���gd||z
g|g|g||z
d|z
gg|f}||fS�NT)�exactr)�expjpi�fnegrv)	r`rar��rzr�r�r�sectorr"s	      ���r/r]z_hyp1f1.<locals>.hEs������J�J�s�x�x���x�'>�?���J�J�q�M���1��B��Q�%�!�Q�B��!��q��s�e�a��1��Q��Z��b�S�I�B��7�7�1�:�a�.�1�Q�q�S�'�A�3���a��c�1�Q�3�Z��R�P�B��r�6�Mr1Tr�)rDrO�isintr�isinf�signrR�nanr�_imrlr�r�r�r�)
rr r!r"rVr`�atyperar�r�r]rcr�s
`  `        @r/rwrw5s�����K�J�Q���K�J�Q����w�w�q�y���7�7�1�:�D��q�y�#�)�)�A�,�3�6�6�!�9��>��9�9�Q�<��x�x��{�c�h�h�q�k�=�S�X�X�a�[�=�A�=��w�w��>��7�7�Q�;��	�
����D� �������a���"��M�M�!�a��U��M�>���$�$�Q�'�C�,=�,=�a�,@�S�EV�EV�WX�EY�����
�A��r�
�H�H���H���
�
�1�a�%���!�Q���=�f�=�A��H���$�$�
��
����H�H���H��C�H�H���H�s%�
BF�F�F8�F�F8�8Gc��||||f\}}}}	|j}
|jdd|
z�}d}	|
|z|_|j|	�}|jd�}
|jd�}|jd�}d}||z|z}||jz}|dz|jz}d|z
}||z}||z|z}||z
|z
}|dz
}|jdz
}|j
}|j}|j}|j}	||z}||z||zz|zd|dzz|z||zzz}||||z|
z|zz
z} ||
|z||z|zzz}!|
|||z||zz|z
z|z
zd|z|zz}"||z|"z
}#t|||#��}||#|z
�|krn| |!|#}}}
|dz
}��|||#�z
}$|$|kr	|#S||$z
}||kDr|j����)Nr4�dr7rrrr�)rr=rr�mpfr�rG�nprintrOr
rPr�)%rr`rarbr"rV�_a�_b�_c�_zrWr4�extrar,�e�fr$�abz�ch�c1h�nz�g�abg�cba�z2�tolrGr�rO�maxmag�kch�kakbz�d1�e1�ft�f1rgs%                                     r/�_hyp2f1_gosperr�YsK���Q��1�*�K�B�r�"�R��8�8�D��j�j��C��H�-�G��E�
��%�<���
�K�K��O���G�G�A�J���G�G�A�J���G�G�A�J��
����c�!�e��
����_����s�c�k�k�!��
�q�S��
�b�D����c�!�e����c�!�e��
�q�S���x�x�i�"�n���x�x�������g�g��������B�$�C��q�S�1�Q�3�K��M�Q��!��W�S�[�!�C�%�%8�9�E���1�S�5�!�)�A�+�
�&�B���#��q��s�A�g�
�&�B��A�s�1�u�Q�r�T�z�!�|�$�S�(�)�1�S�5��8�4�B��Q����B����R��)�F��2�a�4�y�3����"�b�!�q�A�
��F�A����B��'���%���
�I�
�\�!�E��w���'�'�'�Yr1c������|\\}}\}}|\\�}	�dk(�r
�j�|z
|z
�dkD}
�j|�xr|dkxs�j|�xr|dk}�j��xrP�dkxrI�j|�xr�|cxkxrdkncxs#�j|�xr�|cxkxrdknc}|
s|r&|s$�j��|z
|z
g�|z
�|z
gd��S�j||�d�jdzz
��j
zS�s�s
|dk(s|dk(rd�zS�jS�j��rZ�dkrU�j|�r
�|cxkrdks,n�j|�r�|cxkrdkrn
�j
Sn�j
St��}
|
dks6�j|�r
|dkr|dk\s�j|�r'|dkr"|dk\r�jdd|||	f||�g�fi|��S�j}	�xjdz
c_	|
d	k\r���fd
�}�j|||gfi|��}nytd�z
�dkr��fd�}�j|||gfi|��}nKt��dz
z�dkr'�j|�|z
���dz
z�d�z
|zz}nt�||��fi|��}|�_	|�S#|�_	wxYw)
NrrT)�_infsignrg�������?i���r7g�������?c�����j�z
}||z
}d�	z}�	g|g�|g|�|z
g|||zg�j|zg|f}�	g|g�|g|�|z
g|||zg�j|z
g|f}||fSr�)�mpq_1)
r`ra�t�abr�r�r�rbrr"s
       ���r/r]z_hyp2f1.<locals>.h�s�����I�I�a�K��a��c���!��2��r�d�Q�B�4�!�R�C��!�A�a�C��1�Q�q�S�'�3�9�9�R�<�.�2�N���r�d�Q�B�4�!�B���1�Q�3��!�A�a�C��#�)�)�B�,��"�M���2�v�
r1g�?c����|z
|z
}�|z
}�|z
}d�	z
}gg�|g||g||gd|z
g|f}|g|g�||z�z
g||g||gd|zg|f}||fSr�r�)
r`rar��ca�cbr�r�r�rbr"s
        ��r/r]z_hyp2f1.<locals>.h�s�����a�C��E���!��2�!�A�#�R�A�a�C�r���a��U�R��G�a��U�Q�q�S�E�2�=���T�A�3��1�Q�3�q�5�	�A�a�5�2�b�'�A�a�C�5�"�D���2�v�
r1)rr��	gammaprodr��epsrRr��absr�rrlr�)rr r!r"rVr`r�rar��ctype�
convergent�finite�zerodiv�abszrWr]rcrbs`  `             @r/rzrz�s��� ��J�Q��
��E��K�J�Q���A�v��V�V�A�a�C��E�]�Q�&�
��)�)�A�,�)�1��6�G�s�y�y��|�/F��Q����)�)�A�,�O�1��6�O��i�i��l�*�q�A�{��{�
M��	�	�!��0L��a��1��/O��
�&�'��=�=�!�Q�q�S��U��a��c�1�Q�3�Z�$�=�G�G��z�z�!�A�a��#�'�'�!�)��,�s�w�w�6�6�
�
��Q��!�q�&��Q�3�J��w�w���y�y��|��Q���I�I�a�L�Q�!�[�q�[��I�I�a�L�Q�!�[�q�[��7�7�N�
��7�7�N��q�6�D��s�{�s�y�y��|��Q��1��:��y�y��|��Q��1��:��s�z�z�!�Q���u� 5��1�a�y�!�N�v�N�N��8�8�D�����B����3�;�
�
��
�
�a�!�A��1�&�1�A���1��X��
�
�
��
�
�a�!�A��1�&�1�A���A�a�C��\�T�
!��
�
�1�a��c�1�a��1��g�.�!�A�#���9�A��s�1�Q�q��4�V�4�A����
�2�I�����s
� B1K�	K$c
�������������� �t��\�}t��\�}t���t���t��}	d}
�D]}�j|�s�|dks�d}
n|	dks|
r	�j��||z��z�fi|��S�dk(rM�j
t��t��z
�}|dkr"�j��dfi|���jzS��fdk(r�t�dz
�dkr��\����\����z�z
}
�j��z
��z
��g��z
��z
d|
g�}d|if�������fd	�	� 	�j� d�jg|jd
�|jdd���}|�j��g���g�zS|	dkr��j��dkr}d�jif������ �fd
�	� |jdd�}	�j� d�jg|jd
�|jdd�|jdd���S����fd�}�j.|��zfi|��S#�j
$r|	dkDs|
r�Y���wxYw#�j
$rY��wxYw#�j
$rd|vr�YnwxYw|jd
�rt!d�	����fd������fd�}�j"dz}�j$}	d�j&z}�xj$dz
c_t)d�D]�}�j+� fd�t)|�D��}�j-� |�jg|||�|jd
�dd��\}}||kr||z}n?|dz}�xj$�j$dzz
c_|dk(s���jd��|�_�S#|�_wxYw) z&
    Evaluates 3F2, 4F3, 5F4, ...
    FrTrg�������?g�������?)rprg�������?c�����z�z
|z}||vr||}n2||dz
}|�|z�z
dz
�|z�z
dz
zz}|||dz
zz}|||<|�	j��|�
�zSr�)r�)r$�_cacher.r�r�r�r�r�r�rr"s    �������r/rz_hypq1fq.<locals>.terms�����2��b���
�A��F�{��1�I���1�Q�3�K���b��d�2�g�a�i�"�Q�$�r�'�!�)�,�,���Q��!��W�����q�	��s�z�z�"�R��!�,�,�,r1r3�strict)r3r�c���t|�}||k7r`��
j|�z�
j|�z}�D]}|�
j||�z}��	D]}|�
j||�z}�|S||vr||S�
|dz
�}|dz
}t	��D]
}|�||zz}�t	��D]
}|�	||zz}�|�z}||z}|||<|Sr�)rr��fac�rfr)�kkr�r$r�r`ra�mr�r r!rr�r�rr"s        �������r/rz_hypq1fq.<locals>.term's�����B��A��B�w�������$�s�w�w�r�{�2���/�A�a�3�6�6�!�B�<�/�a�/��/�A�a�3�6�6�!�B�<�/�a�/����F�{��a�y� ��Q�q�S�	�A��!��A��A�Y�/���c�!�f�Q�h���/��A�Y�/���c�!�f�Q�h���/�
��F�A�
��F�A��F�1�I��Hr1�
sum_methodzr+s+er��)r3r��methodz$Attempting Euler-Maclaurin summationc3�F���K��dgz}t��fd��D��t��fd�|D��z
��j��zz��d�	t���fd��D��t���fd�|D��z
}�dk(r|�j��z
}|���dz
��P�w)Nrc3�F�K�|]}�j|�z����y�wr���loggamma)�.0r`r�k0s  ��r/�	<genexpr>z._hypq1fq.<locals>.log_diffs.<locals>.<genexpr>is�����6�Q�c�l�l�1�R�4�(�6���!c3�F�K�|]}�j|�z����y�wr�r)rrarrs  ��r/rz._hypq1fq.<locals>.log_diffs.<locals>.<genexpr>js�����3�1�C�L�L��2��&�3�rrc3�H�K�|]}�j�|�z����y�wr���psi)rr`rr)rs  ���r/rz._hypq1fq.<locals>.log_diffs.<locals>.<genexpr>ms�����5�A�����!�B�$��5���"c3�H�K�|]}�j�|�z����y�wr�r	)rrarr)rs  ���r/rz._hypq1fq.<locals>.log_diffs.<locals>.<genexpr>ns�����4�A�����!�B�$��4�r)r�log)rr�rcr)r r!rr"s`  @����r/�	log_diffsz_hypq1fq.<locals>.log_diffsgs��������s��B��6�#�6�6��3��3�3�4�68������m�D�
D��A���5��5�5��4��4�4�5����6�������O�A����Q���
�s�BB!c3���K��j�D�cgc]}|��c}�D�cgc]}|��c}�}�j�	|��D]}||z}|���
ycc}wcc}w�wr�)r��	diffs_exp)
rrar`�Cr,rcr r!rrs
      ����r/�hyper_diffsz_hypq1fq.<locals>.hyper_diffsts_������
�
�#�.�Q�q�.�C�0@�q��0@�A�A��]�]�9�R�=�1�
����E����
��/��0@�s�A"�	A�A"�	A
�7A"i�2r;�c3�.�K�|]}�|����y�wr�r�)rr$rs  �r/rz_hypq1fq.<locals>.<genexpr>�s�����?�A��Q��?�s�)r��adiffsr3�error�_fast_abortrr�z*Euler-Maclaurin summation did not convergec
�`��t|d��}t|�d�}g}�j�z}�j�d��}t�dz�D]�}||}|g}|g}	||gzt�dz�D�
cgc]}
|
|k7s�	||
|z
��c}
z}|t��D�
cgc]
}
||
|z
��c}
z}|gt��D�
cgc]
}
|||
z
dz��c}
z}
t�dz�D�
cgc]}
|
|k7s�	d||
z
|z��}}
|j	||	|||
||f���|Scc}
wcc}
wcc}
wcc}
wr�)�listrDr�rEr)�argsr r!�Ts�recz�negzr$�akr�Cpr��Gn�Gd�Fn�Fdrr�r�r"s               ����r/r]z_hypq1fq.<locals>.h�sL����4���8�n���4���8�n��
���w�w�q�y���x�x���x�&���q��s��	5�A��Q��B���A��#��B���t��%��!��*�G�Q��Q��s�1�v�b�y�G�G�B��5��8�4�a��A��r�	�4�4�B���e�A�h�7���C��F��1��7�7�B�',�Q�q�S�z�<�!�Q�!�V�!�C��F�(�2�+�<�B�<��I�I�q�"�b�"�b�"�d�3�4�	5��	��H��4��7��<s$�4
D�?
D�D!�?D&�&
D+�1
D+)rLrr�r�r�r�rrrIrRr��nsumr=r�rD�replacerAr�r�dpsrrN�sumemrl)!rr�r�r r!r"rV�a_types�b_typesr��ispolyr`�Sr.�initialr�rr�r�truncr)�head�tail�errrcr]r�r�r�r�r�rrs!``````                    @@@@@@@r/r~r~�s����
��9�L�C����9�L�C��
�s�)�C�
�s�)�C��q�6�D�
�F�
����9�9�Q�<�A��F��F���
�a�x�6�	��3�:�:�a��G�G�O�S��W�a�J�6�J�J�&	�A�v�
�F�F�3�s�8�C��H�$�%����6��3�9�9�S�#�s�5�f�5����?�?�	�!�u��~�#�a��c�(�T�/����2�b����2��r�E�"�H���-�-��B��r�"�u�R�� 3�R��U�2�b�5��1�4E�F���g�;�		-�		-�	������#�'�'�{�F�J�J�y�4I��z�z�(�D�1��3�A��s�}�}�b��W�b��B�Z�8�8�8��c�z�c�g�g�a�j�A�o��s�w�w�K�	�	�$�Z�Z��g�6�
�	��8�8�D�1�S�W�W�+�v�z�z�)�7L��z�z�(�D�1�!�)�)�#�b�1��3�
3�v� �3�=�=��C��G�.�v�.�.��k� � �	��c�z�V��$�	��Z� � �	��	��@� � �	��*�$���	��
�:�:�i� ��8�9�	�>	�	��g�g��n���x�x��	�����L�E��H�H��N�H��A�Y�
F���x�x�?����?�?���I�I�d�U�C�G�G�,<�#�&�u�-�"�J�J�y�1�� $�	&�&�	��c�
��9��t��A����
�����C�H�H�a�K�'����6��+�+�D�F�F�
F�$�C�H��r�	���C�H�sK�1I�AI.� AJ�I+�*I+�.J�?J�J�J�&C	O�0O�	Oc��������rt��\�}t���ngdc�}�rt��\�}t���ngdc�}|jd�j�|d<	�j||||z��z�fi|��S#�j
$rYnwxYw�j}		|jd�jdz�}
�xjdz
c_d�jif����fd�	��j}td�j�D](}�|�}
||
z
}t|
�|
ks�|c|	�_S	|	�_n#|	�_wxYw||d	zkr�|jd
�}|s~�j��dkr\�tdt���z}�j��dk\rddd
|zd|z�jg}n&ddd|zd|z�jg}nd�jg}|jdi�}����fd�}�j||fddi|��\}}|t|��jzdzkr|S�j
�)Nr��maxterms�	asymp_tolr�r7rc���||vr||S�|dz
�}�D]
}|||dz
zz}��D]
}|||dz
zz}�|�z}||z}|||<|Sr�r�)	r$�cacher�r`rar r!rr"s	     ����r/rz_hyp_borel.<locals>.term�sv����E�z��Q�x���Q�q�S�	�A��(��!��1�Q�3��.�!�(��(��!��1�Q�3��.�!�(�
��F�A�
��F�A��E�!�H��Hr1rrp�contourg�?y@y@@ry��y@��quad_kwargsc�^���j|��j��dgz|�z�zSr�)rvrI)r�r r!rr"s ����r/r�z_hyp_borel.<locals>.g�s/����7�7�A�2�;�s�y�y��c�1�#�g�q��s�;�;�;r1rTr�)rLrr=rr�r�r�rDrr��argrPrR�quad)rr�r�r r!r"rVr)r*rr��sr$r�r7r.r8r��Ir1rs`  ```              @r/rr�sC���
��C�y���W��3�i���2���W�
��C�y���W��3�i���2���W����J����9�F�:��
��s�z�z�!�Q�����S��!�F�v�F�F�����
��
���8�8�D���j�j��c�g�g�a�i�0�����B����S�W�W�+�		�		�
�G�G����3�8�8�$�	�A��Q��A�
��F�A��1�v��}������
	�����4�����A�a�C�x��*�*�Y�'����w�w�q�z�D� ���A�s�1�v��&���7�7�1�:��?� �"�t�Q�h��!��S�W�W�=�G� �#��q�y�!�A�#�s�w�w�?�G�
�c�g�g�,���j�j���3��	<�����!�W�@�D�@�K�@���3��#�a�&����.��"�"��H�

�
�
�s+�"A?�?B�B�!B
E�,E�7E�	E
c	�^���|\\}}\}}|\\}	}
\}}t��}
�j��}�j}|}|jd�xr�j|
�dkD}|ry		�xj|z
c_��fd�}�j	||||	|gdd�jz��}t�fd�|||	|�fD��dk(r�j
|�}||�_S�jd	d	|||
|f|||	|g�fi|��S#�j$rYnwxYw	|�_�=#|�_wxYw)
Nrqrpc���||z|z
|z
}||z}||z}i}�j|d<|dz
|z||zz||zz
|d<d}d}	d}
	|	|vr�d|z
d|zz|dzzd|zz|dzz||zz
||zz||zzd|zd|dzzz
|	zzd|	dzzz}|	|z
|zdz
|	|z
|zdz
z|	|z
dz
z}�j|	z|||	dz
z|||	dz
zz
z||	<||	�|	zz}
t|
�d�jzkrn8|	dkDr&t|
�t|
�zdkr�j�||
z
}|
}
|	dz
}	���j	��|z}�|g|dg||g||gggdf}�g|g||||z
g|||z
||z
g|||z
dz||z
dzg||z
dzgd�zf}�g|g||||z
g|||z
||z
g|||z
dz||z
dzg||zdzgd�zf}|||fS)	Nrrrrp皙�����?r��?r�)rDr�r�r�rv)r�r�r�r��X�A2�B2rb�s1r$�tprev�uu1�uu2�t1r,r�r��T3rr"s                  ��r/r]z_hyp2f2.<locals>.hs�����2��b����A��B��B��B��B��A��7�7�A�a�D��q�D�!�8�B�r�E�>�"�R�%�/�A�a�D��B��A��E���A�:�"#�B�$�q��t�)�B��E�/�!�B�$�"6�r�1�u�"<�R��U�"B�2�b�5�"H��B��"N�PQ�RT�PT�UV�XZ�[\�X\�U]�P]�_`�O`�"`�ab�cd�fg�cg�ag�"g�C�#$�R�4��7�1�9�q��t�B�w�q�y�"9�1�Q�3�q�5�"A�C�#&�7�7�1�9��A�a��c�F�
�3�q��1��v�:�0E�#F�A�a�D��q�T�A���G�^���r�7�S����[�0�!��q�5�S��Z�#�b�'�%9�C�%?�"%�"3�"3�3��b��� "���Q���� ����
�2�
�A��A���1���2�w��2�w�r�"�Q�>�B��"���s�e�R��2�b�5�M�2�b��e�B�r�E�2B�B�r�"�u�Q�w�r�RT�u�UV�w�CW�Y[�\^�Y^�_`�Y`�Xa�bd�ef�bf�f�B��"���s�e�R��2�b�5�M�2�b��e�B�r�E�2B�B�r�"�u�Q�w�r�RT�u�UV�w�CW�Z\�Y\�]_�Y_�`a�Ya�Xb�ce�fg�cg�g�B��r�2�:�%r1Tr��rqr3c3�@�K�|]}�j|����y�wr��r��rr.rs  �r/rz_hyp2f2.<locals>.<genexpr>!s�����E��s�(�(��+�E���rr)	r�rOrr=rlrrr�r�)rr r!r"rVr��a1typer��a2typer��b1typer��b2typer�r�rW�asymp_extraprec�can_use_asymptoticr]rcs`  `                r/r{r{�s[���!$��L�R��,�2�v�!$��L�R��,�2�v��q�6�D��7�7�1�:�D��8�8�D��O�%�j�j��8�8��	������	��
�+	�(
����O�+��&�>�M�M�!�b��B�r�]��PQ�RU�RZ�RZ�PZ�M�[���E�b��B�r�!�_�E�E��J����q�	�A���C�H��3�:�:�a��V�V�V�V�<�r�2�r�2�>N�PQ�\�U[�\�\���$�$�
��
����C�H��t�C�H�s%�/A/D�D�D#�D�D#�#	D,c����|\\}}|\\}}\}	}
t��}�j��}�j}
�xr|dz}|jd�xr-�j|�dkDxr�j	|�d|
zkD}|rw		�xj|z
c_��fd�}�j||||	gdd�jz��}t
�fd	�|||	�fD��dk(r�j|�}||
�_S�jd
d|||
f|||	g�fi|��S#�j$rYnwxYw	|
�_�;#|
�_wxYw)Nrrq�rAc����j||z
|z
�jzz}i}�j|d<d�jd|z|z|zdz
z||z
|z
z||zz�jz
z|d<d||z�j||z
|z
zd|z|z|zdz
zz�jz
dzz�jdd|zdz
z|z|zd||z
|z
zd|dzzd|zz|z|zdz
zzdz
zz|d<d}d}d}d}	||vr�d|dzzd	|zd|zzd|zzdz
|zzd|dzzz||z
dzz
d|z||zdz
zz
�jz}	||z
|z|z
�jz
||z
|z
|z�jz
z||z
|z|z�j
z
z}
�jd|zz|	||dz
z|
||dz
zz
z||<||�d
|zzz}�j|z�jd�|zz|z}�j|z�jd�|zz|z}
t|�d�jzkrn>|dkDr&t|�t|�zd
kr�j�||z
}||
z
}|}|dz
}����j�j|zd�j��zz�|z�j�j|zd�j��zz�|zz}d|z�j�gdd
|g||g|gggdf}�g|g||g||z
||z
g|||z
dz||z
dzggd�zf}||fS)Nrrrpr��r�����i������r@rrA��?)r�rD�mpq_1_4�mpq_3_16�mpq_1_16�mpq_5_2r�r�r�r�r��expjr�r�)r�r�r�rBrbrE�s2r$rFrGrHr%rI�t2r,r�r�rr"s                 ��r/r]z_hyp1f2.<locals>.hIs������R��U�2�X�c�k�k�%9�:�A��A��7�7�A�a�D��c�k�k�1�R�4��7�2�:�a�<�8�"�R�%��(�C�B�r�E�I�#�,�,�V�W�A�a�D��b��e�C�K�K��B��r��$:�A�b�D��G�B�J�q�L�$I�I�#�,�,�V�YZ�Z�Z����c�1�R�4��6�l�2�o�b�&8��2�b�5��8��b��Q��h�r�"�u�n�R�&7��&:�1�&<�=�'>�>?�'@�A�A�A�a�D��B��B��A��E���A�:�#$�Q��T�6�2�b�5��2��:�a��d�?�1�+<�a�*?�#?�!�B��E�'�#I�!#�B���
�$+�-.�r�T�2�b�5��7�^�$<�>A�k�k�$J�C�#$�R�4��7�2�:�c�k�k�#9�A�b�D��G�B�J�s�{�{�<R�"S�!"�2��b����C�K�K�!7�#9�C�#&�7�7�A�a�C�=�#�a��!��f�*�S��1�Q�3��Z�2G�#H�A�a�D��a�D�Q�B�$�q�&�>�1��"�u�u�f�q�[�3�7�7�1�:���+;�;�a�?�� �U�U�A�X�����
�a�R�(8�8�1�<���r�7�S����[�0�!��q�5�S��Z�#�b�'�%9�C�%?�"%�"3�"3�3��b����b��� "���Q���)�*��������!�C�H�H�a�R�L�.�!8�9�"�<����3�6�6�!�8�A�c�h�h��r�l�N�#:�!;�<�R�?�@�A��a�%����!��,�q�$��l�R��H�r�d��B��"�B��"����u�r�"�g�r�"�u�R��U�m��B�r�E�!�G�B�r�E�!�G�,�b�!�A�#�6�B��r�6�Mr1Tr�rKc3�@�K�|]}�j|����y�wr�rMrNs  �r/rz_hyp1f2.<locals>.<genexpr>rs�����B��s�(�(��+�B�rOr�
r�rOrr=r�rlrrr�r�)rr r!r"rVr�rPr�rRr�rSr�r�rWrTrUr]rcs`  `              r/rxrx-sh����M�L�R��!$��L�R��,�2�v��q�6�D��7�7�1�:�D��8�8�D��m�D�!�G�O�%�j�j��8�8�$�	������	�$�	���$��#�d�(�	"���4	�1
����O�+��'"�P�M�M�!�b��B�Z�d�Q�s�x�x�Z�M�X���B�b��B�q�\�B�B�a�G����q�	�A���C�H��3�:�:�a��V�V�V�4�r�2�r�l�A�P��P�P��
�$�$�
��
����C�H��t�C�H�s%�A-D�D/�,D;�.D/�/D;�;	Ec
����|\\}}\}}|\\}	}
\}}\}
}t��}�j��}�xr|dz}�j}|jd�xr-�j|�dkDxr�j	|�d|zkD}|r{		�xj|z
c_��fd�}�j||||	||
gdd�jz��}t
�fd	�|||	||
�fD��d
k(r�j|�}||�_S�jdd|||
||f|||	||
g�fi|��S#�j$rYnwxYw	|�_�?#|�_wxYw)NrrqrWrAc	�8���j||z|z
|z
|z
�jzz}||z}||z|z}||z}||z||zz||zz}	||z|z}
i}�j|d<d|	|z
�jd|z|zdz
z||z
zz�jz
z|d<�j|ddzz�jdd|zdz
z|	|z
zd|
zzdd|dzzd	|zzd
|zz|zdz
z||z
zzdz
zz|d<d}d}
d}d}	||v�r,|d|zz
dz
|d|zz
d|zz
dz
z|d|zz
d|zz
dz
z|d|zz
d|zz
dz
z}d|dz
dzzdd|z|zz|dz
dzzz
dd|dzzd
|z|zzd|	zz|zdz
z|dz
zzd|dzzz
d|z|dzzz
d|	zz
d
|
zz
dd|	z|zdz
z|zz
d|zzdz
}d|dz
dzzdd|z|zd|zz
dzz|dz
zzd|dzz}�jd|zz|||dz
z|||dz
zz
|||dz
zzz||<||�j�d|z�z}�j|z�jd�|zz|z}�j|z�jd�|zz|z}t|�d�jzkrn>|dkDr&t|�t|�zdkr�j�||z
}|
|z
}
|}|dz
}����j�j|zd�j��zz�|z�j�j|zd�j��zz�|
zz}d|z�j�gdd|g|||g||gggdf}�g|g|||||z
g|||z
||z
||z
g|||z
dz||z
dz||z
dzg||z
dzgd�zf}�g|g|||||z
g|||z
||z
||z
g|||z
dz||z
dz||z
dzg||zdzgd�zf}|||fS)NrrrprrY� r�rZr[r���r�ri����r\r@rAr])r�rDr^r_r`rMr�r�r�r�r�rbr�r�)r�r�r�r�r�rBrC�B3�A�B�RrbrErcr$rFrGrH�uu3r%rIrdr,r�r�rJrr"s                          ��r/r]z_hyp2f3.<locals>.h�s�������R��U�2�X�b�[��^�C�K�K�%?�@�A��B��B��B��r��B��2��A��2��b��e��B�r�E�)�A��2��b��A��A��7�7�A�a�D��a�!�e�c�k�k�1�R�4��7�1�9�&=�r�"�u�&E�E����T�U�A�a�D��;�;�q��t�Q�w�.����s�A�b�D��F�|�Q�q�S�?Q�TV�WX�TX�?X��2�b�!�e�8�b��e�+�a��c�1�B�6��:�;�R��U�C�@D�DE�@F�2G�G�A�a�D��B��B��A��E���A�:�#$�Q�q�S�5��7�Q�q��s�U�1�R�4�Z��\�":�A�a��c�E�!�B�$�J�q�L�"I�!"�1�Q�3��q��t��A��#/�C�#$�a��c�A�X�:��1�Q�3�r�6�
�A�a�C�!�8�0C�#C� !�2�a��d�7�2�b�5��7�?�1�Q�3�#6�r�#9�!�#;� <�a��c� B�$C�EG��1��W�$M� "�2��a��d�
�$+�-.�q�S�$1�34�Q�3�$7�9:�A�a�C��F�1�H��a��$H�JK�B�$�$O�OP�$Q�C�$%�a��c�A�X�:�a��Q��r��!�B�$��q��.A�1�Q�3�.G�#G��!�A�$��#N�C�#&�7�7�A�a�C�=�#�a��!��f�*�S��1�Q�3��Z�2G��A�a�PQ�c�F�
�2R�#S�A�a�D��a�D�3�9�9�a�R��a��#8�8��"�u�u�f�q�[�3�7�7�1�:���+;�;�a�?�� �U�U�A�X�����
�a�R�(8�8�1�<���r�7�S����[�0�!��q�5�S��Z�#�b�'�%9�C�%?�"%�"3�"3�3��b����b��� "���Q���)�*��������!�C�H�H�a�R�L�.�!8�9�"�<����3�6�6�!�8�A�c�h�h��r�l�N�#:�!;�<�R�?�@�A��a�%����!��,�q�$��l�R��R�L�2�r�(��B��"�B��"����u�r�"�R��2��&6��2�b�5��B��r�"�u�7M��B�r�E�!�G�B�r�E�!�G�B�r�E�!�G�4�r�"�u�Q�w�i��1��E�B��"����u�r�"�R��2��&6��2�b�5��B��r�"�u�7M��B�r�E�!�G�B�r�E�!�G�B�r�E�!�G�4�r�c�"�f�Q�h�Z��1��E�B��r�2�:�%r1Tr�rKc3�@�K�|]}�j|����y�wr�rMrNs  �r/rz_hyp2f3.<locals>.<genexpr>�s�����H��s�(�(��+�H�rOrjrprf)rr r!r"rVr�rPr�rQr�rRr�rSr��b3typer�r�rTrWrUr]rcs`  `                  r/r|r|s����!$��L�R��,�2�v�/2�,�L�R��,�2�v���V��q�6�D��7�7�1�:�D��m�D�!�G�O��8�8�D�
%�j�j��8�8�=�	������	�=�"%�(�(�4�.�3�t�8�";���:	�7
����O�+��-&�\�M�M�!�b��B�r�"�%5�D�ST�UX�U]�U]�S]�M�^���H�b��B�r�"�Q�5G�H�H�A�M����q�	�A���C�H��3�:�:�a��V�V�V�V�V�D�r�2�r�SU�WY�FZ�\]�h�ag�h�h���$�$�
��
����C�H��t�C�H�s%�A1D-�-D?�<E�>D?�?E�	Ec�D���|\\}}\}}	|j�}	|	jd�j�|	d<�jdd||f||g�fi|	��S#�j$r|jd�r�YnwxYw��fd�}
�j
|
|d|z|z
gfi|��S)Nr3rrrqc�����j|�}d�z}�j||gdd|gg||z
dz|g|g|g|f}�j||gddd|z|z
gg|d|z
g||z
dzgd|z
g|f}||fS)Nr�rr)�sinpir�)r`rar%r�r�r�rr"s      ��r/r]z_hyp2f0.<locals>.h�s�����I�I�a�L��
��T���v�v�a��m�Q�r�!�H�R��1��Q��q�	�1�#�q�c�"�
=�����w�q��n�a��1�Q�3�q�5�\�"�a��!��W�a��c�!�e�W�a��c�U�2�
F���2�v�
r1r)�copyr=rr�r�rl)rr r!r"rVr`r�rar��kwargsbr]s`  `       r/r}r}�s���� ��J�Q��
��E�
��+�+�-��%�k�k�*�c�h�h�?��
���s�z�z�!�Q��u�
��!�u�a�C�7�C�C�����
��:�:�n�%���
����3�=�=��Q��!��A��J�1�&�1�1s�A	A�!A=�<A=Nc�R��������|\}}|\}	}
t|���t|�z�t|	���t|
�z�||z}|	|
z}|D�
cgc]}
�j|
���}}
|D�
cgc]}
�j|
���}}
�j���|�.��krd}��kDrd}��k(r��z�k(rt��dkDrd}nd}|jd�rt	d����|�|dk(r�������fd�}n�������fd�}�j
|||zfi|��Scc}
wcc}
w)Nrrr3zMeijer G m,n,p,q,series =c
����|d�}|�d}g}t��D�]8}�g}||�zg}t��D�cgc]}||k7s�	||||z
��}}|t��D�cgc]}d||z
||z��c}z
}t���D�cgc]
}||||z
��}	}|	t���D�cgc]}d||z
||z��c}z
}	t��D�cgc]}d||z
||z��}
}t��D�cgc]}||k7s�	d||z
||z��}}�
j��z
�z
z��
j�zzz}|j||||	|
||f���;|Scc}wcc}wcc}wcc}wcc}wcc}wr�)rErDr�rr`rarr$�bases�exptsr��gn�gd�hn�hd�hzrr�r+r�r��rr"s             �������r/r]zmeijerg.<locals>.h�s�����R�a��A��Q�R��A��E��1�X�

A������1��a����).�q��<�A�Q�!�V�a��d�1�Q�4�i�<��<��E�!�H�5�q�q��1��v�a��d�{�5�5��).�q���4�A�a��d�1�Q�4�i�4��4��E�!�A�J�7�q�q��1��v�a��d�{�7�7��+0��8�4�a�a��!��f�Q�q�T�k�4��4�+0��8�>�a�q�A�v�a��!��f�Q�q�T�k�>��>��w�w�h�!�A�#�a�%�(�1�s�w�w�q�y�>�9�����e�U�B��B��B�?�@�

A��L��=��5��4��7��4��>s/�
E�
E�E�E!�)E&�E+�4
E0�?E0c
�,��|d�}|�d}g}t��D�]Y}�g}�dk(r
||dz
g}n||dz
�
j��zg}t��D�cgc]}||k7s�	||||z
��}}|t��D�cgc]}d||z
||z��c}z
}t���D�cgc]
}||||z
��}	}|	t���D�cgc]}d||z
||z��c}z
}	t��D�cgc]}d||z
||z��}
}t��D�cgc]}||k7s�	d||z||z
��}}�
j��z
�z
z��
j�zzz}|j||||	|
||f���\|Scc}wcc}wcc}wcc}wcc}wcc}wr�)rErrrDrrzs             �������r/r]zmeijerg.<locals>.hs�����R�a��A��Q�R��A��E��1�X�

A�������6��q�T�!�V�H�E���d�1�f�c�k�k�!�n�4�5�E�).�q��<�A�Q�!�V�a��d�1�Q�4�i�<��<��E�!�H�5�q�q��1��v�a��d�{�5�5��).�q���4�A�a��d�1�Q�4�i�4��4��E�!�A�J�7�q�q��1��v�a��d�{�7�7��+0��8�4�a�a��!��f�Q�q�T�k�4��4�+0��8�>�a�q�A�v�a��!��f�Q�q�T�k�>��>��w�w�h�!�A�#�a�%�(�1�s�w�w�q�y�>�9�����e�U�B��B��B�?�@�

A��L��=��5��4��7��4��>s0�
E8�"
E8�?E=�'F�
F�1F�
F� F)rFrrr�r=rArl)rr r!r"r��seriesrV�an�ap�bm�bqr`ra�_r]r�r+r�r�s`  ``          @@@@r/�meijergr��s8���
�F�B��
�F�B���B��A�	�C��G��A��B��A�	�C��G��A�
�2��A�
�2��A�!"�#�A����Q��#�A�#�!"�#�A����Q��#�A�#����A��A�
�~��q�5�1�&��q�5�1�&���6���s�a�x�C��F�Q�J�����
�z�z�)��
�)�1�Q�q��6�:�
��{�	�	�"	�	�&�3�=�=��A�a�C�*�6�*�*��e	$��#s�D�0D$c	���t|�t|�kDr||}}||}}d�}|j|�rn�|j|�rn}|j|�r||||f\}}}}n`||�sX||z
|dz
z}	||	�std��d|z
|zd|z
||z
|z
zz|j||z
|||z
|z
||	|fi|��zS|j|g|g|gd�d|gi||fi|��S)Nc��t|�dkS)Ng�G�z��?)r�)r*s r/r-zappellf1.<locals>.ok+s���1�v��}�r1rz%Analytic continuation not implemented��m+nr�r+r�)r�rr?�appellf1�hyper2d)
rr`r�r�rbr*�yrVr-�u1s
          r/r�r�$s���1�v��A����!�1���R�B����{�{�1�~��	���R���	���R���!�R��|���1�b�"�
�!�u��A�#��!���B��b�6� �!H�I�I��a�C�B�3�<��1���!��B���/�����Q�q�S��A�b�D��G�A�b��<�V�<�=�
=��3�;�;�q�c�r�d��t�4�u�a�S�k�1�Q�Q�&�Q�Qr1c�F�|j|g|g|gd�|g|gd�||fi|��S)Nr��r�r+�r�)	rr`r�r��c1�c2r*r�rVs	         r/�appellf2r�As>���3�;�;�q�c�r�d��t�4��T�r�d��Q�q�,�$*�,�,r1c�&�|j|�xs|j|�}	|j|�xs|j|�}
|	s(|
st|�t|�kDr||}}||||f\}}}}|j||g||gd�d|gi||fi|��S)Nr�r�)rr�r�)rr�r�r�r�rbr*r�rV�outer_polynomial�inner_polynomials           r/�appellf3r�Gs����{�{�2��9�#�+�+�b�/���{�{�2��9�#�+�+�b�/����s�1�v��A����a�q�A��R��2�+�K�B�r�"�R��3�;�;�R��G��B��0�5�!��+�a��K�F�K�Kr1c�@�|jd||gi|g|gd�||fi|��S)Nr�r�r�)rr`rar�r�r*r�rVs        r/�appellf4r�Qs1���3�;�;��q��e�}�B�4�R�D�&9�!�A�G��G�Gr1c�����j|�}�j|�}�fd�}t|�}t|�}||d�}	||d�}
||d�}||d�}||d�}
||d�}||d�}||d	�}||d�}||d�}||d�}|rtd
|j�dz��|rtd
|j�dz��d}�j}�jd�}d}�j}|jdd
|z�}	�xjdz
c_�j�}	d}d}t|
�}t|�}|	D�cgc]}||z��	}}|D�cgc]}||z��	} }|D])}||z}|j|�|j|��+|D])}||z}|j|�| j|��+|
D]-}|j||z
�| j||z
dz
��/|D];}|dz}|d|zz}|jd|z
|z
�|j||z
��=|D]<}|j|d|zz�|j|d|zzd|zd|zzz��>|D]D}|dz}|jd|z
d|zz
�|j|d|zzd|zd|zzz��F|D]O}|dz}|jd||z
z�|jd||z
dzz�| j||z
dz
��Q�j||||zfd�ji|��}!||!z|z}"t|"�|kr|dz
}nd}|dk\s|snB||"z
}|D]}||z}�	| D]}||z}�	|dz
}||z|z}||kDr�jd����v	|�_|�Scc}wcc}w#|�_wxYw)az
    Sums the generalized 2D hypergeometric series

    .. math ::

        \sum_{m=0}^{\infty} \sum_{n=0}^{\infty}
            \frac{P((a),m,n)}{Q((b),m,n)}
            \frac{x^m y^n} {m! n!}

    where `(a) = (a_1,\ldots,a_r)`, `(b) = (b_1,\ldots,b_s)` and where
    `P` and `Q` are products of rising factorials such as `(a_j)_n` or
    `(a_j)_{m+n}`. `P` and `Q` are specified in the form of dicts, with
    the `m` and `n` dependence as keys and parameter lists as values.
    The supported rising factorials are given in the following table
    (note that only a few are supported in `Q`):

    +------------+-------------------+--------+
    | Key        |  Rising factorial | `Q`    |
    +============+===================+========+
    | ``'m'``    |   `(a_j)_m`       | Yes    |
    +------------+-------------------+--------+
    | ``'n'``    |   `(a_j)_n`       | Yes    |
    +------------+-------------------+--------+
    | ``'m+n'``  |   `(a_j)_{m+n}`   | Yes    |
    +------------+-------------------+--------+
    | ``'m-n'``  |   `(a_j)_{m-n}`   | No     |
    +------------+-------------------+--------+
    | ``'n-m'``  |   `(a_j)_{n-m}`   | No     |
    +------------+-------------------+--------+
    | ``'2m+n'`` |   `(a_j)_{2m+n}`  | No     |
    +------------+-------------------+--------+
    | ``'2m-n'`` |   `(a_j)_{2m-n}`  | No     |
    +------------+-------------------+--------+
    | ``'2n-m'`` |   `(a_j)_{2n-m}`  | No     |
    +------------+-------------------+--------+

    For example, the Appell F1 and F4 functions

    .. math ::

        F_1 = \sum_{m=0}^{\infty} \sum_{n=0}^{\infty}
              \frac{(a)_{m+n} (b)_m (c)_n}{(d)_{m+n}}
              \frac{x^m y^n}{m! n!}

        F_4 = \sum_{m=0}^{\infty} \sum_{n=0}^{\infty}
              \frac{(a)_{m+n} (b)_{m+n}}{(c)_m (d)_{n}}
              \frac{x^m y^n}{m! n!}

    can be represented respectively as

        ``hyper2d({'m+n':[a], 'm':[b], 'n':[c]}, {'m+n':[d]}, x, y)``

        ``hyper2d({'m+n':[a,b]}, {'m':[c], 'n':[d]}, x, y)``

    More generally, :func:`~mpmath.hyper2d` can evaluate any of the 34 distinct
    convergent second-order (generalized Gaussian) hypergeometric
    series enumerated by Horn, as well as the Kampe de Feriet
    function.

    The series is computed by rewriting it so that the inner
    series (i.e. the series containing `n` and `y`) has the form of an
    ordinary generalized hypergeometric series and thereby can be
    evaluated efficiently using :func:`~mpmath.hyper`. If possible,
    manually swapping `x` and `y` and the corresponding parameters
    can sometimes give better results.

    **Examples**

    Two separable cases: a product of two geometric series, and a
    product of two Gaussian hypergeometric functions::

        >>> from mpmath import *
        >>> mp.dps = 25; mp.pretty = True
        >>> x, y = mpf(0.25), mpf(0.5)
        >>> hyper2d({'m':1,'n':1}, {}, x,y)
        2.666666666666666666666667
        >>> 1/(1-x)/(1-y)
        2.666666666666666666666667
        >>> hyper2d({'m':[1,2],'n':[3,4]}, {'m':[5],'n':[6]}, x,y)
        4.164358531238938319669856
        >>> hyp2f1(1,2,5,x)*hyp2f1(3,4,6,y)
        4.164358531238938319669856

    Some more series that can be done in closed form::

        >>> hyper2d({'m':1,'n':1},{'m+n':1},x,y)
        2.013417124712514809623881
        >>> (exp(x)*x-exp(y)*y)/(x-y)
        2.013417124712514809623881

    Six of the 34 Horn functions, G1-G3 and H1-H3::

        >>> from mpmath import *
        >>> mp.dps = 10; mp.pretty = True
        >>> x, y = 0.0625, 0.125
        >>> a1,a2,b1,b2,c1,c2,d = 1.1,-1.2,-1.3,-1.4,1.5,-1.6,1.7
        >>> hyper2d({'m+n':a1,'n-m':b1,'m-n':b2},{},x,y)  # G1
        1.139090746
        >>> nsum(lambda m,n: rf(a1,m+n)*rf(b1,n-m)*rf(b2,m-n)*\
        ...     x**m*y**n/fac(m)/fac(n), [0,inf], [0,inf])
        1.139090746
        >>> hyper2d({'m':a1,'n':a2,'n-m':b1,'m-n':b2},{},x,y)  # G2
        0.9503682696
        >>> nsum(lambda m,n: rf(a1,m)*rf(a2,n)*rf(b1,n-m)*rf(b2,m-n)*\
        ...     x**m*y**n/fac(m)/fac(n), [0,inf], [0,inf])
        0.9503682696
        >>> hyper2d({'2n-m':a1,'2m-n':a2},{},x,y)  # G3
        1.029372029
        >>> nsum(lambda m,n: rf(a1,2*n-m)*rf(a2,2*m-n)*\
        ...     x**m*y**n/fac(m)/fac(n), [0,inf], [0,inf])
        1.029372029
        >>> hyper2d({'m-n':a1,'m+n':b1,'n':c1},{'m':d},x,y)  # H1
        -1.605331256
        >>> nsum(lambda m,n: rf(a1,m-n)*rf(b1,m+n)*rf(c1,n)/rf(d,m)*\
        ...     x**m*y**n/fac(m)/fac(n), [0,inf], [0,inf])
        -1.605331256
        >>> hyper2d({'m-n':a1,'m':b1,'n':[c1,c2]},{'m':d},x,y)  # H2
        -2.35405404
        >>> nsum(lambda m,n: rf(a1,m-n)*rf(b1,m)*rf(c1,n)*rf(c2,n)/rf(d,m)*\
        ...     x**m*y**n/fac(m)/fac(n), [0,inf], [0,inf])
        -2.35405404
        >>> hyper2d({'2m+n':a1,'n':b1},{'m+n':c1},x,y)  # H3
        0.974479074
        >>> nsum(lambda m,n: rf(a1,2*m+n)*rf(b1,n)/rf(c1,m+n)*\
        ...     x**m*y**n/fac(m)/fac(n), [0,inf], [0,inf])
        0.974479074

    **References**

    1. [SrivastavaKarlsson]_
    2. [Weisstein]_ http://mathworld.wolfram.com/HornFunction.html
    3. [Weisstein]_ http://mathworld.wolfram.com/AppellHypergeometricFunction.html

    c���|j|g�}	t|�}|D�cgc]}�j|���c}S#t$r|g}Y�-wxYwcc}wr�)�popr�	TypeErrorrr)�dct�keyrr:rs    �r/�parsezhyper2d.<locals>.parse�sW����w�w�s�B���	���:�D�-1�1�S����C� �1�1���	��6�D�	��1s�?�A�A�
Ar�r+r�zm-nzn-mz2m+nz2m-nz2n-mzunsupported key: %rrr3r;r7rr�rr�r]r5rpzmaxterms exceeded in hyper2d)rr�dictr?�keysrDr�rr=r�rrrIr�r�)#rr`rar*r�rVr�r r!�a_m�a_n�	a_m_add_n�	a_m_sub_n�	a_n_sub_m�
a_2m_add_n�
a_2m_sub_n�
a_2n_sub_m�b_m�b_n�	b_m_add_nr<�outerr��ok_countrr3r��
inner_sign�
outer_sign�inner_a�inner_b�outer_a�outer_b�innerrs#`                                  r/r�r�Vsj���P	���A��A����A��A�2��q�'�C�
�q�'�C�
��3�-�C�
��3�-�C��a���I��a���I��a���I��q�&�!�J��q�&�!�J��q�&�!�J�
��3�-�C�
��3�-�C��a���I��
�0�1�6�6�8�A�;�>�?�
?��
�0�1�6�6�8�A�;�>�?�
?�	�A��G�G�E�����
�A��H��8�8�D��z�z�*�b��g�.�H�>����B����w�w�h����J��J��3�i�G��3�i�G�$'�(�q�q��s�(�G�(�$'�(�q�q��s�(�G�(��
"���a�C�����q�!����q�!�
"�
�
"���a�C�����q�!����q�!�
"�
�
&�����q��s�#����q��s�1�u�%�
&��
%���r�"�
��r�Q�i�'�
����q��s�1�u�%�����r�!�t�$�	
%� �
2�����q��1��u�%�����!�A�#���!��A�a�C��0�1�
2� �
2���r�"�
����q��s�1�Q�3�w�'�����!�A�#���!��A�a�C��0�1�
2�
 �
&���a��
����s�A�a�C�y�)����s�A�a�C��E�{�+����q��s�1�u�%�	
&�
�C�I�I�g�w�
�1��-����-�%+�-�E��5�=�:�-�D��4�y�3���A�
�����1�}�E��
��I�A��(��e�q�j�e�(��(��e�q�j�e�(�
��F�A��A�I��M�E��8�|��'�'�(F�G�G�s�d����
�2�I��o)��(��j���s,�'AO�(O�4O�:O�H:O�
O�	Oc�����	�|j���||z}t|��	t|�}�	|fdk(s�	|fdk(r|j�zS�	|z
dz���	�fd�}|j||fi|��S)a`
    Evaluates the bilateral hypergeometric series

    .. math ::

        \,_AH_B(a_1, \ldots, a_k; b_1, \ldots, b_B; z) =
            \sum_{n=-\infty}^{\infty}
            \frac{(a_1)_n \ldots (a_A)_n}
                 {(b_1)_n \ldots (b_B)_n} \, z^n

    where, for direct convergence, `A = B` and `|z| = 1`, although a
    regularized sum exists more generally by considering the
    bilateral series as a sum of two ordinary hypergeometric
    functions. In order for the series to make sense, none of the
    parameters may be integers.

    **Examples**

    The value of `\,_2H_2` at `z = 1` is given by Dougall's formula::

        >>> from mpmath import *
        >>> mp.dps = 25; mp.pretty = True
        >>> a,b,c,d = 0.5, 1.5, 2.25, 3.25
        >>> bihyper([a,b],[c,d],1)
        -14.49118026212345786148847
        >>> gammaprod([c,d,1-a,1-b,c+d-a-b-1],[c-a,d-a,c-b,d-b])
        -14.49118026212345786148847

    The regularized function `\,_1H_0` can be expressed as the
    sum of one `\,_2F_0` function and one `\,_1F_1` function::

        >>> a = mpf(0.25)
        >>> z = mpf(0.75)
        >>> bihyper([a], [], z)
        (0.2454393389657273841385582 + 0.2454393389657273841385582j)
        >>> hyper([a,1],[],z) + (hyper([1],[1-a],-1/z)-1)
        (0.2454393389657273841385582 + 0.2454393389657273841385582j)
        >>> hyper([a,1],[],z) + hyper([1],[2-a],-1/z)/z/(a-1)
        (0.2454393389657273841385582 + 0.2454393389657273841385582j)

    **References**

    1. [Slater]_ (chapter 6: "Bilateral Series", pp. 180-189)
    2. [Wikipedia]_ http://en.wikipedia.org/wiki/Bilateral_hypergeometric_series

    )rr)rrrc���t|d��}t|�d�}|D�cgc]}d|z
��	}}|D�cgc]}d|z
��	}}d�z�
zg|D�cgc]}d|z
��	c}z|D�cgc]}d|z
��	c}z}dgdgt|�zzdgt|�zz}gggg|dgz|�
f}	||gg|dgz|d�z�
zf}
|	|
fScc}wcc}wcc}wcc}w)Nrr�r)rrF)rr r!ra�aa_sr`�bb_s�rp�rcr�r��negr�r"s           ���r/r]zbihyper.<locals>.h{s�����3�r��7�m���3�q�r�7�m�� �!���!��!��!� �!���!��!��!��C�i�!�m�_�S�1���!��1�
1�#�4F�Q�Q�q�S�4F�
F���T�Q�C��C��L�
 �B�4��C��=�
0��
��R��S�A�3�Y��Q�
.��
��R��T�Q�C�Z���S�y�1�}�
<���2�v�
��
"��!��1��4Fs�B:�B?�C�"C	)rrrFr<rl)
rr r!r"rVrr�r]r�r�s
   `    @@r/�bihyperr�Ds}���`	���A��A�

��)�C��C��A��C��A�	�1�v���1�a�&�E�/��x�x�!�|���Q�3�!�)�C�	��3�=�=��C�*�6�*�*r1)rN)#�
libmp.backendr�	functionsrrr0r@rlrIr�r�r�r�r�r�r�r�ryrurwr�rzr~rr{rxr|r}r�r�r�r�r�r�r�r�r1r/�<module>r�sl��"�+�-2�^���$&�D�E��E�N�(8��(8�T�(��(��)��)��.��.��+��+��1��1��4��4��*��*��4��4������!8��!8�F�!
��!
�F4�l�M��M�^�G/��G/�R�;��;�|�@]��@]�H�MQ��MQ�b�Oi��Oi�b�2��2�(�;+��;+�z�R��R�8�,��,�
�L��L��H��H��d��d�L��@+��@+r1

Youez - 2016 - github.com/yon3zu
LinuXploit