
    5[g                     j   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 jN                  d
k  Z(h dZ)h dZ*h dZ+h dZ,h dZ-h dZ.ddhZ/h dZ0h dZ1h dZ2h dZ3h dZ4dhZ5g Z6g dddZ7ddgg ddZ8 ejr                  e7e8gddg      D ]G  \  Z:Z; ejx                  dti e:de;iZ= ej|                  e=e;      Z?e6j                  e? eA       f       I d  ZBej                  j                  d! eB             ej                  j                  d"d#g      d$               ZEd% ZFd& ZGej                  j                  d' eF             ej                  j                  d(d)g      d*               ZHej                  j                  d! eB             d+        ZId, ZJd- ZKej                  j                  ej                  j                  d. eK             d/               ZMej                  j                  d0e#      d1        ZNej                  j                  d2g d3      d4        ZOd5 ZPd6 ZQd7 ZRd8 ZSd9 ZTd: ZUd; ZVd< ZWd= ZXd> ZYd? ZZd@ Z[dA Z\dB Z]dC Z^dD Z_dE Z`dF ZadG ZbdH ZcdI ZddJ ZedK ZfdL ZgdM ZhdN ZidO ZjdP Zkej                  j                  dQg dR      ej                  j                  dS ek             dT               ZldU ZmdV ZndW Zoej                  j                  dX        Zpej                  j                  dY        ZrdZ Zsd[ Ztd\g gd]d ggd^d ggd_d ggd`d ggdad ggdbd ggdcd ggddd ggdedggdfg gdgg gdhg gdig gdjg gdkg ggZudldmggdng gdodpgggZvej                  j                  dqeuevz         dr        Zwds Zxy)u    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        >
   ksonekstwokappa4vonmises
gausshypernorminvgaussrecipinvgaussvonmises_linerel_breitwignerstudentized_range>   pearson3crystalballpowerlognorm>   r   r   r   r    r#   r$   r&   >   r   genexpon	exponweib	johnsonsbtukeylambdar)   genhyperbolic>   ncfncx2r!   r#   r$   >   r   r   	irwinhall	trapezoidtruncparetolevy_stabler&   >   chi2exponlognormpowerlaw
loguniform
reciprocal>   betar0   argusexponpowgengammagompertzr,   	kstwobign	truncnormr!   
wrapcauchygeninvgausshalfgennormr"   genhalflogistictruncweibull_min>'   ftr/   nctburrfisklevyalphalomaxburr12cauchykappa3r   levy_lmielkeparetobradfordinvgammar    	betaprimer+   	genparetor1   	jf_skew_t	johnsonsu	powernormr2   
foldcauchy
genextreme
halfcauchy
loglaplace
skewcauchy
truncexponr(   r3   r-   r)   r$   r%   >   r   r   r*   r4   r.   r#   r&   >5   rG   rH   chir/   rI   r;   r5   r0   ricer<   gammar   r   rdistcosinedgammaerlangrR   gennormmaxwelldweibullr>   rV   loggammalogisticnakagamir'   r8   skewnormr    rW   rX   r1   rY   r,   r@   rA   r\   r!   r^   r9   r:   rC   rD   r-   halflogisticr"   r.   r$   r%   rF   r&   rv_histogram_instance)r      rr      rs   rs      rt   rt   rt   r   r   r   r   r      ru   ru   ru      rv   rv      rw   	   rw   )abins)r   r   
   TFdensity)r|   c               #   @  K   t         d d  t        z   D ]  \  } }| dk(  r| t        v r3t        j                  | |t        j
                  j                         G| t        v r3t        j                  | |t        j
                  j                         | |f  y w)Nr4   marks)	r   histogram_test_instancesslow_test_cont_basicpytestparammarkslowxslow_test_cont_basicxslowdistnameargs     b/var/www/html/bid-api/venv/lib/python3.12/site-packages/scipy/stats/tests/test_continuous_basic.pycases_test_cont_basicr   m   s     !!'??#}$++,,xFKK4D4DEE..,,xFKK4E4EFFC- @s   BBzdistname,argsni  c                 j   	 t        t        |       }t        j                  j                  d      } |j                  |||d} |j                  | \  }}| dvrt        |||       t        |||        t        |||        t        |||        t        |||        t        |||        t        |||        t        |||        t        |||        t!        |||        t#        |||        d}| dk(  rt%        |j&                  |||       n| dk7  rt%        | |||       d}	|j(                  |j*                  |j&                  |j,                  |j.                  g}
ddd	d
ddd}|j1                  | d      }| dk(  rd}n| dk(  rd}t3        ||||	|
       t5        ||       | dv rt6        rt9        j:                  d       nt=        ||       t?        ||       | dvrtA        |||        |jB                  dk(  rtE        ||       |jF                  jH                  t        jJ                  jH                  k7  r | dk7  rtM        ||t        jJ                         tO        jP                         5 }|jS                  tT        d       |jS                  tT        d       |jS                  tV        d       tY        ||       d d d        t[        ||       t]        ||       t_        |||
       ta        ||       | tb        vrte        ||       | dk7  rtg        |||        y y # t        $ r | }d} Y w xY w# 1 sw Y   vxY w)Nrq    sizerandom_state>   laplace_asymmetric{Gz?rC   r   r               ?g	@333333?      @)weibull_maxrR   rT   r3   r-   rq         ?
