
    YhE                     @   d dl Z d dlZd dlmZmZmZ d dlZd dlmZ d dl	m
Z
 ej        j        Zej        j        Z G d d          Zej        ej        eej        ej        ej        eeej        ej        eeiZ ed           G d	 d
                      Zd Zd ZdS )    N)assert_almost_equalxp_assert_closexp_assert_equal)raises)signalc                   H   e Zd ZdZ eddg          d             Z eddg          d             Zd Z edd	          d
             Z ed          d             Z	 ed          d             Z
 ed          d             Z ed          d             ZdS )TestBSplinesa  Test behaviors of B-splines. Some of the values tested against were
    returned as of SciPy 1.1.0 and are included for regression testing
    purposes. Others (at integer points) are compared to theoretical
    expressions (cf. Unser, Aldroubi, Eden, IEEE TSP 1993, Table 1).Tcupy)cpu_only
exceptionsc                    t           j                            d          }t          t          t
          j        |                    dg          d           |                    dd          }ddd|z  z
  z  }|                    |          }|                    g dg dg d	g d
g dg dg dg dg dg dg dg dg|j	                  }t          t          j        |d          |           d S )N0  r      
         )g7ݿgD @g5P?gB߸"@gyӼ®@gE5`@g&)#@gEf!g<@gSbυ.X gFgg$`#@)gI<M$k@gYvR@gd)B	g,#@gd?H#gk"`	@g)k	gi	V@gH-Q#@gn)Ҕ g3%0?g@	)gAZg9s@g /{Dg٩@gEC@g'UE@gآ?@gւ gzk@ga"@g򿽴~gKG#@)g$9gL|wglT@g?ss}aJ?g7iH^	gt,`?gh@glA1mg(T#g* #gDV#@g@)g'"@ga+ٓgE=a@g:"gIV@g5@gYd\gі1pg"2 g:i@gi^J @g`HQ!)g^!g#Mƽ!@g论H&@gfQq@g%gm̂r!@g
,@gY bgY#@g/-i@g Vqg")gH"@ggʜ6@ggeg&P@gY2#Cg
0@g|0g *!/@g(gqq?g)a)gF!Ígzǫg^geF?g
O0!gV<gtСg%9W@g^[H#@g3Fz?g.5 gsd?)g놎g,
@gA/u@g`9ogebgȀ@gMgZa@g<ˉ8g:ɘ%@gdH	g4g:Ns@)gg-ҿgE- g)in"@gd.xg2~wgďgYr4Q?gkM-:g,@g+ @gVu)g@+_$?gKGj?g[ C$}@g_gD/X@gHgJSʌ@g@ggN+gë	gk8gjZ@)gݵB`#@g+N?g	gg'g=!"F@gƱҟ"g'-:@g#g'gj \B@g~ge@g$q߼dtype)nprandomRandomStater   	TypeErrorr   spline_filterasarrayrandfloat64r   )selfxprngdata_array_realresult_array_reals        r/var/www/tools.fuzzalab.pt/emblema-extractor/venv/lib/python3.11/site-packages/scipy/signal/tests/test_bsplines.pytest_spline_filterzTestBSplines.test_spline_filter   s   i##E**y&.

A3CCC((2r**a/ 112**_55JJ@ @ @@ @ @? ? ?@ @ @@ @ @@ @ @@ @ @? ? ?@ @ @@ @ @A A A@ @ @C#AF IK
I ' $T $TJ 	,_a@@)	+ 	+ 	+ 	+ 	+    c           
         t           j                            d          }|                    dd          |                    dd          dz  z   }ddd|z  z
  z  }|                    |          }|                    g dg dg d	g d
