mirror/kissfft
1
0
Fork 0
mirror of https://github.com/mborgerding/kissfft.git synced 2025-06-04 01:28:23 -04:00
Code Issues Packages Projects Releases Wiki Activity

Normalized identation in new code.

Browse Source
This commit is contained in:
Ralph Tandetzky 2016-04-25 10:17:41 +02:00
parent 1fbc2b6ab4
commit 6e0d8bbcd2
1 changed files with 14 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
kissfft.hh
View File

@ -96,7 +96,7 @@ class kissfft
{ {
const std::size_t N = _nfft; const std::size_t N = _nfft;
if ( N == 0 ) if ( N == 0 )
return; return;
// perform complex FFT // perform complex FFT
transform( reinterpret_cast<const cpx_type*>(src), dst ); transform( reinterpret_cast<const cpx_type*>(src), dst );
@ -111,21 +111,21 @@ class kissfft
const cpx_type twiddle_mul = exp( cpx_type(0, half_phi_inc) ); const cpx_type twiddle_mul = exp( cpx_type(0, half_phi_inc) );
for ( std::size_t k = 1; 2*k < N; ++k ) for ( std::size_t k = 1; 2*k < N; ++k )
{ {
const cpx_type w = 0.5 * cpx_type( const cpx_type w = 0.5 * cpx_type(
dst[k].real() + dst[N-k].real(), dst[k].real() + dst[N-k].real(),
dst[k].imag() - dst[N-k].imag() ); dst[k].imag() - dst[N-k].imag() );
const cpx_type z = 0.5 * cpx_type( const cpx_type z = 0.5 * cpx_type(
dst[k].imag() + dst[N-k].imag(), dst[k].imag() + dst[N-k].imag(),
-dst[k].real() + dst[N-k].real() ); -dst[k].real() + dst[N-k].real() );
const cpx_type twiddle = const cpx_type twiddle =
k % 2 == 0 ? k % 2 == 0 ?
_twiddles[k/2] : _twiddles[k/2] :
_twiddles[k/2] * twiddle_mul; _twiddles[k/2] * twiddle_mul;
dst[ k] = w + twiddle * z; dst[ k] = w + twiddle * z;
dst[N-k] = conj( w - twiddle * z ); dst[N-k] = conj( w - twiddle * z );
} }
if ( N % 2 == 0 ) if ( N % 2 == 0 )
dst[N/2] = conj( dst[N/2] ); dst[N/2] = conj( dst[N/2] );
} }
private: private: