+
    0iOY              	       x   ^ RI HtHtHt ^ RIt^ RIHtHtH	t	H
t
HtHt ^ RIt^ RIHt R tR t]P$                  P'                  R^]P(                  ! R]P$                  P*                  .R	7      .4      R
 4       t]P$                  P'                  R^]P(                  ! R]P$                  P*                  .R	7      .4      R 4       tR tR t ! R R]P4                  4      t ! R R]P4                  4      t ! R R]P4                  4      tR tR tR t R t!. R)Ot"]P$                  P'                  R]"RR.,           4      ]P$                  P'                  R]"4      R 4       4       t#R t$R t%R t&R  t']P$                  P*                  R! 4       t(R" t)]P$                  P*                  R# 4       t*R$ t+R% t,R& t-R' t.R( t/R# )*    )statslinalg	integrateN)assert_almost_equalassert_assert_equalassert_array_almost_equalassert_array_almost_equal_nulpassert_allclose)raisesc            	         \         P                  P                  R 4      p RpV P                  ^ ^V4      pVP	                  4       pVP                  ^R7      p\        P                  ! V4      p\         P                  ! . RO4      pVP                  \         P                  ! VP                  4      r\        V! V4      \        P                  P                  VR,          WxR7      P                  RR7      VP                   ,          RR7       \         P"                  ! R^R	4      p	VP%                  V	4      p
\        P                  P                  WVR7      pV	^,          V	^ ,          ,
          p\'        \         P                  ! W,
          ^,          4      V,          R
8  4       VP)                  V\         P*                  4      pVP)                  \         P*                  ) V4      p\-        VR^R7       \-        VR^R7       \-        VP/                  V\         P*                  4      V^R7       \-        VP/                  \         P*                  ) V4      V^R7       \-        VP1                  V4      V
^,          P                  4       V,          ^R7       \-        VP3                  W4^,          4      W,          P                  4       V,          ^R7       R# )   )ddof      ?NlocscaleaxisvIh%,=rtol  {Gz?decimalg?r   g?NNNN)nprandomdefault_rngnormalmeanstdr   gaussian_kdeasarraydatasetsqrt
covariancer   normpdfsumnlinspaceevaluater   integrate_box_1dinfr   integrate_boxintegrate_kdeintegrate_gaussian)rngn_basesamplexnxnmeanxnstdgkdexxr   r   xskdepdfnormpdf	intervallprob1prob2s                  [/var/www/html/photoedit/myenv/lib/python3.14/site-packages/scipy/stats/tests/test_kdeoth.pytest_kde_1drG   
   s
   
))


(CL	Aq,	'BWWYFFFFNE b!D 
O	$Brwwt7R

r'{9==2=FO 
RC	 B]]2FjjnnR5n9G11IBFFF$q())3d:;!!&"&&1E!!266'62EsA.sA.**6266:E2N**BFF7F;UBO**40))3Q@//qA,,.y8!E    c                  6   \         P                  P                  R 4      p RpV P                  ^ ^V4      pV P                  V4      p\         P                  ! W#R7      p\         P
                  ! \         P                  ! W$,
          ^,          VR7      4      p\        P                  ! W#R7      p\         P                  ! . RO4      pVP                  \         P
                  ! VP                  4      r\        P                  P                  p
\        V! V4      \         P                  ! V
! VR,          WR7      VP                  ,          RR7      RR7       \         P                   ! R^R	4      pVP#                  V4      p\        P                  P                  WVR7      pV^,          V^ ,          ,
          p\%        \         P                  ! W,
          ^,          4      V,          R
8  4       VP'                  V\         P(                  4      pVP'                  \         P(                  ) V4      p\+        VR^R7       \+        VR^R7       \+        VP-                  V\         P(                  4      V^R7       \+        VP-                  \         P(                  ) V4      V^R7       \+        VP/                  V4      V^,          P                  4       V,          ^R7       \+        VP1                  WE^,          4      W,          P                  4       V,          ^R7       R# )r   r   weightsr   Nr   r   r   r   r   r   r   r   r   r!   r"   )r#   r$   r%   r&   averager,   r   r)   r*   r+   r-   r.   r/   r   r0   rK   r2   r3   r   r4   r5   r   r6   r7   r8   )r9   r:   r;   wnr<   r=   r>   r?   r   r   r/   r@   rA   rB   rC   rD   rE   s                    rF   test_kde_1d_weightedrN   1   s+   
))


