o
    6a                    @   s  d Z ddlZddlZddlZddlZddlmZmZ ddl	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 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&m'Z'm(Z(m)Z)m*Z*m+Z+m,Z,m-Z-m.Z.m/Z/m0Z0m1Z1m2Z2m3Z3m4Z4m5Z5m6Z6m7Z7m8Z8m9Z9m:Z:m;Z; G dd dZ<G d	d
 d
Z=G dd dZ>G dd dZ?G dd dZ@G dd dZAG dd dZBG dd dZCG dd dZDG dd dZEG dd dZFG dd dZGG dd  d ZHG d!d" d"ZIdS )#zTests suite for MaskedArray.
Adapted from the original test_ma by Pierre Gerard-Marchant

:author: Pierre Gerard-Marchant
:contact: pierregm_at_uga_dot_edu
:version: $Id: test_extras.py 3473 2007-10-29 15:18:13Z jarrod.millman $

    N)assert_warnssuppress_warnings)assert_assert_array_equalassert_equalassert_almost_equal)arrayarangemaskedMaskedArraymasked_arraygetmaskarrayshapenomaskoneszeroscount)!
atleast_1d
atleast_2d
atleast_3dmr_dotpolyfitcovcorrcoefmedianaverageuniquesetxor1d	setdiff1dunion1dintersect1din1dediff1dapply_over_axesapply_along_axiscompress_ndcompress_rowcolsmask_rowcolsclump_maskedclump_unmaskedflatnotmasked_contiguousnotmasked_contiguousnotmasked_edges
masked_allmasked_all_likeisindiagflatstackvstackc                   @   sL   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 )TestGenericc                 C   sJ  t dtd}tddgddgtd}t|| tddgddgd}t d|d}td	d	gd
d
g|d}t|| t d|d}td	d	gd	d	ggd
d
gd
d
gg|d}t|| tddddgfg}t d|d}tddgddg|d}t|| t d|d}tddgddg|d}t|| t d
|d}tdggdgg|d}t|| d S )N   dtype   maskr8   abfnamesformatsr   r   r9   r9   r6   r6   r<   r>   bar>   bbr>   r9   rC   )r.   floatr   r   npr8   )selftestcontroldt rQ   </usr/lib/python3/dist-packages/numpy/ma/tests/test_extras.pytest_masked_all%   s6   




zTestGeneric.test_masked_allc                 C   s   t ddtfgdffg}t jd|}tt|d t jjj tt|d d t jjj tt	|d d d t|d d j
d t|d d jj
d d S )Nr=   cr9   r9   rC   rQ   rL   r8   objectmar.   r   typecorer   lenr   _fill_valuerM   my_dtype
masked_arrrQ   rQ   rR   "test_masked_all_with_object_nestedC   s   z.TestGeneric.test_masked_all_with_object_nestedc                 C   sr   t dtdffg}t jd|}tt|d t jjj tt	|d d t|d j
d t|d jj
d d S )Nr=   rU   r9   rC   rQ   rV   r]   rQ   rQ   rR   test_masked_all_with_objectO   s   z'TestGeneric.test_masked_all_with_objectc                 C   s   t ddgtd}t|}t ddgddgtd}t|| tddgddgd}t d	d	gd
d
g|d}t|}t ddgd
d
g|d}t|| tddddgfg}t ddgddg|d}t|}t|| d S )Nr9   r6   r7   r:   r<   r=   r>   r?   rB   rC   )
   rb   rE   rF   rH   rJ   )r   rK   r/   r   rL   r8   )rM   baserN   rO   rP   rQ   rQ   rR   test_masked_all_likeY   s   

z TestGeneric.test_masked_all_likec           	      C   s   t ddD ]W}t d| D ]N}tj|td}tj||td}td| }|d| @ dk|_d}||D ]}||j|  7 }q3|t	krMt
|  | q|j |_t
|  | qqd S )Nr9      r6   r7   r   )rangerL   r	   intfullr   r;   datasumr*   r   
compressed)	rM   r>   ijkjar<   sslrQ   rQ   rR   check_clumpm   s   
zTestGeneric.check_clumpc                 C   sV   t td}t|g d< t|}tddtddtddg}t|| | t d S )Nrb   r   r9   r6         	   r      rt   re   ru   )r   rL   r	   r
   r)   slicer   rr   rM   r<   rN   rO   rQ   rQ   rR   test_clump_masked}   s   
zTestGeneric.test_clump_maskedc                 C   sN   t td}t|g d< t|}tddtddg}t|| | t d S )Nrb   rs   rw   rt   re   ru   )r   rL   r	   r
   r*   rx   r   rr   ry   rQ   rQ   rR   test_clump_unmasked   s   
zTestGeneric.test_clump_unmaskedc                 C   s   t d}t|}t|td|jg tjdtd|_t|td|jg t	||dk |dkB |dkB < t|}t|tddtddg t	|d d < t|}t|g  d S )	Nrb   r   r7   rw   ru      rt   rv   )
r	   r+   r   rx   sizerL   r   boolr;   r
   )rM   r<   rN   rQ   rQ   rR   test_flatnotmasked_contiguous   s   z)TestGeneric.test_flatnotmasked_contiguousN)__name__
