84     int k,ncfft;  in kiss_fftr()  local
118     for ( k=1;k <= ncfft/2 ; ++k ) {  in kiss_fftr()
119         fpk    = st->tmpbuf[k];   in kiss_fftr()
120         fpnk.r =   st->tmpbuf[ncfft-k].r;  in kiss_fftr()
121         fpnk.i = - st->tmpbuf[ncfft-k].i;  in kiss_fftr()
127         C_MUL( tw , f2k , st->super_twiddles[k]);  in kiss_fftr()
129         freqdata[k].r = HALF_OF(f1k.r + tw.r);  in kiss_fftr()
130         freqdata[k].i = HALF_OF(f1k.i + tw.i);  in kiss_fftr()
131         freqdata[ncfft-k].r = HALF_OF(f1k.r - tw.r);  in kiss_fftr()
132         freqdata[ncfft-k].i = HALF_OF(tw.i - f1k.i);  in kiss_fftr()
139     int k, ncfft;  in kiss_fftri()  local
151     for (k = 1; k <= ncfft / 2; ++k) {  in kiss_fftri()
153         fk = freqdata[k];  in kiss_fftri()
154         fnkc.r = freqdata[ncfft - k].r;  in kiss_fftri()
155         fnkc.i = -freqdata[ncfft - k].i;  in kiss_fftri()
161         C_MUL (fok, tmp, st->super_twiddles[k]);  in kiss_fftri()
162         C_ADD (st->tmpbuf[k],     fek, fok);  in kiss_fftri()
163         C_SUB (st->tmpbuf[ncfft - k], fek, fok);  in kiss_fftri()
165         st->tmpbuf[ncfft - k].i *= _mm_set1_ps(-1.0);  in kiss_fftri()
167         st->tmpbuf[ncfft - k].i *= -1;  in kiss_fftri()
176    int k,ncfft;  in kiss_fftr2()  local
206    for ( k=1;k <= ncfft/2 ; ++k )  in kiss_fftr2()
237       f2k.r = SHR32(SUB32(EXTEND32(st->tmpbuf[k].r), EXTEND32(st->tmpbuf[ncfft-k].r)),1);  in kiss_fftr2()
238       f2k.i = PSHR32(ADD32(EXTEND32(st->tmpbuf[k].i), EXTEND32(st->tmpbuf[ncfft-k].i)),1);  in kiss_fftr2()
240       f1kr = SHL32(ADD32(EXTEND32(st->tmpbuf[k].r), EXTEND32(st->tmpbuf[ncfft-k].r)),13);  in kiss_fftr2()
241       f1ki = SHL32(SUB32(EXTEND32(st->tmpbuf[k].i), EXTEND32(st->tmpbuf[ncfft-k].i)),13);  in kiss_fftr2()
243 …twr = SHR32(SUB32(MULT16_16(f2k.r,st->super_twiddles[k].r),MULT16_16(f2k.i,st->super_twiddles[k].i…  in kiss_fftr2()
244 …twi = SHR32(ADD32(MULT16_16(f2k.i,st->super_twiddles[k].r),MULT16_16(f2k.r,st->super_twiddles[k].i…  in kiss_fftr2()
247       freqdata[2*k-1] = PSHR32(f1kr + twr, 15);  in kiss_fftr2()
248       freqdata[2*k] = PSHR32(f1ki + twi, 15);  in kiss_fftr2()
249       freqdata[2*(ncfft-k)-1] = PSHR32(f1kr - twr, 15);  in kiss_fftr2()
250       freqdata[2*(ncfft-k)] = PSHR32(twi - f1ki, 15);  in kiss_fftr2()
252       freqdata[2*k-1] = .5f*(f1kr + twr);  in kiss_fftr2()
253       freqdata[2*k] = .5f*(f1ki + twi);  in kiss_fftr2()
254       freqdata[2*(ncfft-k)-1] = .5f*(f1kr - twr);  in kiss_fftr2()
255       freqdata[2*(ncfft-k)] = .5f*(twi - f1ki);  in kiss_fftr2()
264    int k, ncfft;  in kiss_fftri2()  local
276    for (k = 1; k <= ncfft / 2; ++k) {  in kiss_fftri2()
278       fk.r = freqdata[2*k-1];  in kiss_fftri2()
279       fk.i = freqdata[2*k];  in kiss_fftri2()
280       fnkc.r = freqdata[2*(ncfft - k)-1];  in kiss_fftri2()
281       fnkc.i = -freqdata[2*(ncfft - k)];  in kiss_fftri2()
287       C_MUL (fok, tmp, st->super_twiddles[k]);  in kiss_fftri2()
288       C_ADD (st->tmpbuf[k],     fek, fok);  in kiss_fftri2()
289       C_SUB (st->tmpbuf[ncfft - k], fek, fok);  in kiss_fftri2()
291       st->tmpbuf[ncfft - k].i *= _mm_set1_ps(-1.0);  in kiss_fftri2()
293       st->tmpbuf[ncfft - k].i *= -1;  in kiss_fftri2()