a
    hA                     @   sb  d dl Z d dlZd dlmZ d dlZd dlm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mZmZmZmZmZmZmZm Z m!Z! d dl"m#Z# d dl$m%Z% d	Z&e j'd
k Z(g dZ)g dZ*g dZ+g dZ,g dZ-g dg d dg Z.e,e-e. dZ/g dZ0g dg d dg Z1e0e-e1 dZ2h dZ3g Z4g dddZ5ddgg ddZ6e7e5e6gddgD ]@\Z8Z9ej:f i e8de9iZ;ej<e;e9dZ=e4>e=e? f qd d! Z@ejABd"e@ ejABd#d$gd%d& ZCejABd"e@ d'd( ZDd)d* ZEd+d, ZFejAjGejABd-eF d.d/ ZHejABd0e#d1d2 ZIejABd3g d4d5d6 ZJd7d8 ZKd9d: ZLd;d< ZMd=d> ZNd?d@ ZOdAdB ZPdCdD ZQdEdF ZRdGdH ZSdIdJ ZTdKdL ZUdMdN ZVdOdP ZWdQdR ZXdSdT ZYdUdV ZZdWdX Z[dYdZ Z\d[d\ Z]d]d^ Z^d_d` Z_dadb Z`dcdd Zadedf Zbdgdh Zcdidj Zddkdl ZeejABdmg dnejABdoe#dpdq Zfdrds Zgdtdu Zhdvdw Zidxdy ZjejAjkdzd{ Zld|d} Zmd~d Zndg gdd ggdd ggdd ggdd ggdd ggdd ggdd ggdd ggddggdg gdg gdg gdg gdg gdg ggZoddggdg gddgggZpejABdeoep dd Zqdd ZrdS )    N)raises)IntegrationWarning)stats   )check_normalizationcheck_momentcheck_mean_expectcheck_var_expectcheck_skew_expectcheck_kurt_expectcheck_entropycheck_private_entropycheck_entropy_vect_scalecheck_edge_supportcheck_named_argscheck_random_state_propertycheck_meth_dtypecheck_ppf_dtypecheck_cmplx_derivcheck_picklingcheck_rvs_broadcastcheck_freezingcheck_munp_expect)distcont)rv_continuous_frozen   l        )
recipinvgaussvonmiseskappa4vonmises_line
gausshypernorminvgaussgeninvgaussgenhyperbolic	truncnormtruncweibull_min)studentized_rangekstwoksoneZ
wrapcauchygenexpon)r&   r   r   r(   r'   r   r)   )exponpow	exponweibr    r)   halfgennormgompertz	johnsonsb	johnsonsur   r(   r'   	kstwobignmielkencfnctpowerlognorm	powernormr   	trapezoidr   r   levy_stablerv_histogram_instancer&   )argusr*   r+   r    r)   Zgenhalflogisticr,   r-   r.   r   r0   r   r6   Z
truncexponr   r   r&   )alpha	betaprimeZbradfordburrburr12cauchycrystalballffisk
foldcauchy
genextreme	genpareto
halfcauchyinvgamma	jf_skew_tkappa3levylevy_l
loglaplacelomaxr1   nakagamir2   
skewcauchyttukeylambda
invweibullrel_breitwigner)r#   r/   r(   r'   r3   paretor5   r4   pearson3ZMLEMM)r<   r*   r+   r    r)   r,   r-   r.   r/   r   r(   r'   r0   r7   r1   r2   ncx2r4   r5   rdistr   r6   truncparetor   r   r&   )r:   r;   r<   r=   r>   r?   r@   rA   rB   rC   rD   rE   rF   rG   rH   rI   rJ   rK   rL   r1   rM   r2   r3   rN   rO   rY   rQ   )r#   r/   r(   r'   rS   r5   r4   >4   r!   betar3   r$   r8   rR   r2   ZgengammadweibullrB   rT   r;   dgammar&   rF   ricer0   rW   r   r'   rX   r.   rP   r    
reciprocalZ
loguniformrO   erlangchiskewnormr"   r%   r@   rJ   ZcosinerD   loggammar   rM   r9   Zlogisticr,   rG   ZmaxwellZhalflogisticr#   gammagennormrE   Zchi2r(   Zpowerlaw)r      re      rf   rf      rg   rg   rg   r   r   r   r   r      rh   rh   rh      ri   ri      rj   	   rj   )aZbins)r   r   
   TFdensity)rn   c                  c   st   t d d  t D ]^\} }| dkr$qq| tv rDtj| |tjjdV  q| tv rdtj| |tjjdV  q| |fV  qd S )Nr7   marks)	r   histogram_test_instancesdistslowpytestparammarkslow	distxslowxslow)distnamearg r{   c/var/www/html/assistant/venv/lib/python3.9/site-packages/scipy/stats/tests/test_continuous_basic.pycases_test_cont_basic   s    r}   zdistname,argzsn, n_fit_samples)i     c                 C   s  zt t| }W n ty(   | }d} Y n0 tjd}|j|||d}|j| \}}| dvrjt||| t|||  t	|||  t
