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/__pycache__/

Upload File :
current_dir [ Writeable ] document_root [ Writeable ]

 

Command :


[ Back ]     

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

d�cvN����ddlZddlZddlmcmZddlm	Z	dZ
ejejd��e
��Z
ejejd��e
��Zejdejz��Zejdejz��ZdZejd��ZejdzZejdzZejd	zZgd
�Zd�Zd�Zdd�Zdd�Zd�Zdd�Zdd�Z dd�Z!d�Z"d�Z#dd�Z$d�Z%dd�Z&dS)�N)�_derivative���i<������)g�SˆB��g�A�Az?g}<�ٰj_�g#�+�K?g8�8�C�g��J?gl�l�f�gUUUUUU�?c��d|z}tj|��dz|z
tdzz|tjt||z��zzS)N��?r)�np�log�_LOG_2PI�polyval�_STIRLING_COEFFS)�n�rns  �6/usr/lib/python3/dist-packages/scipy/stats/_ksstats.py�_log_nfactorial_div_n_pow_nr\sG��

�Q��B�
�6�!�9�9�Q�;��?�X�a�Z�'�"�r�z�:J�B�q�D�/Q�/Q�*Q�Q�Q�c�.�tj|dd��S)z%clips a probability to range 0<=p<=1.�r)r�clip)�ps r�
_clip_probrfs��
�7�1�c�3���rTc�L�tj|||��}t|��S)z>Selects either the CDF or SF, and then clips to range 0<=p<=1.)r�wherer)�cdfprob�sfprob�cdfrs    r�_select_and_clip_probr ks!��
���g�v�&�&�A��a�=�=�rc�8�|dkrtdd|��S||z}|dkrtdd|��Sttj|����}||z
}d|zdz
}tj||g��}tjd|dz��}d||zz
}	tj|��}
d}|D]"}||
|dz
<||z}|	|dz
xx|zcc<�#td|zdz
d��|zd||zzz
}
d|
z|z|	d<td|��D]}|
d||z
dz�||dz
d�|f<�|	|dd�df<tj	|	d�	��|ddd�f<tj
tj|��d��}|}d}d}|dkr�|dzrtj||��}||z
}tj||��}|dz}tj
||dz
|dz
f��tkr|tz}|tz
}|dz}|dk��||dz
|dz
f}td|dz��D];}||z|z}tj
|��t kr|tz}|tz}�<|dkrtj||��}t|d|z
|��S)
z�Computes the Kolmogorov CDF:  Pr(D_n <= d) using the MTW approach to
    the Durbin matrix algorithm.

    Durbin (1968); Marsaglia, Tsang, Wang (2003). [1], [3].
    rr��?rrr���N)�axis)r �intr�ceil�zeros�arange�empty�max�range�flip�eye�shape�matmul�abs�_EP128�_E128�_EM128�ldexp)r�dr�nd�k�h�m�H�intm�v�w�fac�j�tt�i�Hpwr�nn�expnt�Hexpntrs                    r�
_kolmogn_DMTWrFqs���	�C�x�x�$�S�#�s�3�3�3�	
�Q��B�	�S�y�y�$�S�#�s�3�3�3��B�G�B�K�K���A�	�B��A�	�A���	�A�
��!�Q����A��9�Q��A����D��a�4�i��A�
�����A�

