Lines Matching refs:ncfft
84 int k,ncfft;
91 ncfft = st->substate->nfft;
111 freqdata[ncfft].r = tdc.r - tdc.i;
113 freqdata[ncfft].i = freqdata[0].i = _mm_set1_ps(0);
115 freqdata[ncfft].i = freqdata[0].i = 0;
118 for ( k=1;k <= ncfft/2 ; ++k ) {
120 fpnk.r = st->tmpbuf[ncfft-k].r;
121 fpnk.i = - st->tmpbuf[ncfft-k].i;
131 freqdata[ncfft-k].r = HALF_OF(f1k.r - tw.r);
132 freqdata[ncfft-k].i = HALF_OF(tw.i - f1k.i);
139 int k, ncfft;
145 ncfft = st->substate->nfft;
147 st->tmpbuf[0].r = freqdata[0].r + freqdata[ncfft].r;
148 st->tmpbuf[0].i = freqdata[0].r - freqdata[ncfft].r;
151 for (k = 1; k <= ncfft / 2; ++k) {
154 fnkc.r = freqdata[ncfft - k].r;
155 fnkc.i = -freqdata[ncfft - k].i;
163 C_SUB (st->tmpbuf[ncfft - k], fek, fok);
165 st->tmpbuf[ncfft - k].i *= _mm_set1_ps(-1.0);
167 st->tmpbuf[ncfft - k].i *= -1;
176 int k,ncfft;
184 ncfft = st->substate->nfft;
204 freqdata[2*ncfft-1] = tdc.r - tdc.i;
206 for ( k=1;k <= ncfft/2 ; ++k )
209 fpnk.r = st->tmpbuf[ncfft-k].r;
210 fpnk.i = - st->tmpbuf[ncfft-k].i;
221 freqdata[2*(ncfft-k)-1] = HALF_OF(f1k.r - tw.r);
222 freqdata[2*(ncfft-k)] = HALF_OF(tw.i - f1k.i);
225 /*f1k.r = PSHR32(ADD32(EXTEND32(st->tmpbuf[k].r), EXTEND32(st->tmpbuf[ncfft-k].r)),1);
226 f1k.i = PSHR32(SUB32(EXTEND32(st->tmpbuf[k].i), EXTEND32(st->tmpbuf[ncfft-k].i)),1);
227 f2k.r = PSHR32(SUB32(EXTEND32(st->tmpbuf[k].r), EXTEND32(st->tmpbuf[ncfft-k].r)),1);
228 f2k.i = SHR32(ADD32(EXTEND32(st->tmpbuf[k].i), EXTEND32(st->tmpbuf[ncfft-k].i)),1);
234 freqdata[2*(ncfft-k)-1] = HALF_OF(f1k.r - tw.r);
235 freqdata[2*(ncfft-k)] = HALF_OF(tw.i - f1k.i);
237 f2k.r = SHR32(SUB32(EXTEND32(st->tmpbuf[k].r), EXTEND32(st->tmpbuf[ncfft-k].r)),1);
238 f2k.i = PSHR32(ADD32(EXTEND32(st->tmpbuf[k].i), EXTEND32(st->tmpbuf[ncfft-k].i)),1);
240 f1kr = SHL32(ADD32(EXTEND32(st->tmpbuf[k].r), EXTEND32(st->tmpbuf[ncfft-k].r)),13);
241 f1ki = SHL32(SUB32(EXTEND32(st->tmpbuf[k].i), EXTEND32(st->tmpbuf[ncfft-k].i)),13);
249 freqdata[2*(ncfft-k)-1] = PSHR32(f1kr - twr, 15);
250 freqdata[2*(ncfft-k)] = PSHR32(twi - f1ki, 15);
254 freqdata[2*(ncfft-k)-1] = .5f*(f1kr - twr);
255 freqdata[2*(ncfft-k)] = .5f*(twi - f1ki);
264 int k, ncfft;
270 ncfft = st->substate->nfft;
272 st->tmpbuf[0].r = freqdata[0] + freqdata[2*ncfft-1];
273 st->tmpbuf[0].i = freqdata[0] - freqdata[2*ncfft-1];
276 for (k = 1; k <= ncfft / 2; ++k) {
280 fnkc.r = freqdata[2*(ncfft - k)-1];
281 fnkc.i = -freqdata[2*(ncfft - k)];
289 C_SUB (st->tmpbuf[ncfft - k], fek, fok);
291 st->tmpbuf[ncfft - k].i *= _mm_set1_ps(-1.0);
293 st->tmpbuf[ncfft - k].i *= -1;
