403Webshell
Server IP : 217.160.0.135  /  Your IP : 216.73.217.37
Web Server : Apache
System : Linux www 6.18.52-i1-ampere #1203 SMP Mon Sep 14 18:29:59 CEST 2026 aarch64
User : sws1074145052 ( 1074145052)
PHP Version : 8.3.32
Disable Function : NONE
MySQL : OFF  |  cURL : ON  |  WGET : ON  |  Perl : ON  |  Python : OFF  |  Sudo : OFF  |  Pkexec : OFF
Directory :  /lib/python3/dist-packages/scipy/special/tests/__pycache__/

Upload File :
current_dir [ Writeable ] document_root [ Writeable ]

 

Command :


[ Back ]     

Current File : /lib/python3/dist-packages/scipy/special/tests/__pycache__/test_mpmath.cpython-311.pyc
�

d�c�%��X�dZddlZddlmZmZddlmZddlZddlZddl	m
Z
ddlmZ
ddlmZmZmZmZddlmZmZmZmZmZmZmZmZmZmZddlmZm Z m!Z!m"Z"m#Z#m$Z$	ddl%Z%n#e&$red	��Z%YnwxYwee%d
��d���Z'ee%d��d
���Z(ee%d��d���Z)ee%d��d���Z*ee%d��d���Z+ee%d��d���Z,ee%d��d���Z-ee%d��d���Z.ee%d��d���Z/ee%d��ej0j1d�����Z2ee%d��d���Z3ee%d��d���Z4ee%d��d���Z5dZ6ee%d��d ���Z7ee%d��d!���Z8ee%d��ej0j1d"�����Z9ee%d��ej0j1d#�����Z:ee%d��d$���Z;ee%d��d%���Z<ee%d��ej0j1d&�����Z=ee%d��ej0j1d'�����Z>ee%d��d(���Z?ee%d��d)���Z@ee%d��d*���ZAd+�ZBej0j1ee%d��d,�����ZCej0j1ee%d��d-�����ZDej0j1ee%d��d.�����ZEej0j1ee%d��d/�����ZFeGd0d0�1��ZHej0j1ee%d2��Gd3�d4������ZIdS)5z4
Test SciPy functions versus mpmath, if available.

