o
    EbE                    @   s  d dl Zd dlmZmZ d dlZd dlZd dlZd dl	m
Z
 d dlmZmZmZmZmZmZmZ d dlZd dlmZ d dlZd dlmZ d dlmZmZmZmZmZmZm Z m!Z!m"Z"m#Z# d dlm$Z$m%Z%m&Z&m'Z'm(Z(m)Z)m*Z*m+Z+m,Z,m-Z-m.Z.m/Z/m0Z0m1Z1m2Z2m3Z3m4Z4m5Z5m6Z6m7Z7m8Z8m9Z9m:Z: g d	Z;ej<g d
g dg dg dg dg dgddZ=ee=Z>i Z?dd Z@e@  dd ZAdoddZB					dpddZCdd ZDdd ZEdod d!ZFeed"d#d$d% dfd&d'ZGd(dd)d% d*ddddddddddd+d eed"d#fd,d-ZHeHed.dd/ZIeHed.ddd0ZJeHed.dd.d1ZKeHed.ddd.d2ZLeHe4ZMeHe,dd3 ZNZOeHe:ZPeHe*ZQeHe'ZReHe&ZSeHe)ZTeHe(ZUeHe/ZVeHe0ZWeHe-ZXeHe+dd4ZYeHe9dd4ZZeHe$Z[eHe%dd4Z\eHe8Z]eHe7Z^eHe5Z_G d5d6 d6Z`G d7d8 d8ZaG d9d: d:Zbd;d< Zcec  G d=d> d>ZdG d?d@ d@ZeG dAdB dBZfdCdD ZgG dEdF dFZhdGdH ZiG dIdJ dJZjejkldKg dLdMdN ZmdOdP ZndQdR ZodSdT ZpdUdV ZqdWdX ZrdYdZ Zsd[d\ Ztd]d^ Zud_d` Zvdadb Zwdcdd Zxdedf Zydgdh Zzdidj Z{dkdl Z|dmdn Z}dS )q    N)wrapspartial)norm)verboseassert_assert_array_equalassert_equalassert_almost_equalassert_allclosesuppress_warnings)raises)_distance_pybind)

squareformpdistcdist	num_obs_y
num_obs_dmis_valid_dm
is_valid_y_validate_vector_METRICS_NAMES_METRICS)
braycurtiscanberra	chebyshev	cityblockcorrelationcosinedice	euclideanhammingjaccardjensenshannon	kulsinskikulczynski1mahalanobismatching	minkowskirogerstanimoto
russellrao
seuclideansokalmichenersokalsneathsqeuclideanyule)zcdist-X1.txtzcdist-X2.txtziris.txtzpdist-boolean-inp.txtzpdist-chebyshev-ml-iris.txtzpdist-chebyshev-ml.txtzpdist-cityblock-ml-iris.txtzpdist-cityblock-ml.txtzpdist-correlation-ml-iris.txtzpdist-correlation-ml.txtzpdist-cosine-ml-iris.txtzpdist-cosine-ml.txtzpdist-double-inp.txtzpdist-euclidean-ml-iris.txtzpdist-euclidean-ml.txtzpdist-hamming-ml.txtzpdist-jaccard-ml.txtzpdist-jensenshannon-ml-iris.txtzpdist-jensenshannon-ml.txtzpdist-minkowski-3.2-ml-iris.txtzpdist-minkowski-3.2-ml.txtzpdist-minkowski-5.8-ml-iris.txtzpdist-seuclidean-ml-iris.txtzpdist-seuclidean-ml.txtzpdist-spearman-ml.txtzrandom-bool-data.txtzrandom-double-data.txtzrandom-int-data.txtzrandom-uint-data.txt)r     m         )r/   r   '        )r0   r4   r     4     )r1   r5   r8   r      e  )r2   r6   r9   r;   r     )r3   r7   r:   r<   r=   r   doubledtypec                  C   s   t D ]'} | dddd}tjtjtd| }t|}t	|t
|< |  qtt
d t
d< tt
d t
d< tt
d t
d< tt
d	 t
d	< tt
d
 t
d
< d S )Nz.txt z-mldatapdist-boolean-inprandom-bool-datarandom-double-datarandom-float32-datarandom-int-datarandom-uint-data)
_filenamesreplaceospathjoindirname__file__opennpZloadtxteocloseZbool_float32int_uint)fnnameZfqfnfp rZ   C/usr/lib/python3/dist-packages/scipy/spatial/tests/test_distance.pyload_testing_filesq   s   
r\   c                   C   s   t djdk S )Nr      )rQ   ZintpitemsizerZ   rZ   rZ   r[   	_is_32bit   s   r_   c                    s|   dd  D  |d u rdd  D  d}t dd  D  |dk r9t fdd D s0tdt d j| } |f S )	Nc                 S      g | ]}t |qS rZ   )rQ   
asanyarray.0arZ   rZ   r[   
<listcomp>       z!_chk_asarrays.<locals>.<listcomp>c                 S   s$   g | ]}|j d krt|n|qS    )ndimrQ   Zravelrb   rZ   rZ   r[   re      s    r   c                 s   s    | ]}t |V  qd S N)rQ   
atleast_1drb   rZ   rZ   r[   	<genexpr>       z _chk_asarrays.<locals>.<genexpr>c                 3   s     | ]}|j  d  j kV  qdS )r   N)ri   rb   arraysrZ   r[   rl          z(array ndim must be the same for neg axis)tupleall
ValueErrorrangeri   )ro   axisrZ   rn   r[   _chk_asarrays   s   
rv   FTc
                    sp  t |  d}
|
d d |
d }  |o| }|s#|r#tdd | D }|r@dd | D }t|r@d}d}tdd t| |D } d urJtn|rWt| d	 j  n|  f S |	rct|rkt	|  t
 fd
d| D s{td|rd	k  rtd|rtd	kd	 jjk rt fdd| D }  |rdk
 rd |  f S )Nru   c                 s   s$    | ]}t j|t jjuV  qd S rj   )rQ   magetmasknomaskrb   rZ   rZ   r[   rl      s   " z_chk_weights.<locals>.<genexpr>c                 S   s   g | ]
}t t |qS rZ   )rQ   Zisnansumrb   rZ   rZ   r[   re          z _chk_weights.<locals>.<listcomp>Tc                 s   s(    | ]\}}|rt j|n|V  qd S rj   rQ   ry   Zmasked_invalid)rc   rd   Zhas_nanrZ   rZ   r[   rl      s    r   c                 3   s"    | ]}j |j   fkV  qd S rj   shaperb   )ru   weightsrZ   r[   rl      s     z*weights shape must match arrays along axiszweights cannot be negativec                 3   s     | ]}t j| d V  qdS )rw   N)rQ   takerb   )ru   pos_weightsrZ   r[   rl      rp   rh   )rv   anyrq   ziprQ   ra   onesr   _freq_weights_weight_maskedrr   rs   Znonzerosize)ro   r   ru   force_weightsZsimplify_weightsZpos_onlyZ	neg_checkZ
nan_screenmask_screenddofchkedZhas_nansrZ   )ru   r   r   r[   _chk_weights   sD   
r   c                 C   s2   | d u r| S |  t}| |k rtd|  |S )Nz2frequency (integer count-type) weights required %s)astypeintr   rs   )r   Zint_weightsrZ   rZ   r[   r      s   
r   c                    s    d u rd t |}| D ]1}t j|}|t jju rq|jdkr5t fddt|jD }|j|d}|d|	t
 9 }q|S )Nr   rh   c                 3   s    | ]	}| kr|V  qd S rj   rZ   rc   irw   rZ   r[   rl      s    z!_weight_masked.<locals>.<genexpr>rw   )rQ   ra   ry   rz   r{   ri   rq   rt   r   r   r   )ro   r   ru   rd   Z	axis_maskZnot_axesrZ   rw   r[   r      s   

r   c           	         s   t j|t jd}t j|}dd | d j  t fdd| D s)J dtt| D ],}|	| | }|
 }|| |< t |d| | } fd	d