invweibullr   r   )rs   r%   fails on Linux 32-bit)r@   r   r/   r   r    z The occurrence of roundoff errorzExtremely bad integrandzinvalid valuerA   )4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_32BITr   skipr   r   r   numargscheck_vecentropy	__class___entropyrv_continuousr   nptsuppress_warningsfilterr   RuntimeWarningr   check_retrieving_supportr   r   r   fails_cmplxr   check_ppf_private)r   r   r   distfnrngr   mvrM   locscale_defaultsmethsspec_xxsups                 r   test_cont_basicr   y   s   +)
 ))


'C
&**s"
5C6<<DAq--aC(&#x(h'h'&#x(fc8$VS(+!&#x8VS(+63)X.E**vzz3s;	]	" 	xeS9ZZ

FMM\\E "TC'*,F 	

8S!A<	W	VQ%6>,&&9+,vs#63 44fc8,~~%!!U%8%8%A%AAJ&fc5+>+>?				 C

%'IJ

%'@A

>?3 -	 
! VS)vs#VS%(FC {"&#&;&#x0 ]  +*+~ 
!	 s   L %AL)L&%L&)L2c            
   #     K   t         j                  j                  } t         j                  j                  }t         j                  j	                  d      }t         j                  j	                  d      }t
        d d  t        z   D ]s  \  }}dD ]g  }dD ]^  }|dk(  r%|t        v rt        j                  |||||        .|dk(  r%|t        v rt        j                  |||||       X|dk(  r%|t        v rt        j                  |||||       |dk(  r%|t        v rt        j                  |||||       |dk(  r%|t        v rt        j                  |||||        |dk(  r&|t        v rt        j                  |||||       |dk(  r&|t        v rt        j                  |||||       ,|dk(  r&|t        v rt        j                  |||||       W||||f a j v y w)	NzTest fails and may be slow.reasonz*Test too slow to run to completion (>10m).)MLEMM)TFr   r~   r   )r   r   r   r   r   r   r   slow_fit_mler   xslow_fit_mlexfail_fit_mleskip_fit_mleslow_fit_mmxslow_fit_mmxfail_fit_mmskip_fit_mm)r   r   failr   r   r   methodfix_argss           r   cases_test_cont_basic_fitr      s    ;;DKKE;;#@AD;;#OPD!!'??##F)U?x<'? ,,xfhdSSU?x='@ ,,xfheTTU?x='@ ,,xfhdSSU?x<'? ,,xfhdSST>h+&= ,,xfhdSST>h,&> ,,xfheTTT>h,&> ,,xfhdSST>h+&= ,,xfhdSSVX555 * $ @s   G=G?c                     t        t        j                  t              j                  t                    t        t              t        t              z   t        t              z   k(  sJ t        t
        j                  t              j                  t                    t        t
              t        t              z   t        t              z   k(  sJ y N)lenr   unionr   r   r   r   r        r   test_cont_basic_fit_casesr      s    ##M288FG]!33c,6GGH I H""<066{CDL 11C4DDE F Er   zdistname, arg, method, fix_argsn_fit_samples   c                     	 t        t        |       }t        j                  j                  d      } |j                  |||d}|rt        ||||       y t        ||||       y # t        $ r | }Y _w xY w)Nr   r   )	r   r   r   r   r   r   r   check_fit_args_fixcheck_fit_args)r   r   r   r   r   r   r   r   s           r   test_cont_basic_fitr      sr    ) ))


'C
&**#-C
@C63V4vsC0  s   A$ $A21A2c                 0   	 t        t        |       }t        j                   |j
                  |       sJ t        j                   |j
                  |ddi      sJ t        j                   |j
                  |dd i      sJ y # t        $ r | }d} Y w xY w)Nrq   r   r   )r   r   r   r   isscalarr   )r   r   r   s      r   test_rvs_scalarr     s    +)
 ;;zvzz3'(((;;zvzz30R0111;;zvzz32T2333  +*+s   B BBc                  8    t        t        j                  d       y )N)r   皙?)r   r   r4   r   r   r   &test_levy_stable_random_state_propertyr     s       1 1:>r   c               #     K   t               } dh}dh}t        d d  t        z   D ]  \  }}|dk(  r|t        v r=t	        j
                  ||ddddt        j                  j                  d             Q|| v}||v}||v}t               }t	        j
                  |||||d|       |r|r|rt	        j
                  ||ddddt        j                  j                  g|z           y w)	Nr/   rZ   r4   Ttoo slowr   r~   F)
