403Webshell
Server IP : 217.160.0.135  /  Your IP : 216.73.217.85
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 :  /usr/lib/python3/dist-packages/scipy/integrate/tests/__pycache__/

Upload File :
current_dir [ Writeable ] document_root [ Writeable ]

 

Command :


[ Back ]     

Current File : /usr/lib/python3/dist-packages/scipy/integrate/tests/__pycache__/test_quadrature.cpython-311.pyc
�

d�c�;�� �ddlZddlZddlmZmZmZddlmZmZm	Z	m
Z
mZddlm
Z
mZmZmZmZmZmZmZmZmZmZmZmZddlmZddlmZmZ Gd�d��Z!Gd	�d
��Z"Gd�d��Z#Gd
�d��Z$Gd�d��Z%dS)�N��cos�sin�pi)�assert_equal�assert_almost_equal�assert_allclose�assert_�suppress_warnings)
�
quadrature�romberg�romb�newton_cotes�cumulative_trapezoid�cumtrapz�trapz�	trapezoid�quad�simpson�simps�
fixed_quad�AccuracyWarning)�	_qmc_quad)�stats�specialc��eZdZd�Zd�ZdS)�
TestFixedQuadc�r��d�dd�zz}t�fd�dd����\}}t||d���dS)	N���c���|d�zdz
zS)Nr!r ���x�ns ��G/usr/lib/python3/dist-packages/scipy/integrate/tests/test_quadrature.py�<lambda>z+TestFixedQuad.test_scalar.<locals>.<lambda>s���a�!�A�#��'�l��r�r&��-���q=��rtol)rr	)�self�expected�got�_r&s    @r'�test_scalarzTestFixedQuad.test_scalarsT���
���a��c�7���2�2�2�2�A�q�A�>�>�>���Q���X�E�2�2�2�2�2�2r)c���d}tjdd|z���d�dzz}t�fd�dd|���\}}t||d���dS)	Nrr r!c�"��|�dd�dfzS�Nr#)r%�ps �r'r(z+TestFixedQuad.test_vector.<locals>.<lambda>s���a��1�1�1�d�7��m�r)rr*r+r,)�np�arangerr	)r.r&r/r0r1r6s     @r'�test_vectorzTestFixedQuad.test_vectorsi���
���I�a��1������a�!�e�9���3�3�3�3�Q��Q�?�?�?���Q���X�E�2�2�2�2�2�2r)N)�__name__�
__module__�__qualname__r2r9r#r)r'rrs2������3�3�3�3�3�3�3�3r)rc��eZdZd�Zd�Zd�Zd�Zd�Zd�Zd�Z	d�Z
d	�Zd
�Zd�Z
d�Zd
�Zej�dddg��d���Zd�ZdS)�TestQuadraturec��t�r5)�NotImplementedError)r.r%�a�b�argss     r'rzTestQuadrature.quads��!�!r)c�h�d�}t|dtd��\}}d}t||d���dS)Nc�\�t||z|t|��zz
��tzSr5r�r%r&�zs   r'�myfuncz.TestQuadrature.test_quadrature.<locals>.myfunc$�'���q��s�1�S��V�V�8�|�$�$�R�'�'r)r�r!�������?��ۗ�?���decimal)rrr�r.rH�val�err�	table_vals     r'�test_quadraturezTestQuadrature.test_quadrature"sH��	(�	(�	(��f�a��X�6�6���S�*�	��C��A�6�6�6�6�6�6r)c�l�d�}t|dtdd���\}}d}t||d���dS)Nc�b�dt||z|t|��zz
��ztzS)N��Y)	k�RrrFs   r'rHz3TestQuadrature.test_quadrature_rtol.<locals>.myfunc+s,���#�a��c�!�C��F�F�(�l�+�+�+�B�.�.r)rrJ绽���|�=r,��.�<�R�rrr	rPs     r'�test_quadrature_rtolz#TestQuadrature.test_quadrature_rtol*sM��	/�	/�	/��f�a��X�E�B�B�B���S�1�	���Y�U�3�3�3�3�3�3r)c��d�}d}dD]C}t|dtd|���\}}t||d���t|d	k���DdS)
Nc�\�t||z|t|��zz
��tzSr5rrFs   r'rHz6TestQuadrature.test_quadrature_miniter.<locals>.myfunc3rIr)rL)��4rrJ)�miniterrMrN��?)rrrr
)r.rHrSr`rQrRs      r'�test_quadrature_miniterz&TestQuadrature.test_quadrature_miniter1su��	(�	(�	(�*�	��	�	�G�!�&�!�R��7�K�K�K�H�C����Y��:�:�:�:��C�#�I�����	�	r)c�l�d�}t|dtdd���\}}d}t||d���dS)Nc�b�dt||zdt|��zz
��ztzS)NrWrKrr$s  r'rHz:TestQuadrature.test_quadrature_single_args.<locals>.myfunc<s,���#�a��c�#�c�!�f�f�*�n�-�-�-�b�0�0r)rr!rX�rCr-rYr,rZrPs     r'�test_quadrature_single_argsz*TestQuadrature.test_quadrature_single_args;sM��	1�	1�	1��f�a��!�%�@�@�@���S�1�	���Y�U�3�3�3�3�3�3r)c�d�d�}t|dtd���}d}t||d���dS)Nc�\�t||z|t|��zz
��tzSr5rrFs   r'rHz+TestQuadrature.test_romberg.<locals>.myfuncDrIr)rrJ)rCrLrMrN)r
rr�r.rHrQrSs    r'�test_rombergzTestQuadrature.test_rombergBsG��	(�	(�	(��f�a��(�3�3�3��*�	��C��A�6�6�6�6�6�6r)c�f�d�}t|dtdd���}d}t||d���dS)Nc�b�dt||z|t|��zz
��ztzS)Ng=�`�X�CrrFs   r'rHz0TestQuadrature.test_romberg_rtol.<locals>.myfuncLs,����A�a�C��#�a�&�&��L�)�)�)�"�,�,r)rrJrXreg���3>�Cr,)r
rr	ris    r'�test_romberg_rtolz TestQuadrature.test_romberg_rtolJsI��	-�	-�	-��f�a��(��?�?�?��/�	���Y�U�3�3�3�3�3�3r)c�d�tttjd����d��dS)N��)rrr7r8)r.s r'�	test_rombzTestQuadrature.test_rombRs(���T�"�)�B�-�-�(�(�#�.�.�.�.�.r)c�*�tjd��}tjd|z��}t|��}t	d�|���|�����\}}t||dd���t��5}|�	td��td�|���|���d	�
��}ddd��n#1swxYwYt||dd���dS)Nro皙�����?c�0�tjd|z��S�Nrs�r7r�r%s r'r(z2TestQuadrature.test_romb_gh_3731.<locals>.<lambda>Zs��2�6�#�a�%�=�=�r)g:�0�yE>r)r-�atolzdivmax .4. exceededc�0�tjd|z��Srurvrws r'r(z2TestQuadrature.test_romb_gh_3731.<locals>.<lambda>`s��R�V�C��E�]�]�r)r)�divmaxr+)r7r8rrr�min�maxr	r�filterrr
)r.r%�yrQ�val2rR�sup�val3s        r'�test_romb_gh_3731z TestQuadrature.test_romb_gh_3731Us;���I�f�����F�3�q�5�M�M���1�g�g���0�0�!�%�%�'�'�1�5�5�7�7�C�C�	��c���T��1�5�5�5�5��
 �
 �	P�C��J�J��(=�>�>�>��2�2�A�E�E�G�G�Q�U�U�W�W�Q�O�O�O�D�	P�	P�	P�	P�	P�	P�	P�	P�	P�	P�	P����	P�	P�	P�	P�	��T��A�6�6�6�6�6�6s�AC5�5C9�<C9c�b�ddl}t|jdd��}d}t||d���dS)Nrr g?��k�?rMrN)�mathr
