comparison fft.c @ 10172:eda985c53dba libavcodec

ARM: 10l: fix large FFTs

10171:39a7bdbf0d40

{
    int i, j, m, n;
    float alpha, c1, s1, s2;
    int split_radix = 1;
    int av_unused has_vectors;
    int revtab_shift = 0;

    if (nbits < 2 || nbits > 16)
        goto fail;
    s->nbits = nbits;
    n = 1 << nbits;

    s->fft_permute = ff_fft_permute_neon;
    s->fft_calc    = ff_fft_calc_neon;
    s->imdct_calc  = ff_imdct_calc_neon;
    s->imdct_half  = ff_imdct_half_neon;
    s->mdct_calc   = ff_mdct_calc_neon;
    revtab_shift = 3;
#endif

    if (split_radix) {
        for(j=4; j<=nbits; j++) {
            int m = 1<<j;

                tab[i] = cos(i*freq);
            for(i=1; i<m/4; i++)
                tab[m/2-i] = tab[i];
        }
        for(i=0; i<n; i++)
            s->revtab[-split_radix_permutation(i, n, s->inverse) & (n-1)] =
                i << revtab_shift;
        s->tmp_buf = av_malloc(n * sizeof(FFTComplex));
    } else {
        int np, nblocks, np2, l;
        FFTComplex *q;

10172:eda985c53dba

{
    int i, j, m, n;
    float alpha, c1, s1, s2;
    int split_radix = 1;
    int av_unused has_vectors;

    if (nbits < 2 || nbits > 16)
        goto fail;
    s->nbits = nbits;
    n = 1 << nbits;

    s->fft_permute = ff_fft_permute_neon;
    s->fft_calc    = ff_fft_calc_neon;
    s->imdct_calc  = ff_imdct_calc_neon;
    s->imdct_half  = ff_imdct_half_neon;
    s->mdct_calc   = ff_mdct_calc_neon;
#endif

    if (split_radix) {
        for(j=4; j<=nbits; j++) {
            int m = 1<<j;

                tab[i] = cos(i*freq);
            for(i=1; i<m/4; i++)
                tab[m/2-i] = tab[i];
        }
        for(i=0; i<n; i++)
            s->revtab[-split_radix_permutation(i, n, s->inverse) & (n-1)] = i;
        s->tmp_buf = av_malloc(n * sizeof(FFTComplex));
    } else {
        int np, nblocks, np2, l;
        FFTComplex *q;