S_fL<dZddlmZmZddlmZddlmZddlZddl Z ddl Z ddl m Z ddlmZmZddlmZe ej&j)d gZej&j,Zgd Zd Zd Zed GddZed GddZeddgGddZeddddgGddZeddddgGdd Zy)!zIncludes test functions for fftpack.helper module Copied from fftpack.helper by Pearu Peterson, October 2005 Modified for Array API, 2023 ) next_fast_len_init_nd_shape_and_axes) assert_equal)raisesN)array_api_compatible)xp_assert_close SCIPY_DEVICE)fftskip_if_array_api)  iicDtD]}tt||y)N)_5_smooth_numbersrr)ns ;lib/python3.12/site-packages/scipy/fft/tests/test_helper.pytest_next_fast_lenrs *]1%q)*c|}|dkrJ t|d\}}|dk7rn|}td|dzdD]} t||\}}|dk7r|}|dk(sJd|d|d|y)Nr rr zx=z is not z-smooth, remainder=)divmodrange)xrx_origqrds r_assert_n_smoothr$ s F1uu a|1 6    1ac1 !Q.nums<s Q~ % %( ( &s   FTr)nprandomseedrr$)selfr+rms rrz"TestNextFastLen.test_next_fast_len9si t ) #Aa A Q # a// //a&A Q "  #rc ,tjtjtjtjtj tj g}|D]3}|d}t|}t|tt|5y)Ni90) r-int16int32int64uint16uint32uint64rrint)r0ITYPESityprtestNs rtest_np_integersz TestNextFastLen.test_np_integersHsa((BHHbhh 299biiP 7DU A!!$E  c!f 5 6 7rcddddddddddd d d d d d}|jD]\}}tt|d|y)Nrr r rrrri@ )rr r rrrrr?r@irBrCTitemsrrr0hamsrys rtestnext_fast_len_smallz'TestNextFastLen.testnext_fast_len_smallOsUQ1qQ1$dh JJL 4DAq q$/ 3 4rlz$Hamming Numbers too large for 32-bitT)reasonrstrictc iddddddddddddd d d d d d dddddddddddddddddddddddddddd }|jD]\}}tt|d!|y)"Nihihii=ui+2i,2iB3l l ls'l>9l>9l lfnlfnl>Ukolq. lq. lxQ+t0 l!l!ldopL?lhRlhR2 f;`71Y!l("  X l ) rPl2 f;rQlrRl`71Y!rSl l X rTl! X TrGrIs rtestnext_fast_len_bigz%TestNextFastLen.testnext_fast_len_bigWs' y "/ y  y #0  ,  )9,  =  +<]  _     _    .  !"2  0  "#4 1#40#40#4 2$6$6 2$69 <JJL 4DAq q$/ 3 4rcLtdddk(sJtdddk(sJy)Nr,T)real rDF)targetrW)r)r0s rtest_keyword_argsz!TestNextFastLen.test_keyword_args|s,Rd+r111AE2a777rN)__name__ __module__ __qualname__rr=rLpytestmarkxfailsysmaxsize ValueErrorrUrZr*rrr'r'6sR #74 [[s{{U*D(7 47 4D8rr')cpu_onlycTeZdZdZdZdZdZdZdZdZ dZ d Z d Z d Z d Zy )Test_init_nd_shape_and_axescr|jd}d}d}d}g}t|||\}}||k(sJ||k(sJy)Nrr*asarrayr r0xprshapeaxesshape_expected axes_expected shape_resaxes_ress rtest_py_0d_defaultsz/Test_init_nd_shape_and_axes.test_py_0d_defaultssR JJqM 5aE 8N***=(((rcr|jd}d}d}d}g}t|||\}}||k(sJ||k(sJy)Ng@r*rhrjs rtest_xp_0d_defaultsz/Test_init_nd_shape_and_axes.test_xp_0d_defaultssR JJrN 5aE 8N***=(((rcx|jgd}d}d}d}dg}t|||\}}||k(sJ||k(sJy)Nrr r )r rrhrjs rtest_py_1d_defaultsz/Test_init_nd_shape_and_axes.test_py_1d_defaultssU JJy ! 5aE 8N***=(((rcx|jddd}d}d}d}dg}t|||\}}||k(sJ||k(sJy)Nrr皙?)r)arangerrjs rtest_xp_1d_defaultsz/Test_init_nd_shape_and_axes.test_xp_1d_defaultssY IIaB  5aE 8N***=(((rc|jgdgdg}d}d}d}ddg}t|||\}}||k(sJ||k(sJy)N)rr r r)rrrDr)r rrrrhrjs rtest_py_2d_defaultsz/Test_init_nd_shape_and_axes.test_py_2d_defaultssa JJ $& 'A 5aE 8N***=(((rc|jddd}|j|d}d}d}d}ddg}t|||\}}||k(sJ||k(sJy)Nrrry)rr )rzreshaperrjs rtest_xp_2d_defaultsz/Test_init_nd_shape_and_axes.test_xp_2d_defaultssk IIaB  JJq& !A 5aE 8N***=(((rcz|jgd}d}d}d}gd}t|||\}}||k(sJ||k(sJy)Nrr rr rrrr r rzerosrrjs rtest_xp_5d_defaultsz/Test_init_nd_shape_and_axes.test_xp_5d_defaultssS HH_ %(' 5aE 8N***=(((rc~|jgd}gd}d}d}gd}t|||\}}||k(sJ||k(sJy)Nr)rrrr)rr rrrrrrjs rtest_xp_5d_set_shapez0Test_init_nd_shape_and_axes.test_xp_5d_set_shapesS HH_ %")' 5aE 8N***=(((rc~|jgd}d}gd}d}gd}t|||\}}||k(sJ||k(sJy)Nr)rrr )rr rrrjs rtest_xp_5d_set_axesz/Test_init_nd_shape_and_axes.test_xp_5d_set_axessS HH_ %"! 5aE 8N***=(((rc|jgd}gd}gd}d}gd}t|||\}}||k(sJ||k(sJy)Nr)rrr )rrr )rrr rrjs rtest_xp_5d_set_shape_axesz5Test_init_nd_shape_and_axes.test_xp_5d_set_shape_axessS HH_ %#! 5aE 8N***=(((rch|jd}t|dd\}}|dk(sJ|gdk(sJy)N)r r rr)rrrrlrmrvr)r0rkrrlrms rtest_shape_axes_subsetz2Test_init_nd_shape_and_axes.test_shape_axes_subsets> HH\ "-aytL t !!!y   rc|jd}ttd5t|dddgddggdddttd5t|dgddddttd 5t|ddgdddttd 5t|dd gdddttd 5t|dd d gdddttd 5t|ddgddggddddttd 5t|gdddddttd5t|jgdgddgdddttd5t|d gddddttd5t|d ddddy#1swYxYw#1swYyxYw#1swY]xYw#1swYAxYw#1swY$xYw#1swYxYw#1swYxYw#1swYxYw#1swYxYw#1swYyxYw)Nrz-axes must be a scalar or iterable of integers)matchr r rr)g?g@g@g@z$axes exceeds dimensionality of inputzall axes must be uniquerz.shape must be a scalar or iterable of integerszBwhen given, axes and shape arguments have to be of the same length)rrrrrvz1invalid number of data points \(\[0\]\) specifiedz2invalid number of data points \(\[-2\]\) specified)r assert_raisesrcr)r0rkrs r test_errorsz'Test_init_nd_shape_and_axes.test_errorssQ HHQK :.23 J #AT!Q!Q8H I J:.23 J #AT8H I J:!GI = #AT < =:!GI > #AT = >:!:< @ #ATA ? @:.23 J #Aq!fq!f-=D I J:.23 J #A-=D I J:"<= ? $BHH\$:*31# ? ? :"34 = $AaSt < = :"45 < $ARd ; <  > @ @ J J J J ? ?  = =  <  GG#&G03G=H  HH#&H/2H;>IN)r[r\r]rrrtrwr{r}rrrrrrrr*rrrfrfs> ) ) ) ) ) ) ) ) ) )!, ?rcD|jgdgdgdg}|jgdgdgdg}ttj|d|ttj|d tj|d ttj|d|ttj|d tj|d ttj||ttj||y) N)rrr )r rr)rrr)grr)r rr)rr rrrrmrrr)r0rkfreqsshifteds rtest_axes_keywordzTestFFTShift.test_axes_keywordQs J LAB**mY CD U8'B U3S\\%d5ST gF;UC gA6 gD9 ; U+W5 g.6rc|jddgddgddgg|j}|jddgddgddgg|j}ttj|d|ttj |d|ttj|d |ttj |dg||jddgddgddgg|j}ttj|d|ttj |d||jddgddgddgg|j}ttj|d |ttj |d |ttj|ddg|ttj |ddg|ttj|d |ttj |d |ttj||ttj ||y ) z1 Test 2D input, which has uneven dimension sizes rrr r rrdtyperrrN)rifloat64rr rr)r0rkr shift_dim0 shift_dim1shift_dim_boths rtest_uneven_dimszTestFFTShift.test_uneven_dims\s F F F  ZZ F F F!     U3Z@ jq95A U6 C js;UCZZ F F F!     U3Z@ jq95A F F F%  $  U8.I n6BEJ U!Q8.I nAq6BEJ  U6G n4@%H U+^< n5u=rN)r[r\r]rrrrr*rrrr>s-? 7+>rrarray_api_strictcupyz+fft not yet implemented by array-api-strictz0cupy.fft not yet implemented by array-api-compatceZdZdZy) TestFFTFreqct} |jgd|j|}|jgd|j|}|jdt j d|z|j}t |||jd|jzt j d|j|z|j}t |||jdt j d|z|j}t |||jd|jzt j d|j|z|j}t ||y#t$rB|jgd|j}|jgd|j}YwxYw)N) rrr r rrrrrrdevice) rrr r rrrrrrrrrkr)r rir TypeErrorr fftfreqrpir0rkrrx2rKs rrzTestFFTFreq.test_definitionsm . :!#Frs6 E&*  /?"FKK$;$;G>RG>T%vINPQQ2%vINPQQr