| D } q1| |fS )Nr?   c                    s6   t t j |d}| fddt jD S )Nrw   c                    s"   g | ]}|kr j | nd qS rg   r   r   rd   ru   rZ   r[   re      s    z/_rand_split.<locals>.mytake.<locals>.<listcomp>)rQ   ra   r   reshapert   ri   )rd   Zixru   recordrZ   r   r[   mytake   s   z_rand_split.<locals>.mytaker   c                 3   s    | ]
}|j   kV  qd S rj   r   rb   )ru   n_obsrZ   r[   rl      s    z_rand_split.<locals>.<genexpr>z#data must be aligned on sample axis      ?c              	      s&   g | ]}t j|| d  d qS rw   rQ   appendrb   )ru   r   split_ixrZ   r[   re      s
    
z_rand_split.<locals>.<listcomp>)rQ   arrayfloat64randomRandomStater   rr   rt   r   Zrandintrandr   )	ro   r   ru   	split_perseedZseeded_randr   Zprev_wqrZ   )ru   r   r   r   r[   _rand_split   s    r   h㈵>atolc                 C      | S rj   rZ   xrZ   rZ   r[   <lambda>       r   c           	   	   C   s   || }||}zt ||k r||| W d S W d S  ty+   ||| Y d S  ttfyG   t||D ]\}}t|||d q8Y d S w )N)compare_assert)rQ   r   r   AttributeError	TypeErrorrs   r   _rough_check)	rd   br   keywZcheck_aZcheck_bZa_iZb_irZ   rZ   r[   r      s   r      c                 C   r   rj   rZ   r   rZ   rZ   r[   r     r   r   r   c                    s8   t  	
fdd}|S )zruns fn on its arguments 2 or 3 ways, checks that the results are the same,
       then returns the same thing it would have returned beforec               
      s  | i |}| d 	 }| 	d  }| d }| d t|| ddd}|d d |d |d }} rAdd |D }zt|| } r[||< t|| i |d r|d	 |< t|| i |d |d
 |< zt|| i |d W n ty } z
t||||f|d }~ww rrt|| d\}}	||	d |j< d|	|jd < t|| }
|	|< t||
i |d |D ]}td g|j	 }t|jd | < t|}|| d ||< qt|| }
t||
i |d |D ] }td g|j	 }t|jd | < t|}|| t
j ||< q| dd dkr@
r@t|| }
t||
i |d rodd |D }t|| }
t||
i |d rod |< t||
i |d ~~
~	r fdd|D }t
||d }t|| }||< t||i |d ~~~rˈdkrt|| d\}}t|| }||< t||i |d W |S W |S W |S  ty } zstdj|f  W Y d }~|S W Y d }~|S d }~ww )Nru   T)r   ru   r   r   rx   c                 S   s   g | ]	}t | qS rZ   )rQ   rk   squeezerb   rZ   rZ   r[   re     s    z4_weight_checked.<locals>.wrapped.<locals>.<listcomp>)r   g     @Y@gB`"۹?)r   r   r   e   Z
nan_policyZomitc                 S   s   g | ]}t j|qS rZ   r~   rb   rZ   rZ   r[   re   D  s    c                    s   g | ]
}t j|| d qS r   r   rb   rw   rZ   r[   re   N  r}          @z%s NotImplemented weights: %s)getr   rq   r   	Exceptiontyper   r   sliceri   rQ   nanr   NotImplementedErrorwarningswarn__name__)argskwargsresultro   restr   r   eZ
dud_arraysZdud_weightsZdud_argsrd   ZindexerZ
dup_arraysZdup_weightsZdup_argsZsplit_arraysZsplit_weightsZ
split_args
const_testdefault_axisdud_testdup_testrW   r   ma_safema_very_safen_argsnan_safe	ones_testr   silentr   
split_testr   
weight_argrw   r[   wrapped  s    


z _weight_checked.<locals>.wrapped)r   )rW   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   rZ   r   r[   _weight_checked  s   .Ur   rh   )r   r   )r   r   r   )r   r   r   )r   r   r   r   )r   )r   c                   @   s   e Zd Zdd Zdd Zdd Zejddd	d
dddde	j
gdd Zdd Zdd Zdd Zd%ddZdd Zdd Zdd Zd d! Zd"d# Zd$S )&	TestCdistc                 C   F   g d| _ tjtjtjtjgtjtjtjgtjtjgtjgd| _d S N)rF   rG   rH   rE   rD   )boolrV   r   rT   rnd_eo_namesrQ   rV   rU   rT   r>   valid_upcastsselfrZ   rZ   r[   setup_method~     

zTestCdist.setup_methodc              
   C   s  ddd}g dg dg dg}g dg d	g d
g}dt dd}dgd }tD ][}ttt||fd|i| ttt||fdt|i| ttt||fdd| i| ttt||g|R d|i ttt||g|R dt|i ttt||g|R dd| i q*ttt||| ttt|||g|R   ttt|||fi | ttt|||ddd ttt|||dddd ttt|||ddd ttt|||dd ttt|||d ttt|||dddd tt|||dddd d S )Nrh   r   c                 S      || | S rj   rZ   r   yargkwargkwarg2rZ   rZ   r[   
_my_metric     z3TestCdist.test_cdist_extra_args.<locals>._my_metricr   r         @333333?ffffff@333333@皙@r   皙@)g      @      @g       @)g      @333333@g @)g      @r   r   Q	@   ZN0tV4l1D_p4raMr      metrictest_r   ffffff
@r   r   r   皙?r   r   r   皙@rh   r   )rQ   aranger   assert_raisesr   r   evalr
   )r   r   X1X2r   r   r   rZ   rZ   r[   test_cdist_extra_args  sl   


zTestCdist.test_cdist_extra_argsc                 C   sD   d}t d }t d }t||d}t||d}t|||tdkd d S )NV瞯<cdist-X1cdist-X2r   test_euclideanr   rtolr   rR   wcdist_no_constr
   r   )r   epsr  r  Y1Y2rZ   rZ   r[   #test_cdist_euclidean_random_unicode  s   z-TestCdist.test_cdist_euclidean_random_unicodep皙?      ?r   Gz?r   gffffff@gffffff@c                 C   sN   d}t d }t d }t||d|d}t||d|d}t||d|tdkd	 d S )
NvIh%<=r	  r
  r'   r  test_minkowskir   r   )r   r  r   r  )r   r  r  r  r  r  r  rZ   rZ   r[   test_cdist_minkowski_random  s   z%TestCdist.test_cdist_minkowski_randomc                 C   sb   d}t d }t d }t||d}dd }dt||| ||| j }t|||tdkd	 d S )
N+=r	  r
  r   c                 S   s   t jj| ddddS )Nrh   rw   rx   )rQ   linalgr   r   )XrZ   rZ   r[   norms     z1TestCdist.test_cdist_cosine_random.<locals>.normsrh   r   r  )rR   wcdistrQ   dotTr
   r   )r   r  r  r  r  r  r  rZ   rZ   r[   test_cdist_cosine_random  s   "z"TestCdist.test_cdist_cosine_randomc                 C   s   t dgdgg}t dgdgg}t||dd}t|dt dgt dt dgg t d	d	gd
d	gg}t d	dgdd	gd	dgg}t||dd}t d}t||||gdd| dgg tttd	dggddggdd d S )Nr   r      r%   r                 @      ?r   rx   rh   r   )rQ   r   r   r
   sqrtr  rs   )r   x1x2distrt2rZ   rZ   r[   test_cdist_mahalanobis  s   (

z TestCdist.test_cdist_mahalanobisc                    sV   G dd d  fdd}t j  ggtd}t|||d}d}t||tdkd	 d S )
Nc                   @      e Zd ZdS )z6TestCdist.test_cdist_custom_notdouble.<locals>.myclassNr   
__module____qualname__rZ   rZ   rZ   r[   myclass      r4  c                    (   t | d  rt |d  stddS Nr   zType has been changed+?
isinstancers   r   r   r4  rZ   r[   r        z9TestCdist.test_cdist_custom_notdouble.<locals>._my_metricr?   r&  r8  r   r   )rQ   r   objectr   r   r   )r   r   rB   Zcdist_yright_yrZ   r<  r[   test_cdist_custom_notdouble  s   z%TestCdist.test_cdist_custom_notdoubleHz>c              
   K   s  z't ||fd|i|}t ||fdt|i|}t ||fdd| i|}W nN tyu }	 zB|	j}
tdkr?t|
j t|	 t|
t ||fd|i| t|
t ||fdt|i| t|
t ||fdd| i| W Y d }	~	d S d }	~	ww t|||tdkd t|||tdkd d S Nr   r   r   r  )	r   r  r   	__class__r   printr   r  r
   )r   r  r  r   r  r   y1y2y3r   e_clsrZ   rZ   r[   _check_calling_conventions  s    
,	z$TestCdist._check_calling_conventionsc                 C   s   | j D ]}t| d d dd d df }t| dd dd d df }tD ]i}tdkr0td|d| |dv r9d|vr9q#| ||| |d	krat||gtj	}tj
|d
dd}| j||||d q#|dkrt||gtj	}tt|j}ttj|j}| j||||d q#qd S )Nr%  r   rh   r   	testing:  with: 	   r   r.   r&   r#   r)   r$   r,   r+   r(   r   r*   r   ru   r   Vr%   VI)r   rR   r   r   rE  rJ  rQ   vstackr   r>   var
atleast_2dcovr#  r   r  inv)r   eo_namer  r  r   ZX12rP  rR  rZ   rZ   r[   test_cdist_calling_conventions   s.   
z(TestCdist.test_cdist_calling_conventionsc                 C   sX  d}t d | jd ft d | jd ft d | jd ft d | jd	 fg}tD ]}|D ]~}|d
 d d dd d df }|d
 dd dd d df }z	t|||d}W n: ty } z.|j}	tdkrit|	j t| |d D ]}
|
|}|
|}t	|	t|||d qmW Y d }~q*d }~ww |d D ]}
t|
||
||d}t
|||tdkd qq*q&d S )NrB  rD   r   rH   rV   rG   r   rF   rT   r   r%  r   rh   r   r&  r  )rR   r   r   r   r   rD  r   rE  r   r  r
   )r   r  testsr   testr  r  rF  r   rI  new_typeZX1newZX2newrG  rZ   rZ   r[   test_cdist_dtype_equivalence   s<   

z&TestCdist.test_cdist_dtype_equivalencec                 C   s  d}t d }t d }|jd |jd }}tD ]}t }|dkr$d|d< tj||ftjd}t|||fi |}	t|||fd	|i|}
t|	|
|t	d
kd t
|
|u  tj|d |d ftjd}ttt|||fd	|i| tjd
| d
| ftjdd d d
d d d
f }tj||ftjdd}ttt|||fd	|i| ttt|||fd	|i| tj||ftjd}ttt|||fd	|i| qd S )Nr  r	  r
  r   r'   r  r  r?   outr   r  rh   F)r@   Zorder)rR   r   r   dictrQ   emptyr>   r   r
   r   r   r  rs   int64)r   r  r  r  Zout_rZout_cr   r   out1r  r  out2out3Zout4out5rZ   rZ   r[   test_cdist_out;  sd   zTestCdist.test_cdist_outc           
      C   s   d}t d d d dd d df }t d d d dd d df }| }| }t|| t|| t|jj  t|jj  t|jj t|jj tD ]+}t }|dkrYd|d< t|||fi |}t|||fi |}	t	||	|t
dkd qLd S )	Nr  r	  r   r
  r'   r  r  r  )rR   copyr   r   flagsc_contiguousr   r`  r   r
   r   )