(CL	Aq,	'B	L	!BZZ'FGGBJJ	A~r:;E b-D 
O	$Brwwt7
**..CR
s2g;C5D2N 
RC	 B]]2FjjnnR5n9G11IBFFF$q())3d:;!!&"&&1E!!266'62EsA.sA.**6266:E2N**BFF7F;UBO**40))3Q@//qA,,.y8!ErH   r:   r   )marksc           
         \         P                  P                  R 4      p\         P                  ! RR.4      p\         P                  ! RR.RR..4      pVP	                  W#V R7      P
                  p\        P                  ! V4      p\         P                  ! ^^.^^.^^..4      pVR,          VP                  P
                  ,
          p\        P                  P                  p\        V! VP
                  4      V! WuP                  R7      P                  RR7      VP                  ,          R	R
7       \        P                  P                  p	RR.^ ^ .rWP                  P
                  ,
          WP                  P
                  ,
          r\        VP!                  WVR7      V	! WVP                  VR7      P                  RR7      VP                  ,          RR
7       \         P"                  R^R1R^R13,          w  r\         P$                  ! VP'                  4       VP'                  4       .4      pVP)                  V4      pVP+                  RR4      p\        P                  P                  \         P,                  ! W.4      W#R7      pVP'                  4       ^,          VP'                  4       ^ ,          ,
          p\/        \         P                  ! VV,
          ^,          4      V^,          ,          R8  4       RpRpVP!                  VV^,          .VV.VR7      pVP!                  VV.VV^,          .VR7      p\1        VR^R7       \1        VR^R7       \1        VP3                  V4      V^,          P                  4       V^,          ,          ^R7       \1        VP5                  W#4      VV,          P                  4       V^,          ,          ^R7       R# )r         ?      @       @      @sizeNcovr   r   r   r9   lower_limitrX   r9   gƠ>             @@r   r'   rX   r   }Ô%ITr   r   r    Nr    r!   r"   }Ô%I)r#   r$   r%   arraymultivariate_normalTr   r)   r*   r+   r/   r   r-   r0   r1   cdfr6   mgridvstackravelr3   reshapedstackr   r   r7   r8   )r:   r9   r'   r-   r;   r>   r?   argr/   re   lohilo_hi_xygrid_coordsrA   rB   rC   smalllargerD   rE   s   &                       rF   test_kde_2dru   \   s    ))


(C88S#JDC:Sz23J 
	 	 	 	E	G	GB b!D 
aVaVaV,	-B
Z.4<<>>
)C

#
#
'
'CRTT
C__%))r)2TVV; 
#
#
'
'C"X1vLLNN"B$72s+Cdoo3?CCCLtvvU 88BqIr!Dy()DA))QWWY	23K]];'F^^C%F''++BIIqf,=15 , GG	!qwwy|+IBFFFW$q()Y\:TABEEtAw/%SIEu~tAw/?SIEsA.sA.**40)9a<8!E//A,,.	1=qJrH   c                    \         P                  P                  R 4      p\         P                  ! RR.4      p\         P                  ! RR.RR..4      pVP	                  W#V R7      P
                  pVP                  V 4      p\        P                  ! WER7      p\         P                  ! ^^.^^.^^..4      pVR,          VP                  P
                  ,
          p\        P                  P                  p	\        V! VP
                  4      \         P                  ! V	! WP                  R7      VP                  ,          RR	7      R
R7       \        P                  P                   p
RR.^ ^ .rWP                  P
                  ,
          WP                  P
                  ,
          r\        VP#                  WVR7      \         P                  ! V
