o
    h})                     @   s   d dl Zd dl mZ d dlmZmZmZmZ d dlZd dlm	Z	 d dl
m  mZ d dlmZ G dd dZdd	 Zd
d Zdd Zdd ZdS )    N)array)assert_allcloseassert_array_equalassert_almost_equalsuppress_warnings)raises)signalc                   @   s`   e Zd ZdZdd Zdd Zdd Zdd	 Zd
d Zdd Z	dd Z
dd Zdd Z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).c                 C   s   t jd tttjtdgd t jd t jdd}ddd|   }t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}t	t|d| d S )Ni0  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߼)
nprandomseedr   	TypeErrorbspspline_filterr   randr   )selfdata_array_realresult_array_real r   [/var/www/html/scripts/venv/lib/python3.10/site-packages/scipy/signal/tests/test_bsplines.pytest_spline_filter   s.   %zTestBSplines.test_spline_filterc                 C   s*  t  }|t ttg ddtg d ttg ddtg d ttg ddtg dd  ttg dd	tg d
d  ttg ddtg dd  ttg ddtg dd  tj	d tttj
dddtdgg W d    d S 1 sw   Y  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   L      r&   r   r   g      x@   )r   r      B   r)   r   r   g      ^@i0  ga$?)r   filterDeprecationWarningr   r   bspliner   r   r   r   r   r   supr   r   r   test_bsplineD   s2   



"zTestBSplines.test_bsplinec                 C   s@   t jd ttddd tttdgdtdg d S )Ni0  r   g<~?      ?r   g?)r   r   r   r   r   gauss_spliner   r   )r   r   r   r   test_gauss_splineY   s   "zTestBSplines.test_gauss_splinec                 C   s&   g d}t t|dtg d d S )N)              r4   r!   )硫Y.?gye~?r6   )r   r   r2   r   )r   knotsr   r   r   test_gauss_spline_list^   s   
z#TestBSplines.test_gauss_spline_listc                 C   sb   t  $}|t tjd ttg dt	g dd  W d    d S 1 s*w   Y  d S )Ni0  r   r"   r$   )
r   r+   r,   r   r   r   r   r   cubicr   r.   r   r   r   
test_cubicd   s   
"zTestBSplines.test_cubicc                 C   sb   t jd t }|t ttg dt	g dd  W d    d S 1 s*w   Y  d S )Ni0  r   r   r    )
r   r   r   r   r+   r,   r   r   	quadraticr   r.   r   r   r   test_quadraticl   s   
"zTestBSplines.test_quadraticc                 C   sr   t jd tttdgdg tg d}tttg dd| tg d}tttg d| d S )Ni0  r   r5   )g{߮]?gzU%?gz@g6fw@gven'@r1   r   r!   r#   r(   r   )g|-?g>!<m @g6L@glMO@g>}@)r   r   r   r   r   	cspline1dr   r   )r   c1dc1d0r   r   r   test_cspline1dt   s   zTestBSplines.test_cspline1dc                 C   s|   t jd tttdgdg tttjtg dd tttjtg dd tg d}t	ttg d| d S )Ni0  r   r5   r=   r1   r4   )g)`O?gAk2 @g@go@g5ڕ@)
r   r   r   r   r   	qspline1dr   r   
ValueErrorr   )r   q1d0r   r   r   test_qspline1d   s   zTestBSplines.test_qspline1dc                 C      t jd tttddgdgtdg tttg dg tg  g d}|d |d  }g d}tg d}t|}tg d	}ttj||||d d
| d S )Ni0  r5   r   r1   r   r   
r%   r   r   r   r   r   r!   r#   r(   r   r   &g      g      g      g      g      g      g      g      g       g      r4   g      r5   g      ?r1   g      ?       @g      @g      @g      @      @g      @g      @g      @r$   g      @g      @g      @r    g      !@g      "@g      #@g      $@g      %@g      &@g      '@g      (@g      )@
gw/@~jt@&1@Zd;@K7	@Zd@"~@jt@tV?K7A`@)&rO   SW@rN   7.O@rM   00/@g)-3/@rX   rM   rW   rN   rV   rO   ^:@rP   4Ev@rQ   c4~#.@rR   uyJ9@rS   /~	@rT   Y@@rU   r^   rT   r]   rS   r\   rR   r[   rQ   rZ   rP   rY   rO   rV   dxx0)	r   r   r   r   r   cspline1d_evalr   r   r>   r   xr`   newxycjnewyr   r   r   test_cspline1d_eval      "
 z TestBSplines.test_cspline1d_evalc                 C   rF   )Ni0  r5   r   rG   rH   r   rI   rL   )&rO   TW
@rN   oR@rM   bP@geW/@rm   rM   rl   rN   rk   rO    B@rP   HG(@rQ   nfWΜ=@rR   %d	O@rS   u#	@rT   )C@rU   rs   rT   rr   rS   rq   rR   rp   rQ   ro   rP   rn   rO   rk   r_   )	r   r   r   r   r   qspline1d_evalr   r   rB   rc   r   r   r   test_qspline1d_eval   rj   z TestBSplines.test_qspline1d_evalN)__name__
__module____qualname____doc__r   r0   r3   r8   r:   r<   rA   rE   ri   ru   r   r   r   r   r	      s    1
r	   c               	   C   sH  t g d} t jdd}t|| | dd   tjtdd t|| | dd   W d    n1 s6w   Y  tjtdd t|| dd  |  W d    n1 sXw   Y  tjtdd t|| 	dd	|  W d    n1 szw   Y  tjtdd t|| | 	dd	 W d    d S 1 sw   Y  d S )
Nr1   rJ   rK   rJ   r1      	   r   z
odd lengthmatchr   object too deepr   )
r   r   r   r   r   sepfir2dpytestr   rC   reshapefiltimager   r   r   test_sepfir2d_invalid_filter   s   "r   c                  C   s   t g d} t jdd}tjtdd t|	ddd| |  W d    n1 s,w   Y  tjtdd t|d | |  W d    d S 1 sMw   Y  d S )Nrz      r   r}   r#   zobject of too small depthr   )
r   r   r   r   r   r   rC   r   r   r   r   r   r   r   test_sepfir2d_invalid_image   s   "r   c                  C   s*   t jd t jdd} t| d d S )NiW
G   I   r    )r   r   r   r   r   	cspline2dr   r   r   r   test_cspline2d   s   r   c                  C   s(   t jd t jdd} t|  d S )NiW
r   r   )r   r   r   r   r   	qspline2dr   r   r   r   test_qspline2d   s   r   )numpyr   r   numpy.testingr   r   r   r   r   r   scipy.signal._bsplinesr   	_bsplinesr   scipyr	   r   r   r   r   r   r   r   r   <module>   s    ,