g dg dg dg|j                  }t          t          j        |d          |d           d S )Nr      y              ?r   y      ?      ?r   )y$
ݿ%)y @;;R@yH?k៌"y?!@3⟥@yʱƯ@bȺ@ymտ5a@,y`#@A4p)yF! y韷@@_@yw.X /@ @y4 L &Q!yjٿ_#@!yߟk@a
 ƽ!@ySh?@ӿ&@)y?@	X@@yH:#@~%yXH#r!@y} "`	@)ƿ@y@k	ubyoU@I"$#@yA Q#@,i@)yb@r ߂KqyӔ?Ȥ"yr	"@y%`Zv.y$?s@*  7@yxebeyC@P@)y?M@k9FyθUE@4%@@y[y@#?0yN 5H /@ywk@y"a"@$?y~H)y}?\#@A?y	dy{_( y^p@;^F?yu?aJ?U 0!y\@^	Kyx_?+ߡ)yl.@@9UY@yNЪmxp[H#@yDT#bz?yh #F5 y#@^d?y<@V?y{?"@&@@r   r   gư>)rtol)	r   r   r   r   r   
complex128r   r   r   )r   r   r   data_array_complexresult_array_complexs        r"   test_spline_filter_complexz'TestBSplines.test_spline_filter_complexI   s/   i##E** XXa^^chhq!nnR.??a(:&:!:;ZZ(:;;!zz& & && & && & &% % %+ + +& & &% % %1&6 .0]9  *  <  <@ 	,-?CC,4	9 	9 	9 	9 	9 	9r$   c                     t          j        t          j        dd          d          sJ t	          t          j        |                    dg          d          |                    dg          d           d S )Nr   g<~?      ?r   g?g&.>atol)mathiscloser   gauss_spliner   r   )r   r   s     r"   test_gauss_splinezTestBSplines.test_gauss_splinet   sy    |F/1557HIIIII+BJJt,<,<a@@

J<00t	
 	
 	
 	
 	
 	
r$   z$deliberate: array-likes are acceptednp_onlyreasonc                 |    g d}t          t          j        |d          t          j        g d                     d S )N)              r8      )硫Y.?gye~?r;   )r   r   r2   r   r   )r   r   knotss      r"   test_gauss_spline_listz#TestBSplines.test_gauss_spline_list{   sN     "!!F/q99J'J'J'JKK	
 	
 	
 	
 	
r$   )r   c                    t          t          j        |                    dg                    |                    dg|j                             |                    g d|j                  }t          t          j        |                    g d          d          |           |                    g d|j                  }t          t          j        |                    g d                    |           d S )Nr   r9   r   )g{߮]?gzU%?gz@g6fw@gven'@r-   r   r:         r   )g|-?g>!<m @g6L@glMO@g>}@)r   r   	cspline1dr   r   r   )r   r   c1dc1d0s       r"   test_cspline1dzTestBSplines.test_cspline1d   s   (QC99

B4rz
::	< 	< 	<jj & & &-/Z  9 9 	(4D4D4D)E)EqII3OOOzz ' ' '.0j  : :(4D4D4D)E)EFFMMMMMr$   c                 8   t          t          j        |                    dg                    |                    dg|j                             t          t          t          j        |                    g d          d           t          t          t          j        |                    g d          d           |                    g d|j                  }t          t          j        |                    g d|j                            |           d S )Nr   r9   r   r?   r-   r8   )g)`O?gAk2 @g@go@g5ڕ@)r   r   	qspline1dr   r   r   
ValueErrorr   )r   r   q1d0s      r"   test_qspline1dzTestBSplines.test_qspline1d   s   (QC99

B4rz
::	< 	< 	< 	z6+RZZ8H8H8H-I-I2NNNz6+RZZ8H8H8H-I-I3OOOzz ' ' '.0j  : :RZZ(8(8(8
ZKKLLd	
 	
 	
 	
 	
r$   c           	         t          j        |                    ddg|j                  |                    dg|j                            }t	          ||                    dg|j                             t          j        |                    g d|j                  |                    g |j                            }t          ||                    g |j                             g d}|d         |d         z
  }g d}|                    g d          }t          j        |          }|                    g d	|j                  }t	          t          j        ||                    |          ||d         
          |           t          j        t          d          5  t          j        |                    g |j                  |                    dg|j                             d d d            d S # 1 swxY w Y   d S )Nr9   r   r   r-   r   r   
