U ʹj*C@sddlZddlmZmZmZmZmZmZmZm Z ddl Z ddl m Z ddlmZddlmZddlmZmZmZddlmZmZdddd d d d d ddddddddddddgZe dedededfZddddhZe jjjZeeee edddZ!e j"e e j"d d!d"Z#d_ee ed$d%d&Z$eeed'feed'fed(d)d*Z%e&eeeeee ed+d,d-Z'e&eeeeee e ed.d/d0Z(e&eeeeee ed+d1d2Z)eej*ed`eeeeeed4d5dZ+eej,edaeeeeeed4d6dZ-eej.edbeeeeeed4d7d Z/eej0edceeeeeed4d8dZ1eej2eddeeeeeed4d9d Z3eej4edeeeeeeed4d:dZ5Gd;d<dd?Z7eeed@dAdBZ8e&eeed'feed'fee edCdDdEZ9eej:edfeeeeeeedFdGdZ;eej<edgeeeeeeedFdHdZ=eej>edheeeeeeedFdIdZ?eej@edieeeeeeedFdJdZAGdKdLdLeZBe&eeeeeeBdMdNdOZCeejDedjeeeeeeedFdPdZEeejFedkeeeeeeedFdQd ZGeejHedleeeeeeedFdSdZIeejJedmeeeeeeedFdTdZKeejLedneeeeeeedFdUd ZMeejNedoeeeeeeedFdVdZOeejPedpeeeeeeedFdWd ZQeejRedqeeeeeeedFdXdZSeeeeedYdZd[ZTeejUdreeeed\d]dZVeejWdseeeed\d^dZXdS)tN)IterableListLiteral NamedTupleOptionalSequenceTupleUnion)register_decomposition)DimsType ShapeTypeTensorLikeType)_maybe_convert_to_dtype out_wrapperfftfft2fftnhffthfft2hfftnrfftrfft2rfftnifftifft2ifftnihfftihfft2ihfftnirfftirfft2irfftnfftshift ifftshiftforwardbackwardortho)xnorm signal_numelr$returncshttkfdddkr2|dt|S| rHdkpRdkpR|oRdk}|rd|d|S|S)z3Apply normalization to the un-normalized FFT resultcs dS)NzInvalid normalization mode: r+r(r+N/var/www/html/bwave/python/venv/lib/python3.8/site-packages/torch/_refs/fft.py.z_apply_norm..r&Nr%r$)torch_check _NORM_VALUESmathsqrt)r'r(r)r$ normalizer+r,r- _apply_norm*s r7)dtyperequire_complexr*cCs*|jr |S|jst}|r&t|}|S)z@Helper to promote a dtype to one supported by the FFT primitives) is_complexZis_floating_pointr1Zget_default_dtypeutilsZcorresponding_complex_dtype)r8r9r+r+r-_promote_type_fft9s r<F)tr9r*cCs|j}t||}t||S)zEHelper to promote a tensor to a dtype supported by the FFT primitives)r8r<r)r=r9Zcur_typenew_typer+r+r-_maybe_promote_tensor_fftHs r?.)r'dimssizesr*cCst|t|kstd}|j}dgt|d}tt|D]}||dkrNq<|||||krd}t|d||d}|||||||<|||||kr<|||d||}q<|rt||S|S)z Fixes the shape of x such that x.size(dims[i]) == sizes[i], either by zero-padding, or by slicing x starting from 0. FrTr0)lenAssertionErrorshaperangeZnarrowr1Zconstant_pad_nd)r'r@rAZ must_copyZx_sizesZ pad_amountiZpad_idxr+r+r-_resize_fft_inputQs rI) func_nameinputndimr(r$r*c st|dd}tj|j|ddf}dk r,nd|j|d}t|dkfdd dk rvt|||ddfd }|rt|}t j |||d }t ||||d S) zBCommon code for performing any complex to real FFT (irfft or hfft)Tr9FZ wrap_scalarNrBr0cs ddSNzInvalid number of data points (z ) specifiedr+r+rLr+r-r.xr/z_fft_c2r..)r@rArM last_dim_sizer(r)r$) r?r;canonicalize_dimndimrFr1r2rIconjprimsfft_c2rr7) rJrKrLrMr(r$r@rSoutputr+rQr-_fft_c2rks   r[)rJrKrLrMr(r$onesidedr*c stjj fddttjj|ddf}|dk rNt||ft j ||d}t ||j ||}|rz|St |S)zBCommon code for performing any real to complex FFT (rfft or ihfft)csdjS)Nz0 expects a floating point input tensor, but got r8r+rJrKr+r-r.r/z_fft_r2c..FrONrMr\)r1r2r8r:r?r;rUrVrIrXfft_r2cr7rFrW) rJrKrLrMr(r$r\r@retr+r^r-_fft_r2cs  rbcshtjjfddtjj|ddf}|dk rDt||ftj ||d}t ||j ||S)zCCommon code for performing any complex to complex FFT (fft or ifft)csdjSNz) expects a complex input tensor, but got r]r+r^r+r-r.r/z_fft_c2c..FrONrMr$) r1r2r8r:r;rUrVrIrXfft_c2cr7rF)rJrKrLrMr(r$r@rar+r^r-_fft_c2cs  rfrC)rKrLrMr(r*c Cs6|jjrtd||||ddStd||||dddSdS)NrTr$Fr$r\r8r:rfrbrKrLrMr(r+r+r-rsc Cs6|jjrtd||||ddStd||||dddSdS)NrFrgrhrirjr+r+r-rsc Cstd||||dddS)NrTrhrbrjr+r+r-rscCstd||||ddS)NrFrgr[rjr+r+r-rscCstd||||ddS)NrTrgrlrjr+r+r-rsc Cstd||||dddS)NrFTrhrkrjr+r+r-rsc@s.eZdZUeedfed<eedfed<dS) _ShapeAndDims.rFr@N__name__ __module__ __qualname__rint__annotations__r+r+r+r-rms rm)rKrFrMr*csP|j|j|dk rTt|ts$|f}tj|dd}ttt |t|kdd|dk rt|tsl|f}t|dkpt|t|kddt|tkfdd|dkrt t }t fdd t ||D}n6|dkr t t }t }nt fd d |D}|D]td kfd dq$t ||d S)zTConvert the shape and dim arguments into a canonical form where neither are optionalNFrOcSsdS)NzFFT dims must be uniquer+r+r+r+r-r.r/z6_canonicalize_fft_shape_and_dim_args..cSsdS)Nz=When given, dim and shape arguments must have the same lengthr+r+r+r+r-r.r/csdddS)NzGot shape with z" values but input tensor only has z dimensions.r+r+) input_dimtransform_ndimr+r-r.r/c3s&|]\}}|dkr|n|VqdS)rCNr+).0sd input_sizesr+r- &sz7_canonicalize_fft_shape_and_dim_args..c3s|]}|VqdSNr+rvrxryr+r-r{/srcs ddSrPr+r+rQr+r-r.2r/)rFr@)rVrF isinstancerr;Zcanonicalize_dimsr1r2rDsettuplerGziprm)rKrFrMZret_dimsZ ret_shaper+)rtrzrLrur-$_canonicalize_fft_shape_and_dim_argss>         r)xsr*cCsd}|D] }||9}q|S)zCompute product of a listr0r+)rprodr'r+r+r-_prod7s r) function_namerKrFrMr(r$r*csJtjjfddt||}tj|||d}t||t||dS)zECommon code for n-dimensional complex to complex FFTs (fftn or ifftn)csdjSrcr]r+rrKr+r-r.Jr/z_fftn_c2c..rdrT) r1r2r8r:rIrXrer7r)rrKrFrMr(r$r'rZr+rr- _fftn_c2c?s   r)rKrwrMr(r*cCs0t|||\}}t|dd}td||||ddS)NTrNrrgrr?rrKrwrMr(rFr'r+r+r-rRs cCs0t|||\}}t|dd}td||||ddS)NTrNrFrgrrr+r+r-r_s csftjj fddt||\}}tddt||tj|dd}t ||t |ddS)Ncs djS)Nz2rfftn expects a real-valued input tensor, but got r]r+rKr+r-r.vr/zrfftn..FrNTr_rT) r1r2r8r:rr?rIrXr`r7r)rKrwrMr(rFoutr+rr-rls   cstjj fddt||\}}tt|dkddtddt||tj |dddd }t|d krt |||ddd }t |St |}tj ||dddd }t ||t|dd S) Ncs djS)Nz3ihfftn expects a real-valued input tensor, but got r]r+rr+r-r.r/zihfftn..rcSsdS)Nz'ihfftn must transform at least one axisr+r+r+r+r-r.r/FrNrCTr_r0rTrd)r1r2r8r:rrDr?rIrXr`r7rW conj_physicalrer)rKrwrMr(rFtmpr+rr-rs      c@s6eZdZUeedfed<eedfed<eed<dS)_CanonicalizeC2rReturn.rFrMrSNrnr+r+r+r-rs r)fnamerKrwrMr*cst|||\}}tt|dkfdd|dks@|ddkrXd|j|ddn|dtdkfddt|}dd|d<tt||d S) zCanonicalize shape and dim arguments for n-dimensional c2r transforms, as well as calculating the last_dim_size which is shape[dim[-1]] for the outputrcs dS)Nz! must transform at least one axisr+r+)rr+r-r.r/z:_canonicalize_fft_c2r_shape_and_dim_args..NrCrBr0cs ddSrPr+r+)rSr+r-r.r/)rFrMrS)rr1r2rDrFlistrr)rrKrwrMrFZ shape_listr+)rrSr-(_canonicalize_fft_c2r_shape_and_dim_argss  rcs^td|||\}}}t|dd}t|||}tj|||dt|tfdd|DddS) Nr!TrNrRc3s|]}j|VqdSr|rFr}rr+r-r{szirfftn..Frg)rr?rIrXrYr7r)rKrwrMr(rFrSr+rr-r!s   cCstd|||\}}}t|dd}t|||}t|dkrPtj||ddddn|}t||t|dddd}t|}tj ||dd|d}t|||ddS) NrTrNr0rCrdrgrR) rr?rIrDrXrer7rrrY)rKrwrMr(rFrSrrr+r+r-rs   ( rCcCstjj||||dSN)rwrMr()r1rrrKrwrMr(r+r+r-rscCstjj||||dSr)r1rrrr+r+r-rscCstjj||||dSr)r1rrrr+r+r-rscCstjj||||dSr)r1rr!rr+r+r-r scCstjj||||dSr)r1rrrr+r+r-rscCstjj||||dSr)r1rrrr+r+r-rs)rMr'r*cCs2|dkrtt|jSt|ts&|gSt|SdS)zIConvert Optional[DimsType] to a simple list, defaulting to all dimensionsN)rrGrVr~r)rMr'r+r+r-_default_alldims%s  r)rKrMr*cs*t|}fdd|D}t||S)Ncsg|]}j|dqS)rBrr}rr+r- 2szfftshift..rr1ZrollrKrMr@shiftr+rr-r"/s cs*t|}fdd|D}t||S)Ncsg|]}j|ddqS)r0rBrr}rr+r-r9szifftshift..rrr+rr-r#6s )F)NrCN)NrCN)NrCN)NrCN)NrCN)NrCN)NNN)NNN)NNN)NNN)NNN)NNN)NrN)NrN)NrN)NrN)NrN)NrN)N)N)Yr4typingrrrrrrrr r1Z torch._primsZ_primsrXZtorch._prims_commonZ _prims_commonr;Z torch._decompr r r r Ztorch._prims_common.wrappersrr__all__ZNormTyper3Z_opsZopsZatenrrboolr7r8r<r?rIstrr[rbrfZfft_fftrZfft_ifftrZfft_rfftrZ fft_irfftrZfft_hfftrZ fft_ihfftrrmrrrZfft_fftnrZ fft_ifftnrZ fft_rfftnrZ fft_ihfftnrrrZ fft_irfftnr!Z fft_hfftnrZfft_fft2rZ fft_ifft2rZ fft_rfft2rZ fft_irfft2r Z fft_hfft2rZ fft_ihfft2rrZ fft_fftshiftr"Z fft_ifftshiftr#r+r+r+r-sz(                   6