|||  t|||  t|||  t|||  t|||  t|||  t|||  t|||  d}	| dkrt|j||	| n| dkrt| ||	| d}
|j|j|j|j|jg}ddd	d
ddd}|| d}| dkrdd}n| dkrrd}t||||
| t|| | dv rtrtd n
t|| t|| | dvrt |||  |j!dkrt"|| |j#j$tj%j$kr| dkrt&||tj% t'( >}|)t*d |)t*d |)t+d t,|| W d    n1 s`0    Y  t-|| t.|| t/||| t0|| | t1vrt2|| | dkrt3|||  dD ]N}| t4| vrt5|||d | | | t6| vrt7|||d | | qd S )Nr8   i )sizerandom_state>   Zlaplace_asymmetric{Gz?r"   r   r               ?g	@333333?      @)Zweibull_maxrJ   rS   rY   rP   r8         ?rQ   r   r(   )rf   rR   fails on Linux 32-bit)r0   r'   r2   r   r   z The occurrence of roundoff errorzExtremely bad integrandzinvalid valuer$   rU   )8getattrr   	TypeErrornprandomRandomStatervscheck_sample_meanvar_check_cdf_ppfcheck_sf_isfcheck_cdf_sfcheck_ppf_isf	check_pdfcheck_pdf_logpdfcheck_pdf_logpdf_at_endpointscheck_cdf_logcdfcheck_sf_logsfcheck_ppf_broadcastcheck_distribution_rvscdfpdflogpdflogcdflogsfgetr   r   	_IS_32BITrs   skipr   r   r   numargscheck_vecentropy	__class___entropyZrv_continuousr   nptsuppress_warningsfilterr   RuntimeWarningr   check_retrieving_supportr   r   r   fails_cmplxr   check_ppf_privateskip_fit_testcheck_fit_argsskip_fit_fix_testcheck_fit_args_fix)ry   rz   ZsnZn_fit_samplesdistfnrngr   mvr:   Zlocscale_defaultsZmethsZspec_xxsupmethodr{   r{   r|   test_cont_basic   s    










*





