o
    '&i                     @   sR  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:di 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 ZcejABdig djejABdke#dldm Zddndo Zedpdq Zfdrds Zgdtdu ZhejAjidvdw Zjdxdy Zkdzd{ Zld|g gd}d ggd~d ggdd 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Zmddggdg gddgggZnejABdemen dd Zodd Zpd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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   r1   rL   r2   r3   rM   rN   rX   rP   )r#   r/   r(   r'   rR   r5   r4   >3   ricer'   r    r0   rW   r   rD   
reciprocalrN   r(   rO   gennormrS   r@   Zgengammar%   rF   dgammaZhalflogisticZcosinerV   rI   r9   r"   rQ   r$   r8   r,   r;   loggammar2   r3   Zmaxwellchir   r.   betaskewnormdweibullZchi2r!   ZpowerlawZlogisticZ
loguniformgammaerlangr&   rL   r#   rE   rB   )r      rd      re   re      rf   rf   rf   r   r   r   r   r      rg   rg   rg      rh   rh      ri   	   ri   )aZbins)r   r   
   TFdensity)rm   c                  c   st    t d d  t D ].\} }| dkrq	| tv r"tj| |tjjdV  q	| tv r2t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 rz   RC:\wamp64\www\opt\env\Lib\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w tjd}|j|||d}|j| \}}| dvr5t||| t|||  t	|||  t
|||  t|||  t|||  t|||  t|||  t|||  d}	| dkrtt|j||	| n| dkrt| ||	| d}
|j|j|j|j|jg}ddd	d
ddd}|| d}| dkrd}n| dkrd}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rt$||tj# t%& }|'t(d |'t(d |'t)d t*|| W d    n	1 sw   Y  t+|| t,|| t-||| t.|| | t/vr>t0|| | dkrIt1|||  dD ]'}| t2| vr_t3|||d | | | t4| vrqt5|||d | | qKd S )Nr8   i )sizerandom_state>   Zlaplace_asymmetric{Gz?r"   r   r               ?g	@333333?      @)Zweibull_maxrI   rR   rX   rO   r8         ?rP   r   r(   )re   rQ   fails on Linux 32-bit)r0   r'   r2   r   r   z The occurrence of roundoff errorzExtremely bad integrandzinvalid valuer$   rT   )6getattrr   	TypeErrornprandomRandomStatervscheck_sample_meanvar_check_cdf_ppfcheck_sf_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_32BITrr   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)rx   ry   ZsnZn_fit_samplesdistfnrngr   mvr:   Zlocscale_defaultsZmethsZspec_xxsupmethodrz   rz   r{   test_cont_basic   s   











r   c                 C   sv   zt t| }W n ty   | }d} Y nw t|j| sJ t|j|ddis,J t|j|dd is9J d S )Nr8   r~   rz   )r   r   r   r   Zisscalarr   )rx   ry   r   rz   rz   r{   test_rvs_scalar   s   r   c                   C   s   t tjd d S )N)r   皙?)r   r   r7   rz   rz   rz   r{   &test_levy_stable_random_state_property  s   r   c            	      c   s    t  } dh}dh}td d  t D ]T\}}|dkrq|tv r2tj||ddddtjjdddV  q|| v}||v}||v}t }tj|||||d|dV  |rT|rT|sftj||ddddtjj	g| dV  qd S )	Nr2   r/   r7   Ttoo slow)reasonrn   F)
setr   rp   distxslow_test_momentsrr   rs   rt   rw   listZxfail)	Zfail_normalizationZfail_higherZfail_momentrx   ry   Zcond1Zcond2Zcond3ro   rz   rz   r{   cases_test_moments  s2   	r   z=distname,arg,normalization_ok,higher_ok,moment_ok,is_xfailingc              	   C   sd  zt t| }W n ty   | }d} Y nw t }|td |td |td |r3|t |j|ddi\}}	}
}tjdd9 |rNt	|||  |rwt
||||  t||||	|
|  t||||	|  t||||	||  t|||  W d    n1 sw   Y  t||||	|  |rt||||	|  W d    d S W d    d S 1 sw   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   )rx   ry   Znormalization_okZ	higher_okZ	moment_okZis_xfailingr   r   r   r   skrz   rz   r{   test_moments0  sH   

"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 ]!}	|	d	 fd
|	d   }
|
||	 t|
  |d|	d	  q4|||g t|| |||d d S )N)r    r&   r   r   r   )r9   r;   r\   ra   Z	exponnormr#   r"   r7   r3   r!   rY   r`   Zsemicircularr[   r]   rd   )re   r   re   rf   r   r   d)rr   r   r   r   r   r   zerosonesr   rangeappendinsertextendr   )distZ
