
    Џkh                     >    d dl Zd dlmZ d dlZd dlmZ d Zd Z	d Z
y)    N)assert_allclose)minimizec                     d } dD ]  }t        dd      D ]p  }t        | ddd||d      }|j                  t        j                  d	g            j                  d	d	      }|j                  j                         }t        ||       r  y )
Nc                     | dz  d| dz  z  fS )N       )xs    d/var/www/teggl/fontify/venv/lib/python3.12/site-packages/scipy/optimize/tests/test_lbfgsb_hessinv.pyfz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m 	$F!jR!%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                  j                         }t        j                  |j                  t        j
                  ddg            |j                  t        j
                  ddg            f      }t        |j                  t        j
                  ddg      j                  dd            j                  d      |j                  t        j
                  ddg                   t        ||       t        ||j                  dd       y )Nr   r   r      c                     t        j                  | t        j                  t        j                  j	                        |             S Nr   dotscipylinalginvr
   H0s    r   r   z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_2r8      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                      d } ddgddggfd}t        |ddd	g
      }t        |j                  j                          | |j                               y )Nc                    | j                   | j                  | j                  | j                  f\  }}}}t	        j
                  | j                  d| j                  i}|}t        |      D ]  }|||   d d t        j                  f   ||   t        j                  d d f   z  ||   z  z
  }|||   d d t        j                  f   ||   t        j                  d d f   z  ||   z  z
  }	t	        j                  |t	        j                  ||	            ||   ||   d d t        j                  f   z  ||   t        j                  d d f   z  z   } |S )Ndtype)skykn_corrsrhor   eyeshaper;   r   newaxisr(   )
selfsyr>   r?   I_arrHkiA1A2s
             r   todense_old_implz test_3.<locals>.todense_old_impl.   s/   !WWdggt||TXXE1gs

5$**5w 	MA1am,qtBJJM/BBSVKKB1am,qtBJJM/BBSVKKBBFF2rN+s1v!Q

]8K/K89!RZZ]8K0L MB		M 	r   r   r   r   r$   c                     t        j                  | t        j                  t        j                  j	                        |             S r&   r'   r,   s    r   r   ztest_3.<locals>.f=   r.   r   r   r/   r   r0   )r   r   r   r   )rK   r   r6   r-   s      @r   test_3rM   ,   sZ     a&1a&	B: 1ZRH=GG$$,,.0@AQAQ0RSr   )numpyr   numpy.testingr   scipy.linalgr)   scipy.optimizer   r"   r8   rM   r	   r   r   <module>rR      s"     )  #$@,Tr   