r   c                 C   sv   zt t| }W n ty(   | }d} Y n0 t|j| s>J t|j|ddisXJ t|j|dd isrJ d S )Nr8   r   r{   )r   r   r   r   Zisscalarr   )ry   rz   r   r{   r{   r|   test_rvs_scalar   s    
r   c                   C   s   t tjd d S )N)r   皙?)r   r   r7   r{   r{   r{   r|   &test_levy_stable_random_state_property  s    r   c            	      c   s   t  } dh}dh}td d  t D ]\}}|dkr4q"|tv rbtj||ddddtjjdddV  q"|| v}||v}||v}t }tj|||||d|dV  |r|r|s"tj||ddddtjj	g| dV  q"d S )	Nr2   r/   r7   Ttoo slow)reasonro   F)
setr   rq   distxslow_test_momentsrs   rt   ru   rx   listZxfail)	Zfail_normalizationZfail_higherZfail_momentry   rz   Zcond1Zcond2Zcond3rp   r{   r{   r|   cases_test_moments  s,    	r   z=distname,arg,normalization_ok,higher_ok,moment_ok,is_xfailingc              	   C   sZ  zt t| }W n ty(   | }d} Y n0 t }|td |td |td |rh|t |j|ddi\}}	}
}tjddr |rt	|||  |rt
||||  t||||	|
|  t||||	|  t||||	||  t|||  W d    n1 s0    Y  t||||	|  |r6t||||	|  W d    n1 sL0    Y  d S )	Nr8   z9The integral is probably divergent, or slowly convergent.z#The maximum number of subdivisions.z The algorithm does not converge.momentsmvskignoreall)r   r   r   r   r   r   r   r   errstater   r   r
   r	   r   r   check_loc_scaler   )ry   rz   Znormalization_okZ	higher_okZ	moment_okZis_xfailingr   r   r   r   skr{   r{   r|   test_moments2  s<    

,r   zdist,shape_argsc                 C   s   | dv rt d | dv r(tr(t d | dv }tt| }td}td}|j}g }ddg}t	|D ]B}	|	d	 fd
|	d   }
|
||	 t|
  |d|	d	  qh|||g t|| |||d d S )N)r    r&   r   r   r   )r9   r;   r\   r[   Z	exponnormr#   r"   r7   r3   r!   r]   ra   Zsemicircularrd   rb   re   )rf   r   rf   rg   r   r   d)rs   r   r   r   r   r   zerosonesr   rangeappendinsertextendr   )distZ
shape_argsZ
shape_onlyZdistfunclocscalenargsZallargsZbshaper   Zshpr{   r{   r|   test_rvs_broadcast\  s"    




r   zx,n,sf,cdf,pdf,rtol))gh㈵> ʚ;g$?gm%?g0eD@g+<)g&.!>r   gי?gscl)D>g ^RUU!@vIh%,=)gMb@?r   g_H'9$      ?gL.'r   c           	      C   s   t j}t|| ||| ||| |g}t|||g}tj|||d t	|d |d  d tj|
||g| gdd d S )Nrtolr   r   r   :0yE>)r   r(   r   arraysfr   r   r   assert_allcloseassert_equalisf)	r   nr   r   r   r   ksvalsexpectedr{   r{   r|   test_gh17775_regression  s    (r   c                  C   st  t jd} tjjt dd| d}t |}t	t 
|dkd tjjdt d| d}t |}t	t 
|dkd tjjt dt d| d}t |}t	t 
|dkd tjjt dgdggt d| d}t | }t	t 
|dkd tttjjddgddggddgddggd tttjjg dddd tttjjg d	g d
dgdggd d S )N{   r   r   )r   r   r   r   z0All the values are equal, but they shouldn't be!)re   rf   rg   r   )re   re   )r   r   r   r   )r   r   r   r   re   )rg   )r   r   r   r   normr   r   diffr   assert_r   r   r   Zravelassert_raises
ValueErrorrc   )r   r   r   r{   r{   r|   test_rvs_gh2069_regression  s,    


r   c                  C   s  t j} t| ddtjd t| dgdgtjgd t| dtj dd t| dgtj gdgd t| dgdgtjgd t| ddtjd t| dgtj gdgd t| dtj dd t| dtj dd d S )Nr   r   gA?g|N?)r   r$   r   assert_almost_equalr   r   infZ_cdf)tnr{   r{   r|   test_nomodify_gh9900_regression  s8    r   c                  C   s   t g d} t g d}dd t| |D }ttjd| || dd t| |D }ttjd| || dd t| |D }ttj| | || d	d t| |D }ttj|| || d S )
N)r   re   rf   rg   r   rh   )rj      r       r   0   c                 S   s    g | ]\}}t jd ||qS )ri   r   r^   r   .0_a_br{   r{   r|   
<listcomp>      z4test_broadcast_gh9990_regression.<locals>.<listcomp>ri   c                 S   s    g | ]\}}t jd ||qS r   r   r   r{   r{   r|   r     r   r   c                 S   s    g | ]\}}t j|||qS r{   r   r   r{   r{   r|   r     r   c                 S   s    g | ]\}}t j|||qS r{   r   r   r{   r{   r|   r     r   )r   r   zipr   assert_array_almost_equalr   r^   r   )rl   bansr{   r{   r|    test_broadcast_gh9990_regression  s    r  c                   C   sB   t jjtg ddtg d d tjtg ddd d S )N)      @       @r   r   )g      @r   g      @r  )rl   r  r   r   )r   r$   r   r   r   r   r{   r{   r{   r|    test_broadcast_gh7933_regression  s    r  c                  C   sb   t jddd } t jt d t df }dd t| |D }tj| |t j}t	|| d S )Nre   y             @Y@2   3   c                 S   s"   g | ]\}}t j||tjqS r{   )r   r$   r   r   r   )r   Z_xr   r{   r{   r|   r     r   z*test_gh2002_regression.<locals>.<listcomp>)
r   Zr_r   r  r   r$   r   r   r   r  )r   rl   r   r  r{   r{   r|   test_gh2002_regression  s
    r  c                  C   s(   d} t jdt| g| d gg d S )Ng(\@r   )r   rC   ppfr   r   )cr{   r{   r|   test_gh1320_regression  s    r  c                  C   s   t jd g d} ddt d d  }ddt d d  }tttjtjj| dd\}}t	j