setr   r   xslow_test_momentsr   r   r   r   listxfail)	fail_normalizationfail_higherfail_momentr   r   cond1cond2cond3r   s	            r   cases_test_momentsr   "  s    'K-K!!'??#}$)),,xdD$%+[[%6%6j%6%IK K 22++ ll8S% / 	/ E ,,xdD$&,kk&7&7%85%@B B9 @s   B,C,/=C,z=distname,arg,normalization_ok,higher_ok,moment_ok,is_xfailingc           
         	 t        t        |       }t        j                         5 }|j                  t        d       |j                  t        d       |j                  t        d       |r|j                  t                |j                  |ddi\  }}	}
}t        j                  d      5  |rt        |||        |rJt        ||||        t        ||||	|
|        t        ||||	|        t        ||||	||        t        |||        d d d        t        ||||	|        |rt!        ||||	|        d d d        y # t        $ r | }d} Y 8w xY w# 1 sw Y   FxY w# 1 sw Y   y xY w)	Nrq   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   )r   r   normalization_ok	higher_ok	moment_okis_xfailingr   r   r   r   sks               r   test_momentsr  G  sT   +)
 
			 C

%N	P

%8	:

%5	7 JJ)*!V\\377
1a[[X&#FC:!&#q(;!&#q!QA aH=!&#q!QA!&#x8 ' 	Q84aH57 
!	 	  +*+" '& 
!	 s6   D: BE-AE	(E:E
EE	EE#zdist,shape_argsc                    | dv rt        j                  d       | dv rt        rt        j                  d       | dv }t        t        |       }t        j                  d      }t        j                  d      }|j                  }g }ddg}t        |      D ]P  }	|	d	z   fd
|	dz   z  z   }
|j                  ||	   t        j                  |
      z         |j                  d|	d	z          R |j                  ||g       t        || |||d       y )N)r!   r&   r   r   r   )r<   rW   rg   rk   	exponnormr.   rC   r4   rI   r"   rc   ro   semicircularri   rl   rr   )rs   r   rs   rt   r   r   d)r   r   r   r   r   r   zerosonesr   rangeappendinsertextendr   )dist
shape_args
shape_onlydistfunclocscalenargsallargsbshaper  shps              r   test_rvs_broadcastr$  q  s    22J""y+,  M MJ
 ud#H
((1+CGGFOEEGVF5\1uhq1u%z!}RWWS\12aQ  NNC<  $SIr   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                    t         j                  }t        j                  |j	                  | |      |j                  | |      |j                  | |      g      }t        j                  |||g      }t        j                  |||       t        j                  |d   |d   z   d       t        j                  |j                  ||      g| gd       y )Nrtolr   r   r'  :0yE>)r   r   r   arraysfr   r   r   assert_allcloseassert_equalisf)	r   nr-  r   r   r*  ksvalsexpecteds	            r   test_gh17775_regressionr5    s    $ 
B88RUU1a["&&A,q!=>DxxS#'HhT2T!WtAw&,A!48r   c                     t         j                  j                  d      } t        j                  j                  t        j                  d      d|       }t        j                  |      }t        j                  t        j                  |dk7        d       t        j                  j                  dt        j                  d      |       }t        j                  |      }t        j                  t        j                  |dk7        d       t        j                  j                  t        j                  d      t        j                  d      |       }t        j                  |      }t        j                  t        j                  |dk7        d       t        j                  j                  t        j                  dgdgg      t        j                  d      |       }t        j                  |j                               }t        j                  t        j                  |dk7        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       y )N{   r   r   )r  r  r   r   z0All the values are equal, but they shouldn't be!)rr   rs   rt   r   )rr   rr   )r   r   r   r   )r   r   r   r   rr   )rt   )r   r   r   r   normr   r  diffr   assert_r  r  r,  ravelassert_raises
ValueErrorrd   )r   r3  r  s      r   test_rvs_gh2069_regressionr>    s    ))


$C::>>bhhqk>ED
AKKqAv RS::>>arwwqz>DD
AKKqAv RS::>>bhhqk#>ND
AKKqAv RS::>>bhhaSz2"''!*'*  ,D


