o
    Eb7-                     @   s   d dl Z d dlZd dlmZ d dlmZmZ d dlm	Z
 d dlm  mZ d dlmZ d dlmZ d-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 ).    N)assert_allcloseassert_equal)raises)getfullargspec_no_self)statsFc                 C   sH   t |D ]\}}|rt| | t| | qt| | t| | qd S N)	enumeratema_nptr   getattrnpt)resZ
attributesmaiattr r   @/usr/lib/python3/dist-packages/scipy/stats/tests/common_tests.pycheck_named_results   s
   r   c           
      C   s   | j dg|R  }t|d |dkrd\}}nd\}}| jdd |d}tj|d|||d	d
 | j| \}}| j|g|R  }	t|	d d S )Nr         ?Zncf)h㈵>r   )Hz>r   c                 S   s   dS )N   r   xr   r   r   <lambda>!       z%check_normalization.<locals>.<lambda>)argsT)atolrtolerr_msgverbose)momentr   r   expectsupportcdf)
distfnr   distnameZnorm_momentr   r   Znormalization_expectZ_aZ_bZnormalization_cdfr   r   r   check_normalization   s   
r&   c                 C   s   | j dg|R  }| j dg|R  }t|s#tj||d|d d ntt||dt|   t|sHtj|||  |d|d d d S tt||dt|   d S )	Nr      
   z - 1st momentZdecimalr   z - 1st moment -infinite, m1=%sz - 2ndt momentz - 2nd moment -infinite, m2=%s)r    npZisinfr   assert_almost_equalassert_str)r$   argmvmsgm1m2r   r   r   check_moment*   s    


r4   c                 C   s8   t |r| dd |}tj||d|d d d S d S )Nc                 S   s   | S r   r   r   r   r   r   r   >   r   z#check_mean_expect.<locals>.<lambda>   z - 1st moment (expect)r)   r*   isfiniter!   r   r+   )r$   r.   r/   r1   r2   r   r   r   check_mean_expect<   s   
r8   c                 C   s@   t |r| dd |}tj||||  d|d d d S d S )Nc                 S   s   | |  S r   r   r   r   r   r   r   E   s    z"check_var_expect.<locals>.<lambda>r5   z - 2st moment (expect)r)   r6   )r$   r.   r/   r0   r1   r3   r   r   r   check_var_expectC   s   
r9   c                    sX   t |r"|  fdd|}tj||t |d d|d d d S tt | d S )Nc                       t |   dS )N   r*   powerr   r/   r   r   r   L       z#check_skew_expect.<locals>.<lambda>g      ?r5   z - skewr)   )r*   r7   r!   r   r+   r=   r,   isnan)r$   r.   r/   r0   sr1   Zm3er   r>   r   check_skew_expectJ   s   

rB   c                    sl   t |r%|  fdd|}tj||d t |d dd|d d d S t |s4tt | d S d S )Nc                    r:   )N   r<   r   r>   r   r   r   U   r?   z#check_kurt_expect.<locals>.<lambda>g      @r'   r   z - kurtosis)r   r   r   )	r*   r7   r!   r   r   r=   Zisposinfr,   r@   )r$   r.   r/   r0   kr1   Zm4er   r>   r   check_kurt_expectS   s   


rE   c                 C   s&   | j | }tt| |d  d S )Nztest Entropy is nan)entropyr   r,   r*   r@   )r$   r.   r1   Zentr   r   r   check_entropy\   s   
rG   c                 C   s$   t | j| |j| g|R   d S r   )r   r   Z_entropy)r$   r   Z
superclassr   r   r   check_private_entropya   s   rH   c                    s   t ddgddgg}j d|i} fdd| D }t ||j}t||dd	 g d
}j d|i} fdd|D }t ||j}t||dd	 d S )Nr   r'   r;   rC   scalec                       g | ]
}j  d |iqS rI   rF   .0rA   r.   r$   r   r   
<listcomp>k       z,check_entropy_vect_scale.<locals>.<listcomp>g+=)r   )r   r'   c                    rJ   rK   rL   rM   rO   r   r   rP   r   rQ   )r*   asarrayrF   ZravelZreshapeshaper   )r$   r.   ZscZv_entZs_entr   rO   r   check_entropy_vect_scaleg   s   rU   c                 C   sJ  | j | }t| tjr|d d |d f}t| j|g|R  ddg t| j|g|R  ddg | jdvrZt| j	|g|R  t
j dg t| j|g|R  dt
j g t| jddgg|R  | t| jddgg|R  |d d d  tt
| jddgg|R    tt
| jddgg|R    d S )Nr   r   g        r   )ZskellamZdlaplacer'   )r"   
isinstancer   Zrv_discreter   r   r#   sfnamelogcdfr*   inflogsfppfisfr,   r@   all)r$   r   r   r   r   r   check_edge_supportw   s   