! WVP                  VR7      VP                  ,          RR	7      RR7       \         P$                  R^R1R^R13,          w  pp\         P&                  ! VP)                  4       VP)                  4       .4      pVP+                  V4      pVP-                  RR4      p\        P                  P                  \         P.                  ! VV.4      W#R7      pVP)                  4       ^,          VP)                  4       ^ ,          ,
          p\1        \         P                  ! VV,
          ^,          4      V^,          ,          R8  4       RpRpVP#                  VV^,          .VV.VR7      pVP#                  VV.VV^,          .VR7      p\3        VR^R7       \3        VR^R7       \3        VP5                  V4      V^,          P                  4       V^,          ,          ^R7       \3        VP7                  W#4      VV,          P                  4       V^,          ,          ^R7       R# )r   rQ   rR   rS   rT   rU   rJ   NrW   r   r   r   rY   rZ   gh㈵>r\   r   r]   r   r^   r   r   r_   r!   r`   r"   ra   )r#   r$   RandomStaterb   rc   rd   randr   r)   r*   r+   r/   r   r0   r-   rK   re   r6   rf   rg   rh   r3   ri   rj   r   r   r7   r8   )r:   r9   r'   r-   r;   rM   r>   r?   rk   r/   re   rl   rm   rn   ro   rp   rq   rr   rA   rB   rC   rs   rt   rD   rE   s   &                        rF   test_kde_2d_weightedry      s*    ))


(C88S#JDC:Sz23J 
	 	 	 	E	G	GB	,	B b-D 
aVaVaV,	-B
Z.4<<>>
)C

#
#
'
'CRTT

s3OO,t||;"E 
#
#
'
'C"X1vLLNN"B$72s+
s3T__#F||"$	&	 88BqIr!Dy()DAq))QWWY	23K]];'F^^C%F''++BIIq!f,=15 , GG	!qwwy|+IBFFFW$q()Y\:TABEEtAw/%SIEu~tAw/?SIEsA.sA.**40)9a<8!E//A,,.	1=qJrH   c                     R  p \         P                  P                  R4      p^2pVP                  ^ ^V4      p\        P
                  ! V4      p\        P
                  ! W0R7      p\        P
                  ! W4P                  R7      p\         P                  ! R^^34      pVP                  V4      pVP                  V4      p	\        W4       VP                  V4      p
\        W4       \        \        \        P
                  VRR7       R# )c                t    \         P                  ! V P                  RV P                  ^,           ,          4      # z*Same as default, just check that it works.g      )r#   powerr1   dkde_objs   &rF   scotts_factor0test_kde_bandwidth_method.<locals>.scotts_factor   s$    xx		3		!#455rH   r   	bw_methodwrongstringNr"   r#   r$   r%   r&   r   r)   factorr2   r3   r   assert_raises
ValueErrorr   r9   r:   r;   r>   gkde2gkde3r@   rA   kdepdf2kdepdf3s              rF   test_kde_bandwidth_methodr      s    6 ))


(CL	Aq,	'B b!Dr;Er[[9E	R"	B]]2FnnR G(nnR G(*e00"NrH   c                     R  p \         P                  P                  R4      p^2pVP                  ^ ^V4      p\        P
                  ! V4      p\        P
                  ! W0R7      p\        P
                  ! W4P                  R7      p\         P                  ! R^^34      pVP                  V4      pVP                  V4      p	\        W4       VP                  V4      p
\        W4       \        \        \        P
                  VRR7       R# )c                t    \         P                  ! V P                  RV P                  ^,           ,          4      # r|   )r#   r}   neffr~   r   s   &rF   r   9test_kde_bandwidth_method_weighted.<locals>.scotts_factor   s$    xxc799Q;&788rH   r   r   r   Nr"   r   r   s              rF   "test_kde_bandwidth_method_weightedr      s    9 ))


(CL	Aq,	'B b!Dr;Er[[9E	R"	B]]2FnnR G(nnR G(*e00"NrH   c                   &   a  ] tR tRt o R tRtV tR# )_kde_subclass1i  c                    \         P                  ! V4      V n        V P                  P                  w  V n        V n        V P                  V n        V P                  4        R # N)	r#   
atleast_2dr+   shaper~   r1   r   covariance_factor_compute_covariance)selfr+   s   &&rF   __init___kde_subclass1.__init__  sC    }}W-++!%!3!3  "rH   )r   r~   r+   r1   N)__name__
__module____qualname____firstlineno__r   __static_attributes____classdictcell____classdict__s   @rF   r   r     s     # #rH   r   c                   2   a a ] tR tRt oV 3R ltRtVtV ;t# )_kde_subclass2i  c                H   < V P                   V n        \        SV `  V4       R # r   )r   r   superr   )r   r+   	__class__s   &&rF   r   _kde_subclass2.__init__   s    !%!3!3!rH   )r   )r   r   r   r   r   r   r   __classcell__)r   r   s   @@rF   r   r     s     " "rH   r   c                   &   a  ] tR tRt o R tRtV tR# )_kde_subclass4i%  c                0    R V P                  4       ,          # )r   )silverman_factor)r   s   &rF   r    _kde_subclass4.covariance_factor&  s    T**,,,rH    N)r   r   r   r   r   r   r   r   s   @rF   r   r   %  s     - -rH   r   c                    a \         P                  ! . RO\        R7      p \         P                  ! R^
^2R7      p\        P
                  ! V 4      oS! V4      p\        V 4      pV! V4      p\        W$^
