`i ddlZejdddddZdd Zejd d d d dddZddZejddddddZdZejddddddZd dZ dS)!Nz T spectrumzT outputz* output = log( abs( spectrum ) ); _real_cepstrum_kernel)z -std=c++11)optionsctj|}tj|||}t|}tj|||jS)aa Calculates the real cepstrum of an input sequence x where the cepstrum is defined as the inverse Fourier transform of the log magnitude DFT (spectrum) of a signal. It's primarily used for source/speaker separation in speech signal processing Parameters ---------- x : ndarray Input sequence, if x is a matrix, return cepstrum in direction of axis n : int Size of Fourier Transform; If none, will use length of input array axis: int Direction for cepstrum calculation Returns ------- ceps : ndarray Complex cepstrum result naxis)cupyasarrayfftrifftreal)xrr spectrums u/home/jaya/work/projects/VOICE-AGENT/VIET/agent-env/lib/python3.11/site-packages/cupyx/signal/_acoustics/_cepstrum.py real_cepstrumr$sT* QAx||A|..H$X..H 8==QT= 2 2 77zC spectrum, raw T unwrappedzC output, T ndelayz ndelay = round( unwrapped[center] / M_PI ); const T temp { unwrapped[i] - ( M_PI * ndelay * i / center ) }; output = log( abs( spectrum ) ) + C( 0, temp ); _complex_cepstrum_kernelTzKconst int center { static_cast( 0.5 * ( _ind.size() + 1 ) ) };)r return_tuple loop_prepc6tj|}tj|||}tjtj|}t ||\}}tj|||j}||fS)a Calculates the complex cepstrum of a real valued input sequence x where the cepstrum is defined as the inverse Fourier transform of the log magnitude DFT (spectrum) of a signal. It's primarily used for source/speaker separation in speech signal processing. The input is altered to have zero-phase at pi radians (180 degrees) Parameters ---------- x : ndarray Input sequence, if x is a matrix, return cepstrum in direction of axis n : int Size of Fourier Transform; If none, will use length of input array axis: int Direction for cepstrum calculation Returns ------- ceps : ndarray Complex cepstrum result r)r r r unwrapanglerr r)rrr r unwrapped log_spectrumndelaycepss rcomplex_cepstrumrPs}* QAx||A|..H DJx0011I3HiHHL& 8=== 6 6 ;D <rz(C log_spectrum, int32 ndelay, float64 piz C spectrumz const double wrapped { log_spectrum.imag() + M_PI * ndelay * i / center }; spectrum = exp( C( log_spectrum.real(), wrapped ) ) _inverse_complex_cepstrum_kernelz2const double center { 0.5 * ( _ind.size() + 1 ) };)rrctj|}tj|}t||tj}tj|j}|S)a6Compute the inverse complex cepstrum of a real sequence. ceps : ndarray Real sequence to compute inverse complex cepstrum of. ndelay: int The amount of samples of circular delay added to `x`. Returns ------- x : ndarray The inverse complex cepstrum of the real sequence `ceps`. The inverse complex cepstrum is given by .. math:: x[n] = F^{-1}\left{\exp(F(c[n]))\right} where :math:`c_[n]` is the input signal and :math:`F` and :math:`F_{-1} are respectively the forward and backward Fourier transform. )r r r rpir r)rrrricepss rinverse_complex_cepstrumr#|sT <  D8<<%%L/ fdgNNH HMM( # # (E LrzT cepszT windowz if ( !i ) { window = ceps; } else if ( i < bend ) { window = ceps * 2.0; } else if ( i == bend ) { window = ceps * ( 1 - odd ); } else { window = 0; } _minimum_phase_kernelzconst bool odd { _ind.size() & 1 }; const int bend { static_cast( 0.5 * ( _ind.size() + odd ) ) };c|t|}t||}t|}tjtjtj|j}|S)aCompute the minimum phase reconstruction of a real sequence. x : ndarray Real sequence to compute the minimum phase reconstruction of. n : {None, int}, optional Length of the Fourier transform. Compute the minimum phase reconstruction of a real sequence using the real cepstrum. Returns ------- m : ndarray The minimum phase reconstruction of the real sequence `x`. N)r)lenrr$r r r expr)rrrwindowms r minimum_phaser*sf y FF a D "4 ( (F  dhtx||F334455:A Hr)Nr)N) r ElementwiseKernelrrrrrr#r$r*rrr-s1, /. 88886241!  ""<$:4#9. ' B $ $ $ ./.   0!,      r