r   r  r  r  X1_copyZX2_copyr   r   r  r  rZ   rZ   r[   test_striding_  s&   

zTestCdist.test_stridingc                 C   s   t D ]>}tjdd}tjdd}t }|dkrd|d< t||fd|i|}dd |||fD }~~~tdd	 |D s@J qd S )
N
   r'   r  r  r   c                 S   r`   rZ   )weakrefref)rc   vrZ   rZ   r[   re     rf   z1TestCdist.test_cdist_refcount.<locals>.<listcomp>c                 s   s    | ]}| d u V  qd S rj   rZ   )rc   Zweak_refrZ   rZ   r[   rl     rm   z0TestCdist.test_cdist_refcount.<locals>.<genexpr>)r   rQ   r   r   r`  r   rr   )r   r   r+  r,  r   r^  Z	weak_refsrZ   rZ   r[   test_cdist_refcountz  s   zTestCdist.test_cdist_refcountNrB  )r   r2  r3  r   r  r  pytestmarkparametrizerQ   infr  r$  r/  rA  rJ  rY  r]  rg  rl  rq  rZ   rZ   rZ   r[   r   |  s"    	(


 $r   c                	   @   s  e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Ze	j
jdd Ze	j
jdd Ze	j
jdd Zdd Zdd Zdd Zdd Zdd Zdd Zdd  Zd!d" Zd#d$ Ze	j
jd%d& Ze	j
jd'd( Ze	j
jd)d* Zd+d, Zd-d. Zd/d0 Zd1d2 Ze	j
jd3d4 Ze	j
jd5d6 Z e	j
jd7d8 Z!d9d: Z"d;d< Z#d=d> Z$e	j
jd?d@ Z%e	j
jdAdB Z&e	j
jdCdD Z'e	j
(dEdFdGdHdIdJe)j*gdKdL Z+dMdN Z,dOdP Z-dQdR Z.e	j
jdSdT Z/e	j
jdUdV Z0e	j
jdWdX Z1e	j
jdYdZ Z2e	j
jd[d\ Z3e	j
jd]d^ Z4d_d` Z5dadb Z6dcdd Z7dedf Z8dgdh Z9didj Z:dkdl Z;dmdn Z<dodp Z=dqdr Z>dsdt Z?dudv Z@dwdx ZAdydz ZBd{d| ZCd}d~ ZDdd ZEdd ZFdd ZGdd ZHdd ZIdd ZJdd ZKdd ZLdd ZMdd ZNdd ZOdd ZPdd ZQdd ZRdd ZSdd ZTdd ZUdd ZVdd ZWdd ZXdd ZYdd ZZdd Z[dd Z\dd Z]e	j
jdd Z^dd Z_dd Z`dddZadd Zbdd Zcdd Zddd ZedS )	TestPdistc                 C   r   r   r   r   rZ   rZ   r[   r     r   zTestPdist.setup_methodc              	   C   s  ddd}ddgddgd	dgg}d
t dd}d
gd }tD ]U}ttt|fd|i| ttt|fdt|i| ttt|fdd| i| ttt|g|R d|i ttt|g|R dt|i ttt|g|R dd| i qttt|| ttt||g|R   ttt||fi | ttt||d	dd ttt||ddd	d ttt||dd	d ttt||dd	 ttt||d ttt||dd	dd tt||dddd d S )Nrh   r   c                 S   r   rj   rZ   r   rZ   rZ   r[   r     r   z3TestPdist.test_pdist_extra_args.<locals>._my_metricr   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r  )rQ   r  r   r  r   r   r  r
   )r   r   r  r   r   r   rZ   rZ   r[   test_pdist_extra_args  sR   



zTestPdist.test_pdist_extra_argsc                 C   0   d}t d }t d }t|d}t|||d d S NrB  pdist-double-inppdist-euclideanr   r  rR   wpdist_no_constr
   r   r  r  Y_rightY_test1rZ   rZ   r[   test_pdist_euclidean_random  
   
z%TestPdist.test_pdist_euclidean_randomc                 C   ry  rz  r~  r  rZ   rZ   r[   test_pdist_euclidean_random_u  r  z'TestPdist.test_pdist_euclidean_random_uc                 C   6   d}t td }td }t|d}t|||d d S rz  rQ   rT   rR   r  r
   r  rZ   rZ   r[   #test_pdist_euclidean_random_float32  
   
z-TestPdist.test_pdist_euclidean_random_float32c                 C   ry  )NrB  r{  r|  r  r}  r~  r   r  r  r  Y_test2rZ   rZ   r[    test_pdist_euclidean_random_nonC  r  z*TestPdist.test_pdist_euclidean_random_nonCc                 C   ry  )NrB  irispdist-euclidean-irisr   r}  r~  r  rZ   rZ   r[    test_pdist_euclidean_iris_double  
   
z*TestPdist.test_pdist_euclidean_iris_doublec                 C   <   d}t td }td }t|d}t|||tdkd d S )Nr   r  r  r   r   r  rQ   rT   rR   r  r
   r   r  rZ   rZ   r[   !test_pdist_euclidean_iris_float32  
   
z+TestPdist.test_pdist_euclidean_iris_float32c                 C   ry  )NrB  r  r  r  r}  r~  r  rZ   rZ   r[   test_pdist_euclidean_iris_nonC  
   
z(TestPdist.test_pdist_euclidean_iris_nonCc                 C   ry  )NrB  r{  pdist-seuclideanr*   r}  rR   r   r
   r  rZ   rZ   r[   test_pdist_seuclidean_random  r  z&TestPdist.test_pdist_seuclidean_randomc                 C   sb   d}t td }td }t|d}t|||d t j|ddd}t|d|d	}t|||d d S )
