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

Upload File :
current_dir [ Writeable ] document_root [ Writeable ]

 

Command :


[ Back ]     

Current File : /usr/lib/python3/dist-packages/scipy/integrate/tests/__pycache__/test_bvp.cpython-311.pyc
�

d�c�N����ddlZ	ddlmZn#e$r	ddlmZYnwxYwddlZddlmZmZm	Z	m
Z
ddlmZ
ddlmZddlmZddlmZmZmZmZmZmZd�Zd	�Zd
�Zd�Zd�Zd
�Zd�Zd�Z d�Z!d�Z"d�Z#d�Z$d�Z%d�Z&d�Z'd�Z(d�Z)d�Z*d�Z+d�Z,d�Z-d�Z.d�Z/d�Z0d �Z1d!�Z2d"�Z3d#�Z4d$�Z5d%�Z6d&�Z7d'�Z8d(�Z9d)�Z:d*�Z;d+�Z<d,�Z=d-�Z>d.�Z?d/�Z@d0�ZAd1�ZBd2�ZCd3�ZDd4�ZEd5�ZFd6�ZGd7�ZHdS)8�N)�StringIO)�assert_�assert_array_equal�assert_allclose�assert_equal)�raises)�
coo_matrix)�erf)�modify_mesh�estimate_fun_jac�estimate_bc_jac�compute_jac_indices�construct_global_jac�	solve_bvpc�F�tj|d|df��S�N�r��np�vstack��x�ys  �@/usr/lib/python3/dist-packages/scipy/integrate/tests/test_bvp.py�exp_funrs��
�9�a��d�A�a�D�\�"�"�"�c�r�tjdd|jdf��}d|d<d|d<d|d<d|d<|S)N�r�rrr�rr�rr�rr�r�empty�shape�rr�df_dys   r�exp_fun_jacr(sG���H�a��A�G�A�J�'�(�(�E��E�$�K��E�$�K��E�$�K��E�$�K��Lrc�L�tj|ddz
|df��S�Nrr�r�hstack��ya�ybs  r�exp_bcr0!s#��
�9�b��e�a�i��A��'�(�(�(rc�R�tj|ddz
dz
|df��S)Nrry�?r+r-s  r�exp_bc_complexr2%s'��
�9�b��e�a�i�"�n�b��e�,�-�-�-rc�r�tjddgddgg��}tjddgddgg��}||fSr�r�array�r.r/�dbc_dya�dbc_dybs    r�
exp_bc_jacr9)�Y���h�	
�A��	
�A�����G��h�	
�A��	
�A�����G��G��rc��tj|��tj|dz
��z
dtjd��z
zS)Nrr���)r�exp�rs r�exp_solr?5s4���F�A�2�J�J����A����&�1�r�v�b�z�z�>�:�:rc�`�tj|d|ddz|dzf��S)Nrrrr�rr�ps   r�sl_funrC9s-��
�9�a��d�Q�q�T�1�W�H�q��t�O�,�-�-�-rc���|j\}}tj|d|f��}d|d<d|d<|ddz|d<d|d<tj|d|f��}d|d<d|dz|dz|d<||fS)	Nrrrrr r!r"r<)r%rr$�rrrB�n�mr'�df_dps       r�
sl_fun_jacrI=s����7�D�A�q��H�a��A�Y���E��E�$�K��E�$�K��Q�4��7�(�E�$�K��E�$�K��H�a��A�Y���E��E�$�K��q��t�)�a��d�"�E�$�K��%�<�rc�f�tj|d|d|d|dz
f��Sr*r+�r.r/rBs   r�sl_bcrLLs-��
�9�b��e�R��U�B�q�E�A�a�D�L�1�2�2�2rc��tjd��}d|d<d|d<tjd��}d|d<tjd��}d|d<|||fS)	N)�rrr)rrr!)rNr���)rr�r�zeros)r.r/rBr7r8�dbc_dps      r�	sl_bc_jacrSPs`���h�v���G��G�D�M��G�D�M��h�v���G��G�D�M�
�X�f�
�
�F��F�4�L��G�V�#�#rc�<�tj|d|z��S)Nr)r�sin�rrBs  r�sl_solrW^s��
�6�!�A�$��(���rc�N�tj|d|ddzf��S)Nrr�rrs  r�	emden_funrZbs%��
�9�a��d�Q�q�T�1�W�H�%�&�&�&rc��tjdd|jdf��}d|d<d|d<d|ddzz|d<d|d	<|S)
Nrrrrr ����r!r"r#r&s   r�
emden_fun_jacr^fsS���H�a��A�G�A�J�'�(�(�E��E�$�K��E�$�K��q��t�Q�w�,�E�$�K��E�$�K��Lrc�L�tj|d|ddz
g��S)Nrr�LX�z��?r4r-s  r�emden_bcraos$��
�8�R��U�B�q�E�J�.�/�0�0�0rc�r�tjddgddgg��}tjddgddgg��}||fSr*r4r6s    r�emden_bc_jacrcsr:rc��d|dzdzzdzS)NrrrN�࿩r>s r�	emden_solrgs��
��1��Q��J���rc�*�tj|��S�N�r�
zeros_likers  r�
undefined_funrl�s��
�=����rc�L�tj|d|ddz
g��Sr*r4r-s  r�undefined_bcrn�s#��
�8�R��U�B�q�E�A�I�&�'�'�'rc�P�tj|��}|ddd�|ddd�<|S)Nrrrj)rr�fs   r�big_funrq�s/��
�
�a���A�
�q�t�!�t�W�A�c�c��c�F��Hrc�X�tj|ddd�|ddd�dz
f��S)Nrrr+r-s  r�big_bcrs�s/��
�9�b���1��g�r�#�#�A�#�w��{�+�,�,�,rc�R�tjd|z|jf��}||ddd�<|S)Nr)r�ones�size)rrFrs   r�big_solrw�s/��
���Q����� � �A�
�A�c�c��c�F��Hrc���tj|��}|ddd�|ddd�<|ddz|ddd�z|ddd�<|ddz|ddd�z|ddd�<|S)a� Big version of sl_fun, with two parameters.

    The two differential equations represented by sl_fun are broadcast to the
    number of rows of y, rotating between the parameters p[0] and p[1].
    Here are the differential equations:

        dy[0]/dt = y[1]
        dy[1]/dt = -p[0]**2 * y[0]
        dy[2]/dt = y[3]
        dy[3]/dt = -p[1]**2 * y[2]
        dy[4]/dt = y[5]
        dy[5]/dt = -p[0]**2 * y[4]
        dy[6]/dt = y[7]
        dy[7]/dt = -p[1]**2 * y[6]
        .
        .
        .

    rNrrr]rNrj)rrrBrps    r�big_fun_with_parametersry�s���(	�
�a���A�
�q�t�!�t�W�A�c�c��c�F���t�Q�w�h��3�3�Q�3���A�a�d��d�G���t�Q�w�h��1�4�a�4�� �A�a�d��d�G��Hrc�\�|j\}}tj|||f��}d|td|d��td|d��f<|ddz|td|d��td|d��f<|ddz|td|d��td|d��f<tj|d|f��}d|dz|td|d��z|td|d��df<d|dz|td|d��z|td|d��df<||fS)Nrrrr]rNr<�r%rrQ�rangerEs       r�big_fun_with_parameters_jacr}�s(���7�D�A�q��H�a��A�Y���E�,-�E�%��1�a�.�.�%��1�a�.�.�
(�)�-.�q�T�1�W�H�E�%��1�a�.�.�%��1�a�.�.�
(�)�-.�q�T�1�W�H�E�%��1�a�.�.�%��1�a�.�.�
(�)��H�a��A�Y���E�!�A�a�D�y�1�U�1�a��^�^�+<�<�E�%��1�a�.�.�!�
��!�A�a�D�y�1�U�1�a��^�^�+<�<�E�%��1�a�.�.�!�
���%�<�rc��tj|ddd�|ddd�|d|dz
|d|dz
f��S)NrrrrNr+rKs   r�big_bc_with_parametersr�sG��
�9�b���1��g�r�#�#�A�#�w��1���!���b��e�a��d�l�C�D�D�Drc��|jd}tj|dz|f��}tj|dz|f��}d|t|dz��td|d��f<d|t|dz|��td|d��f<tj|dzdf��}d||df<d||df<d||dzdf<d||dzdf<|||fS)NrrrrOrNr{)r.r/rBrFr7r8rRs       r�big_bc_with_parameters_jacr��s���
����A��h��A��q�z�"�"�G��h��A��q�z�"�"�G�-.�G�E�!�q�&�M�M�5��A�q�>�>�)�*�01�G�E�!�q�&�!���e�A�q�!�n�n�,�-�
�X�q�1�u�a�j�
!�
!�F��F�1�a�4�L��G�A�q�D�M��F�1�q�5�!�8���G�A��E�1�H���G�V�#�#rc��tjtj|d|z��tj|d|z��f��Sr*)rrrUrVs  r�big_sol_with_parametersr��s;��
�9�b�f�Q�q�T�A�X�&�&���q��t�a�x�(8�(8�9�:�:�:rc	�&�d}tj|d||dz|tjdzztjtj|z��zztj|ztjtj|z��zz|zf��S)N���MbP?rr)rr�pi�cosrU)rr�epss   r�	shock_funr��s���
�C�
�9�	�!��
�a��d�(�S�2�5�!�8�^�b�f�R�U�Q�Y�&7�&7�7�
7�
�%�!�)�b�f�R�U�Q�Y�'�'�
'�(�	)�+.�	/����rc�L�tj|ddz|dg��S)Nrrr4r-s  r�shock_bcr��s#��
�8�R��U�Q�Y��1��&�'�'�'rc���d}tjd|z��}tjtj|z��t	||z��t	d|z��zzS)Nr�rr)r�sqrtr�r�r
)rr��ks   r�	shock_solr��sN��
�C�
���C����A�
�6�"�%�!�)���s�1�q�5�z�z�C��A��J�J�6�6�6rc�^�tj|dtj|��g��S)Nr)r�stackrkrs  r�
nonlin_bc_funr��s%��
�8�Q�q�T�2�=��+�+�,�-�-�-rc�^�|\}}|\}}d\}}}}	}
d|z
dz
}|tj|
|z��tj|
|z��z
z}|||zz}
|	|z
dz
}|tj|
|z��tj|
|z��z
z}|||zz
}tj|
|g��S)N)g=
ףp=�?�{�G�z�?�-C��6?��?g33333sC@���?)rr=r5)r.r/�phiA�phipA�phiC�phipC�kappa�ioA�ioC�Vrp�hA�iA�res0�hC�iC�res1s                 r�nonlin_bc_bcr��s����K�D�%��K�D�%�9��E�3��Q��

�T��#��B�	���q��t���r�v�q�b��e�}�}�,�	-�B�
���
��D�
�T��C��B�	���q��t���r�v�q�b��e�}�}�,�	-�B�
��e���D�
�8�T�4�L�!�!�!rc��dd|zz
S)Ng��/��g)���?rfr>s r�
nonlin_bc_solr�s���)�a�-�/�/rc���tjgd�t���}t|tjdg��tjdg����}t	|tjgd�����tjgd�t���}t|tjdgt
���tjgd�����}t	|gd	���dS)
N)rrrN�	��dtyperr)rr�rrNrY�r�)������rrN�r)rrrN)r�r\���r�g��rrrrNr]rYr�)rr5�floatrr�int)r�x_news  r�test_modify_meshr�s���
�����U�+�+�+�A���2�8�Q�C�=�=�"�(�A�3�-�-�8�8�E��u�b�h�'>�'>�'>�?�?�@�@�@�
��"�"�"�%�0�0�0�A���2�8�Q�C�s�3�3�3�R�X�i�i�i�5H�5H�I�I�E��u�I�I�I�J�J�J�J�Jrc��tjddd��}tjd|jdf��}d|d<d|d<tjg��}td�|||��\}}t
||��}t||��t|du��tjdtj	d��}tjd|jdf��}tj
|��|d<tj|��|d<tjdg��}tt|||��\}}t|||��\}}t||��t||��tjddd	��}tjd|jdf��}d
|d<d|d<tjg��}td�|||��\}}t||��}t||��t|du��dS)
NrrrYrr�g{�G�z�?c�"�t||��Sri�rrAs   r�<lambda>z&test_compute_fun_jac.<locals>.<lambda> s��G�A�q�M�M�rr��
r`r�c�"�t||��Sri)rZrAs   rr�z&test_compute_fun_jac.<locals>.<lambda>4s��I�a��O�O�r)r�linspacer$r%r5rr(rrr�rUr�rCrIr^)rrrBr'rH�df_dy_an�df_dp_ans       r�test_compute_fun_jacr�s���
��A�q�!���A�
��!�Q�W�Q�Z��!�!�A��A�a�D��A�a�D�
�����A�#�$A�$A�1�a��K�K�L�E�5��1�a� � �H��E�8�$�$�$��E�T�M����
��A�r�u�a� � �A�
��!�Q�W�Q�Z��!�!�A�
�6�!�9�9�A�a�D�
�6�!�9�9�A�a�D�
��#����A�#�F�A�q�!�4�4�L�E�5�#�A�q�!�,�,��H�h��E�8�$�$�$��E�8�$�$�$�
��A�q�"���A�
��!�Q�W�Q�Z��!�!�A��A�a�D��A�a�D�
�����A�#�$C�$C�Q��1�M�M�L�E�5��Q��"�"�H��E�8�$�$�$��E�T�M�����rc��tjddg��}tjddg��}tjg��}td�|||��\}}}t||��\}}t	||��t	||��t|du��tjddg��}tjddg��}tjdg��}tt|||��\}}}t|||��\}}}t	||��t	||��t	||��tjdd	g��}tjd
dg��}tjg��}td�|||��\}}}t||��\}}t	||��t	||��t|du��dS)
Ng�rr�rNc�"�t||��Sri�r0rKs   rr�z%test_compute_bc_jac.<locals>.<lambda>?s��&��R�.�.�rr�rrO�di���g%@c�"�t||��Sri)rarKs   rr�z%test_compute_bc_jac.<locals>.<lambda>Rs��(�2�r�*�*�r)	rr5r
r9rrrLrSrc)	r.r/rBr7r8rR�
dbc_dya_an�
dbc_dyb_an�	dbc_dp_ans	         r�test_compute_bc_jacr�:s���	��4��)�	�	�B�	��3��(�	�	�B�
�����A�.�(�(�"�b�!� 5� 5��G�W�f�'��B�/�/��J�
��G�Z�(�(�(��G�Z�(�(�(��F�d�N����	��3��(�	�	�B�	��3��)�	�	�B�
��#����A�.�u�b�"�a�@�@��G�W�f�(1�"�b�!�(<�(<�%�J�
�I��G�Z�(�(�(��G�Z�(�(�(��F�I�&�&�&�	��3��*�	�	�B�	��5�$�-�	 �	 �B�
�����A�.�*�*�B��A� 7� 7��G�W�f�)�"�b�1�1��J�
��G�Z�(�(�(��G�Z�(�(�(��F�d�N�����rc
�,�d}d}d}t|||��\}}ttj|��||ff�����}tjgd�gd�gd�gd�gd�gd�gd�gd�gd�gd�g
��}t
||��dS)Nrr])
rrrrrrrrrr)
rrrrrrrrrr)
rrrrrrrrrr)
rrrrrrrrrr)rr	r�	ones_like�toarrayr5r)rFrGr��i�j�s�s_trues       r�test_compute_jac_indicesr�Ys���	�A�	�A�	�A��q�!�Q�'�'�D�A�q��B�L��O�O�a��V�,�-�-�5�5�7�7�A�
�X�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&�&����F��q�&�!�!�!�!�!rc��d}d}d}tddd��\}}tjddd��}tj|��}tjtjtj|z��tjtjtj|z��zf��}tjdg��}t|||��}	|dd�d|zz}
d|dd�dd�f|dd�dd�fzz|dz|	dd�dd�f|	dd�dd�fz
zz
}t|||��\}}
t|
||��\}}t|dd�df|dd�df|��\}}}t|||||||||
||||�
�
}|�
��}d	�}tj||z|z||z|zf��}t|dz
��D]5}||||d��|||z|dz|z�||z|dz|z�f<�6||dzzd
z|ddd�f|ddd�fz
z|d|dz
|zd�df<|||ddd�f|ddd�fzz|dzd
z|ddd�f|ddd�fz
zzz|d|dz
|zd�df<d|d<d|d<d|d
<d|d<t!||d���t#t|||��\}}
t#t|
||��\}}t%t&|dd�df|dd�df|��\}}}t|||||||||
||||�
�
}|�
��}t!||dd���dS)NrrYrrg@rOr��c
���tj|dz|dzzdzdz
d|z|dz|dzzdzdzd|zgd|z|dzz|dz|dzzdzdz
d|z|dzzd|dz|dzzdzz
gg��S)Nr�rrer�r4)�hrBs  r�J_blockz(test_compute_global_jac.<locals>.J_block�s����x�
��T�!�Q�$�Y�r�\�A�
�t�A�v��1��u�Q��T�z�"�}�q�'8�$�q�&�A�
��U�1�a�4�Z��A��a��d��2���)�3�q�5��A��:�q�1�a�4��1��9�R�<�7G�H�
���	rr�)r�r)r�r�)r�r)r�r�绽���|�=��rtolg:�0�yU>�r��atol)rrr��diffrrUr�r�r5rCrIrSrr�rQr|rrr
rL)rFrGr��i_jac�j_jacrr�rrBrp�x_middle�y_middler'rH�df_dy_middle�df_dp_middler7r8rR�Jr��J_truer�s                       r�test_compute_global_jacr�ns���	�A�	�A�	�A�&�q�!�Q�/�/�L�E�5�
��A�q�!���A�
���
�
�A�
�	�2�6�"�%�!�)�$�$�b�e�b�f�R�U�Q�Y�.?�.?�&?�@�A�A�A�
��#����A��q�!�Q���A���"��v��a���H��a����3�B�3��i�!�A�A�A�q�r�r�E�(�*�+�a��c�Q�q�q�q�!�"�"�u�X��!�!�!�S�b�S�&�	�5I�.J�J�H��a��A�&�&�L�E�5�!+�H�h��!B�!B��L�,�(��1�1�1�a�4��!�A�A�A�r�E�(�A�>�>��G�W�f��Q��1�e�U�A�u�l�"�L�'�7�F�	L�	L�A�	�	�	���A�����X�q�1�u�q�y�!�a�%�!�)�,�
-�
-�F�
�1�q�5�\�\�M�M��9@���1��q��t�9L�9L��q�1�u�q�1�u��k�!�1�q�5�1�q�5�A�+�#5�5�6�6�!"�Q��T��!��q��C�R�C��y�1�Q����U�8�/C�!D�F�>�Q��U�a�K��>�2���"#�q�A�a��"��f�I��!�Q�R�R�%��,@�'A�'(�!�t�A�v��1�c�r�c�6��Q�q�!�"�"�u�X�1E�'F�(G�#H�F�1�a�!�e�q�[��?�B����F�4�L��F�4�L��F�5�M��F�6�N��A�v�E�*�*�*�*�#�F�A�q�!�4�4�L�E�5�!1�&�(�H�a�!P�!P��L�,�.�u�a����1��g�q����B��x��K�K��G�W�f��Q��1�e�U�A�u�l�"�L�'�7�F�	L�	L�A�	�	�	���A��A�v�D�t�4�4�4�4�4�4rc	�~�gd�}tjd��}tttt
t||��tjddd��}tjd��}tttt
t||��d�}d�}tjd	|jdf��}ttt||||dg�
��d�}ttt||||��tj	ddgg��}tttt
t|||���dS)
N)rrr�)rrNrrrY)rr]c�"�t||��Srir�rAs   rr�z+test_parameter_validation.<locals>.<lambda>�s��'�!�Q�-�-�rc�"�t||��Srir�rKs   rr�z+test_parameter_validation.<locals>.<lambda>�s��6�"�b�>�>�rr)rBc�*�tjd��S)NrNrPrs  r�wrong_shape_funz2test_parameter_validation.<locals>.wrong_shape_fun�s���x��{�{�r)�S)
rrQ�
assert_raises�
ValueErrorrrr0r�r%r5)rr�fun�bcr�r�s      r�test_parameter_validationr��s�����A�
�����A��*�i��&�!�Q�?�?�?�
��A�q�!���A�
�����A��*�i��&�!�Q�?�?�?�
'�
'�C�	)�	)�B�
��!�Q�W�Q�Z��!�!�A��*�i��b�!�Q�1�#�>�>�>�>�����*�i��"�a��C�C�C�
��1�a�&����A��*�i��&�!�Q�!�D�D�D�D�D�Drc
�.�tjddd��}tjddd��}tjd|jdf��}dtfD�]�}dt
fD�]�}t
tt||||���}t|j
d��t|j��t|j
jd��|�|��}t!|dt#|��d���t||��}|�|d��|z
}|dtj|��zz}	tj|	dzd�	��d
z}
ttj|
dk����ttj|jdk����t!|�|j
��|jdd�
��t!|�|j
d��|jdd�
��������dS)NrrrYr�r��fun_jac�bc_jac��h㈵��>�r���axisr�r�r�r�)rr�rQr%r(r9rrr0r�statusr�successrrv�solrr?�abs�sum�all�
rms_residualsr�yp�r�x_testrr�r�r��sol_test�f_test�r�rel_res�norm_ress           r�test_no_paramsr	�s���
��A�q�!���A�
�[��A�s�
#�
#�F�
��!�Q�W�Q�Z��!�!�A��+�&�O�O���Z�(�	O�	O�F��G�V�Q��7�#)�+�+�+�C�
���Q�'�'�'��C�K� � � �����Q�'�'�'��w�w�v���H��H�Q�K�����t�D�D�D�D��V�X�.�.�F������"�"�V�+�A��1�r�v�f�~�~�-�.�G��v�g�q�j�q�1�1�1�3�6�H��B�F�8�d�?�+�+�,�,�,��B�F�3�,�t�3�4�4�5�5�5��C�G�G�C�E�N�N�C�E��E�J�J�J�J��C�G�G�C�E�1�-�-�s�v�E��N�N�N�N�N�+	O�O�Orc��tjdtjd��}tjdtjd��}tjd|jdf��}dt
fD�]�}dtfD�]�}ttt||dg||���}t|jd��t|j
��t|jjdk��t!|jdgd	�
��|�|��}t!|dt'|dg��d	d	���t||dg��}|�|d��|z
}|dtj|��zz}	tj|	dzd���dz}
ttj|
d
k����ttj|jd
k����t!|�|j��|jdd���t!|�|jd��|jdd�����Ґ��dS)NrrYr�rr��rBr�r�r�rr�r�r�r�r�r�)rr�r�rur%rIrSrrCrLrr�rr�rrvrrBr�rWr�r�r�rrrrs           r�test_with_paramsr�s(��
��A�r�u�a� � �A�
�[��B�E�3�
'�
'�F�
���A�G�A�J�� � �A��*�%�O�O���Y�'�	O�	O�F��F�E�1�a�C�5�'�#)�+�+�+�C�
���Q�'�'�'��C�K� � � ��C�E�J��O�$�$�$��C�E�A�3�T�2�2�2�2��w�w�v���H��H�Q�K�����)<�)<�!%�D�
2�
2�
2�
2��F�H�q�c�2�2�F������"�"�V�+�A��1�r�v�f�~�~�-�.�G��v�g��l��3�3�3�s�:�H��B�F�8�d�?�+�+�,�,�,��B�F�3�,�t�3�4�4�5�5�5��C�G�G�C�E�N�N�C�E��E�J�J�J�J��C�G�G�C�E�1�-�-�s�v�E��N�N�N�N�N�1	O�O�Orc�<�tjddd��}tjddd��}tjd��}d|d<d|d<tjddgdd	gg��}dtfD�]�}dt
fD�]�}t
tt|||||�
��}t|j
d��t|j��t|j
jd��|�|��}t!|dt#|��d���t||��|�|��|zz}|�|d��|z
}	|	dtj|��zz}
tj|
d
zd���dz}ttj|dk����t!|�|j
��|jdd���t!|�|j
d��|jdd���������dS)Nrrr�g�������?r�)rr�r`r�r<)r�r�r�r�r�rr�r�r�r�r�)rr�r$r5r^rcrrZrarr�rr�rrvr�rrg�dotr�r�r�rr)rrrr�r�r�r�rrrrrs            r�test_singular_termr�s
��
��A�q�"���A�
�[��q�#�
&�
&�F�
�����A��A�a�D��A�a�D�
��1�a�&�1�b�'�"�#�#�A��-�(�O�O���\�*�	O�	O�F��I�x��A��G�#)�+�+�+�C�
���Q�'�'�'��C�K� � � �����R�(�(�(��w�w�v���H��H�Q�K��6�):�):��F�F�F�F��v�x�0�0�1�5�5��?�?�V�3K�K�F������"�"�V�+�A��1�r�v�f�~�~�-�.�G��v�g��l��3�3�3�s�:�H��B�F�8�d�?�+�+�,�,�,��C�G�G�C�E�N�N�C�E��E�J�J�J�J��C�G�G�C�E�1�-�-�s�v�E��N�N�N�N�N�'	O�O�Orc
��tjddd��}tjddd��}tjd|jdft���}dt
fD�]�}dtfD�]�}ttt||||���}t|jd��t|j
��|�|��}t|djt#|��d�	��t|djt#|��d�	��t||��}|�|d��|z
}|dtj|��zz}	tjtj|	tj|	��z��d�
��dz}
ttj|
dk����ttj|jdk����t|�|j��|jd
d
���t|�|jd��|jd
d
��������dS)NrrrYr�rr�r�r�r�r�r�r�r�r�)rr�rQr%�complexr(r9rrr2rr�rr�r�r�realr?�imagr�r��conjr�rrrrrs           r�test_complexrs%��	��A�q�!���A�
�[��A�s�
#�
#�F�
��!�Q�W�Q�Z���0�0�0�A��+�&�O�O���Z�(�	O�	O�F��G�^�Q��7�#)�+�+�+�C�
���Q�'�'�'��C�K� � � ��w�w�v���H��H�Q�K�,�g�f�o�o�D�I�I�I�I��H�Q�K�,�g�f�o�o�D�I�I�I�I��V�X�.�.�F������"�"�V�+�A��1�r�v�f�~�~�-�.�G��v�b�g�g����0@�0@�&@�A�A�#$�&�&�&�),�-�H��B�F�8�d�?�+�+�,�,�,��B�F�3�,�t�3�4�4�5�5�5��C�G�G�C�E�N�N�C�E��E�J�J�J�J��C�G�G�C�E�1�-�-�s�v�E��N�N�N�N�N�+	O�O�Orc���tjddd��}tjd|jf��}t	t
t||dd���}t|jd��t|j
��tjddd��}tjd|jf��}t	tt||��}t|jd��t|j
��dS)Nrrrr�rY)�tol�	max_nodes)
rr�rQrvrrr0rr�rr�rlrn)rr�ress   r�
test_failuresr7s���
��A�q�!���A�
��!�Q�V����A�
�G�V�Q��t�q�
A�
A�
A�C����Q�������O����
��A�q�!���A�
��!�Q�V����A�
�M�<��A�
6�
6�C����Q�������O�����rc���d}tjddd��}tjd|z|jf��}t	t
t||��}t|jd��t|j
��|�|��}t|dt||����t||��}|�|d��|z
}|dtj|��zz}tjtj|tj|��z��d���dz}ttj|dk����ttj|jdk����t|�|j��|jd	d	�
��t|�|jd��|jd	d	�
��dS)N�rrrYrr�r�r�r�r�)rr�rQrvrrqrsrr�rr�r�rrwr�r�rrr�rrrr)	rFrrr�rrrrrs	         r�test_big_problemrEs���
�A�
��A�q�!���A�
��!�a�%����!�!�A�
�G�V�Q��
*�
*�C����Q�����C�K�����w�w�q�z�z�H��H�Q�K���A���/�/�/�
�Q��
!�
!�F�����1�
�
���A��1�r�v�f�~�~�%�&�G��v�b�g�g����(8�(8�8�9�9��B�B�B�c�I�H��B�F�8�d�?�#�#�$�$�$��B�F�3�$�t�+�,�,�-�-�-��C�G�G�C�E�N�N�C�E��E�B�B�B�B��C�G�G�C�E�1�%�%�s�v�E��F�F�F�F�F�Frc��d}tjdtjd��}tjdtjd��}tjd|z|jf��}dt
fD�]}dtfD�]}ttt||ddg||���}t|jd��t|j
��t|jddgd	�
��|�|��}t#d|d��D]c}t||t%|ddg��dd	d	���t||dzt%|ddg��dd	d	����dt||ddg��}	|�|d��|	z
}
|
dtj|	��zz}tj|dzd�
��dz}ttj|dk����ttj|jdk����t|�|j��|jdd���t|�|jd��|jdd�������dS)NrrrYr�rr�rrr�r�r]r�r�r�r�)rr�r�rurvr}r�rryrrr�rr�rrBr�r|r�r�r�r�rrrr)
rFrrrr�r�r�r�isolrrrrs
             r� test_big_problem_with_parametersr ]s���
�A�
��A�r�u�a� � �A�
�[��B�E�3�
'�
'�F�
���Q����� � �A��5�6�O�O���7�8�	O�	O�F��3�5K�Q��#�s��W�V�M�M�M�C�
���Q�'�'�'��C�K� � � ��C�E�A�q�6��5�5�5�5��w�w�v���H��a��A���
6�
6������ 7���A�� G� G�� J�%)��6�6�6�6� ����� 2� 7���A�� G� G�� J�%)��6�6�6�6�6�-�V�X��1�v�F�F�F������"�"�V�+�A��1�r�v�f�~�~�-�.�G��v�g��l��3�3�3�s�:�H��B�F�8�d�?�+�+�,�,�,��B�F�3�,�t�3�4�4�5�5�5��C�G�G�C�E�N�N�C�E��E�J�J�J�J��C�G�G�C�E�1�-�-�s�v�E��N�N�N�N�N�7	O�O�Orc��tjddd��}tjddd��}tjd|jf��}t	t
t||��}t|jd��t|j
��t|jjdk��|�|��}t|dt|��dd�	��t||��}|�|d��|z
}|dtj|��zz}tj|dzd�
��dz}ttj|dk����t|�|j��|jd
d
�	��t|�|jd��|jd
d
�	��dS)NrOrrYr�rr�nr�r�r�r�r�r�)rr�rQrvrr�r�rr�rr�rr�rr�r�r�r�rr�	rrrr�rrrrrs	         r�test_shock_layerr$�s���
��B��1���A�
�[��Q��
$�
$�F�
��!�Q�V����A�
�I�x��A�
.�
.�C����Q�����C�K�����C�E�J�������w�w�v���H��H�Q�K��6�!2�!2��D�I�I�I�I�
�v�x�
(�
(�F��������V�#�A��1�r�v�f�~�~�%�&�G��v�g��l��+�+�+�s�2�H��B�F�8�d�?�#�#�$�$�$��C�G�G�C�E�N�N�C�E��E�B�B�B�B��C�G�G�C�E�1�%�%�s�v�E��F�F�F�F�F�Frc�p�tjddd��}|}tjd|jg��}t	t
t||��}t|jd��t|j
��t|jjdk��|�|��}t|dt|��dd���t||��}|�|d��|z
}|dtj|��zz}tj|dzd�	��d
z}ttj|dk����t|�|j��|jdd���t|�|jd��|jdd���dS)
Nrg�������?rYrr�r�r�rr�r�r�r�)rr�rQrvrr�r�rr�rr�rr�rr�r�r�r�rrr#s	         r�test_nonlin_bcr&�sw��
��A�s�A���A�
�F�
��!�Q�V����A�
�M�<��A�
6�
6�C����Q�����C�K�����C�E�J��N�����w�w�v���H��H�Q�K��v�!6�!6�T��M�M�M�M�
�6�8�
,�
,�F��������V�#�A��1�r�v�f�~�~�%�&�G��v�g��l��+�+�+�s�2�H��B�F�8�d�?�#�#�$�$�$��C�G�G�C�E�N�N�C�E��E�B�B�B�B��C�G�G�C�E�1�%�%�s�v�E��F�F�F�F�F�Frc�2�tjddd��}tjd|jdf��}dD]�}tj}t
��t_	ttt|||���}tj�
��}|t_n#|t_wxYwt|j��|dkrt||��|dkrtd|v|��|dkrtd|v|����dS)	NrrrYr)rrr)�verbosez	Solved inzMax residual)
rr�rQr%�sys�stdoutrrrr0�getvaluerr�)rrr(�
old_stdoutr��texts      r�test_verboser.�s
��
��A�q�!���A�
��!�Q�W�Q�Z��!�!�A��2�2���Z�
��Z�Z��
�	$��G�V�Q��7�C�C�C�C��:�&�&�(�(�D�#�C�J�J���C�J�#�#�#�#��������a�<�<���H�d�#�#�#��a�<�<��K�4�'��.�.�.��a�<�<��N�d�*�D�1�1�1��2�2s�!<B*�*B8)Ir)r�ImportError�io�numpyr�
numpy.testingrrrr�pytestrr��scipy.sparser	�
scipy.specialr
�scipy.integrate._bvprrr
rrrrr(r0r2r9r?rCrIrLrSrWrZr^rarcrgrlrnrqrsrwryr}rr�r�r�r�r�r�r�r�r�r�r�r�r�r�r	rrrrrr r$r&r.rfrr�<module>r7s:��
�
�
�
��!�!�!�!�!�!�!����������������������)�)�)�)�)�)�)�)�)�)�)�)�*�*�*�*�*�*�#�#�#�#�#�#�������C�C�C�C�C�C�C�C�C�C�C�C�C�C�C�C�
#�#�#����)�)�)�.�.�.�	�	�	�;�;�;�.�.�.����3�3�3�$�$�$����'�'�'����1�1�1�	�	�	�������(�(�(�
�
�
�-�-�-�
�
�
�
�
�
�6���E�E�E�
$�$�$�$;�;�;�
���(�(�(�7�7�7�.�.�.�
"�"�"�&0�0�0�K�K�K����@���>"�"�"�*25�25�25�jE�E�E�0O�O�O�:O�O�O�BO�O�O�>O�O�O�>���G�G�G�0"O�"O�"O�JG�G�G�0G�G�G�02�2�2�2�2s�
��

Youez - 2016 - github.com/yon3zu
LinuXploit