
    bidD                     ,   d dl Z d dlZd dlmZ d dlmZmZ d dlm	Z
 d dlmZmZ d Ze j                  j!                  d       G d d	             Z G d
 de      Z G d de      Ze j                  j!                  d       G d d             Zd ZddZy)    N)default_rng)quadratic_assignmentOptimizeWarning)_calc_score)assert_equalassert_c                  R   g dg dg dg dg dg dg dg dg d	g d
g dg dg} g dg dg dg dg dg dg dg dg dg dg dg dg}t        j                  |       t        j                  |      }} | j                  d   }t        j                  g d      dg|z  z
  }| ||fS )N)r   Z   
   r   r   r   r   r   r   r   r   r   )r
   r   r      r   r   r   r   r   r   r   r   )r   r   r   r   +   r   r   r   r   r   r   r   )r   r   r   r   r   X   r   r   r   r   r   r   )r   r   r   r   r   r      r   r   r   r   r   )r   r   r   r   r   r   r      r   r   r   r   )r   r   r   r   r   r   r   r      r   r   r   )r   r   r   r   r   r   r   r   r   `   r   r   )r   r   r   r   r   r   r   r   r   r      r   )r   r   r   r   r   r   r   r   r   r   r   %   )r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   )r   $   6   r   ;   H   	   "   O      .   _   )r   r   I   #   r
   :      N   r    ,   r   r   )r   r   r      r   a   r!   B   E   =   r   ?   )r   r    r%   r   ]      r   (   r   0   D   U   )r   r
   r   r+   r   @      r   L   r   r2   r3   )r   r!   r&   r,   r1   r   r   7   &   r   r   r   )r   r"   r!   r   r2   r   r   S   r       8   r   )r   r#   r'   r-   r   r4   r6   r   r$   r,      P   )r   r    r(   r   r3   r5   r    r$   r   r1   '   !   )r   r$   r)   r.   r   r   r7   r,   r1   r   F   V   )r   r   r   r/   r2   r   r8   r9   r;   r=   r      )r   r   r*   r0   r3   r   r   r:   r<   r>   r?   r   r   )   r2   r      r            r   r7      r,   r   )nparrayshape)ABnopt_perms       i/home/cdr/jupyterlab/.venv/lib/python3.12/site-packages/scipy/optimize/tests/test_quadratic_assignment.pychr12crN      s    ......-.-.--	A 	777756577757	A 88A;qA	
Axx?@A37JHa>    z1ignore:The NumPy global RNG was seeded by callingc                   .    e Zd ZdZd Zd Zd Zd Zd Zy)QAPCommonTestsz6
    Base class for `quadratic_assignment` tests.
    c                    t         j                  j                  d      }g dg dg dg dg}g dg dg dg dg}t        ||| j                  |d	d
      }t        |j                  d       t        |j                  t        j                  g d             t        ||| j                  |dd
      }| j                  dk(  rBt        |j                  d       t        |j                  t        j                  g d             nAt        |j                  d       t        |j                  t        j                  g d             t        ||| j                  |dd
       y )N   mX)Ec )r   rA   rB   rE   )r   r   r   rE   )r   r   r   r   )r   r   r   r   )r   rB   rE   rB   )r   rE   r   rE   )r   r   rE   r   Frngmaximizemethodoptionsr   r   rE   rA   r   Tfaqr   )r   rE   rA   r   r-   )r   rA   r   rE   )	rF   randomr   r   rX   r   funcol_indrG   selfrU   rI   rJ   ress        rM   test_accuracy_1zQAPCommonTests.test_accuracy_1;   s   ii##$78
 
 #1a36E+JL 	SWWb!S[["((<"89"1a36D+IK ;;%"%bhh|&<="%bhh|&<=Q$++-0d%C	ErO   c                    t         j                  j                  d      }t        j                  g dg dg dg dg      }t        j                  g dg dg dg d	g      }t	        ||| j
                  |d
d      }| j
                  dk(  rBt        |j                  d       t        |j                  t        j                  g d             nAt        |j                  d       t        |j                  t        j                  g d             t	        ||| j
                  |dd      }t        |j                  d       t        |j                  t        j                  g d             y )NrS   )r   r2   rC   rD   )r2   r   r2   r   )rC   r2   r   rE   )rD   r   rE   r   )r   r   rC   rB   )r   r   r2   rE   )rC   r2   r   r2   )rB   rE   r2   r   FrT   rW   r[      )r   r   rA   rE      rZ   Ti  )rE   rA   r   r   )	rF   r\   r   rG   r   rX   r   r]   r^   r_   s        rM   test_accuracy_2zQAPCommonTests.test_accuracy_2a   s   ii##$78HHl"""$ %
 HHl"""$ %
 #1a36E+JL ;;%#&bhh|&<=#&bhh|&<="1a36D+IK 	SWWc"S[["((<"89rO   c                 @   t         j                  j                  d      }t               \  }}}t	        ||| j
                  d|i      }t        d|j                  cxk  xr dk  nc        t        |j                  t        |||j                               t	        ||| j
                  |dd      }t        d|j                  cxk  xr d	k  nc        t        |j                  t        |||j                               t        j                  g d
      }t        j                  |||   g      j                  }t	        ||| j
                  ||d      }t        d|j                  cxk  xr dk  nc        t        |j                  |   ||          t        j                  t        j                  t        |            |g      j                  }t	        ||| j
                  ||d      }t        |j                  |d d df   j                         t        |j                  d       t        |j                   d       t        j"                  d      }t	        ||| j
                  d|i      }t        |j                   d       t        |j                  d       y )NrS   rU   rW   +  R  TrT   i! iL )rB   rC   r   )partial_matchrU   r   r   )r   r   )rF   r\   r   rN   r   rX   r   r]   r   _scorer^   rG   asarrayTarangelennitempty)	r`   rU   rI   rJ   rL   ra   	seed_costseedrq   s	            rM   test_accuracy_3zQAPCommonTests.test_accuracy_3   s   ii##$781h #1a,13<9(5()SWWfQ3;;78 #1a36D+IK(5()SWWfQ3;;78 HHZ(	zz9hy&9:;=="1a=A#+NP(5()S[[+Xi-@A zz299SV,h78::"1a=A#+NPS[[$q!t*,,/SWWe$SWWa   "5%,13<9SWWa SWWa rO   c                     t               \  }}}t        j                  t              5  t	        ||| j
                  ddi       d d d        y # 1 sw Y   y xY w)Nz	ekki-ekkiTrW   )rN   pytestwarnsr   r   rX   )r`   rI   rJ   rL   s       rM   test_unknown_optionsz#QAPCommonTests.test_unknown_options   sJ    1h\\/* 	> Adkk*5t)<>	> 	> 	>s   AAc           	         t        j                  d      j                  d      }t        j                  d      j                  d      }t        j                  t
        d      5  t         j                  j                  d      }t        ||| j                  |dd       d d d        t        j                  t        d	      5  t         j                  j                  d       t        ||| j                  d