NrB  r{  r  r*   r}  r   rh   rN  rO  )rQ   rT   rR   r   r
   rT  )r   r  r  r  r  rP  r  rZ   rZ   r[   $test_pdist_seuclidean_random_float32  s   
z.TestPdist.test_pdist_seuclidean_random_float32c                 C   ry  )NrB  r{  r  test_seuclideanr}  r  r  rZ   rZ   r[   !test_pdist_seuclidean_random_nonC  r  z+TestPdist.test_pdist_seuclidean_random_nonCc                 C   ry  )NrB  r  pdist-seuclidean-irisr*   r}  r  r  rZ   rZ   r[   test_pdist_seuclidean_iris  r  z$TestPdist.test_pdist_seuclidean_irisc                 C   r  )Nr   r  r  r*   r}  )rQ   rT   rR   r   r
   r  rZ   rZ   r[   "test_pdist_seuclidean_iris_float32  s
   
z,TestPdist.test_pdist_seuclidean_iris_float32c                 C   ry  )NrB  r  r  r  r}  r  r  rZ   rZ   r[   test_pdist_seuclidean_iris_nonC  s
   
z)TestPdist.test_pdist_seuclidean_iris_nonCc                 C   ry  NrB  r{  pdist-cosiner   r}  rR   wpdistr
   r  rZ   rZ   r[   test_pdist_cosine_random$  r  z"TestPdist.test_pdist_cosine_randomc                 C   r  r  rQ   rT   rR   r  r
   r  rZ   rZ   r[    test_pdist_cosine_random_float32+  r  z*TestPdist.test_pdist_cosine_random_float32c                 C   ry  )NrB  r{  r  test_cosiner}  r  r  rZ   rZ   r[   test_pdist_cosine_random_nonC2  r  z'TestPdist.test_pdist_cosine_random_nonCc                 C   ry  )Nr   r  pdist-cosine-irisr   r   r  r  rZ   rZ   r[   test_pdist_cosine_iris:  r  z TestPdist.test_pdist_cosine_irisc                 C   r  )Nr   r  r  r   r   r   r   )rQ   rT   rR   r  r
   r   r  rZ   rZ   r[   test_pdist_cosine_iris_float32B  r  z(TestPdist.test_pdist_cosine_iris_float32c                 C   ry  )Nr   r  r  r  r   r  r  rZ   rZ   r[   test_pdist_cosine_iris_nonCJ  r  z%TestPdist.test_pdist_cosine_iris_nonCc                 C   sD   t t jdd}t ||g}tt|dd dkdd d S )Ni9  [   r   r   z&cosine distance should be non-negative)msg)rQ   absr   r   r   rS  r   r  )r   r   r  rZ   rZ   r[   test_pdist_cosine_boundsR  s
   
z"TestPdist.test_pdist_cosine_boundsc                 C   ry  NrB  r{  pdist-cityblockr   r}  r~  r  rZ   rZ   r[   test_pdist_cityblock_random[  r  z%TestPdist.test_pdist_cityblock_randomc                 C   r  r  r  r  rZ   rZ   r[   #test_pdist_cityblock_random_float32b  r  z-TestPdist.test_pdist_cityblock_random_float32c                 C   ry  )NrB  r{  r  test_cityblockr}  r~  r  rZ   rZ   r[    test_pdist_cityblock_random_nonCi  r  z*TestPdist.test_pdist_cityblock_random_nonCc                 C   ry  )Nr  r  pdist-cityblock-irisr   r}  r~  r  rZ   rZ   r[   test_pdist_cityblock_irisp  r  z#TestPdist.test_pdist_cityblock_irisc                 C   r  )Nr   r  r  r   r   r  r  r  rZ   rZ   r[   !test_pdist_cityblock_iris_float32x  r  z+TestPdist.test_pdist_cityblock_iris_float32c                 C   ry  )Nr  r  r  r  r}  r~  r  rZ   rZ   r[   test_pdist_cityblock_iris_nonC  r  z(TestPdist.test_pdist_cityblock_iris_nonCc                 C   ry  NrB  r{  pdist-correlationr   r}  r  r  rZ   rZ   r[   test_pdist_correlation_random  r  z'TestPdist.test_pdist_correlation_randomc                 C   r  r  r  r  rZ   rZ   r[   %test_pdist_correlation_random_float32  r  z/TestPdist.test_pdist_correlation_random_float32c                 C   ry  )NrB  r{  r  test_correlationr}  r  r  rZ   rZ   r[   "test_pdist_correlation_random_nonC  r  z,TestPdist.test_pdist_correlation_random_nonCc                 C   ry  )NrB  r  pdist-correlation-irisr   r}  r  r  rZ   rZ   r[   test_pdist_correlation_iris  r  z%TestPdist.test_pdist_correlation_irisc                 C   s<   d}t d }tt d }t|d}t|||tdkd d S )NrB  r  r  r   r   r  )rR   rQ   rT   r  r
   r   r  rZ   rZ   r[   #test_pdist_correlation_iris_float32  s
   
z-TestPdist.test_pdist_correlation_iris_float32c                 C   ry  )NrB  r  r  r  r}  r  r  rZ   rZ   r[    test_pdist_correlation_iris_nonC  r  z*TestPdist.test_pdist_correlation_iris_nonCr  r  r  r   r   皙	@c                 C   s<   d}t d }t|d|d}t|d|d}t||d|d d S )Nr  r{  r'   r  r  r   )r   r  r~  )r   r  r  r  r  r  rZ   rZ   r[   test_pdist_minkowski_random_p  s
   z'TestPdist.test_pdist_minkowski_random_pc                 C   4   d}t d }t d }t|ddd}t|||d d S NrB  r{  pdist-minkowski-3.2r'   r  r  r}  r~  r  rZ   rZ   r[   test_pdist_minkowski_random  
   z%TestPdist.test_pdist_minkowski_randomc                 C   :   d}t td }td }t|ddd}t|||d d S r  r  r  rZ   rZ   r[   #test_pdist_minkowski_random_float32  s
   z-TestPdist.test_pdist_minkowski_random_float32c                 C   r  )NrB  r{  r  r  r  r  r}  r~  r  rZ   rZ   r[    test_pdist_minkowski_random_nonC  r  z*TestPdist.test_pdist_minkowski_random_nonCc                 C   r  )NrB  r  pdist-minkowski-3.2-irisr'   r  r  r}  r~  r  rZ   rZ   r[   test_pdist_minkowski_3_2_iris  
   z'TestPdist.test_pdist_minkowski_3_2_irisc                 C   r  )Nr   r  r  r'   r  r  r}  r  r  rZ   rZ   r[   %test_pdist_minkowski_3_2_iris_float32  s
   z/TestPdist.test_pdist_minkowski_3_2_iris_float32c                 C   r  )NrB  r  r  r  r  r  r}  r~  r  rZ   rZ   r[   "test_pdist_minkowski_3_2_iris_nonC  r  z,TestPdist.test_pdist_minkowski_3_2_iris_nonCc                 C   r  )NrB  r  pdist-minkowski-5.8-irisr'   r   r  r}  r~  r  rZ   rZ   r[   test_pdist_minkowski_5_8_iris  r  z'TestPdist.test_pdist_minkowski_5_8_irisc                 C   s@   d}t td }td }t|ddd}t|||tdkd d S )	Nr   r  r  r'   r   r  r   r  r  r  rZ   rZ   r[   %test_pdist_minkowski_5_8_iris_float32  s
   z/TestPdist.test_pdist_minkowski_5_8_iris_float32c                 C   r  )NrB  r  r  r  r   r  r}  r~  r  rZ   rZ   r[   "test_pdist_minkowski_5_8_iris_nonC  r  z,TestPdist.test_pdist_minkowski_5_8_iris_nonCc                 C   s   t g ddd}t|dd}t|dt dt dt dt dt d	g t d
d
gdd
gd
dgdd
gd
dgg}t|dd}t d}t|||||dd| ddd| dg
 tttd