R7       \        V 4      pV! V4      p\        W&^
R7       \        V 4      pV! V 4      p. R	Op	\        W^R7       Sp
V3R lV
n        V
P                  4        V
! V4      p\        W+^
R7       R# )
   dtypenumnulpr   c                     < S P                   # r   )r   )kdes   rF   <lambda>/test_gaussian_kde_subclassing.<locals>.<lambda>H  s	    SZZrH   Nr"            )gn,?gǻe??gW[3?ga )?gKd8;?)r#   rb   floatr2   r   r)   r   r
   r   r   r	   r   r   )x1r@   yskde1y1kde2y2kde4y4
y_expectedkde5y5r   s               @rF   test_gaussian_kde_subclassingr   *  s    	#5	1B	S""	%B 

R
 C	RB "D	bB"23 "D	bB"23 "D	bBMJja8 D/D	bB"23rH   c                     \         P                  ! . RO\        R7      p \         P                  ! R	^
^R7      p. R
Op\        P
                  ! V 4      pVP                  RR7       VP                  RR7       V! V4      p\        W$^R7       R# )r   r   r   r   r   scottr   Nr   r   )g>L\9?gp4=?gH꟫?g8(H	MT?gĲCR?)r#   rb   r   r2   r   r)   set_bandwidthr	   )r   r@   r   r   r   s        rF   $test_gaussian_kde_covariance_cachingr   N  st    	#5	1B	S"!	$B NJ 

R
 C$(	RBja8rH   c                 L   \         P                  ! . RO\        R7      p \         P                  ! R^
^2R7      p\        P
                  ! V 4      pVP                  Vn        VP                  4        V! V4      p\        P
                  ! V RR7      pV! V4      p\        W5^
R7       R# )	zUgly, but people may rely on this.  See scipy pull request 123,
specifically the linked ML thread "Width of the Gaussian in stats.kde".
If it is necessary to break this later on, that is to be discussed on ML.
r   r   	silvermanr   r   Nr   r   )
r#   rb   r   r2   r   r)   r   r   r   r
   )r   r@   r   r   r   r   s         rF   test_gaussian_kde_monkeypatchr   ^  s    
 
#5	1B	S""	%B 

R
 C00C	RB bK8D	bB"23rH   c                     \         P                  ! ^4      p \        P                  ! V 4      p. ROp\	        V! V 4      V^R7       R# )zRegression test for #1181.r   N)B\A?筱ES?g
f?r   r   )r#   aranger   r)   r	   )r   r   r   s      rF   test_kde_integer_inputr   s  s4    	1B


R
 CMJc"gz1=rH   bw_typer   r   r   c                    \        \        V R 4      p VR9   d   TpM#\        \        VR 4      pV'       d	   V! ^4      MR p\        ;QJ d    R W3 4       F  '       g   K   RM	  RM! R W3 4       4      '       d   \        P                  ! 4        \        P
                  ! ^V R7      p\        P
                  ! ^V R7      p\        P                  ! WBVR7      p\        P
                  ! ^V R7      pV! V4      pVP                  \        P                  ! WF\        P                  ! V4      VP                  4      8X  g   Q hR # )Nc              3   (   "   T F  qR J x  K
  	  R # 5ir   r   ).0dts   & rF   	<genexpr>(test_kde_output_dtype.<locals>.<genexpr>  s     
,":s   TFr   )r   rK   )r   r   )getattrr#   anypytestskipr   r   r)   r   result_typefloat64r   )r   r   bwrK   r+   kpointsresults   &&      rF   test_kde_output_dtyper   ~  s     Bt$E(("gt,"WQZ
s
,
,sss
,
,,,ii'Gii'G7'BAYYq&FvYF<<2>>'2::g;N*+((4 4 4 4rH   c                  \   \         P                  P                  R 4      p V P                  R4      p\        P
                  ! V4      pV P                  R4      pRp\        P                  ! \        VR7      ;_uu_ 4        VP                  V4       RRR4       R#   + '       g   i     R# ; i)l	   Fgu	\U.lt&0 z0points have dimension 3, dataset has dimension 2matchN)   
   )   r   )