�C�
������!�a�%���q���	�!�a�%����C������	�Q��U�S�[�!�	�	�a�	�!�A�q�D�&�	(�B�
�2�X���A�b�E�
�1�a�[�[�%�%����!�a�%�!�)��}��!�a�%�&�&�!�)����A�a�a�a��d�G��w�q�q�!�!�!�A�b�!�!�!�e�H�
�6�"�(�1�+�+�a�.�!�!�D�	
�B�
�E�
�F�
�q�&�&�
��6�	��9�T�1�%�%�D��V�O�E��I�a��O�O���!���
�6�!�A��E�1�q�5�L�/�"�"�V�+�+�
��K�A��e�O�F�
�1�W��
�q�&�&�	
�Q��U�A��E�\��A��1�a�!�e�_�_����
��E�A�I��
�6�!�9�9�v���
��K�A��U�N�E��
��z�z��H�Q����� ��C��E�3�/�/�/rc�V�|dkr||z
dz
||zdz
}}not|dzd��\}}|dkr=||dzkr||z
|z
dz
||z|zdz
}}n3|dz
|z
|z
dz
||zdz
|zdz
}}n|dz
|z
dz
||z|zdz
}}t|dzd��t||��fS)z0Compute the endpoints of the interval for row i.rrr)�divmodr*�min)	rAr�ll�ceilf�roundf�j1�j2�ip1div2�ip1mod2s	         r�_pomeranz_compute_j1j2rQ�s����A�v�v���u��q��"�u�*�q�.�B���"�!�a�%��+�+�����a�<�<��!�a�%����R��%��!�+�Q��V�e�^�a�-?�B��� �1��r�)�F�2�Q�6��"��q�8H�5�8P�ST�8T�B����q�[�2�%��)�7�R�<�&�+@�1�+D��B��r�A�v�q�>�>�3�r�1�:�:�%�%rc���||z}ttj|����}d||z
z}t|d|z
��}|dkrdnd}|dkrdnd}d|dzz}	tj|	��}
tj|	��}tj|	��}d|
d<d|d<d|d<d}
||zd|z|zdd|zz
|z}}}td|	��D]>}|
|dz
|z|z|
|<||dz
|z|z||<||dz
|z|z||<�?tj|	g��}tj|	g��}d|d<d\}}td||||��\}}tdd|zdz��D]�}|}||}}||}}|�d��t|||||��\}}|dks|d|zdzkr|
}n	|dzr|n|}||z
dz}|dkr�tj	|||z
||z
|z�|d|���}||z
}||z
dz}||||z�|d|�<dtj
|��cxkrtkrnn|tz}|
tz}
||z|z
}��|||z
}td|dz��D]8}tj|��tkr|tz}|
tz
}
||z}�9|
dkrtj||
��}t!|d|z
|��}|S)	z[Computes Pr(D_n <= d) using the Pomeranz recursion algorithm.

    Pomeranz (1974) [2]
    rrrr"r)rrrN)r%r�floorrIr)r+r'rQ�fill�convolver*r3r1r2r0r4r ) r�xr�trJ�f�grKrL�npwrs�gpower�	twogpower�onem2gpowerrD�g_over_n�two_g_over_n�one_minus_two_g_over_nr9�V0�V1�V0s�V1srMrNrA�k1�pwrs�ln2�conv�
conv_start�conv_len�anss                                 r�_kolmogn_Pomeranzrl�s}��$	
�A��A�	�R�X�a�[�[�	�	�B��q�2�v��A��A�s�Q�w���A��a�%�%�Q�Q�Q�E��s�7�7�a�a��F�
��a��L�E�
�X�e�_�_�F������I��(�5�/�/�K��F�1�I��I�a�L��K��N�
�E�56�q�S�!�A�#�a�%�!�a��c�'�1��2�l�H�
�1�e�_�_�I�I���1�q�5�M�H�,�q�0��q�	� ��Q��'�,�6��:�	�!��$�Q��U�+�.D�D�q�H��A���	��5�'�	�	�B�	��5�'�	�	�B�
�B�q�E��H�C��
#�A�q�"�e�V�
<�
<�F�B��
�1�a�!�e�a�i�
 �
 � � ��
���R�B����S��
�������'��1�b�%��@�@���B���6�6�Q�!�a�%�!�)�^�^��D�D�!"�Q��7�I�I�K�D��2�g��k����7�7��;�r�"�s�(�2��8�c�>�"9�:�D��#��J�G�G�D��b��J��B�w��{�H� ��J��,A�!A�B�B�y��y�M��2�6�"�:�:�&�&�&�&��&�&�&�&�&��f���������(�R�-�C��
�Q��W�+�C�
�1�a�!�e�_�_����
�6�#�;�;�����6�M�C��U�N�E��q����
��z�z��h�s�E�"�"��
��S�3�Y��
4�
4�C��Jrc	��|dkrtdd|���S|dkrtdd|���Stj|��|z}|dz|dz|dz|dzf\}}}}tdz|z}|tkrtdd|���Stj|��}	|}
tdz}d|zd|zz}d|zd	|zz
tzdz}
td
d|zz
zdz}td	d|zz
zd
z}td|zd|zzzdz}td|zd|zz
zdz}d|zd|dzzz
}tjd��}ttj
d|ztjz����}t|dd��D]w}d|zd
z
}|dz|dz|dz}}}tj
|	d|z��}tjd|
||zz||
|zz||zz|||zz||zz||zzg��}||z}||z
}�x||	z}|tz}|tj|d|zd|dzzd|dzzg��z}tjtdz|z��}	tj|dd��}|dz}t"|z}tj|z}|	|z} tj|| z��}!|!ttzd|zzz}!|dxx|!z
cc<tj||z||z
z|z| z��}"|"ttzd|zzz}"|dxx|"z
cc<tj
|dztjt'|����dz��}#||#z}|s|dz}|dxxd
z
cc<t%|��}$|$S)aPComputes the Pelz-Good approximation to Prob(Dn <= x) with 0<=x<=1.

    Start with Li-Chien, Korolyuk approximation:
        Prob(Dn <= x) ~ K0(z) + K1(z)/sqrt(n) + K2(z)/n + K3(z)/n**1.5
    where z = x*sqrt(n).
    Transform each K_(z) using Jacobi theta functions into a form suitable
    for small z.
    Pelz-Good (1976). [6]
    rr�rrrrr	��r���@i�������`i���Zrr#�H�iP�
i������@)r r�sqrt�_PI_SQUARED�_MIN_LOG�exp�_PI_FOUR�_PI_SIXr'r%r&�pir+�power�array�_SQRT2PIr(�_SQRT3�sum�len)%rrVr�z�zsquared�zthree�zfour�zsix�qlog�q�k1a�k1b�k2a�k2b�k2c�k3d�k3c�k3b�k3a�K0to3�maxkr7r9�msquared�mfour�msix�qpower�coeffs�ks�ksquared�sqrt3z�kspi�qpwers�k2extra�k3extra�powers_of_n�Ksums%                                     r�_kolmogn_PelzGoodr�"s��	�C�x�x�$�S�#�3�7�7�7�7��C�x�x�$�S�#�3�7�7�7�7�
���
�
�Q��A�$%�q�D�!�Q�$��1��a��d�$:�!�H�f�e�T��<�!��h�&�D��h���$�S�#�3�7�7�7�7�
��t���A��)�C�
��/�C�
�d�(�Q��Y�
�C��u�9�q�8�|�#�{�
2�Q�
6�C�
�a�!�h�,�&�
'�"�
,�C�
�Q��h��&�
'�"�
,�C�
�c�H�n�s�U�{�2�
3�b�
8�C�
��u��r�D�y�0�
1�A�
5�C�

��*�r�A�q�D�y�
 �C��H�Q�K�K�E��r�w�r�A�v���~�&�&�'�'�D�
�4��B�
�
�	�	��
��E�A�I�� !�1��a��d�A�q�D��%����!�Q��U�#�#����3���X��-���X��-��E�	�9���X��-��E�	�9�C��H�D�F�G�G��	����
�����	�Q�J�E�	�X��E�	�R�X�q�!�e�)�R�!�Q�$�Y��q�"�u��=�
>�
>�>�E�	���|�a��(�*�+�+�A�	��4��B�	�	�B��Q�w�H�
�a�Z�F�
�5�2�:�D�
�(�]�F��f�X��&�'�'�G��{�X�%�s�V�|�4�4�G�	�!�H�H�H���H�H�H��f�f�t�m���
�6��A�F�J�K�K�G��{�X�%�s�T�z�2�2�G�	�!�H�H�H���H�H�H��(�1�s�7�B�I�c�%�j�j�$9�$9�C�$?�@�@�K�	�[��E���
����
�a����A�
�����u�:�:�D��Krc�|�tj|��r|St|��|ks|dkrtjS|dkrt	dd|���S|dkrt	dd|���S||z}|dkr�|dkrt	dd|���S|dkr:tjtjd|dz��d|zzd|zdz
z��}n?tjt|��|tj	d|zdz
��zz��}t	|d|z
|���S||dz
kr dd|z
|zz}t	d|z
||���S|dkr8dtj�||��z}t	d|z
||���S||z}|dkr�|d	kr't||d
���}t	|d|z
|���S|dkr't||d
���}t	|d|z
|���Sdtj�||��z}t	d|z
||���S|s@|dkrdS|d
kr2dtj�||��z}t|��S|dkrd}n7|dkr||dzzdkrt||d
���}nt!||d
���}t	|d|z
|���S)z�Computes the CDF(or SF) for the two-sided Kolmogorov-Smirnov statistic.

    x must be of type float, n of type integer.

    Simard & L'Ecuyer (2011) [7].
    rrrrnr"�rrg�0�q&�?Trg w@g������@g2@i��g�?gffffff�?)r�isnanr%�nanr �prodr(r�rr
�scipy�special�smirnovrFrlrr�)rrVrrW�prob�	nxsquaredrs       r�_kolmognr�us��
�x��{�{����
�1�v�v��{�{�a�1�f�f��v�
��C�x�x�$�S�#�3�7�7�7�7��C�x�x�$�S�#�3�7�7�7�7�	�A��A��C�x�x���8�8�(��c�s�;�;�;�;���8�8��7�2�9�Q��!��,�,��A��6�!�A�#��'�B�C�C�D�D��6�5�a�8�8�1�r�v�a��c�!�e�}�}�;L�L�M�M�D�$�T�3��:�3�?�?�?�?��A��E�z�z��C�!�G�a�<���$�Q��X�t��=�=�=�=��C�x�x��5�=�(�(��A�.�.�.��$�S�4�Z��3�?�?�?�?��A��I��C�x�x��� � � ��A�4�0�0�0�D�(��s�T�z�s�C�C�C�C���>�>�$�Q��t�4�4�4�D�(��s�T�z�s�C�C�C�C��5�=�(�(��A�.�.�.��$�S�4�Z��3�?�?�?�?��$������3������u�}�,�,�Q��2�2�2�D��d�#�#�#��D������	
�f����Q��V��s�*�*���1�$�/�/�/���#�A�q�d�3�3�3�� ��#��-�S�A�A�A�Arc����tj���r�St����ks�dkrtjS|dks|dkrdS�|z}|dkr�|dkrdS�dkr7tjtjd���d�zzd|zdz
z��}nBtjt����dz
tjd|zdz
��zz��}|dz�dzzS|�dz
krdd|z
�dz
zz�zS|dkr(dtj
j�|���zS|dz}t||d�zz
��}t|d|z
��}�fd	�}t|||d
���S)zvComputes the PDF for the two-sided Kolmogorov-Smirnov statistic.

    x must be of type float, n of type integer.
    rrr"rr�rrg�@c�$��t�|��S�N)�kolmogn)�_xrs �r�_kkz_kolmogn_p.<locals>._kk�s����q�"�~�~�rrp)�dx�order)rr�r%r�r�r(r�rr
r��stats�ksone�pdfrIr)rrVrW�prd�deltar�s`     r�
_kolmogn_pr��s����

�x��{�{����
�1�v�v��{�{�a�1�f�f��v�
��C�x�x�1��6�6��q�	�A��A��C�x�x���8�8��3���8�8��'�"�)�A�q�/�/�S�1�W�5��Q����C�D�D�C�C��&�4�Q�7�7�1�Q�3�"�&��Q��QR��BS�BS�:S�S�T�T�C��Q�w��A��~���A��E�z�z��C�!�G��1��%�%��)�)��C�x�x��5�;�$�(�(��A�.�.�.�.�
��K�E���q�3�q�5�y�!�!�E���s�Q�w���E�������s�A�%�q�1�1�1�1rc����tj���r�St����ks�dkrtjS�dkrd�zS|dkrdStjtj���tj��dz��z
�z��}|d�zkr|d�zzdzStj	tj|dz���z��}|dd�zz
kr|Stj���tj���z}t|dd�zz
��}��fd�}tj�|d�z|d���S)	zeComputes the PPF/ISF of kolmogn.

    n of type integer, n>= 1
    p is the CDF, q the SF, p+q=1
    rrrrr|c�*��t�|���z
Sr�)r�)rVrrs ��r�<lambda>z_kolmogni.<locals>.<lambda>�s���8�A�q�>�>�A�%�rg�+����=)�xtol)rr�r%r�r�r
r�r��loggamma�expm1�scu�	_kolmogcir}rI�optimize�brentq)rrr�r�rV�x1�_fs``     r�	_kolmognir��sM����
�x��{�{����
�1�v�v��{�{�a�1�f�f��v�
��A�v�v��1�u���A�v�v��s��F�B�F�1�I�I��
� 6� 6�q��s� ;� ;�;�Q�>�?�?�E���A��~�~���a���1�$�$�	��"�&��3��-�-��/�	"�	"�"�A��A��A��I�~�~���	��q�	�	�"�'�!�*�*�	$�B�	�R��s�1�u��	�	�B�	%�	%�	%�	%�	%�B��>� � ��S��U�B�U� �;�;�;rc�n�tj|||dgdtjtjtjg���}|D]h\}}}}tj|��r||d<�!t|��|krt
d|�����tt|��||���|d<�i|jd}|S)aComputes the CDF for the two-sided Kolmogorov-Smirnov distribution.

    The two-sided Kolmogorov-Smirnov distribution has as its CDF Pr(D_n <= x),
    for a sample of size n drawn from a distribution with CDF F(t), where
    D_n &= sup_t |F_n(t) - F(t)|, and
    F_n(t) is the Empirical Cumulative Distribution Function of the sample.

    Parameters
    ----------
    n : integer, array_like
        the number of samples
    x : float, array_like
        The K-S statistic, float between 0 and 1
    cdf : bool, optional
        whether to compute the CDF(default=true) or the SF.

    Returns
    -------
    cdf : ndarray
        CDF (or SF it cdf is False) at the specified locations.

    The return value has shape the result of numpy broadcasting n and x.
    N)�	op_dtypes.�n is not integral: rnr#)	r�nditer�float64�bool_r�r%�
ValueErrorr��operands)	rrVr�it�_nr��_cdfr��results	         rr�r��s���0
��A�q�#�t�$�"�B�J���"�*�E�
G�
G�
G�B��1�1���B��a�
�8�B�<�<�	��A�c�F���r�7�7�b�=�=��7�2�7�7�8�8�8��#�b�'�'�2�4�0�0�0��#���
�[��_�F��Mrc��tj||dg��}|D]e\}}}tj|��r||d<� t|��|krt	d|�����tt|��|��|d<�f|jd}|S)a�Computes the PDF for the two-sided Kolmogorov-Smirnov distribution.

    Parameters
    ----------
    n : integer, array_like
        the number of samples
    x : float, array_like
        The K-S statistic, float between 0 and 1

    Returns
    -------
    pdf : ndarray
        The PDF at the specified locations

    The return value has shape the result of numpy broadcasting n and x.
    N.r�r#)rr�r�r%r�r�r�)rrVr�r�r�r�r�s       r�kolmognpr�s���"