dgddggdd d S )N)r   r   r   r   rx   rh   r%   r&  r'  r)  r(  r   r   r   r   r   )	rQ   r   r   r   r
   r*  r  rs   r  )r   r   r-  r.  rZ   rZ   r[   test_pdist_mahalanobis  s   (
&
z TestPdist.test_pdist_mahalanobisc                 C   ry  Nr  rC   pdist-hammingr    r}  r  r  rZ   rZ   r[   test_pdist_hamming_random  r  z#TestPdist.test_pdist_hamming_randomc                 C   r  r  r  r  rZ   rZ   r[   !test_pdist_hamming_random_float32  r  z+TestPdist.test_pdist_hamming_random_float32c                 C   ry  Nr  rC   r  Ztest_hammingr}  r  r  rZ   rZ   r[   test_pdist_hamming_random_nonC#  r  z(TestPdist.test_pdist_hamming_random_nonCc                 C   r  r  rQ   r   rR   r  r
   r  rZ   rZ   r[   test_pdist_dhamming_random*  r  z$TestPdist.test_pdist_dhamming_randomc                 C   r  r  r  r  rZ   rZ   r[   "test_pdist_dhamming_random_float321  r  z,TestPdist.test_pdist_dhamming_random_float32c                 C   r  r  r  r  rZ   rZ   r[   test_pdist_dhamming_random_nonC8  r  z)TestPdist.test_pdist_dhamming_random_nonCc                 C   ry  N:0yE>rC   pdist-jaccardr!   r}  r  r  rZ   rZ   r[   test_pdist_jaccard_random?  r  z#TestPdist.test_pdist_jaccard_randomc                 C   r  r  r  r  rZ   rZ   r[   !test_pdist_jaccard_random_float32F  r  z+TestPdist.test_pdist_jaccard_random_float32c                 C   ry  Nr  rC   r  test_jaccardr}  r  r  rZ   rZ   r[   test_pdist_jaccard_random_nonCM  r  z(TestPdist.test_pdist_jaccard_random_nonCc                 C   r  r  r  r  rZ   rZ   r[   test_pdist_djaccard_randomT  r  z$TestPdist.test_pdist_djaccard_randomc                 C   r  r  r  r  rZ   rZ   r[   "test_pdist_djaccard_random_float32[  r  z,TestPdist.test_pdist_djaccard_random_float32c                 C   ,   d}t tdd}ttd||d d S )Nr  r%  r   r!   rm  r}  r   rQ   zerosr
   r   r  YrZ   rZ   r[   test_pdist_djaccard_allzerosb     z&TestPdist.test_pdist_djaccard_allzerosc                 C   r  r  r  r  rZ   rZ   r[   test_pdist_djaccard_random_nonCg  r  z)TestPdist.test_pdist_djaccard_random_nonCc                 C   ry  )Ndy=r{  pdist-jensenshannonr"   r}  r  r  rZ   rZ   r[   test_pdist_jensenshannon_randomn  r  z)TestPdist.test_pdist_jensenshannon_randomc                 C   r  )Nr  r{  r  r"   r   r  rQ   rT   rR   r   r
   r   r  rZ   rZ   r[   'test_pdist_jensenshannon_random_float32u  
   
z1TestPdist.test_pdist_jensenshannon_random_float32c                 C   ry  )Nr  r{  r  test_jensenshannonr}  r  r  rZ   rZ   r[   $test_pdist_jensenshannon_random_nonC|  r  z.TestPdist.test_pdist_jensenshannon_random_nonCc                 C   s<   t  rd}nd}td }td }t|d}t|||d d S )NgLa㧝=g-q=r  pdist-jensenshannon-irisr"   r   )r_   rR   r   r
   r  rZ   rZ   r[   test_pdist_jensenshannon_iris  s   
z'TestPdist.test_pdist_jensenshannon_irisc                 C   r  )Ngư>r  r   r"   r   r  r  r  rZ   rZ   r[   %test_pdist_jensenshannon_iris_float32  r  z/TestPdist.test_pdist_jensenshannon_iris_float32c                 C   ry  )Ng-C6
?r  r   r  r}  r  r  rZ   rZ   r[   "test_pdist_jensenshannon_iris_nonC  r  z,TestPdist.test_pdist_jensenshannon_iris_nonCc                 C   r  )Nr  r  r  rm  r}  r  r  rZ   rZ   r[   !test_pdist_djaccard_allzeros_nonC  r  z+TestPdist.test_pdist_djaccard_allzeros_nonCc                 C   ry  )Nr  r{  pdist-chebyshevr   r}  r  r  rZ   rZ   r[   test_pdist_chebyshev_random  r  z%TestPdist.test_pdist_chebyshev_randomc                 C   r  )NrB  r{  r  r   r   r  r  r  rZ   rZ   r[   #test_pdist_chebyshev_random_float32  r  z-TestPdist.test_pdist_chebyshev_random_float32c                 C   ry  )Nr  r{  r  test_chebyshevr}  r  r  rZ   rZ   r[    test_pdist_chebyshev_random_nonC  r  z*TestPdist.test_pdist_chebyshev_random_nonCc                 C   ry  )Nr  r  pdist-chebyshev-irisr   r}  r  r  rZ   rZ   r[   test_pdist_chebyshev_iris  r  z#TestPdist.test_pdist_chebyshev_irisc                 C   r  )Nr   r  r
  r   r   r  r  r  rZ   rZ   r[   !test_pdist_chebyshev_iris_float32  r  z+TestPdist.test_pdist_chebyshev_iris_float32c                 C   ry  )Nr  r  r
  r  r}  r  r  rZ   rZ   r[   test_pdist_chebyshev_iris_nonC  r  z(TestPdist.test_pdist_chebyshev_iris_nonCc                 C   h   t tg dtg d}t tjg dtdtjg dtd}t|dddd t|dddd d S Nrh   r   rh   rh   r   rh   rh   r   rh   rh   r?   333333?r   绽|=r  r   	wmatchingrQ   r   r   r
   r   mZm2rZ   rZ   r[   test_pdist_matching_mtica1     z$TestPdist.test_pdist_matching_mtica1c                 C   r  Nrh   r   rh   rh   rh   r   r?   UUUUUU?r   r  r  r  r  rZ   rZ   r[   test_pdist_matching_mtica2  r  z$TestPdist.test_pdist_matching_mtica2c                 C   r  r  wjaccardrQ   r   r   r
   r  rZ   rZ   r[   test_pdist_jaccard_mtica1     z#TestPdist.test_pdist_jaccard_mtica1c                 C   r  r  r   r  rZ   rZ   r[   test_pdist_jaccard_mtica2  r#  z#TestPdist.test_pdist_jaccard_mtica2c                 C   x   t tg dtg d}t tjg dtdtjg dtd}tdkr*t| t|dddd t|dddd d S )Nr  r  r?   r   r   r  r  wyulerQ   r   r   r   rE  r
   r  rZ   rZ   r[   test_pdist_yule_mtica1     z TestPdist.test_pdist_yule_mtica1c                 C   r%  )Nr  r  r?   r   r   r  r  r&  r  rZ   rZ   r[   test_pdist_yule_mtica2  r)  z TestPdist.test_pdist_yule_mtica2c                 C   x   t tg dtg d}t tjg dtdtjg dtd}tdkr*t| t|dddd t|dddd d S )	Nr  r  r?   r   g۶m۶m?r   r  r  wdicerQ   r   r   r   rE  r
   r  rZ   rZ   r[   test_pdist_dice_mtica1  r)  z TestPdist.test_pdist_dice_mtica1c                 C   r+  )	Nr  r  r?   r   r)  r   r  r  r,  r  rZ   rZ   r[   test_pdist_dice_mtica2  r)  z TestPdist.test_pdist_dice_mtica2c                 C   r+  	Nr  r  r?   r   g      ?r   r  r  )r,   rQ   r   r   r   rE  r
   r  rZ   rZ   r[   test_pdist_sokalsneath_mtica1  r)  z'TestPdist.test_pdist_sokalsneath_mtica1c                 C   r+  	Nr  r  r?   r   g?r   r  r  )wsokalsneathrQ   r   r   r   rE  r
   r  rZ   rZ   r[   test_pdist_sokalsneath_mtica2   r)  z'TestPdist.test_pdist_sokalsneath_mtica2c                 C   r+  r0  wrogerstanimotorQ   r   r   r   rE  r
   r  rZ   rZ   r[    test_pdist_rogerstanimoto_mtica1*  r)  z*TestPdist.test_pdist_rogerstanimoto_mtica1c                 C   r+  r2  r5  r  rZ   rZ   r[    test_pdist_rogerstanimoto_mtica24  r)  z*TestPdist.test_pdist_rogerstanimoto_mtica2c                 C   r+  )	Nr  r  r?   r   r  r   r  r  wrussellraorQ   r   r   r   rE  r
   r  rZ   rZ   r[   test_pdist_russellrao_mtica1>  r)  z&TestPdist.test_pdist_russellrao_mtica1c                 C   r+  )	Nr  r  r?   r   r  r   r  r  r9  r  rZ   rZ   r[   test_pdist_russellrao_mtica2H  r)  z&TestPdist.test_pdist_russellrao_mtica2c                 C   sN   t d }tdkrt|j|j d}t|d}t|d}t|||tdkd d S )Nr  r   r  r   test_canberrar  )rR   r   rE  r   r@   r  r
   )r   Dr  rF  rG  rZ   rZ   r[   test_pdist_canberra_matchR  s   

z#TestPdist.test_pdist_canberra_matchc                 C   s2   d}t dgdgfd}d}t|||tdkd d S )Nr  r   r   r   gd>@1?r   r  )r  r
   r   )r   r  pdist_yr@  rZ   rZ   r[   test_pdist_canberra_ticket_711\  s   z(TestPdist.test_pdist_canberra_ticket_711c                    sZ   G dd d  fdd}t j  g  ggtd}t||d}d}t||tdkd	 d S )
Nc                   @   r0  )z6TestPdist.test_pdist_custom_notdouble.<locals>.myclassNr1  rZ   rZ   rZ   r[   r4  f  r5  r4  c                    r6  r7  r9  r;  r<  rZ   r[   r   i  r=  z9TestPdist.test_pdist_custom_notdouble.<locals>._my_metricr?   r&  r8  r   r>  )rQ   r   r?  r   r   r   )r   r   rB   r@  r@  rZ   r<  r[   test_pdist_custom_notdoubled  s   z%TestPdist.test_pdist_custom_notdoublerB  c           
   
   K   s  z$t |fd|i|}t |fdt|i|}t |fdd| i|}W nK tyo } z?|j}	tdkr<t|	j t| t|	t |fd|i| t|	t |fdt|i| t|	t |fdd| i| W Y d }~d S d }~ww t|||tdkd t|||tdkd d S rC  )	r   r  r   rD  r   rE  r   r  r
   )
r   r  r   r  r   rF  rG  rH  r   rI  rZ   rZ   r[   rJ  r  s    
*	z$TestPdist._check_calling_conventionsc              	   C   s   | j D ]k}t| d d dd d df }tD ]X}tdkr"td|d| |dv r+d|vr+q| || |dkrJtj|tj	dd	d
}| j|||d q|dkrmt
t|tj	j}ttj|j}| j|||d qqd S )Nr%  r   rK  rL  rM  r   r*   r   rh   rN  rO  r%   rQ  )r   rR   r   r   rE  rJ  rQ   rT  r   r>   rU  rV  r#  r   r  rW  )r   rX  r  r   rP  rR  rZ   rZ   r[   test_pdist_calling_conventions  s(   
z(TestPdist.test_pdist_calling_conventionsc                 C   s*  d}t d | jd ft d | jd ft d | jd ft d | jd	 fg}tD ]l}|D ]g}|d
 d d dd d df }zt||d}W n5 tyw } z)|j}tdkrZt|j t| |d D ]}	|	|}