di       d d d        t        j                  t        d      5  t        ||| j                  ddd       d d d        y # 1 sw Y   xY w# 1 sw Y   VxY w# 1 sw Y   y xY w)Nr   rB   rB   zUse of `RandomState`*matchr   FrT   rW   z The NumPy global RNG was seeded*rV   z!The behavior when the rng option*)rF   rn   reshaperv   rw   DeprecationWarningr\   RandomStater   rX   FutureWarningrs   )r`   rI   rJ   rU   s       rM    test_deprecation_future_warningsz/QAPCommonTests.test_deprecation_future_warnings   s%   IIbM!!&)IIbM!!&)\\,4KL 	J))''*C Adkk14%)HJ	J
 \\-/QR 	>IINN1 Adkk*4e)<>	>
 \\-/RS 	H Adkk12)FH	H 	H	J 	J
	> 	>
	H 	Hs$   $<E;E!EEEE(N)	__name__
__module____qualname____doc__rb   rf   rt   rx   r    rO   rM   rQ   rQ   2   s$     EL:<%!N>HrO   rQ   c                       e Zd ZdZd Zd Zy)TestFAQr[   c                 Z   t         j                  j                  d      }t               \  }}}t	        |      }t        ||ddi      }t        |j                  d       t        ||ddi      }t        d|j                  cxk  xr dk  nc        t        |||d	d
      }t        d|j                  cxk  xr dk  nc        t        j                  ||f      t        |      z  }t        |      }t        ||d|i      }t        d|j                  cxk  xr
 dk         y c        y )NrS   maxiterr2   rY   shuffle_inputTrh   ri   
