← Back
Editing: test_lbfgsb_setulb.cpython-311.pyc
� d�cd � �& � d dl Zd dlmZ d� Zd� ZdS )� N)�_lbfgsbc � � g d�}g d�}g d�}g d�}d}d}d}d}t j g d �� � } t j g d �� � } t j g d �� � } t j g d �� � } t j | |� � r|} | }n]t j | |� � r|} | }nCt j | |� � r|} | }n)t j | |� � r|} | }nt d � � �t j | � � t j |� � fS )z8simplified objective func to test lbfgsb bound violation�g� �?g� �?gleffff�?g�%I�$I�?gbs�E]�?)� �?� r r r )r r gD��W���?r gB�|m��?)r r g�ҥ�?r g�0��Ԕ?g ��ܥ+�@g�5s� �@g2�a �@g`���_ �@)gF���g���� �?g����w�Q�g�����b@gf漍B@)g8���ۿg MpB#�?g���8B92@g�柇P�?@g ���3�)g��8Fgƌ��J�?g���%�G��g;S�9H@g���MG@)g4</�g �dZ���?g�c�4�y?g(�,��7H@g ���@z<Simplified objective function not defined at requested point)�np�array�allclose� ValueError�copy)�x�x0�x1�x2�x3�f0�f1�f2�f3�g0�g1�g2�g3�f�gs �I/usr/lib/python3/dist-packages/scipy/optimize/tests/test_lbfgsb_setulb.py�objfunr s� � � � � �B� #� "� "�B� � � �B� � � �B� �B� �B� �B� �B� �� &� &� &� '� '�B� �� (� (� (� )� )�B� �� &� &� &� '� '�B� �� &� &� &� '� '�B� �{�1�b��� "� �� ��� ��Q�� � � "� �� ��� ��Q�� � � "� �� ��� ��Q�� � � "� �� ���� !�"� "� "� �G�A�J�J���� � � #�#� c �` � d} d}d}d}d}d}t j | fd� � }t j | t j � � }t j | t j � � }t j g d�� � } t j | � � } t j d t j � � } t j | t j � � } t j j j } t j d|z | z d| z z d |z |z z d |z z t j � � }t j d | z | � � }t j d d� � }t j d d� � }t j d| � � }t j d| � � }t j dt j � � }d|dd�< t d� � D ]s}t | � � \ } } t j || |||| | |||||||||||� � | |k � � � r| |k � � � s J d� � ��tdS )zdtest if setulb() violates bounds checks for violation due to floating point rounding error � � g �cAg�h㈵��>� ���� r r � � � � �S60� �, � s STARTN� z9_lbfgsb.setulb() stepped to a point outside of the bounds)r �full�zeros�float64�onesr r r �types�intvar�dtype�ranger �setulb�all)�n�m�factr�pgtol�maxls�iprint�nbd�low_bnd� upper_bndr r r r � fortran_int�wa�iwa�task�csave�lsave�isave�dsave�n_iters r �test_setulb_floatroundrJ E s/ � � �A� �A� �E� �E� �E� �F� �'�1�$�� � �C��h�q�"�*�%�%�G����2�:�&�&�I� �� � � � � �B� ��� � �A� ���b�j�!�!�A� ���B�J� � �A��-�&�,�K� ��!�A�#�a�%�!�A�#�+��1��Q��&��1��,�b�j� 9� 9�B� �(�1�Q�3� � $� $�C� �8�A�u� � �D� �H�Q�� � �E� �H�Q� � $� $�E� �H�R�� %� %�E� �H�R��� $� $�E��D����G���(�(� I� I���a�y�y���1���q�!�W�i��a��E� �b�#�t�V�U�E� �e�U� ,� ,� ,� �Y��#�#�%�%� I�1��<�*<�*<�*>�*>� I� I� G� I� I� I� I�I� Ir )�numpyr �scipy.optimizer r rJ � r r �<module>rN sQ �� � � � � "� "� "� "� "� "�=$� =$� =$�@/I� /I� /I� /I� /Ir
Save File
Cancel