tKgdZddlZddlZddlmZddlmZm Z m Z m Z m Z m Z mZ d dddZej ed Zd e_ej edZd e_ d ddd Zej ed Zd e_ej edZde_ d dddZej ed Zde_ej edZde_d!dddZd!dddZ d"dddZej ed Zde_ej edZde_ d"dddZej ed Z de _ej edZ!de!_ d"dddZ"ej e"d Z#de#_ej e"dZ$de$_d#dZ%ej e%d Z&de&_ej e%dZ'de'_y)$z( Discrete Fourier Transforms - basic.py N) pypocketfft) _asfarray_init_nd_shape_and_axes _datacopied _fix_shape _fix_shape_1d_normalization_workersF)planc| tdt|}|xs t||}t||}t |}|t |||\}} |xs| }n0|j |dkrd|j |d} t| |r|jjdk(r|nd} tj||f||| |S)z@ Return discrete Fourier transform of real or complex sequence. NFPassing a precomputed plan is not yet supported by scipy.fft functionsrinvalid number of data points ( ) specifiedc) NotImplementedErrorrrr r r shape ValueErrordtypekindpfftc2c) forwardxnaxisnorm overwrite_xworkersr tmpcopiedmessageouts ^/home/alanp/www/video.onchill/myenv/lib/python3.12/site-packages/scipy/fft/_pocketfft/basic.pyrr s !#EF F A,C4S!!4K $ (DwG}#CD1 V!+V 41 3CIIdO3DKP!!#))..C"73TC 88C$'4g >>TfftifftcX| tdt|}t||}t|}t j |s t d|t|||\}} n.|j|dkrtd|j|dtj||f||d|S)z8 Discrete Fourier transform of a real sequence. Nrx must be a real sequencerrr) rrr r np isrealobj TypeErrorr rrrr2c rrrrrrrr r _s r$r-r-(s  !#EF F A,C $ (DwG << 344}sAt,Q 41 :399T?:K;WXX 88C$'4w ??r%rfftihfftc N| tdt|}t||}t|}t j |r|dz}|)|j |dz dz}|dkr%td|dt||dzdz|\}} tj||f|||d|S)zG Return inverse discrete Fourier transform of real sequence x. NrrzInvalid number of data points (r) rrr r r*r+rrr rc2rr.s r$r5r5Fs  !#EF F A,C $ (DwG ||CCi y YYt_q A % q5>qcMN NsQTQJ5Q 88C$!WdD' BBr%hfftirfftc<| tdt||||||S)z@ 2-D discrete Fourier transform of a Hermitian sequence r)rhfftnrsaxesrrrr s r$hfft2r=hs4  !#EF F AtT; 88r%c<| tdt||||||S)zH 2-D discrete inverse Fourier transform of a Hermitian sequence r)rihfftnr:s r$ihfft2r@ss4  !#EF F !QdK 99r%c^| tdt|}t|||\} }|xs t||}t |}t |dk(r|St || |\}} |xs| }t||}|r|jjdk(r|nd} tj||||| |S)z= Return multidimensional discrete Fourier transform. Nrrr) rrrrr lenrr rrrr) rrr;r<rrrr r rr!r#s r$c2cnrC~s  !#EF F A,C)#q$7KE44S!!4KwG 4yA~S%.KC'K $ (D#))..C"73TC 88Cwc7 ;;r%fftnifftncF| tdt|}tj|s t dt |||\} }t || |\}} t||}t|}t|dk(r tdtj||||d|S)z@Return multidimensional discrete Fourier transform of real inputNrr)r#at least 1 axis must be transformed) rrr*r+r,rrr r rBrrr-) rrr;r<rrrr r rr/s r$r2cnrHs !#EF F A,C << 344)#q$7KE4 UD )FC $ (DwG 4yA~>?? 88CwdG <?? KEWWT"X&*a/b $ (DwGRyHra1$E"I :c5$/ 0FC 88Cx$g FFr%r9irfftnc t|}|xs t||}t||}td}|jj dk(r t d|t|||\}}|xs|}n.|j|dkrtd|j|d|r|nd} tj||f|||| |S)z=FFT of a real sequence, returning fftpack half complex formatNrr)rrr) rrr r rrr,r rrr r2r_fftpack) rrrrrrr rr!r#s r$rSrSs A,C4S!!4K $ (DtnG yy~~344}#CD1 V!+V 41 :399T?:K;WXX34C   C$'7D#w OOr% rfft_fftpack irfft_fftpack)NrKNFN)N)rKNFN)NNNFN)NrKNF)(__doc__numpyr* functoolsrrhelperrrrrr r r rpartialr&__name__r'r-r0r1r5r6r7r=r@rCrDrErHrIr?rPr9rQrSrTrUr%r$r_sX!=B?"?.iT" ye$ =B@"@0yd#  #u%=BC"C8yd#  #u%99::@E<#<4yt$  $&?D=#=.  $%   4 '@E G# GF  $%   4 'P*!y  d3 & ! !!+u5 ( r%