randomized)rU   P0r   )rF   r\   r   rN   ro   r   r   rp   r   r]   onesfloat_doubly_stochastic)r`   rU   rI   rJ   rL   rK   ra   Ks           rM   test_optionszTestFAQ.test_options   s    ii##$781hF #1a)Q@SWWa  #1a/41HI(5() #1aL1QR(5() GGQFOeAh&q!"1a$;(5()()rO   c           	      
   t        j                  d      }|}t        j                  t        d      5  t        ||ddi       d d d        t        j                  t        d      5  t        ||dd	i       d d d        t        j                  t        d
      5  t        ||dd	i       d d d        t        j                  t              5  t        ||ddi       d d d        t        j                  t        d      5  t        t        j                  d      t        j                  d      dt        j                  d      i       d d d        g dg dg dg}t        j                  t        d      5  t        t        j                  d      t        j                  d      d|i       d d d        y # 1 sw Y   wxY w# 1 sw Y   PxY w# 1 sw Y   )xY w# 1 sw Y   xY w# 1 sw Y   xY w# 1 sw Y   y xY w)NrE   zInvalid 'P0' parameterr{   r   r\   r   z$'maxiter' must be a positive integerr   z'tol' must be a positive floattolg      ?z1`P0` matrix must have shape m' x m', where m'=n-mrB   rA   rA   )g?皙?333333?)r   g333333?r   )r   r   gffffff?z%`P0` matrix must be doubly stochasticrA   )rF   identityrv   raises
ValueErrorr   	TypeErrorr   )r`   rI   rJ   r   s       rM   test_specific_input_validationz&TestFAQ.test_specific_input_validation   s   KKN ]]:-EF 	A Ah/?@	A]]"HJ 	@ A	2?	@ ]]:-MN 	< Ar{;	<
 ]]9% 	A A	3/?@	A ]]IK 	 !AArwwv/	  ]]"IK 	 AAq		 	5	A 	A	@ 	@	< 	<
	A 	A	 		 	sH   F9'GGG A
G-97G99GGG G*-G69HN)r   r   r   rX   r   r   r   rO   rM   r   r      s    F*2'rO   r   c                   "    e Zd ZdZd Zd Zd Zy)Test2opt2optc                    d}t        d      }|j                  ||f      }|j                  ||f      }t        ||| j                  d|i      }t        d      }|j                  ||f      }|j                  ||f      }t        ||| j                  d|i      }t	        |j
                  |j
                         y )N   i<2sizerU   rW   )r   r\   r   rX   r   rp   )r`   rK   rU   rI   rJ   res1res2s          rM   test_deterministiczTest2opt.test_deterministic
  s    (#JJQFJ#JJQFJ##AquclS(#JJQFJ#JJQFJ##AquclSTXXtxx(rO   c           	         d}t         j                  j                  d      }|j                  ||f      }|j                  ||f      }t        ||| j                  d|i      }t        j
                  t        j                  d      |j                  g      j                  }t        ||| j                  ||d      }ddg}t        j
                  t        j                  d      |   |j                  |   g      j                  }	t        ||| j                  |||	d	      }
t        |j                  ||d
z   z  dz  k7         t        |j                  ||d
z   z  dz         t        |
j                  |dz
  |d
z
  z  dz         y )Nr2   rS   r   rU   rW   )rU   partial_guessrE   rB   )rU   r   rj   r   )rF   r\   r   r   rX   rG   rn   r^   rm   r   rp   r   )r`   rK   rU   rI   rJ   r   guessr   fixr|   res3s              rM   test_partial_guesszTest2opt.test_partial_guess  sS   ii##$78JJQFJ#JJQFJ##Aq-2CL:"))A,5688#Aq47%,PR!f"))A,s+T\\#->?@BB#Aq47%>C-EF 	AqsGAI%&TXXq!A#wqy)TXX!ac{1}-rO   c                    t         }t        j                  t        d      5  t	        t        j                  d      t        j                  d      | j                  d |dd      i       d d d        t        j                  t        d      5  t	        t        j                  d      t        j                  d      | j                  dt        dd      i       d d d        t        j                  t        d	      5  t	        t        j                  d      t        j                  d      | j                  dt
        j                  j                  ddd      i       d d d        t        j                  t        d
      5  t	        t        j                  d      t        j                  d      | j                  ddt        dd      z  i       d d d        t        j                  t        d      5  t	        t        j                  d      t        j                  d      | j                  ddt        dd      z  i       d d d        t        j                  t        d      5  t	        t        j                  d      t        j                  d      | j                  dt        j                  d      i       d d d        y # 1 sw Y   ZxY w# 1 sw Y   xY w# 1 sw Y   xY w# 1 sw Y   xY w# 1 sw Y   xY w# 1 sw Y   y xY w)Nz0`partial_guess` can have only as many entries asr{   rA   r   r2   rE   rW   z%`partial_guess` must have two columnsz%`partial_guess` must have exactly twoz%`partial_guess` must contain only posr   z0`partial_guess` entries must be less than numberrB   z-`partial_guess` column entries must be uniquerE   rE   )_range_matrixrv   r   r   r   rF   r   rX   r\   randr   )r`   _rms     rM   r   z'Test2opt.test_specific_input_validation-  sR   ]]HJ 	G !QQ(,*93q!9)EG	G ]]"IK 	 AAt{{(-1*=>	 ]]"IK 	 AAt{{("))..Aq*AB	 ]]"IK 	 AAt{{("}Q/B*BC	 ]]HJ 	 !AAt{{(!mAq.A*AB	 ]]EG 	 !AAt{{("''&/:	 	I	G 	G	 		 		 		 		 	sM   A	KAK=A!KAK*3AK7%ALKKK'*K47L LN)r   r   r   rX   r   r   r   r   rO   rM   r   r     s    F).*-rO   r   c                       e Zd Zd Zy)TestQAPOncec                    t        d      }t        j                  t        d      5  t	        |j                  d      |j                  d             d d d        t        j                  t        d      5  t	        |j                  d      |j                  d             d d d        t        j                  t        d      5  t	        |j                  d      |j                  d             d d d        t        j                  t        d	      5  t	        |j                  d      |j                  d
             d d d        t        }t        j                  t        d      5  t	        t        j                  d      t        j                  d      d |dd      i       d d d        t        j                  t        d      5  t	        t        j                  d      t        j                  d      dt        dd      i       d d d        t        j                  t        d      5  t	        t        j                  d      t        j                  d      dt        j
                  j                  ddd      i       d d d        t        j                  t        d      5  t	        t        j                  d      t        j                  d      ddt        dd      z  i       d d d        t        j                  t        d      5  t	        t        j                  d      t        j                  d      ddt        dd      z  i       d d d        t        j                  t        d      5  t	        t        j                  d      t        j                  d      dt        j                  d      i       d d d        y # 1 sw Y   txY w# 1 sw Y   3xY w# 1 sw Y   xY w# 1 sw Y   xY w# 1 sw Y   WxY w# 1 sw Y    xY w# 1 sw Y   xY w# 1 sw Y   :xY w# 1 sw Y   xY w# 1 sw Y   y xY w)Ninn z`A` must be squarer{   )rA   rB   r   z`B` must be squarez!`A` and `B` must have exactly two)rA   rA   rA   z*`A` and `B` matrices must be of equal sizerz   z.`partial_match` can have only as many seeds asrA   rj   r2   rE   r   z%`partial_match` must have two columnsz%`partial_match` must have exactly twoz%`partial_match` must contain only posr   z0`partial_match` entries must be less than numberrB   z-`partial_match` column entries must be uniquer   )r   rv   r   r   r   r\   r   rF   r   r   r   )r`   rU   r   s      rM   test_common_input_validationz(TestQAPOnce.test_common_input_validationa  sY   (#]]:-AB 	 

