81             int m = _stageRemainder[stage];  in kf_work()  local
83             cpx_type * Fout_end = Fout + p*m;  in kf_work()
85             if (m==1) {  in kf_work()
98                 }while( (Fout += m) != Fout_end );  in kf_work()
105                 case 2: kf_bfly2(Fout,fstride,m); break;  in kf_work()
106                 case 3: kf_bfly3(Fout,fstride,m); break;  in kf_work()
107                 case 4: kf_bfly4(Fout,fstride,m); break;  in kf_work()
108                 case 5: kf_bfly5(Fout,fstride,m); break;  in kf_work()
109                 default: kf_bfly_generic(Fout,fstride,m,p); break;  in kf_work()
123         void kf_bfly2( cpx_type * Fout, const size_t fstride, int m)  in kf_bfly2()  argument
125             for (int k=0;k<m;++k) {  in kf_bfly2()
126                 cpx_type t = Fout[m+k] * _traits.twiddle(k*fstride);  in kf_bfly2()
127                 Fout[m+k] = Fout[k] - t;  in kf_bfly2()
132         void kf_bfly4( cpx_type * Fout, const size_t fstride, const size_t m)  in kf_bfly4()  argument
136             for (size_t k=0;k<m;++k) {  in kf_bfly4()
137                 scratch[0] = Fout[k+m] * _traits.twiddle(k*fstride);  in kf_bfly4()
138                 scratch[1] = Fout[k+2*m] * _traits.twiddle(k*fstride*2);  in kf_bfly4()
139                 scratch[2] = Fout[k+3*m] * _traits.twiddle(k*fstride*3);  in kf_bfly4()
147                 Fout[k+2*m]  = Fout[k] - scratch[3];  in kf_bfly4()
149                 Fout[k+m] = scratch[5] + scratch[4];  in kf_bfly4()
150                 Fout[k+3*m] = scratch[5] - scratch[4];  in kf_bfly4()
154         void kf_bfly3( cpx_type * Fout, const size_t fstride, const size_t m)  in kf_bfly3()  argument
156             size_t k=m;  in kf_bfly3()
157             const size_t m2 = 2*m;  in kf_bfly3()
161             epi3 = _twiddles[fstride*m];  in kf_bfly3()
166                 C_FIXDIV(*Fout,3); C_FIXDIV(Fout[m],3); C_FIXDIV(Fout[m2],3);  in kf_bfly3()
168                 C_MUL(scratch[1],Fout[m] , *tw1);  in kf_bfly3()
176 …Fout[m] = cpx_type( Fout->real() - HALF_OF(scratch[3].real() ) , Fout->imag() - HALF_OF(scratch[3]…  in kf_bfly3()
182 …   Fout[m2] = cpx_type(  Fout[m].real() + scratch[0].imag() , Fout[m].imag() - scratch[0].real() );  in kf_bfly3()
184                 C_ADDTO( Fout[m] , cpx_type( -scratch[0].imag(),scratch[0].real() ) );  in kf_bfly3()
189         void kf_bfly5( cpx_type * Fout, const size_t fstride, const size_t m)  in kf_bfly5()  argument
197             ya = twiddles[fstride*m];  in kf_bfly5()
198             yb = twiddles[fstride*2*m];  in kf_bfly5()
201             Fout1=Fout0+m;  in kf_bfly5()
202             Fout2=Fout0+2*m;  in kf_bfly5()
203             Fout3=Fout0+3*m;  in kf_bfly5()
204             Fout4=Fout0+4*m;  in kf_bfly5()
207             for ( u=0; u<m; ++u ) {  in kf_bfly5()
259                 int m,  in kf_bfly_generic()  argument
269             for ( u=0; u<m; ++u ) {  in kf_bfly_generic()
274                     k += m;  in kf_bfly_generic()
287                     k += m;  in kf_bfly_generic()