| 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/optimize/_lsq/__pycache__/ |
Upload File : |
�
d�cK � �L � d Z ddlZddlmZmZ ddlmZ ddlm Z m
Z
d� Z d d�Z dS )
z)Bounded-variable least-squares algorithm.� N)�norm�lstsq)�OptimizeResult� )�print_header_linear�print_iteration_linearc �z � | |z }|dk }t j | | � � ||< t j |� � S )z0Compute the maximum violation of KKT conditions.r )�np�abs�max)�g�on_bound�g_kkt�free_sets �:/usr/lib/python3/dist-packages/scipy/optimize/_lsq/bvls.py�compute_kkt_optimalityr s9 � �
��L�E��1�}�H��f�Q�x�[�)�)�E�(�O�
�6�%�=�=�� c
�r
� | j \ } }
|� � � } t j |
� � } | |k }
||
| |
< d| |
< | |k }
||
| |
< d| |
< | dk }| }t j |� � \ }| � | � � |z
}dt j ||� � z }|}| j � |� � }d }d }d}|dk rt � � |j dk �r�|dk r#t || � � }t |||||� � |dz
}| | � � � }| d d �|f }|| � | |z � � z
}t |||�� � d }||| k } ||| k } | | z } t j | � � r || } || | | < d|| < d| | < t j | � � r || } || | | < d|| < d| | < || } || | | < | � | � � |z
}dt j ||� � z } || z
}| }| j � |� � }t | | |z
� � }t j | � � r
|| }nn |j dk ���|�|
}||z
}d }!t || � � }t ||� � D �]O}|dk rt |||||� � ||k rd}!|!� �n(t j || z � � }"d| |"< | dk }| }t j |� � \ }| | }#|#� � � }|| }$|| }%| d d �|f }|| � | |z � � z
}t |||�� � d }t j ||$k � � \ } t j ||%k � � \ } t j | | f� � } | j dk r�t j |$| |#| z
|%| |#| z
f� � || |#| z
z }&t j |&� � }'| |' }(|&|' })|#d|)z
z }#|#|)|z z
}#|#| |< |'| j k r d| ||( < nd| ||( < n|}#|#| |< n���t |#|z
� � }| � | � � |z
}dt j ||� � z } || z
}|||z k rd}!| }| j � |� � }t || � � }��Q|!�d}!t% | |||| |dz |!|�� � S ) N���r r g �?� )�rcondT)�x�fun�cost�
optimality� active_mask�nit�status� initial_cost)�shape�copyr
�zeros�nonzero�dot�Tr �sizer r r �anyr �range�argmax�hstack�argminr )*�A�b�x_lsq�lb�ub�tol�max_iter�verboser �m�nr r �maskr �
active_set�rr r r
� cost_change� step_norm� iterationr �
x_free_old�A_free�b_free�z�lbv�ubv�v�ind�cost_new�termination_status� move_to_free�x_free�lb_free�ub_free�alphas�i�i_free�alphas* r �bvlsrN s� � � �7�D�A�q�
�
�
� � �A��x��{�{�H� ��7�D���h�A�d�G��H�T�N� ��7�D���h�A�d�G��H�T�N��1�}�H���J��
�8�$�$�I�H� ���a���1� �A� ����1� � � �D��L� �����
�
�A��K��I��I��!�|�|�� � � � �-�!�
�
� �a�<�<�/��8�<�<�J� "�9�d�K��#-�
/�
/�
/� �Q�� ��x�[�%�%�'�'�
��1�1�1�h�;� ���Q�U�U�1�z�>�*�*�*�� �&�&�� .� .� .�q� 1���"�X�,� ���"�X�,� �� �#�I��
�6�#�;�;� ��3�-�C���W�A�c�F� "�J�s�O� �H�S�M�
�6�#�;�;� ��3�-�C���W�A�c�F� "�J�s�O� �H�S�M���r�l���A�2���#��
�E�E�!�H�H�q�L������1� � �%���X�o� ���
�C�G�G�A�J�J����8��z�1�2�2� �
�6�!�9�9� � ���|�H�H� �U �-�!�
�
�X ���� � ��H� �� (��8�4�4�J��9�h�/�/� B9� B9� � �a�<�<� "�9�d�K�#,�j�
:�
:�
:� �� � �!"� � � )� �E��y��X� �.�.� �!"�� � �' � �1�}�H�"��J��
�8�,�,�I�H��x�[�F� ��� � �J���l�G���l�G��q�q�q�(�{�^�F�����q�:�~�.�.�.�F��f�f�E�2�2�2�1�5�A��:�a�'�k�*�*�D�C��:�a�'�k�*�*�D�C�� �3��*�%�%�A��v��z�z����C�L�6�#�;�.��C�L�6�#�;�.�$0� 1� 1�45�a�D�6�!�9�4D�F�� �I�f�%�%���1��� �q� ���!�e�)�#���%�!�)�#�� $��(� ��s�x�<�<�13�H�X�f� -�.�.�12�H�X�f� -�.�.��� $��(� ��O' �R ��*�,�-�-� �
�E�E�!�H�H�q�L������1� � �%���X�o� � ��t�� #� #�!"� ���
�C�G�G�A�J�J��+�A�x�8�8�
�
��!� �� �
���*�(� ��M�"4�!� #� #� #� #r )N)
�__doc__�numpyr
� numpy.linalgr r �scipy.optimizer �commonr r r rN � r r �<module>rU s� �� /� /� � � � � $� $� $� $� $� $� $� $� )� )� )� )� )� )� ?� ?� ?� ?� ?� ?� ?� ?�� � �f#� f#� f#� f#� f#� f#r