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/sparse/linalg/tests/__pycache__/

Upload File :
current_dir [ Writeable ] document_root [ Writeable ]

 

Command :


[ Back ]     

Current File : /lib/python3/dist-packages/scipy/sparse/linalg/tests/__pycache__/test_matfuncs.cpython-311.pyc
�

d�c S���dZddlZddlZddlmZmZmZmZddlm	Z	ddl
mZmZm
Z
mZmZmZddlmZmZddlmZddlmZmZmZmZmZdd	lmZdd
lmZddl m!Z!m"Z"ddlZ#ddl$Z#d�Z%d
�Z&Gd�d��Z'Gd�d��Z(dS)z3 Test functions for scipy.linalg._matfuncs module

�N)�array�eye�exp�random)�matrix_power)�assert_allclose�assert_�assert_array_almost_equal�assert_equal�assert_array_almost_equal_nulp�suppress_warnings)�
csc_matrix�SparseEfficiencyWarning)r)�expm�_expm�ProductOperator�MatrixPowerOperator�_onenorm_matrix_power_nnm)�matrix)�logm)�	factorial�binomc���|t|��ks|dkrtd���t|��}|t|��ks|dkrtd���t|��}t||��\}}tjd||z��}|tjd|��z}tj|g||z
z|��tj|g|z||z
��zS)aj
    A helper function for testing matrix functions.

    Parameters
    ----------
    n : integer greater than 1
        Order of the square matrix to be returned.
    p : non-negative integer
        Power of the matrix.

    Returns
    -------
    out : ndarray representing a square matrix
        A Forsythe matrix of order n, raised to the power p.

    �z#n must be an integer greater than 1rz p must be a non-negative integer�$@)�int�
ValueError�divmod�np�power�diag)�n�p�a�b�large�smalls      �I/usr/lib/python3/dist-packages/scipy/sparse/linalg/tests/test_matfuncs.py�_burkardt_13_powerr)s���$	�C��F�F�{�{�a�!�e�e��>�?�?�?��A���A��C��F�F�{�{�a�!�e�e��;�<�<�<��A���A��!�Q�<�<�D�A�q��H�T�A�2�a�4� � �E��B�H�T�A�2�&�&�&�E�
�7�E�7�A�a�C�=�!�$�$�r�w��w�q�y�!�A�#�'>�'>�>�>�c��tj�d��tdd��D]�}td��D]�}tj�||f��}tj�||��}t
||��}tj�|d��}t||������dS)N����)	rr�seed�range�linalgrr�normr)r"r#�M�Mp�observed�expecteds      r(�test_onenorm_matrix_power_nnmr7<s����I�N�N�4����
�1�a�[�[�0�0���q���	0�	0�A��	� � �!�Q��(�(�A���'�'��1�-�-�B�0��A�6�6�H��y�~�~�b�!�,�,�H��H�h�/�/�/�/�	0�0�0r*c���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�Zd�Zd�Zd�Zd�Zd�Z dS) �TestExpMc�v�tddgddgg��}tt|��ddgddgg��dS�N�rr-)rr
r��selfr$s  r(�test_zero_ndarrayzTestExpM.test_zero_ndarrayHsC���B�q�6�1�Q�%�.�!�!��!�$�q�'�'�A�a�5�!�A��-�8�8�8�8�8r*c��tddgddgg��}tt|�����ddgddgg��dSr;)rr
r�toarrayr=s  r(�test_zero_sparsezTestExpM.test_zero_sparseLsM����A���!�u�~�&�&��!�$�q�'�'�/�/�"3�"3�a��U�A�a�5�M�B�B�B�B�Br*c�v�tddgddgg��}tt|��ddgddgg��dSr;)rr
rr=s  r(�test_zero_matrixzTestExpM.test_zero_matrixPsC���R��F�A�a�5�>�"�"��!�$�q�'�'�A�a�5�!�A��-�8�8�8�8�8r*c�X�ttjdgg����}ttd��|��ttdgg��|��ttt	dgg����|��tttjdgg����|��tttdgg����j|��ttjdgg����}ttd��|��ttdgg��|��ttt	dgg����|��tttdgg����j|��dS)Nr-))r-��?))rF)rrrrrr�A)r>rG�Bs   r(�test_misc_typeszTestExpM.test_misc_typesTsR�����A�3�%���!�!����W�
�
�q�)�)�)���q�c�U���Q�'�'�'���V�a�S�E�]�]�+�+�Q�/�/�/���R�X��s�e�_�_�-�-�q�1�1�1���Z�!���.�.�/�/�1�1�5�5�5����B�4�&�!�!�"�"����X����*�*�*���r�d�V���a�(�(�(���V�b�T�F�^�^�,�,�a�0�0�0���Z�"���/�/�0�0�2�A�6�6�6�6�6r*c�h�tgd�gd�gd�gt���}tjd��}tjd��}t	j|d|zd|d|zz
zgd	|d
||z
zgd	d	|ggt���}t
|�����}t||��dS)N)r-�r)rr-r.)rrr��dtyper-rrK�rr.)	r�float�mathrrrrrAr
)r>rG�e1�e2r6r5s      r(�test_bidiagonal_sparsezTestExpM.test_bidiagonal_sparseas�����I�I��I�I��I�I��$�
%�
%�
%���X�a�[�[��
�X�a�[�[���8�
��2��r�2��"��9�~�&�
��A�r�B�w�K� �
��2�J�� %�&�&�&����7�7�?�?�$�$��!�(�H�5�5�5�5�5r*c��tjtjfD]_}dD]Z}|td|���z}t	|��}t||���td|���z}t
||d����[�`dS�N��{�G�z�?皙�����?g�?r-�
rKrL�d��nulp)r�float32�float64rrrr�r>rM�scalerGr5r6s      r(�test_padecases_dtype_floatz#TestExpM.test_padecases_dtype_floatos����j�"�*�-�	M�	M�E�2�
M�
M���C���/�/�/�/����7�7���u�E�2�2�2�S��%�5H�5H�5H�H��.�x���L�L�L�L�L�	
M�	M�	Mr*c��tjtjfD]_}dD]Z}|td|���z}t	|��}t||���td|���z}t
||d����[�`dSrU)r�	complex64�
complex128rrrrr_s      r(�test_padecases_dtype_complexz%TestExpM.test_padecases_dtype_complexws����l�B�M�2�	M�	M�E�2�
M�
M���C���/�/�/�/����7�7���u�E�2�2�2�S��%�5H�5H�5H�H��.�x���L�L�L�L�L�	
M�	M�	Mr*c���tj}dD]�}|tdd|d���z}t||���t	d|���z}t��5}|�td��t|d����	��}t|d	����	��}ddd��n#1swxYwYt||d
���t||d
�����dS)NrVrK�csc�rM�formatrL�=Changing the sparsity structure of a csc_matrix is expensive.T)�use_exact_onenormFrZr[)rr^�speyerrr
�filterrrrAr)r>rMr`r$�e�sup�
exact_onenorm�inexact_onenorms        r(�!test_padecases_dtype_sparse_floatz*TestExpM.test_padecases_dtype_sparse_floats]���
��.�		I�		I�E���a��%��>�>�>�>�A��E��'�'�'�#�a�u�*=�*=�*=�=�A�"�$�$�
N���
�
�2�Z�\�\�\� %�a�4� @� @� @� H� H� J� J�
�"'��U�"C�"C�"C�"K�"K�"M�"M��	
N�
N�
N�
N�
N�
N�
N�
N�
N�
N�
N����
N�
N�
N�
N�

+�=�!�#�F�F�F�F�*�?�A�C�H�H�H�H�H�		I�		Is�A"C�C
	�
C
	c�z�tj}dD]�}|tdd|d���z}t|��t	d|���z}t��5}|�td��tt|���
��|d���ddd��n#1swxYwY��dS)	NrVrKrgrhrLrjrZr[)rrdrlrrr
rmrrrrA)r>rMr`r$rnros      r(�#test_padecases_dtype_sparse_complexz,TestExpM.test_padecases_dtype_sparse_complex�s���
��.�	O�	O�E���a��%��>�>�>�>�A��E�
�
�S��%�0�0�0�0�A�"�$�$�
O���
�
�2�Z�\�\�\�.�t�A�w�w���/@�/@�!�#�N�N�N�N�
O�
O�
O�
O�
O�
O�
O�
O�
O�
O�
O����
O�
O�
O�
O��	O�	Os�A
B/�/B3	�6B3	c
���tjd��tjtjfD]�}tdd��D]�}dD]�}t
|��tj||��|zz�|��}tj	|��r|dtj||��z|zz}ttt|����|��������dS)Nr,r-rY)�-C��6?���MbP?rWrXr-rgY@rF)
rr/rr^rdr0r�rand�astype�iscomplexobjr
rr)r>rMr"r`rGs     r(�test_logm_consistencyzTestExpM.test_logm_consistency�s�����D�����j�"�-�0�	@�	@�E��1�b�\�\�
@�
@��B�@�@�E��Q���&�+�a��"3�"3�e�";�;�C�C�E�J�J�A���q�)�)�?���V�[��A�%6�%6� 6�� >�>��-�d�4��7�7�m�m�Q�?�?�?�?�@�
@�	@�	@r*c��tjgd�gd�gd�gd�g��}tt|��td|z����dS)N)���r-r-r-)r-r}r-r-)r-r-r}r-)r-r-r-r}��?)rrrr�r>�Qs  r(�test_integer_matrixzTestExpM.test_integer_matrix�s\���H��M�M��M�M��M�M��M�M�	�
�
��
	��Q����c�A�g���/�/�/�/�/r*c�D�tjgd�gd�gd�gd�gtj���}tt	|��t	d|z����t|��}tt	|��jt	d|z��j��dS)N)i���i�rr)ri���ih�)rivi����r�rrrrrLr~)rr�int16rrrrGrs  r(�test_integer_matrix_2zTestExpM.test_integer_matrix_2�s����H�'�'�'�)�)�)�'�'�'�"�l�l�$�,.�8�
5�
5�
5��	��Q����c�A�g���/�/�/��q�M�M����Q���	�4��a��=�=�?�3�3�3�3�3r*c�8�tjgd�gd�gd�gd�gt���}tjgd�gd�gd�gd	�gt���}tt	|��|d
���tjd��d
}|���}||d<t��5}|�	td��t	|��}ddd��n#1swxYwYd
}d|z}ttj||||�����dS)N)g3d����?�L�@r�r�)rgRal!�A�?r�r�)rrg���QI��?r�)rrrg�^���?rL)g��<��g�Gg�x�@g9��=��g��[d��B)r<g��'�V7�gZ�����@g��2p�)r<r<g��ѤT �g\Ge��E�@)r<r<r<gP�dž!��rv)�rtolr,g��FF�g<)r-rzIll-conditioned.*rZ)r��atol)
rrrOrrrr/�copyr
rm�RuntimeWarningr	�allclose)	r>rG�A_logm�tiny�A_logm_perturbedro�A_expm_logm_perturbedr�r�s	         r(�test_triangularity_perturbationz(TestExpM.test_triangularity_perturbation�s���

�H�&�&�&�$�$�$�!�!�!� � � �	"�
�
�
�
����
@�
@�
@�
@�
@�
@�
@�
@�
@�
@�
@�
@�A��	�	�	��	��V���a�d�3�3�3�3�	��D������!�;�;�=�=��!%����
�
 �
 �	;�C��J�J�~�':�;�;�;�$(�)9�$:�$:�!�	;�	;�	;�	;�	;�	;�	;�	;�	;�	;�	;����	;�	;�	;�	;����T�z���B�K� 5�q�t�$�O�O�O�O�P�P�P�P�Ps�++C"�"C&�)C&c��tjd��}tjd��}tjddgddggt���}tj|dgd|ggt���}t	|��}t||��dS)Nr-rrrL�rrrrOrr)r>�exp1�exp2rG�desired�actuals      r(�test_burkardt_1zTestExpM.test_burkardt_1�s���6�v�a�y�y���v�a�y�y���H�
��F�
��F���
�
�
���(�
�1�I�
��I��������a�������(�(�(�(�(r*c���tjddgddggt���}tjddgddggt���}t|��}t	||��dS)	Nr-rKrrLgk�Q�C@gb��]IG@gd��]IG@g��J�[K@�rrrOrr�r>rGr�r�s    r(�test_burkardt_2zTestExpM.test_burkardt_2s���
�H�
��F�
��F���
�
�
���(�
�!3�4�
�!3�4��������a�������(�(�(�(�(r*c��tjd��}tjd��}tjddgddggt���}tjdd|zzdd|zzz
tjd��d|zzgdtjd��zd|zzd	d|zzdd|zzzggt���}t|��}t
||��dS)
Nr-�'ri���i���rL�&i������)rrrrO�expm1rr)r>r��exp39rGr�r�s      r(�test_burkardt_3zTestExpM.test_burkardt_3s����v�a�y�y����r�
�
���H�
��F�
�#�J���
�
�
���(��B�t�G��q�"�U�(�|�+���#����"�T�'�*�
,��2�8�C�=�=� �B�t�G�,��B�t�G��r�2�e�8�}�,�
.�	�������a�������(�(�(�(�(r*c��tjddgddggt���}tjddgdd	ggt���}tjdd
gddggt���}tjd
dgt���}tj|tj|��z|��}t|��}t
||��dS)Ni����i�����rLrKr-�r�����g�?i��r�)rrrO�dotrrr)r>rG�U�V�wr�r�s       r(�test_burkardt_4zTestExpM.test_burkardt_4%s���
�H�
�"�I�
�"�I���
�
�
��
�H�q�!�f�q�!�f�%�U�3�3�3���H�q�$�i�"�c��+�5�9�9�9���H�c�2�Y�e�,�,�,���&��R�V�A�Y�Y���*�*���a�������(�(�(�(�(r*c���tjgd�gd�gd�gd�gt���}tjgd�gd�gd�gd	�gt���}t|��}t	||��dS)
N)r�rr)rrr�r)rrrr�r�rL)r-r���$)rr-r�r�)rrr-r��rrrr-r�r�s    r(�test_burkardt_5zTestExpM.test_burkardt_54s���

�H��L�L��L�L��L�L��L�L�	�
�
�
�
���(��N�N��M�M��L�L��L�L�	�
������a�������(�(�(�(�(r*c���tjd��}tjddgddggt���}tj||gd|ggt���}t	|��}t||��dS�Nr-rrLr�)r>r�rGr�r�s     r(�test_burkardt_6zTestExpM.test_burkardt_6Hs����v�a�y�y���H�
��F�
��F���
�
�
���(�
�4�L�
��I��������a�������(�(�(�(�(r*c�$�tjd��}tjd��}tjd|zdgdd|z
ggt���}tj||gd|ggt���}t|��}t
||��dSr�)rr�spacingrrOrr)r>r��epsrGr�r�s      r(�test_burkardt_7zTestExpM.test_burkardt_7Xs���
�v�a�y�y���j��m�m���H�
��W�a�L�
��C��L���
�
�
���(�
�4�L�
��I��������a�������(�(�(�(�(r*c��tjd��}tjd��}tjgd�gd�gd�gt���}tjd|z|z
d|zd|zz
d	|zd	|zz
gd
|z|zd
|zd|zzd|zd	|zzgd|zd|zd|zggt���dz}t	|��}t||��dS)
Nr��)��r�)���r�i����)r�r�r�rL�
r.ri����r�g�?r�)r>�exp4�exp16rGr�r�s      r(�test_burkardt_8zTestExpM.test_burkardt_8js���v�a�y�y����r�
�
���H��K�K��L�L��J�J���	
�
�
��
�(�
��X��_�b��h��4��/��5��1�T�6�1A�B�
��X��_�b��h��4��/��E��A�d�F�1B�C�
��X�r�%�x��5��)���	���#�	#��
�a�������(�(�(�(�(r*c���tjgd�gd�gd�gd�gt���}tjgd�gd�gd�gd	�gt���}t|��}t	||��dS)
N)r-rrr)rKr-r-r)rKrr-r)rKrKrKr-rL)�f��a�%�@�������@�X�2�1�@g鷯g)�@)�^�Iچ@�|�Pk܂@g�'���@r�)��v����@gB�f�g;�@r�r�)gD�l�{3�@r�r�r�r�r�s    r(�test_burkardt_9zTestExpM.test_burkardt_9{s���
�H��L�L��L�L��L�L��L�L�	�
�
�
�
���(�4�4�4�4�4�4�4�4�4�4�4�4�	�
������a�������(�(�(�(�(r*c�H�tjgd�gd�gd�gt���}tt	t
j�|����d��tjgd�gd�gd�gt���}t|��}t||��dS)	N)r�rr)r-r�r-)r-r-r�rL)rKrKr�)gl$�^�{b@���
|�f@g��.��Q@)g�\"N��_@r�g��
|�V@)g�\"N��_@g��L�ud@g��s��[@�	rrrOr�sorted�scipyr1�eigvalsrr�s    r(�test_burkardt_10zTestExpM.test_burkardt_10�s���
�H��I�I��I�I��I�I���	
�
�
��
	��u�|�3�3�A�6�6�7�7��C�C�C��(�E�E�E�E�E�E�E�E�E���	����
�a�������(�(�(�(�(r*c�.�tjgd�gd�gd�gt���}ttj�|��d��tjgd�gd�gd�gt���}t|��}t||��dS)	N)g��!�=@�l����?�Q����P�)r�g���t�9@猽誉\!@)r�r�g��0�n2A@rL)���()g���Jކ3Cg�;%s1P���MSUc[�)g�;%s1P�gcx
M�jC��O�}vC)r�r�g���˂C)rrrOrr�r1�eigvalshrr�s    r(�test_burkardt_11zTestExpM.test_burkardt_11�s���
�H�G�G�G�F�F�F�F�F�F���	
�
�
��
	���-�-�a�0�0�,�?�?�?��(�)�)�)�'�'�'�'�'�'�
��
�
�
���a�������(�(�(�(�(r*c�H�tjgd�gd�gd�gt���}tt	t
j�|����d��tjgd�gd�gd�gt���}t|��}t||��dS)	N)i}����r�)iz����8�6)i}����9�4rL)i��r�r�)g��)��'��g�1�3V��?g�-R��R�?)g�ǣ���g�׃4V��?g�3���?)gWN�x`��gp������?g�}m��R�?r�r�s    r(�test_burkardt_12zTestExpM.test_burkardt_12�s���
�H��N�N��N�N��N�N���	
�
�
��
	��u�|�3�3�A�6�6�7�7��G�G�G��(�H�H�H�G�G�G�G�G�G���	����
�a�������(�(�(�(�(r*c��tdd��}gd�gd�gd�gd�g}t||��dD�])}tdtt	jd|z������}t	j||ft�	��}t||z��D]�}t||��}tt	j
|��d
��tt	j|��t	jdt	j||z��|z����||t|��zz
}��tt|d����}t||����+dS)Nr�r-)rr-rr)rrr-rr�)rvrrr)rrKr�rYr�rLrrY)r)r�maxrr�ceil�zerosrOr0r�minr �floorrr)	r>�	A4_actual�
A4_desiredr"�kr�r#�Apr�s	         r(�test_burkardt_13zTestExpM.test_burkardt_13�sW��'�q�!�,�,�	�"�l�l�"�l�l�"�l�l�%�o�o�'�
�	�	�:�.�.�.��	-�	-�A��A�s�2�7�2�a�4�=�=�)�)�*�*�A��h��1�v�U�3�3�3�G��1�Q�3�Z�Z�
-�
-��'��1�-�-���R�V�B�Z�Z��+�+�+����r�
�
�B�H�R�"�(�1�Q�3�-�-���9I�,J�,J�K�K�K��2�	�!���,�,����,�Q��2�2�3�3�F��F�G�,�,�,�,�	-�	-r*c���tjgd�gd�gd�gt���}tjgd�gd�gd�gt���}t|��}t	||��dS)N)rg:�0�yE>r)g���"D��r}g _�B)g������P@rg������P�rL)g� �.��?g�zD�v>g1�O�?)g���v�U�g԰g�L��g�@��sDQ�)g��F~��?gg�=�>gc�E����?r�r�s    r(�test_burkardt_14zTestExpM.test_burkardt_14�s���
�H��L�L�(�(�(�"�"�"���	
�
�
��
�(�H�H�H�G�G�G�H�H�H���	����
�a�������(�(�(�(�(r*c	�`�dD�])}tddd��D�]}|tj|dd��z}tj|dk��rn�tjtjd|dz��d��|z}t|��}|}t
tj|dz��dd�dftj|dz��ddd�f��|ddd�fz|dd�dfz}dt|�����z}t|||�	������+dS)
N)r~rwg���ư>r�PrKr�gY��n�r-g�vIh�%<=)r�)
r0r�arange�anyr!rr�absr�r)	r>r`r"�scrGrH�gotr6r�s	         r(�test_pascalzTestExpM.test_pascal�sH��'�
	:�
	:�E��1�b�!�_�_�
:�
:���b�i��2�r�2�2�2���6�"�v�+�&�&���E��G�B�I�a��Q��/�/��4�4�u�<����G�G���� ���1�q�5�!1�!1�!�!�!�D�&�!9�!#��1�q�5�!1�!1�$�q�q�q�&�!9�;�;�=?��Q�Q�Q��Z�H�JL�Q�Q�Q�t�V�*�U���s�8�}�}�0�0�2�2�2����X�D�9�9�9�9�9��
	:�
	:r*c�p�tjd��}d|d<t|��}t��5}|�t
d��|�td��ttj|����}ddd��n#1swxYwYt||��dS)N)��r�r-)r�rzthe matrix subclass.*)	rr�rr
rm�DeprecationWarning�PendingDeprecationWarningrr)r>rG�B0rorHs     r(�test_matrix_inputzTestExpM.test_matrix_inputs����H�Z� � ����$��
�!�W�W��
�
 �
 �	#�C��J�J�)�+B�C�C�C��J�J�0�2I�J�J�J��R�Y�q�\�\�"�"�A�	#�	#�	#�	#�	#�	#�	#�	#�	#�	#�	#����	#�	#�	#�	#�	��2�����s�AB�B�"Bc
���tjgd�gd�gd�gd�gd�gd�gd�g��}t|��}td|z��}|}td��D]}||z}�t	||��dS)N)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<r<i����)rrrr0r)r>�L�E0�E1�E2�js      r(�test_exp_sinch_overflowz TestExpM.test_exp_sinch_overflows����H�;�;�;�;�;�;�;�;�;�9�9�9�9�9�9�9�9�9�9�9�9�
;�
<�
<���1�"�X�X��
�&�1�*�
�
��
���r���	�	�A��b��B�B���B�����r*N)!�__name__�
__module__�__qualname__r?rBrDrIrSrarerrrtr{r�r�r�r�r�r�r�r�r�r�r�r�r�r�r�r�r�r�rr	�r*r(r9r9Gs�������9�9�9�C�C�C�9�9�9�7�7�7�6�6�6�M�M�M�M�M�M�I�I�I�	O�	O�	O�	@�	@�	@�0�0�0�	4�	4�	4�"Q�"Q�"Q�H&)�&)�&)�P)�)�)�)�)�)�*
)�
)�
)�)�)�)�()�)�)� )�)�)�$)�)�)�")�)�)�$)�)�)�")�)�)�4)�)�)�$-�-�-�>)�)�)� :�:�:�&	�	�	� � � � � r*r9c��eZdZd�Zd�ZdS)�
TestOperatorsc��tjd��d}d}d}t|��D�]V}tj�||��}tj�||��}tj�||��}tj�||��}t|||��}	t
|	�|��|�|���|���|����t
|	j	�|��|�|���|��j	�|������XdS)Nr,r.rrY)
rr/r0r�randnrr�matmatr��T)
r>r"r��nsamples�irGrH�C�D�ops
          r(�test_product_operatorz#TestOperators.test_product_operator,s)����D����
��
�����x���	H�	H�A��	����1�%�%�A��	����1�%�%�A��	����1�%�%�A��	����1�%�%�A� ��A�q�)�)�B��B�I�I�a�L�L�!�%�%��(�(�,�,�q�/�/�*=�*=�a�*@�*@�A�A�A��B�D�K�K��N�N�Q�U�U�1�X�X�\�\�!�_�_�,?�,C�,C�A�,F�,F�G�G�G�G�	H�	Hr*c�&�tjd��d}d}d}d}t|��D]�}tj�||��}tj�||��}t||��}t
|�|��t||���	|����t
|j
�|��t||��j
�	|������dS)Nr,r.rrKrY)rr/r0rrrrrrr�r)	r>r"r�r#rrrGrHrs	         r(�test_matrix_power_operatorz(TestOperators.test_matrix_power_operator:s�����D����
��
��
�����x���	I�	I�A��	����1�%�%�A��	����1�%�%�A�$�Q��*�*�B��B�I�I�a�L�L�,�q�!�*<�*<�*@�*@��*C�*C�D�D�D��B�D�K�K��N�N�L��A�,>�,>�,@�,D�,D�Q�,G�,G�H�H�H�H�	I�	Ir*N)r
rrrrr
r*r(rr*s:������H�H�H�I�I�I�I�Ir*r))�__doc__rP�numpyrrrrr�numpy.linalgr�
numpy.testingrr	r
rrr
�scipy.sparserr�scipy.sparse._constructrl�scipy.sparse.linalg._matfuncsrrrrr�scipy.sparse._sputilsr�scipy.linalgr�
scipy.specialrrr��scipy.sparse.linalgr)r7r9rr
r*r(�<module>r's������������)�)�)�)�)�)�)�)�)�)�)�)�%�%�%�%�%�%�;�;�;�;�;�;�;�;�;�;�;�;�;�;�;�;�=�<�<�<�<�<�<�<�0�0�0�0�0�0�#�#�#�#�#�#�#�#�#�#�#�#�#�#�)�(�(�(�(�(�������*�*�*�*�*�*�*�*���������?�?�?�@0�0�0�` �` �` �` �` �` �` �` �FI�I�I�I�I�I�I�I�I�Ir*

Youez - 2016 - github.com/yon3zu
LinuXploit