r#   r$   r%   standard_normalr   r)   r   r   r   logpdf)r9   r;   r>   r@   msgs        rF   test_pdf_logpdf_validationr     sw    
))

 F
GC			W	%Bb!D			W	%B
<C	z	-	-B 
.	-	-	-s   >BB+	c                  B   \         P                  P                  ^4      p ^2pV P                  ^ ^V4      p\        P
                  ! V4      p\         P                  ! R^^4      pVP                  V4      pVP                  V4      p\        WV^R7       \         P                  ! V4      pVP                  V4      p\        Wx^R7       \        P
                  ! V4      p\         P                  ! VP                  V4      4      pVP                  V4      p\        WV^R7       R# )r   r   N)r#   r$   r%   r&   r   r)   r2   r3   r/   r   logr   )	r9   r:   r;   r>   r@   r/   pdf2r   logpdf2s	            rF   test_pdf_logpdfr    s    
))


"CL	Aq,	'B b!D	S"b	!B
--
C88B<D2.VVC[Fkk"oG4 b!D
&&r"
#C;;r?D2.rH   c                     \         P                  P                  ^4      p ^2pV P                  ^ ^V4      pV P                  V4      p\        P
                  ! W#R7      p\         P                  ! R^^4      pVP                  V4      pVP                  V4      p\        Wg^R7       \         P                  ! V4      pVP                  V4      p	\        W^R7       \         P                  P                  R4      p \        P
                  ! WPP                  \        V4      4      R7      p\         P                  ! VP                  V4      4      pVP                  V4      p\        Wg^R7       R# )r   rJ   r   l   qX? Nr   )r#   r$   r%   r&   r   r)   r2   r3   r/   r   r   r   len)
r9   r:   r;   rM   r>   r@   r/   r   r   r   s
             rF   test_pdf_logpdf_weightedr    s	   
))


"CL	Aq,	'B	L	!B b-D	S"b	!B
--
C88B<D2.VVC[Fkk"oG4 ))



+Cb**SW*=>D
&&r"
#C;;r?D2.rH   c                    a	a
a \         P                  P                  R 4      p ^2p^
pV P                  W!3R7      pV P                  V^3R7      p\         P                  ! . RO4      o	\
        P                  ! V4      o
S
P                  S	4      pVP                  VS	,          4      pV	V
3R loV3R lpV! V4      p\        WhRR7       R# )   K b1d& rU   c                    <a  VVV 3R  lp\         P                  ! V\        P                  ) \        P                  4      ^ ,          # )c                 t   < \         P                  ! V .SS,          34      p SP                  V 4      ^ ,          # r   )r#   concatenater/   )rp   
dimensionsr   points   &rF   f<test_marginal_1_axis.<locals>.marginal_pdf_single.<locals>.f  s/    U:%678A771:a= rH   )r   quadr#   r5   r  r  r  r   s   f rF   marginal_pdf_single1test_marginal_1_axis.<locals>.marginal_pdf_single  s,    	! ~~a"&&"&&1!44rH   c                 6   < \         P                  ! S^ V R7      # )r   r   arrr#   apply_along_axis)r   r  s   &rF   marginal_pdf*test_marginal_1_axis.<locals>.marginal_pdf  s    ""#6QFKKrH   ư>r   N)	r   r   r   r   r      r      	   
r#   r$   r%   r&   rb   r   r)   marginalr/   r   )r9   n_datan_dimr+   r   r  r/   r  refr  r   r  s            @@@rF   test_marginal_1_axisr#    s    
))

 3