""(&*r`   c                    s  t | j}t|jd u  t|jd u  t|j  tt|jt|k |j	d t
|  }| jr?| jdd }nd}tt
|| jk tt
|t
|k tfdd|D }ttt| |d d  d d  i }	 |	r|	   i  fdd|D }
t||
 d vrt| jdg R i | jdgR   |	sd	d
i tt| jfi  d S )N,  c                    s   g | ]
}|g R  qS r   r   rN   meth)
shape_argsr   r   r   rP      rQ   z$check_named_args.<locals>.<listcomp>c                    s"   g | ]}|g R i qS r   r   rd   )arD   r   r   r   rP      s   " nr   Zkaboom*   )_getfullargspecZ_parse_argsr   r,   ZvarargsZvarkw
kwonlyargslistdefaultsr   lenZshapesreplacesplitZnumargsr*   r_   r7   updatepopZassert_array_equalkeysr   r    assert_raises	TypeErrorr#   )r$   r   rf   rm   methsZ	signatureZshape_argnamesZshapes_Zvalsnamesr0   r   )rg   rD   rf   r   r   check_named_args   s4   
 
rx   c           	      C   s   | j }tjd d | _ | j|ddi}d| _ | j|ddi}t|| tjd| _ | j|ddi}t|| ttjdrOtj	d}| j|d|d d| _ | j 
 }| j|dtjdd}t|| t| j 
 | || _ d S )N  size   default_rngr   )rz   random_stater'   )r}   r*   randomseedrvsr   r   ZRandomStatehasattrr|   Z	get_state)	r$   r   rndmr0r1Zr2rngZ
orig_stateZr3r   r   r   check_random_state_property   s&   

r   c                    s   g d}| j |g|R    fddtjtjtjtjfD }|D ])}| j|  || j|k || jk @  }|D ]}||g|R  }t	
|jtjk q5q d S )N      ?      ?      ?c                       g | ]}  |qS r   ZastyperN   tpZx0r   r   rP          z$check_meth_dtype.<locals>.<listcomp>)r]   r*   int_float16float32float64	_argcheckrg   br   r,   dtypefloat_)r$   r.   rv   q0x_castr   re   valr   r   r   check_meth_dtype   s   

r   c                    sn   t g d  fddt jt jt jfD }|D ]}| j| jfD ]}||g|R  }t|j	t j
k q qd S )Nr   c                    r   r   r   r   r   r   r   rP      r   z#check_ppf_dtype.<locals>.<listcomp>)r*   rS   r   r   r   r]   r^   r   r,   r   r   )r$   r.   Zq_castqre   r   r   r   r   check_ppf_dtype   s   r   c                    sn  dd }| j g dg|R    fddtjtjtjtjfD }|D ]}| j|  || j|k || jk @  }| j	|g|R  | j
|g|R  | j|g|R  }}}t|| j
|g|R  |dd t|| j|g|R  || dd t|| j|g|R  | dd t|| j|g|R  | | dd t|| j|g|R  || j	|g|R  | j	|g|R   dd q"d S )Nc                 W   s,   t |}d}| ||d  g|R  | jS )Ng|=y              ?)r*   rS   imag)fr   r.   hr   r   r   deriv   s   
z check_cmplx_deriv.<locals>.deriv)r   gRQ?r   c                    r   r   r   r   r   r   r   rP      r   z%check_cmplx_deriv.<locals>.<listcomp>r   r   )r]   r*   r   r   r   r   r   rg   r   pdfr#   rX   r   rZ   r\   Zlogpdf)r$   r.   r   r   r   r   r#   rX   r   r   r   check_cmplx_deriv   s$   

:" $$r   c                 C   sR  | j }d| _ | j|ddi t| }| j|ddi}t|}|j|ddi}t|| | jdg|R  |jdg|R  g}t|d |d  t| j|d g|R  |j|d g|R   | | }t|}	t|	}|jdd}|jdd}t|| t	| dr| j
}
t|
}t|}|
j|j  krdksJ  J || _ d S )	Nry   rz   r{   r   r   r   )rz   fit)r}   r   pickledumpsloadsr   r   r]   r#   r   r   __name__)r$   r   r   rA   r   Z	unpickledr   ZmediansZfrozen_distZpklZfit_functionZpickled_fit_functionZunpickled_fit_functionr   r   r   check_pickling  s2   

$




 
r   c                 C   sZ   t | tjrddd}nddi}| |i |}|j| | jks!J |j| | jks+J d S )Nr   r'   )locrI   r   )rW   r   Zrv_continuousrg   r   )r$   r   ZlocscaleZrvr   r   r   check_freezing<  s   r   c           	         sl   t jd  j| }t|j|d|  |s4t j fdd|d}t jd || }t||dd d S d S )N{   z%s: rvs failed to broadcastc                     s
    j |  S r   )r   )allargsdistfuncr   r   r   N  s   
 z%check_rvs_broadcast.<locals>.<lambda>)ZotypesgvIh%<=r   )r*   r~   r   r   r   rT   Z	vectorizer   )	r   r%   r   rT   Z
shape_onlyZotypeZsampler   Zexpectedr   r   r   check_rvs_broadcastI  s   
r   )F)'r   Znumpyr*   Znumpy.testingZtestingr   r   r   Zpytestr   rt   Znumpy.ma.testutilsr   Z	testutilsr	   Zscipy._lib._utilr   rj   Zscipyr   r   r&   r4   r8   r9   rB   rE   rG   rH   rU   r`   rx   r   r   r   r   r   r   r   r   r   r   r   <module>   s6    
		')	+