t	|t|
|d q^W Y d }~q*d }~ww |d D ]}	t|	||d}t
|||tdkd q|q*q&d S )NrB  rD   r   rH   rV   rG   r   rF   rT   r   r%  r   r&  rh   r  )rR   r   r   r   r   rD  r   rE  r   r  r
   )r   r  rZ  r   r[  r  rF  r   rI  r\  r  rG  rZ   rZ   r[   test_pdist_dtype_equivalence  s8   
	z&TestPdist.test_pdist_dtype_equivalencec                 C   sF  d}t d d d dd d df }t|jd |jd d  d }tD ]~}t }|dkr/d|d	< tj|tjd
}t||fi |}t||fd|i|}t	|||d t
||u  tj|d tjd
}	ttt||fd|	i| tjd| tjd
d d d }
ttt||fd|
i| tj|tjd
}ttt||fd|i| q"d S )Nr  rF   r%  r   r   rh   r'   r  r  r?   r^  r}  r   )rR   r   r   r   r`  rQ   ra  r>   r   r
   r   r  rs   rb  )r   r  r  Zout_sizer   r   rc  r  r  rd  re  rf  rZ   rZ   r[   test_pdist_out  s&    zTestPdist.test_pdist_outc                 C   s   d}t d d d dd d df }| }t|jj  t|jj tD ])}t }|dkr0d|d< t||fi |}t||fi |}t|||t	dkd q#d S )	Nr  rF   r%  r   r'   r  r  r  )
rR   rh  r   ri  rj  r   r`  r   r
   r   )r   r  r  ZX_copyr   r   r  r  rZ   rZ   r[   rl    s   zTestPdist.test_stridingNrr  )fr   r2  r3  r   rx  r  r  r  r  rs  rt  Zslowr  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  ru  rQ   rv  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/  r1  r4  r7  r8  r;  r<  r?  rA  rB  rJ  rC  rD  rE  rl  rZ   rZ   rZ   r[   rw    s    	$


		


	


	









		










	
rw  c                   @   L   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	dd Z
dS )TestSomeDistanceFunctionsc                 C   s,   t g d}t g d}||fg| _d S Nr   r   r   r   )rQ   r   cases)r   r   r   rZ   rZ   r[   r     s   z&TestSomeDistanceFunctions.setup_methodc           	      C   s   | j D ]4\}}t||dd}t|d t||dd}t|d t||dd}t|d t||dd}t|d	 qtd
dg}tddg}tt||t|d|d d S )Nrh   r  r         ?gqaz@r   gw@r  gR"6@i`  i  i^  i  uint16)rJ  r'   r	   rQ   r   r   r   )	r   r   r   dist1Zdist1p5dist2Zdist0p25rd   r   rZ   rZ   r[   r    s   



z(TestSomeDistanceFunctions.test_minkowskic                 C   s.   | j D ]\}}t||}t|td qd S )Nr%  )rJ  
weuclideanr	   rQ   r*  r   r   r   r-  rZ   rZ   r[   r    s   
z(TestSomeDistanceFunctions.test_euclideanc                 C   s(   | j D ]\}}t||}t|d qd S )Nr   )rJ  wsqeuclideanr	   rP  rZ   rZ   r[   test_sqeuclidean	  s   
z*TestSomeDistanceFunctions.test_sqeuclideanc              	   C   s@   | j D ]\}}t||}t|ddtdtd    qd S )Nr   g      2@      )rJ  wcosiner	   rQ   r*  rP  rZ   rZ   r[   r    s   
$z%TestSomeDistanceFunctions.test_cosinec                 C   s`   t g d}t g d}| jD ]\}}t||}t|dt ||t|t|    qd S )N)      r   r   )UUUUUUrW  gUUUUUU@r   )rQ   r   rJ  wcorrelationr	   r"  r   )r   ZxmZymr   r   r-  rZ   rZ   r[   r    s   
(z*TestSomeDistanceFunctions.test_correlationc                 C   sT   t g d}t g d}t||}d|  kr%dt t jj ks(J  J d S )N)%r'  r'  r'  r'  r'  r'         r'  r'  r'  rY  rY  rY  r'  rY  r'  rY  r'  r'  rV  rY  r'  r   r'  r'  rY  r'  r'  rY  r'  rY  rY  rY  rY  rY  rY  r'  )%r   r   r   r   r   r   rV  r   r   r   rV  rV  rV  r   rV  r   rV  r   r   r'  rV  r   r   r   r   rV  r   r   rV  r   rV  rV  rV  rV  rV  rV  r   r   rm  )rQ   r   r   Zfinfor   r  rP  rZ   rZ   r[   test_correlation_positive  s   
.z3TestSomeDistanceFunctions.test_correlation_positivec                 C   sh   t g d}t g d}t g dg dg dg}| jD ]\}}t|||}t|t d qd S )Nr   rI  )r   r   r'  )r   r   r   )r'  r   r   g      @)rQ   r   rJ  r%   r	   r*  )r   r   r   vir-  rZ   rZ   r[   test_mahalanobis%  s   z*TestSomeDistanceFunctions.test_mahalanobisN)r   r2  r3  r   r  r  rR  r  r  rZ  r\  rZ   rZ   rZ   r[   rG    s    rG  c                     sZ   dd } t  td t  fddt D }| |d< d}t| f|t |< d S )Nc                 S   sh   t g d}t g d}|d d t jf }|d d t jf }|j}|j}||f||f||fg| _d S rH  )rQ   r   newaxisr#  rJ  )r   r   r   Zx31Zy31Zx13Zy13rZ   rZ   r[   r   1  s   z-construct_squeeze_tests.<locals>.setup_methodz?.*distance metrics ignoring length-1 dimensions is deprecated.*c                    s&   i | ]}| d r|t |qS )r   )
startswithgetattr)rc   rX   basesuprZ   r[   
<dictcomp>B  s    z+construct_squeeze_tests.<locals>.<dictcomp>r   ZTestDistanceFunctionsSqueeze)r   filterDeprecationWarningrG  dirr   globals)r   ZattrsrX   rZ   r`  r[   construct_squeeze_tests.  s   rh  c                   @   sR   e Zd ZejejejejegZ	dd Z
dd Zdd Zdd Zd	d
 Zdd ZdS )TestSquareFormc                 C      | j D ]}| | qd S rj   )checked_dtypescheck_squareform_matrixr   r@   rZ   rZ   r[   test_squareform_matrixR     
z%TestSquareForm.test_squareform_matrixc                 C   rj  rj   )rk  check_squareform_vectorrm  rZ   rZ   r[   test_squareform_vectorV  ro  z%TestSquareForm.test_squareform_vectorc                 C   s   t jd|d}t|}t|jd t|j| t jd|d}t|}t|jd t|j| t jddgddgg|d}t|}t|jd t|j| t|t jdg|d d S )N)r   r   r?   r   rh   rh   r   g@rg   )rQ   r  r   r   r   r@   r   r   )r   r@   AZrArZ   rZ   r[   rl  Z  s   z&TestSquareForm.check_squareform_matrixc                 C   s   t jd|d}t|}t|jd t|j| t|dgg t jdg|d}t|}t|jd t|j| t|t jddgddgg|d d S )Nrr  r?   rs  r   g @)r   r   )rQ   r  r   r   r   r@   r   r   )r   r@   rp  ZrvrZ   rZ   r[   rp  k  s   $z&TestSquareForm.check_squareform_vectorc                 C      t ddD ]}| | qd S )Nr   r%  )rt   check_squareform_multi_matrixr   nrZ   rZ   r[   test_squareform_multi_matrixx     z+TestSquareForm.test_squareform_multi_matrixc           
      C   s   t j|d}t|}tt|jd t|}t|}|j}d}tdkr-t	|j|j|j tt|d tt|jd t|d |d  t