r   r   r   r:   r@   rA      r   &g      g      g      g      g      g      g      g      g       g      r8   g      r9         ?r-   g      ?       @      @      @g      @      @g      @g      @      @      @g      @g      @g      @       @g      !@g      "@g      #@g      $@g      %@g      &@g      '@g      (@g      )@
gw/@~jt@&1@Zd;@K7	@Zd@"~@jt@tV?K7A`@)&r^   SW@r]   7.O@r\   00/@g)-3/@rg   r\   rf   r]   re   r^   ^:@r_   4Ev@r`   c4~#.@ra   uyJ9@rb   /~	@rc   Y@@rd   rm   rc   rl   rb   rk   ra   rj   r`   ri   r_   rh   r^   re   dxx0+Spline coefficients 'cj' must not be empty.match)
r   cspline1d_evalr   r   r   r   rB   pytestr   rH   )	r   r   rxro   newxycjnewys	            r"   test_cspline1d_evalz TestBSplines.test_cspline1d_eval   sx   !"**b!WBJ*"G"G!zz2$bjzAAC C2::rd"*:==>>>!"**ZZZrz*"J"J!zz"BJz??A A2::b
:;;<<<---qTAaD[   JJ & & & ' 'a  zz : : : BD  M M 	!"bjj&6&62!A$GGG	
 	
 	
 ]:"OQ Q Q 	G 	G!"**Rrz*"B"B"$**cU"**"E"EG G G	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	Gs   )AHHHc                    t          t          j        |                    ddg          |                    dg                    |                    dg                     t	          t          j        |                    g d          |                    g                     |                    g                      g d}|d         |d         z
  }g d}|                    g d          }t          j        |          }|                    g d|j        	          }t          j        ||                    ||j        	          ||d         
          }t          ||           t          j        t          d          5  t          j        |                    g |j        	          |                    dg|j        	                     d d d            d S # 1 swxY w Y   d S )Nr9   r   rL   rM   r   rR   r[   )&r^   TW
@r]   oR@r\   bP@geW/@r   r\   r   r]   r~   r^    B@r_   HG(@r`   nfWΜ=@ra   %d	O@rb   u#	@rc   )C@rd   r   rc   r   rb   r   ra   r   r`   r   r_   r   r^   r~   r   rn   rq   rr   )
r   r   qspline1d_evalr   r   rG   r   ru   r   rH   )	r   r   rw   ro   rx   ry   rz   r{   rv   s	            r"   test_qspline1d_evalz TestBSplines.test_qspline1d_eval   sV   -bjj"a.A.A2::rdCSCSTT

B4((	
 	
 	
 	-bjj.D.DbjjQSnnUU

2	
 	
 	
 .--qTAaD[   JJ & & & ' 'a  zz : : : BD  M M !

4rz
22rad
 
 
 	4   ]:!NP P P 	G 	G!"**Rrz*"B"B"$**cU"**"E"EG G G	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	G 	Gs   AGG #G N)__name__
__module____qualname____doc__skip_xp_backendsr#   r+   r3   r=   rE   rJ   r|   r    r$   r"   r	   r	      so       H H
 t999/+ /+ :9/+b t999(9 (9 :9(9T
 
 
 d+QRRR
 
 SR
 t$$$	N 	N %$	N t$$$

 

 %$

 t$$$G G %$GB t$$$G G %$G G Gr$   r	   T)r5   c                   l   e Zd Zd Zd Zej                            dej	        e
ej        eej        eg          d             Z edd          ej                            dej	        e
ej        eej        eg          d                         Z edd          ej                            d	
          d                         Z edd          ej                            d
          ej                            dej	        e
ej        eej        eg          d                                     ZdS )TestSepfir2dc           	         |                     g d          }t          j                            dd          }|                     |          }t	          j        |||dd                     t          j        t          d          5  t	          j        |||dd                     d d d            n# 1 swxY w Y   t          j        t          d          5  t	          j        ||dd          |           d d d            n# 1 swxY w Y   t          j        t          d          5  t	          j        ||	                    |d	          |           d d d            n# 1 swxY w Y   t          j        t          d          5  t	          j        |||	                    |d	                     d d d            d S # 1 swxY w Y   d S )
Nr-   rT   rW   rT   r-   r&   	   r   z
odd lengthrr   r   object too deep)r   rP   )
r   r   r   r   r   sepfir2dru   r   rH   reshaper   r   filtimages       r"   test_sepfir2d_invalid_filterz)TestSepfir2d.test_sepfir2d_invalid_filter   s   zz33344	q!$$

