a
    /'a+                     @   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   =lib/python3.9/site-packages/scipy/stats/tests/common_tests.pycheck_named_results   s    r   c           
      C   s   | j dg|R  }t|d |dkr0d\}}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|sFtj||d|d d ntt||dt|   t|stj|||  |d|d d n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   s4   t |r0| dd |}tj||d|d d 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 |r8| dd |}tj||||  d|d d d S )Nc                 S   s   | |  S r   r   r   r   r   r   r   E   r   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                    sV   t |rB|  fdd|}tj||t |d d|d d ntt | d S )Nc                    s   t |   dS )N   r*   powerr   r/   r   r   r   L   r   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    
r@   c                    sf   t |rH|  fdd|}tj||d t |d dd|d d nt |sbtt | d S )Nc                    s   t |   dS )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    

rC   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    
rE   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    rF   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:   rA   scalec                    s   g | ]}j  d |iqS rG   rD   .0r?   r.   r$   r   r   
<listcomp>k   r   z,check_entropy_vect_scale.<locals>.<listcomp>g+=)r   )r   r'   c                    s   g | ]}j  d |iqS rH   rI   rJ   rL   r   r   rM   r   r   )r*   asarrayrD   ZravelZreshapeshaper   )r$   r.   ZscZv_entZs_entr   rL   r   check_entropy_vect_scaleg   s    rQ   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rt| 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   }	 |	rx|	   i  fdd|D }
t||
 d vrt| jdg R i | jdgR   qd	d
i tt| jfi  d S )N,  c                    s   g | ]}|g R  qS r   r   rK   meth)
shape_argsr   r   r   rM      r   z$check_named_args.<locals>.<listcomp>c                    s"   g | ]}|g R i qS r   r   r`   )arB   r   r   r   rM      r   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   rb   ri   methsZ	signatureZshape_argnamesZshapes_Zvalsnamesr0   r   )rc   rB   rb   r   r   check_named_args   s2    
 rt   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rtj	d}| j|d|d d| _ | j 
 }| j|dtjdd}t|| t| j 
 | || _ d S )N  size   default_rngr   )rv   random_stater'   )ry   r*   randomseedrvsr   r   ZRandomStatehasattrrx   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 ]R}| j|  || j|k || jk @  }|D ]&}||g|R  }t	
|jtjk qjq@d S )N      ?      ?      ?c                    s   g | ]}  |qS r   ZastyperK   tpZx0r   r   rM      r   z$check_meth_dtype.<locals>.<listcomp>)rY   r*   int_float16float32float64	_argcheckrc   br   r,   dtypefloat_)r$   r.   rr   q0x_castr   ra   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 ]8}| j| jfD ]&}||g|R  }t|j	t j
k q@q0d S )Nr   c                    s   g | ]}  |qS r   r   r   r   r   r   rM      r   z#check_ppf_dtype.<locals>.<listcomp>)r*   rO   r   r   r   rY   rZ   r   r,   r   r   )r$   r.   Zq_castqra   r   r   r   r   check_ppf_dtype   s    r   c                    sp  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 qDd S )Nc                 W   s,   t |}d}| ||d  g|R  | jS )Ng|=y              ?)r*   rO   imag)fr   r.   hr   r   r   deriv   s    
z check_cmplx_deriv.<locals>.deriv)r   gRQ?r   c                    s   g | ]}  |qS r   r   r   r   r   r   rM      r   z%check_cmplx_deriv.<locals>.<listcomp>r   r   )rY   r*   r   r   r   r   r   rc   r   pdfr#   rT   r   rV   rX   Zlogpdf)r$   r.   r   r   r   r   r#   rT   r   r   r   check_cmplx_deriv   s"    


:" $$r   c           
      C   s  | 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|| || _ d S )Nru   rv   rw   r   r   r   )rv   )	ry   r|   pickledumpsloadsr   r   rY   r#   )
r$   r   r~   r?   r   Z	unpickledr   ZmediansZfrozen_distZpklr   r   r   check_pickling  s(    

$

r   c                 C   sZ   t | tjrddd}nddi}| |i |}|j| | jksBJ |j| | jksVJ d S )Nr   r'   )locrG   r   )rS   r   Zrv_continuousrc   r   )r$   r   Zlocscalervr   r   r   check_freezing5  s    r   c           	         sh   t jd  j| }t|j|d|  |sdt j fdd|d}t jd || }t||dd d S )N{   z%s: rvs failed to broadcastc                     s
    j |  S r   )r|   )allargsdistfuncr   r   r   G  r   z%check_rvs_broadcast.<locals>.<lambda>)ZotypesgvIh%<=r   )r*   rz   r{   r|   r   rP   Z	vectorizer   )	r   r%   r   rP   Z
shape_onlyZotypeZsampler|   Zexpectedr   r   r   check_rvs_broadcastB  s    
r   )F)'r   Znumpyr*   Znumpy.testingZtestingr   r   r   Zpytestr   rp   Znumpy.ma.testutilsr   Z	testutilsr	   Zscipy._lib._utilr   rf   Zscipyr   r   r&   r4   r8   r9   r@   rC   rE   rF   rQ   r\   rt   r   r   r   r   r   r   r   r   r   r   r   <module>   s4   
		')	$