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

Upload File :
current_dir [ Writeable ] document_root [ Writeable ]

 

Command :


[ Back ]     

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

d�c�y���ddlZddlZddlmZddlZddlmZddlm	Z	ddlm
Z
mZddlmZm
Z
ddlmZmZmZddlmZd	�ZGd
�d��Zd�Zej�d
���d���Zd�Zd�Zd�ZGd�d��ZdS)�N)�assert_array_almost_equal)�raises)�solve_sylvester)�solve_continuous_lyapunov�solve_discrete_lyapunov)�solve_continuous_are�solve_discrete_are)�
block_diag�solve�LinAlgError)�matrixc�`�tj�tj�tj�t
����d|��}t
j|��5}t|�	����cddd��S#1swxYwYdS)zo
    Load npz data file under data/
    Returns a copy of the data, rather than keeping the npz file open.
    �dataN)
�os�path�join�abspath�dirname�__file__�np�load�dict�items)�name�filename�fs   �A/usr/lib/python3/dist-packages/scipy/linalg/tests/test_solvers.py�
_load_datars���
�w�|�|�B�G�O�O�B�G�O�O�H�,E�,E�F�F�"�D�*�*�H�	���	�	��a��A�G�G�I�I�����������������������s�5!B#�#B'�*B'c��eZdZejddgddgg��ejddgddgg��fejd	d
gddgg��ejd
dgdd
gg��fejdd
gddgg��ejd
dgdd
gg��fejd	d
gddgg��ejd
d
gdd
gg��fejgd�gd�gd�gd�gd�g��ejgd�gd�gd�gd�gd�g��fejgd�gd�gd�gd �gd!�gd"�gd#�gd$�gd%�gd&�gd'�g��ejd��fed(dgd)d*gg��ed(dg��jed(dg��jjzfed(dgd)d*gg��ejed(dg��jed(dg��jjz��fgZd+�Z	d,�Z
d0d.�Zd/�Zd-S)1�TestSolveLyapunov�����	�
����?�?�@�@��@�@��@@�����?�@��)r#r%�r!r$)r!r"r#�r$)r$�r5r5r#)r!r3r"r�)r3r#r#r!r3)r"r$r!rr!)r$r!rr"r)r!rr#rr#)rr"rr!r)r!rr#rr$)y�������?y����K�?yˡE��?y㥛� ��?y����Mb�?y)\��(�?yZd;�O��?y�l��?yy�&1��?y�~j�t�?y{�G�z�?)��?�r8r8r8r8r8r8r8r8r8)r8r7r8r8r8r8r8r8r8r8r8)r8r8r7r8r8r8r8r8r8r8r8)r8r8r8r7r8r8r8r8r8r8r8)r8r8r8r8r7r8r8r8r8r8r8)r8r8r8r8r8r7r8r8r8r8r8)r8r8r8r8r8r8r7r8r8r8r8)r8r8r8r8r8r8r8r7r8r8r8)r8r8r8r8r8r8r8r8r7r8r8)r8r8r8r8r8r8r8r8r8r7r8rg����c�"�tjd��}tjd��}ttt
||��ttt
||��ttt
|tjd����dS)N�r#r"r#r")r�ones�eye�
assert_raises�
ValueErrorr)�self�nsq�sqs   r�$test_continuous_squareness_and_shapez6TestSolveLyapunov.test_continuous_squareness_and_shapeXsi���g�f�o�o��
�V�A�Y�Y���j�";�S�"�E�E�E��j�";�R��E�E�E��j�";�R������K�K�K�K�K�c���t||��}ttj||��tj||��������z|��dS�N)rrr�dot�conj�	transpose)r@�a�q�xs    r�check_continuous_casez'TestSolveLyapunov.check_continuous_case_sf��%�a��+�+��!��&��A�,�,����1�6�6�8�8�3E�3E�3G�3G�)H�)H�H�!�	M�	M�	M�	M�	MrDNc���t|||���}ttjtj||��|��������|z
d|z��dS)N��methodr2)rrrrGrHrI)r@rJrKrPrLs     r�check_discrete_casez%TestSolveLyapunov.check_discrete_casedso��#�A�q��8�8�8��!��f�R�V�A�q�\�\�1�6�6�8�8�+=�+=�+?�+?�@�@�1�D�d�1�f�	N�	N�	N�	N�	NrDc�2�|jD]�}|�|d|d��|�|d|d��|�|d|dd���|�|d|dd�����dS)Nrr!�directrO�bilinear)�casesrMrQ�r@�cases  r�
test_caseszTestSolveLyapunov.test_casesis����J�	J�	J�D��&�&�t�A�w��Q��8�8�8��$�$�T�!�W�d�1�g�6�6�6��$�$�T�!�W�d�1�g�h�$�G�G�G��$�$�T�!�W�d�1�g�j�$�I�I�I�I�		J�	JrDrF)
�__name__�
__module__�__qualname__r�arrayr=r
�TrUrCrMrQrX�rDrr r sP������
���A�q�6�A�q�6�"�	#�	#�	���A�r�7�R��H�%�	&�	&�	(�
���F�C�=�8�S�/�2�	3�	3�	���F�F�#�g�s�^�4�	5�	5�	7�
���C��:��S�z�*�	+�	+�	���F�F�#�g�s�^�4�	5�	5�	7�
���F�C�=�8�S�/�2�	3�	3�	���C��:��c�{�+�	,�	,�	.�
���?�?�?�O�O�O�_�_�_�"�?�?�O�O�O�5�
6�
6�	���?�?�?�O�O�O�_�_�_�"�?�?�O�O�O�5�
6�
6�	7�

������������������������������������=  � 
!� 
!�B
������C!	�F
��!�Q��$���$�	%�	%�
�&�!�Q��.�.�
�V�V�Q��F�^�^�-�/�
/�	2�
��!�Q��$���$�	%�	%�
�"�(�6�6�1�a�&�>�>�#�f�f�a��V�n�n�&6�&8�8�
9�
9�	<�q:
�E�xL�L�L�M�M�M�
N�N�N�N�
J�J�J�J�JrDr c�~�td��}td��}td��}td��}td��}tjdgd��tjdgdgg��t	dd	��ddftjd
dgdd
gg��tjdgdgg��tjddgddgg��ddftjgd�gd�gd�gd�g��tjddgddgddgddgg��tjgd�gd�gd�gd�g��tjd��dftjgd �gd!�gd"�gd#�gd$�gd%�gd&�gd'�g��tjgd(�gd)�g��jd*ztjgd+�gd,�gd-�gd.�gd/�gd0�gd1�gd2�g��tjd��dftjgd3�gd4�gd5�gd6�gd7�gd8�gd9�gd:�gd;�g	��tjgd<�gd=�gd>�gd?�gd@�gdA�gdA�gdA�gdA�g	��tjd��tjd��df|dB|dC|dD|dEdftjddgddFgg��tjdGgdgg��tjdH��ddIft	dJdK��tjdLdMgd*dNgg��tjdOdPgdPdQgg��tjdH��t	dGd��zdftjddRgddgg��tjdgdgg��tjd��ddftjdSdgddSgg��tjd��tjd��tjd��dftjddgd
dgg��tjdgdgg��tjdTdUgdUdFgg��ddftjgdV�gdW�gdX�g��dztjd��tjgdY�gdZ�gd[�g���tjgd\������tjgdY�gdZ�gd[�g����dztjd��dRzd]ftjgd^�gd_�gd`�gda�g��tjgdb�g��jtjgdc���ddftjgdd�gde�gdf�gdg�g��tjdh��tjdi��ddf|dB|dC|dD|dEdftjdjdjd�k��tjdjdj��dFzztj	t	dtj
dl��d����ztjdjdjd�k��ztjdj��tjdj��tjdj��dftjtjdm��d��tjtjdnd����tjdnd��tjdnd��jzdd]f|dB|dC|dD|dEdf|dB|dC|dD|dEd]f|dB|dC|dD|dEd]fg}do}dp�}t|��D]\}}	||	||���dS)qNzcarex_6_data.npzzcarex_15_data.npzzcarex_18_data.npzzcarex_19_data.npzzcarex_20_data.npzr0r!rr*r$r#����r9r%r5�@�rr!rr)rg=
ףp=��g�(\���?g��Q��)rg�� �rh��gףp=
��gq=
ףp@)g�� �rh�?g/n��R�g�G�z��g�z�G�ʿg
ףp=
�?g��������gffffff�g����Mb��g���Q��?)g/�$�@�7�A`��@�j�t��?�Zd;�O��?)rdg/�$�@��|?5^�?�y�&1��?)rergg�x�&1�?��&1��?)rfrhrigsh��|?�?r")g��ʡE��g!�rh���?rrrrrr)�NbX9��?g7�A`�����O��n�?rrrrr)rrjg��S��rkrrrr)rrrjg�S㥛��g`��"��?rrr)rrr�/�$��?g=
ףp=��g� �rh��?rr)rrrrrlg��K7���g�E����?r)rrrrrrlg^�I+��g�t�V�?)rrrrrrrlgJ+���)g���Q�@rbg������B@�p=
ף@g�z�G�@g
ףp=
@rmr1)g
ףp=
�gR���Q�gffffff�g��(\��g��(\�
�g��(\��g{�G�z�g�G�z�����MbP?)r0�roro��?roro皙�����?)ror0rororqrororo)roror0rororproro)rororor0rorororo)rprqrororqrororo)rororprororqroro)rororororororqro)rqrororororororq)	g�~j�t�g{�G�z@rorog-����rororogףp=
��?)	g��/�$ֿgZd;�O��?rorog���Q��rororog
ףp=
�?)	g�A`��g���Mb�.@g-���F�rog�������rororogq=
ףp@)	g7�A`��5�g!�rh��A@gV-�տg�(\���g������1�rororog��(\�@)	g+�N�gy�&1�X@gT㥛Ġ�g�(\���?g/�$�J�rororogffffff4@)	rrrr��W@�ffffffb�rog������J@ro)	rrrrrrrrsrr)	rrrrr皙����)@ro皙����?�r)	rrrrrtrorog������2@ru)�{�G�z�?g�I+���g���Sÿ)g�~j�t�h?g/�$���ro)g;�O��n�?gh��|?5��ro)g�~j�t��?g#�~j�Ŀro)g�� �rh�?g{�G�zܿro)rororo�A�B�Q�R�����ư>�r"r"zBad residual accuracyg��������g{�G�z��rqrorv�di�i'���.Ag����?i�������)g�ZA���>Ag�)r�g`�VA���>�)ror�g�SA)r0r{r{)r{r0r{)r{r{r0)r|r!rzBad Residual Accuracy)rg�������?rr)rrg�G�z�?r)rg���g�a�g��A)rrrg��.�)rrrr)r!rr!r)���ư�r!rr)r9r�rr)rrr|r!)rrr9r|)r$r!)r$r$�@��k)�>r�)��)�r(�
r�r'r5Nr3r6r�r�Nr%r�r�r�Nr(NNc��|\}}}}}|rtj|���t||||��}|�|��|���j�|��z|z}|�|��}	||	�t
tj|��|	���j����z}t|tj
|��|���dS)z3Checks if 0 = XA + A'X - XB(R)^{-1} B'X + Q is true��reason��decimalN)�pytest�xfailrrGrHr]rr�
atleast_2dr�
zeros_like)
rW�decrJ�brK�r�knownfailurerL�res�out_facts
          r�
_test_factoryz0test_solve_continuous_are.<locals>._test_factory&s���#'� ��1�a��L��	.��L��-�-�-�-� ��A�q�!�,�,���e�e�A�h�h����������*�*�*�Q�.���5�5��8�8���x�|�|�E�"�-��"2�"2�H�M�M�O�O�4E�F�F�G�G�G��!�#�r�}�S�'9�'9�3�G�G�G�G�G�GrD)
rr�diagr\r
r=r]r<rG�rot90�zeros�flipud�	enumerate)
�mat6�mat15�mat18�mat19�mat20rU�min_decimalr��indrWs
          r�test_solve_continuous_arer�qs����(�)�)�D��*�+�+�E��*�+�+�E��*�+�+�E��*�+�+�E�
��"��q�	�	�	��A�3���*�	�	�	�B��	�	�	
�	
�		�
��A�q�6�D�$�<�(�	)�	)�	��A�3���+�	�	�	��A�q�6�A�r�7�#�	$�	$�	
�	
�		�
��<�<�<�*�*�*�+�+�+�1�1�1�3�
4�
4�
��A�q�6�D�$�<�%���4��)�D�	E�	E�	��/�/�/�/�/�/�/�/�/�/�/�/�1�
2�
2�
�����	
�
	�
��4�4�4�8�8�8�8�8�8�8�8�8�7�7�7�8�8�8�8�8�8�4�4�4�6�
7�
7�
��D�D�D�K�K�K�M�
�
���
�
��;�;�;�;�;�;�;�;�;�;�;�;�;�;�;�;�;�;�;�;�;�;�;�;�=�
>�
>�
�����	
�)	�,
��=�=�=�=�=�=�B�B�B�H�H�H�H�H�H�9�9�9�2�2�2�5�5�5�>�>�>�@�	
A�	
A�
��*�*�*�)�)�)�)�)�)�)�)�)�)�)�)�(�(�(�(�(�(�(�(�(�(�(�(�*�
+�
+�
�����	�����	
�+	�.
�c��D��I�t�C�y�$�s�)�T�:�	��A�q�6�A�s�8�$�	%�	%�	��D�6�A�3�-�	 �	 �	�����	�	 �		"�
�D�%�	 �	 �	��E�5�>�E�5�>�2�	3�	3�	��C��;��u�
�.�	/�	/�	�����:�d�A�.�.�	.�	
�		�
��A�s�8�a��V�$�	%�	%�	��A�3���+�	�	�	�����	�	
�		�
��I�q�>�B�	�?�3�	4�	4�	�����	�����	�����	
�		�
��A�r�7�Q��F�#�	$�	$�	��A�3���*�	�	�	��C��9�r�3�i�(�	)�	)�	�	
�		�
��,�,�,�2�2�2�,�,�,�.�
/�
/�12�
3�
�����	��>�>�>�>�>�>�>�>�>�B�	C�	C�	G�	G������'�'�
)�
)�),���H�n�n�n�n�n�n�n�n�n�E�F�F�*H�*H�JK�
L�
�����S��	 �	"�
��>�>�>�#�#�#�3�3�3�"�?�?�$�
%�
%�
��>�>�>�"�	#�	#�	%�	�����	�	�	�	
�	�
��#�#�#�$�$�$�"�?�?�#�#�#�%�
&�
&�
�����	�����	�	
�	�
�s��U�3�Z��s��U�3�Z��>�	���B�"�	�	�	���r�2���� 4�	4�r�x��A�r�x��1�1�1�5�5�87�87�
7�9;���B�!�9L�9L�9L�
M�	�����	�����	�����	
�	�
�������!�	$�	$�	��2�6�"�a�=�=�	!�	!�	���A������A����	(�	
�	 �		"�
�s��U�3�Z��s��U�3�Z��>�	�s��U�3�Z��s��U�3�Z�	 �	"�
�s��U�3�Z��s��U�3�Z�	 �	"�Cc
�E�X>�K�
H�
H�
H��u�%�%�.�.�	��T��
�d�K��,�-�-�-�-�.�.rDz;test fails, see https://github.com/scipy/scipy/issues/14517r�c��tjddgddgg��tjdgdgg��tjddgddgg��tjdgg��dftjddgddgg��tjdgdgg��tjddgddgg��tjdgg��dftjd	dgddgg��tjddgdd	gg��tjdd
gddgg��tjddgdd	gg��dftjgd
�gd�gd�g��tjdgdgdgg��tjd	��tjd��dftjdd	gddgg��tjdgdgg��tjddgddgg��tjdgg��dftjddgddgg��tjddgddgg��tjddgddgg��tjddgdd	gg��dftjddgddgg��tjdgdgg��tjddgddgg��tjdgg��dftjddgddgg��tjddgddgg��tjd d gd d!gg��d"ztjdd	gd	dgg��dftjddgddgg��tjdgdgg��tjddgddgg��tjdgg��dftjgd#�gd$�gd%�gd&�g��tjd'dgd(d)gd*d+gd,d(gg��tjgd-�gd.�gd/�gd0�g��tjd��dftjgd1�gd2�gd3�gd4�g��tjd5d6gd7d8gd9d:gd;d<gg��tjd��d=ztjd��dftjgd>�gd?�gd@�gdA�g��tjgdB�gdC�gdD�gdE�g��tjgdF�gdG�gdH�gdI�g��tjd��dftjgdJ�gdK�gdL�gdM�gdN�g��d=ztjdOdPgdQdRgdSdTgdUdVgdWdXgg��d=ztjdY��tjd��dftjtjd��tjddgd�Z����tjtjd��tjdgdgdgg����tjgd[�gd[�gd\�gd]�gd^�gd\�g��tjd	dgddgg��dfd_tjgd`�gda�gdb�gdc�gdd�gde�gdf�gdg�gdh�g	��ztjgdi�gdj�gdk�gdl�gdm�gdn�gdo�gdp�gdq�g	��d_ztjgdr���tjd	��dftjddsgddgg��tjdgdgg��tjd��tjdgg��dtftjgdu�gdv�gdw�g��dxztjd	��dstjd	��zdstjd	��zdyftjgdz�gd{�gd|�gd}�g��tjd~gdgdgdgg��tjgd���tjd�gg��dftjd�d�Z��tjtjd�d����tjd���tjdgg��dfg}d�}d��}t
|��D]\}}||||���tjtjd�����}d|d�<tjgd��gd��g��j	}tj
|d���tjgd����z}tjd�d(g��}ttt||||��dS)�Nr"y�?�ry��r!r9y����?r#r)y�?�@)g�X=X�I�?gPŽ���?g�'9ԭf�?)g��ҍK�?g��h�2�?g�N�ցZ�?)gS����?g�s �,�?g^L��^��?g�=S���?g�>%�S�?g���FB�?r$r`rar%r5g�~�:p�?g2w-!��?g5^�I�@g�x�&1�gZd;�O�@g{�G�zt?g{�G�z�?gUUUUUU�?���r6gF]t�E�?)�V-����?g��x�&�?rr)g��x�&��r�rr)rrr�g/�$���?)rrg/�$��ÿr�gF%u�k?rqgǺ���F�g{�G�z�?g�<,Ԛ�}?gǺ���f�)g�Q���?rr��Zd;Ͽ)rg+����?�=
ףp=�?r)rr�g+���?r)r�rrg�S㥛��?)gP��n��?g�/��t��g����M?g[A�+�Q�)g��ϷK�?gH�`��?g$��(�[��g��u���?)gY��9�}�g�_"�:��?g!�'�>տgC�K��?)g��_�L@ge�X���g������g+��пg� T�fn?g`[���H?gQ�n�H��g	v���>g�}�k�,�g�UW�;�g{�/L��	�gO�a�eRN?rv)g333333�g�������g�������g�>�Q���)r0g333333�?g�������?g��G�)3@)ror0g������?g{�\f@)rororogg�R@��)r0r2r2r2)ror0r2r2)roror0r2)rororor0)r"r!r#r5)r!r"r"r3)r#r"r5r')r5r3r'�)g5^�I�W@gY�8��m�?g��&S�?g&S��:�?g�~j�t�?)g�&1�lD@gj�t��D@g/�$�0@g�1w-!�@ge�`TR'�?)g�v��o(@g�~j�tS:@gP��nB@g\��(�/@gj�t��(@)gm��{r@g���K�)@g/�$�5;@g��x�&q5@g�rh��|D@)gNbX9��?g������?g�?g�\�C��@gR���W@g��JY�8�?g�HP���g����H@g$(~���g�/�$@g
ףp=
�g�e�c]�@gffffff�gx$(~��?g�/L�
F�r3r�)r!r!rrrr)rrrrrr)rrrr!r9r)rrrr9r!rrn)	g�����0�@g�`@g�G�z.'@g]�Fx�?g\��(�B�g#J{�/L�?rg+���@g��"��@)	g33333#S@g33333�@gq=
ףp)@g"��u��?g�G�zD�g�#�����?rg��Q�@g��"��@)	g������_�gXv@i1gH�z�G�?g33333�Y�g/�$���?rg33333�'@g���Q�2@)	g�v�g33333σ@g
ףp=�R@g������@gq�g�MbX9@rg���Q�?@g)\��hH@)	i@���g�������@g�����`�g��Sc�g{�G��Q@g;�O��n@rg�z�G!U@g�����y_@)	g33333Ä�g=
ףp=\@g)\��8V�g�E����g�G�zU@g333333+@rg�����	b@gffffffY@)	g33333�y�i�g{�G�ZK�g+����g��(\�P@g{�G�(@g���z6�?g�G�z�X@g{�G�jQ@)	g�����|f�g33333�r@g�G�z�7�g��(\��g�Q��KN@g)\��(6@rg����̼j@g��Q��A@)	g������u�g33333#�@g{�G�F�gm����gffffffZ@g\��(�3@rg33333ck@gfffff�j@)g
ףp=
@g c�ZB>�g�Q���T�)gT㥛� �?gd;�O��gףp=
��)g��ʡE��?g=
ףp=*�g��n��!@)g�A`��"@gq=
ף�A�g������8@)g�G�z�$@g0W�g33333�P@)gZd;��@g�Q���N�g�Q��+C@)g7�A`��@g
ףp=jB�g
ףp=J4@)g��K7��?g�G�z/�gsh��|�@)g�v��/@gH�z�G>�g�z�Ga-@)	�2rrrr�rrrrrz*Presumed issue with OpenBLAS, see gh-16926)�r&���)r&r����)r�r�r6g�q�q�?z!Issue with OpenBLAS, see gh-16926)gG������?rrr)r!rrrrc)rrr!rg:�0�yE>)rrrr!g�?r~)r(r�r�r�r�r��r�r�r�r�r�r�r�r(r"r3r5r&c��|\}}}}}|rtj|���t||||��}|���j�|�|����|z
|z}||���j�|�|�����t
||���j�|�|����z|���j���|�|������z}t|tj	|��|���dS��5Checks if X = A'XA-(A'XB)(R+B'XB)^-1(B'XA)+Q) is truer�r�N)
r�r�r	rHr]rGrrrr�)	rWr�rJr�rKr�r�rLr�s	         rr�z.test_solve_discrete_are.<locals>._test_factorys&��#'� ��1�a��L��	.��L��-�-�-�-��q�!�Q��*�*���f�f�h�h�j�n�n�Q�U�U�1�X�X�&�&��*�Q�.���q�v�v�x�x�z�~�~�a�e�e�A�h�h�'�'�+�+��!�A�F�F�H�H�J�N�N�1�5�5��8�8�4�4�4�a�f�f�h�h�j�A�A�E�E�a�e�e�A�h�h�O�O���	��	"�#�r�}�S�'9�'9�3�G�G�G�G�G�GrD�r#r#)rr!)r!r!r)rrr!r�)r4r9gffffff��i����)rr\r=�kronr�r�r��triur<r]�	full_liker>rr)	rUr�r�r�rWrwrxryrzs	         r�test_solve_discrete_arer�5s�
��$
��A�t�9�q�#�h�'�	(�	(�	��A�3���*�	�	�	��A�q�6�A�q�6�"�	#�	#�	��A�3�%���	
�		�
��A�q�6�A�r�7�#�	$�	$�	��C�5�2�$�-�	 �	 �	��A�q�6�A�q�6�"�	#�	#�	��A�3�%���	
�		�
��A�q�6�A�r�7�#�	$�	$�	��A�q�6�A�q�6�"�	#�	#�	��A�t�9�t�Q�i�(�	)�	)�	��A�s�8�b�!�W�%�	&�	&�	
�		�
��7�7�7�6�6�6�6�6�6�8�
9�
9�
��J�<�*��
�|�<�	=�	=�	�����	�����	
�
	�
��A�q�6�D�$�<�(�	)�	)�	��A�3���+�	�	�	��A�q�6�A�q�6�"�	#�	#�	��A�3�%���	
�		�
��F�A�;��F��,�	-�	-�	��E�5�>�G�U�#3�4�	5�	5�	��E�1�:��4�y�)�	*�	*�	��C��8�a��V�$�	%�	%�	
�		�
��A�r�7�Q��F�#�	$�	$�	��A�3���*�	�	�	��A�q�6�A�q�6�"�	#�	#�	��A�3�%���	
�		�
��A�q�6�A�r�7�#�	$�	$�	��A�q�6�A�q�6�"�	#�	#�	��B��8�b�!�W�%�	&�	&�$�	/�	��A�q�6�A�q�6�"�	#�	#�	
�		�
��A�q�6�A�q�6�"�	#�	#�	��A�3���*�	�	�	��A�q�6�A�q�6�"�	#�	#�	��A�3�%���	
�		�
��'�'�'�'�'�'�'�'�'�'�'�'�)�
*�
*�
��F�F�#��F�#��F�#��F�#�%�
&�
&�
��'�'�'�'�'�'�'�'�'�(�(�(�*�
+�
+�
�����	
�
	�
��=�=�=�=�=�=�=�=�=�=�=�=�?�
@�
@�
��I�y�)��y�)��z�*��y�)�+�
,�
,�
�����4��	�����	
�
	�
��C�C�C�?�?�?�?�?�?�@�@�@�B�
C�
C�
��*�*�*�)�)�)�(�(�(�'�'�'�)�
*�
*�
��<�<�<��<�<� �=�=�!�>�>�#�
$�
$�
�����	
�
	�
��<�<�<�>�>�>�@�@�@�?�?�?�<�<�<�	>�
?�
?�BF�	
F�

��F�G�$��G�$��G�$��G�$��G�$�	&�
'�
'�*.�	
.�

�����	�����	
�	�
�������B�G�Q��F�a�0�0�0�	1�	1�	�������B�H�q�c�A�3���_�5�5�	6�	6�	��%�%�%�%�%�%�%�%�%�&�&�&�&�&�&�%�%�%�'�
(�
(�
��A�q�6�A�q�6�"�	#�	#�	
�		�
���
F�
F�
F�
F�
F�
F�
D�
D�
D�
G�
G�
G�
G�
G�
G�
H�
H�
H�
J�
J�
J�
H�
H�
H�
H�
H�
H�
J�	K�	K�	
K�
��.�.�.�-�-�-�-�-�-�.�.�.�/�/�/�.�.�.�.�.�.�-�-�-�.�.�.�0�
1�
1�49�
9�
��.�.�.�	/�	/�	�����	
�+	�.
��A�s�8�a��V�$�	%�	%�	��A�3���*�	�	�	�����	��A�3�%���	5�		7�
��<�<�<��,�,��+�+��
 �
 �#&�
(�
�����	�r�v�a�y�y��	�r�v�a�y�y��	,�
	.�
��'�'�'��,�,��,�,��,�,� �
!�
!�
��E�7�Q�C�!��q�c�*�	+�	+�	�����	�	�	��D�6�(�	�	�	
�	�
���q�	�	�	�	��2�6�#�q�>�>�	"�	"�	�����	��A�3�%���	
�		�q}
�E�N4�K�H�H�H��u�%�%�.�.�	��T��
�d�K��,�-�-�-�-�	������� � �A��A�d�G�
��)�)�)�Y�Y�Y�'�(�(�*�A�
��Q����b�g�m�m�m�4�4�4�A�
���c�
���A��+�3�Q��1�a�@�@�@�@�@rDc��tjgd�gd�gd�g��tjddgddgdd	gg��tjd
��tjd��tjgd�gd
�gd�g��tjd��dftjgd�gd�gd�g��tjddgddgdd	gg��tjd
��tjd��tjgd�gd
�gd�g��tjd��dfg}d}d�}t|��D]\}}||||���dS)N�g?� ?��?g[,�Y�?g��7h�?�g;)�ƣ�?g7j*0	<�?gɚ&�ꢡ?�g���ڠݸ?g*���[K�?g�fSdc�?�6��Wtk�?��]�p��?�;��?��/�LX�?�	�Gr�?�m!#�{��?r#r"�gGY�����?g X.e��?g6��u��?�g�����?g�p���?g����v�?�g�<�d&�?g������?g���2��?r;)r&r&c���|\}}}}}}}|rtj|���t||||||��}	|���j�|	�|����|���j�|	�|����z|z}
|���j�|	���|��|z}|
|�t
tj|��|���j����z}
t|
tj
|
��|���dSr�)r�r�rrHr]rGrrr�rr�)rWr�rJr�rKr��e�sr�rLr�r�s            rr�z<test_solve_generalized_continuous_are.<locals>._test_factory=s��)-�&��1�a��A�q�,��	.��L��-�-�-�-� ��A�q�!�Q��2�2���f�f�h�h�j�n�n�Q�U�U�1�X�X�&�&������������a���)A�)A�A�A�E���6�6�8�8�:�>�>�!�$�$�(�(��+�+�a�/���x�|�|�E�"�-��"2�"2�H�M�M�O�O�4E�F�F�G�G�G��!�#�r�}�S�'9�'9�3�G�G�G�G�G�GrD)rr\r=r�r<r�)rUr�r�r�rWs     r�%test_solve_generalized_continuous_arer�s���
��=�=�=�=�=�=�=�=�=�?�
@�
@�
��L�,�/� �,�/� �,�/�1�
2�
2�
�����	�����	��=�=�=�=�=�=�=�=�=�?�
@�
@�
��&�	�	�	
�	�
��=�=�=�=�=�=�=�=�=�?�
@�
@�
��L�,�/� �,�/� �,�/�1�
2�
2�
�����	�����	��=�=�=�=�=�=�=�=�=�?�
@�
@�
�����	
�	�!
�E�>�K�
H�
H�
H��u�%�%�.�.�	��T��
�d�K��,�-�-�-�-�.�.rDc��td��}tjgd�gd�gd�g��tjddgddgd	d
gg��tjd��tjd��tjgd
�gd�gd�g��tjd��dftjgd�gd�gd�g��tjddgddgd	d
gg��tjd��tjd��tjgd
�gd�gd�g��tjd��df|d|d|d|dd|ddfg}d}d�}t
|��D]\}}||||���dS)Nzgendare_20170120_data.npzr�r�r�r�r�r�r�r�r�r#r"r�r�r�r;rwrxryrz�S)r'r'r�c�$�|\}}}}}}}|rtj|���t||||||��}	|�tj|jd��}|�tj|��}|���j�	|	�	|����|���j�	|	�	|����z
|z}
|
|���j�	|	�	|����|z�	t||���j�	|	�	|����z|���j�	|	�	|����|���jz����z}
t|
tj|
��|���dS)r�r�Nrr�)r�r�r	rr=�shaper�rHr]rGrr)rWr�rJr�rKr�r�r�r�rLr�s           rr�z:test_solve_generalized_discrete_are.<locals>._test_factoryzs���)-�&��1�a��A�q�,��	.��L��-�-�-�-��q�!�Q��1�a�0�0���9���q�w�q�z�"�"�A��9��
�a� � �A��f�f�h�h�j�n�n�Q�U�U�1�X�X�&�&������������a���)A�)A�A�A�E��������
���q�u�u�Q�x�x�(�(�1�,�1�1��!�A�F�F�H�H�J�N�N�1�5�5��8�8�4�4�4��6�6�8�8�:�>�>�!�%�%��(�(�3�3�a�f�f�h�h�j�@�����	��
	"�#�r�}�S�'9�'9�3�G�G�G�G�G�GrD)rrr\r=r�r<r�)�mat20170120rUr�r�r�rWs      r�#test_solve_generalized_discrete_arer�Ms%���8�9�9�K�

��=�=�=�=�=�=�=�=�=�?�
@�
@�
��L�,�/� �,�/� �,�/�1�
2�
2�
�����	�����	��=�=�=�=�=�=�=�=�=�?�
@�
@�
��&�	�	�	
�	�
��=�=�=�=�=�=�=�=�=�?�
@�
@�
��L�,�/� �,�/� �,�/�1�
2�
2�
�����	�����	��=�=�=�=�=�=�=�=�=�?�
@�
@�
�����	
�	�
�S�	�	�S�	�	�S�	�	�S�	�	
�	�S�	�	
�
	�?&
�E�P�K�H�H�H�&�u�%�%�.�.�	��T��
�d�K��,�-�-�-�-�.�.rDc�$�d�}d�}d�}d�}d�}dS)Nc
�B�tjd��}tjd��}ttfD]g}tt||ddd��tt||||d��tt|||||��tt||||||���hdS)Nr;r#r!�rr<r=rr	r>r?�rArBrLs   r�test_square_shapez1test_are_validate_args.<locals>.test_square_shape�s����g�f�o�o��
�V�A�Y�Y��&�(:�;�	>�	>�A��*�a��a��A�6�6�6��*�a��R��a�8�8�8��*�a��R��S�9�9�9��*�a��R��R��=�=�=�=�		>�	>rDc���tjd��}tjd��}ttfD]�}tt|||dd��tt|||||||��tt|||tjd��|��tt||||tjd����tt|||||tjd������dS)Nr;r$r!r#r�r�s   r�test_compatible_sizesz5test_are_validate_args.<locals>.test_compatible_sizes�s����g�f�o�o��
�V�A�Y�Y��&�(:�;�	D�	D�A��*�a��S�!�Q�7�7�7��*�a��R��R��S�A�A�A��*�a��R������B�?�?�?��*�a��R��R�V�A�Y�Y�?�?�?��*�a��R��R������C�C�C�C�	D�	DrDc	��tjd���dd��}tjd��}tt
fD]4}t
t|||||��t
t|||||���5dS)Nr%r#)r�arange�reshaper=rr	r>r?)�nsym�symrLs   r�
test_symmetryz-test_are_validate_args.<locals>.test_symmetry�s~���y��|�|�#�#�A�q�)�)���f�Q�i�i��&�(:�;�	>�	>�A��*�a��c�4��=�=�=��*�a��c�3��=�=�=�=�	>�	>rDc
�
�tjdd��}|dxxdzcc<tjd��}ttfD]}tt||||||���ttt||||��dS)Nr�g��mBr}r!r#)r�fullr=rr	r>r?)�singrBrLs   r�test_singularityz0test_are_validate_args.<locals>.test_singularity�s����w�v�t�$�$���T�
�
�
�a��
�
�
�
�V�A�Y�Y��&�(:�;�	?�	?�A��*�a��R��R��>�>�>�>��j�"6��B��D�I�I�I�I�IrDc��tjdtj��}tjd��}tt
fD]�}t
t|||||��t
t|||||��t
t|||||��t
t|||||��t
t||||||��t
t|||||||����dS)Nr}r")rr��nanr=rr	r>r?)�nmrBrLs   r�test_finitenessz/test_are_validate_args.<locals>.test_finiteness�s���
�W�V�R�V�
$�
$��
�V�A�Y�Y��&�(:�;�	A�	A�A��*�a��R��R�8�8�8��*�a��R��R�8�8�8��*�a��R��R�8�8�8��*�a��R��R�8�8�8��*�a��R��R��<�<�<��*�a��R��R��R�@�@�@�@�
	A�	ArDr^)r�r�r�r�r�s     r�test_are_validate_argsr��sb��>�>�>�D�D�D�>�>�>�J�J�J�	A�	A�	A�	A�	ArDc�`�eZdZejddgddgg��ejddgddgg��ejdd	gd