5!!tT!""X... ]:\::: 	3 	3OE4abb222	3 	3 	3 	3 	3 	3 	3 	3 	3 	3 	3 	3 	3 	3 	3]:\::: 	3 	3OE48T222	3 	3 	3 	3 	3 	3 	3 	3 	3 	3 	3 	3 	3 	3 	3 ]:->??? 	D 	DOE2::dG#<#<dCCC	D 	D 	D 	D 	D 	D 	D 	D 	D 	D 	D 	D 	D 	D 	D]:->??? 	D 	DOE4D')B)BCCC	D 	D 	D 	D 	D 	D 	D 	D 	D 	D 	D 	D 	D 	D 	D 	D 	D 	DsH   B11B58B5DD	D(+EE#&E#+F==GGc                     |                     g d          }t          j                            dd          }|                     |          }t	          j        t          d          5  t          j        |	                    |d          ||           d d d            n# 1 swxY w Y   t	          j        t          d          5  t          j        |dd d f         ||           d d d            d S # 1 swxY w Y   d S )Nr      r   rr   )r@   r@   r@   zobject of too small depthr   )
r   r   r   r   ru   r   rH   r   r   r   r   s       r"   test_sepfir2d_invalid_imagez(TestSepfir2d.test_sepfir2d_invalid_image   s|   zz33344	q!$$