||d	d
 t	j
|| |d	d
 d S )Ni  )r   r   r   r   r   皙?re   rf   r   rV   r   rg   )decimal)r   r   seedsqrtsupertyper   uniformfitr   r   )r   rl   r  r   r   r{   r{   r|   test_method_of_moments&  s     r  c                 C   s,   t | rt||  t |r(t|| d S N)r   isfinitecheck_sample_meancheck_sample_var)popmeanpopvarsampler{   r{   r|   r   2  s    


r   c                 C   s   t | |j}|dksJ d S )Nr   )r   Zttest_1sampZpvalue)r   r  Zprobr{   r{   r|   r  9  s    r  c                 C   sH   t j| fdd dd}|j}|j|j }}||  kr>|ksDn J d S )Nc                 S   s   | j d|dS )Nr   )Zddofaxis)var)r   r!  r{   r{   r|   <lambda>E  r   z"check_sample_var.<locals>.<lambda>gףp=
?)Zconfidence_level)r   Z	bootstrapZconfidence_intervallowhigh)r   r  resZconfr$  r%  r{   r{   r|   r  ?  s    r  c                 C   s>   g d}t j| j| j|g|R  g|R  |t|d d d S )N)gMbP?r   g+?z - cdf-ppf roundtripr  err_msg)r   r   r   r  DECIMAL)r   rz   msgvaluesr{   r{   r|   r   M  s    "r   c                 C   s>   t j| j| jg dg|R  g|R  g dt|d d d S )Nr   r   ?z - sf-isf roundtripr'  )r   r   r   r   r)  r   rz   r*  r{   r{   r|   r   T  s
    &
r   c                 C   sB   t j| jddgg|R  d| jddgg|R   t|d d d S )Nr   r-  r   z - cdf-sf relationshipr'  )r   r   r   r   r)  r.  r{   r{   r|   r   Z  s    r   c                 C   sH   t ddg}tj| j|g|R  | jd| g|R  t|d d d S )Nr   r-  r   z - ppf-isf relationshipr'  )r   r   r   r   r   r  r)  )r   rz   r*  pr{   r{   r|   r   a  s    (r   c                 C   s   | j dg|R  }d}| j|g|R  }|dk s8|dkrR|d }| j|g|R  }| j|| g|R  | j|| g|R   | d }|d7 }tj||t|d d S )	Nr   gư>g-C6?     @r   r  z - cdf-pdf relationshipr'  )r  r   r   r   r   r)  )r   rz   r*  medianepsZpdfvZcdfdiffr{   r{   r|   r   h  s    r   c                 C   s   t g d}| j|g|R  }|t | }| j|g|R  }| j|g|R  }||dkt |@  }|t | }|d7 }tjt ||d|d d S )N)r  r   皙?r   333333?ffffff?皙?r    - logpdf-log(pdf) relationshipri   r'  	r   r   r  r  r   r   r   r   logr   argsr*  pointsr   r   r   r{   r{   r|   r   z  s    r   c                 C   s   t ddg}| j|g|R  }|t | }| j|g|R  }| j|g|R  }||dkt |@  }|t | }|d7 }tjt ||d|d d S )Nr   r   r7  ri   r'  r8  r:  r{   r{   r|   r     s    r   c                 C   s   t g d}| j|g|R  }|t | }| j|g|R  }| j|g|R  }||dk }|t | }|d7 }tjt ||d|d d S )N)	        r  r   r3  r   r4  r5  r6  r   r   z - logsf-log(sf) relationshipri   r'  )	r   r   r  r  r   r   r   r   r9  )r   r;  r*  r<  r   r   r   r{   r{   r|   r     s    r   c                 C   s   t g d}| j|g|R  }|t | }| j|g|R  }| j|g|R  }||dk }|t | }|d7 }tjt ||d|d d S )N)	r   r  r   r3  r   r4  r5  r6  r   r   z - logcdf-log(cdf) relationshipri   r'  )	r   r   r  r  r   r   r   r   r9  )r   r;  r*  r<  r   r   r   r{   r{   r|   r     s    r   c                    sj   d g   }|r" fdd|D }| j dg|R  }| j dg|R  }|d7 }tj||g  d|d d S )Nr   c                    s   g | ]}t |g  qS r{   )r   r   )r   _Znum_repeatsr{   r|   r     r   z'check_ppf_broadcast.<locals>.<listcomp>r   z - ppf multipleri   r'  )r  r   r   )r   rz   r*  r;  r1  Zmediansr{   r?  r|   r     s    r   c                 C   st   t j|| |dd\}}||k rpt j| | |dd\}}t||kdt| d t| d t| d t|  d S )Ni  )r;  NzD = z	; pval = z
; alpha = z
args = )r   Zkstestr   r   str)r   r;  r:   r   DZpvalr{   r{   r|   r     s     r   c                 C   s   t | j| | j|  d S r  )r   r   Z
vecentropyr   )r   r;  r{   r{   r|   r     s    r   c           	      C   s`   t ddgt ddg }}| j|||d\}}t|| | | t|| | | d S )N      $@g      4@r   r   )r   r   r   r   r   )	r   rz   r   r   r*  r   r   mtvtr{   r{   r|   r     s    r   c                 C   s>   | j tg dg|R  }ttt| |d  d S )Nr,  zppf private is nan)Z_ppfr   r   r   r   anyisnan)r   rz   r*  Zppfsr{   r{   r|   r     s    r   c                 C   sH   d\}}| j | }| j |||d}tt|| | t| d S )N)r   re   rD  )supportr   r   r   r   )r   r;  r   r   ZsuppZsupp_loc_scaler{   r{   r|   r     s    
r   c              	   C   s   t jddr t J}|jtdd |jtdd | j||d}| j|d|d}W d    n1 sf0    Y  W d    n1 s0    Y  tt|d	t| k tt|d	t| k d S )
Nr   r   !The shape parameter of the erlangcategorymessagezfloating point number truncatedr  Zpowell)Z	optimizerr   re   	r   r   r   r   r   r   r  r   len)r   rz   r   r   r   r   vals2r{   r{   r|   r     s    Lr   c           
   	   C   s  t jdd t }|jtdd | j|d|d}| j|d|d}tt|d	t| k t|d
 dk t|d dk tt|d	t| k t|dkr| j||d |d}tt|d	t| k t|d |d k t|dkrD| j||d |d}tt|d	t| k t|d |d k t|d	kr| j||d	 |d}	tt|	d	t| k t|	d	 |d	 k W d    n1 s0    Y  W d    n1 s0    Y  d S )Nr   r   rJ  rK  r   )Zflocr   r   )Zfscaler   re   r	  r   )Zf0r   )f1r   )f2r   rN  )
r   rz   r   r   r   r   rP  vals3vals4vals5r{   r{   r|   r     s,    r   r   )r   r   r   r   r   r   r  r   zdistname, argsc           	         s   t t|}t ||  |dkr0| dr0ddg}nddg}dd |D }d	dg}d
dg} |g|R ||d}tj| fddt|g|||R  D ddd d S )NrQ   r9  r   re   r   r  c                 S   s   g | ]}|gd  qS )re   r{   )r   rl   r{   r{   r|   r     r   z+test_methods_with_lists.<locals>.<listcomp>r   r   g)\(?rD  c                    s   g | ]} | qS r{   r{   )r   r   r@   r{   r|   r   !  r   g+=r   )r   Zatol)r   r   
startswithr   r   r  )	r   ry   r;  r   r   Zshape2r   r   resultr{   rV  r|   test_methods_with_lists  s    


$rY  c                  C   s>   t jddd} t| tsJ t jdd}t|ts:J d S )Nr   r   rg   rj   )r   r<   moment
isinstancefloatrA   )vals0vals1r{   r{   r|   (test_burr_fisk_moment_gh13234_regression%  s    r_  c               	   C   s  t jjdtg ddd} tg d}t| | t jjdtg ddd}ttjtjtjg}t|| t jjdtg dg dd}ttjdtjg}t|| t jjdd	d
d}tj}t|| t||j	sJ t jjdg dg dd}tdtjtjg}	t||	 t jjdg dg dd}
tdtjdg}t|
| t jjdg dg dd}ttjtjtjg}t|| t j
jddd	d	d}tj}t|| t||j	sJ t j
jdg dd	d	d}ttjtjtjg}t|| t j
jdg dg dg dd}tdtjtjg}tj||dd t jdddgddg}tddg}tj||dd dd	ddtjg}dd	ddtjg}dd	tdg}d d	ddtjg}t|d!}t|d"}t|d#}t|}t jjd||||d$}t||||\}}}}t|D ]v}tjd%d%d&P |d	 }t jjd|| || || || d$}tj|| | W d    n1 s0    Y  q4d S )'Nr   r   re   rf   orderr   r   r   r  r  r   )r   r   r  re   r   )r   r   re   )rf   re  rC  )r   r   r   )r   r	  g      Y@g      9@r0  )g      r	  g      Yrb  Zdfr   r   )r   r=  r  )r   g      r=  g ﲙ"@r   r   r   r   g      :@g      @gvIh%<=gg@g )r   r   r   r   )r   r   r   )r   r   )rb  rl   r  r   r   r   )invaliddivide)r   r   rZ  r   r   r   r   nanr[  r   r`   r   pir  ZreshaperZ   Zbroadcast_arraysZndenumerater   testing)r]  	expected0r^  	expected1rP  	expected2rS  	expected3rT  Z	expected4rU  Z	expected5Zvals6Z	expected6Zvals7Z	expected7Zvals8Z	expected8Zvals9Z	expected9Zvals10Z
expected10rl   r  r   r   Zvals11ir   r{   r{   r|   *test_moments_with_array_gh12192_regression-  sp    
rr  c                  C   sn  t jjdtg ddggd} tg dg}t| | | j|jksLJ t jjdtdgdgdggg dd}tg dg dg d	g}t|| |j|jksJ t jjdg dd
dd}tg d}tj	||dd |j|jksJ t jjddgdgdggg dg dd}ttj
tj
dgtj
tj
dgtj
tj
dgg}tj	||dd |j|jksjJ d S )Nr   r`  ra  rc  re   rf   )r   r   r   )r  r  r  )r  r  r  r=  r   rg  )g e3E?g> ?g e3E?r   r   r  r  )r=  r   r  )g      r=  r  g@3F@gw,W.
@gf3&@)r   r   rZ  r   r   r   r   shaper`   r   rj  )r]  rm  r^  rn  rP  ro  rS  rp  r{   r{   r|   /test_broadcasting_in_moments_gh12192_regression  s.     rt  c                     sJ   g d d} t  fdd| D }t tjj | d}t|| d S )N)r   r         @      @      @r   c                    s   g | ]  fd dD qS )c                    s   g | ]}t jj | d qS r   )r   rH   )r   rs  rZ  r{   r|   r     s   z8test_kappa3_array_gh13582.<locals>.<listcomp>.<listcomp>r{   r   Zshapesrz  r|   r     s   
z-test_kappa3_array_gh13582.<locals>.<listcomp>ry  )r   r   r   rH   r   r   r   r&  Zres2r{   r|  r|   test_kappa3_array_gh13582  s    r~  c                     s  t g d t g dd} t  fdd| D }t tjj | d}t|| t g d t g dt  fd	d| D }t tjj | d}t|| t g d
 t g dd d d f t tjj | d}|jdksJ d S )N)r   ru  rv  rw  rd  )r   r   g      r   rv  r   c                    s&   g | ]  fd dt dD qS )c                    s&   g | ]}t jj  | | d qS rx  r   r   r   rq  hr   rZ  r{   r|   r     s   8test_kappa4_array_gh13582.<locals>.<listcomp>.<listcomp>r   r   r{  r  r   rz  r|   r     s   z-test_kappa4_array_gh13582.<locals>.<listcomp>ry  )r         пr  r   r   r   )r   r   r   gUUUUUUտr   r   c                    s&   g | ]  fd dt dD qS )c                    s&   g | ]}t jj  | | d qS rx  r  r  r  r{   r|   r     s   r  rh   r  r{  r  rz  r|   r     s   )r   r   r   )r   r   r   r   )rg   rg   rf   )r   r   r   r   r   r   rs  r}  r{   r  r|   test_kappa4_array_gh13582  s&    r  c                  C   s   d} t jt| d t j W d    n1 s20    Y  t jt| d t j W d    n1 sj0    Y  dtj_t }t|tsJ t	tjd d S )Nz.'rv_continuous_frozen' object has no attribute)matchZherringpmf)