6"

6"	
 ]]:-AB 	 

6"

6"	 ]]"EG 	 

9%

9%	 ]]BD 	 !

6"

6"	 ]]FH 	G !QQ*93q!9)EG	G ]]"IK 	 AA(-1*=>	 ]]"IK 	 AA("))..Aq*AB	 ]]"IK 	 AA("}Q/B*BC	 ]]HJ 	 !AA(!mAq.A*AB	 ]]EG 	 !AA("''&/:	 	}	 	
	 		 		 		G 	G	 		 		 		 		 	s|   +O5+O+O+O+%>O8AP*AP#AP
AP,1A
P8OOO(+O58PPPP),P58QN)r   r   r   r   r   rO   rM   r   r   ]  s    GrO   r   c                     t        j                  | |f      }t        |      D ]  }t        j                  |       |d d |f<     |S )N)rF   zerosrangern   )abmatis       rM   r   r     sA    
((Aq6
C1X !IIaLAqD	!JrO   c                    d}d| j                  d      z  }d| |z  z  }| }t        |      D ]  }t        j                  |j                  d      dz
        |k  j	                         r<t        j                  |j                  d      dz
        |k  j	                         r |S d|| z  z  }d| |z  z  }|d d d f   | z  |z  } |S )Ni  r   r   )axis)sumr   rF   absall)Pr   max_itercrP_epsits          rM   r   r     s     H	AEEqEMA	QUAEHo #VVEII1I%)*S0557		q	)A-.499; L	 QKQK!T'
Q"# LrO   )gMbP?)rv   numpyrF   numpy.randomr   scipy.optimizer   r   scipy.optimize._qapr   rk   numpy.testingr   r   rN   markfilterwarningsrQ   r   r   r   r   r   r   rO   rM   <module>r      s      $ @ 5 /"J OPKH KH QKH\Cn CLS~ Sl OPJ J QJZrO   