shape_argsZ
shape_onlyZdistfunclocscalenargsZallargsZbshaper   Zshprz   rz   r{   test_rvs_broadcastZ  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expectedrz   rz   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!)rd   re   rf   r   )rd   rd   )r   r   r   r   )r   r   r   r   rd   )rf   )r   r   r   r   normr   r   diffr   assert_r   r   r   Zravelassert_raises
ValueErrorrb   )r   r   r   rz   rz   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)tnrz   rz   r{   test_nomodify_gh9900_regression  s     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   rd   re   rf   r   rg   )ri      r       r   0   c                 S       g | ]\}}t jd ||qS )rh   r   rZ   r   .0_a_brz   rz   r{   
<listcomp>       z4test_broadcast_gh9990_regression.<locals>.<listcomp>rh   c                 S   r   r   r   r   rz   rz   r{   r     r   r   c                 S   s    g | ]\}}t j|||qS rz   r   r   rz   rz   r{   r     r   c                 S   s    g | ]\}}t j|||qS rz   r   r   rz   rz   r{   r     r   )r   r   zipr   assert_array_almost_equalr   rZ   r   )rk   bansrz   rz   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  )rk   r  r   r   )r   r$   r   r   r   r   rz   rz   rz   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 )Nrd   y             @Y@2   3   c                 S   s"   g | ]\}}t j||tjqS rz   )r   r$   r   r   r   )r   Z_xr   rz   rz   r{   r   
  s   " z*test_gh2002_regression.<locals>.<listcomp>)
