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 :  /lib/python3/dist-packages/scipy/stats/tests/__pycache__/

Upload File :
current_dir [ Writeable ] document_root [ Writeable ]

 

Command :


[ Back ]     

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

d�ckO��x�ddlmZmZmZddlZddlmZmZm	Z	m
Z
mZmZddl
Z
ddl
mZd�Zd�Ze
jjd���Ze
jjd���Zd	�Zd
�ZGd�dej��ZGd
�dej��ZGd�dej��Zd�Zd�Zd�Zd�Zgd�Z e
j�!de ddgz��e
j�!de ��d�����Z"d�Z#d�Z$d�Z%d�Z&e
jj'd���Z(d �Z)e
jj'd!���Z*d"�Z+d#�Z,d$�Z-d%�Z.d&�Z/dS)'�)�stats�linalg�	integrateN)�assert_almost_equal�assert_�assert_equal�assert_array_almost_equal�assert_array_almost_equal_nulp�assert_allclose)�raisesc��tj�d��d}tj�|��}|���}|�d���}t
j|��}tjddd��}|�	|��}tj
�|||���}|d|d	z
}ttj
||z
d
z��|zdk��|�|tj��}	|�tj|��}
t!|	dd�
��t!|
dd�
��t!|�|tj��|	d�
��t!|�tj|��|
d�
��t!|�|��|d
z�
��|zd
�
��t!|�||d
z��||z�
��|zd
�
��dS)N�����)�ddof�������loc�scaler��{�G�z�?��?��decimal�
)�np�random�seed�randn�mean�stdr�gaussian_kde�linspace�evaluate�norm�pdfr�sum�integrate_box_1d�infr�
integrate_box�
integrate_kde�integrate_gaussian)�n_basesample�xn�xnmean�xnstd�gkde�xs�kdepdf�normpdf�	intervall�prob1�prob2s           �?/usr/lib/python3/dist-packages/scipy/stats/tests/test_kdeoth.py�test_kde_1dr;	s���I�N�N�7�����L�	�����	&�	&�B�
�W�W�Y�Y�F��F�F��F�N�N�E���b�!�!�D�
��R��#�	�	�B�
�]�]�2�
�
�F��j�n�n�R�V�5�n�9�9�G��1���1��
�I��B�F�F�W�$�q�(�)�)�)�3�d�:�;�;�;��!�!�&�"�&�1�1�E��!�!�2�6�'�6�2�2�E���s�A�.�.�.�.���s�A�.�.�.�.���*�*�6�2�6�:�:�E�2�N�N�N�N���*�*�B�F�7�F�;�;�U�B�O�O�O�O���*�*�4�0�0������)�)�)�3�Q�@�@�@�@���/�/���q��A�A����,�,�.�.�y�8�!�E�E�E�E�E�E�c��tj�d��d}tj�|��}tj�|��}tj||���}tjtj||z
dz|�����}tj||���}tj	ddd��}|�
|��}tj�|||���}|d	|d
z
}	ttj||z
dz��|	zdk��|�|tj��}
|�tj|��}t#|
dd	�
��t#|dd	�
��t#|�|tj��|
d�
��t#|�tj|��|d�
��t#|�|��|dz���|	zd�
��t#|�||dz��||z���|	zd�
��dS)Nrr��weightsrrrrrrrrrrr)rrr r!�rand�average�sqrtrr$r%r&r'r(rr)r*r+rr,r-r.)r/r0�wnr1r2r3r4r5r6r7r8r9s            r:�test_kde_1d_weightedrD(sR���I�N�N�7�����L�	�����	&�	&�B�	�����	%�	%�B�
�Z��B�
'�
'�
'�F��G�B�J��6�	�A�~�r�:�:�:�;�;�E���b�"�-�-�-�D�
��R��#�	�	�B�
�]�]�2�
�
�F��j�n�n�R�V�5�n�9�9�G��1���1��
�I��B�F�F�W�$�q�(�)�)�)�3�d�:�;�;�;��!�!�&�"�&�1�1�E��!�!�2�6�'�6�2�2�E���s�A�.�.�.�.���s�A�.�.�.�.���*�*�6�2�6�:�:�E�2�N�N�N�N���*�*�B�F�7�F�;�;�U�B�O�O�O�O���*�*�4�0�0������)�)�)�3�Q�@�@�@�@���/�/���q��A�A����,�,�.�.�y�8�!�E�E�E�E�E�Er<c�,�tj�d��d}tjddg��}tjddgddgg��}tj�|||���j}t
j|��}tjdd	d
�dd	d
�f\}}tj	|�
��|�
��g��}|�|��}|�dd��}tj�
tj||g��||���}	|�
��d|�
��d
z
}
ttj||	z
dz��|
dzzdk��d}d}|�||dg||g��}
|�||g||dg��}t%|
dd���t%|dd���t%|�|��|dz���|
dzzd���t%|�||��||	z���|
dzzd���dS)Nrr��?�@�@�@��sizerr�@@�r"�covrrrr�}Ô%�I���}Ô%�I�Trr)rrr �array�multivariate_normal�Trr$�mgrid�vstack�ravelr&�reshaper(�dstackrr)r,rr-r.)r/r"�
covariancer0r3�x�y�grid_coordsr5r6r7�small�larger8r9s               r:�test_kde_2dr_Hs}���I�N�N�7�����L�
�8�S�#�J���D���C��:��S�z�2�3�3�J�
��	&�	&�t�Z�l�	&�	K�	K�	M�B���b�!�!�D��8�B�q��I�r�!�D�y�(�)�D�A�q��)�Q�W�W�Y�Y����	�	�2�3�3�K�
�]�]�;�
'�
'�F�
�^�^�C��
%�
%�F��'�+�+�B�I�q�!�f�,=�,=�D�j�+�Y�Y�G����	�	�!��q�w�w�y�y��|�+�I��B�F�F�W�$�q�(�)�)�Y��\�:�T�A�B�B�B��E��E�����t�A�w�/�%���@�@�E�����u�~��t�A�w�/?�@�@�E���s�A�.�.�.�.���s�A�.�.�.�.���*�*�4�0�0������)�)�9�a�<�8�!�E�E�E�E���/�/��j�A�A����,�,�.�.�	�1��=�q�J�J�J�J�J�Jr<c�n�tj�d��d}tjddg��}tjddgddgg��}tj�|||���j}tj�|��}tj||���}tj	d	d
d�d	d
d�f\}}tj
|���|���g��}|�|��}	|	�
dd��}	tj�tj||g��||���}
|���d
|���dz
}t!tj|	|
z
dz��|dzzdk��d}d}
|�||d
g|
|
g��}|�||g|
|d
g��}t'|dd
���t'|dd
���t'|�|��|	dz���|dzzd���t'|�||��|	|
z���|dzzd���dS)NrrrFrGrHrIrJr>rrrLrMrrrrrOrPrr)rrr rQrRrSr@rr$rTrUrVr&rWr(rXrr)r,rr-r.)r/r"rYr0rCr3rZr[r\r5r6r7r]r^r8r9s                r:�test_kde_2d_weightedraos����I�N�N�7�����L�
�8�S�#�J���D���C��:��S�z�2�3�3�J�
��	&�	&�t�Z�l�	&�	K�	K�	M�B�	�����	%�	%�B���b�"�-�-�-�D��8�B�q��I�r�!�D�y�(�)�D�A�q��)�Q�W�W�Y�Y����	�	�2�3�3�K�
�]�]�;�
'�
'�F�
�^�^�C��
%�
%�F��'�+�+�B�I�q�!�f�,=�,=�D�j�+�Y�Y�G����	�	�!��q�w�w�y�y��|�+�I��B�F�F�W�$�q�(�)�)�Y��\�:�T�A�B�B�B��E��E�����t�A�w�/�%���@�@�E�����u�~��t�A�w�/?�@�@�E���s�A�.�.�.�.���s�A�.�.�.�.���*�*�4�0�0������)�)�9�a�<�8�!�E�E�E�E���/�/��j�A�A����,�,�.�.�	�1��=�q�J�J�J�J�J�Jr<c�D�d�}tj�d��d}tj�|��}t	j|��}t	j||���}t	j||j���}tjddd��}|�|��}|�|��}t||��|�|��}	t||	��tttj|d���dS)	Nc�L�tj|jd|jdzz��S�z*Same as default, just check that it works.g��)r�power�n�d��kde_objs r:�
scotts_factorz0test_kde_bandwidth_method.<locals>.scotts_factor�s!���x��	�3��	�!��#4�5�5�5r<r�2��	bw_methodrr�3�wrongstring�rrr r!rr$�factorr%r&r�
assert_raises�
ValueError�
rkr/r0r3�gkde2�gkde3r4r5�kdepdf2�kdepdf3s
          r:�test_kde_bandwidth_methodrz�s��6�6�6��I�N�N�7�����L�	�����	&�	&�B���b�!�!�D���r�]�;�;�;�E���r�T�[�9�9�9�E�	��R��"�	�	�B�
�]�]�2�
�
�F��n�n�R� � �G����(�(�(��n�n�R� � �G����(�(�(��*�e�0�"�
�N�N�N�N�N�Nr<c�D�d�}tj�d��d}tj�|��}t	j|��}t	j||���}t	j||j���}tjddd��}|�|��}|�|��}t||��|�|��}	t||	��tttj|d���dS)	Nc�L�tj|jd|jdzz��Srd)rrf�neffrhris r:rkz9test_kde_bandwidth_method_weighted.<locals>.scotts_factor�s!���x���c�7�9�Q�;�&7�8�8�8r<rrlrmrrrorprqrus
          r:�"test_kde_bandwidth_method_weightedr~�s��9�9�9��I�N�N�7�����L�	�����	&�	&�B���b�!�!�D���r�]�;�;�;�E���r�T�[�9�9�9�E�	��R��"�	�	�B�
�]�]�2�
�
�F��n�n�R� � �G����(�(�(��n�n�R� � �G����(�(�(��*�e�0�"�
�N�N�N�N�N�Nr<c��eZdZd�ZdS)�_kde_subclass1c��tj|��|_|jj\|_|_|j|_|���dS�N)	r�
atleast_2d�dataset�shaperhrgrk�covariance_factor�_compute_covariance)�selfr�s  r:�__init__z_kde_subclass1.__init__�sH���}�W�-�-�����+������!%�!3���� � �"�"�"�"�"r<N)�__name__�
__module__�__qualname__r��r<r:r�r��s#������#�#�#�#�#r<r�c���eZdZ�fd�Z�xZS)�_kde_subclass2c�b��|j|_t���|��dSr�)rkr��superr�)r�r��	__class__s  �r:r�z_kde_subclass2.__init__�s,���!%�!3���
������!�!�!�!�!r<)r�r�r�r��
__classcell__)r�s@r:r�r��s8�������"�"�"�"�"�"�"�"�"r<r�c��eZdZd�ZdS)�_kde_subclass4c�0�d|���zS)Nr)�silverman_factor)r�s r:r�z _kde_subclass4.covariance_factor�s���T�*�*�,�,�,�,r<N)r�r�r�r�r�r<r:r�r��s#������-�-�-�-�-r<r�c�8��tjgd�t���}tjddd���}t	j|����|��}t
|��}||��}t||d���t|��}||��}t||d���t|��}||��}gd�}	t|	|d	�
���}
�fd�|
_|
���|
|��}t||d���dS)N�r���rre���dtype����
rl��num��nulp)g�n�,�?gǻe�?�?gW�[3�?g�a��)�?gKd8�;�?�rc����jSr�)rr)�kdes�r:�<lambda>z/test_gaussian_kde_subclassing.<locals>.<lambda>s	���S�Z�r<)
rrQ�floatr%rr$r�r
r�r�r	r�r�)
�x1r4�ys�kde1�y1�kde2�y2�kde4�y4�
y_expected�kde5�y5r�s
            @r:�test_gaussian_kde_subclassingr��sC���	��#�#�#�5�	1�	1�	1�B�	��S�"�"�	%�	%�	%�B��
�R�
 �
 �C�	��R���B��"���D�	
��b���B�"�2�r��3�3�3�3��"���D�	
��b���B�"�2�r��3�3�3�3��"���D�	
��b���B�M�M�M�J��j�"�a�8�8�8�8��D�/�/�/�/�D��������	
��b���B�"�2�r��3�3�3�3�3�3r<c�0�tjgd�t���}tjddd���}gd�}t	j|��}|�d�	��|�d
�	��||��}t||d���dS)
Nr�r�r�r�r�r�)g>L\�9�?gp4=�?g�H��꟫?g8(H	MT�?gIJ�CR�?rrm�scottrr)rrQr�r%rr$�
set_bandwidthr	)r�r4r�r�r�s     r:�$test_gaussian_kde_covariance_cachingr�s���	��#�#�#�5�	1�	1�	1�B�	��S�"�!�	$�	$�	$�B�N�M�M�J��
�R�
 �
 �C������$�$�$������(�(�(�	��R���B��j�"�a�8�8�8�8�8�8r<c�R�tjgd�t���}tjddd���}t	j|��}|j|_|���||��}t	j|d���}||��}t||d�	��d
S)z�Ugly, but people may rely on this.  See scipy pull request 123,
    specifically the linked ML thread "Width of the Gaussian in stats.kde".
    If it is necessary to break this later on, that is to be discussed on ML.
    r�r�r�r�rlr��	silvermanrmr�N)
rrQr�r%rr$r�r�r�r
)r�r4r�r�r�r�s      r:�test_gaussian_kde_monkeypatchr�s���

��#�#�#�5�	1�	1�	1�B�	��S�"�"�	%�	%�	%�B��
�R�
 �
 �C��0�C��������	��R���B���b�K�8�8�8�D�	
��b���B�"�2�r��3�3�3�3�3�3r<c��tjd��}tj|��}gd�}t	||��|d���dS)zRegression test for #1181.r�)�B��\A�?筱�ES�?g�
f���?r�r�r�rN)r�arangerr$r	)r�r�r�s   r:�test_kde_integer_inputr�+sM��	��1���B�
�
�R�
 �
 �C�M�M�M�J��c�c�"�g�g�z�1�=�=�=�=�=�=r<)�float32�float64�float96�float128�int32�int64�bw_typer�r�r�c� �tt|d��}|dvr|}n%tt|d��}|r|d��nd}td�||fD����rtj��tjd|���}tjd|���}t
j|||���}tjd|���}||��}|jtj	||tj
|��|j��ksJ�dS)N)r�r��c3�K�|]}|duV��	dSr�r�)�.0�dts  r:�	<genexpr>z(test_kde_output_dtype.<locals>.<genexpr>As&����
,�
,�"�2��:�
,�
,�
,�
,�
,�
,r<r�r�)rnr?)�getattrr�any�pytest�skipr�rr$r��result_typer�rr)r�r��bwr?r��k�points�results        r:�test_kde_output_dtyper�5s#��
�B��t�$�$�E��(�(�(�
����"�g�t�,�,��"�
,�W�W�Q�Z�Z�Z���
�
,�
,����
,�
,�
,�,�,����
�
�
��i���'�'�'�G��i���'�'�'�G�
��7�b�'�B�B�B�A�
�Y�q��
&�
&�
&�F�
�Q�v�Y�Y�F��<�2�>�'�6�2�:�g�;N�;N�*+�(�4�4�4�4�4�4�4�4r<c�V�tj�d��}|�d��}t	j|��}|�d��}d}t
jt|���5|�	|��ddd��dS#1swxYwYdS)Nl	�F�g�u	\�U.�lt&0)rr�)r�r�z0points have dimension 3, dataset has dimension 2��match)
rr�default_rng�standard_normalrr$r�rrt�logpdf)�rngr0r3r4�msgs     r:�test_pdf_logpdf_validationr�Ns���
�)�
�
� F�
G�
G�C�	�	�	�W�	%�	%�B���b�!�!�D�	�	�	�W�	%�	%�B�
<�C�	��z��	-�	-�	-������B������������������������s�;B�B"�%B"c��tj�d��d}tj�|��}t	j|��}tjddd��}|�|��}|�|��}t||d���tj
|��}|�|��}t||d���t	j|��}tj
|�|����}|�|��}t||d���dS)Nrrl����r)rrr r!rr$r%r&r(r�logr�)r/r0r3r4r(�pdf2r��logpdf2s        r:�test_pdf_logpdfr�Ys���I�N�N�1�����L�	�����	&�	&�B���b�!�!�D�	��S�"�b�	!�	!�B�
�-�-��
�
�C��8�8�B�<�<�D���T�2�.�.�.�.�
�V�C�[�[�F��k�k�"�o�o�G�����4�4�4�4���b�!�!�D�
�&����r�"�"�
#�
#�C��;�;�r�?�?�D���T�2�.�.�.�.�.�.r<c�&�tj�d��d}tj�|��}tj�|��}tj||���}tjddd��}|�|��}|�	|��}t||d���tj|��}|�|��}t||d���tj|tj�t|�������}tj|�|����}|�|��}t||d���dS)Nrrlr>r�r�r�r)rrr r!r@rr$r%r&r(rr�r��len)	r/r0rCr3r4r(r�r�r�s	         r:�test_pdf_logpdf_weightedr�qsA���I�N�N�1�����L�	�����	&�	&�B�	�����	%�	%�B���b�"�-�-�-�D�	��S�"�b�	!�	!�B�
�-�-��
�
�C��8�8�B�<�<�D���T�2�.�.�.�.�
�V�C�[�[�F��k�k�"�o�o�G�����4�4�4�4���b�"�)�.�.��R���*A�*A�B�B�B�D�
�&����r�"�"�
#�
#�C��;�;�r�?�?�D���T�2�.�.�.�.�.�.r<c��	�
��tj�d��}d}d}|�||f���}|�|df���}tjgd����	tj|���
�
��	��}|�|�	��}�	�
fd���fd�}||��}t||d	�
��dS)N�Kb1d�&rlr�rJr�)	rrr�rer�r�r��	c�v������fd�}tj|tjtj��dS)Nc�t��tj|g��f��}��|��dS�Nr)r�concatenater()rZ�
dimensionsr��points ���r:�fz<test_marginal_1_axis.<locals>.marginal_pdf_single.<locals>.f�s2�������U�:�%6�7�8�8�A��7�7�1�:�:�a�=� r<r)r�quadrr+�r�r�r�r�s` ��r:�marginal_pdf_singlez1test_marginal_1_axis.<locals>.marginal_pdf_single�sH����	!�	!�	!�	!�	!�	!�	!��~�a�"�&��"�&�1�1�!�4�4r<c�2��tj�d|���S)Nr��axis�arr�r�apply_along_axis)r�r�s �r:�marginal_pdfz*test_marginal_1_axis.<locals>.marginal_pdf�s����"�#6�Q�F�K�K�K�Kr<���ư>��rtol�
rrr��normalrQrr$�marginalr(r)r��n_data�n_dimr�r�rr(r�refr�r�r�s         @@@r:�test_marginal_1_axisr
�s�����
�)�
�
� 3�
4�
4�C�
�F��E��j�j�u�f�o�j�.�.�G�
�Z�Z�e�Q�Z�Z�
(�
(�F���5�5�5�6�6�J�
�
�W�
%�
%�C��|�|�J�'�'�H�
�,�,�v�j�)�
*�
*�C�5�5�5�5�5�5�L�L�L�L�L��,�v�
�
�C��C��4�(�(�(�(�(�(r<c��	�
�tj�d��}d}d}|�||f���}|�|df���}tjddg���	tj|���
�
��	��}|�|�	��}�	�
fd�}||��}t||d�	��dS)
Nr��rerJr�rc�>����fd�}tj|d|���S)Nc�������fd�}tj|tjtjtjtj��dS)Nc�����\}}tj||||g��}��|��dSr�)rrQr()r[rZ�w�zr�r�r�s    ���r:r�zRtest_marginal_2_axis.<locals>.marginal_pdf.<locals>.marginal_pdf_single.<locals>.f�s<����Z�(���1��H�a��A�q�\�*�*���w�w�q�z�z�!�}�$r<r)r�dblquadrr+r�s` ��r:r�zGtest_marginal_2_axis.<locals>.marginal_pdf.<locals>.marginal_pdf_single�sS����
%�
%�
%�
%�
%�
%�
%��$�Q������"�&��"�&�I�I�!�L�Lr<rr�r�)r�r�r�r�s  ��r:rz*test_marginal_2_axis.<locals>.marginal_pdf�sB���	M�	M�	M�	M�	M�	M��"�#6�Q�F�K�K�K�Kr<rrr)r�rrr�r�rr(rr	r�r�s         @@r:�test_marginal_2_axisr�s�����
�)�
�
� 3�
4�
4�C�
�F�
�E��j�j�u�f�o�j�.�.�G�
�Z�Z�e�Q�Z�Z�
(�
(�F���1�a�&�!�!�J�
�
�W�
%�
%�C��|�|�J�'�'�H�
�,�,�v�j�)�
*�
*�C�L�L�L�L�L�L��,�v�
�
�C��C��4�(�(�(�(�(�(r<c��tj�d��}d}d}|�||f���}|�|df���}t	j|��}ddg}|�|��}|�||��}ddg}	|�|	��}
|
�||	��}t||��d	}tj
t|�
��5|�ddg��ddd��n#1swxYwYd}tj
t|�
��5|�gd
���ddd��n#1swxYwYd}tj
t|�
��5|�gd���ddd��dS#1swxYwYdS)Nr�rrerJr����r���z,Elements of `dimensions` must be integers...r�g@z,All elements of `dimensions` must be unique.)rrrz;Dimensions \[-5  6\] are invalid for a distribution in 4...)rr�r�)rrr�rrr$rr(rr�rrt)
r�rrr�r�r��dimensions1�	marginal1�pdf1�dimensions2�	marginal2r��messages
             r:�test_marginal_ivr�sw��
�)�
�
� 3�
4�
4�C�
�F�
�E��j�j�u�f�o�j�.�.�G�
�Z�Z�e�Q�Z�Z�
(�
(�F�
�
�W�
%�
%�C��q�'�K����[�)�)�I��=�=���,�-�-�D��b�'�K����[�)�)�I��=�=���,�-�-�D���t����=�G�	��z��	1�	1�	1������a��X����������������������=�G�	��z��	1�	1�	1� � ����Y�Y�Y���� � � � � � � � � � � ���� � � � �N�G�	��z��	1�	1�	1�!�!����Z�Z�Z� � � �!�!�!�!�!�!�!�!�!�!�!�!����!�!�!�!�!�!s6�=D!�!D%�(D%�	E-�-E1�4E1�F:�:F>�F>c���tj�d��d}d�tj�fd�t	d|��D����}tj|��}|�tjd|����}tj	�
tj|d��d��tj	�
tj|d��d��dS)Nri�	i�c�T��g|]$}tj����|z��%Sr�)rrr!)r�rg�	n_sampless  �r:�
<listcomp>z(test_logpdf_overflow.<locals>.<listcomp>�s8�������2�9�?�?�9�-�-��3���r<rF)
rrr rQ�rangerr$r�r��testingr�isneginf�isnan)�n_dimensionsr0r3r�r s    @r:�test_logpdf_overflowr'�s�����I�N�N�1�����L��I�	������U�	�<�>�>����
�
�B���b�!�!�D�
�[�[���1�l�3�3�
4�
4�F��J���B�K��q�	�2�2�E�:�:�:��J���B�H�V�A�Y�/�/��7�7�7�7�7r<c�r�tj�d��tj�d���}tj�gd�|j���}|���}tjtj	|��|���t||dd���dS)N�90�drJ)rFg$@r*r>�+����=��atolr)rrr �	lognormal�choicerK�copyrr$�log10r)�valsr?�orig_weightss   r:�test_weights_intactr4�s����I�N�N�5����
�9���C��(�(�D��i���/�/�/�d�i��@�@�G��<�<�>�>�L�	��r�x��~�~�w�7�7�7�7��G�\��E�B�B�B�B�B�Br<c�J�tj�d��gd�}gd�}tj||���}tj|tj|�����}gd�}t
|�|��|�|��dd���dS)Nr))g�������?g+@g5@g�R@g�X@)rrrer��r>)g333333�?��Xr+r,)rrr rr$r�rr&)�valuesr?�pdf_i�pdf_fr0s     r:�test_weights_integerr<s����I�N�N�5����
*�
*�
*�F����G���v�w�7�7�7�E���v�r�z�'�/B�/B�C�C�C�E�	���B��E�N�N�2�&�&��N�N�2�&�&�U��@�@�@�@�@�@r<c�x�d�}tj�d��d}tj�|��}tj�|��}tj|��}||��tj||���}||��tjddg��}tjddgddgg��}tj�|||�	��j	}tj|��}	||	��tj||���}
||
��dS)
Nc���d}|�|��}|�|��}ttt||d���d}|�||���}|�||���}t||d���tj�d���}|�||���}tj�d���}|�||���}t||d���ttjd��r8tj�d��}|�||���dSdS)	N��g�vIh�%<=)r-i?)r �r�i�)	�resamplers�AssertionErrorrrr�RandomState�hasattrr�)�
gkde_trail�n_sample�samp1�samp2r �rstate1�rstate2r�s        r:�
test_seed_subz test_seed.<locals>.test_seed_subsd�����#�#�H�-�-���#�#�H�-�-����O�U�E��	
�	
�	
�	
����#�#�H�4�#�8�8���#�#�H�4�#�8�8����u�5�1�1�1�1��)�'�'�S�'�1�1���#�#�H�7�#�;�;���)�'�'�S�'�1�1���#�#�H�7�#�;�;����u�5�1�1�1�1��2�9�m�,�,�	4��)�'�'��-�-�C�����s��3�3�3�3�3�	4�	4r<rrr>rFrGrHrIrJ)
rrr r@r!rr$rQrRrS)rKr/rC�xn_1d�gkde_1d�gkde_1d_weightedr"rY�xn_2d�gkde_2d�gkde_2d_weighteds           r:�	test_seedrRs6��4�4�4�4�I�N�N�7�����L�	�����	%�	%�B��I�O�O�L�)�)�E�� ��'�'�G��M�'�����)�%��<�<�<���M�"�#�#�#��8�S�#�J���D���C��:��S�z�2�3�3�J��I�)�)�$�
��)�N�N�P�E�� ��'�'�G��M�'�����)�%��<�<�<���M�"�#�#�#�#�#r<c���tj�d��}tjgd���}tjgd�gd�gd�g��}|�||d��}	tj|j��dS#tj	$rTd}ttj	|���5tj|j��ddd��YdS#1swxYwYYdSwxYw)	NlSD�G.�)rr��)rer�r)r�r�r)rrr*i�z:The data appears to lie in a lower-dimensional subspace...r�)rrr�rQrRrr$rSr�LinAlgErrorrs)r��mu�sigma�datar�s     r:�%test_singular_data_covariance_gh10205rYCs>��
�)�
�
� 3�
4�
4�C�	��+�+�+�	�	�B��H�j�j�j�+�+�+�{�{�{�;�<�<�E��"�"�2�u�d�3�3�D�'�
��4�6�"�"�"�"�"����'�'�'�J��
�6�-�S�
9�
9�
9�	'�	'���t�v�&�&�&�	'�	'�	'�	'�	'�	'�	'�	'�	'�	'�	'�	'�	'����	'�	'�	'�	'�	'�	'�	'�'���s6�+B�,C)�2C�C)�C	�C)�!C	�"C)�(C)c�P�tj�d��}|�tjd��tjd��d���}d}t
jt|���5tj
|��ddd��dS#1swxYwYdS)Nl��&�8*cr�r�rJz9Number of dimensions is greater than number of samples...r�)rrr�rR�zeros�eyer�rrtrr$)r��rvsrs   r:�)test_fewer_points_than_dimensions_gh17436r^Rs���
�)�
�
� 3�
4�
4�C�

�
!�
!�"�(�1�+�+�r�v�a�y�y�q�
!�
A�
A�C�I�G�	��z��	1�	1�	1� � �
��3���� � � � � � � � � � � � ���� � � � � � s�9B�B�"B)0�scipyrrr�numpyr�
numpy.testingrrrr	r
rr�rrsr;rD�mark�slowr_rarzr~r$r�r�r�r�r�r�r��_ftypes�parametrizer�r�r�r�r
�xslowrrr'r4r<rRrYr^r�r<r:�<module>rgsM��*�*�*�*�*�*�*�*�*�*�����P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�P�
�
�
�
�*�*�*�*�*�*�E�E�E�>E�E�E�@���#J�#J���#J�L���$J�$J���$J�NO�O�O�4O�O�O�<#�#�#�#�#�U�'�#�#�#�"�"�"�"�"�U�'�"�"�"�-�-�-�-�-�U�'�-�-�-�
!4�!4�!4�H
9�
9�
9� 4�4�4�*>�>�>�J�
I�
I�������G�w��.D�$D�E�E������'�*�*�4�4�+�*�F�E�4�.���/�/�/�0/�/�/�2)�)�)�6���)�)���)�8"!�"!�"!�J���8�8���8�"C�C�C�
@�
@�
@�/$�/$�/$�d'�'�'�
 �
 �
 �
 �
 r<

Youez - 2016 - github.com/yon3zu
LinuXploit