4CFEjjuoj.GZZeQZZ(F56J


W
%C||J'H
,,vj)
*C5L v
CC4(rH   c                    a	a
 \         P                  P                  R 4      p ^p^pV P                  W!3R7      pV P                  V^3R7      p\         P                  ! ^^.4      o	\
        P                  ! V4      o
S
P                  S	4      pVP                  VS	,          4      pV	V
3R lpV! V4      p\        WhRR7       R# )r  rU   c                 D   < VV3R  lp\         P                  ! V^ V R7      # )c                    <a  VVV 3R  lp\         P                  ! V\        P                  ) \        P                  \        P                  ) \        P                  4      ^ ,          # )c                 z   < SS,          w  r#\         P                  ! WW.4      pSP                  V4      ^ ,          # r	  )r#   rb   r/   )rq   rp   wzr  r   r  s   &&  rF   r  Rtest_marginal_2_axis.<locals>.marginal_pdf.<locals>.marginal_pdf_single.<locals>.f  s3    Z(HHaA\*wwqz!}$rH   )r   dblquadr#   r5   r  s   f rF   r  Gtest_marginal_2_axis.<locals>.marginal_pdf.<locals>.marginal_pdf_single  s<    % $$Q"&&"&&I!LLrH   r  r  )r   r  r  r   s   & rF   r  *test_marginal_2_axis.<locals>.marginal_pdf  s!    	M ""#6QFKKrH   r  r   Nr  )r9   r   r!  r+   r   r  r/   r  r"  r  r   s            @@rF   test_marginal_2_axisr.    s    
))

 3
4CFEjjuoj.GZZeQZZ(F1a&!J


W
%C||J'H
,,vj)
*CL v
CC4(rH   c                     \         P                  P                  R 4      p ^p^pV P                  W!3R7      pV P                  V^3R7      p\        P
                  ! V4      pR^.pVP                  V4      pVP                  WF,          4      p^R	.p	VP                  V	4      p
V
P                  WI,          4      p\        W4       Rp\        P                  ! \        VR7      ;_uu_ 4        VP                  ^R.4       RRR4       Rp\        P                  ! \        VR7      ;_uu_ 4        VP                  . R
O4       RRR4       Rp\        P                  ! \        VR7      ;_uu_ 4        VP                  . RO4       RRR4       R#   + '       g   i     L; i  + '       g   i     Lg; i  + '       g   i     R# ; i)r  rU   z,Elements of `dimensions` must be integers...r   g      @Nz,All elements of `dimensions` must be unique.z;Dimensions \[-5  6\] are invalid for a distribution in 4...r!   )r   r   r   )r   r   r  )r#   r$   r%   r&   r   r)   r  r/   r   r   r   r   )r9   r   r!  r+   r   r   dimensions1	marginal1pdf1dimensions2	marginal2r   messages                rF   test_marginal_ivr7    sR   
))

 3