d|d D ]+}t
|d |d D ]}	||	krmt|||	f ||  |d7 }qWt|||	f d qWqLd S )N   rh   r   r   r   )rQ   r   r   r  r   lenr   r   r   rE  rt   )
r   rx  r  r  rt  ZYrskr   jrZ   rZ   r[   rv  |  s(   
z,TestSquareForm.check_squareform_multi_matrixN)r   r2  r3  rQ   r   rT   int32int8r   rk  rn  rq  rl  rp  ry  rv  rZ   rZ   rZ   r[   ri  O  s    ri  c                   @   d   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	dd Z
dd Zdd Zdd ZdS )TestNumObsYc                 C   s8   t ddD ]}tj|d}t|}tt|| qd S )Nr   rm  r{  )rt   rQ   r   r   r  r   r   )r   rx  r  r  rZ   rZ   r[   test_num_obs_y_multi_matrix  s
   z'TestNumObsY.test_num_obs_y_multi_matrixc                 C   s   t t| jd d S Nrh   )r  rs   check_yr   rZ   rZ   r[   test_num_obs_y_1     zTestNumObsY.test_num_obs_y_1c                 C      t | d d S Nr   r   r  r   rZ   rZ   r[   test_num_obs_y_2  r  zTestNumObsY.test_num_obs_y_2c                 C   r  )Nr   r  r   rZ   rZ   r[   test_num_obs_y_3     zTestNumObsY.test_num_obs_y_3c                 C   r  Nr{  r  r   rZ   rZ   r[   test_num_obs_y_4  r  zTestNumObsY.test_num_obs_y_4c                 C   ru  )Nr%     )rt   minit)r   r   rZ   rZ   r[   test_num_obs_y_5_10  rz  zTestNumObsY.test_num_obs_y_5_10c                 C   X   t g }tddD ]}|||d  d  q	tddD ]}||vr)tt| j| qd S Nr   r  rh   r%  i   setrt   addr  rs   bad_yr   rd   rx  r   rZ   rZ   r[   test_num_obs_y_2_100  s   z TestNumObsY.test_num_obs_y_2_100c                 C   s   t | | d S rj   r  rw  rZ   rZ   r[   r    r  zTestNumObsY.minitc                 C   s   t j|}t|S rj   )rQ   r   r   r   r   rx  r   rZ   rZ   r[   r    s   zTestNumObsY.bad_yc                 C      t | ||kS rj   )r   make_yrw  rZ   rZ   r[   r    r  zTestNumObsY.check_yc                 C   s   t j||d  d S Nrh   r   rQ   r   r   rw  rZ   rZ   r[   r    r   zTestNumObsY.make_yN)r   r2  r3  r  r  r  r  r  r  r  r  r  r  r  rZ   rZ   rZ   r[   r    s    
r  c                   @   rF  )TestNumObsDMc                 C   sV   t ddD ]#}tj|d}t|}t|}tdkr!t|j|j t	t
|| qd S )Nrh   rm  r{  r   )rt   rQ   r   r   r  r   r   rE  r   r   r   )r   rx  r  r  rt  rZ   rZ   r[   test_num_obs_dm_multi_matrix  s   z)TestNumObsDM.test_num_obs_dm_multi_matrixc                 C   r  )Nr   r   check_Dr   rZ   rZ   r[   test_num_obs_dm_0     zTestNumObsDM.test_num_obs_dm_0c                 C   r  r  r  r   rZ   rZ   r[   test_num_obs_dm_1  r  zTestNumObsDM.test_num_obs_dm_1c                 C   r  r  r  r   rZ   rZ   r[   test_num_obs_dm_2  r  zTestNumObsDM.test_num_obs_dm_2c                 C   r  r  r  r   rZ   rZ   r[   test_num_obs_dm_3  r  zTestNumObsDM.test_num_obs_dm_3c                 C   r  r  r  r   rZ   rZ   r[   test_num_obs_dm_4  r  zTestNumObsDM.test_num_obs_dm_4c                 C   r  rj   )r   make_Drw  rZ   rZ   r[   r    r  zTestNumObsDM.check_Dc                 C   s   t j||S rj   r  rw  rZ   rZ   r[   r    s   zTestNumObsDM.make_DN)r   r2  r3  r  r  r  r  r  r  r  r  rZ   rZ   rZ   r[   r    s    	r  c                 C      t | ddS NT)throw)r   )r>  rZ   rZ   r[   is_valid_dm_throw  r   r  c                   @   st   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	dd Z
dd Zdd Zdd Zdd Zdd ZdS )TestIsValidDMc                 C       t jdt jd}ttt| d S )Nr%  r?   rQ   r  r>   r  rs   r  r   r>  rZ   rZ   r[   $test_is_valid_dm_improper_shape_1D_E     z2TestIsValidDM.test_is_valid_dm_improper_shape_1D_Ec                 C   "   t jdt jd}tt|d d S )Nr  r?   FrQ   r  r>   r   r   r  rZ   rZ   r[   $test_is_valid_dm_improper_shape_1D_F     z2TestIsValidDM.test_is_valid_dm_improper_shape_1D_Fc                 C   r  Nr   r   r   r?   r  r  rZ   rZ   r[   $test_is_valid_dm_improper_shape_3D_E  r  z2TestIsValidDM.test_is_valid_dm_improper_shape_3D_Ec                 C   r  Nr  r?   Fr  r  rZ   rZ   r[   $test_is_valid_dm_improper_shape_3D_F  r  z2TestIsValidDM.test_is_valid_dm_improper_shape_3D_Fc                 C   s@   t jd}t|}tddD ]}d|||f< qttt| d S )Nrm  r   r%  r   )rQ   r   r   r   rt   r  rs   r  r   r   r>  r   rZ   rZ   r[   #test_is_valid_dm_nonzero_diagonal_E  s
   z1TestIsValidDM.test_is_valid_dm_nonzero_diagonal_Ec                 C   sB   t jd}t|}tddD ]}d|||f< qtt|d d S )Nrm  r   r%  r   F)rQ   r   r   r   rt   r   r   r  rZ   rZ   r[   #test_is_valid_dm_nonzero_diagonal_F  s
   z1TestIsValidDM.test_is_valid_dm_nonzero_diagonal_Fc                 C   s4   t jd}t|}|d d |d< ttt| d S )Nrm  r   rh   rh   rh   r   )rQ   r   r   r   r  rs   r  r   r   r>  rZ   rZ   r[   test_is_valid_dm_asymmetric_E  s   z+TestIsValidDM.test_is_valid_dm_asymmetric_Ec                 C   s6   t jd}t|}|d d |d< tt|d d S )Nrm  r  rh   r  FrQ   r   r   r   r   r   r  rZ   rZ   r[   test_is_valid_dm_asymmetric_F  s   z+TestIsValidDM.test_is_valid_dm_asymmetric_Fc                 C   r  )Nrs  r?   Tr  r  rZ   rZ   r[   test_is_valid_dm_correct_1_by_1  r  z-TestIsValidDM.test_is_valid_dm_correct_1_by_1c                 C   &   t jd}t|}tt|d d S )Nrh   Tr  r  rZ   rZ   r[   test_is_valid_dm_correct_2_by_2     z-TestIsValidDM.test_is_valid_dm_correct_2_by_2c                 C   r  Nr   Tr  r  rZ   rZ   r[   test_is_valid_dm_correct_3_by_3"  r  z-TestIsValidDM.test_is_valid_dm_correct_3_by_3c                 C   r  )N   Tr  r  rZ   rZ   r[   test_is_valid_dm_correct_4_by_4'  r  z-TestIsValidDM.test_is_valid_dm_correct_4_by_4c                 C   r  )Nrm  Tr  r  rZ   rZ   r[   test_is_valid_dm_correct_5_by_5,  r  z-TestIsValidDM.test_is_valid_dm_correct_5_by_5N)r   r2  r3  r  r  r  r  r  r  r  r  r  r  r  r  r  rZ   rZ   rZ   r[   r    s    r  c                 C   r  r  )r   )r   rZ   rZ   r[   is_valid_y_throw2  r   r  c                   @   r  )TestIsValidYc                 C   r  )Nr   r   r?   rQ   r  r>   r  rs   r  r   r   rZ   rZ   r[   #test_is_valid_y_improper_shape_2D_E;  r  z0TestIsValidY.test_is_valid_y_improper_shape_2D_Ec                 C   r  )Nr  r?   FrQ   r  r>   r   r   r  rZ   rZ   r[   #test_is_valid_y_improper_shape_2D_F?  r  z0TestIsValidY.test_is_valid_y_improper_shape_2D_Fc                 C   r  r  r  r  rZ   rZ   r[   #test_is_valid_y_improper_shape_3D_EC  r  z0TestIsValidY.test_is_valid_y_improper_shape_3D_Ec                 C   r  r  r  r  rZ   rZ   r[   #test_is_valid_y_improper_shape_3D_FG  r  z0TestIsValidY.test_is_valid_y_improper_shape_3D_Fc                 C      |  d}tt|d d S )Nr   Tcorrect_n_by_nr   r   r  rZ   rZ   r[   test_is_valid_y_correct_2_by_2K     
z+TestIsValidY.test_is_valid_y_correct_2_by_2c                 C   r  r  r  r  rZ   rZ   r[   test_is_valid_y_correct_3_by_3O  r  z+TestIsValidY.test_is_valid_y_correct_3_by_3c                 C   r  )Nr{  Tr  r  rZ   rZ   r[   test_is_valid_y_correct_4_by_4S  r  z+TestIsValidY.test_is_valid_y_correct_4_by_4c                 C   r  )Nr%  Tr  r  rZ   rZ   r[   test_is_valid_y_correct_5_by_5W  r  z+TestIsValidY.test_is_valid_y_correct_5_by_5c                 C   r  r  r  r  rZ   rZ   r[   test_is_valid_y_2_100[  s   z"TestIsValidY.test_is_valid_y_2_100c                 C   s   t j|}t|ddS r  )rQ   r   r   r   r  rZ   rZ   r[   r  c  s   zTestIsValidY.bad_yc                 C   s   t j||d  d }|S r  r  r  rZ   rZ   r[   r  g  s   zTestIsValidY.correct_n_by_nN)r   r2  r3  r  r  r  r  r  r  r  r  r  r  r  rZ   rZ   rZ   r[   r  6  s    r  r  )g      $g      r'  c                 C   s:   t ttddgddg|  t ttddgddg| ddg d S )Nrh   r   r   r{  )r  rs   r'   r  rZ   rZ   r[   
test_bad_pl  s   "r  c                   C   s   t ttg dg d d S )N)FFF)r  rs   r,   rZ   rZ   rZ   r[   test_sokalsneath_all_falses  s   r  c                   C   s4   t tg dg dd t tg dg dd d S )Nrh   r   r   r   r{  r  rh   rh   rh   r   r   rh   r   rh   r   r   )r   	wcanberrarZ   rZ   rZ   r[   r=  x  s   r=  c                   C   s<   t tg dg dddd t tg dg dddd d S )	Nr  r  gUUUUUU?   Zdecimalr  r  r)  )r	   wbraycurtisrZ   rZ   rZ   r[   test_braycurtis~  s    r  c                  C   s  t g d} t g d}tt| |ddd tt| |t ddd tjtdd# tt| t j	d d f |t j	d d f t ddd W d    n1 sRw   Y  tjtdd  tt| t j	d d f |t j	d d f ddd W d    n1 sw   Y  tjtdd  tt| d d t j	f |d d t j	f ddd W d    n1 sw   Y  t 
d	d
d
}ttt|| ttt|| t jd}|d}|d}t||}t||}t|d
 |dd d S )N)rh   rh   rh   )r   r   r   r   rS  r  r   *ignoring length-1 dimensions is deprecatedmatchr{  r   iIrm  )rQ   r   r	   rQ  rO  r*  rs  warnsre  r]  r  r   r  rs   r   r   r   )r+  r,  r   Zrsr   d1d2rZ   rZ   r[   test_euclideans  sD   &
&&



r  c                  C   s"   g d} g d}t tt| | d S )N)r   r   rh   r  )r  rs   whammingr;  rZ   rZ   r[   test_hamming_unequal_length  s   r  c                  C   s<   t jg ddd} t jg ddd}d}tt| || d S )N)eggsspamr  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  z|S4r?   )r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  g?)rQ   r   r
   r  )rd   r   ZdesiredrZ   rZ   r[   test_hamming_string_array  s   

r  c                  C   sv   t g dg dg} t| ddd d}t| | ddd d}t| ddd}t| | ddd}t||dd t||dd d S )N)
QUUT@      Y@r  r  g      B@g      N@g     V@g     b@g      8@g      H@r'   rh   )r   r  r   )r   r  r  r}  )rQ   r   r   r   r
   )Zarr_inZp0Zc0Zp1Zc1rZ   rZ   r[   test_minkowski_w  s   
r  c                  C   s<  g d} g d}t jt jt jt jfD ]}tt j| |dt j||d}tt |j	t j
 qt jt jt jt jfD ]-}tdgt jdg|d}tt jdg|ddg}t|| t|t t |jd  q8t jt jt jt jg}dD ]}tt |r|tt | qr|D ]}tt j| |dt j||d}t|j	| qd S )Nr  )r{  r%  r  r?   r   rx   r   )Zfloat16Zfloat128)rQ   r  Zint16r  rb  rQ  asarrayr   Z
