
    \h                        S SK Jr  S SKJr  S SKJr  S SKJr  S SKJ	r	  S SK
Jr  S SKJr  S SKJrJrJrJrJrJrJrJrJrJr  S S	KJrJrJrJrJrJrJ r   S S
K!J"r"  \	" SSS/0S9r#\S 5       r$\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S j5       r&\$RK                  \5      S\4S  j5       r&\$RK                  \5      S\4S! j5       r&\$RK                  \5      S\4S" j5       r&\$RK                  \ 5      S\ 4S# j5       r&\$RK                  \"5      S\"4S$ j5       r&g%)&    )singledispatch)Dummy)exp)lambdify)import_module)DiscreteDistributionHandmade)SingleContinuousDistribution)
ChiSquaredDistributionExponentialDistributionGammaDistributionLogNormalDistributionNormalDistributionParetoDistributionUniformDistributionBetaDistributionStudentTDistributionCauchyDistribution)GeometricDistributionLogarithmicDistributionNegativeBinomialDistributionPoissonDistributionSkellamDistributionYuleSimonDistributionZetaDistribution)SingleFiniteDistributionscipyfromliststats)import_kwargsc                     g N distsizeseeds      Y/var/www/auris/envauris/lib/python3.13/site-packages/sympy/stats/sampling/sample_scipy.pydo_sample_scipyr(      s        r$   c                 B  ^ SS K n[        S5      n[        X@R                  U5      SS/5      m " U4S jSUR                  R
                  5      nU" [        U R                  R                  5      [        U R                  R                  5      SS9nUR                  XS9$ )	Nr   znumpyr   c                   "   > \ rS rSrU 4S jrSrg)_.<locals>.scipy_pdf#   c                    > T" U5      $ r!   r"   )r$   xhandmade_pdfs     r'   _pdf_.<locals>.scipy_pdf._pdf$       ?"r)   r"   N)__name__
__module____qualname____firstlineno__r3   __static_attributes__)r2   s   r'   	scipy_pdfr.   #       	# 	#r)   r;   abnamer%   random_state)scipy.statsr   r   pdfr   rv_continuousfloatset_inf_suprvs)r$   r%   r&   r   r+   r;   scipy_rvr2   s          @r'   _rL      s     c
AAxx{Wg,>?L#EKK-- # 5/ /kCH<<T<55r)   c                 x    [         R                  R                  R                  [	        U R
                  5      XS9$ N)dfr%   rB   )r   r   chi2rJ   rF   kr#   s      r'   rL   rL   ,   s,     ;;5=tOOr)   c                 ~    [         R                  R                  R                  S[	        U R
                  5      -  XS9$ )N   )scaler%   rB   )r   r   exponrJ   rF   rater#   s      r'   rL   rL   2   s2     ;;  q5+;';$ ZZr)   c                     [         R                  R                  R                  [	        U R
                  5      [	        U R                  5      XS9$ )N)r>   rT   r%   rB   )r   r   gammarJ   rF   rQ   thetar#   s      r'   rL   rL   8   s8     ;;  5=djj8IPT hhr)   c                     [         R                  R                  R                  [	        [        U R                  5      5      [	        U R                  5      XS9$ )N)rT   sr%   rB   )r   r   lognormrJ   rF   r   meanstdr#   s      r'   rL   rL   >   s<     ;;""s499~)>%/X\"ppr)   c                     [         R                  R                  R                  [	        U R
                  5      [	        U R                  5      XS9$ N)locrT   r%   rB   )r   r   normrJ   rF   r]   r^   r#   s      r'   rL   rL   D   s6    ;;E$))$4E$((ORVjjr)   c                     [         R                  R                  R                  [	        U R
                  5      [	        U R                  5      XS9$ )N)r?   rT   r%   rB   )r   r   paretorJ   rF   alphaxmr#   s      r'   rL   rL   I   s8     ;;!!E$**$5U477^RV!jjr)   c                 x    [         R                  R                  R                  [	        U R
                  5      XS9$ rN   )r   r   trJ   rF   nur#   s      r'   rL   rL   O   s(    ;;==dggTMMr)   c                     [         R                  R                  R                  [	        U R
                  5      [	        U R                  U R
                  -
  5      XS9$ r`   )r   r   uniformrJ   rF   leftrightr#   s      r'   rL   rL   T   sE     ;;""uTYY'7uTZZRVR[R[E[?\cg"{{r)   c                     [         R                  R                  R                  [	        U R
                  5      [	        U R                  5      XS9$ )N)r>   r?   r%   rB   )r   r   betarJ   rF   re   r#   s      r'   rL   rL   Z   s8     ;;%

"3uTYY7Gdffr)   c                     [         R                  R                  R                  [	        U R
                  5      [	        U R                  5      XS9$ r`   )r   r   cauchyrJ   rF   x0rX   r#   s      r'   rL   rL   `   s6    ;;!!eDGGnE$**<MTX!llr)   c                   ^ SSK Jn  [        S5      n[        X@R	                  U5      SS/5      m " U4S jSU5      nU" [        U R                  R                  5      [        U R                  R                  5      SS9nUR                  XS	9$ )