rr)r.r��valmath�expected_vals    r'�test_non_dtypezTestQuadrature.test_non_dtypecs>�������$�(�A�q�)�)��*���G�\�1�=�=�=�=�=�=r)c��d}t|d��\}}t||tjddg��z��t	||dzdz��d}t|d��\}}t	||tjgd���zdz��t	||dzd	z��d}t|d��\}}t	||tjgd
���zdz��t	||dzdz��d
}t|d��\}}t	||tjgd���zdz��t	||dzdz��dS)z5Test the first few degrees, for evenly spaced points.r ��?��(@r!)ra�@ra�@r^g��@)ra�@r�ra� @gP�@r)�@�@@r�r�r�g�V@rMg�=AN)rrr7�arrayr)r.r&�wts�errcoffs    r'�test_newton_cotesz TestQuadrature.test_newton_cotesjsq��
��#�A�q�)�)���W��S�!�B�H�c�3�Z�0�0�0�1�1�1��G�a��d�U�4�Z�0�0�0�
��#�A�q�)�)���W��C��2�8�O�O�O�#<�#<�!<�S�!@�A�A�A��G�a��d�U�6�\�2�2�2�
��#�A�q�)�)���W��C��2�8�,@�,@�,@�#A�#A�!A�#�!E�F�F�F��G�a��d�U�6�\�2�2�2�
��#�A�q�)�)���W��C��2�8�,H�,H�,H�#I�#I�!I�$�!N�O�O�O��G�a��d�U�9�_�5�5�5�5�5r)c�V�tjgd���}|dz}t|��\}}d}tj||��}t	||��tjgd���}|dz}t|��\}}d}tj||��}t	||��dS)z9Test newton_cotes with points that are not evenly spaced.)�g�?�@r!gUUUUUU@)r�gffffff�?g������@r��"@N)r7r�r�dotr)r.r%r~r�r��exact_integral�numeric_integrals       r'�test_newton_cotes2z!TestQuadrature.test_newton_cotes2�s���
�H�_�_�_�%�%��
�q�D��#�A�����W����6�#�q�>�>���,�n�=�=�=��H�)�)�)�*�*��
�q�D��#�A�����W����6�#�q�>�>���,�n�=�=�=�=�=r)c
��tjd��}tt|��d��tt|d���d��tt|tjddd�����d	��tjd��}d
|z}tt||d���d
��tt||d���d��tt||d���d��tjdg��}tj|d
��}tt||d���d��tt||d���d��tjgd���}tj|d
��}tt||d���d��tt||d���d��tjgd�gd�gd�g��}tj|d
��}gd�}gd�}tt||d���|��tt||d���|��tjgd�gd�gd�g��}tj|d
��}gd�}gd�}tt||d���|��tt||d���|��dS)Nrorpr���dx�@rrrw� r!�avg)r%�eveng�+@�firstg�+@�last�r�)r%�axisr����)r�r�r�r�)r r!r�)r�r�r�r�)��e@r�r�)r r��r�)r�ga@g�@g�@)r�r�g �@)r7r8rr�linspacer��power)r.r~r%�	zero_axis�default_axiss     r'�test_simpsonzTestQuadrature.test_simpson�s����I�b�M�M���W�Q�Z�Z��%�%�%��W�Q�3�'�'�'��,�,�,��W�Q�"�+�a��B�"7�"7�8�8�8�"�=�=�=��I�a�L�L��
�q�D���W�Q�!�%�0�0�0�&�9�9�9��W�Q�!�'�2�2�2�E�:�:�:��W�Q�!�&�1�1�1�2�6�6�6�
�H�a�S�M�M���H�Q��N�N���W�Q�!�!�,�,�,�c�2�2�2��W�Q�!�"�-�-�-�s�3�3�3��H�\�\�\�"�"���H�Q��N�N���W�Q�!�!�,�,�,�c�2�2�2��W�Q�!�"�-�-�-�s�3�3�3��H�l�l�l�L�L�L�,�,�,�?�@�@���H�Q��N�N��(�(�(�	�/�/�/���W�Q�!�!�,�,�,�i�8�8�8��W�Q�!�"�-�-�-�|�<�<�<��H�l�l�l�L�L�L�.�.�.�A�B�B���H�Q��N�N��0�0�0�	�1�1�1���W�Q�!�!�,�,�,�i�8�8�8��W�Q�!�"�-�-�-�|�<�<�<�<�<r)�droplastFTc��tjgd�gd�g��}|r|dd�dd�f}t|d���}ttj|tj���d���}t	||��dS)N)r!r!rrr�r����r^)rrr!r��
����r�r�)r�)�dtype)r7r�r�float64r)r.r�r~�resultr/s     r'�test_simpson_2d_integer_no_xz+TestQuadrature.test_simpson_2d_integer_no_x�s���
�H�/�/�/�3�3�3�5�
6�
6���	��!�!�!�S�b�S�&�	�A����$�$�$���2�8�A�R�Z�8�8�8�r�B�B�B���V�X�&�&�&�&�&r)c	��tjd��}d|z}tt||dd���t	||dd�����dS)Nrr!r�r�)r%r�r�)r7r8rrr�r.r~r%s   r'�
test_simpszTestQuadrature.test_simps�s[���I�a�L�L��
�q�D���W�Q�!��'�:�:�:��1��c��8�8�8�	:�	:�	:�	:�	:r)N)r:r;r<rrTr[rbrfrjrmrqr�r�r�r�r��pytest�mark�parametrizer�r�r#r)r'r>r>s������"�"�"�7�7�7�4�4�4����4�4�4�7�7�7�4�4�4�/�/�/�7�7�7�>�>�>�6�6�6�,>�>�>�"#=�#=�#=�J�[���Z�%���7�7�	'�	'�8�7�	'�:�:�:�:�:r)r>c�,�eZdZd�Zd�Zd�Zd�Zd�ZdS)�TestCumulative_trapezoidc���tjddd���}|}t||d���}gd�}t||��t||d���}t||dd���dS)	Nr�r!r^��numr��initial)r������r�r�r �r7r�rr	)r.r%r~�y_int�
y_expecteds     r'�test_1dz TestCumulative_trapezoid.test_1d�s~���K��A�1�%�%�%��
��$�Q��1�5�5�5��.�.�.�
���z�*�*�*�$�Q��4�8�8�8����z�!�"�"�~�.�.�.�.�.r)c���tjd���ddd��}|}t||d���}tjgd�gd�ggd	�gd
�ggd�gd�gg��}t||��gd
�}t
gd�|��D]U\}}t||d|���}t|jd��t||d|���}t|j|���VdS)N�r�r!rrr�)r�r�r��@)r�r��$@��0@)r��!@�2@g�<@)r��)@g:@g@D@)r�r�gA@g@J@)r���4@gE@g P@))r!r!r)r�r r)r�r!r��rr r!g������@)r�r�)r�r!r)	r7r8�reshaperr�r	�zipr�shape)r.r%r~r�r��shapesr�r�s        r'�test_y_nd_x_ndz'TestCumulative_trapezoid.test_y_nd_x_nd�s4���I�i� � �(�(��A�q�1�1��
��$�Q��1�5�5�5���X� 2� 2� 2� 4� 4� 4� 6� 4� 4� 4� 5� 5� 5� 7� 5� 5� 5� 5� 5� 5� 7�	8�9�9�
�	��z�*�*�*�3�2�2���y�y�y�&�1�1�	-�	-�K�D�%�(��A�t�$�G�G�G�E����i�0�0�0�(��A�t�$�G�G�G�E����e�,�,�,�,�		-�	-r)c	��tjd���ddd��}tjd��dz}tjgd�gd�ggd�gd�gg��tjgd	�ggd
�ggd�gg��tjgd�gd
�ggd�gd�ggd�gd�gg��f}t	gd�|��D];\}}t||d|j|�|d���}t||���<dS)Nr�r�r!r)r��@r�r�)r�r�r��&@)gD@gF@gH@gJ@)gL@gN@gP@gQ@)r�r�r�r�)r�r�r�g*@)r�g3@g4@�5@)r�r�g�1@)r�r�g�J@)r�g�B@g`V@)r�g�J@g`_@)r�g@Q@g0d@)r�g@U@g�h@r�)r%r�r�)r7r8r�r�r�rr�r	)r.r~r%�ys_expectedr�r�r�s       r'�test_y_nd_x_1dz'TestCumulative_trapezoid.test_y_nd_x_1d�s~���I�i� � �(�(��A�q�1�1���I�a�L�L�!�O��
�H�'�'�'�)�)�)�+�+�+�+�+�+�+�-�.�
/�
/�
�H�'�'�'�(�+�+�+�,�+�+�+�,�.�
/�
/�
�H����'�'�'�)�'�'�'�)�)�)�+�)�)�)�)�)�)�+�	,�
-�
-�
.��!$�I�I�I�{� ;� ;�	/�	/��D�*�(��a�������.?�d�15�7�7�7�E��E�:�.�.�.�.�	/�	/r)c�Z�tjddd���}t|��}gd�}t||��t|d���}gd�}t||��t|d	�
��}gd�}t||��t|d	d���}gd
�}t||��dS)Nr�r!r^r�)r�r�r�r��G�z��?r�)r�r�r�r�r�r�r�)����r�r�)r�r�)r�r�r�r�r�r�)r.r~r�r�s    r'�test_x_nonez$TestCumulative_trapezoid.test_x_nones����K��A�1�%�%�%��$�Q�'�'��*�*�*�
���z�*�*�*�$�Q��5�5�5��0�0�0�
���z�*�*�*�$�Q�1�-�-�-��*�*�*�
���z�*�*�*�$�Q�1�d�;�;�;��0�0�0�
���z�*�*�*�*�*r)c
���tjd���ddd��}|}tt	||ddd���t||ddd���d�	��dS)
Nr�r�r!rr�r)r�r�r�g�+����=r,)r7r8r�r	rr)r.r%r~s   r'�
test_cumtrapzz&TestCumulative_trapezoid.test_cumtrapzsx���I�i� � �(�(��A�q�1�1��
���,�Q��c��1�M�M�M� ��A�#�A�q�A�A�A�"�	$�	$�	$�	$�	$�	$r)N)r:r;r<r�r�r�r�r�r#r)r'r�r��s_������/�/�/�-�-�-�*/�/�/�0+�+�+�&$�$�$�$�$r)r�c��eZdZdZd�Zd�ZdS)�
TestTrapezoidz[This function is tested in NumPy more extensive, just do some
    basic due diligence here.c
�~�tjd��}tt|��d��tt|d���d��tt|tjddd�����d	��tjd��}d
|z}tt||d���d
��dS)Nrorpr�r�r�rrrwr�r!g�������?)r%r�g+@)r7r8rrr�r�s   r'�test_trapezoidzTestTrapezoid.test_trapezoid!s����I�b�M�M���Y�q�\�\�3�'�'�'��Y�q�S�)�)�)�2�.�.�.��Y�q�B�K��1�b�$9�$9�:�:�:�B�?�?�?��I�a�L�L��
�q�D���Y�q�A�#�.�.�.��5�5�5�5�5r)c	��tjd��}d|z}tt||dd���t	||dd�����dS)Nrr!r�r)r%r�r�)r7r8rrrr�s   r'�
test_trapzzTestTrapezoid.test_trapz+s[���I�a�L�L��
�q�D���Y�q�A�#�A�6�6�6��1��c��2�2�2�	4�	4�	4�	4�	4r)N)r:r;r<�__doc__r�r�r#r)r'r�r�s<������!�!�6�6�6�4�4�4�4�4r)r�c��eZdZd�Zddejd��fd�Zej�	dddg��ej�	ddd	g��d
�����Z
ej�	dddgd
d
gd
dgdd
gg��d���Zej�	dddg��d���Zd�Z
dS)�TestQMCQuadc	��d}tjt|���5tdddgddg��ddd��n#1swxYwYd}tjt|���5td�ddgddg��ddd��n#1swxYwYd�}d	}tjt|���5t|ddgddg��ddd��n#1swxYwYd
}tjt|���5td�ddgddgd�
��ddd��n#1swxYwYd}tjt|���5td�ddgddgd���ddd��n#1swxYwYd}tjt|���5td�ddgddgd���ddd��n#1swxYwYd}tjt|���5td�ddgddgtj�	d�����ddd��n#1swxYwYd}tjt|���5td�ddgddgd���ddd��dS#1swxYwYdS)Nz`func` must be callable.��matchza duckrr z/`func` must evaluate the integrand at points...c��dS�Nr r#r#r)r'r(z3TestQMCQuad.test_input_validation.<locals>.<lambda>;s��Q�r)c�D�|jdksJ�tj|��Sr�)�ndimr7�sumrws r'�funcz/TestQMCQuad.test_input_validation.<locals>.func=s���6�Q�;�;�;�;��6�!�9�9�r)z8Exception encountered when attempting vectorized call...z`n_points` must be an integer.c��dSr�r#rws r'r(z3TestQMCQuad.test_input_validation.<locals>.<lambda>F���q�r)g�@)�n_pointsz!`n_estimates` must be an integer.c��dSr�r#rws r'r(z3TestQMCQuad.test_input_validation.<locals>.<lambda>Jr�r)r�)�n_estimatesz8`qrng` must be an instance of scipy.stats.qmc.QMCEngine.c��dSr�r#rws r'r(z3TestQMCQuad.test_input_validation.<locals>.<lambda>Nr�r))�qrngz8`qrng` must be initialized with dimensionality equal to c��dSr�r#rws r'r(z3TestQMCQuad.test_input_validation.<locals>.<lambda>Rr�r)z,`log` must be boolean \(`True` or `False`\).c��dSr�r#rws r'r(z3TestQMCQuad.test_input_validation.<locals>.<lambda>Vr�r)r���log)
r��raises�	TypeError�qmc_quad�
ValueError�warns�UserWarningr�qmc�Sobol)r.�messager�s   r'�test_input_validationz!TestQMCQuad.test_input_validation4s���,��
�]�9�G�
4�
4�
4�	/�	/��X��1�v��1�v�.�.�.�	/�	/�	/�	/�	/�	/�	/�	/�	/�	/�	/����	/�	/�	/�	/�D��
�]�:�W�
5�
5�
5�	0�	0��Y�Y��A���A��/�/�/�	0�	0�	0�	0�	0�	0�	0�	0�	0�	0�	0����	0�	0�	0�	0�	�	�	�M��
�\�+�W�
5�
5�
5�	+�	+��T�A�q�6�A�q�6�*�*�*�	+�	+�	+�	+�	+�	+�	+�	+�	+�	+�	+����	+�	+�	+�	+�3��
�]�9�G�
4�
4�
4�	C�	C��[�[�1�a�&�1�a�&�6�B�B�B�B�	C�	C�	C�	C�	C�	C�	C�	C�	C�	C�	C����	C�	C�	C�	C�6��
�]�9�G�
4�
4�
4�	C�	C��[�[�1�a�&�1�a�&�c�B�B�B�B�	C�	C�	C�	C�	C�	C�	C�	C�	C�	C�	C����	C�	C�	C�	C�M��
�]�9�G�
4�
4�
4�	A�	A��[�[�1�a�&�1�a�&�x�@�@�@�@�	A�	A�	A�	A�	A�	A�	A�	A�	A�	A�	A����	A�	A�	A�	A�M��
�]�:�W�
5�
5�
5�	K�	K��[�[�1�a�&�1�a�&�u�y���q�7I�7I�J�J�J�J�	K�	K�	K�	K�	K�	K�	K�	K�	K�	K�	K����	K�	K�	K�	K�B��
�]�9�G�
4�
4�
4�	:�	:��[�[�1�a�&�1�a�&�b�9�9�9�9�	:�	:�	:�	:�	:�	:�	:�	:�	:�	:�	:�	:����	:�	:�	:�	:�	:�	:s��A�A�A�(B�B�B�6C�C�C�D%�%D)�,D)�
E2�2E6�9E6�F?�?G�G�'6H)�)H-�0H-�I7�7I;�>I;�r�r!c
��
���d}tj|���tj|���
�
�fd��tj�d��}t
j�||���}tj|��}tj|��|z}t�||||��
f|���}	t
j
�|��
|���}
tj
|dz
d��|	jz}t|	j|
|�	��tj|��|	jzd
ksJ�tj�d��}t
j�||���}t�fd�||||��
fd|�
��}ttj|j��|	j��tj|j��tj|��d
krtjnd
ksJ�dS)Nr!c�F��tj�|����Sr5)r�multivariate_normal�pdf)r%�cov�means ��r'r�z$TestQMCQuad.basic_test.<locals>.func^s����,�0�0��D�#�>�>�>r)l�^^@E9�?)�seed)r�rrCr)�lower_limitr gףp=
��?)rxrc�2��tj�|���Sr5)r7r)rCr�s �r'r(z(TestQMCQuad.basic_test.<locals>.<lambda>ns�����t�t�T�{�(;�(;�r)T)r�rrCrr)r7�zeros�eye�random�default_rngrrr�onesr
r�cdf�sc�stdtrit�standard_errorr	�integral�prod�exp�imagr)r.r�r�signsr��rngrrArB�res�refrx�logresrr�rs             @@@r'�
basic_testzTestQMCQuad.basic_testXs���������x��~�~���f�T�l�l��	?�	?�	?�	?�	?�	?��i�#�#�$7�8�8���y���t�#��.�.���H�T�N�N���G�D�M�M�E�!���t�Q��H�#.�d�C�[�t�M�M�M���'�+�+�A�t�S�a�+�H�H���z�+�a�-��/�/�#�2D�D�����c��5�5�5�5��w�u�~�~�c�l�*�Q�.�.�.�.��i�#�#�$7�8�8���y���t�#��.�.���;�;�;�;�Q��#+�� $�c�{��4�A�A�A��	���v��/�/���>�>�>��w�v��'�'�R�W�U�^�^�a�5G�5G�B�E�E�Q�O�O�O�O�O�Or)r�irr�c�2�|�||��dSr5�r.)r.r�rs   r'�
test_basiczTestQMCQuad.test_basicts��	
����+�.�.�.�.�.r)r)r r�c�2�|�|���dS)N)r)r0)r.r)s  r'�	test_signzTestQMCQuad.test_signys�����e��$�$�$�$�$r)rFTc���d}tjt|���5td�ddgddg|���}ddd��n#1swxYwY|j|r
t
jndksJ�|jdksJ�dS)Nz-A lower limit was equal to an upper limit, sor�c��dSr�r#rws r'r(z'TestQMCQuad.test_zero.<locals>.<lambda>�s��Q�r)rr r)r�rr
r
r%r7�infr$)r.rrr+s    r'�	test_zerozTestQMCQuad.test_zero}s���A��
�\�+�W�
5�
5�
5�	A�	A��;�;��A���A��C�@�@�@�C�	A�	A�	A�	A�	A�	A�	A�	A�	A�	A�	A����	A�	A�	A�	A��|�3� 5�����A�6�6�6�6��!�Q�&�&�&�&�&�&s�A�A�
Ac���d�}t|dd��}tj�dd���tj�dd���z
}t	|j|d��dS)Nc�D�tj�|d���S)Nr!��scale)r�normrrws r'r�z-TestQMCQuad.test_flexible_input.<locals>.func�s���:�>�>�!�1�>�-�-�-r)rr r!r:g{�G�z�?)r
rr<r!r	r%)r.r�r+r,s    r'�test_flexible_inputzTestQMCQuad.test_flexible_input�sl��	.�	.�	.��t�Q��"�"���j�n�n�Q�a�n�(�(�5�:�>�>�!�1�>�+E�+E�E�����c�4�0�0�0�0�0r)N)r:r;r<rr7r r.r�r�r�r1r3r7r=r#r)r'r�r�3s+������":�":�":�H#'�A�W�R�W�Q�Z�Z�P�P�P�P�8�[���Z�$���7�7��[���]�Q��G�4�4�/�/�5�4�8�7�/��[���W��1�v��B�x�"�a��1�b�'�&J�K�K�%�%�L�K�%��[���U�U�D�M�2�2�'�'�3�2�'�1�1�1�1�1r)r�)&r��numpyr7rrr�
numpy.testingrrr	r
r�scipy.integraterr
rrrrrrrrrrr�scipy.integrate._quadraturerr
�scipyrrr"rr>r�r�r�r#r)r'�<module>rCsL��
�
�
�
���������������7�7�7�7�7�7�7�7�7�7�7�7�7�7�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�>�=�=�=�=�=�&�&�&�&�&�&�&�&�
3�
3�
3�
3�
3�
3�
3�
3� i:�i:�i:�i:�i:�i:�i:�i:�XQ$�Q$�Q$�Q$�Q$�Q$�Q$�Q$�h4�4�4�4�4�4�4�4�*Z1�Z1�Z1�Z1�Z1�Z1�Z1�Z1�Z1�Z1r)

Youez - 2016 - github.com/yon3zu
LinuXploit