issubdtyper@   ZfloatingZuint8rL  Zuint32Zuint64r   r   ZiinfomaxrT   Z	complex64Z
complex128hasattrr   r_  )r   r   r@   dr  r  ZdtypesrZ   rZ   r[   test_sqeuclidean_dtypes  s&   

r  c                  C   sN   g d} g d}dd | D }dd |D }t | |}t ||}t|| d S )N)TTF)TFTc                 S      g | ]}t |qS rZ   r   rc   r   rZ   rZ   r[   re         z&test_sokalmichener.<locals>.<listcomp>c                 S   r  rZ   r  r  rZ   rZ   r[   re     r  )r+   r   )r  r   r   r   rM  rN  rZ   rZ   r[   test_sokalmichener  s   

r	  c            	      C   s   d} d}d}d}d||   || d||     }t |d tddgddgddgd}t || g d	}g d
}dD ]}t t|||gd q9d S )Nr'  g?r   r   gX$I?rh   r   r   )FFTTTFFTTTTTTFTFFFTT)TTTFFTTTFTTTTTFFFTTT)g?r  r   g      4@r  )r	   r+   )	ZntfZnftZnttZnffZexpectedZactualZa1Za2r   rZ   rZ   r[   test_sokalmichener_with_weight  s    

r  c                  C   sX   t g dg dg dg dg} |  }tD ]}t| | | t| | t| | qd S )Nr   r   r   g3333336@gL7@g333333F@)rQ   r   rh  r   r   r   r   )r  rk  r   rZ   rZ   r[   test_modifies_input  s   

r  c               
   C   sf  t g dg dg dg dg} t d}tD ]}tt t| | |d W d    n1 s1w   Y  tt t| |d W d    n1 sLw   Y  dD ]\}|di}|d	kra|d
ksq|dkri|dksq|dkrr|dkrrqStt t| | |fi | W d    n1 sw   Y  tt t| |fi | W d    n1 sw   Y  qSqd S )Nr   r   r   r  r   r   )r  rP  rR  ZfoorP  r*   rR  r%   r  r'   )	rQ   r   r  r   rs  r   r   r   r   )r  r   r   r   r   rZ   rZ   r[   test_Xdist_deprecated_args  sD   

r  c               
   C   s   t d d d dd d df } t| jd }|d d d  |d d d< tD ]&}|dv r,q%|t|d| fD ]}ttt| ||d ttt	| | ||d q5q%d S )NrF   r%  r   rh   )r*   r%   r"   r   r
  )
rR   rQ   r   r   r   r  r  rs   r   r   )r  r   r   r  rZ   rZ   r[   test_Xdist_non_negative_weights8  s   r  c                  C   sD  g d} t | }t||  t | tjd}t||  t|jtj dg} t | }t|jd t||  d} tjt	dd t | }W d    n1 sJw   Y  t|jd t|| g t
dddd} tjt	dd t | }W d    n1 s{w   Y  t|jd t|| dd d df  dd	gd
dgg} ttt |  d S )Nr  r?   rh   r  r  r%  rx   r   r   r   r{  )r   r   rQ   r   r   r@   ri   rs  r  re  r  r   r  rs   r;  rZ   rZ   r[   test__validate_vectorE  s8   




r  c                  C   sp   t jdtd} t| d | d }|dksJ t| d}t|dg t| d d | d d d}t|dgg d S )N)r   r  r?   r   r'  r.   rh   )rQ   r   r   r'  r   r   r   )r   r  rZ   rZ   r[   test_yule_all_samed  s   
r  c                  C   sJ  t tg dg ddd t tddgddgd t tg dg dd t tddggddggd	d
ddg t tddggddggdd
dg t tddggddggd	ddddgg t tddggddggddddgg tg dg dg dg} tg dg dg dg}t t| |d	d
g d t t| |dd
g d d S )N)r   r'  r'  )r'  r   r'  r   r   r'  r)  g)>d?rK  r   rw   rh   gД?T)ru   Zkeepdims)rh   r   r   r{  )r%  r     r]   )	   rm        )   rS  r  r  )            )            )g'Y?gTyCЇ?gFw$?g\BZ?)g;0%7/?g)o?gku?)r	   r"   rQ   r   )rd   r   rZ   rZ   r[   r  q  sJ   

r  rj   )	NNFTFFFFN)~os.pathrK   	functoolsr   r   rn  ZnumpyrQ   r   Znumpy.linalgr   Znumpy.testingr   r   r   r   r	   r
   r   rs  r   r  Zscipy.spatial.distanceZscipyZscipy.spatialr   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.   rI   r   Z_tdistZ_ytdistrR   r\   r_   rv   r   r   r   r   r   r   r!  r  r  r  r6  r  r  r'  r-  Z
wcityblockZ
wchebyshevrU  rX  Z
wkulsinskiZwkulczynski1r!  rO  rQ  r  r  r3  Zwsokalmichenerr:  r   rw  rG  rh  ri  r  r  r  r  r  r  rt  ru  r  r  r=  r  r  r  r  r  r  r	  r  r  r  r  r  r  r  rZ   rZ   rZ   r[   <module>   s   "$0d
 

0	




`        bFC4#E6
%