Nr   rv_discreter+   r,   r   c                   "   > \ rS rSrU 4S jrSrg)_.<locals>.scipy_pmfn   c                    > T" U5      $ r!   r"   )r$   r1   handmade_pmfs     r'   _pmf_.<locals>.scipy_pmf._pmfo   r5   r)   r"   N)r6   r7   r8   r9   r{   r:   )rz   s   r'   	scipy_pmfrw   n   r<   r)   r}   r=   rA   )
rC   ru   r   r   rD   rF   rG   rH   rI   rJ   )r$   r%   r&   ru   r+   r}   rK   rz   s          @r'   rL   rL   g   sr    'c
AAxx{Wg,>?L#K # 5/53G)+H<<T<55r)   c                 x    [         R                  R                  R                  [	        U R
                  5      XS9$ N)pr%   rB   )r   r   geomrJ   rF   r   r#   s      r'   rL   rL   w   *    ;;%-dNNr)   c                 x    [         R                  R                  R                  [	        U R
                  5      XS9$ r   )r   r   logserrJ   rF   r   r#   s      r'   rL   rL   |   s*    ;;!!E$&&M!PPr)   c                     [         R                  R                  R                  [	        U R
                  5      [	        U R                  5      XS9$ )N)nr   r%   rB   )r   r   nbinomrJ   rF   rr   r#   s      r'   rL   rL      s4    ;;!!E$&&MU466]!aar)   c                 x    [         R                  R                  R                  [	        U R
                  5      XS9$ )N)mur%   rB   )r   r   poissonrJ   rF   lamdar#   s      r'   rL   rL      s+    ;;""eDJJ&7d"VVr)   c                     [         R                  R                  R                  [	        U R
                  5      [	        U R                  5      XS9$ )N)mu1mu2r%   rB   )r   r   skellamrJ   rF   r   r   r#   s      r'   rL   rL      s5    ;;""uTXXE$((ORV"jjr)   c                 x    [         R                  R                  R                  [	        U R
                  5      XS9$ )N)re   r%   rB   )r   r   	yulesimonrJ   rF   rhor#   s      r'   rL   rL      s*    ;;  $$5?$YYr)   c                 x    [         R                  R                  R                  [	        U R
                  5      XS9$ )N)r>   r%   rB   )r   r   zipfrJ   rF   r[   r#   s      r'   rL   rL      r   r)   c                     SSK Jn  U R                  n/ / peUR                  5        H9  u  pxUR	                  [        U5      5        UR	                  [        U5      5        M;     U" SXV4S9n	U	R                  XS9$ )Nr   rt   rK   )r@   valuesrA   )rC   ru   dictitemsappendintrF   rJ   )
r$   r%   r&   ru   density_r1   yrQ   vrK   s
             r'   rL   rL      sk     (yyHrq 	Q	q ! 
A6:H<<T<55r)   N)'	functoolsr   sympy.core.symbolr   &sympy.functions.elementary.exponentialr   sympy.utilities.lambdifyr   sympy.externalr   sympy.statsr   sympy.stats.crvr	   sympy.stats.crv_typesr
   r   r   r   r   r   r   r   r   r   sympy.stats.drv_typesr   r   r   r   r   r   r   sympy.stats.frvr   r   r(   registerrL   r"   r)   r'   <module>r      s   $ # 6 - ( 4 8- - -V V V 4 	gj'-CD   676( 6 86  01P" P 2P
 12[# [ 3[
 +,i i -i
 /0q! q 1q
 ,-k k .k ,-k k .k
 ./N  N 0N -.| | /|
 *+g g ,g
 ,-m m .m 676( 6 86 /0O! O 1O 12Q# Q 3Q 67b( b 8b -.W W /W -.k k /k /0Z! Z 1Z *+O O ,O 23
6$ 
6 4
6r)   