rs   r   AttributeErrorr   r   r  Zlogpmfr[  r   delattr)rM  Zfrozen_normr{   r{   r|   test_frozen_attributes  s    ((r  c                  C   s   t jd} t j d }}ttj||tj	| | jdd| jdd }}| jdddk }d||< || }tj||}t|| tj	| t||  tj||  ||   d S )Nr   )rf   rf   )r   r   )
r   r   Zdefault_rngr   r   r   r   ra   r   r   )r   r   rl   maskZx_normr&  r{   r{   r|   test_skewnorm_pdf_gh16038  s    r  r   r   r   r   r   r   r   r  r   rZ  Zentropyexpectr1  meanZstdr"  intervalgffffff?rI  r   mvcasec                 C   sd   | \}}t t |}|| }| tv r8t|tjs`J n(t|d tjsLJ t|d tjs`J d S )Nr   r   )r   r   r   
scalar_outr[  r   number)r  method_namer;  r   r&  r{   r{   r|   test_scalar_for_scalar  s    r  c                  C   s   t jg d} t| d tjs$J t| d tjs8J t jg d} t| d tjs\J t| d tjspJ t jdg d} t| tjsJ d S )Nr`  r   r   r   )r   r   r  r[  r   r  Zfit_loc_scaleZnnlf)r&  r{   r{   r|   test_scalar_for_scalar2  s    r  )ssysnumpyr   Znumpy.testingrl  r   rs   r   r   Zscipy.integrater   	itertoolsZscipyr   Zcommon_testsr   r   r   r	   r
   r   r   r   r   r   r   r   r   r   r   r   r   r   r   Zscipy.stats._distr_paramsr   Z!scipy.stats._distn_infrastructurer   r)  maxsizer   rr   rw   r   Zskip_fit_test_mleZslow_fit_test_mmZfail_fit_test_mmr   Zskip_fit_fix_test_mleZfail_fit_fix_test_mmr   r   rq   Zcase1Zcase2productr  rn   Z	histogramZ_histZrv_histogramZ_rv_histr   tupler}   ru   Zparametrizer   r   r   r   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   rY  r_  rr  rt  r~  rx   r  r  r  r  Zscalars_outr  r  r{   r{   r{   r|   <module>   s   T
		

	^
%&
5

$
	T

 