dgg��fejgd�gd
�gd�gd�g��ejgd�gd�gd�gd�g��ejgd�gd�gd�gd�g��fejddgddgg��ejddgddgg��ejdd gd!dgg��fejd"dgddgg��ejddgddgg��ejdd gd!dgg��fejddgddgg��ejddgddgg��ejdd gd!dgg��fejddgddgg��ejddgddgg��ejddgddgg��fejgd#�gd$�gd%�g��ejdd&gddgg��ejddgd&dgddgg��fejgd'�gd$�gd%�g��ejdd&gdd(gg��ejddgd&d)gd*d+gg��fgZd,�Zd-�Zd.�Zd/S)0�TestSolveSylvesterr!r"rr$r3r5r4r%r&r'r()r0rrr)rr0r*ro)rrr1r�)rrr"r3)r*rrr0)rr0roro)rrr0r9)rrr!r!)rr0rr)rrr0r)rrrr0r)r*r+r,r2r�r1rbr-r.r/r0)r4r!r5)r#r3r6)r$r%r"r#)r4r�y@@y@�y@y@y@@c��t|||��}ttj||��tj||��z|��dSrF)rrrrG�r@rJr��crLs     r�
check_casezTestSolveSylvester.check_case�sB���A�q�!�$�$��!�"�&��A�,�,����1���"=�q�A�A�A�A�ArDc�l�|jD]+}|�|d|d|d���,dS)Nrr!r")rUr�rVs  rrXzTestSolveSylvester.test_cases�sB���J�	7�	7�D��O�O�D��G�T�!�W�d�1�g�6�6�6�6�	7�	7rDc�L�tjddgddgg��}tjdgg��}tjddg���dd��}t|||��}t	|tjddg���dd����dS)Nr0ror*r9r!)rr\r�rrr�s     r�test_trivialzTestSolveSylvester.test_trivial�s����H�s�C�j�3��*�-�.�.���H�s�e�W�����H�c�3�Z� � �(�(��Q�/�/���A�q�!�$�$��!�!�R�X�s�C�j�%9�%9�%A�%A�"�a�%H�%H�I�I�I�I�IrDN)	rYrZr[rr\rUr�rXr�r^rDrr�r��s�������
���A�q�6�A�q�6�"�	#�	#�	���A�q�6�A�q�6�"�	#�	#�	���A�r�7�R��H�%�	&�	&�	(�

���>�>�>�%�%�%�"�?�?��<�<�!�
"�
"�
���#�#�#�%�%�%�"�?�?��<�<�!�
"�
"�
���>�>�>�!�>�>�!�>�>�!�>�>�#�
$�
$�	%�
���F�C�=�8�S�/�2�	3�	3�	���D�"�:��S�z�*�	+�	+�	���F�F�#�g�s�^�4�	5�	5�	7�
���C��:��S�z�*�	+�	+�	���D�!�9�s�C�j�)�	*�	*�	���F�F�#�g�s�^�4�	5�	5�	7�
���F�C�=�8�S�/�2�	3�	3�	���D�!�9�s�C�j�)�	*�	*�	���F�F�#�g�s�^�4�	5�	5�	7�
���F�C�=�8�S�/�2�	3�	3�	���D�!�9�s�C�j�)�	*�	*�	���C��:��c�{�+�	,�	,�	.�
���9�9�9�i�i�i����3�	4�	4�	���A�q�6�A�q�6�"�	#�	#�	���A�q�6�A�q�6�A�q�6�*�	+�	+�	-�
���=�=�=�)�)�)�Y�Y�Y�7�	8�	8�	���A�q�6�A�t�9�%�	&�	&�	���A�r�7�Q��G�b�$�Z�0�	1�	1�	3�Q+
�E�ZB�B�B�7�7�7�J�J�J�J�JrDr�)r�numpyr�
numpy.testingrr�rr>�scipy.linalgrrrrr	r
rr�scipy.sparse._sputilsr
rr r��mark�skipr�r�r�r�r�r^rDr�<module>r�s���	�	�	�	�����3�3�3�3�3�3�
�
�
�
�*�*�*�*�*�*�(�(�(�(�(�(�K�K�K�K�K�K�K�K�A�A�A�A�A�A�A�A�7�7�7�7�7�7�7�7�7�7�(�(�(�(�(�(����TJ�TJ�TJ�TJ�TJ�TJ�TJ�TJ�nB.�B.�B.�H����V��W�W�bA�bA�X�W�bA�J/.�/.�/.�dA.�A.�A.�H.A�.A�.A�b<J�<J�<J�<J�<J�<J�<J�<J�<J�<JrD

Youez - 2016 - github.com/yon3zu
LinuXploit