
    -i                     >    S SK rS SKJr  S SKrS SKJr  S rS r	S r
g)    N)assert_allclose)minimizec                     S n S H  n[        SS5       Hm  n[        U SSSXS.S9nUR                  [        R                  " S	/5      5      R                  S	S	5      nUR                  R                  5       n[        XE5        Mo     M     g )
Nc                     U S-  SU S-  -  4$ )N       )xs    [/var/www/html/venv/lib/python3.13/site-packages/scipy/optimize/tests/test_lbfgsb_hessinv.pyftest_1.<locals>.f   s    !tQq!tV|    )g:0yE>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mF!jR!%8:F 1#/77!<B((*BB# $ %r   c            	      x  ^ SS/SS//mU4S jn [        U SSS/S	9n[        U S
SS/S	9nUR                  R                  5       n[        R                  " UR                  [        R
                  " SS/5      5      UR                  [        R
                  " SS/5      5      45      n[        UR                  [        R
                  " SS/5      R                  SS5      5      R                  S5      UR                  [        R
                  " SS/5      5      5        [        X45        [        X2R                  SSS9  g )Nr   r   r      c                    > [         R                  " U [         R                  " [        R                  R	                  T5      U 5      5      $ Nr   dotscipylinalginvr
   H0s    r   r   test_2.<locals>.f   .    vva 0 0 4a899r   r   
   r   r   r   r   BFGSg{Gz?gQ?)rtolatol)r   r   r   r   vstackr   r   r   )r   result1result2r!   r"   r.   s        @r   test_2r:      s   a&1a&	B: 1ZRH=G1VR9G				!	!	#B	1a&)*1a&)*, 
-B 1a&)11!A67??C1a&)*, BB((t$?r   c                     ^ S n SS/SS//mU4S jn[        USSS	/S
9n[        UR                  R                  5       U " UR                  5      5        g )Nc                    U R                   U R                  U R                  U R                  4u  pp4[        R
                  " U R                  SU R                  06nUn[        U5       H  nXQU   S S 2[        R                  4   X'   [        R                  S S 24   -  XG   -  -
  nXRU   S S 2[        R                  4   X   [        R                  S S 24   -  XG   -  -
  n	[        R                  " U[        R                  " Xi5      5      XG   X   S S 2[        R                  4   -  X   [        R                  S S 24   -  -   nM     U$ )Ndtype)skykn_corrsrhor   eyeshaper=   r   newaxisr)   )
selfsyr@   rA   I_arrHkiA1A2s
             r   todense_old_impl test_3.<locals>.todense_old_impl.   s   !WWdggt||TXXEg

5$**5wA1am,qtBJJM/BBSVKKB1am,qtBJJM/BBSVKKBBFF2N+svQ

]8K/K89RZZ]8K0L MB	   	r   r   r   r   r%   c                    > [         R                  " U [         R                  " [        R                  R	                  T5      U 5      5      $ r'   r(   r-   s    r   r   test_3.<locals>.f=   r0   r   r   r1   r   r2   )r   r   r   r   )rM   r   r8   r.   s      @r   test_3rQ   ,   sZ     a&1a&	B: 1ZRH=GG$$,,.0@AQAQ0RSr   )numpyr   numpy.testingr   scipy.linalgr*   scipy.optimizer   r#   r:   rQ   r	   r   r   <module>rV      s"     )  #$@,Tr   