r   Zr_r   r   r   r$   r   r   r   r   )r   rk   r   r  rz   rz   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   )crz   rz   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   皙?rd   re   r   rU   r   rf   )decimal)r   r   seedsqrtsupertyper   uniformfitr   r   )r   rk   r  r   r   rz   rz   r{   test_method_of_moments  s    r  c                 C   s0   t | r
t||  t |rt|| d S d S N)r   isfinitecheck_sample_meancheck_sample_var)popmeanpopvarsamplerz   rz   r{   r   '  s
   


r   c                 C   s   t | |j}|dksJ d S )Nr   )r   Zttest_1sampZpvalue)r  r  Zprobrz   rz   r{   r  .  s   r  c                 C   sJ   t j| fdd dd}|j}|j|j}}||  kr |ks#J  J d S )Nc                 S   s   | j d|dS )Nr   )Zddofaxis)var)r   r  rz   rz   r{   <lambda>:  s    z"check_sample_var.<locals>.<lambda>gףp=
?)Zconfidence_level)r   Z	bootstrapZconfidence_intervallowhigh)r  r  resconfr"  r#  rz   rz   r{   r  4  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   ry   msgvaluesrz   rz   r{   r   B  s   "
r   c                 C   s|   t j| j| jg dg|R  g|R  g dt|d d t j| jddgg|R  d| jddgg|R   t|d d d S )Nr   r   ?z - sf-isf roundtripr&  r   r,  r   z - cdf-sf relationship)r   r   r   r   r(  r   )r   ry   r)  rz   rz   r{   r   I  s   &

r   c                 C   s   | j dg|R  }d}| j|g|R  }|dk s|dkr)|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   ry   r)  medianZepsZpdfvZcdfdiffrz   rz   r{   r   S  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) relationshiprh   r&  	r   r   r  r  r   r   r   r   logr   argsr)  pointsr   r   r   rz   rz   r{   r   e     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   r3  rh   r&  r4  r6  rz   rz   r{   r   r  r9  r   c                 C      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/  r   r0  r1  r2  r   r   z - logsf-log(sf) relationshiprh   r&  )	r   r   r  r  r   r   r   r   r5  )r   r7  r)  r8  r   r   r   rz   rz   r{   r        r   c                 C   r:  )N)	r   r  r   r/  r   r0  r1  r2  r   r   z - logcdf-log(cdf) relationshiprh   r&  )	r   r   r  r  r   r   r   r   r5  )r   r7  r)  r8  r   r   r   rz   rz   r{   r     r<  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 rz   )r   r   )r   _Znum_repeatsrz   r{   r     s    z'check_ppf_broadcast.<locals>.<listcomp>r   z - ppf multiplerh   r&  )r  r   r   )r   ry   r)  r7  r.  Zmediansrz   r>  r{   r     s   r   c                 C   sx   t j|| |dd\}}||k r:t j| | |dd\}}t||kdt| d t| d t| d t|  d S d S )Ni  )r7  NzD = z	; pval = z
; alpha = z
args = )r   Zkstestr   r   str)r   r7  r:   r   DZpvalrz   rz   r{   r     s    
r   c                 C   s   t | j| | j|  d S r  )r   r   Z
vecentropyr   )r   r7  rz   rz   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   ry   r   r   r)  r   r   mtvtrz   rz   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   ry   r)  Zppfsrz   rz   r{   r     s   "r   c                 C   sH   d\}}| j | }| j |||d}tt|| | t| d S )N)r   rd   rC  )supportr   r   r   r   )r   r7  r   r   ZsuppZsupp_loc_scalerz   rz   r{   r     s   
r   c              	   C   s   t jdd9 t %}|jtdd |jtdd | j||d}| j|d|d}W d    n1 s3w   Y  W d    n1 sBw   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   rd   	r   r   r   r   r   r   r  r   len)r   ry   r   r   r   r   vals2rz   rz   r{   r     s    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rw| 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 t|d	kr| j||d	 |d}	tt|	d	t| k t|	d	 |d	 k W d    n1 sw   Y  W d    d S W d    d S 1 sw   Y  d S )Nr   r   rI  rJ  r   )Zflocr   r   )Zfscaler   rd   r  r   )Zf0r   )f1r   )f2r   rM  )
r   ry   r   r   r   r   rO  vals3vals4vals5rz   rz   r{   r     s0   Pr   r   )r   r   r   r   r   r   r  r   zdistname, argsc           	         s   t t|}t ||  |dkr| drd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 )NrP   r5  r   rd   r   r  c                 S   s   g | ]}|gd  qS )rd   rz   )r   rk   rz   rz   r{   r     s    z+test_methods_with_lists.<locals>.<listcomp>r   r   g)\(?rC  c                    s   g | ]} | qS rz   rz   )r   r   r@   rz   r{   r     s    g+=r   )r   Zatol)r   r   
startswithr   r   r   )	r   rx   r7  r   r   Zshape2r   r   resultrz   rU  r{   test_methods_with_lists  s   


$
rX  c                  C   s>   t jddd} t| tsJ t jdd}t|tsJ d S )Nr   r   rf   ri   )r   r<   moment
isinstancefloatrA   )vals0vals1rz   rz   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	stJ 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 ];}tjd%d%d&( |d	 }t jjd|| || || || d$}tj|| | W d    n	1 sw   Y  qd S )'Nr   r   rd   re   orderr   r   r   r  r  r   )r   r   r  rd   r   )r   r   rd   )re   rd  rB  )r   r   r   )r   r  g      Y@g      9@r-  )g      r  g      Yra  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   )ra  rk   r  r   r   r   )invaliddivide)r   r   rY  r   r   r   r   nanrZ  r   r^   r   pir  Zreshaper_   Zbroadcast_arraysZndenumerater   testing)r\  	expected0r]  	expected1rO  	expected2rR  	expected3rS  Z	expected4rT  Z	expected5Zvals6Z	expected6Zvals7Z	expected7Zvals8Z	expected8Zvals9Z	expected9Zvals10Z
expected10rk   r  r   r   Zvals11ir   rz   rz   r{   *test_moments_with_array_gh12192_regression  sv   
rq  c                  C   sl  t jjdtg ddggd} tg dg}t| | | j|jks&J 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sVJ t jjdg dd
dd}tg d}tj	||dd |j|jksyJ 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sJ d S )Nr   r_  r`  rb  rd   re   )r   r   r   )r  r  r  )r  r  r  r;  r   rf  )g e3E?g> ?g e3E?r   r   r  r  )r;  r   r  )g      r;  r  g@3F@gw,W.
@gf3&@)r   r   rY  r   r   r   r   shaper^   r   ri  )r\  rl  r]  rm  rO  rn  rR  ro  rz   rz   r{   /test_broadcasting_in_moments_gh12192_regressionl  s.    rs  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   rG   )r   rr  rY  rz   r{   r     s    z8test_kappa3_array_gh13582.<locals>.<listcomp>.<listcomp>rz   r   Zshapesry  r{   r     s
    

z-test_kappa3_array_gh13582.<locals>.<listcomp>rx  )r   r   r   rG   r   r   r   r$  Zres2rz   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sJ d S )N)r   rt  ru  rv  rc  )r   r   g      r   ru  r   c                    &   g | ]  fd dt dD qS )c                    &   g | ]}t jj  | | d qS rw  r   r   r   rp  hr   rY  rz   r{   r         8test_kappa4_array_gh13582.<locals>.<listcomp>.<listcomp>r   r   rz  r  r   ry  r{   r     
    
z-test_kappa4_array_gh13582.<locals>.<listcomp>rx  )r         пr  r   r   r   )r   r   r   gUUUUUUտr   r   c                    r~  )c                    r  rw  r  r  r  rz   r{   r     r  r  rg   r  rz  r  ry  r{   r     r  )r   r   r   )r   r   r   r   )rf   rf   re   )r   r   r   r   r   r   rr  r|  rz   r  r{   test_kappa4_array_gh13582  s&   r  c                  C   s   d} t jt| d t j W d    n1 sw   Y  t jt| d t j W d    n1 s5w   Y  dtj_t }t|tsIJ t	tjd d S )Nz.'rv_continuous_frozen' object has no attribute)matchZherringpmf)
rr   r   AttributeErrorr   r   r  ZlogpmfrZ  r   delattr)rL  Zfrozen_normrz   rz   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   )re   re   )r~   r   )
r   r   Zdefault_rngr   r   r   r   r`   r   r   )r   r   rk   maskZx_normr$  rz   rz   r{   test_skewnorm_pdf_gh16038  s   ,r  r   r   r   r   r   r   r   r  r   rY  Zentropyexpectr.  meanZstdr   intervalgffffff?rH  r   mvcasec                 C   sf   | \}}t t |}|| }| tv rt|tjsJ d S t|d tjs'J t|d tjs1J d S )Nr   r   )r   r   r   
scalar_outrZ  r   number)r  method_namer7  r   r$  rz   rz   r{   test_scalar_for_scalar  s   r  c                  C   s   t jg d} t| d tjsJ t| d tjsJ t jg d} t| d tjs.J t| d tjs8J t jdg d} t| tjsIJ d S )Nr_  r   r   r   )r   r   r  rZ  r   r  Zfit_loc_scaleZnnlf)r$  rz   rz   r{   test_scalar_for_scalar2  s   r  rz   )qsysnumpyr   Znumpy.testingrk  r   rr   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   rq   rv   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   rp   Zcase1Zcase2productr  rm   Z	histogramZ_histZrv_histogramZ_rv_histr   tupler|   rt   Zparametrizer   r   r   r   ru   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   rX  r^  rq  rs  r}  rw   r  r  r  r  Zscalars_outr  r  rz   rz   rz   r{   <module>   s    T
		

	\
%&
5



	T

 