�N)�assert_�assert_allclose)�pi)�_pep440)�
MissingModule�
check_version�FuncData�assert_func_equal)
�Arg�FixedArg�
ComplexArg�IntArg�assert_mpmath_equal�nonfunctional_tooslow�
trace_args�time_limited�exception_to_nan�
inf_to_nan)�_sinpi�_cospi�_lgam1p�_lanczos_sum_expg_scaled�_log1pmx�	_igam_fac�mpmathz0.10c
���g}tjddd��D]x}tjddtjzd��D]R}|tjd|z��z}|�|t
tj|����f���S�ytj	|tj
���}ttj
|dd�����dS)	Ni������
r���?��dtype�)�np�logspace�linspacer�exp�append�complexr�ei�array�complex_r	�sc�expi�check)�dataset�r�p�zs    �A/usr/lib/python3/dist-packages/scipy/special/tests/test_mpmath.py�test_expi_complexr5#s����G�
�[��a��
$�
$�7�7����Q��"�%���,�,�	7�	7�A��"�&��A��,�,��A��N�N�A�w�v�y��|�|�4�4�5�6�6�6�6�	7��h�w�b�k�2�2�2�G��R�W�g�q�!�$�$�*�*�,�,�,�,�,�z0.19c
��g}dD]}tjddd��D]f}tjd��5|�||ttj||����f��ddd��n#1swxYwY�g��tj|��}ttj
|ddd�	�����dS)
N)�2�3r����d�rr#r�vIh�%<=��rtol)r$r%r�workdpsr(�float�expint�asarrayr	r-�expnr/)r0�n�xs   r4�test_expn_large_nrH3s!���G�
�C�C����Q��3�'�'�	C�	C�A����$�$�
C�
C�����1�e�F�M�!�Q�,?�,?�&@�&@�A�B�B�B�
C�
C�
C�
C�
C�
C�
C�
C�
C�
C�
C����
C�
C�
C�
C��	C��j��!�!�G��R�W�g�v�q�u�5�5�5�;�;�=�=�=�=�=s�9A8�8A<�?A<c
��g}gd�}|D]z}|D]u}|D]p}|d|zz}tjd��5ttj||����}ddd��n#1swxYwY|�|||f���q�v�{tj|��}td�|ddd������dS)	N)g�X�g#�����?g#@g�X@r �xc�6�tj|j|��S�N�r-�hyp0f1�real��vr3s  r4�<lambda>z$test_hyp0f1_gh5764.<locals>.<lambda>Us��"�)�A�F�A�.�.�r6r=rr>r?)	rrAr)rPr(r$r+r	r/)r0�axisrSrG�yr3�ress       r4�test_hyp0f1_gh5764rXDs4���G�.�.�.�D�
�	,�	,���	,�	,�A��
,�
,����1��H���^�C�(�(�7�7�!�&�-��1�"5�"5�6�6�C�7�7�7�7�7�7�7�7�7�7�7����7�7�7�7�����1�c�{�+�+�+�+�
,�	,��h�w���G��
.�
.��������������s�#A	�A"
�%A"
c���tjddd��}tj|d��}tjd�|D����}t||�t��d���dS)N���rKc�8�g|]}tj|d����S)rK�rrP)�.0rSs  r4�
<listcomp>z'test_hyp0f1_gh_1609.<locals>.<listcomp>^s$��5�5�5�Q�6�=��C�(�(�5�5�5r6��-���q=r?)r$r&r-rPr+r�astyperB)�vv�af�mfs   r4�test_hyp0f1_gh_1609rfYsm��
��S�#�r�	"�	"�B�	��2�s�	�	�B�	��5�5�"�5�5�5�	6�	6�B��B��	�	�%�(�(�u�5�5�5�5�5�5r6z1.1.0c�~�g}tjdd��D]�}tjddd��D]�}|}|�||dt	tj||d����f��||zdz
}|�||dt	tj||d����f������dD]A}dD]<}|�||dt	tj||d����f���=�Btj|��}ttj|dd	d
d����
��dS)
N�����rr#)�%����rJrrKr#r)��rJrrKr#��?r�@�rr#r��V瞯�<��-���a=�r@�atol)r$�aranger&r(rBr�hyperur+r	r-r/)r0rF�b�as    r4�test_hyperu_around_0rzfs`���G�
�Y�r�1�
�
�E�E����R��B�'�'�	E�	E�A���A��N�N�A�q�!�U�6�=��A�q�+A�+A�%B�%B�C�D�D�D���Q���
�A��N�N�A�q�!�U�6�=��A�q�+A�+A�%B�%B�C�D�D�D�D�		E�0�E�E��5�	E�	E�A��N�N�A�q�!�U�6�=��A�q�+A�+A�%B�%B�C�D�D�D�D�	E��h�w���G��R�Y���A�E��F�F�F�L�L�N�N�N�N�Nr6�1.0.0c�T��ddg}|ttjgd�d�����z
}d�|D��}td�����fd	�|D��}t	j|tj�
��}ttj	|ddd�
���
��dS)N)r���r}�ffffff�?)r���rr~)rr#gffffff�����r:)�repeatc�n�g|]2\}}}}||k�
t|��|k� |dk�&|dk�,||||f��3S)rr�)�round)r_ryrx�crGs     r4r`z.test_hyp2f1_strange_points.<locals>.<listcomp>�sW�����#��A�q�!���6�6�e�A�h�h�!�m�m��A���!�u�*�*�
�A�q�!��2<�*�*r6T)�	eliminatec
�T��g|]$}|ttj|i�����fz��%S��rBr�hyp2f1)r_r2�kws  �r4r`z.test_hyp2f1_strange_points.<locals>.<listcomp>�s8���B�B�B��q�E�&�-��1�b�1�1�2�2�4�4�B�B�Br6r!�rr#rrq绽���|�=r?)�list�	itertools�product�dictr$r+�float_r	r-r�r/)�ptsr0r�s  @r4�test_hyp2f1_strange_pointsr�}s����	���C��4�	�!�"5�"5�"5�a�@�@�@�A�A�A�C���'*����C�
��	�	�	�B�B�B�B�B�c�B�B�B�G��h�w�b�i�0�0�0�G��R�Y���A�E�:�:�:�@�@�B�B�B�B�Br6z0.13c�"�gd�}d�|D��}tj|tj���}tjd���5t	t
j|ddd�	�����ddd��dS#1swxYwYdS)
N)�r#rrqr)�UUUUUU�?gUUUUUU�?g�������?g�?)��?rKg�?gX�<ݚ�?)rr���rq)rr�rrq)rrlrlrqr�)g��x�&�?r}rhg333333�?)r�r�r�+����?)r�r�rr})rrqri�G�z��?)rnrJrqr�)rrog
�r�)���g*4�]92@g|��O�%@gC��u��?)����rkr�)r�r��%@r�)r}rr#��?)r}rr#rm)r��
rir�)r�r�rirm)rKg�p�rngm˟��?c�L�g|]!}|ttj|���fz��"Sr�r�)r_r2s  r4r`z0test_hyp2f1_real_some_points.<locals>.<listcomp>�s/��<�<�<�1�q�E�&�-��+�,�,�.�.�<�<�<r6r!�ignore��invalidr�r:r�r?)r$r+r��errstater	r-r�r/)r�r0s  r4�test_hyp2f1_real_some_pointsr��s������C�.=�<��<�<�<�G��h�w�b�i�0�0�0�G�	��X�	&�	&�	&�G�G����G�Y���>�>�>�D�D�F�F�F�G�G�G�G�G�G�G�G�G�G�G�G����G�G�G�G�G�Gs�1B�B�B�0.14c����gd�}d���fd�|D��}tj|tj���}ttj|ddd������dS)	N))�p)r9r)���rgwJ��)r�|���r�r�)rr�rkr�c�l�t|t��rt|d��|dzS|S)Nrr#)�
isinstance�tuplerB�rGs r4�fevz&test_hyp2f1_some_points_2.<locals>.fev�s2���a����	���1��;�;��1��%�%��Hr6c���g|]<}tt�|����ttj|���fz��=Sr�)r��maprBrr�)r_r2r�s  �r4r`z-test_hyp2f1_some_points_2.<locals>.<listcomp>�s@���M�M�M�A�u�S��a�[�[�!�!�U�6�=�!�+<�%=�%=�$?�?�M�M�Mr6r!r�r:r�r?)r$r+r�r	r-r�r/)r�r0r�s  @r4�test_hyp2f1_some_points_2r��s�������C����N�M�M�M��M�M�M�G��h�w�b�i�0�0�0�G��R�Y���A�E�:�:�:�@�@�B�B�B�B�Br6c���g}dD]a}dD]\}dD]W}dD]R}	ttj||||����}n#t$rY�4wxYw|�|||||f���S�X�]�btj|t
j���}tjd���5ttj|dd	d
d������ddd��dS#1swxYwYdS)
N)r�rh��������g������?rir)g�r}r#g������@)r�r�rigffffff4@)r�g)\��(��G�z��rg333333�?gffffff�?rnrr!r�r�r�r:��&�.>T�r@�ignore_inf_sign)rBrr��	Exceptionr(r$r+r�r�r	r-r/)r0ryrxr�r3rSs      r4�test_hyp2f1_real_somer��sp���G�
(�4�4��#�	4�	4�A�(�
4�
4��C�4�4�A�!�!�&�-��1�a��";�";�<�<����$�!�!�!� ��!�����N�N�A�q�!�Q��?�3�3�3�3�4�
4�	4��h�w�b�i�0�0�0�G�	��X�	&�	&�	&�/�/����G�Y���!%�	'�	'�	'�',�u�w�w�w�/�/�/�/�/�/�/�/�/�/�/�/����/�/�/�/�/�/s!�$9�
A�A�2C�C�"Cz0.12c
���d}tj|dftj��}tj�d��tj�d|��|dd�df<tj�d|��|dd�df<tj�d|��|dd�df<dtj�|��zdz
|dd�df<|dd�dfxxd	tj�d|��zzcc<|dd�dfxxd	tj�d|��zzcc<|dd�dfxxd	tj�d|��zzcc<|D]�}tj	d
krlt|dd������t|d��kr,t|dd������|d<ttj
t|dd������|d<��tt j
|ddd
������dS)N��ri��rnrr#rrqr}r�r:r�r�r?)r$�zerosr��random�seed�pareto�rand�randintr�__version__�abs�maxrBr�r�r	r-r/)�npointsr0�dss   r4�test_hyp2f1_real_randomr��s%���G��h���|�R�Y�/�/�G��I�N�N�4�����I�$�$�S�'�2�2�G�A�A�A�q�D�M��I�$�$�S�'�2�2�G�A�A�A�q�D�M��I�$�$�S�'�2�2�G�A�A�A�q�D�M��b�i�n�n�W�-�-�-��1�G�A�A�A�q�D�M��A�A�A�q�D�M�M�M�b�2�9�,�,�Q��8�8�8�8�M�M�M��A�A�A�q�D�M�M�M�b�2�9�,�,�Q��8�8�8�8�M�M�M��A�A�A�q�D�M�M�M�b�2�9�,�,�Q��8�8�8�8�M�M�M��5�5�����&�&��2�b�q�b�6�{�{��� � �3�r�!�u�:�:�-�-��B�r��r�F�����)�)��1���f�m�U�2�b�q�b�6�]�]�3�4�4��1����R�Y���q�t�<�<�<�B�B�D�D�D�D�Dr6c�h�tjjtjj}}	dtj_t	jt	jddd��t	jddd����\}}t	jt	jddd��t	jddd����\}}tj|�	��|�	��fdtj|�	��|�	��fzz}ttjd	�|d
d���ttj
d
�|d
d���||ctj_tj_dS#||ctj_tj_wxYw)N�Fr�r#��i����g�������?r c�D�ttj|����SrN)r)r�erfr�s r4rTz"test_erf_complex.<locals>.<lambda>s��G�F�J�q�M�M�,B�,B�r6Fr>��
vectorizedr@c�D�ttj|����SrN)r)r�erfcr�s r4rTz"test_erf_complex.<locals>.<lambda>s��W�V�[��^�^�-D�-D�r6)r�mp�dps�precr$�meshgridr&r%�r_�ravelr
r-r�r�)�old_dps�old_prec�x1�y1�x2�y2�pointss       r4�test_erf_complexr��sV���	�
�v�y�~�X�G�:���	�
���R�[��a��4�4�b�k�#�q�"�6M�6M�N�N���B���R�[��b�"�5�5�r�{�3��B�7O�7O�P�P���B���r�x�x�z�z�"�(�(�*�*�,�-��2�5������R�X�X�Z�Z�9O�3P�0P�P���"�&�"B�"B�F�%*��	8�	8�	8�	8��"�'�#D�#D�f�%*��	8�	8�	8�	8�)0��%��	�
�v�y�~�~�~����%��	�
�v�y�~�9�9�9�9s�EF�%F1z0.15c�"��g}dD]}|�dd|fdd|fdd|fdd|fdd|fdd|fdd|fdd|fdd|fdd	|fd
d|fdd	|fd
d|fdd	|fd
d|fdd|fd
d|fdd|fd
d|fdd|fd
d|fdd|fd
d|fdd|fd
d|fdd|fg����d���fd�|D��}tj|tj���}d�}tjd���5t||dd
dd������ddd��dS#1swxYwYdS)N)r�gm�������ư>gL7�A`�?r#r#r}rg333333�?g333333��g�������rr�rq�rir�rhr�i����g������ @g������&@g������ �g������&�c�t�|t|��kr|dkr
|dkrdSdStj|||��S�Nr#r)�intr�legenp)�nu�murGs   r4�mplegenpztest_lpmv.<locals>.mplegenp-s@��
��R���=�=�Q�!�V�V��Q�w�w��q��q��}�R��Q�'�'�'r6c	�V��g|]%}|�|d|d|d��fz��&S)r#rrr�)r_r2r�s  �r4r`ztest_lpmv.<locals>.<listcomp>6s;���>�>�>�Q�q�H�H�Q�q�T�1�Q�4��1��.�.�0�0�>�>�>r6r!c�^�tj|�t��||��SrN)r-�lpmvrbr�)r�r�rGs   r4�evfztest_lpmv.<locals>.evf9s ���w�r�y�y��~�~�r�1�-�-�-r6r�r�rpr��+����=rt)�extendr$r+r�r�r	r/)r�rGr0r�r�s    @r4�	test_lpmvr�s$���
�C�
,�����
�
�
��1�I�
��A�J�
��A�J�
��Q�K�
��Q�K�
��a�L�
��a�L�
��q�M�
��A�J�
��A�J�
��1�I�
��A�J�
��A�J�
��Q�K�
��A�J�
��Q�K�
��Q�K�
��a�L�
��Q�K�
��a�L�
��a�L�
��q�M�
��a�L�
��q�M�
��q�M�
���N�5
�	�	�	�	�:(�(�(�?�>�>�>�#�>�>�>�G��h�w�b�i�0�0�0�G�.�.�.�
��X�	&�	&�	&�K�K���g�w���E�B�B�B�H�H�J�J�J�K�K�K�K�K�K�K�K�K�K�K�K����K�K�K�K�K�Ks�(D�D�Dc�r�tj�d��tjtjddd��tjddd��tjddd��tjdd	d��tjdd	d��d
zdddd
df
}|}tjtj|dd�df|ddd�f�����dd��j	}tjjtjj
}}	dtj_ttjd�|ddd���ttjd�|dd���||ctj_tj_
dS#||ctj_tj_
wxYw)Nr��8���r;r:r�rr}r#r�rKgffffff�r�g33333Y�g3333���r�c�F�ttj||����SrN)rBr�beta�ryrxs  r4rTztest_beta.<locals>.<lambda>Ws��u�V�[��A�->�->�'?�'?�r6Fr�T)r�r@r�c
��ttjttj||��������SrN)rBr�logr�r�r�s  r4rTztest_beta.<locals>.<lambda>_s,���v�z�#�f�k�!�Q�.?�.?�*@�*@�A�A�B�B�r6r�)r$r�r�r�r%rvr+�broadcast_arrays�reshape�Trr�r�r�r
r-r��betaln)rxry�abr�r�s     r4�	test_betarDs����I�N�N�4����
��b�k�$��Q�'�'��k�#�r�1�%�%��k�"�a��#�#��i��R��#�#��i��R��#�#�c�)��$��F�H�-�	.�A�	
�A�	��"�%�a����$��i��4����6��;�;�	<�	<�	D�	D�Q��	K�	K�	M�B��	�
�v�y�~�X�G�:���	�
��"�'�?�?��%*�$�*.�	0�	0�	0�	0�	��I�B�B����	�	�	�	�)0��%��	�
�v�y�~�~�~����%��	�
�v�y�~�9�9�9�9s
�AF�%F6g�������?c���ttjgd���z}tjddtjzd��}tj||��\}}|tjd|z��z}tjd|zd|zf���}d�|D��}tj|��}ttj|ddd�	�����dS)
N)g��������g{�G�z��rg{�G�z�?g�������?rrrjr r#c�T�g|]%}|ttj|����f��&Sr��r)r�loggamma�r_�z0s  r4r`z3test_loggamma_taylor_transition.<locals>.<listcomp>y�/��>�>�>�b��G�F�O�B�/�/�0�0�1�>�>�>r6�vIh�%,=r?)
�LOGGAMMA_TAYLOR_RADIUSr$r+r&rr�r'r��flattenr	r-rr/�r1�theta�dzr3r0s     r4�test_loggamma_taylor_transitionrns���
	���*E�*E�*E�!F�!F�F�A��K��1�R�U�7�B�'�'�E��{�1�e�$�$�H�A�u�	
�2�6�"�U�(���	�B�
��a�"�f�a�"�f�n��%�%�'�'�A�>�>�A�>�>�>�G��h�w���G��R�[�'�1�a�e�4�4�4�:�:�<�<�<�<�<r6c��tjdtjt��d��}tjddtjzd��}tj||��\}}|tjd|z��z}tjd|zd|zf�	��}d�|D��}tj
|��}ttj
|ddd	�
�����dS)Ni�rrrrjr r#c�T�g|]%}|ttj|����f��&Sr�rrs  r4r`z(test_loggamma_taylor.<locals>.<listcomp>�rr6r	r?)r$r%�log10r
r&rr�r'r�rr+r	r-rr/rs     r4�test_loggamma_taylorrs���	��C���"8�9�9�2�>�>�A��K��1�R�U�7�B�'�'�E��{�1�e�$�$�H�A�u�	
�2�6�"�U�(���	�B�
��a�"�f�a�"�f�n��%�%�'�'�A�>�>�A�>�>�>�G��h�w���G��R�[�'�1�a�e�4�4�4�:�:�<�<�<�<�<r6c��tjtjddd��dtjddd��f}|���}tj||��\}}|d|zz}tjddd���ddd��}|tj|f|jz��z�	��}tjd��5d	�|D��}ddd��n#1swxYwYtj|��}ttj|ddd
������dS)Nr}��rqrr iV���r#r<c�T�g|]%}|ttj|����f��&Sr�)r)r�rgammars  r4r`z%test_rgamma_zeros.<locals>.<listcomp>�s/��@�@�@��B���
�b� 1� 1�2�2�3�@�@�@r6rar?)r$r�r%�copyr�rvr��dstack�sizerrrAr+r	r-rr/)�dx�dyrr�r3r0s      r4�test_rgamma_zerosr�sz��
����B��Q�'�'�'��B�K��R��,C�,C�C�	D�B�	�����B�
�[��R�
 �
 �F�B��	�b��e��B��I�a��r�"�"�*�*�1�a��4�4�E�	���B�5���+�,�,�	,�5�5�7�7�A�	���	�	�A�A�@�@�a�@�@�@��A�A�A�A�A�A�A�A�A�A�A����A�A�A�A��h�w���G��R�Y���A�E�2�2�2�8�8�:�:�:�:�:s�$
C=�=D�Dc	�R�tjtjd��}t|��g}tjtjd��}|�t|����tj|��}t
jdtjddd��dtjddd��df}|�	��}tj
||��\}}|d	|zz}|tj|f|jz��z�
��}tjd
��5d�|D��}ddd��n#1swxYwYtj|��}tt j|ddd
������dS)NrnrJg���Q�οr}��rrg���Q��?r rc�T�g|]%}|ttj|����f��&Sr��r)r�digammars  r4r`z&test_digamma_roots.<locals>.<listcomp>��/��A�A�A��B����r� 2� 2�3�3�4�A�A�Ar6r#r�r?)r�findrootr"rBr(r$r+r�r%rr�rrrrAr	r-r/)�root�rootsrrrr3r0s       r4�test_digamma_rootsr'�s����?�6�>�3�/�/�D�
�4�[�[�M�E��?�6�>�4�0�0�D�	�L�L��t�������H�U�O�O�E�
��u�r�{�2�s�B�/�/�/��B�K��R��4L�4L�d�R�	S�B�	�����B�
�[��R�
 �
 �F�B��	�b��e��B�	���B�5���+�,�,�	,�5�5�7�7�A�	���	�	�B�B�A�A�q�A�A�A��B�B�B�B�B�B�B�B�B�B�B����B�B�B�B��h�w���G��R�Z��!�Q�U�3�3�3�9�9�;�;�;�;�;s�?
E�E�Ec�D��ttj���tjddd��}tjtjddd��dtjddd��f}tj||��\}}|d|zz���}tjd��5�fd	�|D��}ddd��n#1swxYwYtj	|��}ttj|dd
d������dS)
N�,���r<rr�rir �(c�D��g|]}|t�|����f��Sr�)r))r_rr"s  �r4r`z(test_digamma_negreal.<locals>.<listcomp>�s.���:�:�:�"�B�������,�,�-�:�:�:r6r#r>r?)
rrr"r$r%r�r�rrArDr	r-r/)rGrVr3r0r"s    @r4�test_digamma_negrealr-�sD����v�~�.�.�G�	��S�#�s�	#�	#�#�A�
��r�{�1�b�!�$�$�$�a���R��A�)>�)>�>�?�A�
�;�q�!���D�A�q�	
�R��T������A�	���	�	�;�;�:�:�:�:��:�:�:��;�;�;�;�;�;�;�;�;�;�;����;�;�;�;��j��!�!�G��R�Z��!�Q�U�3�3�3�9�9�;�;�;�;�;s�6C�C�Cc��tjddd��}tjgd���}tj||��\}}|d|zz���}tjd��5d�|D��}ddd��n#1swxYwYtj|��}ttj
|dd	d
������dS)Nr)r*r<)gffffff�g�������g������@gffffff@r rc�T�g|]%}|ttj|����f��&Sr�r!rs  r4r`z)test_digamma_boundary.<locals>.<listcomp>�r#r6rr#r>r?)r$r%r+r�rrrArDr	r-r"r/�rGrVr3r0s    r4�test_digamma_boundaryr1�s#��

��S�#�s�	#�	#�#�A�
��'�'�'�(�(�A�
�;�q�!���D�A�q�	
�R��T������A�	���	�	�B�B�A�A�q�A�A�A��B�B�B�B�B�B�B�B�B�B�B����B�B�B�B��j��!�!�G��R�Z��!�Q�U�3�3�3�9�9�;�;�;�;�;s�4
B
�
B�Bc��d}tjd|zd|zd��}|���}tj||��\}}|���|���}}tjd��5d�t||��D��}ddd��n#1swxYwYtj|��}ttj|ddd�	�����dS)
NrjrKrr8r<c�b�g|],\}}||ttj||d�����f��-S)T�rx�regularized)rBr�gammainc)r_�a0�x0s   r4r`z*test_gammainc_boundary.<locals>.<listcomp>�sL��,�,�,��r�2���E�&�/�"���"M�"M�"M�N�N�O�,�,�,r6r=rar?)
r$r&rr�rrrA�zipr+r	r-r6r/)�smallryrGr0s    r4�test_gammainc_boundaryr;�s(��
�E�
��C��I�q��w��+�+�A�	�����A�
�;�q�!���D�A�q��9�9�;�;��	�	���q�A�	���	�	�,�,�,�,�!$�Q����,�,�,��,�,�,�,�,�,�,�,�,�,�,����,�,�,�,��h�w���G��R�[�'�6�1�5�9�9�9�?�?�A�A�A�A�As�B.�.B2�5B2c���d��tjdd��}tjddtz��}dtj|tjd|z����z���}tj�fd�|D����}ttj	|ddd	�
���
��dS)Nc�L�ttjdd|z
����S�Nrr#)r)r�polylog�r3s r4�spencez"test_spence_circle.<locals>.spences ���v�~�a��Q��/�/�0�0�0r6rKrnrrr#r c�*��g|]}|�|��f��Sr�r�)r_rrAs  �r4r`z&test_spence_circle.<locals>.<listcomp>
s&���7�7�7�r�2�v�v�b�z�z�*�7�7�7r6r�r?)r$r&r�outerr'rrDr	r-rAr/)r1r
r3r0rAs    @r4�test_spence_circlerDs����1�1�1�	��C����A��K��1�R�4� � �E�	
�R�X�a����5��)�)�
*�
*�	*�3�3�5�5�A��j�7�7�7�7�Q�7�7�7�8�8�G��R�Y���A�E�2�2�2�8�8�:�:�:�:�:r6c�x�tjt��j}tjtjddd��dtjddd��f}|���}tj||��\}}|d|zz}tjddd���	ddd��}|tj
|f|jz��z���}tj
d	�|D����}tt|ddd
|z������dS)Nrrrqr ���r<r#r}c�T�g|]%}|ttj|����f��&Sr�)r)r�sinpirs  r4r`z$test_sinpi_zeros.<locals>.<listcomp>�?��'�'�'� ��w�v�|�B�'7�'7�8�8�9�'�'�'r6rr?)r$�finforB�epsr�r%rr�rvr�rrrrDr	rr/�rKrrrr�r3r0s       r4�test_sinpi_zerosrMs��
�(�5�/�/�
�C�	����A�s�A�&�&�&��2�;�s�A�q�+A�+A�A�	B�B�	�����B�
�[��R�
 �
 �F�B��	�b��e��B��I�d�C��#�#�+�+�A�q�"�5�5�E�	���B�5���+�,�,�	,�5�5�7�7�A��j�'�'�$%�'�'�'�(�(�G��V�W�a���3��/�/�/�5�5�7�7�7�7�7r6c�~�tjt��j}tjtjddd��dtjddd��f}|���}tj||��\}}|d|zz}tjddd��dz�	ddd	��}|tj
|f|jz��z���}tj
d
�|D����}tt|ddd|z������dS)
Nrrrqr rFr<r#rKr}c�T�g|]%}|ttj|����f��&Sr�)r)r�cospirs  r4r`z$test_cospi_zeros.<locals>.<listcomp>-rIr6rr?)r$rJrBrKr�r%rr�rvr�rrrrDr	rr/rLs       r4�test_cospi_zerosrQ$s$��
�(�5�/�/�
�C�	����A�s�A�&�&�&��2�;�s�A�q�+A�+A�A�	B�B�	�����B�
�[��R�
 �
 �F�B��	�b��e��B�
�Y�t�S�!�
$�
$�s�
*�3�3�A�q�"�=�=�E�	���B�5���+�,�,�	,�5�5�7�7�A��j�'�'�$%�'�'�'�(�(�G�
�V�W�a���3��/�/�/�5�5�7�7�7�7�7r6c
���d�}d�}tjddd��}tjtjddd��dtjddd��f}g}|D]O}t	tj|����}|D])}||z}|�|||||��f���*�Ptj|��}t||d	d
d����
��dS)
Nc�8�tj||��dS�Nr�r-�ellipj��u�ms  r4�dnz"test_dn_quarter_period.<locals>.dn9s���y��A���q�!�!r6c�J�ttjd||�����S�NrZrW)rBr�ellipfunrWs  r4�	mpmath_dnz)test_dn_quarter_period.<locals>.mpmath_dn<s!���V�_�T�Q�!�4�4�4�5�5�5r6rr#rjr}rrr=rr�r?)r$r&r�r%rBr�ellipkr(rDr	r/)	rZr^rY�dur0�m0�u0�du0r2s	         r4�test_dn_quarter_periodrd7s��"�"�"�6�6�6�	��A�q�"���A�	����B��R�(�(�(�!�R�[��b�"�-E�-E�E�	F�B��G��6�6��
�6�=��$�$�
%�
%���	6�	6�C��S��A��N�N�A�r�9�9�Q��#3�#3�4�5�5�5�5�	6��j��!�!�G��R��&�!�%�0�0�0�6�6�8�8�8�8�8r6c�@�tj|��5tj|��}tj|jtjz
dtjzz��}tjtj|��|��}ddd��n#1swxYwY|SrT)rrA�mpc�ceil�imagr�lambertwr')r3r��unwindrWs    r4�_mpmath_wrightomegarkPs���	���	�	�5�5��J�q�M�M����a�f�v�y�0�1�V�Y�;�?�@�@���o�f�j��m�m�V�4�4��5�5�5�5�5�5�5�5�5�5�5����5�5�5�5��Js�A2B�B�Bc��tjddd��}tjtjtj��g}tjtjtj��g}tjtjtj��g}tjtjtj��g}td��D]�}|�tj|dtj����|�tj|dtj����|�tj|dtj����|�tj|dtj������tj||||f��}tj||��\}}|d|zz�	��}tj
d�|D����}ttj
|ddd	�
�����dS)Nrr�r8r}r c	�L�g|]!}|tt|d����f��"S�rm�r)rkrs  r4r`z+test_wrightomega_branch.<locals>.<listcomp>i�@��'�'�'� ��w�':�2�r�'B�'B�C�C�D�'�'�'r6r#�:�0�yE>r?)r$r%�	nextafterr�inf�ranger(�hstackr�rrDr	r-�wrightomegar/)	rG�picut_above�picut_below�npicut_above�npicut_below�irVr3r0s	         r4�test_wrightomega_branchr}Xs���
��R��B�	�	��A��<���r�v�.�.�/�K��<�����w�/�/�0�K��L�"�%����0�0�1�L��L�"�%��"�&��1�1�2�L�
�2�Y�Y�E�E�����2�<��B����@�@�A�A�A����2�<��B��"�&��A�A�B�B�B����B�L��b�)9�2�6�B�B�C�C�C����B�L��b�)9�B�F�7�C�C�D�D�D�D�
�	�;��\�<�H�I�I�A�
�;�q�!���D�A�q�	
�R��T������A��j�'�'�$%�'�'�'�(�(�G�
�R�^�W�a���6�6�6�<�<�>�>�>�>�>r6c�t�tjdd��}tjddtjz��}tj||��\}}|d|zz���}tjd�|D����}t
tj|ddd����	��dS)	Nrr#rr c	�L�g|]!}|tt|d����f��"Srorprs  r4r`z,test_wrightomega_region1.<locals>.<listcomp>xrqr6rrrr?�
r$r&rr�rrDr	r-rwr/r0s    r4�test_wrightomega_region1r�os���	��B����A�
��A�q���w���A�
�;�q�!���D�A�q�	
�R��T������A��j�'�'�$%�'�'�'�(�(�G�
�R�^�W�a���7�7�7�=�=�?�?�?�?�?r6c�t�tjdd��}tjdtjzd��}tj||��\}}|d|zz���}tjd�|D����}t
tj|ddd����	��dS)	Nrr#r}r c	�L�g|]!}|tt|d����f��"Srorprs  r4r`z,test_wrightomega_region2.<locals>.<listcomp>�rqr6rrrr?r�r0s    r4�test_wrightomega_region2r�~s���	��B����A�
��B�r�u�H�b�!�!�A�
�;�q�!���D�A�q�	
�R��T������A��j�'�'�$%�'�'�'�(�(�G�
�R�^�W�a���7�7�7�=�=�?�?�?�?�?r6c�^�tjddd��tjddd��}}tj||��\}}|d|zz���}tjd�|D����}ttj|ddd������dS)	Nr}r#rmr c�T�g|]%}|ttj|����f��&Sr�)r)rrirs  r4r`z(test_lambertw_smallz.<locals>.<listcomp>�s?��'�'�'� ��w�v��r�':�':�;�;�<�'�'�'r6rr>r?)	r$r&r�rrDr	r-rir/r0s    r4�test_lambertw_smallzr��s����;�r�1�b�!�!�2�;�r�1�b�#9�#9�q�A�
�;�q�!���D�A�q�	
�R��T������A��j�'�'�$%�'�'�'�(�(�G�
�R�[�'�1�a�e�4�4�4�:�:�<�<�<�<�<r6r;��maxprec�maxtermsz0.17c���eZdZd�Zd�Zd�Zd�Zd�Zd�Zd�Z	d�Z
d	�Zd
�Zd�Z
d�Zd
�Zd�Zd�Zd�Zd�Zd�Zd�Zd�Zd�Zd�Zd�Zd�Zd�Zej�dd���d���Zd�Z ej�dd���d���Z!d �Z"d!�Z#d"�Z$d#�Z%d$�Z&d%�Z'd&�Z(d'�Z)d(�Z*d)�Z+d*�Z,d+�Z-d,�Z.d-�Z/d.�Z0d/�Z1d0�Z2d1�Z3d2�Z4d3�Z5d4�Z6d5�Z7d6�Z8d7�Z9d8�Z:d9�Z;d:�Z<d;�Z=d<�Z>d=�Z?d>�Z@d?�ZAd@�ZBdA�ZCdB�ZDdC�ZEdD�ZFdE�ZGdF�ZHdG�ZIdH�ZJdI�ZKdJ�ZLdK�ZMdL�ZNdM�ZOej�d�N��dO���ZPdP�ZQej�d�N��dQ���ZReSdR���ZTdS�ZUdT�ZVej�ddU���dV���ZWdW�ZXdX�ZYdY�ZZej�d�N��dZ���Z[eSd[���Z\ej�d�N��d\���Z]ej�^d]��d^���Z_d_�Z`d`�Zaej�d�N��da���Zbdb�Zcdc�Zddd�ZeeSde���Zfdf�Zgej�^dg��dh���Zhdi�ZieSdj���Zjdk�Zkdl�Zldm�Zmdn�Znej�ddo���dp���ZoeSdq���Zpdr�Zqds�Zrej�d�N��dt���Zsej�ddu���dv���Ztdw�Zuej�ddx���dy���Zvdz�Zwd{�Zxej�d|�}��d~���Zyd�Zzd��Z{d��Z|d��Z}d��Z~d��Zd��Z�d��Z�d��Z�d��Z�d��Z�d��Z�d��Z�d��Z�d��Z�d��Z�d��Z�d��Z�d��Z�d��Z�d��Z�d��Z�d��Z�d��Z�d��Z�ej�dd����d����Z�d�S)��TestSystematicc��td�tjtdd��gd���td�tjtdd��g��dS)	Nc�6�tj|��dS�Nr�r-�airyr@s r4rTz,TestSystematic.test_airyai.<locals>.<lambda>����b�g�a�j�j��m�r6��ח���חA��h㈵��>r?c�6�tj|��dSr�r�r@s r4rTz,TestSystematic.test_airyai.<locals>.<lambda>�r�r6�@���@�@)rr�airyair��selfs r4�test_airyaizTestSystematic.test_airyai�sl���3�3�"�M� ��s�^�^�,�!%�	'�	'�	'�	'�	�3�3�"�M� ��s�^�^�,�	.�	.�	.�	.�	.r6c�X�td�tjt��g��dS)Nc�6�tj|��dSr�r�r@s r4rTz4TestSystematic.test_airyai_complex.<locals>.<lambda>�r�r6)rrr�r
r�s r4�test_airyai_complexz"TestSystematic.test_airyai_complex�s2���3�3�"�M�'�\�\�N�	,�	,�	,�	,�	,r6c��td�d�tdd��gd���td�d�td	d
��g��dS)Nc�6�tj|��dS�Nr#r�r@s r4rTz2TestSystematic.test_airyai_prime.<locals>.<lambda>�r�r6c�.�tj|d���S�Nr#)�
derivative�rr�r@s r4rTz2TestSystematic.test_airyai_prime.<locals>.<lambda>����"�M�!��:�:�:�r6r�r�r�r?c�6�tj|��dSr�r�r@s r4rTz2TestSystematic.test_airyai_prime.<locals>.<lambda>�r�r6c�.�tj|d���Sr�r�r@s r4rTz2TestSystematic.test_airyai_prime.<locals>.<lambda>�r�r6r�r��rrr�s r4�test_airyai_primez TestSystematic.test_airyai_prime��v���3�3�6;�6;� ��s�^�^�,�!%�	'�	'�	'�	'�	�3�3�6;�6;� ��s�^�^�,�	.�	.�	.�	.�	.r6c�F�td�d�t��g��dS)Nc�6�tj|��dSr�r�r@s r4rTz:TestSystematic.test_airyai_prime_complex.<locals>.<lambda>�r�r6c�.�tj|d���Sr�r�r@s r4rTz:TestSystematic.test_airyai_prime_complex.<locals>.<lambda>�r�r6�rr
r�s r4�test_airyai_prime_complexz(TestSystematic.test_airyai_prime_complex��7���3�3�6;�6;�'�\�\�N�	,�	,�	,�	,�	,r6c��td�d�tdd��gd���td�d�td	d
��g��dS)Nc�6�tj|��dSrTr�r@s r4rTz,TestSystematic.test_airybi.<locals>.<lambda>�r�r6c�*�tj|��SrN�r�airybir@s r4rTz,TestSystematic.test_airybi.<locals>.<lambda>����"�M�!�,�,�r6r�r�r�r?c�6�tj|��dSrTr�r@s r4rTz,TestSystematic.test_airybi.<locals>.<lambda>�r�r6c�*�tj|��SrNr�r@s r4rTz,TestSystematic.test_airybi.<locals>.<lambda>�r�r6r�r�r�r�s r4�test_airybizTestSystematic.test_airybi�sv���3�3�6-�6-� ��s�^�^�,�!%�	'�	'�	'�	'�	�3�3�6-�6-� ��s�^�^�,�	.�	.�	.�	.�	.r6c�F�td�d�t��g��dS)Nc�6�tj|��dSrTr�r@s r4rTz4TestSystematic.test_airybi_complex.<locals>.<lambda>�r�r6c�*�tj|��SrNr�r@s r4rTz4TestSystematic.test_airybi_complex.<locals>.<lambda>�r�r6r�r�s r4�test_airybi_complexz"TestSystematic.test_airybi_complex�s7���3�3�6-�6-�'�\�\�N�	,�	,�	,�	,�	,r6c��td�d�tdd��gd���td�d�td	d
��g��dS)Nc�6�tj|��dS�Nrqr�r@s r4rTz2TestSystematic.test_airybi_prime.<locals>.<lambda>�r�r6c�.�tj|d���Sr�r�r@s r4rTz2TestSystematic.test_airybi_prime.<locals>.<lambda>�r�r6r�r�r�r?c�6�tj|��dSr�r�r@s r4rTz2TestSystematic.test_airybi_prime.<locals>.<lambda>�r�r6c�.�tj|d���Sr�r�r@s r4rTz2TestSystematic.test_airybi_prime.<locals>.<lambda>�r�r6r�r�r�r�s r4�test_airybi_primez TestSystematic.test_airybi_prime�r�r6c�F�td�d�t��g��dS)Nc�6�tj|��dSr�r�r@s r4rTz:TestSystematic.test_airybi_prime_complex.<locals>.<lambda>�r�r6c�.�tj|d���Sr�r�r@s r4rTz:TestSystematic.test_airybi_prime_complex.<locals>.<lambda>�r�r6r�r�s r4�test_airybi_prime_complexz(TestSystematic.test_airybi_prime_complex�r�r6c�v�ttjtd���t	dd��g��dS)Nc�2�tjd|fit��Sr�)r�bei�HYPERKWr@s r4rTz)TestSystematic.test_bei.<locals>.<lambda>����v�z�!�Q�7R�7R�'�7R�7R�r6r�r�)rr-r�rrr�s r4�test_beizTestSystematic.test_bei��?���B�F�,�-R�-R�S�S� ��s�^�^�,�	.�	.�	.�	.�	.r6c�v�ttjtd���t	dd��g��dS)Nc�2�tjd|fit��Sr�)r�berr�r@s r4rTz)TestSystematic.test_ber.<locals>.<lambda>�r�r6r�r�)rr-r�rrr�s r4�test_berzTestSystematic.test_ber�r�r6c�P�td�d�tdd��gdd���dS)Nc�j�tjt|����t|��SrN)r-�	bernoullir��rFs r4rTz/TestSystematic.test_bernoulli.<locals>.<lambda>�s!��b�l�3�q�6�6�&:�&:�3�q�6�6�&B�r6c�^�ttjt|������SrN)rBrr�r�r�s r4rTz/TestSystematic.test_bernoulli.<locals>.<lambda>�s��e�F�,<�S��V�V�,D�,D�&E�&E�r6ri�2r��r@rF)rrr�s r4�test_bernoullizTestSystematic.test_bernoulli�sD���B�B�E�E�#�A�u�-�-�.�!%��	0�	0�	0�	0�	0�	0r6c��ttjtd���t	dd��t	��gd���dS)Nc�2�tj||fit��SrN�r�besselir�rRs  r4rTz-TestSystematic.test_besseli.<locals>.<lambda>����&�.��A�:Y�:Y�QX�:Y�:Y�r6�}Ô%�I���}Ô%�I�TgJD������ru)rr-�ivrrr�s r4�test_besselizTestSystematic.test_besseli�sO���B�E�,�-Y�-Y�Z�Z� ���/�/����7�!'�	)�	)�	)�	)�	)�	)r6c�~�td�td���tdd��t��g��dS)Nc�6�tj|j|��SrN)r-r�rQrRs  r4rTz5TestSystematic.test_besseli_complex.<locals>.<lambda>������q�v�q�)9�)9�r6c�2�tj||fit��SrNr�rRs  r4rTz5TestSystematic.test_besseli_complex.<locals>.<lambda>�r�r6r�r��rrrr
r�s r4�test_besseli_complexz#TestSystematic.test_besseli_complex��N���9�9�,�-Y�-Y�Z�Z� ���/�/����>�	@�	@�	@�	@�	@r6c	�,�ttjtd���t	dd��t	dd��gd���ttjtd���t	dd��t	d	d
��gdd���dS)
Nc�2�tj||fit��SrN�r�besseljr�rRs  r4rTz-TestSystematic.test_besselj.<locals>.<lambda>r�r6r�r�r�r�T�r�c�2�tj||fit��SrNr�rRs  r4rTz-TestSystematic.test_besselj.<locals>.<lambda>r�r6r�r�r�)r�r@)rr-�jvrrr�s r4�test_besseljzTestSystematic.test_besseljs����B�E�,�-Y�-Y�Z�Z� ���/�/��T�3���@�,0�	2�	2�	2�	2�	�B�E�,�-Y�-Y�Z�Z� ���/�/��T�3���@�,0�!%�		'�	'�	'�	'�	'�	'r6c�z�td�td���t��t��g��dS)Nc�6�tj|j|��SrN)r-r�rQrRs  r4rTz5TestSystematic.test_besselj_complex.<locals>.<lambda>r�r6c�2�tj||fit��SrNr�rRs  r4rTz5TestSystematic.test_besselj_complex.<locals>.<lambda>r�r6r�r�s r4�test_besselj_complexz#TestSystematic.test_besselj_complex
sC���9�9�,�-Y�-Y�Z�Z� �U�U�J�L�L�1�	3�	3�	3�	3�	3r6c	��ttjtjtdd��tdtj��gdd���dS�Nr�r;rFra��nan_okr@)rr-�kvr�besselkrr$rtr�s r4�test_besselkzTestSystematic.test_besselksK���B�E�"�N� ��s�^�^�S��B�F�^�^�<�#(�u�	6�	6�	6�	6�	6�	6r6c	��ttjtjtdd��t
dtj��gdd���dSr�)	rr-�knrr�rrr$rtr�s r4�test_besselk_intzTestSystematic.test_besselk_intsM���B�E�"�N�#�D�#�.�.��A�r�v���?�#(�u�	6�	6�	6�	6�	6�	6r6c�~�td�td���tdd��t��g��dS)Nc�6�tj|j|��SrN)r-r�rQrRs  r4rTz5TestSystematic.test_besselk_complex.<locals>.<lambda>r�r6c�2�tj||fit��SrN)rr�r�rRs  r4rTz5TestSystematic.test_besselk_complex.<locals>.<lambda> r�r6r�r�r�r�s r4�test_besselk_complexz#TestSystematic.test_besselk_complexr�r6c	��d�}ttjt|��t	dd��t	dd��gd���dS)Nc��ttj||fit����}t	|��dkr!t
jtj|��z}t	|��dkr|dkrt
jS|S)N�قnQ:Br)	rBr�besselyr�r�r$rt�sign�nan�rSrGr1s   r4�	mpbesselyz.TestSystematic.test_bessely.<locals>.mpbessely$si���f�n�Q��5�5�W�5�5�6�6�A��1�v�v��~�~��F�R�W�Q�Z�Z�'���1�v�v��{�{�q�A�v�v��v�
��Hr6r�r�r�r��r�)rr-�yvrr�r�rs  r4�test_besselyzTestSystematic.test_bessely#s`��	�	�	�	�B�E�,�Y�7�7� ���/�/��T�3���@�"�	$�	$�	$�	$�	$�	$r6c��d�}td�t|��t��t��gd���dS)Nc��ttj||fit����}t	|��dkrNtjd���5t
jtj|��z}ddd��n#1swxYwY|S)Nr
r�r�)	r)rrr�r�r$r�rtrrs   r4rz6TestSystematic.test_bessely_complex.<locals>.mpbessely3s������q�!�7�7�w�7�7�8�8�A��1�v�v��~�~��[��2�2�2�,�,��������+�A�,�,�,�,�,�,�,�,�,�,�,����,�,�,�,��Hs�"A<�<B�Bc�6�tj|j|��SrN)r-rrQrRs  r4rTz5TestSystematic.test_bessely_complex.<locals>.<lambda>:r�r6i�:r�r�rs  r4�test_bessely_complexz#TestSystematic.test_bessely_complex2sX��	�	�	�	�9�9�,�Y�7�7� �U�U�J�L�L�1�#�	%�	%�	%�	%�	%�	%r6c	��d�}td�t|��tdd��tdd��g��dS)Nc��ttj||����}t|��dkr|dkrtjS|Sr�)rBrrr�r$r
rs   r4rz2TestSystematic.test_bessely_int.<locals>.mpbessely@s>���f�n�Q��*�*�+�+�A��1�v�v��{�{�q�A�v�v��v�
��Hr6c�F�tjt|��|��SrN)r-�ynr�rRs  r4rTz1TestSystematic.test_bessely_int.<locals>.<lambda>Fs����s�1�v�v�q�)9�)9�r6r���r�r��rrrrrs  r4�test_bessely_intzTestSystematic.test_bessely_int?s_��	�	�	�	�9�9�,�Y�7�7�#�E�4�0�0�#�d�C�.�.�A�	C�	C�	C�	C�	Cr6c�����g�d	�fd�	�ttj�fd�t��t��gdd���ttj�t	j���ddd���dS)
NFc�@��|dks|dkrtjS|dks|dkrett||z����dzdkr?|r=��t|��t|��f��tjStj||��S)Ng��m�rr#)r$r
r�rBr(rr�)ryrx�nonzero�
bad_pointss   �r4r�z&TestSystematic.test_beta.<locals>.betaMs�����5�y�y�A��I�I��v�
��A����Q���S��q�1�u���%6�%6��%:�q�$@�$@��"��%�%�u�Q�x�x��q���&:�;�;�;��6�M��;�q�!�$�$�$r6c� ���||d���S�NT)r!r�)ryrxr�s  �r4rTz*TestSystematic.test_beta.<locals>.<lambda>]s�����a��D�)A�)A�)A�r6r�T�r�r��dy���=)r�r�ru�F)rr-r�rr$r+)r�r"r�s @@r4rzTestSystematic.test_betaJs������
�
	%�
	%�
	%�
	%�
	%�
	%�	�B�G�A�A�A�A� �U�U�C�E�E�N� #�,0�		2�	2�	2�	2�	�B�G� ��H�Z�0�0� #�,0�!&�	(�	(�	(�	(�	(�	(r6c���ttjt��t	d�����t��t��t��g��dS)Nc�4�tj||d|d���S)NrT)r5)r�betainc�ryrxrGs   r4rTz-TestSystematic.test_betainc.<locals>.<lambda>ks ��F�N�[\�^_�ab�de�sw�Lx�Lx�Lx�r6)rr-r*rrrr�s r4�test_betainczTestSystematic.test_betaincisV���B�J�*�L�N�N�+;�<x�<x�+y�+y�z�z� �U�U�C�E�E�3�5�5�1�	3�	3�	3�	3�	3r6c�����g�d�fd�	�ttj�fd�t��t��gd���ttj�t	j���dd���dS)	NFc	���t|��dt|��dzzkrtjS||kr�tt||z
��tjt||z
����z
��dkr?|r=��t|��t|��f��tjSt
j||��S)Nr�r#rr)r�r$r
rBr�r(r�binomial)rF�kr!r"s   �r4r/z+TestSystematic.test_binom.<locals>.binomialqs�����1�v�v��S��V�V�a�Z�(�(�(��v�
��1�u�u��U�1�Q�3�Z�Z�"�(�5��1��:�:�*>�*>�>�?�?�%�G�G��"��%�%�u�Q�x�x��q���&:�;�;�;��6�M��?�1�a�(�(�(r6c� ���||d���Sr$r�)rFr0r/s  �r4rTz+TestSystematic.test_binom.<locals>.<lambda>�s�����!�Q��)E�)E�)E�r6r��r�r��r�rur')rr-�binomrr$r+)r�r"r/s @@r4�
test_binomzTestSystematic.test_binomns������
�
	)�
	)�
	)�
	)�
	)�
	)�	�B�H�E�E�E�E� �U�U�C�E�E�N� #�	%�	%�	%�	%�
	�B�H�$��H�Z�0�0� #�!&�		(�	(�	(�	(�	(�	(r6c�~�td�td���t��t��gd���dS)Nc�F�tjt|��|��SrN)r-�eval_chebytr��rFrGs  r4rTz0TestSystematic.test_chebyt_int.<locals>.<lambda>�������A����)B�)B�r6c�2�tj||fit��SrN)r�chebytr�r9s  r4rTz0TestSystematic.test_chebyt_int.<locals>.<lambda>����&�-��1�:X�:X�PW�:X�:X�r6r8r2rr�s r4�test_chebyt_intzTestSystematic.test_chebyt_int��J���B�B�,�-X�-X�Y�Y�#�X�X�s�u�u�-�2�	7�	7�	7�	7�	7�	7r6Fz'some cases in hyp2f1 not fully accurate)�run�reasonc�z�ttjd�tdd��t��gd���dS)Nc�x�t��ttj����||fit��SrN)rrrr<r�r9s  r4rTz,TestSystematic.test_chebyt.<locals>.<lambda>��7��)X�����8H���8W�8W�)X�)X�YZ�\]�)i�)i�ah�)i�)i�r6����e�'r�)rr-r8rr�s r4�test_chebytzTestSystematic.test_chebyt�sC���B�N�i�i� ��s�^�^�S�U�U�3�u�	>�	>�	>�	>�	>�	>r6c�~�td�td���t��t��gd���dS)Nc�F�tjt|��|��SrN)r-�eval_chebyur�r9s  r4rTz0TestSystematic.test_chebyu_int.<locals>.<lambda>�r:r6c�2�tj||fit��SrN)r�chebyur�r9s  r4rTz0TestSystematic.test_chebyu_int.<locals>.<lambda>�r=r6r8r2rr�s r4�test_chebyu_intzTestSystematic.test_chebyu_int�r?r6c�v�ttjd�tdd��t��g��dS)Nc�x�t��ttj����||fit��SrN)rrrrMr�r9s  r4rTz,TestSystematic.test_chebyu.<locals>.<lambda>�rDr6rErF)rr-rKrr�s r4�test_chebyuzTestSystematic.test_chebyu�s<���B�N�i�i� ��s�^�^�S�U�U�3�	5�	5�	5�	5�	5r6c��d�}t|tjt��g��t|tjt	gd���g��dS)Nc�6�tj|��dSr��r-�shichir�s r4�chiz$TestSystematic.test_chi.<locals>.chi�����9�Q�<�<��?�"r6�g���U@�Xg�V@)rrrVrr�r�rVs  r4�test_chizTestSystematic.test_chi��W��	#�	#�	#��C���c�e�e�W�5�5�5��C���h�7Q�7Q�7Q�.R�.R�-S�T�T�T�T�Tr6c���d�}t|tjtt	t
jd��t	t
jd����gd���dS)Nc�6�tj|��dSr�rTr@s r4rVz,TestSystematic.test_chi_complex.<locals>.chi�rWr6r�r�rar?)rrrVr
r)r$rtrZs  r4�test_chi_complexzTestSystematic.test_chi_complex��i��	#�	#�	#�	�C�"�J�'������(>�(>����PS�@T�@T�U�U�V�!&�	(�	(�	(�	(�	(�	(r6c�`�d�}t|tjtdd��g��dS)Nc�6�tj|��dSr��r-�sicir�s r4�ciz"TestSystematic.test_ci.<locals>.ci�����7�1�:�:�a�=� r6r�r�)rrrer�r�res  r4�test_cizTestSystematic.test_ci�sD��	!�	!�	!�	�B�"�I� ��s�^�^�,�	.�	.�	.�	.�	.r6c���d�}t|tjtt	dt
j��t	dt
j����gd���dS)Nc�6�tj|��dSr�rcr@s r4rez*TestSystematic.test_ci_complex.<locals>.ci�rfr6r�r�rrr?)rrrer
r)r$rtrgs  r4�test_ci_complexzTestSystematic.test_ci_complex�sh��	!�	!�	!�	�B�"�I�'���r�v�g�(>�(>���R�V�@T�@T�U�U�V�!%�	'�	'�	'�	'�	'�	'r6c��tjt��j}t	t
tjt��gdd|z���dS�NFrr�)	r$rJrBrKrrrrPr�r�rKs  r4�
test_cospizTestSystematic.test_cospi�sK���h�u�o�o�!���F�"�L� �U�U�G�E��#��	?�	?�	?�	?�	?�	?r6c�f�tttjt	��gdd���dS)NFr>r�)rrrrPr
r�s r4�test_cospi_complexz!TestSystematic.test_cospi_complex�s=���F�"�L�'�\�\�N�5�u�	F�	F�	F�	F�	F�	Fr6c��ttjttj��t��gdd���dS)Nrar8�r@r�)rr-r"rrrr�s r4�test_digammazTestSystematic.test_digamma�sA���B�J�,�V�^�<�<� �U�U�G�%�R�	9�	9�	9�	9�	9�	9r6c��d�}ttjttj��t��gdd|���dS)Nc�|�tj|jdktj|j��dkzdd��S)Nrg�Q����?FT)r$�whererQr�rhr@s r4�param_filterz9TestSystematic.test_digamma_complex.<locals>.param_filter�s0���8�Q�V�a�Z�B�F�1�6�N�N�T�,A�B�E�4�P�P�Pr6r>r+�r@r�rx)rr-r"rrr
)r�rxs  r4�test_digamma_complexz#TestSystematic.test_digamma_complex�sX��	Q�	Q�	Q�	�B�J�,�V�^�<�<�'�\�\�N��B�)5�	7�	7�	7�	7�	7�	7r6c�n�ttjtjt��gd���dS)Nr�r?)rr-�exp1r�e1rr�s r4�test_e1zTestSystematic.test_e1�s7���B�G�"�I� �U�U�G�%�	1�	1�	1�	1�	1�	1r6c
��ttjtjtt
tjd��t
tjd����gd���ttjtjtj	ddd��dd�dftj
dtjd	d
d��tjd	d
d��fd
zz���d���ttjtjtj	ddd��dzd���dS)Nr�r�r&r?����r8�rr�r�=r�r i���rGy)
rr-r|rr}r
r)r$rtr&r�r%r�r�s r4�test_e1_complexzTestSystematic.test_e1_complex�s#���B�G�"�I�'������(>�(>����PS�@T�@T�U�U�V�!&�	(�	(�	(�	(�	�B�G�"�I��[��b�#�6�6�q�q�q�$�w�?��U�1�b�k�"�a��&<�&<�$&�K��A�r�$:�$:�#:�$;�<�<>�?�?�@E����!&�	(�	(�	(�	(�	�B�G�"�I��[��c�5�9�9�B�>�!&�	(�	(�	(�	(�	(�	(r6c��ttjd�tt	jt	jtj��j��t	jt	jtj��j�����g��ttjd�t	j	dddddtj
g��d	�
��tt	jtjtj
������ttjtj
��dk��dS)Nc�d�|dkrtj|��|zntjd��S�Nrz1.0�r�expm1�mpfr�s r4rTz,TestSystematic.test_exprel.<locals>.<lambda>��+��1��6�6�f�l�1�o�o�a�&7�&7�v�z�RW�GX�GX�r6r�c�d�|dkrtj|��|zntjd��Sr�r�r�s r4rTz,TestSystematic.test_exprel.<locals>.<lambda>�r�r6rag�����W�:rg��mBg���yCx�Dr&r?)
rr-�exprelrr$r�rJ�doubler�r+rtr�isinfr�s r4�test_exprelzTestSystematic.test_exprel�s����B�I�X�X� �B�F�2�8�B�I�+>�+>�+B�$C�$C�#C�r�v�b�h�WY�W`�Na�Na�Ne�Gf�Gf�g�g�g�h�	j�	j�	j�	�B�I�X�X��H�e�U�A�t�T�2�6�%J�K�K�RW�	Y�	Y�	Y�	Y�	�����2�6�*�*�+�+�,�,�,���	�2�6�'�"�"�a�'�(�(�(�(�(r6c���ttjtjt	ttjd��ttjd����g��dS)N��c���cA)rr-r�rr
r)r$rtr�s r4�test_expm1_complexz!TestSystematic.test_expm1_complex�sW���B�H�"�L�'������(>�(>����PS�@T�@T�U�U�V�	X�	X�	X�	X�	Xr6c�\�ttjd�t��gd���dS)Nc�0�tj|dz��Sr��rr�r�s r4rTz3TestSystematic.test_log1p_complex.<locals>.<lambda>s��f�j��1��o�o�r6�<r2)rr-�log1pr
r�s r4�test_log1p_complexz!TestSystematic.test_log1p_complexs8���B�H�5�5�'�\�\�N��	4�	4�	4�	4�	4�	4r6c�T�ttd�t��gdd���dS)Nc�6�tj|dz��|z
Sr�r�r�s r4rTz-TestSystematic.test_log1pmx.<locals>.<lambda>s��f�j��Q��&7�&7�!�&;�r6r�r�)r�r@)rrrr�s r4�test_log1pmxzTestSystematic.test_log1pmx
s8���H�;�;� �U�U�G��%�	9�	9�	9�	9�	9�	9r6c�n�ttjtjt��gd���dS)Nr&r?)rr-r.rr*rr�s r4�test_eizTestSystematic.test_eis8���B�G�"�I� �U�U�G�!&�	(�	(�	(�	(�	(�	(r6c���ttjtjtt
tjd��t
tjd����gd���dS)Nr�r�r�r?)	rr-r.rr*r
r)r$rtr�s r4�test_ei_complexzTestSystematic.test_ei_complexsZ���B�G�"�I�'������(>�(>����PS�@T�@T�U�U�V�!%�	'�	'�	'�	'�	'�	'r6c�n�ttjtjt	d���g��dS)Nr��rx)rr-�elliperrr�s r4�test_ellipezTestSystematic.test_ellipes4���B�I�"�M� �3�Z�Z�Z�L�	*�	*�	*�	*�	*r6c��ttjtjtdd��td���g��dS)Nr�r�r�r��rr-�	ellipeincrr�rr�s r4�test_ellipeinczTestSystematic.test_ellipeinc!s?���B�L�"�M� ��s�^�^�S�3�Z�Z�Z�8�	:�	:�	:�	:�	:r6c��ttjtjt��t��g��dSrNr�r�s r4�test_ellipeinc_largephiz&TestSystematic.test_ellipeinc_largephi&�6���B�L�"�M� �U�U�C�E�E�N�	,�	,�	,�	,�	,r6c��ttjtjtdd��t��g��dS�Nr�r��rr-�	ellipkincr�ellipfrr�s r4�test_ellipfzTestSystematic.test_ellipf+s;���B�L�"�M� ��s�^�^�S�U�U�3�	5�	5�	5�	5�	5r6c��ttjtjt��t��g��dSrNr�r�s r4�test_ellipf_largephiz#TestSystematic.test_ellipf_largephi0r�r6c���ttjtjt	d���g��ttjd�t	d���gd���dS)Nr�r�c�0�tjd|z
��Sr�)rr_)rYs r4rTz,TestSystematic.test_ellipk.<locals>.<lambda>:s��f�m�A��E�&:�&:�r6g)ryr�r2)rr-r_rr�ellipkm1r�s r4�test_ellipkzTestSystematic.test_ellipk5sg���B�I�"�M� �3�Z�Z�Z�L�	*�	*�	*�	�B�K�:�:� �3�Z�Z�Z�L� #�	%�	%�	%�	%�	%�	%r6c��d�}ttj|tdd��td���gd���dS)Nc�.�tjd||��Sr��r�ellippi��phirYs  r4r�z0TestSystematic.test_ellipkinc.<locals>.ellipkinc?����>�!�S�!�,�,�,r6r�r�r�r�Tr��rr-r�r�r�r�s  r4�test_ellipkinczTestSystematic.test_ellipkinc>sT��	-�	-�	-��B�L�%� ��s�^�^�S�3�Z�Z�Z�8�,0�	2�	2�	2�	2�	2�	2r6c�~�d�}ttj|t��td���gd���dS)Nc�.�tjd||��Sr�r�r�s  r4r�z9TestSystematic.test_ellipkinc_largephi.<locals>.ellipkincGr�r6r�r�Tr�r�r�s  r4�test_ellipkinc_largephiz&TestSystematic.test_ellipkinc_largephiFsP��	-�	-�	-��B�L�%� �U�U�C�#�J�J�J�/�,0�	2�	2�	2�	2�	2�	2r6c	�r�d�}td�|tdd��tdd���gd�	��dS)
Nc�@�|dkrdStjd||���S)Nr�snrW�rr]rWs  r4r�z+TestSystematic.test_ellipfun_sn.<locals>.snOs'���A�v�v��q���t�q�A�6�6�6�6r6c�8�tj||��dSr�rUrWs  r4rTz1TestSystematic.test_ellipfun_sn.<locals>.<lambda>Y�����1�a����);�r6���.����.Arr#r�rrr?r�)r�r�s  r4�test_ellipfun_snzTestSystematic.test_ellipfun_snNsZ��	7�	7�	7�	�;�;�� ��s�^�^�S�1��]�]�]�;�!%�	'�	'�	'�	'�	'�	'r6c	�n�td�d�tdd��tdd���gd�	��dS)
Nc�8�tj||��dSr�rUrWs  r4rTz1TestSystematic.test_ellipfun_cn.<locals>.<lambda>`r�r6c�0�tjd||���S)N�cnrWr�rWs  r4rTz1TestSystematic.test_ellipfun_cn.<locals>.<lambda>a�������a�)H�)H�)H�r6r�r�rr#r�rrr?r�r�s r4�test_ellipfun_cnzTestSystematic.test_ellipfun_cn^�L���;�;�H�H� ��s�^�^�S�1��]�]�]�;�!%�	'�	'�	'�	'�	'�	'r6c	�n�td�d�tdd��tdd���gd�	��dS)
Nc�8�tj||��dSrTrUrWs  r4rTz1TestSystematic.test_ellipfun_dn.<locals>.<lambda>gr�r6c�0�tjd||���Sr\r�rWs  r4rTz1TestSystematic.test_ellipfun_dn.<locals>.<lambda>hr�r6r�r�rr#r�rrr?r�r�s r4�test_ellipfun_dnzTestSystematic.test_ellipfun_dner�r6c�X�ttjd�t��g��dS)Nc�*�tj|��SrN�rr�r@s r4rTz)TestSystematic.test_erf.<locals>.<lambda>n���f�j��m�m�r6)rr-r�rr�s r4�test_erfzTestSystematic.test_erfls1���B�F�3�3� �U�U�G�	%�	%�	%�	%�	%r6c�\�ttjd�t��gd���dS)Nc�*�tj|��SrNr�r@s r4rTz1TestSystematic.test_erf_complex.<locals>.<lambda>sr�r6r;r�)rr-r�r
r�s r4r�zTestSystematic.test_erf_complexqs8���B�F�3�3�'�\�\�N�c�	3�	3�	3�	3�	3�	3r6c�v�ttjtd���t	��gd���dS)Nc�*�tj|��SrN�rr�r@s r4rTz*TestSystematic.test_erfc.<locals>.<lambda>x���v�{�1�~�~�r6r>r?)rr-r�rrr�s r4�	test_erfczTestSystematic.test_erfcvsA���B�G�,�-E�-E�F�F� �U�U�G�%�	1�	1�	1�	1�	1�	1r6c�v�ttjtd���t	��gd���dS)Nc�*�tj|��SrNr�r@s r4rTz2TestSystematic.test_erfc_complex.<locals>.<lambda>}r�r6r;r�)rr-r�rr
r�s r4�test_erfc_complexz TestSystematic.test_erfc_complex{sA���B�G�,�-E�-E�F�F�'�\�\�N�c�	3�	3�	3�	3�	3�	3r6c�n�ttjtjt	��gd���dS�Nr;r�)rr-�erfirrr�s r4�	test_erfizTestSystematic.test_erfi�s7���B�G�"�K� �U�U�G�s�	,�	,�	,�	,�	,�	,r6c�n�ttjtjt	��gd���dSr�)rr-r�rr
r�s r4�test_erfi_complexz TestSystematic.test_erfi_complex�s7���B�G�"�K�'�\�\�N�c�	3�	3�	3�	3�	3�	3r6c�v�ttjtd���t	��gd���dS)Nc�*�tj|��SrN)r�ncdfr@s r4rTz*TestSystematic.test_ndtr.<locals>.<lambda>�r�r6r;r�)rr-�ndtrrrr�s r4�	test_ndtrzTestSystematic.test_ndtr�sA���B�G�,�-E�-E�F�F� �U�U�G�s�	,�	,�	,�	,�	,�	,r6c��ttjd�tt	dd��t	dd�����gd���dS)Nc�\�tj|tjd��z��dzS�Ng@)rr�r$�sqrtr@s r4rTz2TestSystematic.test_ndtr_complex.<locals>.<lambda>�s#��f�k�1�"�R�W�R�[�[�.�&A�&A�"�&D�r6����rGr�r�r�)rr-r�r
r)r�s r4�test_ndtr_complexz TestSystematic.test_ndtr_complex�s^���B�G�D�D�'�'�&�&�*A�*A�W�U�TY�EZ�EZ�[�[�[�\�`c�	e�	e�	e�	e�	e�	er6c�z�ttjtd���t	��gddd���dS)Nc�N�tjtj|����SrN)rr�r�r@s r4rTz.TestSystematic.test_log_ndtr.<locals>.<lambda>�s��v�z�&�+�a�.�.�7Q�7Q�r6iXr)r>�rFr�r@)rr-�log_ndtrrrr�s r4�
test_log_ndtrzTestSystematic.test_log_ndtr�sK���B�K�,�-Q�-Q�R�R� �U�U�G�s��%�	A�	A�	A�	A�	A�	Ar6c��ttjtd���t	tdd��tdd�����gdd�	��dS)
Nc��tjtj|tjd��z��dz��Sr�)rr�r�r$r�r@s r4rTz6TestSystematic.test_log_ndtr_complex.<locals>.<lambda>�s3��v�z�&�+�q�b�QS�QX�Y[�Q\�Q\�n�B]�B]�^`�B`�7a�7a�r6r�rFrGr<r�r;r)�rFr�)rr-r�rr
r)r�s r4�test_log_ndtr_complexz$TestSystematic.test_log_ndtr_complex�sr���B�K�,�-a�-a�b�b�'�'�&�$�*?�*?�*1�%��*=�*=�?�?�?�@�CF�C�	Q�	Q�	Q�	Q�	Q�	Qr6c�`�td�tjtdd��gd���dS)Nc�6�tj|��dS)Nr})r-�eulerr�s r4rTz.TestSystematic.test_eulernum.<locals>.<lambda>�s��b�h�q�k�k�"�o�r6r#rGr�)rr�eulernumrr�s r4�
test_eulernumzTestSystematic.test_eulernum�s@���5�5�"�O�#�A�u�-�-�.�%�	9�	9�	9�	9�	9�	9r6c	��ttjtjtdd��t
dtj��gdd���dS)Nrr;r>�rs)	rr-rErrCrrr$rtr�s r4�test_expintzTestSystematic.test_expint�sK���B�G�"�M�#�A�s�^�^�S��B�F�^�^�<�!&�C�	1�	1�	1�	1�	1�	1r6c�\�d�}t|tjt��g��dS)Nc�6�tj|��dSr��r-�fresnelr�s r4�fresnelsz.TestSystematic.test_fresnels.<locals>.fresnels�����:�a�=�=��#�#r6)rrrr)r�rs  r4�
test_fresnelszTestSystematic.test_fresnels��=��	$�	$�	$��H�"�O� �U�U�G�	%�	%�	%�	%�	%r6c�\�d�}t|tjt��g��dS)Nc�6�tj|��dSr�rr�s r4�fresnelcz.TestSystematic.test_fresnelc.<locals>.fresnelc�rr6)rrrr)r�rs  r4�
test_fresnelczTestSystematic.test_fresnelc�rr6c��ttjttj��t��g��dSrN)rr-�gammarrrr�s r4�
test_gammazTestSystematic.test_gamma�s8���B�H�,�V�\�:�:� �U�U�G�	%�	%�	%�	%�	%r6c��ttjttj��t��gd���dS�Nrsr?)rr-rrrr
r�s r4�test_gamma_complexz!TestSystematic.test_gamma_complex�s?���B�H�,�V�\�:�:�'�\�\�N��	8�	8�	8�	8�	8�	8r6c	��ttjd�tddd���tdd��gdd���dS)Nc�0�tj||d���S)NTr4�rr6)r3rxs  r4rTz.TestSystematic.test_gammainc.<locals>.<lambda>�������a�T�)R�)R�)R�r6r���@F��inclusive_ar&r�)rr-r6rr�s r4�
test_gammainczTestSystematic.test_gammainc�sQ���B�K�R�R� ��C�U�;�;�;�S��C�[�[�I�#(�u�	6�	6�	6�	6�	6�	6r6c	��ttjd�tddd���tdd��gdd���dS)Nc�0�tj||d���S)NT)ryr5r)r3rys  r4rTz/TestSystematic.test_gammaincc.<locals>.<lambda>�rr6rr Fr!r&r�)rr-�	gammainccrr�s r4�test_gammaincczTestSystematic.test_gammaincc�sQ���B�L�R�R� ��C�U�;�;�;�S��C�[�[�I�#(�u�	6�	6�	6�	6�	6�	6r6c�v�d�}ttjt|��t	��g��dS)Nc�4�tj|��jSrN�rrrQr@s r4�fz&TestSystematic.test_gammaln.<locals>.f�s���?�1�%�%�*�*r6)rr-�gammalnrr)r�r+s  r4�test_gammalnzTestSystematic.test_gammaln�s<��	+�	+�	+�	�B�J�(8��(;�(;�c�e�e�W�E�E�E�E�Er6)r@c��ttjttj��t
dd��t
��t
��g��dSr�)rr-�eval_gegenbauerrr�
gegenbauerrr�s r4�test_gegenbauerzTestSystematic.test_gegenbauer�sK���B�.�,�V�->�?�?� ��s�^�^�S�U�U�C�E�E�:�	<�	<�	<�	<�	<r6c
�l�d�}d�}t|t|��tdd��tdd��t��gdddd	�
��t|t|��tdd��t��t	tjddd
����gdd���dS)Nc�2�t|��dkrtjS|dkrd}n%|dkr	d|z|z}ntj|||��}t|��dkr�|dkr�t|��t
t|����krdtj||tjd��z|��}t|��tjd��krtjd��}t|��d	krtjS|S)
Nr�rr�r#rr}z1e-50z0.0�\��4@I�w)	r�r$r
rr0rBr�r�rt�rFryrGr1s    r4r0z6TestSystematic.test_gegenbauer_int.<locals>.gegenbauer�s����1�v�v��~�~��v�
��A�v�v�����a����a�C��E����%�a��A�.�.���Q�x�x�1�}�}��R���E�!�H�H��E�!�H�H�
�
�,E�,E��%�a��V�Z��-@�-@�)@�!�D�D���q�6�6�F�J�w�/�/�/�/��
�5�)�)�A��1�v�v��~�~��v�
��Hr6c��tjt|��||��}t|��dkrtjS|S)Nr4)r-r/r�r�r$rtr5s    r4�
sc_gegenbauerz9TestSystematic.test_gegenbauer_int.<locals>.sc_gegenbauer�s7���"�3�q�6�6�1�a�0�0�A��1�v�v��~�~��v�
��Hr6rr<ge��ge��Ai@�Tr�)rFr�r�r@r*���rr%)rrrrrr$r%)r�r0r7s   r4�test_gegenbauer_intz"TestSystematic.test_gegenbauer_int�s���	�	�	�6	�	�	�	�M�,�Z�8�8�#�A�s�^�^�S��s�^�^�S�U�U�C�#��,0�t�		=�	=�	=�	=�	�M�,�Z�8�8�#�A�s�^�^�S�U�U�H�R�[��b�RT�=U�=U�4V�4V�W� #�,0�		2�	2�	2�	2�	2�	2r6c��td�ttj��t	dd��t��t
��g��dS)Nc�R�tjt|��|j|��SrN)r-r/r�rQ�rFryrGs   r4rTz8TestSystematic.test_gegenbauer_complex.<locals>.<lambda>s��B�,>�s�1�v�v�q�v�q�,Q�,Q�r6rr<)rrrr0rrr
r�s r4�test_gegenbauer_complexz&TestSystematic.test_gegenbauer_complex
sQ���Q�Q�,�V�->�?�?�#�A�s�^�^�S�U�U�J�L�L�A�	C�	C�	C�	C�	Cr6c��td�ttj��t	dd��t	��t��g��dS)Nc�B�tj|j|j|��SrN)r-r/rQr<s   r4rTz@TestSystematic.test_gegenbauer_complex_general.<locals>.<lambda>s��B�,>�q�v�q�v�q�,Q�,Q�r6r�r�)rrrr0rr
r�s r4�test_gegenbauer_complex_generalz.TestSystematic.test_gegenbauer_complex_generalsQ���Q�Q�,�V�->�?�?� ��s�^�^�S�U�U�J�L�L�A�	C�	C�	C�	C�	Cr6c��ttjtd���t	dd��t	��g��dS)Nc�2�tj||fit��SrN)r�hankel1r��rSrGs  r4rTz-TestSystematic.test_hankel1.<locals>.<lambda>s'��&�.��A�;U�;U�LS�;U�;U�r6�@��x���@��x�D)rr-rCrrr�s r4�test_hankel1zTestSystematic.test_hankel1sS���B�J�,�.U�.U�V�V� ���-�-�s�u�u�5�	7�	7�	7�	7�	7r6c��ttjtd���t	dd��t	��g��dS)Nc�2�tj||fit��SrN)r�hankel2r�rDs  r4rTz-TestSystematic.test_hankel2.<locals>.<lambda>!r�r6rErF)rr-rJrrr�s r4�test_hankel2zTestSystematic.test_hankel2sG���B�J�,�-Y�-Y�Z�Z� ���-�-�s�u�u�5�	7�	7�	7�	7�	7r6z%issues at intermediately large ordersc��td�ttj��t	dd��t��g��dS)Nc�F�tjt|��|��SrN)r-�eval_hermiter�r9s  r4rTz-TestSystematic.test_hermite.<locals>.<lambda>&������Q����)C�)C�r6rrG)rrr�hermiterrr�s r4�test_hermitezTestSystematic.test_hermite$sG���C�C�,�V�^�<�<�#�A�u�-�-�s�u�u�5�	7�	7�	7�	7�	7r6c	���tdd����ttj�fd�t	dd��t	dd��gd	�
��dS)Nr�i�r�c�*��tj||fi���SrNr^)ryrG�KWs  �r4rTz,TestSystematic.test_hyp0f1.<locals>.<lambda>1s�����q�!�)B�)B�r�)B�)B�r6r�r�rgj�@rr�)r�rr-rPr)r�rTs @r4�test_hyp0f1zTestSystematic.test_hyp0f1,sa���
�#��
-�
-�
-���B�I�B�B�B�B� ��s�^�^�S��C�[�[�9�"�	$�	$�	$�	$�	$�	$r6c��td�td���tdd��tt	dd��t	dd����g��dS)Nc�6�tj|j|��SrNrO)ryr3s  r4rTz4TestSystematic.test_hyp0f1_complex.<locals>.<lambda>:s����1�6�1�)=�)=�r6c�2�tj||fit��SrN)rrPr��ryrGs  r4rTz4TestSystematic.test_hyp0f1_complex.<locals>.<lambda>;r=r6r�ri����rL)rrrr
r)r�s r4�test_hyp0f1_complexz"TestSystematic.test_hyp0f1_complex9sh���=�=�,�-X�-X�Y�Y� ��b�\�\�:�g�d�D�6I�6I�7�SV�X[�K\�K\�+]�+]�^�	`�	`�	`�	`�	`r6c
��d�}ttj|tdd��tddd���tdd��gdd���dS)	Nc�h�	tj|||��S#t$rtjcYSwxYwrN)r�hyp1f1�ZeroDivisionErrorr$rtr+s   r4�
mpmath_hyp1f1z1TestSystematic.test_hyp1f1.<locals>.mpmath_hyp1f1BsA��
��}�Q��1�-�-�-��$�
�
�
��v�
�
�
�
���s��1�1r�r8r#Fr!r�)rFr�)rr-r]r)r�r_s  r4�test_hyp1f1zTestSystematic.test_hyp1f1Ask��	�	�	�	��I��
��b�\�\�3�q�"�%�8�8�8�#�c�2�,�,�G���	
�	
�	
�	
�	
�	
r6c	��ttd���td���tdd��tdd��t	��gd���dS)Nc�B�tj|j|j|��SrN)r-r]rQr+s   r4rTz4TestSystematic.test_hyp1f1_complex.<locals>.<lambda>Rs��r�y�����QR�7S�7S�r6c�4�tj|||fit��SrN)rr]r�r+s   r4rTz4TestSystematic.test_hyp1f1_complex.<locals>.<lambda>S���V�]�1�a�QR�=^�=^�V]�=^�=^�r6r�r�i�r�)rrrrr
r�s r4�test_hyp1f1_complexz"TestSystematic.test_hyp1f1_complexPsa���J�'S�'S�T�T�,�-^�-^�_�_� ��s�^�^�S��s�^�^�Z�\�\�J�"�	$�	$�	$�	$�	$�	$r6c
��td�td���tdd��tdd��tdd��t��gd���dS)Nc�N�tj|j|j|j|��SrN)r-r�rQ�ryrxr�rGs    r4rTz4TestSystematic.test_hyp2f1_complex.<locals>.<lambda>Zs��r�y�������QR�/S�/S�r6c�6�tj||||fit��SrN)rr�r�rhs    r4rTz4TestSystematic.test_hyp2f1_complex.<locals>.<lambda>[�"���
�a�QR�TU�WX�@d�@d�\c�@d�@d�r6gY�gY@rr�r�r�s r4�test_hyp2f1_complexz"TestSystematic.test_hyp2f1_complexWse��	�S�S�,�-d�-d�e�e� ��s�^�^�S��s�^�^�S��s�^�^�Z�\�\�Z� �	"�	"�	"�	"�	"�	"r6c��ttjtd���t	��t	��t	��g��dS)Nc�4�tj|||fit��SrN)rrwr�r+s   r4rTz,TestSystematic.test_hyperu.<locals>.<lambda>brdr6)rr-rwrrr�s r4�test_hyperuzTestSystematic.test_hyperu_sG���B�I�,�-^�-^�_�_� �U�U�C�E�E�3�5�5�1�	3�	3�	3�	3�	3r6z:mpmath issue gh-342: unsupported operand mpz, long for powc	�|�d�}tt|tddd���tdd��gd���dS)Nc��tj||��tj|��ztj|��zSrN)r�powerr'rrYs  r4�mp_igam_facz1TestSystematic.test_igam_fac.<locals>.mp_igam_facgs0���<��1�%�%�f�j�!��n�n�4�V�\�!�_�_�D�Dr6rg�ļ�BFr!r�r?)rrr)r�rrs  r4�
test_igam_faczTestSystematic.test_igam_faces^��	E�	E�	E�	�I�'� ��D�e�<�<�<�c�!�T�l�l�K�!&�	(�	(�	(�	(�	(�	(r6c���ttjtjt	dd��g��ttjtjt	dd��gd���dS�Nr�r�r�r�r�r?)rr-�j0rrr�s r4�test_j0zTestSystematic.test_j0osj��	�B�E�"�I� ��s�^�^�,�	.�	.�	.�	�B�E�"�I� ��s�^�^�,�!%�	'�	'�	'�	'�	'�	'r6c���ttjtjt	dd��g��ttjtjt	dd��gd���dSru)rr-�j1rrr�s r4�test_j1zTestSystematic.test_j1}sh���B�E�"�I� ��s�^�^�,�	.�	.�	.�	�B�E�"�I� ��s�^�^�,�!%�	'�	'�	'�	'�	'�	'r6c	�h�ttjtd���t	��t	��t	��t	��g��td�td���t��t	��t	��t	��g��dS)Nc�6�tj||||fit��SrN�r�jacobir�rhs    r4rTz,TestSystematic.test_jacobi.<locals>.<lambda>�rjr6c�J�tjt|��|||��SrN�r-�eval_jacobir�)rFrxr�rGs    r4rTz,TestSystematic.test_jacobi.<locals>.<lambda>����r�~�c�!�f�f�a��A�/N�/N�r6c�6�tj||||fit��SrNr}rhs    r4rTz,TestSystematic.test_jacobi.<locals>.<lambda>�rjr6)rr-r�rrrr�s r4�test_jacobizTestSystematic.test_jacobi�s����B�N�,�-d�-d�e�e� �U�U�C�E�E�3�5�5�#�%�%�8�	:�	:�	:�	�N�N�,�-d�-d�e�e�#�X�X�s�u�u�c�e�e�S�U�U�;�	=�	=�	=�	=�	=r6c	���d��td��fd�t��t��t��t��gdd���dS)Nc�@�|dkrdStj||||��S)Nrr�)rr~�rFryrxrGs    r4r~z.TestSystematic.test_jacobi_int.<locals>.jacobi�s&���A�v�v��s��=��A�q�!�,�,�,r6c�J�tjt|��|||��SrNr�r�s    r4rTz0TestSystematic.test_jacobi_int.<locals>.<lambda>�r�r6c�@��t���||||fit��SrN)rr�)rFryrxrGr~s    �r4rTz0TestSystematic.test_jacobi_int.<locals>.<lambda>�s+���/G�/?��/G�/G��1�a�QR�/^�/^�V]�/^�/^�r6� Nr8r�rrr)r�r~s @r4�test_jacobi_intzTestSystematic.test_jacobi_int�sh���	-�	-�	-�
	�N�N�^�^�^�^�#�X�X�s�u�u�c�e�e�S�U�U�;�#��	-�	-�	-�	-�	-�	-r6c�d�d�}ttj|tdd��gd���dS)Nc�p�|dkrtdzSttj��d|fit��S)Nrr:)rrr�keir�r�s r4r�z$TestSystematic.test_kei.<locals>.kei�s;���A�v�v��s�1�u��/�#�F�J�/�/��1�@�@��@�@�@r6�ꌠ9Y>)��ꌠ9Y>)Frr�)rr-r�r)r�r�s  r4�test_keizTestSystematic.test_kei�sP��	A�	A�	A�
	�B�F�� ���-�-�.�$�	8�	8�	8�	8�	8�	8r6c�z�ttjtd���t	dd��gd���dS)Nc�2�tjd|fit��Sr�)r�kerr�r�s r4rTz)TestSystematic.test_ker.<locals>.<lambda>�r�r6r�r�rr�)rr-r�rrr�s r4�test_kerzTestSystematic.test_ker�sH���B�F�,�-R�-R�S�S� ���-�-�.�$�	8�	8�	8�	8�	8�	8r6c��tttj��d�t	��t	��g��dS)Nc�N�ttj��||fit��SrN�rr�laguerrer�r9s  r4rTz.TestSystematic.test_laguerre.<locals>.<lambda>��'��)J�)9�&�/�)J�)J�1�a�)[�)[�SZ�)[�)[�r6)rrr-�
eval_laguerrerr�s r4�
test_laguerrezTestSystematic.test_laguerre�s@���J�r�'7�8�8�[�[� �U�U�C�E�E�N�	,�	,�	,�	,�	,r6c�d�td�d�t��t��gd���dS)Nc�F�tjt|��|��SrN)r-r�r�r9s  r4rTz2TestSystematic.test_laguerre_int.<locals>.<lambda>�����)9�#�a�&�&�!�)D�)D�r6c�N�ttj��||fit��SrNr�r9s  r4rTz2TestSystematic.test_laguerre_int.<locals>.<lambda>�r�r6r�r�r�r�s r4�test_laguerre_intz TestSystematic.test_laguerre_int�sA���D�D�[�[�#�X�X�s�u�u�-��	8�	8�	8�	8�	8�	8r6zsee gh-3551 for bad pointsc	��td�d�ttjtj��t	dd��gdd���dS)Nc�P�tj|t|j����SrN)r-rir�rQ�rGr0s  r4rTz3TestSystematic.test_lambertw_real.<locals>.<lambda>�s����Q��A�F���)D�)D�r6c�P�tj|t|j����SrN)rrir�rQr�s  r4rTz3TestSystematic.test_lambertw_real.<locals>.<lambda>�s�����C���K�K�)H�)H�r6rrr>F�r@r�)rr
r$rtrr�s r4�test_lambertw_realz!TestSystematic.test_lambertw_real�sT���D�D�H�H�'������8�8�&��B�-�-�H�!&�u�	6�	6�	6�	6�	6�	6r6c	����d}tjd���d���fd�}t|tjtd|d���gd�	��dS)
NgG�a�se@r#g��E@c� ��tjd���5|�zdz
�z|dz
z}|tjkr|t|��z}n+|�zdz
�zd|dz
zz}|t|��z}||z}ddd��n#1swxYwY|S)Nr�)�overrK)r$r�rtr)rG�facrW�e�gs   ��r4rz:TestSystematic.test_lanczos_sum_expg_scaled.<locals>.gamma�s������(�+�+�+�
�
��A����Q��!�c�'�2���"�&�=�=��6�q�9�9�9�C�C���E�C�K��?�c�1�s�7�m�<�C��6�q�9�9�9�C��3�J�C�
�
�
�
�
�
�
�
�
�
�
����
�
�
�
��Js�A B�B�
BrFr!r>r?)r$r'rrrr)r��maxgammarr�r�s   @@r4�test_lanczos_sum_expg_scaledz+TestSystematic.test_lanczos_sum_expg_scaled�s{����&���F�1�I�I��)��		�		�		�		�		�		�	�E�"�L� ��H�%�@�@�@�A�!&�	(�	(�	(�	(�	(�	(r6c��ttjtjt��t��g��dSrN)rr-�
eval_legendrer�legendrerr�s r4�
test_legendrezTestSystematic.test_legendre�s7���B�,�"�O� �U�U�C�E�E�N�	,�	,�	,�	,�	,r6c���td�d�t��t��gd���td�d�t��tt	jddd	����g��dS)
Nc�F�tjt|��|��SrN�r-r�r�r9s  r4rTz2TestSystematic.test_legendre_int.<locals>.<lambda>�r�r6c�N�ttj��||fit��SrN�rrr�r�r9s  r4rTz2TestSystematic.test_legendre_int.<locals>.<lambda>�r�r6r�r�c�F�tjt|��|��SrNr�r9s  r4rTz2TestSystematic.test_legendre_int.<locals>.<lambda>�r�r6c�N�ttj��||fit��SrNr�r9s  r4rTz2TestSystematic.test_legendre_int.<locals>.<lambda>�r�r6r*r8rj)rrrrr$r%r�s r4�test_legendre_intz TestSystematic.test_legendre_int�s����D�D�[�[�#�X�X�s�u�u�-�#�	%�	%�	%�	%�	�D�D�[�[�#�X�X�x���C��R�0H�0H�'I�'I�J�	L�	L�	L�	L�	Lr6c
��d�}d�}d�}t||tdd��tdd��t��g��t||tdd��tdd��tdd��gd�	��dS)
Nc��	tj|||��dd}n#t$rtjcYSwxYwt|��dkr&tjtj|j��z}|S)Nr�r}r}�)�#���v)	r-�lpmn�
ValueErrorr$r
r�rtrrQ�rFrYr3rSs    r4�lpnmz(TestSystematic.test_legenp.<locals>.lpnm�sv��
��G�A�q�!�$�$�Q�'��.�����
�
�
��v�
�
�
�
�����1�v�v��~�~��F�R�W�Q�V�_�_�,���Hs�"%�>�>c��tj|||��}t|��dkr&tjtj|j��z}|S)Nr�)r-r�r�r$rtrrQr�s    r4�lpnm_2z*TestSystematic.test_legenp.<locals>.lpnm_2�s?�����1�a� � �A��1�v�v��~�~��F�R�W�Q�V�_�_�,���Hr6c���|dks|dkrNt|��|kr;|dkr3|dkr|dz
}tjtj|��|��SdSt	|��dkrt
jSt	|��dkrdnd}ttj��||||���}t	|��dkr&tj	tj|j
��z}|S)	Nr#r}rrrrrq��typer�)r�rrqrr�r$r
rr�rtrQ)rFrYr3�typrSs     r4r�z*TestSystematic.test_legenp.<locals>.legenp�s����Q���!�r�'�'�s�1�v�v��{�{���6�6��1�u�u��B��F��!�<���A����:�:�:��1��1�v�v��~�~��v�
��1�v�v��z�z�!�!�q�C�/� ���/�/��1�a�c�B�B�B�A��1�v�v��~�~��J���Q�V�!4�!4�4���Hr6rFr<r}r#r�r�r�)r�r�r�r�s    r4�test_legenpzTestSystematic.test_legenp�s���	�	�	�	�	�	�	�	�	�0	�D�"�#�D�#�.�.��t�S�0A�0A�3�5�5�I�	K�	K�	K�	�F�"�#�D�#�.�.��D�#����B��
�
�K�!&�	(�	(�	(�	(�	(�	(r6c	�J�d�}d�}tjgd���}tjgd���}|dd�dfd|ddd�fzz���}t||t	gd���t	gd���t	|��gdd�	��dS)
Nc��	tj|j|j|d���ddS#t$rtjcYSwxYw)Nrr�rr��r-�clpmnrQr�r$r
�rFrYr3s   r4�clpnmz3TestSystematic.test_legenp_complex_2.<locals>.clpnm�T��
��x�������:�:�:�1�=�e�D�D���
�
�
��v�
�
�
�
�����-0�A	�A	c���t|��dkrtjSttj��t
|j��t
|j��|d���S�Nrrrr��r�r$r
rrr�r�rQr�s   r4r�z4TestSystematic.test_legenp_complex_2.<locals>.legenp�M���1�v�v��~�~��v�
�2�#�F�M�2�2�3�q�v�;�;��A�F���Q�UV�W�W�W�Wr6�	rr�rJrr�rKr�rjg@�@�r�rJrKg�������?r �rr}rr#rrr�r�r��r$r+r�rr�r�r�r�rGrVr3s      r4�test_legenp_complex_2z$TestSystematic.test_legenp_complex_2����	�	�	�	X�	X�	X�
�H�C�C�C�D�D���H�+�+�+�,�,��
�q�q�q��v�Y��A�d�1�1�1�f�I��
%�,�,�.�.���E�"�%�&;�&;�&;�<�<�h�G\�G\�G\�>]�>]�_g�hi�_j�_j�k�!%�!�		#�	#�	#�	#�	#�	#r6c	�J�d�}d�}tjgd���}tjgd���}|dd�dfd|ddd�fzz���}t||t	gd���t	gd���t	|��gdd�	��dS)
Nc��	tj|j|j|d���ddS#t$rtjcYSwxYw)Nrqr�rr�r�r�s   r4r�z3TestSystematic.test_legenp_complex_3.<locals>.clpnm/r�r�c���t|��dkrtjSttj��t
|j��t
|j��|d���S)Nrrrqr�r�r�s   r4r�z4TestSystematic.test_legenp_complex_3.<locals>.legenp5r�r6r�r�r r�r�r�r�r�r�s      r4�test_legenp_complex_3z$TestSystematic.test_legenp_complex_3.r�r6z*apparently picks wrong function at |z| > 1c	��d�}d�}t||tdd��tdd��t��g��dS)Nc�F�tj|||��ddS�Nrr�)r-�lqmnr�s   r4�lqnmz(TestSystematic.test_legenq.<locals>.lqnmHs ���7�1�a��#�#�A�&�u�-�-r6c��t|��dkrtjSttj��|||d���Sr�)r�r$r
rr�legenqr�s   r4r�z*TestSystematic.test_legenq.<locals>.legenqKs<���1�v�v��~�~��v�
�2�#�F�M�2�2�1�a���C�C�C�Cr6rr<r��r�r�r�s   r4�test_legenqzTestSystematic.test_legenqFsi��	.�	.�	.�	D�	D�	D�	�D�"�#�A�s�^�^�V�A�s�^�^�S�U�U�C�	E�	E�	E�	E�	Er6c	��d�}d�}t||tdd��tdd��t��gd���dS)Nc��tjt|j��t|j��|��ddSr�)r-r�r�rQr�s   r4r�z0TestSystematic.test_legenq_complex.<locals>.lqnmWs0���7�3�q�v�;�;��A�F���Q�7�7��:�5�A�Ar6c���t|��dkrtjSttj��t
|j��t
|j��|d���Sr�)r�r$r
rrr�r�rQr�s   r4r�z2TestSystematic.test_legenq_complex.<locals>.legenqZr�r6rr<r�)rrr
r�s   r4�test_legenq_complexz"TestSystematic.test_legenq_complexUso��	B�	B�	B�	X�	X�	X�	�D�"�#�A�s�^�^�V�A�s�^�^�Z�\�\�J�!�	#�	#�	#�	#�	#�	#r6c�`�d�}d�}tt|t��gdd|���dS)Nc�h�tjtj|��|k|dkzdd��S)NrFT)r$rw�floorr�s r4rxz0TestSystematic.test_lgam1p.<locals>.param_filterfs,���8�R�X�a�[�[�A�-�!�q�&�9�5�$�G�G�Gr6c�:�tjd|z��jSr�r*r@s r4�	mp_lgam1pz-TestSystematic.test_lgam1p.<locals>.mp_lgam1pjs���?�1�q�5�)�)�.�.r6r>r<ry)rrr)r�rxr�s   r4�test_lgam1pzTestSystematic.test_lgam1pes[��	H�	H�	H�	/�	/�	/�	�G�%� �U�U�G�%�S�)5�	7�	7�	7�	7�	7�	7r6c�d�d�}ttj|t��gddd���dS)Nc��	tj|��}n4#t$r'ttjtj��}YnwxYw|SrN)rrr�r)r$r
)r3rWs  r4�mpmath_loggammaz5TestSystematic.test_loggamma.<locals>.mpmath_loggammatsO��
.��o�a�(�(�����
.�
.�
.��b�f�b�f�-�-����
.�����Js��.A�AFr	)r��distinguish_nan_and_infr@)rr-rr
)r�r�s  r4�
test_loggammazTestSystematic.test_loggammassQ��	�	�	�	�B�K�+�'�\�\�N�5�49��	G�	G�	G�	G�	G�	Gr6c�~�d�}t|td���t��t��g��dS)Nc�8�tj||��dSr�)r-�pbdvrDs  r4�pcfdz&TestSystematic.test_pcfd.<locals>.pcfd�����7�1�a�=�=��#�#r6c�2�tj||fit��SrN)rr�r�rDs  r4rTz*TestSystematic.test_pcfd.<locals>.<lambda>�s��&�+�a��:V�:V�g�:V�:V�r6)rrr)r�r�s  r4�	test_pcfdzTestSystematic.test_pcfd�sM��	$�	$�	$��D�,�-V�-V�W�W� �U�U�C�E�E�N�	,�	,�	,�	,�	,r6zHit's not the same as the mpmath function --- maybe different definition?c�h�d�}t|d�t��t��gd���dS)Nc�8�tj||��dSr�)r-�pbvvrDs  r4�pcfvz&TestSystematic.test_pcfv.<locals>.pcfv�r�r6c�x�t��ttj����||fit��SrN)rrrrr�rDs  r4rTz*TestSystematic.test_pcfv.<locals>.<lambda>�s7��)V�����8H���8U�8U�)V�)V�WX�Z[�)g�)g�_f�)g�)g�r6rr�r�)r�rs  r4�	test_pcfvzTestSystematic.test_pcfv�sK��	$�	$�	$��D�g�g� �U�U�C�E�E�N�d�	4�	4�	4�	4�	4�	4r6c	���d�}d�}d�}t|tjtdd��tdd��gdd���t||tdd��tdd��gd	d���dS)
Nc�8�tj||��dSr��r-�pbwarYs  r4�pcfwz&TestSystematic.test_pcfw.<locals>.pcfw�r�r6c�8�tj||��dSr�r
rYs  r4�dpcfwz'TestSystematic.test_pcfw.<locals>.dpcfw�r�r6c�F�tjtj||fd��S)Nr=)r�diffrrYs  r4�mpmath_dpcfwz.TestSystematic.test_pcfw.<locals>.mpmath_dpcfw�s���;�v�{�Q��F�F�;�;�;r6rhrig:�0�yU>r<r�g��&�.!>)rrrr)r�rrrs    r4�	test_pcfwzTestSystematic.test_pcfw�s���	$�	$�	$�	$�	$�	$�	<�	<�	<�
	�D�"�K� ��Q�Z�Z��R����4�4�3�	H�	H�	H�	H�	�E�(� ��Q�Z�Z��R����4�4�3�	H�	H�	H�	H�	H�	Hr6zCissues at large arguments (atol OK, rtol not) and <eps-close to z=0c���ttjt��t	t
j����t
dd��t��g��dS)Nrr)rr-�	polygammarrrrrr�s r4�test_polygammazTestSystematic.test_polygamma�sS���B�L�*�L�N�N�+;�F�<L�+M�+M�N�N�#�A�t�_�_�c�e�e�4�	6�	6�	6�	6�	6r6c��ttjtjt	dt
j��gddd���dS)Ni��rFT)rFr�r�)rr-rrrr$rtr�s r4�test_rgammazTestSystematic.test_rgamma�sH����I��M�
����
�
� ��� �
	
�	
�	
�	
�	
�	
r6c��ttjttj��t��gd���dSr)rr-rrrr
r�s r4�test_rgamma_complexz"TestSystematic.test_rgamma_complex�s?���B�I�,�V�]�;�;�'�\�\�N��	8�	8�	8�	8�	8�	8r6zNsee gh-3551 for bad points on 32 bit systems and gh-8095 for another bad point)rAc��tjtj��tjd��kr
tj}nd�}t
tj|t��t��gd���dS)Nr{c��t||z��t||z��kr?t||z��dkr)tj|��}t||z��|z
}tj||��Sr�)rBr�rr��rf)ryrYs  r4�mppochz&TestSystematic.test_rf.<locals>.mppoch�sg����Q��<�<�3�q�1�u�:�:�-�-�%��A��,�,�!�2C�2C��
�1�
�
�A��A��E�
�
�Q��A��y��A���&r6r�r2)
r�parserr��Versionrrr-�pochr)r�rs  r4�test_rfzTestSystematic.test_rf�sz���=��+�,�,����0H�0H�H�H��Y�F�F�
'�
'�
'�	�B�G�"� �U�U�C�E�E�N� #�	%�	%�	%�	%�	%�	%r6c��tjt��j}t	t
tjt��gdd|z���dSrm)	r$rJrBrKrrrrHrrns  r4�
test_sinpizTestSystematic.test_sinpi�sJ���h�u�o�o�!���F�F�L� �U�U�G�E��#��	?�	?�	?�	?�	?�	?r6c�f�tttjt	��gdd���dS)NFg�+����=r�)rrrrHr
r�s r4�test_sinpi_complexz!TestSystematic.test_sinpi_complex�s<���F�F�L�'�\�\�N�5�u�	F�	F�	F�	F�	F�	Fr6c��d�}t|tjt��g��t|tjt	gd���g��dS)Nc�6�tj|��dSr�rTr�s r4�shiz$TestSystematic.test_shi.<locals>.shi�rWr6rX)rrr(rr�r�r(s  r4�test_shizTestSystematic.test_shi�r\r6c���d�}t|tjtt	t
jd��t	t
jd����gd���dS)Nc�6�tj|��dSr�rTr@s r4r(z,TestSystematic.test_shi_complex.<locals>.shi�rWr6r�r�rar?)rrr(r
r)r$rtr)s  r4�test_shi_complexzTestSystematic.test_shi_complex�r`r6c�\�d�}t|tjt��g��dS)Nc�6�tj|��dSr�rcr�s r4�siz"TestSystematic.test_si.<locals>.si�rfr6)rrr0r�r�r0s  r4�test_sizTestSystematic.test_si�s1��	!�	!�	!��B��	�C�E�E�7�3�3�3�3�3r6c���d�}t|tjtt	dt
j��t	dt
j����gd���dS)Nc�6�tj|��dSr�rcr@s r4r0z*TestSystematic.test_si_complex.<locals>.si�rfr6r�r�rar?)rrr0r
r)r$rtr1s  r4�test_si_complexzTestSystematic.test_si_complex�sh��	!�	!�	!�	�B�"�I�'���r�v�g�(>�(>���R�V�@T�@T�U�U�V�!&�	(�	(�	(�	(�	(�	(r6c��d�}ttjt|��t	dt
j��gd���dS)Nc�2�tjdd|z
��Sr>�rr?r�s r4�dilogz)TestSystematic.test_spence.<locals>.dilog�����>�!�Q��U�+�+�+r6rr�r?)rr-rArrr$rt�r�r9s  r4�test_spencezTestSystematic.test_spence�sU��	,�	,�	,�	�B�I�,�U�3�3� ��B�F�^�^�,�5�	:�	:�	:�	:�	:�	:r6c�z�d�}ttjt|��t	��gd���dS)Nc�2�tjdd|z
��Sr>r8r@s r4r9z1TestSystematic.test_spence_complex.<locals>.dilogr:r6r�r?)rr-rArr
r;s  r4�test_spence_complexz"TestSystematic.test_spence_complex�sL��	,�	,�	,��B�I�,�U�3�3�'�\�\�N��	8�	8�	8�	8�	8�	8r6c���d�}t|tjtdd��tdd��t	dt
���t	ddt
z���gddd�	��dS)
Nc�T�||krtjStj||||��SrN)r$r
r-�sph_harm)�lrYr
r�s    r4�	spherharmz0TestSystematic.test_spherharm.<locals>.spherharms(���1�u�u��v�
��;�q�!�S�%�0�0�0r6rr<r�rrriprZ)rurFr�)rrrDrrr)r�rDs  r4�test_spherharmzTestSystematic.test_spherharms~��	1�	1�	1�	�I�"�,�#�A�s�^�^�V�A�s�^�^� �1��^�^�^�S�1��"��-=�-=�-=�?�!%�� #�	%�	%�	%�	%�	%�	%r6c	��ttjttj��t
dd��t
dd��gd���dS)N���r r��&�.>r?)rr-�struverr�struvehrr�s r4�test_struvehzTestSystematic.test_struvehsO���B�I�,�V�^�<�<� ��s�^�^�S��C�[�[�9�!&�	(�	(�	(�	(�	(�	(r6c	��d�}ttjt|��t	dd��t	dd��gdd���dS)Nc�*�|dkry||krrt|��dkr_tjj}	dtj_tj||��|tj_S#|tj_wxYwtj||��S)Nrrr))r�rr�r��struvel)rSr3r�s   r4�
mp_struvelz/TestSystematic.test_struvel.<locals>.mp_struvels{���1�u�u��a�R���C��F�F�T�M�M� �)�-��,�$'�F�I�M�!�>�!�Q�/�/�$+�F�I�M�M��G�F�I�M�+�+�+�+��>�!�Q�'�'�'s�%A*�*A=rGr rrHTr�)rr-�	modstruverr)r�rOs  r4�test_struvelzTestSystematic.test_struvelsa��		(�		(�		(�	�B�L�,�Z�8�8� ��s�^�^�S��C�[�[�9�!&�,0�		2�	2�	2�	2�	2�	2r6c�h�d�}ttj|tdd��gddd���dS)Nc�t�tjtj|��tjd����S)Nz-0.5)rrir'r�r�s r4�mpmath_wrightomega_realzETestSystematic.test_wrightomega_real.<locals>.mpmath_wrightomega_real+s&���?�6�:�a�=�=�&�*�V�2D�2D�E�E�Er6r�gP����KDg�+�����<rF)r@rur�)rr-rwr)r�rTs  r4�test_wrightomega_realz$TestSystematic.test_wrightomega_real*sX��	F�	F�	F�	��N�#�
���
�
�����
	
�	
�	
�	
�	
�	
r6c�^�ttjd�t��gdd���dS)Nc�"�t|d��S)Nrm)rkr@s r4rTz1TestSystematic.test_wrightomega.<locals>.<lambda><s��&9�!�R�&@�&@�r6r�Fr�)rr-rwr
r�s r4�test_wrightomegazTestSystematic.test_wrightomega:s@���B�N�@�@�'�\�\�N��u�	F�	F�	F�	F�	F�	Fr6c	��ttjttj��tddd���tdd���g��dS)Nr#� _�BF)ryrxr"r�ryr")rr-�zetarrrr�s r4�test_hurwitz_zetaz TestSystematic.test_hurwitz_zeta?sV���B�G�,�V�[�9�9� �1��%�@�@�@� �1�%�8�8�8�:�	;�	;�	;�	;�	;r6c�b�ttjd�tdd��gdd���dS)Nc�N�|dkrtj|��ntjSr��rr\rtr�s r4rTz2TestSystematic.test_riemann_zeta.<locals>.<lambda>Hs���Q���f�k�!�n�n�n�F�J�r6rFr<Frsr�)rr-r\rr�s r4�test_riemann_zetaz TestSystematic.test_riemann_zetaEsB����G�>�>�
��s�^�^����	
�	
�	
�	
�	
�	
r6c�d�ttjd�tdd��gddd���dS)Nc�T�|dkrtj|��dz
ntjSr�r`r�s r4rTz+TestSystematic.test_zetac.<locals>.<lambda>Ps*��*+�q�&�&�(.�{�1�~�~��'9�'9�6<�j�r6rFr<F�-rs)r�r�r@)rr-�zetacrr�s r4�
test_zetaczTestSystematic.test_zetacNsH���B�H�B�B� ��s�^�^�,�#(�b�u�		>�	>�	>�	>�	>�	>r6c��d�}ttjt|��t	dd���t	��gddd���dS)	Nc��tj�|��}tj�|��}|dkrtj�|��Stj�||��|zSr��rr�r�r��powm1)rG�lmbdas  r4�	mp_boxcoxz-TestSystematic.test_boxcox.<locals>.mp_boxcoxWs`���	�
�
�a� � �A��I�M�M�%�(�(�E���z�z��y�}�}�Q�'�'�'��y���q�%�0�0�5�8�8r6rFr[r;r�r>r�)rr-�boxcoxrr)r�rls  r4�test_boxcoxzTestSystematic.test_boxcoxUse��	9�	9�	9�	�B�I�,�Y�7�7� �1�%�8�8�8�#�%�%�@�!� "�!&�	(�	(�	(�	(�	(�	(r6c��d�}ttjt|��t	dd���t	��gddd���dS)	Nc�X�tj�|��}tj�|��}tj�d��}|dkr"tj�||z��Stj�||z|��|zSr�ri)rGrk�ones   r4�mp_boxcox1pz1TestSystematic.test_boxcox1p.<locals>.mp_boxcox1phs{���	�
�
�a� � �A��I�M�M�%�(�(�E��)�-�-��"�"�C���z�z��y�}�}�S�1�W�-�-�-��y���s�Q�w��6�6��>�>r6r}Fr[r;r�r>r�)rr-�boxcox1prr)r�rrs  r4�
test_boxcox1pzTestSystematic.test_boxcox1pfse��	?�	?�	?�	�B�K�,�[�9�9� �2�5�9�9�9�3�5�5�A�!� "�!&�	(�	(�	(�	(�	(�	(r6c	��d�}td�t|��tdd��tdd��gd���dS)	Nc��tj|��}tj|tjd��dzz|��tjd|ztjz��z}|jdkr|jS|S�Nr#rr�r�	mpmathifyr�r�r�rrhrQ�rFr3�arg�outs    r4�mp_spherical_jnz9TestSystematic.test_spherical_jn.<locals>.mp_spherical_jny�g���"�1�%�%�C��>�!�f�j��m�m�A�o�"5�s�;�;��;�q��u�V�Y��/�/�0�C��x�1�}�}��x���
r6c�F�tjt|��|��SrN)r-�spherical_jnr��rFr3s  r4rTz2TestSystematic.test_spherical_jn.<locals>.<lambda>�rOr6rr;r�r�r)r2r�r�r}s  r4�test_spherical_jnz TestSystematic.test_spherical_jnxs`��	�	�	�	�C�C�,�_�=�=�#�A�s�^�^�S��s�^�^�<� #�	%�	%�	%�	%�	%�	%r6c��d�}td�t|��tdd��t��g��dS)Nc��tj|��}tj|tjd��dzz|��tjd|ztjz��z}|jdkr|jS|Srwrxrzs    r4r}zATestSystematic.test_spherical_jn_complex.<locals>.mp_spherical_jn�r~r6c�P�tjt|j��|��SrN)r-r�r�rQr�s  r4rTz:TestSystematic.test_spherical_jn_complex.<locals>.<lambda>�������Q�V���a�)H�)H�r6rr;�rrrr
r�s  r4�test_spherical_jn_complexz(TestSystematic.test_spherical_jn_complex��T��	�	�	�	�H�H�,�_�=�=�#�A�s�^�^�Z�\�\�:�	<�	<�	<�	<�	<r6c	��d�}td�t|��tdd��tdd��gd���dS)	Nc��tj|��}tj|tjd��dzz|��tjd|ztjz��z}|jdkr|jS|Srw�rryrr�r�rrhrQrzs    r4�mp_spherical_ynz9TestSystematic.test_spherical_yn.<locals>.mp_spherical_yn�r~r6c�F�tjt|��|��SrN)r-�spherical_ynr�r�s  r4rTz2TestSystematic.test_spherical_yn.<locals>.<lambda>�rOr6rr;g _��rZr<r2r�r�r�s  r4�test_spherical_ynz TestSystematic.test_spherical_yn�sb��	�	�	�	�C�C�,�_�=�=�#�A�s�^�^�S���-=�-=�>� #�	%�	%�	%�	%�	%�	%r6c��d�}td�t|��tdd��t��g��dS)Nc��tj|��}tj|tjd��dzz|��tjd|ztjz��z}|jdkr|jS|Srwr�rzs    r4r�zATestSystematic.test_spherical_yn_complex.<locals>.mp_spherical_yn�r~r6c�P�tjt|j��|��SrN)r-r�r�rQr�s  r4rTz:TestSystematic.test_spherical_yn_complex.<locals>.<lambda>�r�r6rr;r�r�s  r4�test_spherical_yn_complexz(TestSystematic.test_spherical_yn_complex�r�r6c��d�}td�t|��tdd��t��gdd���dS)Nc��tj|��}tj|tjd��dzz|��tjd|ztjz��z}|jdkr|jS|Srw�rryr�r�r�rrhrQrzs    r4�mp_spherical_inz9TestSystematic.test_spherical_in.<locals>.mp_spherical_in�r~r6c�F�tjt|��|��SrN)r-�spherical_inr�r�s  r4rTz2TestSystematic.test_spherical_in.<locals>.<lambda>�rOr6rr;g���訰6r3r�r�r�s  r4�test_spherical_inz TestSystematic.test_spherical_in�s^��	�	�	�	�C�C�,�_�=�=�#�A�s�^�^�S�U�U�3� #�*�	6�	6�	6�	6�	6�	6r6c��d�}td�t|��tdd��t��g��dS)Nc��tj|��}tj|tjd��dzz|��tjd|ztjz��z}|jdkr|jS|Srwr�rzs    r4r�zATestSystematic.test_spherical_in_complex.<locals>.mp_spherical_in�r~r6c�P�tjt|j��|��SrN)r-r�r�rQr�s  r4rTz:TestSystematic.test_spherical_in_complex.<locals>.<lambda>�r�r6rr;r�r�s  r4�test_spherical_in_complexz(TestSystematic.test_spherical_in_complex�r�r6c��d�}td�t|��tdd��t��gd���dS)Nc��tj|tjd��dzz|��tjtjdtj|��zz��z}tj|��jdkr|jS|Srw)rr�r�r�rryrhrQ)rFr3r|s   r4�mp_spherical_knz9TestSystematic.test_spherical_kn.<locals>.mp_spherical_kn�ss���>�!�f�j��m�m�A�o�"5�q�9�9��;�v�y�!�F�,<�Q�,?�,?�*?�@�A�A�B�C����"�"�'�1�,�,��x���
r6c�F�tjt|��|��SrN)r-�spherical_knr�r�s  r4rTz2TestSystematic.test_spherical_kn.<locals>.<lambda>�rOr6rrZr<r2r�r�r�s  r4�test_spherical_knz TestSystematic.test_spherical_kn�s\��	�	�	�	�C�C�,�_�=�=�#�A�s�^�^�S�U�U�3� #�	%�	%�	%�	%�	%�	%r6z.Accuracy issues near z = -1 inherited from kv.c��d�}td�t|��tdd��t��gd���dS)Nc��tj|��}tj|tjd��dzz|��tjd|ztjz��z}|jdkr|jS|Srw)rryr�r�r�rrhrQrzs    r4r�zATestSystematic.test_spherical_kn_complex.<locals>.mp_spherical_kn�r~r6c�P�tjt|j��|��SrN)r-r�r�rQr�s  r4rTz:TestSystematic.test_spherical_kn_complex.<locals>.<lambda>�r�r6rr;r2r�r�s  r4�test_spherical_kn_complexz(TestSystematic.test_spherical_kn_complex�s\��	�	�	�	�H�H�,�_�=�=�#�A�s�^�^�Z�\�\�:� #�	%�	%�	%�	%�	%�	%r6N)��__name__�
__module__�__qualname__r�r�r�r�r�r�r�r�r�r�r�r�r�r�r�rrrrrrrr,r5r>�pytest�mark�xfailrHrNrQr[r_rhrkrorqrtrzr~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�rrr	rrrrr#r'r-r1r9r=rr@rGrKrQrUrZr`rerkrn�xfail_on_32bitrsrwrzr�r�r�r�r�r�r�r�r�r�r�r�r�r�r�r�r�rrrrrrr!r#r%r*r-r2r5r<r?rErKrQrUrXr]rarfrnrtr�r�r�r�r�r�r�r�r�r6r4r�r��s�
������.�.�.�,�,�,�
.�.�.�,�,�,�
.�.�.�,�,�,�
.�.�.�,�,�,�
.�.�.�
.�.�.�
0�0�0�)�)�)�@�@�@�
'�'�'�3�3�3�
6�6�6�6�6�6�@�@�@�

$�
$�
$�%�%�%�	C�	C�	C�(�(�(�>3�3�3�
(�(�(�:7�7�7�
�[���5�)R��S�S�>�>�T�S�>�
7�7�7�
�[���5�)R��S�S�5�5�T�S�5�
U�U�U�(�(�(�.�.�.�'�'�'�?�?�?�F�F�F�
9�9�9�
	7�	7�	7�1�1�1�
(�(�(�&)�)�)�X�X�X�4�4�4�
9�9�9�
(�(�(�'�'�'�*�*�*�
:�:�:�
,�,�,�
5�5�5�
,�,�,�
%�%�%�2�2�2�2�2�2�'�'�'� '�'�'�'�'�'�%�%�%�
3�3�3�
1�1�1�
3�3�3�
,�,�,�
3�3�3�
,�,�,�
e�e�e�
A�A�A�
Q�Q�Q�9�9�9�
1�1�1�%�%�%�%�%�%�%�%�%�
8�8�8�
6�6�6�6�6�6�F�F�F��[���5��!�!�<�<�"�!�<�
.2�.2�.2�`�[���5��!�!�C�C�"�!�C�
�C�C���C�
7�7�7�7�7�7�
�[���5�)P��Q�Q�7�7�R�Q�7�$�$�$�`�`�`�

�

�

��[���5��!�!�$�$�"�!�$��"�"���"��[���5��!�!�3�3�"�!�3�
�[��� \�]�]�(�(�^�]�(�'�'�'�'�'�'��[���5��!�!�=�=�"�!�=�
-�
-�
-�8�8�8�8�8�8�
�,�,���,�
8�8�8�
�[��� <�=�=�6�6�>�=�6�(�(�(�*�,�,���,�
	L�	L�	L�1(�1(�1(�f#�#�#�0#�#�#�0�[���5�)U��V�V�E�E�W�V�E��
#�
#���
#�7�7�7�G�G�G��[���5��!�!�,�,�"�!�,��[���5�)s��t�t�4�4�u�t�4�H�H�H�(�[���5�)n��o�o�6�6�p�o�6�

�
�
�8�8�8�
�[��� '��)�)�%�%�)�)�%�&?�?�?�
F�F�F�U�U�U�(�(�(�4�4�4�
(�(�(�:�:�:�8�8�8�
%�
%�
%�(�(�(�2�2�2�$
�
�
� F�F�F�
;�;�;�
�
�
�>�>�>�(�(�(�"(�(�(�$
%�
%�
%�<�<�<�
%�
%�
%�<�<�<�
6�
6�
6�<�<�<�%�%�%��[���5�)Y��Z�Z�
%�
%�[�Z�
%�
%�
%r6r�)J�__doc__�numpyr$�
numpy.testingrrrr�r��
scipy._libr�
scipy.special�specialr-�scipy.special._testutilsrrr	r
�scipy.special._mptestutilsrrr
rrrrrrr�scipy.special._ufuncsrrrrrrr�ImportErrorr5rHrXrfrzr�r�r�r�r��slowr�r�r�rr
rrrr'r-r1r;rDrMrQrdrkr}r�r�r�r�r�r�r�r6r4�<module>r�sa��������2�2�2�2�2�2�2�2�������
�
�
�
���������������������������������������������������������������������%��M�M�M�M���%�%�%�
�]�8�
$�
$�F�F�F�%������v�v���-�-���-���v�v���	>�	>���	>� ��v�v���!�!���!�(��v�v���6�6���6���v�w���O�O� ��O�,��v�w���C�C� ��C�"��v�v���G�G���G�>��v�v���C�C���C�*��v�v���/�/���/�$��v�v������E�E�����E�8��v�v���:�:���:�*��v�v���0K�0K���0K�n��v�v���:�:���:�L����v�v���
=�
=���
=� ��v�v���=�=���=�&��v�v������;�;�����;�.��v�v������<�<�����<�*��v�v���<�<���<�&��v�v���
<�
<���
<�(��v�v������B�B�����B�&��v�v������;�;�����;�&��v�v���
8�
8���
8���v�v���8�8���8�$��v�v���9�9���9�0��������v�v���?�?�����?�*�����v�v���
@�
@�����
@������v�v���
@�
@�����
@�"�����v�v���=�=�����=��$�s�S�
)�
)�
)�������v�v���D%�D%�D%�D%�D%�D%�D%�����D%�D%�D%s�A#�#A6�5A6

Youez - 2016 - github.com/yon3zu
LinuXploit