o
    EbT!                     @   s  d dl mZ d dlZd dlmZmZ d dlZd dlm	Z	m
Z
mZ d dlmZm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ejdg ddd Zejd g d!g d"g d"gg d!g d#g d#gg d"g d$gg d%g d&gd'd gd(d gd)d gggd*d+ Zd,d- Zd.d/ Zed0d'gg d1g d2g d3gg d4fg d5g d6g d7gg d8fg d9g d:g d;gg d<fg d=g d>g d?gg d@fdAe dBe dBge dBe dBd'ge dBdCe dBggg d<fgZ!ejdDe!dEdF Z"dS )G    )productN)assert_array_equalassert_equal)
csr_matrix
coo_matrixdiags)maximum_bipartite_matching"min_weight_full_bipartite_matchingc                  C   sN   t t tddgddgg} t|  W d    d S 1 s w   Y  d S )Nr      )pytestraises	TypeErrornparrayr   )graph r   J/usr/lib/python3/dist-packages/scipy/sparse/csgraph/tests/test_matching.py5test_maximum_bipartite_matching_raises_on_dense_input   s   
"r   c                  C   sB   t d} t| dd}t| dd}tg }t|| t|| d S )Nr   r   rowZ	perm_typecolumnr   r   r   r   r   r   xyexpected_matchingr   r   r   +test_maximum_bipartite_matching_empty_graph   s   

r   c                  C   sH   t d} t| dd}t| dd}ttg | ttddg| d S )N   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   4test_maximum_bipartite_matching_empty_left_partition   s
   r#   c                  C   sH   t d} t| dd}t| dd}ttg d| ttg | d S )Nr      r   r   r   )r    r    r    r!   r"   r   r   r   5test_maximum_bipartite_matching_empty_right_partition$   s
   r&   c                  C   sL   t d} t| dd}t| dd}ttddg| ttddg| d S )Nr   r   r   r   r   r    r!   r"   r   r   r   3test_maximum_bipartite_matching_graph_with_no_edges,   s
   r(   c                  C   sR   t ddgddgg} t| dd}t| dd}tddg}t|| t|| d S )Nr
   r   r   r   r   r   r   r   r   r   >test_maximum_bipartite_matching_graph_that_causes_augmentation4   s   
r)   c                  C   s^   t ddgddgddgg} t| dd}t| dd}ttg d| ttddg| d S )Nr
   r   r   r   r   r   r    r
   r   r!   r"   r   r   r   Atest_maximum_bipartite_matching_graph_with_more_rows_than_columns?   s
   r+   c                  C   sX   t g dg dg} t| dd}t| dd}ttddg| ttg d| d S )	N)r
   r
   r   r   r   r
   r   r   r   r   r   r*   r!   r"   r   r   r   Atest_maximum_bipartite_matching_graph_with_more_columns_than_rowsG   s
   r-   c                  C   sh   ddg} ddg}g d}t | ||fdd}t|dd}t|dd}tddg}t|| t|| d S )	Nr   r
   )r   r
   r   r'   shaper   r   r   r   )dataindicesindptrr   r   r   r   r   r   r   =test_maximum_bipartite_matching_explicit_zeros_count_as_edgesO   s   
r3   c                  C   s   t jdtd} g d}g d}t| ||fdd}t|dd}t|d	d}|d
k dks.J |d
k dks8J tt|jd |D ]\}}|d
krR|||f sRJ qBt|t|jd D ]\}}|d
krm|||f smJ q]d S )N2   Zdtype)2                  r9   r%      
   r;      r6   r7      r;   r?         r9   r%      r%   r?   r@   r6   r7   r8   r9   r:   r;   r9   r%   r<   r=   r;   r>   r6   r7   r?   r;   r?   r@   rA   r9   r%   rB   r%   r?   r@   )r   r;      r=   r=   rB   rA   r9   r9   r:             #   rG   (   -   /   rJ   0   r4   )rA   rD   r.   r   r   r   r    r?   r   r
   )	r   onesintr   r   sumzipranger/   )r0   r1   r2   r   r   r   uvr   r   r   5test_maximum_bipartite_matching_feasibility_of_result[   s"   rS   c                  C   s\  t jd tt dddd} t jd}t jd}t d}|}t jdtd}t|||ff	 }|}t d}t jdtd}	t|	||ff	 }
||  |
 }t
|dd}t d}|}t jdtd}t|||ff	 }|| }t
|d	d}|}t d}t jdtd}	t|	||ff	 }
||
 }tt| dkd
 tt| dkd
 d S )N*   rD   r   Zcsr)Zoffsetsformatr5   r   r   r   F)r   ZrandomZseedr   rL   ZpermutationZarangerM   r   Ztocsrr   r   anyZdiagonal)AZ	rand_permZ
rand_perm2ZRrowZRcolZRdataZRmatZCrowZCcolZCdataZCmatBZpermZC1Zperm2ZC2r   r   r   Ftest_matching_large_random_graph_with_one_edge_incident_to_each_vertexq   s6   



rY   znum_rows,num_cols)r   r   r$   c                 C   s<   t || f}t|\}}t|dksJ t|dksJ d S )Nr   )r   r	   len)Znum_rowsZnum_colsbiadjacency_matrixrow_indcol_indr   r   r   +test_min_weight_full_matching_trivial_graph   s   r^   r[   )r
   r
   r
   )r
   r   r   r,   )r   r   r   )r   r
   r   )r   r   r   r
   r   r;   c                 C   s<   t t tt|  W d    d S 1 sw   Y  d S N)r   r   
ValueErrorr	   r   r[   r   r   r   1test_min_weight_full_matching_infeasible_problems   s   	"rb   c                  C   s@   t t td} t|  W d    d S 1 sw   Y  d S )N))r   r   r%   )r   r
   r
   )r   r   r%   )r   ZwarnsUserWarningr   r	   ra   r   r   r   !test_explicit_zero_causes_warning   s   
"rd   c           	   	   C   s   |\}}|dk}||| }|t | }| ||d\}}t|t | t|t |||f   |j}| ||d\}}t|t | tt |t t |||f   d S )Nr    )maximize)r   r   r   sortZflattenT)	ZsolverZ
array_typesign	test_caseZcost_matrixZexpected_costre   r\   r]   r   r   r    linear_sum_assignment_assertions   s$   

rj   r    )     rk   )rk     X  ),     ro   )rl   rk   ro   )rk   rl   rk   r
   )rk   rm   rn   r   )ro   rp   ro   r%   )rl   r   ro   )r=   r=   r<   )	   r<   r
   )rq   rC      )r=   r
   rC   )r=   r=   r<   r6   )rq   r<   r
   r
   )rq   rC   rr   r=   )r=   r
   rr   r=   infrC   zsign,test_casec                 C   s   t tt| | d S r_   )rj   r	   r   )rh   ri   r   r   r   *test_min_weight_full_matching_small_inputs   s   rt   )#	itertoolsr   Znumpyr   Znumpy.testingr   r   r   Zscipy.sparser   r   r   Zscipy.sparse.csgraphr   r	   r   r   r#   r&   r(   r)   r+   r-   r3   rS   rY   ZmarkZparametrizer^   rb   rd   rj   floatZ linear_sum_assignment_test_casesrt   r   r   r   r   <module>   s    	'

"