|
|
|
@@ -26,31 +26,6 @@ |
|
|
|
|
|
|
|
#include "dsputil_ppc.h" |
|
|
|
#include "util_altivec.h" |
|
|
|
/* |
|
|
|
those three macros are from libavcodec/fft.c |
|
|
|
and are required for the reference C code |
|
|
|
*/ |
|
|
|
/* butter fly op */ |
|
|
|
#define BF(pre, pim, qre, qim, pre1, pim1, qre1, qim1) \ |
|
|
|
{\ |
|
|
|
FFTSample ax, ay, bx, by;\ |
|
|
|
bx=pre1;\ |
|
|
|
by=pim1;\ |
|
|
|
ax=qre1;\ |
|
|
|
ay=qim1;\ |
|
|
|
pre = (bx + ax);\ |
|
|
|
pim = (by + ay);\ |
|
|
|
qre = (bx - ax);\ |
|
|
|
qim = (by - ay);\ |
|
|
|
} |
|
|
|
#define MUL16(a,b) ((a) * (b)) |
|
|
|
#define CMUL(pre, pim, are, aim, bre, bim) \ |
|
|
|
{\ |
|
|
|
pre = (MUL16(are, bre) - MUL16(aim, bim));\ |
|
|
|
pim = (MUL16(are, bim) + MUL16(bre, aim));\ |
|
|
|
} |
|
|
|
|
|
|
|
|
|
|
|
/** |
|
|
|
* Do a complex FFT with the parameters defined in ff_fft_init(). The |
|
|
|
* input data must be permuted before with s->revtab table. No |
|
|
|
|