4CFEjjuoj.GZZeQZZ(F


W
%C q'K[)I==,-Db'K[)I==,-D =G	z	1	1aX 
2 =G	z	1	1Y 
2 NG	z	1	1Z  
2	1 
2	1
 
2	1
 
2	1	1s$   *F,F.F2F	F/	2G	c            
      H   \         P                  P                  ^4      p RpRp\         P                  ! \	        ^ V4       Uu. uF  q0P                  ^ ^V4      V,           NK  	  up4      p\        P                  ! V4      pVP                  \         P                  ! ^ V4      4      p\         P                  P                  \         P                  ! V^ ,          4      R4       \         P                  P                  \         P                  ! V^ ,          4      R4       R# u upi )r   i	  i  FN)r#   r$   r%   rb   ranger&   r   r)   r   r   testingr   isneginfisnan)r9   n_dimensions	n_samplesr1   r;   r>   r   s          rF   test_logpdf_overflowr?  1  s     ))


"CLI	e	<?  ?::aI.!44 ?  
B b!D[[1l34FJJBKKq	2E:JJBHHVAY/7s   "Dc                  :   \         P                  P                  R 4      p V P                  ^dR7      pV P	                  . ROVP
                  R7      pVP                  4       p\        P                  ! \         P                  ! V4      VR7       \        W#RRR7       R# )i90  rU   rJ   +=atolr   N)rQ   g      $@d   )r#   r$   r%   	lognormalchoicerV   copyr   r)   log10r   )r9   valsrK   orig_weightss       rF   test_weights_intactrK  C  sl    
))


&C==c="Djj)		j:G<<>L	rxx~w7GEBrH   c                      . ROp . ROp\         P                  ! WR7      p\         P                  ! V \        P                  ! V4      R7      p. ROp\	        VP                  V4      VP                  V4      RRR7       R# )皙?rJ   rA  rB  N)rM  g      +@g      5@g     R@g     X@)r   r   r   r     )g333333?   X   )r   r)   r#   r   r   r3   )valuesrK   pdf_ipdf_fr;   s        rF   test_weights_integerrT  N  s`    *FGv7Evrzz'/BCE	BENN2&NN2&U@rH   c                  ,   R  p \         P                  P                  R4      pRpVP                  V4      pVP                  ^ ^V4      p\        P
                  ! V4      pV ! V4       \        P
                  ! WCR7      pV ! V4       \         P                  ! RR.4      p\         P                  ! RR.RR..4      pVP                  WxVR7      P                  p	\        P
                  ! V	4      p
V ! V
4       \        P
                  ! WR7      pV ! V4       R	# )
c                 j   ^pV P                  V4      pV P                  V4      p\        \        \        W#RR7       RpV P                  WR7      pV P                  WR7      p\        W#RR7       \        P
                  P                  ^R7      pV P                  WR7      p\        P
                  P                  ^R7      pV P                  WR7      p\        W#RR7       \        \        P
                  R4      '       d4   \        P
                  P                  R4      pV P                  WR7       R# R# )   gvIh%<=)rC  i?  )seedr%   i  N)	resampler   AssertionErrorr   r#   r$   rw   hasattrr%   )
gkde_trailn_samplesamp1samp2rX  rstate1rstate2r9   s   &       rF   test_seed_sub test_seed.<locals>.test_seed_sub\  s
   ##H-##H-OU	
 ##H#8##H#851))''S'1##H#;))''S'1##H#;51 299m,,))''-C3 -rH   r   r   rJ   rQ   rR   rS   rT   rU   N)	r#   r$   r%   r&   r   r)   rb   rc   rd   )rb  r9   r:   rM   xn_1dgkde_1dgkde_1d_weightedr'   r-   xn_2dgkde_2dgkde_2d_weighteds               rF   	test_seedrj  Z  s    44 ))


(CL	L	!BJJq!\*E  'G'))%<"# 88S#JDC:Sz23J##D<#HJJE  'G'))%<"#rH   c                     \         P                  P                  R 4      p \         P                  ! . RO4      p\         P                  ! . RO. RO. RO.4      pV P	                  WR4      p \
        P                  ! VP                  4       R#   \        P                   dh    Rp\        \        P                  TR7      ;_uu_ 4        \
        P                  ! TP                  4       RRR4        R#   + '       g   i      R# ; ii ; i)	l   SDG. i  z:The data appears to lie in a lower-dimensional subspace...r   N)r   r      )r   r   r   )r      r   )r   r   rD  )r#   r$   r%   rb   rc   r   r)   rd   r   LinAlgErrorr   )r9   musigmadatar   s        rF   %test_singular_data_covariance_gh10205rr    s     ))

 3
4C	+	BHHj+{;<E""2d3D'466" 'J6--S99tvv& :999's*   * B :D!C2'D2D=DDc                  n   \         P                  P                  R 4      p V P                  \         P                  ! ^4      \         P
                  ! ^4      ^R7      pRp\        P                  ! \        VR7      ;_uu_ 4        \        P                  ! V4       RRR4       R#   + '       g   i     R# ; i)l   &8*c rU   z9Number of dimensions is greater than number of samples...r   N)r#   r$   r%   rc   zeroseyer   r   r   r   r)   )r9   rvsr6  s      rF   )test_fewer_points_than_dimensions_gh17436rw    ss     ))

 3
4C

!
!"((1+rvvayq
!
ACIG	z	1	13 
2	1	1	1s   B##B4	)float32r   float96float128int32int64)0scipyr   r   r   numpyr#   numpy.testingr   r   r   r	   r
   r   r   r   r   rG   rN   markparametrizeparamxslowru   ry   r   r   r)   r   r   r   r   r   r   r   _ftypesr   r   r  r  r#  r.  r7  r?  rK  rT  rj  rr  rw  r   rH   rF   <module>r     s   * * L L  *$EN(EV "LLV[[5F5F4GH7J7Jt "LLV[[5F5F4GH:J:JzO4O<#U'' #"U'' "-U'' -
!4H9 4*> J Gw.D$DE'*4 + F4./0/4)6 ) )8"!J 8 8"C	@/$d'
 rH   