5!! ]:->??? 	F 	FOBJJui88$EEE	F 	F 	F 	F 	F 	F 	F 	F 	F 	F 	F 	F 	F 	F 	F ]:-HIII 	5 	5OE!QQQ$Kt444	5 	5 	5 	5 	5 	5 	5 	5 	5 	5 	5 	5 	5 	5 	5 	5 	5 	5s$   (+BB#&B#!C33C7:C7dtypc                    t          j        g dg dg dg dg|          }g d}dg}t          j        |||          }t          |         }t          j        g dg dg dg dg|          }t          ||d           t          j        |||          }t          j        g dg dg dg dg|          }t          ||d           d S )	N)r   r   r:   r:   r   r   r   rS   r   rS   r   )rU   rW   rX   rX   rW   rU   gؗҜ<r.   )rT   rW   rY   rY   rW   rT   )r   arrayr   r   sepfir_dtype_mapr   r   )	r   r   r   ah1h2resultdtexpecteds	            r"   test_simplezTestSepfir2d.test_simple  sP   
 H((((((((((((* 267 7 7 ]]SB++d#:============? GIJ J J 	u5555B++:7777777777779 ACD D D 	u555555r$   TzTODO: convert this testr4   c                 L   t          j        g dg dg dg dg          }g ddg}}t          j        |d d d d df         ||          }t          j        |d d d d df                                         ||          }t          ||d           |j        |j        k    sJ d S )N)	r   r   r:   r:   r   r   r   r   r:   r   r   r   V瞯<r.   )r   r   r   r   copyr   r   )r   r   r   r   r   r   result_stridedresult_contigs           r"   test_stridedzTestSepfir2d.test_strided"  s    
 H1110000000002 3 3 B111ccc6B;;!!!SSqS&	(8(8"bAAEBBBB#}':::::::r$   z#XXX: filt.size > image.shape: flaky)r6   c           
         t          j        g d          }t           j                            dd          }t          j        g dg dg dg dg          }t          t          j        |||d d d                   |           d S )N)r-   rT   rW   rT   r-   rV   rT   r@   )g!SB@g¿3=>@gs`[C@gsLhE@)gͨ*C@g5A@g;/K;E@g<f2E@)gE@g>#jD@gB`G@g]n0a#F@)gٕzH@g|FC@gǜgKE@gCF@r:   )r   r   r   r   r   r   r   r   )r   r   r   r   r   s        r"   test_sepfir2d_strided_2z$TestSepfir2d.test_sepfir2d_strided_21  s    
 x;;;<<	q!$$:KKK L L L L L L L L LN O O 	tT##A#Y??JJJJJr$   z*XXX: flaky. pointers OOB on some platformsc                 J   t          j        g d|          }t          j        g dg dg dg dg dg|          }g dg d	g d
g dg dg}t          j        |          }t          j        |||d d d                   }t          ||d           |j        t          |         k    sJ g dg dg dg dg dg}t          j        |          }t          j        ||d d d         |d d d                   }t          ||d           |j        t          |         k    sJ d S )N)r   r   r@   r   r   r:   r   r   )r   r:   r   r   r   )r   r   r:   r:   r:   )r   r   r:   r   r:   )r   r:   r   r   r:   )r:   r:   r   r   r   )g     ^@g     @Y@g     V@      a@g     _@)     `@g     @_@g     _@      c@g      d@)r   g      a@g     b@g     @d@g      f@)r   g      _@g     `@g     b@g     `b@)g     e@g     c@r   g     d@g     a@r:   r   r.   )g      6@g     A@     D@g      ?@     G@)g      ;@     C@      H@r   g     K@)g     @@g      E@g     H@g     J@g     M@)r   g      F@r   g      B@r   )g     P@g      O@r   g      A@g      G@)r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   s          r"   test_sepfir2d_strided_3z$TestSepfir2d.test_sepfir2d_strided_3?  s    x---T:::
OOO+OO+OO+OO+OO	- 59	: : : 322222222222222	4
 :h''d33Q3i88u5555|/55555---------------	/
 :h''SSqS	4!9==u5555|/5555555r$   N)r   r   r   r   r   ru   markparametrizer   uint8intfloat32float	complex64complexr   r   r   xfailr   r   r   r$   r"   r   r      s       D D D&
5 
5 
5 [V	3
E2<A 6 6 6. d+DEEE[V	3
E2<A 	; 	;  FE	; d+DEEE[CDD
K 
K ED FE
K d+DEEE[JKK[V	3
E2<A 6 6  LK FE
6 6 6r$   r   c                     t           j                            d          }|                    dd          }t	          j        |d           d S )NiW
G   I   rZ   )r   r   r   r   r   	cspline2dr   r   r   s      r"   test_cspline2dr   d  sC    
)

	
*
*CHHRE
UC     r$   c                     t           j                            d          }|                    dd          }t	          j        |           d S )NiW
r   r   )r   r   r   r   r   	qspline2dr   s      r"   test_qspline2dr   j  sA    
)

	
*
*CHHRE
Ur$   )r0   numpyr   scipy._lib._array_apir   r   r   ru   r   scipyr   r   r   xfail_xp_backendsr	   r   r   r   r   r   r   r   r   r   r   r   r   r$   r"   <module>r      sr                             ;/ K1 MG MG MG MG MG MG MG MGb Hbj#rzJ
E5L",B 
 $y6 y6 y6 y6 y6 y6 y6  y6x! ! !    r$   