xref: /third_party/ffmpeg/libavutil/tx_template.c

65     TXSample *tab = TX_TAB(ff_tx_tab_ ##len);                      \
68         *tab++ = RESCALE(cos(i*freq));                             \
70     *tab = 0;                                                      \
180     const TXSample *tab = TX_TAB(ff_tx_tab_53);
192     mtmp[0] = (int64_t)tab[0] * tmp[0].re;
193     mtmp[1] = (int64_t)tab[1] * tmp[0].im;
194     mtmp[2] = (int64_t)tab[2] * tmp[1].re;
195     mtmp[3] = (int64_t)tab[2] * tmp[1].im;
201     tmp[0].re = tab[0] * tmp[0].re;
202     tmp[0].im = tab[1] * tmp[0].im;
203     tmp[1].re = tab[2] * tmp[1].re;
204     tmp[1].im = tab[2] * tmp[1].im;
217     const TXSample *tab = TX_TAB(ff_tx_tab_53);                     \
227     SMUL(t[4].re, t[0].re, tab[4], tab[6], t[2].re, t[0].re);       \
228     SMUL(t[4].im, t[0].im, tab[4], tab[6], t[2].im, t[0].im);       \
229     CMUL(t[5].re, t[1].re, tab[5], tab[7], t[3].re, t[1].re);       \
230     CMUL(t[5].im, t[1].im, tab[5], tab[7], t[3].im, t[1].im);       \
256     const TXComplex *tab = (const TXComplex *)TX_TAB(ff_tx_tab_7);
272     mtmp[ 0] = ((int64_t)tab[0].re)*t[0].re - ((int64_t)tab[2].re)*t[4].re;
273     mtmp[ 1] = ((int64_t)tab[0].re)*t[4].re - ((int64_t)tab[1].re)*t[0].re;
274     mtmp[ 2] = ((int64_t)tab[0].re)*t[2].re - ((int64_t)tab[2].re)*t[0].re;
275     mtmp[ 3] = ((int64_t)tab[0].re)*t[0].im - ((int64_t)tab[1].re)*t[2].im;
276     mtmp[ 4] = ((int64_t)tab[0].re)*t[4].im - ((int64_t)tab[1].re)*t[0].im;
277     mtmp[ 5] = ((int64_t)tab[0].re)*t[2].im - ((int64_t)tab[2].re)*t[0].im;
279     mtmp[ 6] = ((int64_t)tab[2].im)*t[1].im + ((int64_t)tab[1].im)*t[5].im;
280     mtmp[ 7] = ((int64_t)tab[0].im)*t[5].im + ((int64_t)tab[2].im)*t[3].im;
281     mtmp[ 8] = ((int64_t)tab[2].im)*t[5].im + ((int64_t)tab[1].im)*t[3].im;
282     mtmp[ 9] = ((int64_t)tab[0].im)*t[1].re + ((int64_t)tab[1].im)*t[3].re;
283     mtmp[10] = ((int64_t)tab[2].im)*t[3].re + ((int64_t)tab[0].im)*t[5].re;
284     mtmp[11] = ((int64_t)tab[2].im)*t[1].re + ((int64_t)tab[1].im)*t[5].re;
286     z[0].re = (int32_t)(mtmp[ 0] - ((int64_t)tab[1].re)*t[2].re + 0x40000000 >> 31);
287     z[1].re = (int32_t)(mtmp[ 1] - ((int64_t)tab[2].re)*t[2].re + 0x40000000 >> 31);
288     z[2].re = (int32_t)(mtmp[ 2] - ((int64_t)tab[1].re)*t[4].re + 0x40000000 >> 31);
289     z[0].im = (int32_t)(mtmp[ 3] - ((int64_t)tab[2].re)*t[4].im + 0x40000000 >> 31);
290     z[1].im = (int32_t)(mtmp[ 4] - ((int64_t)tab[2].re)*t[2].im + 0x40000000 >> 31);
291     z[2].im = (int32_t)(mtmp[ 5] - ((int64_t)tab[1].re)*t[4].im + 0x40000000 >> 31);
293     t[0].re = (int32_t)(mtmp[ 6] - ((int64_t)tab[0].im)*t[3].im + 0x40000000 >> 31);
294     t[2].re = (int32_t)(mtmp[ 7] - ((int64_t)tab[1].im)*t[1].im + 0x40000000 >> 31);
295     t[4].re = (int32_t)(mtmp[ 8] + ((int64_t)tab[0].im)*t[1].im + 0x40000000 >> 31);
296     t[0].im = (int32_t)(mtmp[ 9] + ((int64_t)tab[2].im)*t[5].re + 0x40000000 >> 31);
297     t[2].im = (int32_t)(mtmp[10] - ((int64_t)tab[1].im)*t[1].re + 0x40000000 >> 31);
298     t[4].im = (int32_t)(mtmp[11] - ((int64_t)tab[0].im)*t[3].re + 0x40000000 >> 31);
300     z[0].re = tab[0].re*t[0].re - tab[2].re*t[4].re - tab[1].re*t[2].re;
301     z[1].re = tab[0].re*t[4].re - tab[1].re*t[0].re - tab[2].re*t[2].re;
302     z[2].re = tab[0].re*t[2].re - tab[2].re*t[0].re - tab[1].re*t[4].re;
303     z[0].im = tab[0].re*t[0].im - tab[1].re*t[2].im - tab[2].re*t[4].im;
304     z[1].im = tab[0].re*t[4].im - tab[1].re*t[0].im - tab[2].re*t[2].im;
305     z[2].im = tab[0].re*t[2].im - tab[2].re*t[0].im - tab[1].re*t[4].im;
310     t[0].re = tab[2].im*t[1].im + tab[1].im*t[5].im - tab[0].im*t[3].im;
311     t[2].re = tab[0].im*t[5].im + tab[2].im*t[3].im - tab[1].im*t[1].im;
312     t[4].re = tab[2].im*t[5].im + tab[1].im*t[3].im + tab[0].im*t[1].im;
313     t[0].im = tab[0].im*t[1].re + tab[1].im*t[3].re + tab[2].im*t[5].re;
314     t[2].im = tab[2].im*t[3].re + tab[0].im*t[5].re - tab[1].im*t[1].re;
315     t[4].im = tab[2].im*t[1].re + tab[1].im*t[5].re - tab[0].im*t[3].re;
342     const TXComplex *tab = (const TXComplex *)TX_TAB(ff_tx_tab_9);
379     y[3].re = (int32_t)(((int64_t)tab[0].im)*mtmp[0] + 0x40000000 >> 31);
380     y[3].im = (int32_t)(((int64_t)tab[0].im)*mtmp[1] + 0x40000000 >> 31);
382     mtmp[0] = (int32_t)(((int64_t)tab[0].re)*z[1].re + 0x40000000 >> 31);
383     mtmp[1] = (int32_t)(((int64_t)tab[0].re)*z[1].im + 0x40000000 >> 31);
384     mtmp[2] = (int32_t)(((int64_t)tab[0].re)*t[4].re + 0x40000000 >> 31);
385     mtmp[3] = (int32_t)(((int64_t)tab[0].re)*t[4].im + 0x40000000 >> 31);
392     mtmp[0] = ((int64_t)tab[1].re)*w[0].re;
393     mtmp[1] = ((int64_t)tab[1].re)*w[0].im;
394     mtmp[2] = ((int64_t)tab[2].im)*w[0].re;
395     mtmp[3] = ((int64_t)tab[2].im)*w[0].im;
396     mtmp[4] = ((int64_t)tab[1].im)*w[2].re;
397     mtmp[5] = ((int64_t)tab[1].im)*w[2].im;
398     mtmp[6] = ((int64_t)tab[2].re)*w[2].re;
399     mtmp[7] = ((int64_t)tab[2].re)*w[2].im;
401     x[1].re = (int32_t)(mtmp[0] + ((int64_t)tab[2].im)*w[1].re + 0x40000000 >> 31);
402     x[1].im = (int32_t)(mtmp[1] + ((int64_t)tab[2].im)*w[1].im + 0x40000000 >> 31);
403     x[2].re = (int32_t)(mtmp[2] - ((int64_t)tab[3].re)*w[1].re + 0x40000000 >> 31);
404     x[2].im = (int32_t)(mtmp[3] - ((int64_t)tab[3].re)*w[1].im + 0x40000000 >> 31);
405     y[1].re = (int32_t)(mtmp[4] + ((int64_t)tab[2].re)*w[3].re + 0x40000000 >> 31);
406     y[1].im = (int32_t)(mtmp[5] + ((int64_t)tab[2].re)*w[3].im + 0x40000000 >> 31);
407     y[2].re = (int32_t)(mtmp[6] - ((int64_t)tab[3].im)*w[3].re + 0x40000000 >> 31);
408     y[2].im = (int32_t)(mtmp[7] - ((int64_t)tab[3].im)*w[3].im + 0x40000000 >> 31);
410     y[0].re = (int32_t)(((int64_t)tab[0].im)*t[5].re + 0x40000000 >> 31);
411     y[0].im = (int32_t)(((int64_t)tab[0].im)*t[5].im + 0x40000000 >> 31);
414     y[3].re = tab[0].im*(t[1].re - t[3].re + t[7].re);
415     y[3].im = tab[0].im*(t[1].im - t[3].im + t[7].im);
417     x[3].re = z[0].re  + tab[0].re*z[1].re;
418     x[3].im = z[0].im  + tab[0].re*z[1].im;
419     z[0].re = in[0].re + tab[0].re*t[4].re;
420     z[0].im = in[0].im + tab[0].re*t[4].im;
422     x[1].re = tab[1].re*w[0].re + tab[2].im*w[1].re;
423     x[1].im = tab[1].re*w[0].im + tab[2].im*w[1].im;
424     x[2].re = tab[2].im*w[0].re - tab[3].re*w[1].re;
425     x[2].im = tab[2].im*w[0].im - tab[3].re*w[1].im;
426     y[1].re = tab[1].im*w[2].re + tab[2].re*w[3].re;
427     y[1].im = tab[1].im*w[2].im + tab[2].re*w[3].im;
428     y[2].re = tab[2].re*w[2].re - tab[3].im*w[3].re;
429     y[2].im = tab[2].re*w[2].im - tab[3].im*w[3].im;
431     y[0].re = tab[0].im*t[5].re;
432     y[0].im = tab[0].im*t[5].im;
1289     TXSample *tab;
1300     tab = (TXSample *)s->exp;
1306     *tab++ = RESCALE((inv ? 0.5 : 1.0) * m);
1307     *tab++ = RESCALE(inv ? 0.5*m : 1.0);
1308     *tab++ = RESCALE( m);
1309     *tab++ = RESCALE(-m);
1311     *tab++ = RESCALE( (0.5 - 0.0) * m);
1312     *tab++ = RESCALE( (0.0 - 0.5) * m);
1313     *tab++ = RESCALE( (0.5 - inv) * m);
1314     *tab++ = RESCALE(-(0.5 - inv) * m);
1317         *tab++ = RESCALE(cos(i*f));
1319         *tab++ = RESCALE(cos(i*f) * (inv ? +1.0 : -1.0));

Indexes created Thu Nov 07 10:32:03 CST 2024