
    Ki                     x   S SK r S SKJr  S SKrS SKrS SKJr  S SKJr  S SK	J
r
JrJrJrJrJr  S SKJr  S SKJr  S SKJrJrJrJrJrJrJrJrJrJrJrJ r J!r!J"r"J#r#J$r$J%r%J&r&J'r'J(r(J)r)J*r*J+r+J,r,J-r-J.r.J/r/J0r0J1r1  S S	K2J3r3  S S
K4J5r5J6r6J7r7J8r8J9r9  S SK:J;r;J<r<J=r=J>r>J?r?  S SK@JArAJBrBJCrCJDrDJErE  S SKFJGrGJHrH  S rI\R                  R                  S\9" 5       \6S9\R                  R                  S/ SQ5      S 5       5       rL\R                  R                  S\B5      \R                  R                  S\C5      \R                  R                  S\A5      \R                  R                  S\D5      S 5       5       5       5       rM\?" \NS9\R                  R                  S\5      S 5       5       rOS rP\R                  R                  S\*\./5      S 5       rQS rR\R                  " SS 5      R                  S!S"S#9SS$.rUS% rV\R                  R                  S&\*S'0 4\*\\U4\*S(\U4\.S)S*S04\.\VS+S,04/5      \R                  R                  S-\R                  \R                  \Y/5      S. 5       5       rZ\R                  R                  S\9" 5       \6S9\R                  R                  S&\*S'0 4\*S/0 4\.S)S*S04\.\VS+S,04\.S0S+S,04/5      S1 5       5       r[S2 r\\R                  R                  S/ S3Q5      \R                  R                  S\D5      S4 5       5       r]\R                  R                  S\9" 5       \6S9\R                  R                  S/ S5Q5      S6 5       5       r^S7 r_S8 r`\R                  R                  S9\R                  " 5       5      \R                  R                  S\D5      S: 5       5       rbS; rc\R                  R                  \R                  R                  S<\E5      \R                  R                  S\D5      S= 5       5       5       reS> rfS? rgS@ rh\R                  R                  SASB SC SD /\D V s/ s H	  n U 4SE jPM     sn -   \E Vs/ s H	  nU4SF jPM     sn-   5      SG 5       ri\R                  R                  SHSI \jSJ4SK \jSL4SM \jSN4SO \kSP4SQ \kSR4SS \kST4/5      SU 5       rlSSV jrm\R                  R                  SSW5      SX 5       rn\R                  R                  SSW5      SY 5       ro\R                  R                  SZ5      S[ 5       rq\R                  R                  S\\R                  /\D-   S]/\D Vs/ s H  o"R                  PM     sn-   S9\R                  R                  S^\R                  /\D-   S]/\D Vs/ s H  o"R                  PM     sn-   S9S_ 5       5       rt\R                  R                  S^\R                  /\D-   S]/\D Vs/ s H  o"R                  PM     sn-   S9S` 5       ru\R                  R                  SaSbS"/5      Sc 5       rvSd rw\R                  R                  S\\R                  /\D-   S]/\D Vs/ s H  o"R                  PM     sn-   S9\R                  R                  S^\R                  /\D-   S]/\D Vs/ s H  o"R                  PM     sn-   S9Se 5       5       rx\R                  R                  S\\R                  /\D-   S]/\D Vs/ s H  o"R                  PM     sn-   S9Sf 5       ry\R                  R                  Sg/ ShQ5      \R                  R                  S\\R                  /\D-   S]/\D Vs/ s H  o"R                  PM     sn-   S9\R                  R                  S^\R                  /\D-   S]/\D Vs/ s H  o"R                  PM     sn-   S9Si 5       5       5       rz\R                  R                  Sg/ ShQ5      \R                  R                  S\\R                  /\D-   S]/\D Vs/ s H  o"R                  PM     sn-   S9Sj 5       5       r{\R                  R                  Sk\R                  SlSm4\R                  " \R                  SnSo\R                  R                  SpSq9Sr9/5      \R                  R                  SsSSt/5      Su 5       5       r~\R                  R                  SvSbS"/5      Sw 5       r\R                  R                  Sx\R                  " \GR                   S //5      \R                  " S \GR                   * //5      /5      \R                  R                  Sy\R                  " \GR                   S //5      \R                  " S \GR                   * //5      S/5      Sz 5       5       r\R                  R                  S{\R                  " S S/SS //5      \GR                  " S|5      \GR                  4\R                  " S S/S\GR                  //5      \GR                  " S|5      \GR                  4\R                  " \GR                  S/S\GR                  //5      \GR                  \GR                  4\R                  " \GR                  S/\GR                  S //5      \GR                  " S|5      \GR                  4\R                  " S \GR                  /S\GR                  //5      \GR                  " S|5      \GR                  4\R                  " S S/SS //5      \GR                  " S|5      S}4\R                  " S S/SS}//5      \GR                  " S|5      S}4\R                  " S}S/SS}//5      \GR                  S}4\R                  " S}S/S}S //5      \GR                  " S|5      S}4\R                  " S S}/SS}//5      \GR                  " S|5      S}4/
5      S~ 5       r\R                  R                  S\GR                  S}/5      S 5       r\R                  R                  S\GR                  S}/5      S 5       r\R                  R                  S\GR                  S}/5      S 5       rS rS rS rS rS rS r\R                  R                  S\#\/\0\"\1\45      S 5       r\R                  R                  S\#\/\0\"\1\45      \R                  R                  S\D5      S 5       5       rS rS rS r\R                  R                  SS\#4S\4/5      \R                  R                  S\D5      S 5       5       r\R                  R                  S\D5      S 5       rS rS rS rS r\R                  R                  S\D5      S 5       rS rS rS r\R                  R                  SSS|/5      \R                  R                  SSS/5      \R                  R                  S\*\-/5      S 5       5       5       r\R                  R                  SSS/5      S 5       r\R                  R                  S/ SQ5      \R                  R                  SSbS"/SS/S9S 5       5       r\R                  R                  S\*\+\,/5      S 5       rS r\R                  R                  S/ SQS/ SQ/ SQ/ SQ/4/ SQSS/SS/SS/SS//4/5      S 5       rS r\R                  R                  \R                  R                  S\D5      S 5       5       rgs  sn f s  snf s  snf s  snf s  snf s  snf s  snf s  snf s  snf s  snf s  snf )    N)GeneratorType)linalg)issparse)cdist	cityblockcosine	minkowskipdist
squareform)config_context)DataConversionWarning)PAIRED_DISTANCESPAIRWISE_BOOLEAN_FUNCTIONSPAIRWISE_DISTANCE_FUNCTIONSPAIRWISE_KERNEL_FUNCTIONS_euclidean_distances_upcastadditive_chi2_kernelcheck_paired_arrayscheck_pairwise_arrayschi2_kernelcosine_distancescosine_similarityeuclidean_distanceshaversine_distanceslaplacian_kernellinear_kernelmanhattan_distancesnan_euclidean_distancespaired_cosine_distancespaired_distancespaired_euclidean_distancespaired_manhattan_distancespairwise_distancespairwise_distances_argminpairwise_distances_argmin_minpairwise_distances_chunkedpairwise_kernelspolynomial_kernel
rbf_kernelsigmoid_kernel)	normalize)_convert_to_numpy_get_namespace_device_dtype_idsget_namespacexpx)yield_namespace_device_dtype_combinations)_array_api_for_testsassert_allcloseassert_almost_equalassert_array_equalignore_warnings)BSR_CONTAINERSCOO_CONTAINERSCSC_CONTAINERSCSR_CONTAINERSDOK_CONTAINERS)Paralleldelayedc                 R	   [         R                  R                  S5      nUR                  S5      R	                  U SS9n[        USS9n[        U5      n[        X45        UR                  UR                  s=:X  a  U :X  d   e   eUR                  S5      R	                  U SS9n[        X%SS9n[        X%5      n[        X45        UR                  UR                  s=:X  a  U :X  d   e   eUR                  S5      R	                  U SS9nUR                  S5      R	                  U SS9n[         R                  US'   [         R                  US'   [        XgS	S9n[        Xg5      n	[        X5        UR                  U	R                  s=:X  a  U :X  d   e   e[        U V
Vs/ s H  n
[        U
 Vs/ s H  oPM     sn5      PM      snn
5      n[        U V
Vs/ s H  n
[        U
 Vs/ s H  oPM     sn5      PM      snn
5      n[        XSS9n[        X45        UR                  UR                  s=:X  a  U :X  d   e   eUR                  S
5      R	                  U SS9nUS S 2S4   S-
  S-  [         R                  -  S-  US S 2S4'   US S 2S4   S-
  S-  [         R                  -  US S 2S4'   [        USS9n[        U5      n[        X45        UR                  S5      R	                  U SS9nUS S 2S4   S-
  S-  [         R                  -  S-  US S 2S4'   US S 2S4   S-
  S-  [         R                  -  US S 2S4'   [        X%SS9n[        X%5      n[        X45        [        USS9n[        U[        S9nUR                  S   UR                  S   :X  d   eUR                  S   UR                  S   :X  d   e[        X45        [        X%SS9n[        X%[        S9nUR                  S   UR                  S   :X  d   eUR                  S   UR                  S   :X  d   e[        X45        [        X%SS9n[        X%[         S9nUR                  S   UR                  S   :X  d   eUR                  S   UR                  S   :X  d   e[        X45        g s  snf s  snn
f s  snf s  snn
f )Nr         Fcopy	euclideanmetric   r@   r   r   nan_euclideanr?   rG         ?rG      	haversine)rG   rG   r   	manhattanr   )nprandomRandomStaterandom_sampleastyper#   r   r2   dtypenanr   tuplepir   r   shaper   )global_dtyperngXSS2YX_maskedY_maskedS_masked	S2_maskedrowvX_tuplesY_tupless                 e/var/www/html/dynamic-report/venv/lib/python3.13/site-packages/sklearn/metrics/tests/test_pairwise.py&test_pairwise_distances_for_dense_datarh   K   sS   
))


"C 	&!((E(BA1[1A	Q	BA77bhh.,..... 	&!((E(BA14A	Q	"BA77bhh.,.....   (//5/IH  (//5/IHVVHTNVVHTN!(_MH';IH(>>Y__<<<<<< ;#e,1Q,-;<H;#e,1Q,-;<H	H{	CBA77bhh.,.....
 	&!((E(BAAw}!BEE)A-AadGAw}!BEE)AadG1[1A	Q	BA 	&!((E(BAAw}!BEE)A-AadGAw}!BEE)AadG14A	Q	"BA
 	1[1A	Ai	0B771:###771:###A 	14A	A	3B771:###771:###A
 	11A	A	0B771:###771:###A_ -;,;s0   R
-R8R
R#
$R/R#
R
R#
z#array_namespace, device, dtype_name)idsrE   )r   rC   rN   c                 r   [        X5      n[        R                  R                  S5      nUR	                  S5      R                  USS9nUR	                  S5      R                  USS9nUR                  XaS9nUR                  XqS9n	[        SS9   [        XS9n
[        XS	9n[        U
5      S   R                  UR                  :X  d   eU
R                  UR                  :X  d   eU
R                  UR                  :X  d   e[        XcS9n[        X5        [        XUS
9n
[        XS	9n[        U
5      S   R                  UR                  :X  d   eU
R                  UR                  :X  d   eU
R                  UR                  :X  d   e[        XgUS
9n[        X5        S S S 5        g ! , (       d  f       g = f)Nr   r>   FrA   deviceTarray_api_dispatchrD   xpr^   rE   )r1   rO   rP   rQ   rR   rS   asarrayr   r#   r,   r.   __name__rl   rT   r2   )array_namespacerl   
dtype_namerE   rp   rZ   X_npY_npX_xpY_xpD_xpD_xp_npD_nps                rg   !test_pairwise_distances_array_apir}      s    
o	6B
))


"CV$++JU+CDV$++JU+CD::d:*D::d:*D	4	0!$6#D0T"1%.."++==={{dkk)))zzTZZ'''!$6& "$v>#D0T"1%.."++==={{dkk)))zzTZZ'''!$v>&' 
1	0	0s   DF((
F6coo_containercsc_containerbsr_containercsr_containerc                    [         R                  R                  S5      nUR                  S5      R	                  USS9nUR                  S5      R	                  USS9nU" U5      nU" U5      n	[        XSS9n
[        X5      n[        X5        U
R                  UR                  s=:X  a  U:X  d   e   e[        XSS9n
[        X5      n[        X5        U
R                  UR                  s=:X  a  U:X  d   e   e[        X" U5      S	S9n
[        U" U5      U " U5      5      n[        X5        U[         R                  :X  a'  U
R                  UR                  s=:X  a  U:X  d   e   eOJ[        R                  " [        5         U
R                  UR                  s=:X  a  U:X  d   e   e S S S 5        [        Xg5      n[        X5        U[         R                  :X  a'  U
R                  UR                  s=:X  a  U:X  d   e   eOJ[        R                  " [        5         U
R                  UR                  s=:X  a  U:X  d   e   e S S S 5        S
S0n[        Xg4SS0UD6n
[        Xg4S[        0UD6n[        X5        S
S0n[        U4SS0UD6n
[        U4S[        0UD6n[        X5        [        R                  " [         5         [        USS9  S S S 5        [        R                  " [         5         [        XiSS9  S S S 5        g ! , (       d  f       GNc= f! , (       d  f       N= f! , (       d  f       N_= f! , (       d  f       g = f)Nr   r>   FrA   rF   rC   rD   r   rN   p       @rE   r	   )rO   rP   rQ   rR   rS   r#   r   r2   rT   r   r   float64pytestraisesAssertionErrorr	   	TypeError)r~   r   r   r   rY   rZ   r[   r^   X_sparseY_sparser\   r]   kwdss                rg   'test_pairwise_distances_for_sparse_datar      s    ))


"C&!((E(BA&!((E(BA QHQH8kBA	X	0BA77bhh.,.....8h?A	(	-BA77bhh.,.....8]1%5kJA	]1-}Q/?	@BArzz!ww"((2l22222
 ]]>*77bhh6,66666 + 
Q	"BArzz!ww"((2l22222
 ]]>*77bhh6,66666 + :D1<<t<A	A	;	;d	;BA :D19[9D9A	A	8i	84	8BA 
y	!8K8 
"	y	!1{; 
"	!; +* +*  
"	!	!	!s0   'L'LL)2L:
L
L&)
L7:
Mcategoryc                 l   [         R                  R                  S5      nUR                  SS5      nUR	                  5       nSUS   -
  US'   [        [        S9   US 4 HA  n[        X$U S9n[        R                  " USS9  [         R                  " US:g  5      S:X  a  MA   e   S S S 5        S	U -  n[        R                  " [        US
9   [        X S9  S S S 5        [        R                  " [        US
9   [        UR                  [        5      X0S9  S S S 5        [        R                   " 5          [        R"                  " S[        5        [        UR                  [        5      U S9  S S S 5        g ! , (       d  f       N= f! , (       d  f       N= f! , (       d  f       N= f! , (       d  f       g = f)Nr   r?   r@   rL   rH   r   rD   )
fill_valuez+Data was converted to boolean for metric %smatchrq   error)rO   rP   rQ   randnrB   r5   r   r#   r/   
nan_to_numsumr   warnsrS   boolwarningscatch_warningssimplefilter)rE   rZ   r[   r^   Zresmsgs          rg   test_pairwise_boolean_distancer     sC    ))


"C		!QA	A!D'kAdG 
"7	8TA$Q&9CNN31-66#(#q(((  
9 8&
@C	+3	71, 
8 
+3	7188D>Q> 
8 
	 	 	"g'<=188D>&9 
#	"# 
9	8 
8	7 
8	7 
#	"s7   AE2 E2

F5F09F%2
F 
F
F"%
F3c                     [         R                  R                  S5      n U R                  SS5      n[        R
                  " 5          [        R                  " S[        5        [        USS9  S S S 5        g ! , (       d  f       g = f)Nr   r?   r@   r   r	   rD   )	rO   rP   rQ   r   r   r   r   r   r#   )rZ   r[   s     rg   test_no_data_conversion_warningr   (  sX    
))


"C		!QA		 	 	"g'<=1[1 
#	"	"s   &A66
Bfuncc                    [         R                  " [        SS9   U " [        R                  " S5      SS9  S S S 5        [         R                  " [        SS9   U " [        R                  " S5      [        R                  " S5      SS9  S S S 5        [         R                  " [        SS9   U " [        R                  " S5      [        R                  " S5      SS9  S S S 5        [        R                  " S5      nU " USS9nXL d   e[        R                  " S5      nU " U[        R                  " S	5      SS9nXL d   eU " [        R
                  " S
//SS9SS9nSUR                  R                  :X  d   eU " S//SS9n[        U[        R                  5      (       d   eg ! , (       d  f       GNq= f! , (       d  f       GN1= f! , (       d  f       N= f)Nz.* shape .*r   )r?      precomputedrD   )r@   r@   )r@   r   r?   r?   )r   r   rL   intrT   f      ?)
r   r   
ValueErrorrO   zerosarrayrT   kind
isinstancendarray)r   r\   r]   s      rg   test_pairwise_precomputedr   1  sZ    
z	7RXXfm4 
8 
z	7RXXfrxx/F 
8 
z	7RXXfrxx/F 
8 	A	a	&B7N7
A	a&!-	8B7N7 	RXXse5)-@A!'',, 	seW]+Aa$$$$1 
8	7 
8	7 
8	7s#   F 1F2)1G 
F/2
G
Gc                      [         R                  " [        SS9   [        [        R
                  " SS5      SS9  S S S 5        g ! , (       d  f       g = f)Nz.* non-negative values.*r   r   r   rD   )r   r   r   r#   rO   full     rg   &test_pairwise_precomputed_non_negativer   O  s3    	z)C	D27762.}E 
E	D	Ds    A
ArL   r?   doubleFrA   )wr   c           	         [        X5      u  p4[        [        R                  " U SUS9[        R                  " USUS940 UD6nUR                  S:  a0  UR                  US5      nUR                  S:X  a  [        US   5      $ [        S5      e)NrG   )ndimrp   r   )r   )rL   z6can only convert an array of size 1 to a Python scalar)	r.   r)   r/   
atleast_ndr   reshaperX   floatr   )xyr   rp   _KK_flats          rg   callable_rbf_kernelr   X  s    !EBqqR(#..r*J	NR	A
 	vvzAu%<<4##
M
NNr   zfunc, metric, kwdsrC   r	   
polynomialdegreegamma皙?rT   c                 j   [         R                  R                  S5      n[         R                  " SUR	                  S5      -  US9n[         R                  " SUR	                  S5      -  US9nU " U4USS.UD6nU " U4USS.UD6n[        Xx5        U " XV4USS.UD6nU " XV4USS.UD6n[        Xx5        g )	Nr   r?   r>   r   r   r@   rL   rE   n_jobsrG   )rO   rP   rQ   r   rR   r2   )	r   rE   r   rT   rZ   r[   r^   r\   r]   s	            rg   test_pairwise_parallelr   g  s    ( ))


"C
S&&v..e<A
S&&v..e<AQ0va040A	a	1q	1D	1BAQ3&3d3A	a	46!	4t	4BAr   rN   	laplacianc                    [        X45      n[        R                  R                  S5      n[        R                  " SUR                  S5      -  US9n[        R                  " SUR                  S5      -  US9n	UR                  XS9n
UR                  XS9n[        SS9   S	 GH   nUc  S OUnUc  S OU	n	U " X4US
S.UD6n[        XS9n[        U5      S   R                  UR                  :X  d   eUR                  U
R                  :X  d   eUR                  U
R                  :X  d   eU " X4USS.UD6n[        XS9n[        U5      S   R                  UR                  :X  d   eUR                  U
R                  :X  d   eUR                  U
R                  :X  d   eU " U4USS.UD6n[        UU5        [        UU5        GM#     S S S 5        g ! , (       d  f       g = f)Nr   r?   r>   r   r   rk   Trm   )Nznot nonerL   r   ro   rG   )r1   rO   rP   rQ   r   rR   rr   r   r,   r.   rs   rl   rT   r2   )r   rE   r   rt   rl   ru   rp   rZ   rv   rw   rx   ry   y_val	n_job1_xpn_job1_xp_np	n_job2_xpn_job2_xp_np	n_job2_nps                     rg    test_pairwise_parallel_array_apir     s   $ 
o	6B
))


"C88A))&11DD88A))&11DD::d:*D::d:*D	4	0'E =4dD =4dDTIqIDII,Y>L +A.772;;FFF##t{{222??djj000TIqIDII,Y>L +A.772;;FFF##t{{222??djj000TC&CdCIL,7L)4' ( 
1	0	0s    D(G
Gc                  0    [        S//S S9S   S:X  d   eg )Nr   c                     g)Nr?   r   r   r   s     rg   <lambda>9test_pairwise_callable_nonstrict_metric.<locals>.<lambda>  s    1r   rD   rH   r?   )r#   r   r   rg   'test_pairwise_callable_nonstrict_metricr     s"     ugn=dCqHHHr   )rbfr   sigmoidr   linearchi2additive_chi2c                    [         R                  R                  S5      nUR                  S5      nUR                  S5      n[        U    n[        X0S9nU" U5      n[        Xg5        [        X4U S9nU" X4S9n[        Xg5        [        U VV	s/ s H  n[        U V	s/ s H  oPM     sn	5      PM      sn	n5      n
[        U VV	s/ s H  n[        U V	s/ s H  oPM     sn	5      PM      sn	n5      n[        XU S9n[        Xg5        U" U5      nU" U5      nU S;   a  g [        XU S9n[        Xg5        g s  sn	f s  sn	nf s  sn	f s  sn	nf )Nr   r>   rF   rD   rq   r^   )r   r   )rO   rP   rQ   rR   r   r'   r2   rV   )rE   r   rZ   r[   r^   functionK1K2rc   rd   re   rf   r   r   s                 rg   test_pairwise_kernelsr     s-    ))


"C&!A&!A(0H	!	+B	!BB	!	0B	!	BB;#e,1Q,-;<H;#e,1Q,-;<H	(V	<BB QHQH**	(v	>BB -;,;s0   D7
D2&D7
E
D=E
2D7
=E
)r   r   r   r   r   r   r   c                    [        X5      n[        R                  R                  S5      nSUR	                  S5      -  nUR                  USS9nSUR	                  S5      -  nUR                  USS9nUR                  XbS9nUR                  XrS9n	[        SS	9   [        XS
9n
[        XS9n[        U
5      S   R                  UR                  :X  d   eU
R                  UR                  :X  d   eU
R                  UR                  :X  d   e[        X`S
9n[        X5        [        XU S9n
[        XS9n[        U
5      S   R                  UR                  :X  d   eU
R                  UR                  :X  d   eU
R                  UR                  :X  d   e[        XgU S9n[        X5        S S S 5        g ! , (       d  f       g = f)Nr   
   r>   FrA   rF   rk   Trm   rD   ro   rq   )r1   rO   rP   rQ   rR   rS   rr   r   r'   r,   r.   rs   rl   rT   r2   )rE   rt   rl   ru   rp   rZ   rv   rw   rx   ry   K_xpK_xp_npK_nps                rg   test_pairwise_kernels_array_apir     s    
o	6B
))


"C!!&))D;;z;.D!!&))D;;z;.D::d:*D::d:*D	4	04#D0T"1%.."++==={{dkk)))zzTZZ'''4&  V<#D0T"1%.."++==={{dkk)))zzTZZ'''V<&' 
1	0	0s   DF22
G c                  6   [         R                  R                  S5      n U R                  S5      nU R                  S5      n[        nSS0n[        U4X#S.UD6n[        U4SU0UD6n[        XV5        [        U4XS.UD6n[        U4SU0UD6n[        XV5        g )Nr   r>   rF   r   r   rq   r^   )rO   rP   rQ   rR   r   r'   r)   r2   )rZ   r[   r^   rE   r   r   r   s          rg   test_pairwise_kernels_callabler     s     ))


"C&!A&!A FS>D	!	8q	84	8B	A	#	#d	#BB 
!	8q	84	8B	A	#	#d	#BBr   c                  \   [         R                  R                  S5      n U R                  S5      nU R                  S5      n[	        XSS9nSSS.n[        X4SS	S
.UD6n[        X55        [        R                  " [        5         [        X4SS0UD6  S S S 5        g ! , (       d  f       g = f)Nr   r>   rF   r   r   z:))r   blablar   T)rE   filter_paramsrE   )
rO   rP   rQ   rR   r)   r'   r2   r   r   r   )rZ   r[   r^   r   paramsr   s         rg   "test_pairwise_kernels_filter_paramr   "  s    
))


"C&!A&!A1s#Ad+F	!	KuD	KF	KBA	y	!6e6v6 
"	!	!s   B
B+zmetric, funcc                 v   [         R                  R                  S5      nUR                  S5      nUR                  S5      n[	        XEU S9nU" XE5      n[        Xg5        U" U" U5      U" U5      5      n[        Xh5        U [        ;   a1  [        U    " XE5      n	[         R                  " U	5      n	[        X5        g g )Nr   r>   rD   )rO   rP   rQ   rR   r    r2   r   diag)
rE   r   r   rZ   r[   r^   r\   r]   S3	distancess
             rg   test_paired_distancesr   /  s     ))


"C&!A&!Af-A	aBA	mAa 0	1BA,, 07=	GGI&		% -r   c                    [         R                  R                  S5      nUR                  S5      R	                  U SS9nUR                  S5      R	                  U SS9n[        X#SS9n[        X#S S9n[        XE5        UR                  S5      n[        R                  " [        5         [        X#5        S S S 5        g ! , (       d  f       g = f)	Nr   r>   FrA   rN   rD   c                 L    [         R                  " X-
  5      R                  SS9$ )Nr   axis)rO   absr   r   s     rg   r   0test_paired_distances_callable.<locals>.<lambda>P  s    BFF15M4E4E14E4Mr   r   )
rO   rP   rQ   rR   rS   r    r2   r   r   r   )rY   rZ   r[   r^   r\   r]   s         rg   test_paired_distances_callabler  F  s     ))


"C&!((E(BA&!((E(BAk2A	!'M	NBA 	&!A	z	" 
#	"	"s   )B>>
Cdok_containerc                 l   [         R                  " S/S//US9n[         R                  " S/S//US9nU " U5      nU" XBS9nSS/nSS/nSS/n	[        X4SS	9u  p[        X4SS	9n[	        X5        [	        X5        [	        X5        [        XVSS	9u  p[        XVSS	9n[	        X5        [	        X5        [	        X5        [        U5      [         R                  :X  d   e[        U5      [         R                  :X  d   e[        X4S
S	9u  p[        X4SSS0S9u  nn[        X4S
S	9n[        X4SSS0S9n[	        X5        [	        UU	5        [	        X5        [	        X5        [	        UU5        [	        UU5        [        X4SS	9u  p[        X4SS	9n[	        X5        [	        X5        [	        X5        [        XVSS	9u  p[        XVSS	9n[	        X5        [	        X5        [	        X5        [        X4[        SS0S9u  p[	        X5        [	        X5        [        X4SSS0S9u  p[	        X5        [	        X5        [         R                  R                  S5      nUR                  SS5      nUR                  SS5      n[        X4SS	9nUR                  SS9nUU[        [        U5      5      4   n[        X4SSS9u  nn[	        UUSS9  [	        UUSS9  [        X4SS9u  nn[        XCSS9u  nn[	        UU5        [        UU5        [        X3SS9u  nn[        X3SS9u  nn[	        UU5        [        UU5        [        X4SS9n[        XCSS9n[        UU5        [        X3SS9n[        X3SS9n[        UU5        [        X45      n[        [         R                   " U5      [         R                   " U5      5      n[        UU5        g )Nr   rL   r   r   rG   r@   rC   rD   sqeuclideansquaredT)rE   metric_kwargsrN   r   r	   a      o   r   )r   rE   Hz>rtol)rO   rr   r%   r$   r2   typer   r	   rP   rQ   r   r#   argminrangelenr4   asfortranarray)r  r   rY   r[   r^   XspYspexpected_idxexpected_valsexpected_vals_sqidxvalsidx2idxspvalsspidxsp2vals2idx3idx4rZ   distdist_orig_inddist_orig_valdist_chunked_inddist_chunked_valargmin_0dist_0argmin_1dist_1argmin_C_contiguousargmin_F_contiguouss                                  rg   "test_pairwise_distances_argmin_minr-  \  s   
 	

QC!:\2A


RD1#;l3A

C

.Cq6LFM1v .a;GIC$Q+>DC&D'D(1#;OME&sDFE(F)F*;"**$$$<2::%%% .a=IIC/	[D0AKD% %Q-@D$	[D0AD D+E+,C&D'D,'D,' .a;GIC$Q+>DC&D'D(1#;OME&sDFE(F)F* .	YsAhIC C&D( .	[aIC C&D( ))


"C		"cA		#sAa;7DKKQK'Mc-.@(AABM)F	1[*&& M#3$?M#3$? 5QBHf4QBHfFF#x*4QBHf4QBHfFF#x* )A6H(A6Hx*(A6H(A6Hx* 4A93
!b//2 *,?@r   c                     U S S 2S S24   $ )Nd   r   r"  starts     rg   _reduce_funcr2    s    4C4=r   c                    [         R                  R                  S5      nUR                  S5      R	                  U SS9n[        U5      S S 2S S24   n[        US [        SS9n[        U[        5      (       d   e[        U5      n[        U5      S:  d   eUS   R                  UR                  :X  d   e[        [         R                  " U5      US	S
9  g )Nr   )i  r@   FrA   r/        >reduce_funcworking_memoryrL   r  atol)rO   rP   rQ   rR   rS   r#   r&   r2  r   r   listr  rT   r2   vstack)rY   rZ   r[   r\   S_chunkss        rg   &test_pairwise_distances_chunked_reducer=    s    
))


"C(#**<e*DA1a#g&A)	4\&H h....H~Hx=1A;''' BIIh'6r   c                 2   [         R                  R                  S5      nUR                  S5      R	                  U SS9n[        US S SS9n[        U[        5      (       d   e[        U5      n[        U5      S:  d   e[        S	 U 5       5      (       d   eg )
Nr   r   r@   FrA   c                     g Nr   r0  s     rg   r   =test_pairwise_distances_chunked_reduce_none.<locals>.<lambda>  s    r   r4  r5  rL   c              3   (   #    U  H  oS L v   M
     g 7frA  r   ).0chunks     rg   	<genexpr>>test_pairwise_distances_chunked_reduce_none.<locals>.<genexpr>  s     3(}(s   )rO   rP   rQ   rR   rS   r&   r   r   r:  r  all)rY   rZ   r[   r<  s       rg   +test_pairwise_distances_chunked_reduce_nonerI    s    
))


"C'")),U)CA)	45fH h....H~Hx=13(33333r   good_reducec                     [        U 5      $ rA  r:  Dr1  s     rg   r   r     s    ar   c                 .    [         R                  " U 5      $ rA  )rO   r   rM  s     rg   r   r     s    !r   c                 .    [        U 5      [        U 5      4$ rA  rL  rM  s     rg   r   r     s    $q'47+r   c                     U" U 5      $ rA  r   )rN  r1  scipy_csr_types      rg   r   r     s	    q8Ir   c                 R    U" U 5      [         R                  " U 5      [        U 5      4$ rA  )rO   r   r:  )rN  r1  scipy_dok_types      rg   r   r     s!    1HHQKG9
r   c                 ~    [         R                  " S5      R                  SS5      n[        US U SS9n[	        U5        g )Nr   r   rL   @   r5  )rO   aranger   r&   next)rJ  r[   r<  s      rg   ,test_pairwise_distances_chunked_reduce_validrY    s:    * 			"b!$A)	4[H 	Nr   )
bad_reduceerr_typemessagec                 6    [         R                  " X SS  /5      $ Nr   rO   concatenaterN  ss     rg   r   r     s    bcF4r   zlength 11\..* input: 10\.c                 :    U [         R                  " X SS  /5      4$ r^  r_  ra  s     rg   r   r     s    !R^^Q"#K89r   z!length \(10, 11\)\..* input: 10\.c                     U S S U 4$ )N	   r   ra  s     rg   r   r   #  s    q!uajr   z length \(9, 10\)\..* input: 10\.c                     g)N   r   ra  s     rg   r   r   %  s    r   z2returned 7\. Expected sequence\(s\) of length 10\.c                     g)N)rg     r   ra  s     rg   r   r   *  s    r   z9returned \(7, 8\)\. Expected sequence\(s\) of length 10\.c                 2    [         R                  " S5      S4$ )Nr   re  )rO   rW  ra  s     rg   r   r   /  s    "))B-+r   z-, 9\)\. Expected sequence\(s\) of length 10\.c                     [         R                  " S5      R                  SS5      R                  U SS9n[	        US USS9n[
        R                  " X#S9   [        U5        S S S 5        g ! , (       d  f       g = f)	Nr   r   rL   FrA   rV  r5  r   )rO   rW  r   rS   r&   r   r   rX  )rY   rZ  r[  r\  r[   r<  s         rg   .test_pairwise_distances_chunked_reduce_invalidrl    sc    D 			"b!$++Lu+EA)	4ZH 
x	/X 
0	/	/s   A**
A8c                 2   [        XX#S9n[        U[        5      (       d   e[        U5      nUc  U OUn[	        U5      S-  S-  nU H#  nUR
                  nU[        X&5      S-  ::  a  M#   e   [        R                  " U5      n[        XUS9n	[        XYSS9  g )Nr7  rE   ri  g      >i   rD   r  r8  )r&   r   r   r:  r  nbytesmaxrO   r;  r#   r2   )
r[   r^   r7  rE   genblockwise_distancesmin_block_mibblockmemory_usedr\   s
             rg    check_pairwise_distances_chunkedrv  @  s    
$Q.
XCc=))))s)YAAFQJ'M$llc.@5HHHH % ))$781/A'6r   )rC   l2r  c                 .   [         R                  R                  S5      nUR                  SSS9R	                  USS9n[        [        USU S95      n[        U5      S:  d   e[        [         R                  " [         R                  " U5      5      SS	S
9  g )Nr     r       _BsizescaleFrA   rL   rn  绽|=r  )rO   rP   rQ   normalrS   r:  r&   r  r2   r   r;  )rE   rY   rZ   r[   chunkss        rg   (test_pairwise_distances_chunked_diagonalr  P  sz    
))


"C


$
/66|%6PA,QqPQFv;??BGGBIIf-.>r   c                     [         R                  R                  S5      nUR                  SSS9R	                  USS9n[        X0SS9n[        [         R                  " U5      SS	S
9  g )Nr   ry  r{  r|  FrA   rG   r   r  r8  )rO   rP   rQ   r  rS   r#   r2   r   )rE   rY   rZ   r[   r   s        rg   )test_parallel_pairwise_distances_diagonalr  Y  sY    
))


"C


$
/66|%6PA"1A>IBGGI&6r   z0ignore:Could not adhere to working_memory configc                    [         R                  R                  S5      nUR                  S5      R	                  U SS9n[        US SSS9  [        SS5       H  n[        US S	U-  SS9  M     [        UR                  5       S SSS9  UR                  S
5      R	                  U SS9n[        X$SSS9  [        UR                  5       UR                  5       SSS9  [        X$SSS9  [        X$SSS9  [        U5      n[        USSS9n[        U[        5      (       d   e[        U5      UL d   e[        R                  " [        5         [        U5        S S S 5        g ! , (       d  f       g = f)Nr   )   r@   FrA   rL   rC   rn  irG   )r/  r@   i'  r   r4  r   )rO   rP   rQ   rR   rS   rv  r  tolistr#   r&   r   r   rX  r   r   StopIteration)rY   rZ   r[   powerr^   rN  rq  s          rg   test_pairwise_distances_chunkedr  a  sG    ))


"C(#**<e*DA$QQ{SsA(tAuH[	
 
 %	
D; 	(#**<e*DA$Q!KP$	
AHHJq %Q%T %Q!KP 	1A
$Qvm
TCc=))))9>>	}	%S	 
&	%	%s   E
E)x_array_constrdensey_array_constrc                 d    U " S//5      nU" S/S//5      n[        X#5      n[        USS//5        g Nr   rL   rG   r   r   )r   r2   )r  r  r[   r^   rN  s        rg   %test_euclidean_distances_known_resultr    s>     	uAaSz"AA!AAc
|$r   c                 D   [         R                  R                  S5      nUR                  S5      R	                  U SS9nUR                  S5      R	                  U SS9nUR	                  [         R
                  5      S-  R                  SS9R                  SS	5      nUR	                  [         R
                  5      S-  R                  SS9R                  SS	5      nU" U5      n[        X45      n[        X4US
9n[        X4US9n	[        X4XVS9n
[        X5        [        X5        [        X5        [        UU[         R                  " U5      [         R                  " U5      S9n[        R                  " [        5         [        X5        S S S 5        g ! , (       d  f       g = f)Nr   r   r   FrA      r   rG   rL   r   r   X_norm_squaredY_norm_squaredr  r  )rO   rP   rQ   rR   rS   r   r   r   r   r2   
zeros_liker   r   r   )rY   r  rZ   r[   r^   	X_norm_sq	Y_norm_sqD1D2D3D4wrong_Ds               rg   #test_euclidean_distances_with_normsr    sZ    ))


"C(#**<e*DA(#**<e*DA "**%*//Q/7??2FI"**%*//Q/7??2FIqA	Q	"B	Q)	<B	Q)	<B	Q)	VBBBB "		}}Y/}}Y/	G 
~	&$ 
'	&	&s   <F
F	symmetricTc                 (   [         R                  R                  U 5      nUR                  S5      nU(       a  UOUR                  S5      nUR	                  [         R
                  5      S-  R                  SS9R                  SS5      nUR	                  [         R
                  5      S-  R                  SS9R                  SS5      n[        X45      n[        X4US9n[        X4US9n	[        X4XVS	9n
[        X5        [        X5        [        X5        g )
Nr  r  rG   rL   r   r   r  r  r  )
rO   rP   rQ   rR   rS   float32r   r   r   r2   )global_random_seedr  rZ   r[   r^   r  r  r  r  r  r  s              rg   &test_euclidean_distances_float32_normsr    s     ))

 2
3C(#AC--h7A"**%*//Q/7??2FI"**%*//Q/7??2FI	Q	"B	Q)	<B	Q)	<B	Q)	VBBBBr   c            	         [         R                  R                  S5      n U R                  S5      nU R                  S5      nUS-  R	                  SS9nUS-  R	                  SS9n[        XX4S9n[        UUUR                  SS5      UR                  SS5      S9n[        UUUR                  SS5      UR                  SS5      S9n[        Xe5        [        Xu5        [        R                  " [        S	S
9   [        XUS S S9  S S S 5        [        R                  " [        SS
9   [        XUS S S9  S S S 5        g ! , (       d  f       N>= f! , (       d  f       g = f)Nr   r  r  rG   rL   r   r  r   zIncompatible dimensions for Xr   r?   r  zIncompatible dimensions for Yr  )rO   rP   rQ   rR   r   r   r   r2   r   r   r   )rZ   r[   r^   r  r  r  r  r  s           rg   $test_euclidean_distances_norm_shapesr    sC   
))


"C(#A(#AdZZQZ'NdZZQZ'N		^
B 
		%--b!4%--b!4	
B 
		%--a4%--a4	
B BB	z)H	IA1CD 
J	z)H	IA1CD 
J	I 
J	I	I	Is   7D=&E=
E
Ec                 ^   [         R                  R                  S5      nUR                  S5      R	                  U SS9nSXDS:  '   UR                  S5      R	                  U SS9nSXUS:  '   [        XE5      nU" U5      nU" U5      n[        XE5      n[        XvSS9  UR                  U :X  d   eg )	Nr   r/  r   FrA   皙?r  ư>r  )	rO   rP   rQ   rR   rS   r   r   r2   rT   )rY   r  r  rZ   r[   r^   expectedr   s           rg   test_euclidean_distancesr    s     ))


"C)$++Lu+EAA#gJ(#**<e*DAA#gJQ{HqAqA#A)I Id3??l***r   c                    [         R                  R                  S5      nUR                  S5      R	                  U SS9nSX3S:  '   [        [        U5      5      nU" U5      n[        U5      n[        XTSS9  UR                  U :X  d   eg )Nr   r  FrA   r  r  r  )
rO   rP   rQ   rR   rS   r   r
   r   r2   rT   )rY   r  rZ   r[   r  r   s         rg   test_euclidean_distances_symr    s     ))


"C)$++Lu+EAA#gJ%(#HqA#A&I Id3??l***r   
batch_size)Nr?   rg  e   c                    [         R                  R                  S5      nUR                  S5      R	                  [         R
                  5      nSXDS:  '   UR                  S5      R	                  [         R
                  5      nSXUS:  '   [        XE5      nU" U5      nU" U5      n[        XEU S9n[         R                  " [         R                  " US5      5      n[        XvSS9  g )Nr   r  r  r  r^   r  r  r  )rO   rP   rQ   rR   rS   r  r   r   sqrtmaximumr2   )r  r  r  rZ   r[   r^   r  r   s           rg   test_euclidean_distances_upcastr  &  s     ))


"C)$++BJJ7AA#gJ(#**2::6AA#gJQ{HqAqA+AzJI

9a01I Id3r   c                 b   [         R                  R                  S5      nUR                  S5      R	                  [         R
                  5      nSX3S:  '   [        [        U5      5      nU" U5      n[        X3U S9n[         R                  " [         R                  " US5      5      n[        XTSS9  g )Nr   r  r  r  r  r  )rO   rP   rQ   rR   rS   r  r   r
   r   r  r  r2   )r  r  rZ   r[   r  r   s         rg   #test_euclidean_distances_upcast_symr  E  s     ))


"C)$++BJJ7AA#gJ%(#HqA+AzJI

9a01I Id3r   zdtype, eps, rtolg-C6?h㈵>g:0yE>gGz?z failing due to lack of precision)reason)marksdimi@B c                     [         R                  " S/U-  /U S9n[         R                  " SU-   /U-  /U S9n[        XE5      n[        XE5      n[	        XgSS9  g )Nr   r   r  r  )rO   r   r   r   r2   )rT   epsr  r  r[   r^   r   r  s           rg   'test_euclidean_distances_extreme_valuesr  \  sW      	3%#+e,A
39+#$E2A#A)IQ{HId3r   r  c                     [         R                  R                  S5      nUR                  SS5      nUR                  SS5      n[	        X#U S9n[        X#U S9n[        XE5        g )N9  r   r@   )r^   r  )rO   rP   rQ   r   r   r   r2   )r  rZ   r[   r^   normal_distancenan_distances         rg   8test_nan_euclidean_distances_equal_to_euclidean_distancer  u  sV     ))


%C		!QA		!QA)!'BO*17CLO2r   r[   r^   c                     [         R                  " [        5       n[        XS9  S S S 5        SnU[	        WR
                  5      :X  d   eg ! , (       d  f       N,= f)Nr   zBInput contains infinity or a value too large for dtype('float64').)r   r   r   r   strvalue)r[   r^   excinfoexp_msgs       rg   ,test_nan_euclidean_distances_infinite_valuesr    sG     
z	"g' 
# SGc'--((((	 
#	"s   
A
AzX, X_diag, missing_valuerG   r   c                    [         R                  " SU/US//5      n[        XS9n[        X45        [        U SUS9n[        US-  U5        [        X US9n[        X65        [        X R	                  5       US9n[        X75        g )N        r   missing_valuesTr  r  rG   )rO   r   r   r2   rB   )r[   X_diagmissing_valueexp_distr"  dist_sqdist_twodist_two_copys           rg    test_nan_euclidean_distances_2x2r    s{      xx#v45H"1CDH#%amTGHaK)&qMJHH'+AvvxVMH,r   r  c                 4   [         R                  " X /SS//5      n[         R                  " [         R                  [         R                  /[         R                  S//5      n[        XS9n[	        X#5        [        XR                  5       U S9n[	        X#5        g )Nr   rL   r  )rO   r   rU   r   r2   rB   )r  r[   r  r"  s       rg   )test_nan_euclidean_distances_complete_nanr    sp    
=01a&9:Axx"&&"&&)BFFA;78H"1CDH#"1ffh}MDH#r   c           	      P   [         R                  " SU SSS/U SSSU /SX U S//5      n[         R                  " U SSU S/X SSS/X U SS//5      n[        XU S9n[        X!U S9n[        X4R                  5        [        [        US S	 US S	 S
U S9S//5        [        [        US	S US	S SU S9[         R                  " S5      //5        [        XS9n[        XU S9n[        XR                  5       U S9n[        XV5        [        Xg5        [        XS
S9n[        XSS9n	[        X5        g )Nr   g      @g      @r   g      @g      @g      @r  rL   Tr  g      D@rG   Fg      9@rA   )rO   r   r   r3   Tr2   r  rB   )
r  r[   r^   r  r  r  r  D5D6D7s
             rg   'test_nan_euclidean_distances_not_trivalr    sh   
-c3/Cc=9-sC	
	A 	CmS93S9=#sC	
	A 
!m	DB	 m	DBDD! bqE1Ra5$}	
 4	45	 aFAaFE-	
 '';
<	=>	 
!	AB	 m	DB	 FFH]	KBBB 
!D	1B	 E	2BBr   c                     [         R                  " SSU S/SU SU //5      n[        XSS9n[         R                  " US:  5      (       d   e[        XSS9n[	        US	5        g )
NgzG^g     @@gB@g|ï@T)r  r  r   Fr  )rO   r   r   rH  r2   )r  r[   dist_squaredr"  s       rg   7test_nan_euclidean_distances_one_feature_match_positiver    sr    
 	e]E2mZ?	
	A +	L 66,!#$$$$"1ERDD#r   c                     [         R                  R                  S5      n [         R                  " U R	                  S5      5      n[         R
                  " X/5      n[        U5      n[        USS/SS//SS9  [         R                  " US:  5      (       d   e[         R                  " US:*  5      (       d   e[        U[         R                  " U5         SS/5        [         R
                  " X* /5      n[        U5      n[         R                  " US:  5      (       d   e[         R                  " US:*  5      (       d   e[        USS/SS//5        [         R                  " U R	                  SS5      5      n[        U5      n[        U[         R                  " U5         S/UR                  S	   -  5        [         R                  " US:  5      (       d   e[         R                  " US:*  5      (       d   eg )
Nr  i  r  r  r8  r   rz  i  r   )rO   rP   rQ   r   randr;  r   r2   rH  diag_indices_fromrX   )rZ   r   XArN  XBr  r[   s          rg   test_cosine_distancesr    s   
))


%C
sxx}A	A6	BAAc
S#J/e<66!s(66!s(Ab**1-.c
;	Ar7	B	"	B66")66")B#sc3Z01 	sxxd#$AAAb**1-.
0BC66!s(66!s(r   c                     S n [         R                  R                  S5      nUR                  S5      nUR                  S5      n[         R                  " U VVs/ s H  oC Vs/ s H
  oP" XE5      PM     snPM     snn5      n[        X#5      n[        Xg5        UR                  S5      nSn[        R                  " [        US9   [        U5        S S S 5        g s  snf s  snnf ! , (       d  f       g = f)Nc                 f   US   U S   -
  nUS   U S   -
  n[         R                  " US-  5      S-  [         R                  " U S   5      [         R                  " US   5      -  [         R                  " US-  5      S-  -  -   nS[         R                  " [         R                  " U5      5      -  nU$ )Nr   rL   rG   )rO   sincosarcsinr  )r   r   diff_latdiff_lonacs         rg   slow_haversine_distances:test_haversine_distances.<locals>.slow_haversine_distances!  s    Q4!A$;Q4!A$;FF8a< A%FF1Q4L266!A$<'"&&A*>!*CC
 		"''!*%%r   r   rJ   )r   rG   )r   r   z-Haversine distance only valid in 2 dimensionsr   )
rO   rP   rQ   rR   r   r   r2   r   r   r   )	r  rZ   r[   r^   r   r   r  r  err_msgs	            rg   test_haversine_distancesr    s     ))


"C&!A'"A	J1A>Aq,Q2A>J	KB	Q	"BB'"A=G	z	1A 
2	1 ?J 
2	1s$   	C
#C4C
C%C
%
C3c                  N    S/S//n S/S//n[        X5      n[        USS/5        g r  )r!   r2   r[   r^   rN  s      rg   test_paired_euclidean_distancesr  :  4    
qc
A
qc
A"1(AASz"r   c                  N    S/S//n S/S//n[        X5      n[        USS/5        g r  )r"   r2   r  s      rg   test_paired_manhattan_distancesr  B  r  r   c                  N    S/S//n S/S//n[        X5      n[        USS/5        g )Nr   rL   rG   rK   )r   r2   r  s      rg   test_paired_cosine_distancesr  J  s4    
qc
A
qc
A%AASz"r   c                     [         R                  R                  S5      n U R                  S5      nU R                  S5      n[	        X5      nSn[        XUS9nUR                  [        :X  d   e[        U5       Hr  u  pg[        U5       H^  u  p[         R                  " Xy-
  S-  Xy-   -  5      * n
[         R                  " XJ-  5      n[        X6U4   U
5        [        XVU4   U5        M`     Mt     [        U5      n[        [         R                  " U5      S5        [         R                  " US:  5      (       d   e[         R                  " U[         R                  " [         R                  " U5      5      -
  S:  5      (       d   eU R                  S5      R                  [         R                   5      nU R                  S5      R                  [         R                   5      n[        X5      nUR                  [         R                   :X  d   eU R                  S5      R                  [         R"                  5      n[        X5      n[         R$                  " U5      R                  5       (       d   eUR                  [        :X  d   eSS	/S
S//nSS/SS//n[        X5      nUS   US   :  d   eUS   US   :  d   e[&        R(                  " [*        5         [        SS//5        S S S 5        [&        R(                  " [*        5         [        SS//SS//5        S S S 5        [&        R(                  " [*        5         [        SS//SS//5        S S S 5        [&        R(                  " [*        5         [        SS/// SQ/5        S S S 5        g ! , (       d  f       N= f! , (       d  f       N= f! , (       d  f       Nf= f! , (       d  f       g = f)Nr   r>   r?  r   r   rG   rL   g333333?gffffff?r   g?rH   )r   rL   )rL   rL   )rL   r   r   )皙?r   g333333?)rO   rP   rQ   rR   r   r   rT   r   	enumerater   expr3   r4   r   rH  rS   r  int32isfiniter   r   r   )rZ   r[   r^   K_addr   r   ir   jr   r   chi2_exps               rg   test_chi_square_kernelr	  R  s    
))


"C&!A'"A &EEA&A77e!aLDAFFAEa<15122Dvvel+HdT2Q$2	 !  	AArwwqz1%66!a%===66!bggbggaj))A-....&!((4A'"))"**5AAA77bjj    	'"))"((3AAA;;q>77e sc1XA
Q#sAAAT7QtWT7QtW 
z	"aWI 
#	z	"aWIRz* 
#	z	"aVHBxj) 
# 
z	"aVH/0 
#	" 
#	"	"	"	"	" 
#	"s0   N	N >N13O
N 
N.1
N?
Okernelc                     [         R                  R                  S5      nUR                  S5      nU " X"5      n[	        X3R
                  S5        g )Nr   r>      )rO   rP   rQ   rR   r2   r  )r
  rZ   r[   r   s       rg   test_kernel_symmetryr    s?     ))


"C&!AqAAssBr   c                     [         R                  R                  S5      nUR                  S5      nU" U5      nU " X35      nU " XD5      n[	        XV5        g Nr   r>   )rO   rP   rQ   rR   r2   )r
  r   rZ   r[   r   r   r   s          rg   test_kernel_sparser    sM     ))


"C&!AQHqA		#BAr   c            	         [         R                  R                  S5      n U R                  S5      n[	        X5      n[        UR                  S S S2   U Vs/ s H  n[        R                  " U5      S-  PM     sn5        g s  snf )Nr   r>      rG   )	rO   rP   rQ   rR   r   r2   flatr   norm)rZ   r[   r   r   s       rg   test_linear_kernelr    se    
))


"C&!AaAAFF3Q3Kq!Aq!&++a.A"5q!AB!As   #B
c                      [         R                  R                  S5      n U R                  S5      n[	        X5      n[        UR                  S S S2   [         R                  " S5      5        g )Nr   r>   r  r?   )rO   rP   rQ   rR   r)   r2   r  onesrZ   r[   r   s      rg   test_rbf_kernelr    sN    
))


"C&!A1AAFF3Q3K,r   c                     [         R                  R                  S5      n U R                  S5      n[	        X5      n[        [         R                  " U5      [         R                  " S5      5        [         R                  " US:  5      (       d   e[         R                  " U[         R                  " [         R                  " U5      5      -
  S:  5      (       d   eg )Nr   r>   r?   rL   )	rO   rP   rQ   rR   r   r2   r   r  rH  r  s      rg   test_laplacian_kernelr    s    
))


"C&!AABGGAJ
+ 66!a%===66!bggbggaj))A-....r   zmetric, pairwise_funcr   r   c                    [         R                  R                  S5      nUR                  S5      nUR                  S5      nU" U5      nU" U5      nU" XgSS9n[	        U5      (       d   eU" XESS9n	[	        U	5      (       a   e[        UR                  5       U	5        [        XEU S9n
[        UR                  5       U
5        g )Nr   r>   r   F)dense_outputTrq   )rO   rP   rQ   rR   r   r2   toarrayr'   )rE   pairwise_funcr   rZ   r[   r^   XcsrYcsrr   r   K3s              rg   &test_pairwise_similarity_sparse_outputr#    s     ))


"C&!A&!ADD 
t	6BB<<< 
q$	/B||BJJL"% 
!	0BBJJL"%r   c                 @   [         R                  R                  S5      nUR                  S5      nUR                  S5      nU " U5      nU " U5      nUS 4X#4US 4XE44 H=  u  pg[	        XgSS9n[        U5      nUb  [        U5      n[	        XgSS9n	[        X5        M?     g )Nr   r>   r   r   rq   r   )rO   rP   rQ   rR   r'   r+   r2   )
r   rZ   r[   r^   r   r!  X_Y_r   r   s
             rg   test_cosine_similarityr'    s     ))


"C&!A&!ADDt9qftTlTLA bx8r]>2Bbx8 Br   c                      [         R                  " [         R                  " S5      S5      n [        U S 5      u  pXL d   e[	        X5        g N(   r?   ri  )rO   resizerW  r   r4   )r  
XA_checked
XB_checkeds      rg   test_check_dense_matricesr/    s?     
299R=&	)B22t<J###r&r   c                     [         R                  " [         R                  " S5      S5      n [         R                  " [         R                  " S5      S5      n[        X5      u  p#[	        X5        [	        X5        [         R                  " [         R                  " S5      S5      n[        X5      u  p#[	        X5        [	        X5        g )Nr*  r+      r@   ri  )rO   r,  rW  r   r4   r   r  r  r-  r.  s       rg   test_check_XB_returnedr4    s     
299R=&	)B	299R=&	)B22:Jr&r&	299R=&	)B08Jr&r&r   c                     [         R                  " [         R                  " S5      S5      n [         R                  " [         R                  " S5      S5      n[        R                  " [
        5         [        X5        S S S 5        [         R                  " [         R                  " S5      S5      n[        R                  " [
        5         [        X5        S S S 5        g ! , (       d  f       Nh= f! , (       d  f       g = f)N-   )r?   re  r1  r2  $   )r@   re  )rO   r,  rW  r   r   r   r   r   r  r  s     rg   test_check_different_dimensionsr9    s    	299R=&	)B	299R=&	)B	z	"b% 
# 
299U#V	,B	z	"B# 
#	"	 
#	" 
#	"s   1C
C0
C-0
C>c                  0   [         R                  " S5      R                  SS5      n [         R                  " S5      R                  SS5      n[        R                  " [
        5         [        X5        S S S 5        [         R                  " S5      R                  SS5      n [         R                  " S5      R                  SS5      n[        R                  " [
        5         [        X5        S S S 5        g ! , (       d  f       N= f! , (       d  f       g = f)Nr6  re  r?   r1  r@   ri  )rO   rW  r   r   r   r   r   r8  s     rg   test_check_invalid_dimensionsr;     s     
2		q!	$B	2		q!	$B	z	"b% 
#	2		q!	$B	2		q!	$B	z	"b% 
#	"	 
#	" 
#	"s   'C6!D6
D
Dc                 r   [         R                  R                  S5      nUR                  S5      nU " U5      nUR                  S5      nU " U5      n[	        X55      u  pg[        U5      (       d   e[        X6-
  5      R                  5       S:X  d   e[        U5      (       d   e[        XW-
  5      R                  5       S:X  d   e[	        X35      u  ph[        U5      (       d   e[        X6-
  5      R                  5       S:X  d   e[        U5      (       d   e[        X-
  5      R                  5       S:X  d   eg r  )rO   rP   rQ   rR   r   r   r   r   )	r   rZ   r  	XA_sparser  	XB_sparser-  r.  XA_2_checkeds	            rg   test_check_sparse_arraysr@  .  s(    ))


"C			6	"Bb!I			6	"Bb!I29HJ Jy%&**,111Jy%&**,1114YJJJy%&**,111L!!!!|()--/1444r   c                     U R                   n[        U5      S:  a  [        S U  5       5      $ [        S U  5       5      $ )NrL   c              3   8   #    U  H  n[        U5      v   M     g 7frA  )tuplify)rD  rc   s     rg   rF  tuplify.<locals>.<genexpr>J  s     /QcWS\\Qs   c              3   $   #    U  H  ov   M     g 7frA  r   )rD  rs     rg   rF  rD  M  s     N1Qs   )rX   r  rV   )r[   rb  s     rg   rC  rC  E  s8    	A
1vz/Q/// NN""r   c                      [         R                  R                  S5      n U R                  S5      n[	        U5      nU R                  S5      n[	        U5      n[        X$5      u  pV[        X%5        [        XF5        g r  )rO   rP   rQ   rR   rC  r   r4   )rZ   r  	XA_tuplesr  	XB_tuplesr-  r.  s          rg   test_check_tuple_inputrJ  P  se    
))


"C			6	"BI			6	"BI29HJy-y-r   c                  L   [         R                  " [         R                  " S5      S5      R                  [         R                  5      n [         R                  " [         R                  " S5      S5      R                  [         R                  5      n[        U S 5      u  p#UR                  [         R                  :X  d   e[        X5      u  p#UR                  [         R                  :X  d   eUR                  [         R                  :X  d   e[        U R                  [        5      U5      u  p#UR                  [        :X  d   eUR                  [        :X  d   e[        XR                  [        5      5      u  p#UR                  [        :X  d   eUR                  [        :X  d   eg r)  )rO   r,  rW  rS   r  r   rT   r   r3  s       rg   test_check_preserve_typerL  \  s8   	299R=&	)	0	0	<B	299R=&	)	0	0	<B22t<Jrzz))) 32:Jrzz)))rzz))) 3299U3CRHJu$$$u$$$ 32yy7GHJu$$$u$$$r   r   
seuclideanmahalanobisdist_functionc                 ,   [        SS9   [        R                  R                  S5      nUR	                  S5      n[        [        XAS95      n[        R                  " [        U" XAU S95      5      n[        Xe5        S S S 5        g ! , (       d  f       g = f)Nr   )r7  r   r  rD   r   )
r   rO   rP   rQ   rR   r   r
   r;  rV   r2   )r   rE   rO  rZ   r[   expected_distr"  s          rg   +test_pairwise_distances_data_derived_paramsrR  t  sm     
s	+ii##A&i("5#:;yy}QfMNO, 
,	+	+s   A2B
Bc                    [         R                  R                  S5      nUR                  S5      nUR                  S5      n[        R
                  " [        SU  S3S9   [        X#U S9  S S S 5        g ! , (       d  f       g = f)Nr   r  z+The '(V|VI)' parameter is required for the z metricr   rD   )rO   rP   rQ   rR   r   r   r   r#   )rE   rZ   r[   r^   s       rg   1test_pairwise_distances_data_derived_params_errorrT    so     ))


"C)$A)$A	<VHGL
 	1/	
 
 
s   A33
B)
braycurtiscanberra	chebyshevcorrelationhammingrN  r	   rM  r  r   r   rC   y_is_xzY is Xz
Y is not Xc           	         [         R                  R                  S5      nUR                  S5      R	                  USS9n0 nU(       a  Un[        [        X@S95      nOUR                  S5      R	                  USS9n[        XFU S9nU S:X  a=  S[         R                  " [         R                  " XF/5      SS[         R                  S	90nOdU S
:X  a^  S[         R                  R                  [         R                  " [         R                  " XF/5      R                  5      5      R                  0n[        XF4SU 0UD6n[!        X5        g )Nr   r>   FrA   rD   rM  VrL   )r   ddofrT   rN  VIrE   )rO   rP   rQ   rR   rS   r   r
   r   varr;  r   r   invcovr  r#   r2   )	rE   rY   rZ  rZ   r[   r   r^   rQ  r"  s	            rg   )test_numeric_pairwise_distances_datatypesrb    s   6 ))


"C&!((E(BAF"5#:;f%,,\,Fa62\!266"))QF"3!!2::VWF}$BIIMM"&&A61B1D1D*EFHHIFa<6<V<DD(r   pairwise_distances_funcc                     SS/S[         R                  /SS/SS//nU " XSS9n[         R                  " U5      R                  5       (       a   eg)	z8Check that `nan_euclidean` is lenient with `nan` values.r   rL   rG   r   r?   rI   rD   N)rO   rU   isnanany)rc  r[   outputs      rg   test_nan_euclidean_supportrh    sS     Q!RVVq!fq!f-A$Q/BFxx##%%%%%r   c                     [         R                  [         R                  /[         R                  [         R                  /[         R                  [         R                  //n [        X SS9nSS/SS/SS//n[        X"SS9n[        X5        g)zqCheck that the behavior of constant input is the same in the case of
full of nan vector and full of zero vector.
rI   rD   r   N)rO   rU   r$   r2   )X_nan
argmin_nanX_constargmin_consts       rg   (test_nan_euclidean_constant_input_argminrn    so    
 ffbff/"&&"&&1ABE*5PJ1v1v1v&G,WoVLJ-r   zX,Y,expected_distance)r  ababc)r  r   r   )r   r  r   )r   r   r  r  ro  r  r   r   c                 4    S n[        XUS9n[        XB5        g)z8Check pairwise_distances with lists of strings as input.c                 X    [         R                  " [        U 5      [        U5      -
  5      $ rA  )rO   r   r  r   s     rg   dummy_string_similarityLtest_pairwise_dist_custom_metric_for_string.<locals>.dummy_string_similarity  s    vvc!fs1vo&&r   )r[   r^   rE   N)r#   r2   )r[   r^   expected_distancers  actual_distances        rg   +test_pairwise_dist_custom_metric_for_stringrw    s    $' )1:QROO7r   c                      S n [         R                  " / SQ/ SQ/ SQ/[        S9n[         R                  " / SQ/ SQ/ SQ/5      n[        XS	9n[	        X25        g
)zbCheck that pairwise_distances does not convert boolean input to float
when using a custom metric.
c                 R    SX-  R                  5       X-  R                  5       -  -
  $ )NrL   )r   )v1v2s     rg   dummy_bool_distBtest_pairwise_dist_custom_metric_for_bool.<locals>.dummy_bool_dist  s"    BG==?bg]]_444r   )rL   r   r   r   )rL   r   rL   r   )rL   rL   rL   rL   r   )r  rK         ?)rK   r  rK   )r~  rK   r  )r[   rE   N)rO   r   r   r#   r2   )r|  r[   ru  rv  s       rg   )test_pairwise_dist_custom_metric_for_boolr    sQ    
5 	,l;4HA	
 )1EOO7r   c                     U " [         R                  " S5      5      /nU " [         R                  " S5      5      /n[        SSS9" S [        X5       5       5        g )Nr   rG   r   )r   
max_nbytesc              3   P   #    U  H  u  p[        [        5      " X5      v   M     g 7frA  )r<   r   )rD  m1m2s      rg   rF  9test_sparse_manhattan_readonly_dataset.<locals>.<genexpr>  s$      %;T#$R,,;Ts   $&)rO   r  r;   zip)r   	matrices1	matrices2s      rg   &test_sparse_manhattan_readonly_datasetr    sQ     rwwv/0Irwwv/0I A!$ %;>y;T% r   )rC   )r   typesr   numpyrO   r   r   scipy.sparser   scipy.spatial.distancer   r   r   r	   r
   r   sklearnr   sklearn.exceptionsr   sklearn.metrics.pairwiser   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*   sklearn.preprocessingr+   sklearn.utils._array_apir,   r-   r.   r/   r0   sklearn.utils._testingr1   r2   r3   r4   r5   sklearn.utils.fixesr6   r7   r8   r9   r:   sklearn.utils.parallelr;   r<   rh   markparametrizer}   r   DeprecationWarningr   r   r   r   rW  rS   _minkowski_kwdsr   r   r  r   r   r   r   r   r   r   r   itemsr   r  thread_unsafer-  r2  r=  rI  rY  r   r   rl  rv  r  r  filterwarningsr  r   rs   r  r  r  r  r  r  r  r  paramxfailr  r  infr  r  rU   r  r  r  r  r  r  r  r  r  r	  r  r  r  r  r  r#  r'  r/  r4  r9  r;  r@  rC  rJ  rL  rR  rT  rb  rh  rn  rw  r  r  )rR  rT  	containers   000rg   <module>r     s~        !  # 4       > ,    5L^ )-/'  
 #GH' I'B .9.9.9.9>< : : : :><F ,-#=>: ? .:<2 "46F!GH% I%:F 		!Q..xe.D1MO 	["-	
 	

 
<(A7	.#?$ 2::rzz3"?@ A%& )-/'  
 	["-	["-	<(A7	.#?	;#7	5	5BI T .9 :	
> )-/'  
 T'	'D&
7 )9)?)?)AB.9& : C&*, .9.9vA : : vAr7"
4  $+ -,N )7I, - -N )7 	

 -()( ) 5(	
 :0	

 
!*.QRA	
  H	
 ,<	
->?>7  #EF? G? #EF7 G7 NO! P!H XXJ		HI''HH  
 XXJ		HI''HH  
%% XXJ		HI''HH  
%
%@ tUm4 5 ED XXJ		HI''HH  
 XXJ		HI''HH  
++* XXJ		HI''HH  
+
+$ '89XXJ		HI''HH  
 XXJ		HI''HH  
4 :4( '89XXJ		HI''HH  
4 :4" 	T4 JJ++##+M#N		
 G-	4 .	4 T5M23 33 rxx"&&!61rvvg,8PQRrxx"&&!61rvvg,8PRVWX) Y S) 	Aq6Aq6"	#RWWQZ8	Aq6Arvv;'	("''!*bff=	BFFA;BFF,	-rvvrvv>	BFFA;,	-rwwqz266B	Arvv;BFF,	-rwwqz266B	Aq6Aq6"	#RWWQZ4	Aq6Ar7#	$bggaj"5	B7QG$	%rvvr2	B7RG$	%rwwqz26	Ar7QG$	%rwwqz26--  2662,7	$ 8	$ 2662,71 81h 2662,7 8(<6###31l 
 
  
 .9 :
C-	/ +< => .9& :	
&* .9  : ('' 	$& .95 :5,#	.%0 Aq6*L-#@A(*DE
- B +

- L-#@A0 B0 " D%=x6NO) P#$)> 24QR&	&.  o?	
 $K3Z#sc3Z0	
8882 .9 : i H I
 I Ip I
 I4 I0 I
 I4 IsB   /||
;|||||#|(|-|2?|7