��A�q�$�<�	 �	 �B��)�)�	��B��
�8�B�<�<�	��A�c�F���r�7�7�b�=�=��7�2�7�7�8�8�8��C��G�G�R�(�(��#���
�[��_�F��Mrc�J�tj|||dg��}|D]z\}}}}tj|��r||d<�!t|��|krt	d|�����|r|d|z
fnd|z
|f\}}	tt|��||	��|d<�{|jd}
|
S)a�Computes the PPF(or ISF) for the two-sided Kolmogorov-Smirnov distribution.

    Parameters
    ----------
    n : integer, array_like
        the number of samples
    q : float, array_like
        Probabilities, float between 0 and 1
    cdf : bool, optional
        whether to compute the PPF(default=true) or the ISF.

    Returns
    -------
    ppf : ndarray
        PPF (or ISF if cdf is False) at the specified locations

    The return value has shape the result of numpy broadcasting n and x.
    N.r�rr#)rr�r�r%r�r�r�)rr�rr�r��_qr�r��_pcdf�_psfr�s           r�kolmognir�7s���&
��A�q�#�t�$�	%�	%�B��1�1���B��a�
�8�B�<�<�	��A�c�F���r�7�7�b�=�=��7�2�7�7�8�8�8�$(�8�r�1�R�4�j�j�q��t�R�j���t��3�r�7�7�E�4�0�0��#���
�[��_�F��Mr)T)'�numpyr�
scipy.specialr��scipy.special._ufuncsr��_ufuncsr��scipy._lib._finite_differencesrr2r4�
longdoubler1r3r}r�r�r
rrr�r~r�r�rrrr rFrQrlr�r�r�r�r�r�r��rr�<module>r�s��H��������#�#�#�#�#�#�#�#�#�6�6�6�6�6�6���	���-�"�-��"�"�E�	*�	*��	���-�"�-��"�"�U�F�	+�	+���2�7�1�r�u�9�����2�6�!�b�e�)������	��������e�q�j��
�5�A�:��
�%�1�*��I�I�I��
R�R�R� � � �
����H0�H0�H0�H0�V&�&�&�$Q�Q�Q�Q�hP�P�P�P�f9B�9B�9B�9B�x'2�'2�'2�T<�<�<�4"�"�"�"�J���:�����r

Youez - 2016 - github.com/yon3zu
LinuXploit