__module____qualname__rS   r`   ra   rd   rr   rz   r{   r   rQ   rQ   rQ   rR   r4   #   s    

r4   c                   @   s<   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd ZdS )TestAveragec                 C   s$  t g dg dd}tdt|dd tdt|g dd t|g dd	d
\}}td| t|dk t|d d < tt|ddjd	g t g dg dd}|dd}t|d d df< tt|ddddg tt|ddjd d	g tddgt|dd t|dd	d\}}t|ddg d S )N)              ?       @      @)TFFFr;   r   r   axis)r   r   r   r   weightsTr   returned      @r6   r9   r   r   r   r   )r   r   r   r   r
   r;   reshape)rM   ottresultwtsrQ   rQ   rR   test_testAverage1   s    
zTestAverage.test_testAverage1c                 C   sB  g d}g dg dg}t dtjd}tt|ddd tt|d|dd tt dtjdd	t d g}tt|d tjt dd
 d  tt|ddt dd
 d	  tt|ddt|ddt|ddd	 g tt|d |dd tt|d|dg d tt|ddt|ddt|ddd	 g td}g d}g dg dg}t	d}g d}	ttt
||ddd ttt
||ddd ttt
||ddjdg ttt
||	ddd tttt
||ddd t
||}
tt|
d d tt|
ddg d tt|
ddddg tt|
d|dg d d S )Nr   r9   r9   r9   r9   r   r9   r   r   r   r   r9   rt   r7   r   r         @)r   r   r   r         (@r9   r   g
@)r   r   r   r   r         $@)r   r   r9   r9   r   r   )r   r9   r9   r9   r9   r9   Tr   )r   r        X@r   r   g      @      @)r   r   r   r   r   r   )r	   rL   float_r   r   r   addreducer   r   r   r;   r   )rM   w1w2xym1m2m3m4m5zrQ   rQ   rR   test_testAverage2   sF   & 
zTestAverage.test_testAverage2c                 C   sf  t d}t dd }t||g||ggddd\}}tt|t| t|j|j ttddddgdd\}}tt|t| ttddd	\}}tt|t| ttdtddd
\}}tt|t| tddgddggt}t|ddgddgg}t|dd}	t|	ddg t|dd}
t|
ddg t|d d}
t|
d t|dd}
t|
ddg d S )Nrt   rw   r9   Tr   )r6   r6   rw   r   r   r   r   )r   r   r6      Fr         ?r   r   g@      ?r   )r	   r   r   r   r   r   rK   r   )rM   r<   r=   r1r   r2r   a2da2dma2daa2dmarQ   rQ   rR   test_testAverage3   s*   
zTestAverage.test_testAverage3c                 C   s@   t tddg}t|d t tg dg dd}t|d d S )Nr9   r6   r   )r9   r6   rw   r   )FFTTr   )r   r   r   rM   r<   rQ   rQ   rR   test_onintegers_with_mask   s   
z%TestAverage.test_onintegers_with_maskc                 C   s  t jg dg dgtd}tg dg dg|d}t|}t | }t|j|j t|j|j t|dd}t|jddt|jddd	  }t|j|j t|j|j t|d
d}t|jd
dt|jd
dd	  }t|j|j t|j|j t g dg dg}	t||	d}
t j| |	|  d}t|
j|j t|
j|j t||	dd}t|j|	ddt|j|	ddd	  }t|j|j t|j|j t||	d
d}t|j|	d
dt|j|	d
dd	  }t|j|j t|j|j d S )N)r   r   r   r9   r   r   r9   r   r   r   r7   )r   y      ?       @y      @      @y      @      @y      @       @)y              "@              ?y       @      @y      @      @y      @      @r   r   r   r   r9   )r   r   r   r   r   )r   r   r   r   r   r   )r   r   )	rL   r   r~   r   r   rk   r   realimag)rM   r;   r<   avexpectedav0	expected0av1	expected1r   wavwav0wav1rQ   rQ   rR   test_complex   sR   
  
zTestAverage.test_complexc                 C   s   t jjt dddg dg dg dgd}tg ddd}tg dg dd}t|d|dd	}t g d
}t|| t|d|dd	}t g d}t|| d S )Nrv   rw   r9   r   r   r   r   r   r   )r|         Fr   r   )      @g     @g     @)r   g0N@g0N@)rL   rX   r   r	   r   r   r   r   )rM   r<   weights_unmaskedweights_maskedavg_unmaskedexpected_unmasked
avg_maskedexpected_maskedrQ   rQ   rR   test_masked_weights'  s   
zTestAverage.test_masked_weightsN)	r   r   r   r   r   r   r   r   r   rQ   rQ   rQ   rR   r      s    ",r   c                   @   $   e Zd Zdd Zdd Zdd ZdS )TestConcatenatorc                 C   sz   t td tg d td}g d}t||d}t|dd|f }tt|t t |g d t |jt|dd|f  d S )Nr9   r6   rw   r   r|   rt   r|   )r9   r   r   r   r   r   r   )r9   r9   r9   r9   r9   r   r   r9   r9   r9   r9   r9   )	r   r   r   r   r   r   
isinstancer   r;   )rM   r=   mdrT   rQ   rQ   rR   test_1d<  s   zTestConcatenator.test_1dc                 C   s0  t jdd}t jdd}t t jddd}t t jddd}t||d}t||d}td||f }t|jdk t|d d d df | t|d d dd f | t|j	t j
d||f  t||f }t|jdk t|d dd d f | t|dd d d f | t|j	t j
||f  d S )Nr|   r   r   1)r|   rb   )rb   r|   )rL   randomrandround_r   r   r   r   r   r;   r_)rM   a_1a_2m_1m_2b_1b_2r   rQ   rQ   rR   test_2dG  s    zTestConcatenator.test_2dc                 C   sp   t tjjdf }t|jddg t|jd d t ddgtjjf }t|jg d t|jd d ddg d S )Nr9   TFr6   )FFT)r   rL   rX   r
   r   r;   ri   )rM   actualrQ   rQ   rR   test_masked_constant[  s   z%TestConcatenator.test_masked_constantN)r   r   r   r   r   r   rQ   rQ   rQ   rR   r   9  s    r   c                   @      e Zd Zdd Zdd ZdS )TestNotMaskedc                 C   s  t tdddg dg dg dg dg dgd}t|d }t|d	d
g t|d	}t|d	 ddg t|d ddg t|d}t|d	 ddg t|d ddg t|jd }t|d	d
g t|jd	}t|d	 ddg t|d ddg t|jd}t|d	 ddg t|d ddg t|d< t|d	}t|d	 ddg t|d ddg t|d}t|d	 ddg t|d ddg d S )N   r|   )r   r   r9   r   r   )r   r   r   r9   r9   )r9   r9   r   r   r   )r   r   r   r   r   )r9   r9   r9   r   r   r   r      r   r9   r6   rw   r   r9   )rw   rw   rw   r   r   )r   r   r6   r   rw   )r   r6   r   r   r   )r   r   r   r   r   )r9   r9   r6   r   r   )r   r9   r6   r   )r   r   r6   rw   )r   r6   r   r   )r   rL   r	   r   r-   r   ri   r
   rM   ri   rN   rQ   rQ   rR   
test_edgesh  s>   




zTestNotMasked.test_edgesc                 C   s,  t tdddg dg dg dgd}t|d }t|tdd	d td
dd tddd g t|d}t|tddd tddd gtddd tddd gtddd tddd gtddd tddd gtddd gtddd gg tddd gg t|d}t|tdd	d gg tddd tddd gg d S )Nr   rw   ru   )r   r   r   r   r9   r9   r9   r9   )r9   r9   r9   r9   r9   r9   r9   r9   )r   r   r   r   r   r   r9   r   r   r   r            r9   r6   rt   re   )r   rL   r	   r   r,   r   rx   )rM   r<   tmprQ   rQ   rR   test_contiguous  s:   






zTestNotMasked.test_contiguousN)r   r   r   r   r   rQ   rQ   rQ   rR   r   e  s    "r   c                   @   s|   e Zd Zdd Zdd Zdd Zejdg dejd	d
ge	j
jdfe	j
jdfgdd Zdd Zdd Zdd ZdS )TestCompressFunctionsc              	   C   s(  t ttdddd}t dt}d|d< t||d}t|}t	|g d	g d
g dgg dg dg dgg t|d}t	|g dg dg dg dgg dg dg dg dgg t|d}t	|g dg dg dgg dg dg dgg dg dg dgg t|d}t|d}t|d}t	|| t	|| t	|| t|d}t	|g d	g d g d
g dgg d!g d"g d#g d$gg dg d%g dg dgg t|d&}t|d'}t|d(}t	|| t	|| t	|| t|d)}t	|g dg dg dgg dg dg dgg t|d*}t	|| t|d+}t	|g d	g d
g dgg d!g d#g d$gg dg dg dgg t|d,}t|d-}t|d.}t	|| t	|| t	|| t|d/}t	|g d	g d g d
g dgg dg d%g dg dgg t|d0}t	|| d S )1N<   rw   r   r|   rw   r   r|   Tr9   r9   r9   r   )r   r6   rw   r   )rb            )            )(   *   +   ,   )2   4   5   6   )7   9   :   ;   r   r   )r|   rt   re   ru   rv   )rb      r   r   r   )r   r   r   r   r   )r   )   r  r  r  )-   .   /   0   1   )r  3   r  r  r  )r  8   r	  r
  r  r9   )      r   r   r   )   r       !   "   )#   $   %   &   '   rU   r   )r   r6   )r|   re   ru   rv   )r  r   r   r   )r      r      )r  r  r  r  )r  r  r  r  )r  r  r  r  r5   r   )r   )r   r9   )r   r   r9   r6   )r   r6   )r9   r   )r   r   r   r6   )r   r   )
rL   r   listrf   r   r   astyper~   r&   r   )rM   r   r   r<   a2a3a4rQ   rQ   rR   test_compress_nd  s   




































	z&TestCompressFunctions.test_compress_ndc                 C   s  t tdddg dg dg dgd}tt|ddgdd	gg tt|d
g dg dg tt|dddgddgdd	gg t |jg dg dg dgd}tt|d
dgdd	gg tt|d
g dg dg tt|dd
dgddgdd	gg t |jg dg dg dgd}tt|d	gg tt|d
g dg tt|ddgdgd	gg t |jg dg dg dgd}tt|jd
 tt|d
jd
 tt|djd
 d S )Nrv   rw   r   r   r   r   r|   re   ru   r   r   )rt   re   ru   r9   r6   r   r9   r   rt   r   r9   r6   r   r   r9   )r   rL   r	   r   r   r'   _datar}   rM   r   rQ   rQ   rR   test_compress_rowcols  s$   " "  z+TestCompressFunctions.test_compress_rowcolsc                 C   sH  t tdddg dg dg dgd}tt|jg dg dg dg tt|djg dg dg dg tt|djg dg dg dg t |jg dg d	g dgd}tt|jg d	g dg d	g tt|djg dg dg dg tt|djg d	g d	g d	g t |jg dg d	g dgd}tt|jg dg dg d
g tt|djg dg dg dg tt|djg d
g d
g d
g t |jg dg d	g dgd}tt|	 t
u  tt|d	 t
u  tt|d	 t
u  tt|j	  tt|dj	  tt|dj	  d S )Nrv   rw   r   r   r   r   r   r9   r*  )r9   r9   r   r,  )r   rL   r	   r   r   r(   r;   r-  r   allr
   r.  rQ   rQ   rR   test_mask_rowcols1  sN   
 
 
 z'TestCompressFunctions.test_mask_rowcolsr   )Nr   r9   funcrowcols_axisr   r9   c                 C   sv   t tdddg dg dg dgd}tt |||d}t|t|| W d    d S 1 s4w   Y  d S )Nrv   rw   r   r   r   r   )r   rL   r	   r   r   DeprecationWarningr   r(   )rM   r   r2  r3  r   resrQ   rQ   rR   #test_mask_row_cols_axis_deprecationQ  s   
"z9TestCompressFunctions.test_mask_row_cols_axis_deprecationc                 C   s  t dd}g d}t||ddd}t||ddd}t||dd}t|jddgdd	gg t||dd}t|jg d
g dg dg t||dd}t|t |d	|d	 t||dd}t|t |d	|d	 g d}t||ddd}t||ddd}t||dd}t|jd	dgddgg t||dd}t|jg dg dg d
g t||dd}t|t |d	|d	 t|t|| t||dd}t|t |d	|d	 g d}t||ddd}t||ddd}t||}t|jt t||}t|jt t|g dddd}t|g dddd}t||dd}t|jddgd	d	gg t||dd}t|t |d	|d	 t||dd}t|jg dg dg dg t||dd}t|t |d	|d	 t|g dddd}t|g dddd}t||dd}t|jd	d	gddgg t||}t|t |d	|d	 t||dd}t|jg dg dg dg t||dd}t|t |d	|d	 t|g dddd}t|g dddd}t||dd}t|jdd	gddgg t||dd}t|t |d	|d	 t||dd}t|jg dg d
g dg t||dd}t|t |d	|d	 d S )Nr9   re   )r9   r   r   r   r   r   r   r6   rw   T)strictr   r   r   Fr   r   r   r   r   r9   r,  )r   r   r   r   r   r   )r   r   r9   r   r   r   )	rL   r	   r   r   r   r   r;   filledr   )rM   nr   r<   r=   rT   rQ   rQ   rR   test_dot]  sz   


"zTestCompressFunctions.test_dotc                 C   sn   t d}t|}ttt||tu  ttt||tu  ttt||tu  ttt||tu  d S )Nrw   )rL   eyer   r   rY   r   r   rM   r<   r=   rQ   rQ   rR   test_dot_returns_maskedarray  s   
z2TestCompressFunctions.test_dot_returns_maskedarrayc                 C   sD   t td}t td}t|||d}t||u  t|| d S )Nrw   )rw   rw   out)r   rL   r<  r   r   r   r   )rM   r<   r@  r5  rQ   rQ   rR   test_dot_out  s
   z"TestCompressFunctions.test_dot_outN)r   r   r   r)  r/  r1  pytestmarkparametrizerL   rX   	mask_rows	mask_colsr6  r;  r>  rA  rQ   rQ   rQ   rR   r     s    p 	F	r   c                   @   r   )TestApplyAlongAxisc                 C   s@   t dddd}dd }t|d|}t|ddgdd	gg d S )
Nr   r6   rw   c                 S   s   | d S Nr9   rQ   )r=   rQ   rQ   rR   myfunc  s   z*TestApplyAlongAxis.test_3d.<locals>.myfuncr9   r   re   rb   r	   r   r%   r   rM   r<   rI  xarQ   rQ   rR   test_3d  s   zTestApplyAlongAxis.test_3dc                 C   sF   t dddd}ddd}t|d|dd}t|dd	gd
dgg d S )Nr   r6   rw   r   c                 S   s   | d|  S rH  rQ   )r=   offsetrQ   rQ   rR   rI    s   z1TestApplyAlongAxis.test_3d_kwargs.<locals>.myfuncr9   )rN  r|   ru   r  )r   rJ  rK  rQ   rQ   rR   test_3d_kwargs  s   
z!TestApplyAlongAxis.test_3d_kwargsN)r   r   r   rM  rO  rQ   rQ   rQ   rR   rG        
rG  c                   @   s   e Zd Zdd ZdS )TestApplyOverAxesc                 C   s   t dddd}ttj|ddg}tdgdgdggg}t|| t||d t	< ttj|ddg}td	gd
gdggg}t|| d S )Nr   r6   rw   r   r   r   \   |   r   r  )
r	   r   r$   rL   rj   r   r   r
   r%  r~   )rM   r<   rN   ctrlrQ   rQ   rR   
test_basic  s   
zTestApplyOverAxes.test_basicN)r   r   r   rU  rQ   rQ   rQ   rR   rQ    s    rQ  c                   @   s   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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 ).
TestMedianc                 C   s4   t jjt jt jgt jt jggdd}t|t j d S )Nr   r   )rL   rX   r   infr   rM   rrQ   rQ   rR   test_pytype  s   $zTestMedian.test_pytypec                 C   s   t jjt jt jt jgt jt jggdd}t|t j t jjt jt jt jgt jt jggd d}t|t j t jjt jjt jt jgt jt jggdddd}t|jd t jjt jjt jt jgt jt jggddd d}t|jd d S )Nr   r   Tr   )rL   rX   r   r   rW  r   r;   rX  rQ   rQ   rR   test_inf  s8   



zTestMedian.test_infc                 C   sl  t d}tt j|d ttt j|tu td}tt j|d ttt j|tu d}tt j|d ttt j|tu t d	dd}tt jj|dd	t j|dd	 tt jj|d
d	t j|d
d	 tt jj|d
d	tu t d	dd}tt jj|dd	t j|dd	 tt jj|d
d	t j|d
d	 tt jj|d
d	tu d S )Nrv   r   ru   g      @r|   r   H   r   r   r9   g      R@)
rL   r	   r   rX   r   r   rY   r   rf   r   r.  rQ   rQ   rR   test_non_masked  s"   
    zTestMedian.test_non_maskedc                 C   s  t tddgd dgd  d}ttj|d ttj|jdd tttj|t	u t td	
d
ddgd dgd  d}ttj|d ttj|jdd tttj|t	u tjj|ddd}t|ddg t|jdd tt|t	u  dS )z5test the examples given in the docstring of ma.medianru   r   r   r9   r   r   rQ   shape mismatchrb   r6   r|   rt   r   r   Tr   overwrite_inputr   r   r5   N)r   rL   r	   r   rX   r   r   r   rY   r   r   )rM   r   ma_xrQ   rQ   rR   test_docstring_examples  s   "*z"TestMedian.test_docstring_examplesc                 C   s   d}t dD ]p}dD ]k}td||d}tt | |ddg}|D ] \}}ztjj|||d W q! tyA   t|||||f w t|d  |gddg}|D ]#\}}ztjj|||d W n
 tj	yj   Y qQw t|||||f q
qd S )	Nz5mask = %s, ndim = %s, axis = %s, overwrite_input = %sr|   )FTr9   )ndminr;   FTr_  )
rf   r   	itertoolsproductrL   rX   r   	ExceptionAssertionError	AxisError)rM   msgrc  r;   r   argsr   overrQ   rQ   rR   test_axis_argument_errors  s,   z$TestMedian.test_axis_argument_errorsc                 C   sD   t ddd}ttj|d t ddd}ttj|tjj d S )Nr9   Fr   T)r   r   rL   rX   r   r
   r.  rQ   rQ   rR   test_masked_0d*  s   zTestMedian.test_masked_0dc                 C   s  t tddd}ttj|tjj ttj|jdd tt	tj|tjj
ju  t tddd}ttj|d ttj|jdd tt	tj|tu t tdg dd}ttj|d	 ttj|jdd tt	tj|tu t tdg d
d}ttj|d ttj|jdd tt	tj|tu t tdg dd}ttj|d ttj|jdd tt	tj|tu t tdg dd}ttj|d ttj|jdd tt	tj|tu t tdg dd}ttj|d	 ttj|jdd tt	tj|tu t tdg dd}ttj|d	 ttj|jdd tt	tj|tu d S )Nr|   Tr   rQ   r^  Fr   r   r   )r   r9   r9   r9   r9   r   )r   r9   r9   r   r   r   r   rt   r   r   )r   rL   r	   r   rX   r   r
   r   r   rY   rZ   MaskedConstantr   r.  rQ   rQ   rR   test_masked_1d1  s@   zTestMedian.test_masked_1dc                 C   sB   t tjtg dg ddjtjtg dg ddj d S )Nr9   r6   rw   r   r   r*  )r   rL   rX   r   r   r   rM   rQ   rQ   rR   test_1d_shape_consistencyW  s   z$TestMedian.test_1d_shape_consistencyc                 C   s   d\}}t tdd|}t |d d< |dd < t tj||ftd}|d d  |d d df< tt|}td|D ]}tj	
| || |d d |f< q<tt|d d df d tt|d tt|dd	t| tt|jdd	t| d S )
N)e   r        r   rb   r7   r   r9   r   )r   rL   linspacer
   emptyrK   r	   r[   rf   r   shuffler   r   r   T)rM   r:  pr   r   idxrl   rQ   rQ   rR   r   [  s   zTestMedian.test_2dc                 C   s   t tddd}t |d d< |dd < tt|d tttj	|t
u tt|ddg d tttj	j|ddt
u  tt|d	dg d
 tttj	j|d	dt
u  tt|d	djg d d S )Nr  rb   rw         -@r   r   )      +@r}  g      /@r9   )
r   r   r   rb   r   r   r   r   r   r   )
r9   r9   r9   r   r   r   r   r9   r9   r9   )r   rL   r	   r   r
   r   r   r   rY   rX   r   r;   r.  rQ   rQ   rR   test_2d_waxisk  s   zTestMedian.test_2d_waxisc                 C   s   t jdddd}t||d dk< tt|dddgdd	gddgd
d	gg d|_tt|dddgddgddgg t jdddd}t||d dk< tt|dddgddgddgg d S )Nr   rw   r   r6   r   r   rv   rt   r   r   )r   rw   r6   c   rb   r  r   r   r|   ru   r   r   )rL   rX   r	   r   r
   r   r   r   r.  rQ   rQ   rR   rM  w  s   ("&zTestMedian.test_3dc                 C   sL   t tddd}t |d d< |dd < tt|ddt|dd d S )Nr  rb   rw   r|  r   r   r9   )r   rL   r	   r   r
   r   r   r.  rQ   rQ   rR   test_neg_axis  s   zTestMedian.test_neg_axisc                 C   s   dD ]?}t t|}t |d d< |dd < t td}t||d}|dkr.t|d nt|d t||u  tt|t	u  qd S )	N)r        >@r   g      ?@rw   r|  rQ   r?  r  r}  g      .@)
r   rL   r	   r
   r   r   r   r   rY   r   )rM   vr   r@  rY  rQ   rQ   rR   test_out_1d  s   
zTestMedian.test_out_1dc                 C   s  dD ]}}t t|dd}t |d d< |dd < t td}t|d|d}|dkrMt d	gd g d
 d	gd  dgd dgd  dgd  d}nt d	gd g d d	gd  dgd dgd  dgd  d}t|| t||u  tt	|t
u  qd S )N)r   g      D@r  r  rb   r   rw   r|  r9   )r   r@  r  r   )rb   r   r   r   TFr   r   )r~  g     1@g     5@g     9@)r   rL   r	   r   r
   r   r   r   r   rY   r   )rM   r  r   r@  rY  erQ   rQ   rR   test_out  s    
zTestMedian.test_outc                 C   sF   ddgddgddgg}t j|d}ddg}tt jj|dd| d S )Nr   r   r   r   r   )rL   rX   masked_equalr   r   )rM   ri   r_   r   rQ   rQ   rR   $test_single_non_masked_value_on_axis  s   z/TestMedian.test_single_non_masked_value_on_axisc                 C   s  dt jdtdfD ]}t jdt jdgg dg}||_t jj|d d}tt 	| t
|t j t jj| dd}tt 	| t
|t j t jj|dd}tt|t t
|dt jdg t jj|dd}tt|t t
|t jd	g t jj|d
d}tt|t t
|t jd	g q	t jdt jdgg dg}t jj|d d d	f< t
t jj|d dt j t
t jj|dddt jdg t
t jj|ddt jdg d S )NFrt   r7   r9   rw   rp  r   r   r6   r   r   )rL   r   r~   rX   r   nanr;   r   r   isscalarr   ravelr   rY   r   r
   )rM   r;   dmrY  rQ   rQ   rR   test_nan  s.    zTestMedian.test_nanc                 C   s   t jt d}t jt d}t j|d< t jj|d< tt jj|d|d| t jt d}tt jj|d|d| t jt d	}tt jj||d| d S )
N)r   )rw   r   )r6   r9   rD   r   r?  rw   r9   rQ   )	rL   rX   r   r   r   r  r
   r   r   )rM   or   rQ   rQ   rR   test_out_nan  s   
zTestMedian.test_out_nanc                 C   s  t jt jdtd}t jj|d d d< t j|d< tt j|t j tt jj|ddt j t jt jdtd	ddd}t |j
d dk|_| }t j|d	< t j|d
< tt j|t j tt t j| t jj|dd}t j|d< t j|d< tt j|d| t jj|dd}t j|d< t j|d< tt j|d| t jj|dd}t j|d< t j|d< tt j|d| d S )Nr   r7   rw   r6   r   r   r   r9   rp  )r9   r9   r6   r6   rw   r"  r9   rw   r#  )rL   rX   r   r	   rK   r
   r  r   r   r   r}   r;   copyr   r  r   )rM   r<   aorigr=   rQ   rQ   rR   test_nan_behavior  s0   
 







zTestMedian.test_nan_behaviorc                 C   s   t jg dg dgt jd}t jj||dkd}tt jj|ddd tt jj|ddjd tt jj|d	d|d	  tt j|d d S )
N)rw   rw      r7   rw   r   r9   r   r  Fr   )rL   r   uint8rX   r   r   r   r;   r   rQ   rQ   rR   test_ambigous_fill  s   zTestMedian.test_ambigous_fillc                 C   s  t jt j fD ]v}t |t jgt jt jgg}t jj|t |d}tt jj|dd|t jg tt jj|dd|t jg tt j|| t t jt j|gt jt j|gg}t jj|t |d}t	t jj|dd| t	t jj|ddj
d t	t jj|dd|d  t	t j|| t ||g||gg}tt j|| tt jj|dd| tt jj|dd| t j|d| dgdt jt jd	gd
t jt j|ggt jd}t jj|t |d}|dkrtt jj|dddd| dg tt j|d ntt jj|dddd| dg tt j|d tt jj|dddd|g tddD ]P}tddD ]G}t t jg| |g|  gd }t jj|t |d}tt j|| tt jj|dd| tt jj|ddt jg| |g|   q4q-qd S )Nr   r   r   r9   Fre   iru  r|   r   r7   r         @r   g      @g      $g      "g      rt  rb   r6   )rL   rW  r   r  rX   r   isnanr   r   r   r;   float32rf   )rM   rW  r<   rl   rm   rQ   rQ   rR   test_special  sR   "
  "zTestMedian.test_specialc                 C   s  t jt jg td}t !}|t tt j	|t j
 t|jd jtu  W d    n1 s2w   Y  t jt jg tdd}t (}|t tddt tt j	|t j
 t|jd jtu  W d    n1 sqw   Y  t jt jg tdd}tt jj	|dd| tt jj	|d	d| t jt jt j
tdd}tjd
d$}tddt tt jj	|dd| t|d jtu  W d    d S 1 sw   Y  d S )Nr7   r   rw   )r8   rc  always r6   r   r9   T)record)rL   rX   r   r   rK   r   r  RuntimeWarningr   r   r  r   logcategorywarningsfilterwarningsr   catch_warnings)rM   r<   wr=   rQ   rQ   rR   
test_empty.  s,   

"zTestMedian.test_emptyc                 C   sX   t jt d}ttt j|tt	 t j
|d< ttt j|tt	 d S )Nr  r6   )rL   rX   r   r	   r   rY   r   r%  rW   rK   r  )rM   r  rQ   rQ   rR   test_objectK  s   
 zTestMedian.test_objectN)r   r   r   rZ  r[  r]  rb  rl  rm  ro  rr  r   r  rM  r  r  r  r  r  r  r  r  r  r  r  rQ   rQ   rQ   rR   rV    s.    &	#	+rV  c                   @   s4   e Zd Zdd Zdd Zdd Zdd Zd	d
 ZdS )TestCovc                 C   s   t tjd| _d S Nr   )r   rL   r   r   ri   rq  rQ   rQ   rR   setupT  s   zTestCov.setupc                 C   sZ   | j }tt|t| ttj|ddt|dd ttj|dddt|ddd d S )NFrowvarTr  bias)ri   r   rL   r   r.  rQ   rQ   rR   test_1d_without_missingW  s   zTestCov.test_1d_without_missingc                 C   sb   | j dd}tt|t| ttj|ddt|dd ttj|dddt|ddd d S )Nrw   r   Fr  Tr  )ri   r   r   rL   r   r.  rQ   rQ   rR   test_2d_without_missing_  s   zTestCov.test_2d_without_missingc                 C   s>  | j }t|d< || 8 }| }tt|t| ttj|ddt|dd ttj|dddt|ddd zt|dd W n	 tyJ   Y nw |dd }tt||d d d t||d d d  ttj||d d d ddt||d d d dd ttj||d d d dddt||d d d ddd d S )Nr   Fr  Tr  allow_maskedr9   )ri   r
   meanrk   r   rL   r   
ValueError)rM   r   nxrQ   rQ   rR   test_1d_with_missingg  s,   ,zTestCov.test_1d_with_missingc                 C   s6  | j }t|d< |dd}tt|t}t||j	}||
dd d d f  d}tt|t||jd d  |d   tt|ddtj|dd|jd  |  t|j	|}||
d d}tt|d	d
tj|d	d
|jd d  |d   tt|d	ddtj|d	dd|jd  |  d S )Nr   rw   r   r9   r   r   Tr  Fr  r  )ri   r
   r   rL   logical_notr   r%  rg   r   ry  r  r9  r   r   r   )rM   r   validfracxfrQ   rQ   rR   test_2d_with_missing  s8    zTestCov.test_2d_with_missingN)r   r   r   r  r  r  r  r  rQ   rQ   rQ   rR   r  R  s    r  c                   @   sD   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S )TestCorrcoefc                 C   s(   t tjd| _t tjd| _d S r  )r   rL   r   r   ri   data2rq  rQ   rQ   rR   r    s   zTestCorrcoef.setupc                 C   s   | j | j}}t|}t||}t P}td ttt|dd |	td t
tj|ddt|dd t
t|dd| t
t||dd| t
t|dd| t
t||dd| W d    d S 1 shw   Y  d S )Nr  r   ddofbias and ddof have no effectr   rw   ri   r  rL   r   r   r  simplefilterr   r4  filterr   )rM   r   r   r   	expected2suprQ   rQ   rR   	test_ddof  s   

"zTestCorrcoef.test_ddofc                 C   s   | j | j}}t|}t 7}td ttt||dd ttt||dd ttt|dd |	td t
t|dd| W d    d S 1 sIw   Y  d S )Nr  TFr  r  r9   r  )rM   r   r   r   r  rQ   rQ   rR   	test_bias  s   

"zTestCorrcoef.test_biasc                 C   s   | j }tt|t| ttj|ddt|dd t }|td ttj|dddt|ddd W d    d S 1 s@w   Y  d S )NFr  r  Tr  )ri   r   rL   r   r   r  r4  rM   r   r  rQ   rQ   rR   r    s   
"z$TestCorrcoef.test_1d_without_missingc                 C   s   | j dd}tt|t| ttj|ddt|dd t }|td ttj|dddt|ddd W d    d S 1 sDw   Y  d S )Nrw   r   Fr  r  Tr  )ri   r   r   rL   r   r   r  r4  r  rQ   rQ   rR   r    s   
"z$TestCorrcoef.test_2d_without_missingc              	   C   s  | j }t|d< || 8 }| }tt|t| ttj|ddt|dd t }|t	d ttj|dddt|ddd W d    n1 sMw   Y  zt|dd W n	 t
yc   Y nw |dd }tt||d d d t||d d d  ttj||d d d ddt||d d d dd t ?}|t	d tt||d d d t||d d d dd	 tt||d d d t||d d d d
d W d    d S 1 sw   Y  d S )Nr   Fr  r  Tr  r  r9   r  r6   r  )ri   r
   r  rk   r   rL   r   r   r  r4  r  )rM   r   r  r  rQ   rQ   rR   r    sB   
,"z!TestCorrcoef.test_1d_with_missingc                 C   s$  | j }t|d< |dd}t|}t|}t|d dd df |d dd df  t Z}|td tt|ddd dd df |d dd df  tt|ddd dd df |d dd df  tt|ddd dd df |d dd df  W d    d S 1 sw   Y  d S )	Nr   rw   r   r  r   r  r9   r  )	ri   r
   r   r   rL   r   r   r  r4  )rM   r   rN   rO   r  rQ   rQ   rR   r    s$   
*"z!TestCorrcoef.test_2d_with_missingN)
r   r   r   r  r  r  r  r  r  r  rQ   rQ   rQ   rR   r    s    r  c                   @   r   )TestPolynomialc                 C   s  t jd}t jddd}tt||dt ||d |t}t|d< |t}t |d< |d< t||d d df dd	d
\}}}}}t j|dd  |dd df 	 dd	d
\}}	}
}}t
|||||f||	|
||fD ]	\}}t|| qpt||d d df dd	d
\}}}}}t j|dd |dddf dd	d
\}}	}
}}t
|||||f||	|
||fD ]	\}}t|| qt||dd	d
\}}}}}t j|dd |ddd d f dd	d
\}}	}
}}t
|||||f||	|
||fD ]	\}}t|| qt jdd }| }|dd }|dd }|dd }t||dd	|d\}}}}}t j||dd	|d\}}	}
}}t|| t
|||||f||	|
||fD ]
\}}t|| qMd S )Nrb   r  r   r6   rw   r   rB   r   r   Trh   r9   )rh   r  )rL   r   r   r   r   r   viewr   r
   rk   zipr  r   )rM   r   r   CRKSDrT   rY  rn   rp   r   r<   a_r  woxsyswsrQ   rQ   rR   test_polyfit  s@   

&"&&0&4&
&zTestPolynomial.test_polyfitc                 C   s   t jd}t jddd}t j|d< t j|d< |t}|t}t|d< t|d< t||ddd	\}}}}}t j|d
d |d
dd d f ddd	\}}	}
}}t	|||||f||	|
||fD ]	\}}t
|| qbd S )Nrb   r  r   r6   r   r  rw   Tr  r9   )rL   r   r   r   r  r  r   r
   r   r  r   )rM   r   r   r  r  r  r  r  rT   rY  rn   rp   r   r<   r  rQ   rQ   rR   test_polyfit_with_masked_NaNs)  s   



4&z,TestPolynomial.test_polyfit_with_masked_NaNsN)r   r   r   r  r  rQ   rQ   rQ   rR   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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 )"TestArraySetOpsc                 C   sn   g d}t |ddd}tt|d t t|d tg dg dd t|d g d	 t|d
 g d d S )Nr9   r9   r9   r6   r6   rw   Treturn_indexreturn_inverser   rp  r   r   r9   )r   rw   r|   r6   )r   r   r   r9   r9   r6   )r   r   r   r   r   r   r   rQ   rQ   rR   test_unique_onlist<  s   z"TestArraySetOps.test_unique_onlistc                 C   s   t g dg dd}t|ddd}t|d t g dg dd t|d	 g d
 t|d g d d|_t g dg ddd}t|ddd}t|d t g dg dd t|d	 g d
 t|d g d d S )Nr  )r   r   r9   r   r9   r   r   Tr  r   )r9   r6   rw   r   r   r   r   r9   r9   )r   rw   r|   r6   r6   )r   r   rw   r9   rw   r6   rw   )ri   r;   
fill_value)r   r   r   r  r   rQ   rQ   rR   test_unique_onmaskedarrayE  s   z)TestArraySetOps.test_unique_onmaskedarrayc                 C   s   t g ddd}t|ddd}t|d t dgdgd t|d dg t|d g d t}t|ddd}t|d t t t|d dg t|d dg d S )	Nr   Tr   r  r   r9   r6   r   )r   r   r   r
   r   rQ   rQ   rR   test_unique_allmaskedU  s   z%TestArraySetOps.test_unique_allmaskedc                 C   sd   t tdg dd}tg dg dd}t|}t|| t|d|d t|j|j d S )Nr|   r9   r   r   r   r9   r   )r9   r9   r9   r   )r9   r   r   r9   r   )r   rL   r	   r   r#   r   r9  r;   )rM   r   rO   rN   rQ   rQ   rR   test_ediff1dd  s   
zTestArraySetOps.test_ediff1dc                 C      t tdg dd}t|td}tg dg dd}t|| t|d|d t|j|j t|g dd}tg d	g d
d}t|| t|d|d t|j|j d S )Nr|   r  r   )to_begin)r   r9   r9   r9   r   )r9   r9   r   r   r9   r   rp  )r9   r6   rw   r9   r9   r9   r   )r   r   r   r9   r   r   r9   	r   rL   r	   r#   r
   r   r   r9  r;   rM   r   rN   rO   rQ   rQ   rR   test_ediff1d_tobeginm     

z$TestArraySetOps.test_ediff1d_tobeginc                 C   r  )Nr|   r  r   )to_end)r9   r9   r9   r   r   )r9   r   r   r9   r9   r   rp  )r9   r9   r9   r   r9   r6   rw   )r9   r   r   r9   r   r   r   r  r  rQ   rQ   rR   test_ediff1d_toend|  r  z"TestArraySetOps.test_ediff1d_toendc                 C   s   t tdg dd}t|ttd}tg dg dd}t|| t|d|d t|j|j t|g dtd}tg d	g d
d}t|| t|d|d t|j|j d S )Nr|   r  r   r  r  )r   r9   r9   r9   r   r   )r9   r9   r   r   r9   r9   r   rp  )r   r9   r9   r9   r   r9   r6   rw   )r9   r9   r   r   r9   r   r   r   r  r  rQ   rQ   rR   test_ediff1d_tobegin_toend  s   

z*TestArraySetOps.test_ediff1d_tobegin_toendc                 C   s   t d}t|}tg dg dd}t|| tt|t t|d|d t|j	|j	 t|t
t
d}tg dg dd}tt|t t|d|d t|j	|j	 d S )	Nr|   )r9   r9   r9   r9   )r   r   r   r   r   r   r  r   r   )rL   r	   r#   r   r   r   r   r   r9  r;   r
   r  rQ   rQ   rR   test_ediff1d_ndarray  s   

z$TestArraySetOps.test_ediff1d_ndarrayc                 C   sT   t g dg dd}t g dg dd}t||}t g dg dd}t|| d S )N)r9   rw   rw   rw   r  r   )rw   r9   r9   r9   )r9   rw   r   r,  )r   r!   r   )rM   r   r   rN   rO   rQ   rQ   rR   test_intersect1d  s
   
z TestArraySetOps.test_intersect1dc                 C   s$  t g dg dd}t g dg dd}t||}t|t g d t g dg dd}g d}t||}t|t g dg d	d t g d
}t g d}t||}tt|t t|g d t g dg dd}t g dg d	d}t||}tt|t t|g d tg tg g  d S )Nr9   r6   r|   re   r   r   r   r   r   r9   r   r9   r6   rw   r   r|   r   r8  )rw   r   re   )r9   r6   rw   r   r|   )rw   r   re   r   r  rp  )rt   r|   r   r   )r9   ru   r6   rw   )r   r9   r   r   )rt   r|   r   ru   )r   r   r   r   r   r   r   rM   r<   r=   rN   rQ   rQ   rR   test_setxor1d  s&   



zTestArraySetOps.test_setxor1dc                 C   s   t dg d}t g d}d|d< t||d}tg dg dd}tdtd	}d
|d< d
|d< d
|d< t||}tt|t	 t
|| t |||j  |j @ }t
|| d S )Nr   )r6   rw   r   r9   )r9   r6   r   r   )	r   rb   r  r  r9   rw   r  r   r  )	r   r9   r   r9   r   r9   r   r9   r   )ri   r;   r7   Tr   r,  )r   r6   rw   )rL   r	   r   r   r   r~   r0   r   r   r   r   r;   )rM   r<   r;   r=   ecrT   r   rQ   rQ   rR   	test_isin  s    

zTestArraySetOps.test_isinc                 C   s   t g dg dd}t g dg dd}t||}t|g d t g dg dd}t g dg d	d}t||}t|g d
 tg tg g  d S )Nr  r  r   r  r8  )TTTFTr|   r|   r6   r9   r   r9   r|   r   r,  )TTFTT)r   r"   r   r   r  rQ   rQ   rR   	test_in1d  s   

zTestArraySetOps.test_in1dc                 C   s   t g dg dd}t g dg dd}ttt||t||dd t g dg dd}t g d	g d
d}ttt||t||dd tg tg g dd d S )Nr  r  r   r  r8  T)invertr  r  r,  )r   r   rL   r  r"   r   r=  rQ   rQ   rR   test_in1d_invert  s     z TestArraySetOps.test_in1d_invertc           	      C   s   t g dg dd}t g dg dd}t||}t g dg dd}t|| t g dg dgg d	g d
gd}t g dg dd}t g dg dd}t||}t|| tg tg g  d S )N)r9   r6   r|   re   r|   r   r8  r   r  )r9   r6   rw   r   r|   re   r   )r   r   r   r   r   r   r9   r+  r   r   r,  r   r  )r   r9   r6   rw   r   r|   )r   r    r   r   )	rM   r<   r=   rN   rO   r   r   ezr   rQ   rQ   rR   test_union1d  s   

$

zTestArraySetOps.test_union1dc                 C   s   t g dg dd}t g d}t||}t|t g dg dd td}td}tt||t dd	g t g tjg d}tt|g jtj d S )
N)rt   r|   r   re   re   r9   r6   r9   )r   r   r   r   r   r   r   r9   r   )r6   r   rw   rw   r6   r9   r|   )rt   re   r   r,  rb   ru   rv   )r   r   r   r	   rL   uint32r8   r  rQ   rQ   rR   test_setdiff1d  s   
zTestArraySetOps.test_setdiff1dc                 C   s8   t g d}t g d}tt||t dg d S )N)r<   r=   rT   )r<   r=   rp   rT   )rL   r   r   r   r=  rQ   rQ   rR   test_setdiff1d_char_array  s   z)TestArraySetOps.test_setdiff1d_char_arrayN)r   r   r   r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  rQ   rQ   rQ   rR   r  :  s"    		r  c                   @   r   )TestShapeBasec                 C   sn   t g dg dd}t|}t|jd t|jj|jj t|jd t|jj|jj t|jj|jj d S )Nr+  r*  r   r  r  )r   r   r   r   r;   ri   r=  rQ   rQ   rR   test_atleast_2d$  s   zTestShapeBase.test_atleast_2dc                 C   s  t d}t|jd t|jj|j t|jj|j t dd}|D ]}t|jd t|jj|j t|jj|j q!td}t|jd t|jj|j t|jj|j tdd}|D ]}t|jd t|jj|j t|jj|j q[td}t|jd t|jj|j t|jj|j tdd}|D ]}t|jd t|jj|j t|jj|j qtd}t|jd t|jj|jj d S )Nr   rU   r   rC   r   )r   r   r   r;   ri   r   r   r1   )rM   r=   r<   rQ   rQ   rR   test_shape_scalar.  s<   


zTestShapeBase.test_shape_scalarN)r   r   r   r  r  rQ   rQ   rQ   rR   r  "  rP  r  c                   @   r   )	TestStackc                 C   s   t g dg dd}t g dg dd}t||gdd}t|jd t|j|d j t|j|d	 j t||g}t|j|j t|j|j t||gd	d}t|jd
 t|j|d d df j t|j|d d d	f j d S )Nr+  r*  r   rv   ru   re   r   r   r   r  r9   rw   r6   r   r2   r   r   r   r;   r3   ri   rM   r<   r=   rT   r   rQ   rQ   rR   test_stack_1d\  s   zTestStack.test_stack_1dc                 C   s   t g ddd}t g ddd}t||gdd}t|jd t|j|d j t|j|d	 j t||g}t|j|j t|j|j t||gd	d}t|jd
 t|j|d d df j t|j|d d d	f j d S )Nr+  Tr   r	  Fr   r   r  r9   r
  r  r  rQ   rQ   rR   test_stack_masksn  s   zTestStack.test_stack_masksc           
      C   s  d}t jdd|}t jdd|}t jdd|t}t jdd|t}t||d}t||d}t||gdd}d| }	t|j|	 t	|j
|d j
 t	|j
|d j
 t||gd	d}|d }	t|j|	 t	|j
|d
 j
 t	|j
|d j
 d}t jdd|}t jdd|}t jdd|t}t jdd|t}t||d}t||d}t||gdd}d| }	t|j|	 t	|j
|d j
 t	|j
|d j
 t||gd	d}|d }	t|j|	 t	|j
|d
 j
 t	|j
|d j
 d S )Nr
  r   rb   r6   r   r   r5   r9   r   ).r   ).r9   )rw   r6   r   r|   )rL   r   randintr%  r~   r   r2   r   r   r   r;   )
rM   shpd1d2r   r   a1r&  rT   c_shprQ   rQ   rR   test_stack_nd  sD   zTestStack.test_stack_ndN)r   r   r   r  r  r  rQ   rQ   rQ   rR   r  Z  s    r  )J__doc__r  rd  rB  numpyrL   numpy.testingr   r   numpy.ma.testutilsr   r   r   r   numpy.ma.corer   r	   r
   r   r   r   r   r   r   r   r   numpy.ma.extrasr   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/   r0   r1   r2   r3   r4   r   r   r   r   rG  rQ  rV  r  r  r  r  r  r  rQ   rQ   rQ   rR   <module>   s<   4
  ,F    }Ck: i8