AKKqAv RS*ejjnn1v1v.>q6Aq6"A'*ekkoo|Q6J*ekkoo|\3*d$r   c                      t         j                  } t        j                  | j	                  ddt
        j                        d       t        j                  | j                  dgdgt
        j                  g      d       t        j                  | j	                  dt
        j                   d      d       t        j                  | j                  dgt
        j                   gdg      d       t        j                  | j                  dgdgt
        j                  g      d       t        j                  | j	                  ddt
        j                        d       t        j                  | j                  dgt
        j                   gdg      d       t        j                  | j	                  dt
        j                   d      d       t        j                  | j	                  dt
        j                   d      d       y )Nr   r   gA?g|N?)r   rA   r   assert_almost_equalr   r   inf_cdf)tns    r   test_nomodify_gh9900_regressionrE    s{    
B BFF1a0.0BGGQC!rvvh7.0 BFF2w2/1BGGRDBFF7)aS9/1 BGGQC!rvvh7.0BFF1a0.0 BGGRDBFF7)aS9/1BFF1rvvgq1BFF2w2/1r   c            	          t        j                  g d      } t        j                  g d      }t        | |      D cg c]&  \  }}t        j                  j                  d||      ( }}}t        j                  t        j                  j                  d| |      |       t        | |      D cg c]&  \  }}t        j                  j                  d||      ( }}}t        j                  t        j                  j                  d| |      |       t        | |      D cg c]&  \  }}t        j                  j                  |||      ( }}}t        j                  t        j                  j                  | | |      |       t        | |      D cg c]&  \  }}t        j                  j                  |||      ( }}}t        j                  t        j                  j                  || |      |       y c c}}w c c}}w c c}}w c c}}w )N)r   rr   rs   rt   r   ru   )rw      r       r   0   rv   r   )r   r,  zipr   r:   r   r   assert_array_almost_equal)ry   b_a_banss        r    test_broadcast_gh9990_regressionrP    s    	#$A
&'A9<Qq
Bvr252r*C
B!!%"2"2"6"6q!Q"?E9<Qq
Bvr252r*C
B!!%"2"2"6"6q!Q"?E:=a(
C(B5B+(C
C!!%"2"2"6"6q!Q"?E:=a(
C(B5B+(C
C!!%"2"2"6"6q!Q"?E C C D Ds   +G80+G>"+H+H
c            	          t         j                  j                  t        j                  g d      dt        j                  g d      z
  dz  t        j
                  t        j                  g d      d       y )N)      @       @r'  r   )g      @r   g      @rR  )ry   rL  r  r  )r   rA   r   r   r,  rB  r   r   r    test_broadcast_gh7933_regressionrT    sQ    	OO
!/**c
1
&&HH_%  r   c            	         t         j                  ddd   } t         j                  t        j                  d       t        j                  d      f   }t        | |      D cg c]4  \  }}t        j
                  j                  ||t         j                        6 }}}t        j
                  j                  | |t         j                        }t        j                  ||       y c c}}w )Nrr   y             @Y@2   3   )
r   r_r  rJ  r   rA   r   rB  r   rK  )r   ry   _xrM  r4  rO  s         r   test_gh2002_regressionr[  %  s     	b4iA
rwwr{lBGGBK'(ABEa)L)B##BBFF3)HL
//

aBFF
+C!!#x0 Ms   "9C$c                  |    d} t         j                  j                  dt        j                  | g| dz   gg             y )Ng(\@r   )r   r]   ppfr   r,  )cs    r   test_gh1320_regressionr_  /  s3    A	bhha#gY'789r   c                     t         j                  j                  d       g d} ddt        j                  d      z  dz  z
  }ddt        j                  d      z  dz  z   }t        t        t        j                        t        j                  #  | d      \  }}t        j                  ||d	
       t        j                  ||z   |d	
       y )Ni  )r   r   r   r   r   皙?rr   rs   r   r   r   rt   )decimal)r   r   seedsqrtsupertyper   uniformfitr   rA  )r   ry   rL  r  r  s        r   test_method_of_momentsrj  ;  s    IINN4Aa
l1nAa
l1nAtEMM*EMM>q>NJCCA.CIq!4r   c                     t        j                  |       rt        ||        t        j                  |      rt        ||       y y r   )r   isfinitecheck_sample_meancheck_sample_var)popmeanpopvarsamples      r   r   r   G  s3    	{{7&'*	{{6( r   c                 R    t        j                  | |      j                  }|dkD  sJ y )Nr   )r   ttest_1samppvalue)rq  ro  probs      r   rm  rm  N  s&    VW-44D$;;r   c                     t        j                  | fd d      }|j                  }|j                  |j                  }}||cxk  r|k  sJ  J y )Nc                 (    | j                  d|      S )Nr   )ddofaxis)var)r   ry  s     r   <lambda>z"check_sample_var.<locals>.<lambda>Z  s    140r   gףp=
?)confidence_level)r   	bootstrapconfidence_intervallowhigh)rq  rp  resconfr  r  s         r   rn  rn  T  sT     //		0C
 ""D$))C& D     r   c                     g d}t        j                   | j                   | j                  |g| g| |t        |dz          y )N)gMbP?r   g+?z - cdf-ppf roundtriprc  err_msg)r   rA  r   r]  DECIMAL)r   r   msgvaluess       r   r   r   b  sF     FJFJJzvzz&'?3'?F#F"GS2>34r   c                     t        j                   | j                   | j                  g dg| g| g dt        |dz          y )Nr   r   ?z - sf-isf roundtripr  )r   rA  r-  r0  r  r   r   r  s      r   r   r   i  sB    IFIIjfjj&G3&GN#N+Wc1G23r   c                     t        j                   | j                  ddgg| d | j                  ddgg| z
  t        |dz          y )Nr   r  r'  z - cdf-sf relationshipr  )r   rA  r   r-  r  r  s      r   r   r   o  sO    JFJJSz8C8)&))S#J"="==$+S4656r   c                     t        j                  ddg      }t        j                   | j                  |g|  | j
                  d|z
  g| t        |dz          y )Nr   r  r   z - ppf-isf relationshipr  )r   r,  r   rA  r0  r]  r  )r   r   r  ps       r   r   r   v  sW    
#sAJFJJq/3/AaC1F#1F$+S5667r   c                 B    | j                   dg| }d} | j                  |g| }|dk  s|dkD  r|dz   } | j                  |g| } | j                  ||z   g|  | j                  ||z
  g| z
  |z  dz  }|dz  }t        j                  ||t
        |       y )	Nr   gư>g-C6?     @r   rS  z - cdf-pdf relationshipr  )r]  r   r   r   rA  r  )r   r   r  medianepspdfvcdfdiffs          r   r   r   }  s    VZZ"c"F
C6::f#s#Dt #vzz&'3'vzz&3,--vzz&3,--./23367G $$CD'7CHr   c                    t        j                  g d      } | j                  |g| }|t        j                  |         } | j                  |g| } | j
                  |g| }||dk7  t        j                  |      z     }|t        j                  |         }|dz  }t        j                  t        j                  |      |d|       y )N)ra  r   皙?r   333333?ffffff?皙?r    - logpdf-log(pdf) relationshiprv   r  	r   r,  r]  rl  r   r   r   rA  logr   argsr  pointsr3  r   r   s          r   r   r     s    XX9:F6::f$t$DD!"D
&**T
!D
!CV]]4'$'F
sax2;;s++
,CBKK'(F,,CBFF3KCHr   c                    t        j                  ddg      } | j                  |g| }|t        j                  |         } | j                  |g| } | j
                  |g| }||dk7  t        j                  |      z     }|t        j                  |         }|dz  }t        j                  t        j                  |      |d|       y )Nr   r   r  rv   r  r  r  s          r   r   r     s    XXq!fF6::f$t$DD!"D
&**T
!D
!CV]]4'$'F
sax2;;s++
,CBKK'(F,,CBFF3KCHr   c                 v   t        j                  g d      } | j                  |g| }|t        j                  |         } | j                  |g| } | j
                  |g| }||dk7     }|t        j                  |         }|dz  }t        j                  t        j                  |      |d|       y )N)	        ra  r   r  r   r  r  r  r'  r   z - logsf-log(sf) relationshiprv   r  )	r   r,  r]  rl  r-  r   r   rA  r  )r   r  r  r  r3  r-  r   s          r   r   r     s    XXCDF6::f$t$DD!"D	4	$	BFLL%%E	B!GB"++e$%E**CBFF2Jq#Fr   c                 v   t        j                  g d      } | j                  |g| }|t        j                  |         } | j                  |g| } | j
                  |g| }||dk7     }|t        j                  |         }|dz  }t        j                  t        j                  |      |d|       y )N)	r   ra  r   r  r   r  r  r  r'  r   z - logcdf-log(cdf) relationshiprv   r  )	r   r,  r]  rl  r   r   r   rA  r  )r   r  r  r  r3  r   r   s          r   r   r     s    XXABF6::f$t$DD!"D
&**T
!D
!CV]]4'$'F
cQh-CBKK'(F,,CBFF3KCHr   c                     d}g |z  }|r&|D cg c]  }t        j                  |g|z         }} | j                  dg| } | j                  dg| }|dz  }t        j                  ||g|z  d|       y c c}w )Nr   r   z - ppf multiplerv   r  )r   r,  r]  r   rA  )r   r   r  num_repeatsr  _r  medianss           r   r   r     s    KD
589S!{*+S9VZZ"c"Ffjj$t$GCGfX%;QPST :s    A9c                 $   t        j                  || |d      \  }}||k  rot        j                  | | |d      \  }}t        j                  ||kD  dt	        |      z   dz   t	        |      z   dz   t	        |      z   dz   t	        |      z          y y )Ni  )r  NzD = z	; pval = z
; alpha = z
args = )r   kstestr   r:  str)r  r  rM   r   Dpvals         r   r   r     s     ll3448GAtu ,,tT=4D5L&3q6/K"?#d)"K #!#&u:#.0;#<>A$i#H 	I 	r   c                 d    t        j                   | j                  |  | j                  |        y r   )r   r/  
vecentropyr   )r   r  s     r   r   r     s+    &V&&-v/EFr   c                     t        j                  ddg      t        j                  ddg      }} | j                  |||d\  }}t        j                  ||z  |z   |       t        j                  ||z  |z  |       y )N      $@g      4@r  r  )r   r,  r   r   r.  )	r   r   r   r   r  r  r  mtvts	            r   r  r    so     4,'4,)?CV\\3Cu5FB%#r*%r*r   c                      | j                   t        j                  g d      g| }t        j                  t        j
                  t        j                  |             |dz          y )Nr  zppf private is nan)_ppfr   r,  r   r:  anyisnan)r   r   r  ppfss       r   r   r     sH    6;;rxx0737DKKBFF288D>**C2F,FGr   c                     d\  }} | j                   | } | j                   |||d}t        j                  t        j                  |      |z  |z   t        j                  |             y )N)r   rr   r  )supportr   rA  r   r,  )r   r  r  r  suppsupp_loc_scales         r   r   r     s]    JC6>>4 D#V^^Ts%@NBHHTN5036HH^46r   c                 
   t        j                  d      5  t        j                         5 }|j	                  t
        d       |j	                  t
        d       | j                  ||      }| j                  |d|      }d d d        d d d        t        j                  t              d	t        |      z   k(         t        j                  t              d	t        |      z   k(         y # 1 sw Y   lxY w# 1 sw Y   pxY w)
Nr  r  !The shape parameter of the erlangcategorymessagezfloating point number truncatedrb  powell)	optimizerr   rr   	r   r  r   r   r   r   ri  r:  r   )r   r   r   r   r   r3  vals2s          r   r   r     s    		"C$9$9$;s

N> 	 	@

N< 	 	>zz#fz-

3(6
B %<	" KKD	Qs3xZ'(KKE
aCj() %<$;	"	"s#   C9AC-C9-C6	2C99Dc           	         t        j                  d      5  t        j                         5 }|j	                  t
        d       | j                  |d|      }| j                  |d|      }t        j                  t        |      d	t        |      z   k(         t        j                  |d
   dk(         t        j                  |d   dk(         t        j                  t        |      d	t        |      z   k(         t        |      dkD  rb| j                  ||d   |      }t        j                  t        |      d	t        |      z   k(         t        j                  |d   |d   k(         t        |      dkD  rb| j                  ||d   |      }t        j                  t        |      d	t        |      z   k(         t        j                  |d   |d   k(         t        |      d	kD  rb| j                  ||d	   |      }	t        j                  t        |	      d	t        |      z   k(         t        j                  |	d	   |d	   k(         d d d        d d d        y # 1 sw Y   xY w# 1 sw Y   y xY w)Nr  r  r  r  r   )flocr   r   )fscaler   rr   rV  r@  )f0r   )f1r   )f2r   r  )
r   r   r   r   r   r3  r  vals3vals4vals5s
             r   r   r     s   		"C$9$9$;s

N> 	 	@ zz#Afz5

3q
8CI3s8+,DHM"E"IN#CJ!CH*,-s8a<JJss1vfJ=EKKE
aCj01KKaCF*+s8a<JJss1vfJ=EKKE
aCj01KKaCF*+s8a<JJss1vfJ=EKKE
aCj01KKaCF*++ %<	"	"$;$;	"	"s#   I)H II)I&	"I))I2c               #      K   t         D ]F  \  } }| t        v r3t        j                  | |t        j                  j
                         A| |f H y w)Nr~   )r   slow_with_listsr   r   r   r   r   s     r   cases_test_methods_with_listsr  $  sC     !#&,,xFKK4D4DEEC-	 "s   AAr   )r   r   r   r   r-  r   r]  r0  zdistname, argsc           	      N   t        t        |      }t        ||       }|dk(  r| j                  d      rddg}nddg}|D cg c]  }|gdz  
 }}ddg}dd	g}	 ||g|||	d
}
t        j                  |
t        |g|||	 D cg c]  } || 	 c}dd       y c c}w c c}w )Nr   r  r   rr   r   ra  r   r   g)\(?r  g+=r&  )r*  atol)r   r   
startswithr   r.  rJ  )r   r   r  r  rG   r   ry   shape2r  r  resultr   s               r   test_methods_with_listsr  ,  s     5(#DfA<F$5$5e$<!H#J!"Tqc!eTF"c(CIEq060s%0F(+A(C(C(CU(CD(C1A(CD"0	 #
 Es    BB"
c                      t         j                  j                  ddd      } t        | t              sJ t         j
                  j                  dd      }t        |t              sJ y )Nr   r   rt   rw   )r   rJ   moment
isinstancefloatrK   )vals0vals1s     r   (test_burr_fisk_moment_gh13234_regressionr  B  sP    JJaA&EeU###JJa#EeU###r   c            
         t         j                  j                  dt        j                  g d      d      } t        j                  g d      }t        j                  | |       t         j                  j                  dt        j                  g d      d      }t        j                  t        j                  t        j                  t        j                  g      }t        j                  ||       t         j                  j                  dt        j                  g d      g d      }t        j                  t        j                  dt        j                  g      }t        j                  ||       t         j                  j                  dd	d
      }t        j                  }t        j                  ||       t        ||j                        sJ t         j                  j                  dg dg d      }t        j                  dt        j                  t        j                  g      }	t        j                  ||	       t         j                  j                  dg dg d      }
t        j                  dt        j                  dg      }t        j                  |
|       t         j                  j                  dg dg d      }t        j                  t        j                  t        j                  t        j                  g      }t        j                  ||       t         j                  j                  ddd	d	      }t        j                  }t        j                  ||       t        ||j                        sJ t         j                  j                  dg dd	d	      }t        j                  t        j                  t        j                  t        j                  g      }t        j                  ||       t         j                  j                  dg dg dg d      }t        j                  dt        j                  t        j                  g      }t        j                  ||d       t         j                  j                  dddgddg      }t        j                  ddg      }t        j                  ||d       dd	ddt        j                  g}dd	ddt        j                  g}dd	t        j                  d      g}d d	ddt        j                  g}t        j                  |      j                  d!      }t        j                  |      j                  d"      }t        j                  |      j                  d#      }t        j                  |      }t         j                  j                  d||||$      }t        j                   ||||      \  }}}}t        j"                  |      D ]z  }t        j$                  d%d%&      5  |d	   }t         j                  j                  d||   ||   ||   ||   $      }t        j&                  j                  ||   |       d d d        | y # 1 sw Y   xY w)'Nr   r   rr   rs   orderr  r  r'  rS  rR  r@  )r   r   rS  rr   r   )r   r   rr   )rs   r  r  )r   r   r   )r   rV  g      Y@g      9@r  )g      rV  g      Yr  dfr  r  )r'  r  rS  )r'  g      r  g ﲙ"@r+  r)  r   r'  g      :@g      @gvIh%<=gg@g )r@  r   r   r   )r@  r   r   )r@  r   )r  ry   rL  r  r  r  )invaliddivide)r   r8  r  r   r,  r   r/  nanr  r   rb   r.  pire  reshaper;   broadcast_arraysndenumerater  testing)r  	expected0r  	expected1r  	expected2r  	expected3r  	expected4r  	expected5vals6	expected6vals7	expected7vals8	expected8vals9	expected9vals10
expected10ry   rL  r  r  vals11ir4  s                                r   *test_moments_with_array_gh12192_regressionr  J  s   JJA288I+>aHE&IUI& JJA288I+>bIE"&&"&&"&&12IUI& JJA288I+>$.  0E"&&"bff-.IUI& JJA1B7EIUI&eY00111 JJA9KHE#rvvrvv./IUI& JJA9NKE#rvvv./IUI& JJA9<LME"&&"&&"&&12IUI& II1;EIUI&eY00111 II1CE"&&"&&"&&12IUI& II1#0  2E*bffbff56Iyt4 ZZq2r(RH5F3+&J
7 
q!S"%% A	q!S"%% ABGGAJ
C1abee$E
M*A
J'A
((3-


(CHHUOEZZQ!qcGF**1ae<Aq#u^^A[[(;!Azz((qAaDAaD-0V58 ) EHJJ##F1Ix8 <; ;;s   AY::Z	c                     t         j                  j                  dt        j                  g d      dgg      } t        j                  g dg      }t        j                  | |       | j                  |j                  k(  sJ t         j                  j                  dt        j                  dgdgdgg      g d      }t        j                  g dg dg d	g      }t        j                  ||       |j                  |j                  k(  sJ t         j                  j                  dg dd
d      }t        j                  g d      }t        j                  ||d       |j                  |j                  k(  sJ t         j                  j                  ddgdgdggg dg d      }t        j                  t        j                  t        j                  dgt        j                  t        j                  dgt        j                  t        j                  dgg      }t        j                  ||d       |j                  |j                  k(  sJ y )Nr   r  r  r  rr   rs   )r'  r'  r'  )rS  rS  rS  )rR  rR  rR  r  r'  r  )g e3E?g> ?g e3E?r+  r)  rS  rR  )r  r'  rS  )g      r  rR  g@3F@gw,W.
@gf3&@)r   r8  r  r   r,  r   r/  shaperb   r.  r  )r  r  r  r  r  r  r  r  s           r   /test_broadcasting_in_moments_gh12192_regressionr    s   JJA288I+>seLE,(IUI&;;)//)))JJA288aS1#sO+D$-  /E,lCDIUI&;;)//)))II12RHE=>Iyt4;;)//)))II12$rd);#0  2E266266:666266:666266:68 9I yt4;;)//)))r   c                  ^   g d} d}t        j                  |D cg c]1  }| D cg c]#  }t        j                  j                  ||      % c}3 c}}      }t        j                  t        j                  j                  | |            }t	        j
                  ||       y c c}w c c}}w )N)r   r         @      @      @r  r  )r   r,  r   rQ   r   r.  )shapesr  r  r  r  res2s         r   test_kappa3_array_gh13582r    s     'FG
((6=?6=F( &u \\''v'> &(6=? @C88ELL&&vw&?@DT"( ?s   	B)
(B$B)
$B)
c                     t        j                  g d      } t        j                  g d      }d}t        j                  |D cg c]A  }t        d      D cg c]*  }t        j                  j                  | |   ||   |      , c}C c}}      }t        j                  t        j                  j                  | ||            }t        j                  ||       t        j                  g d      } t        j                  g d      }t        j                  |D cg c]A  }t        d      D cg c]*  }t        j                  j                  | |   ||   |      , c}C c}}      }t        j                  t        j                  j                  | ||            }t        j                  ||       t        j                  g d	      } t        j                  g d
      d d d f   }t        j                  t        j                  j                  | ||            }|j                  dk(  sJ y c c}w c c}}w c c}w c c}}w )N)r   r
  r  r  r  )r   r   g      r   r  r  r   r  )r@        пr  r   r@  r   )r   r   r   gUUUUUUտr@  r   ru   )r@  r   r   )r@  r   r   r   )rt   rt   rs   )r   r,  r  r   r   r   r.  r  )hr  r  r  r  r  r  s          r   test_kappa4_array_gh13582r    s   
*+A
()AG
((4;=4;&E!H&$q \\''!adF'C$&4;= >C88ELL&&q!W&=>DT" 	+,A
)*A
((4;=4;&E!H&$q \\''!adF'C$&4;= >C88ELL&&q!W&=>DT" 	A
!"1d7+A88ELL&&q!W&=>D::"""#& =& =s0   H:
/H5H:
 I
2/I !I
5H:
 I
c                     d} t        j                  t        |       5  t        j                         j
                   d d d        t        j                  t        |       5  t        j                         j                   d d d        dt        j                  _        t        j                         }t        |t              sJ t        t        j                  d       y # 1 sw Y   xY w# 1 sw Y   kxY w)Nz.'rv_continuous_frozen' object has no attribute)matchherringpmf)
r   r   AttributeErrorr   r8  r  logpmfr  r   delattr)r  frozen_norms     r   test_frozen_attributesr    s     ?G	~W	5

 
6	~W	5

 
6EJJN**,Kk#7888EJJ 
6	5	5	5s   C C)C&)C2c                     t         j                  j                  d      } t         j                   d}}t	        j
                  t        j                  j                  ||      t        j                  j                  |             | j                  d      | j                  d      }}| j                  d      dk  }d||<   ||   }t        j                  j                  ||      }t	        j
                  ||   t        j                  j                  |             t	        j
                  ||    t        j                  j                  ||    ||                 y )Nr   )rs   rs   )r   r   )
r   r   default_rngrB  r   r/  r   ro   r   r8  )r   r   ry   maskx_normr  s         r   test_skewnorm_pdf_gh16038r"    s   
))


"CFF7AqAU^^''1-uzz~~a/@A::6:"CJJFJ$;qA::6:"S(DAdGtWF
..

Q
"CSY

v 67S$Z!3!3AteHah!GHr   r   r   r   r   r   r-  r   r]  r0  r  entropyexpectr  meanstdrz  intervalgffffff?r  r   mvcasec                    | \  }}t        t        j                         |      } || }| t        v rt	        |t
        j                        sJ y t	        |d   t
        j                        sJ t	        |d   t
        j                        sJ y )Nr   r   )r   r   r8  
scalar_outr  r   number)r)  method_namer  r   r  s        r   test_scalar_for_scalarr.    su     KUZZ\;/F
$-Cz#ryy)))#a&")),,,#a&")),,,r   c                     t         j                  j                  g d      } t        | d   t        j
                        sJ t        | d   t        j
                        sJ t         j                  j                  g d      } t        | d   t        j
                        sJ t        | d   t        j
                        sJ t         j                  j                  dg d      } t        | t        j
                        sJ y )Nr  r   r   r   )r   r8  ri  r  r   r,  fit_loc_scalennlf)r  s    r   test_scalar_for_scalar2r2    s    
**..
#Cc!fbii(((c!fbii(((
**
"
"9
-Cc!fbii(((c!fbii(((
**//&)
,Cc299%%%r   r   )ysysnumpyr   numpy.testingr  r   r   r   r<  scipy.integrater   	itertoolsscipyr   common_testsr   r   r	   r
   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   scipy.stats._distr_paramsr   !scipy.stats._distn_infrastructurer   r  maxsizer   r   r   r   r   r   r   r   r   r   r   r   r   r  r   case1case2productr)  r|   	histogram_histrv_histogram_rv_histr  tupler   r   parametrizer   r   r   r   r   r   r   r   r  r$  r5  r>  rE  rP  rT  r[  r_  rj  r   rm  rn  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+  scalars_outr.  r2  r   r   r   <module>rG     s   
    * .  
? 
? 
? 
? 
? 
? / B	 [[5 	<  C : 9OK23R2K5? ''  *34	6!fj)&Y&&u~e}EMD'BLL1411E!u!!%9H##Xuw$78 F	  )>)@Au%R1 & BR1j"6JF :2463%01 161 )>)@A
4 B
4?"BJ  '+-/#6/ #6L *H5&J 6&Jj .?	@9	@9$$:!1HF,1:	5)!4367I$
I
I
G
I
UI&G+H6
*,2   $A B)+H+JK0 LB0&$Q9h*4 # # # #0
I bkEA3<(QC51#,!oqc{WqcNUQCLqclXsOi_xnnvrlUBK%F
 TF#i_w6GH k!9:
- ;
-	&r   