§
    µ³Ùgq  ã                   ó:   — d dl Zd dlmZ d dlZd dlmZ d„ Zd„ Z	dS )é    N)Úassert_allclose)Úminimizec                  ó4  — d„ } dD ]‘}t          dd¦  «        D ]~}t          | ddd||dœ¬¦  «        }|                     t          j        d	g¦  «        ¦  «                             d	d	¦  «        }|j                             ¦   «         }t          ||¦  «         ŒŒ’d S )
Nc                 ó   — | dz  d| dz  z  fS )Né   é   © )Úxs    ú^/granary/edev/quantum/lib/python3.11/site-packages/scipy/optimize/tests/test_lbfgsb_hessinv.pyÚfztest_1.<locals>.f   s   € Ø�!‰t�Q�q˜!‘t‘Vˆ|Ðó    )g:Œ0âŽyE>gê-�™—q=g#B’¡œÇ;é   é#   TúL-BFGS-B)ÚgtolÚmaxcor)ÚfunÚjacÚmethodÚx0Úoptionsé   )Úranger   Úhess_invÚnpÚarrayÚreshapeÚtodenser   )r   r   r   ÚresultÚH1ÚH2s         r   Útest_1r"      sÂ   € ðð ð ð %ð $ð $ˆÝ˜B ‘m”mð 	$ð 	$ˆFÝ !¨°jÀRØ!%°Ð8Ð8ð:ñ :ô :ˆFð —’¥¤¨1¨#¡¤Ñ/Ô/×7Ò7¸¸!Ñ<Ô<ˆBØ”×(Ò(Ñ*Ô*ˆBå˜B Ñ#Ô#Ð#Ð#ð	$ð$ð $r   c            	      óÆ  ‡— ddgddggŠˆfd„} t          | dddg¬	¦  «        }t          | d
ddg¬	¦  «        }|j                             ¦   «         }t          j        |                     t          j        ddg¦  «        ¦  «        |                     t          j        ddg¦  «        ¦  «        f¦  «        }t          |                     t          j        ddg¦  «                             dd¦  «        ¦  «                             d¦  «        |                     t          j        ddg¦  «        ¦  «        ¦  «         t          ||¦  «         t          ||j        dd¬¦  «         d S )Nr   r   r   é   c                 óŽ   •— t          j        | t          j        t          j                             ‰¦  «        | ¦  «        ¦  «        S )N)r   ÚdotÚscipyÚlinalgÚinv)r
   ÚH0s    €r   r   ztest_2.<locals>.f   s1   ø€ ÝŒv�a�œ¥¤× 0Ò 0°Ñ 4Ô 4°aÑ8Ô8Ñ9Ô9Ð9r   r   é
   r   )r   r   r   ÚBFGSéÿÿÿÿg{®Gáz„?g¸…ëQ¸ž?)ÚrtolÚatol)r   r   r   r   Úvstackr   r   r   )r   Úresult1Úresult2r    r!   r*   s        @r   Útest_2r3      sj  ø€ Øˆaˆ&�1�a�&Ð	€Bð:ð :ð :ð :ð :õ ˜1 Z°R¸°HÐ=Ñ=Ô=€GÝ˜1 V°°R°Ð9Ñ9Ô9€Gà	Ô	×	!Ò	!Ñ	#Ô	#€Bå	ŒØ×Ò�œ 1 a &Ñ)Ô)Ñ*Ô*Ø×Ò�œ 1 a &Ñ)Ô)Ñ*Ô*ð,ñ 
-ô 
-€Bõ Ø×Ò�œ 1 a &Ñ)Ô)×1Ò1°!°AÑ6Ô6Ñ7Ô7×?Ò?ÀÑCÔCØ×Ò�œ 1 a &Ñ)Ô)Ñ*Ô*ñ,ô ,ð ,õ �B˜ÑÔÐÝ�B˜Ô(¨t¸$Ð?Ñ?Ô?Ð?Ð?Ð?r   )
Únumpyr   Únumpy.testingr   Úscipy.linalgr'   Úscipy.optimizer   r"   r3   r	   r   r   ú<module>r8      sn   ðØ Ð Ð Ð Ø )Ð )Ð )Ð )Ð )Ð )Ø Ð Ð Ð Ø #Ð #Ð #Ð #Ð #Ð #ð$ð $ð $ð@ð @ð @ð @ð @r   