void fft_calc_altivec(FFTContext *s, FFTComplex *z)

{

POWERPC_TBL_DECLARE(altivec_fft_num, s->nbits >= 6);

#ifdef ALTIVEC_USE_REFERENCE_C_CODE

    int ln = s->nbits;

    int	j, np, np2;

    int	nblocks, nloops;

    register FFTComplex *p, *q;

    FFTComplex *exptab = s->exptab;

    int l;

    FFTSample tmp_re, tmp_im;

    

POWERPC_TBL_START_COUNT(altivec_fft_num, s->nbits >= 6);

 

    np = 1 << ln;



    /* pass 0 */



    p=&z[0];

    j=(np >> 1);

    do {

        BF(p[0].re, p[0].im, p[1].re, p[1].im, 

           p[0].re, p[0].im, p[1].re, p[1].im);

        p+=2;

    } while (--j != 0);



    /* pass 1 */



    

    p=&z[0];

    j=np >> 2;

    if (s->inverse) {

        do {

            BF(p[0].re, p[0].im, p[2].re, p[2].im, 

               p[0].re, p[0].im, p[2].re, p[2].im);

            BF(p[1].re, p[1].im, p[3].re, p[3].im, 

               p[1].re, p[1].im, -p[3].im, p[3].re);

            p+=4;

        } while (--j != 0);

    } else {

        do {

            BF(p[0].re, p[0].im, p[2].re, p[2].im, 

               p[0].re, p[0].im, p[2].re, p[2].im);

            BF(p[1].re, p[1].im, p[3].re, p[3].im, 

               p[1].re, p[1].im, p[3].im, -p[3].re);

            p+=4;

        } while (--j != 0);

    }

    /* pass 2 .. ln-1 */



    nblocks = np >> 3;

    nloops = 1 << 2;

    np2 = np >> 1;

    do {

        p = z;

        q = z + nloops;

        for (j = 0; j < nblocks; ++j) {

            BF(p->re, p->im, q->re, q->im,

               p->re, p->im, q->re, q->im);

            

            p++;

            q++;

            for(l = nblocks; l < np2; l += nblocks) {

                CMUL(tmp_re, tmp_im, exptab[l].re, exptab[l].im, q->re, q->im);

                BF(p->re, p->im, q->re, q->im,

                   p->re, p->im, tmp_re, tmp_im);

                p++;

                q++;

            }



            p += nloops;

            q += nloops;

        }

        nblocks = nblocks >> 1;

        nloops = nloops << 1;

    } while (nblocks != 0);



POWERPC_TBL_STOP_COUNT(altivec_fft_num, s->nbits >= 6);



#else /* ALTIVEC_USE_REFERENCE_C_CODE */

#ifdef CONFIG_DARWIN

    register const vector float vczero = (const vector float)(0.);

#else

    register const vector float vczero = (const vector float){0.,0.,0.,0.};

#endif

    

    int ln = s->nbits;

    int	j, np, np2;

    int	nblocks, nloops;

    register FFTComplex *p, *q;

    FFTComplex *cptr, *cptr1;

    int k;



POWERPC_TBL_START_COUNT(altivec_fft_num, s->nbits >= 6);



    np = 1 << ln;



    {

        vector float *r, a, b, a1, c1, c2;



        r = (vector float *)&z[0];



        c1 = vcii(p,p,n,n);

        

        if (s->inverse)

            {

                c2 = vcii(p,p,n,p);

            }

        else

            {

                c2 = vcii(p,p,p,n);

            }

        

        j = (np >> 2);

        do {

            a = vec_ld(0, r);

            a1 = vec_ld(sizeof(vector float), r);

            

            b = vec_perm(a,a,vcprmle(1,0,3,2));

            a = vec_madd(a,c1,b);

            /* do the pass 0 butterfly */

            

            b = vec_perm(a1,a1,vcprmle(1,0,3,2));

            b = vec_madd(a1,c1,b);

            /* do the pass 0 butterfly */

            

            /* multiply third by -i */

            b = vec_perm(b,b,vcprmle(2,3,1,0));

            

            /* do the pass 1 butterfly */

            vec_st(vec_madd(b,c2,a), 0, r);

            vec_st(vec_nmsub(b,c2,a), sizeof(vector float), r);

            

            r += 2;

        } while (--j != 0);

    }

    /* pass 2 .. ln-1 */



    nblocks = np >> 3;

    nloops = 1 << 2;

    np2 = np >> 1;



    cptr1 = s->exptab1;

    do {

        p = z;

        q = z + nloops;

        j = nblocks;

        do {

            cptr = cptr1;

            k = nloops >> 1;

            do {

                vector float a,b,c,t1;



                a = vec_ld(0, (float*)p);

                b = vec_ld(0, (float*)q);

                

                /* complex mul */

                c = vec_ld(0, (float*)cptr);

                /*  cre*re cim*re */

                t1 = vec_madd(c, vec_perm(b,b,vcprmle(2,2,0,0)),vczero);

                c = vec_ld(sizeof(vector float), (float*)cptr);

                /*  -cim*im cre*im */

                b = vec_madd(c, vec_perm(b,b,vcprmle(3,3,1,1)),t1);

                

                /* butterfly */

                vec_st(vec_add(a,b), 0, (float*)p);

                vec_st(vec_sub(a,b), 0, (float*)q);

                

                p += 2;

                q += 2;

                cptr += 4;

            } while (--k);

            

            p += nloops;

            q += nloops;

        } while (--j);

        cptr1 += nloops * 2;

        nblocks = nblocks >> 1;

        nloops = nloops << 1;

    } while (nblocks != 0);



POWERPC_TBL_STOP_COUNT(altivec_fft_num, s->nbits >= 6);



#endif /* ALTIVEC_USE_REFERENCE_C_CODE */

}

void fft_calc_altivec(FFTContext *s, FFTComplex *z)

{

POWERPC_TBL_DECLARE(altivec_fft_num, s->nbits >= 6);

#ifdef ALTIVEC_USE_REFERENCE_C_CODE

    int ln = s->nbits;

    int	j, np, np2;

    int	nblocks, nloops;

    register FFTComplex *p, *q;

    FFTComplex *exptab = s->exptab;

    int l;

    FFTSample tmp_re, tmp_im;

    

POWERPC_TBL_START_COUNT(altivec_fft_num, s->nbits >= 6);

 

    np = 1 << ln;



    



    p=&z[0];

    j=(np >> 1);

    do {

        BF(p[0].re, p[0].im, p[1].re, p[1].im, 

           p[0].re, p[0].im, p[1].re, p[1].im);

        p+=2;

    } while (--j != 0);



    



    

    p=&z[0];

    j=np >> 2;

    if (s->inverse) {

        do {

            BF(p[0].re, p[0].im, p[2].re, p[2].im, 

               p[0].re, p[0].im, p[2].re, p[2].im);

            BF(p[1].re, p[1].im, p[3].re, p[3].im, 

               p[1].re, p[1].im, -p[3].im, p[3].re);

            p+=4;

        } while (--j != 0);

    } else {

        do {

            BF(p[0].re, p[0].im, p[2].re, p[2].im, 

               p[0].re, p[0].im, p[2].re, p[2].im);

            BF(p[1].re, p[1].im, p[3].re, p[3].im, 

               p[1].re, p[1].im, p[3].im, -p[3].re);

            p+=4;

        } while (--j != 0);

    }

    



    nblocks = np >> 3;

    nloops = 1 << 2;

    np2 = np >> 1;

    do {

        p = z;

        q = z + nloops;

        for (j = 0; j < nblocks; ++j) {

            BF(p->re, p->im, q->re, q->im,

               p->re, p->im, q->re, q->im);

            

            p++;

            q++;

            for(l = nblocks; l < np2; l += nblocks) {

                CMUL(tmp_re, tmp_im, exptab[l].re, exptab[l].im, q->re, q->im);

                BF(p->re, p->im, q->re, q->im,

                   p->re, p->im, tmp_re, tmp_im);

                p++;

                q++;

            }



            p += nloops;

            q += nloops;

        }

        nblocks = nblocks >> 1;

        nloops = nloops << 1;

    } while (nblocks != 0);



POWERPC_TBL_STOP_COUNT(altivec_fft_num, s->nbits >= 6);



#else 

#ifdef CONFIG_DARWIN

    register const vector float vczero = (const vector float)(0.);

#else

    register const vector float vczero = (const vector float){0.,0.,0.,0.};

#endif

    

    int ln = s->nbits;

    int	j, np, np2;

    int	nblocks, nloops;

    register FFTComplex *p, *q;

    FFTComplex *cptr, *cptr1;

    int k;



POWERPC_TBL_START_COUNT(altivec_fft_num, s->nbits >= 6);



    np = 1 << ln;



    {

        vector float *r, a, b, a1, c1, c2;



        r = (vector float *)&z[0];



        c1 = vcii(p,p,n,n);

        

        if (s->inverse)

            {

                c2 = vcii(p,p,n,p);

            }

        else

            {

                c2 = vcii(p,p,p,n);

            }

        

        j = (np >> 2);

        do {

            a = vec_ld(0, r);

            a1 = vec_ld(sizeof(vector float), r);

            

            b = vec_perm(a,a,vcprmle(1,0,3,2));

            a = vec_madd(a,c1,b);

            

            

            b = vec_perm(a1,a1,vcprmle(1,0,3,2));

            b = vec_madd(a1,c1,b);

            

            

            

            b = vec_perm(b,b,vcprmle(2,3,1,0));

            

            

            vec_st(vec_madd(b,c2,a), 0, r);

            vec_st(vec_nmsub(b,c2,a), sizeof(vector float), r);

            

            r += 2;

        } while (--j != 0);

    }

    



    nblocks = np >> 3;

    nloops = 1 << 2;

    np2 = np >> 1;



    cptr1 = s->exptab1;

    do {

        p = z;

        q = z + nloops;

        j = nblocks;

        do {

            cptr = cptr1;

            k = nloops >> 1;

            do {

                vector float a,b,c,t1;



                a = vec_ld(0, (float*)p);

                b = vec_ld(0, (float*)q);

                

                

                c = vec_ld(0, (float*)cptr);

                

                t1 = vec_madd(c, vec_perm(b,b,vcprmle(2,2,0,0)),vczero);

                c = vec_ld(sizeof(vector float), (float*)cptr);

                

                b = vec_madd(c, vec_perm(b,b,vcprmle(3,3,1,1)),t1);

                

                

                vec_st(vec_add(a,b), 0, (float*)p);

                vec_st(vec_sub(a,b), 0, (float*)q);

                

                p += 2;

                q += 2;

                cptr += 4;

            } while (--k);

            

            p += nloops;

            q += nloops;

        } while (--j);

        cptr1 += nloops * 2;

        nblocks = nblocks >> 1;

        nloops = nloops << 1;

    } while (nblocks != 0);



POWERPC_TBL_STOP_COUNT(altivec_fft_num, s->nbits >= 6);



#endif 

}

void FUNC_0(FFTContext *VAR_0, FFTComplex *VAR_1)

{

POWERPC_TBL_DECLARE(altivec_fft_num, VAR_0->nbits >= 6);

#ifdef ALTIVEC_USE_REFERENCE_C_CODE

    int VAR_3 = VAR_0->nbits;

    int	VAR_4, VAR_5, VAR_6;

    int	VAR_7, VAR_8;

    register FFTComplex *VAR_9, *q;

    FFTComplex *exptab = VAR_0->exptab;

    int l;

    FFTSample tmp_re, tmp_im;

    

POWERPC_TBL_START_COUNT(altivec_fft_num, VAR_0->nbits >= 6);

 

    VAR_5 = 1 << VAR_3;



    



    VAR_9=&VAR_1[0];

    VAR_4=(VAR_5 >> 1);

    do {

        BF(VAR_9[0].re, VAR_9[0].im, VAR_9[1].re, VAR_9[1].im, 

           VAR_9[0].re, VAR_9[0].im, VAR_9[1].re, VAR_9[1].im);

        VAR_9+=2;

    } while (--VAR_4 != 0);



    



    

    VAR_9=&VAR_1[0];

    VAR_4=VAR_5 >> 2;

    if (VAR_0->inverse) {

        do {

            BF(VAR_9[0].re, VAR_9[0].im, VAR_9[2].re, VAR_9[2].im, 

               VAR_9[0].re, VAR_9[0].im, VAR_9[2].re, VAR_9[2].im);

            BF(VAR_9[1].re, VAR_9[1].im, VAR_9[3].re, VAR_9[3].im, 

               VAR_9[1].re, VAR_9[1].im, -VAR_9[3].im, VAR_9[3].re);

            VAR_9+=4;

        } while (--VAR_4 != 0);

    } else {

        do {

            BF(VAR_9[0].re, VAR_9[0].im, VAR_9[2].re, VAR_9[2].im, 

               VAR_9[0].re, VAR_9[0].im, VAR_9[2].re, VAR_9[2].im);

            BF(VAR_9[1].re, VAR_9[1].im, VAR_9[3].re, VAR_9[3].im, 

               VAR_9[1].re, VAR_9[1].im, VAR_9[3].im, -VAR_9[3].re);

            VAR_9+=4;

        } while (--VAR_4 != 0);

    }

    



    VAR_7 = VAR_5 >> 3;

    VAR_8 = 1 << 2;

    VAR_6 = VAR_5 >> 1;

    do {

        VAR_9 = VAR_1;

        q = VAR_1 + VAR_8;

        for (VAR_4 = 0; VAR_4 < VAR_7; ++VAR_4) {

            BF(VAR_9->re, VAR_9->im, q->re, q->im,

               VAR_9->re, VAR_9->im, q->re, q->im);

            

            VAR_9++;

            q++;

            for(l = VAR_7; l < VAR_6; l += VAR_7) {

                CMUL(tmp_re, tmp_im, exptab[l].re, exptab[l].im, q->re, q->im);

                BF(VAR_9->re, VAR_9->im, q->re, q->im,

                   VAR_9->re, VAR_9->im, tmp_re, tmp_im);

                VAR_9++;

                q++;

            }



            VAR_9 += VAR_8;

            q += VAR_8;

        }

        VAR_7 = VAR_7 >> 1;

        VAR_8 = VAR_8 << 1;

    } while (VAR_7 != 0);



POWERPC_TBL_STOP_COUNT(altivec_fft_num, VAR_0->nbits >= 6);



#else 

#ifdef CONFIG_DARWIN

    register const vector float VAR_2 = (const vector float)(0.);

#else

    register const vector float VAR_2 = (const vector float){0.,0.,0.,0.};

#endif

    

    int VAR_3 = VAR_0->nbits;

    int	VAR_4, VAR_5, VAR_6;

    int	VAR_7, VAR_8;

    register FFTComplex *VAR_9, *q;

    FFTComplex *cptr, *cptr1;

    int VAR_10;



POWERPC_TBL_START_COUNT(altivec_fft_num, VAR_0->nbits >= 6);



    VAR_5 = 1 << VAR_3;



    {

        vector float *r, a, b, a1, c1, c2;



        r = (vector float *)&VAR_1[0];



        c1 = vcii(VAR_9,VAR_9,n,n);

        

        if (VAR_0->inverse)

            {

                c2 = vcii(VAR_9,VAR_9,n,VAR_9);

            }

        else

            {

                c2 = vcii(VAR_9,VAR_9,VAR_9,n);

            }

        

        VAR_4 = (VAR_5 >> 2);

        do {

            a = vec_ld(0, r);

            a1 = vec_ld(sizeof(vector float), r);

            

            b = vec_perm(a,a,vcprmle(1,0,3,2));

            a = vec_madd(a,c1,b);

            

            

            b = vec_perm(a1,a1,vcprmle(1,0,3,2));

            b = vec_madd(a1,c1,b);

            

            

            

            b = vec_perm(b,b,vcprmle(2,3,1,0));

            

            

            vec_st(vec_madd(b,c2,a), 0, r);

            vec_st(vec_nmsub(b,c2,a), sizeof(vector float), r);

            

            r += 2;

        } while (--VAR_4 != 0);

    }

    



    VAR_7 = VAR_5 >> 3;

    VAR_8 = 1 << 2;

    VAR_6 = VAR_5 >> 1;



    cptr1 = VAR_0->exptab1;

    do {

        VAR_9 = VAR_1;

        q = VAR_1 + VAR_8;

        VAR_4 = VAR_7;

        do {

            cptr = cptr1;

            VAR_10 = VAR_8 >> 1;

            do {

                vector float a,b,c,t1;



                a = vec_ld(0, (float*)VAR_9);

                b = vec_ld(0, (float*)q);

                

                

                c = vec_ld(0, (float*)cptr);

                

                t1 = vec_madd(c, vec_perm(b,b,vcprmle(2,2,0,0)),VAR_2);

                c = vec_ld(sizeof(vector float), (float*)cptr);

                

                b = vec_madd(c, vec_perm(b,b,vcprmle(3,3,1,1)),t1);

                

                

                vec_st(vec_add(a,b), 0, (float*)VAR_9);

                vec_st(vec_sub(a,b), 0, (float*)q);

                

                VAR_9 += 2;

                q += 2;

                cptr += 4;

            } while (--VAR_10);

            

            VAR_9 += VAR_8;

            q += VAR_8;

        } while (--VAR_4);

        cptr1 += VAR_8 * 2;

        VAR_7 = VAR_7 >> 1;

        VAR_8 = VAR_8 << 1;

    } while (VAR_7 != 0);



POWERPC_TBL_STOP_COUNT(altivec_fft_num, VAR_0->nbits >= 6);



#endif 

}

void FUNC_0(FFTContext *VAR_0, FFTComplex *VAR_1)
{ POWERPC_TBL_DECLARE(altivec_fft_num, VAR_0->nbits >= 6); #ifdef ALTIVEC_USE_REFERENCE_C_CODE
int VAR_3 = VAR_0->nbits; int	VAR_4, VAR_5, VAR_6; int	VAR_7, VAR_8; register FFTComplex *VAR_9, *q; FFTComplex *exptab = VAR_0->exptab; int l; FFTSample tmp_re, tmp_im; POWERPC_TBL_START_COUNT(altivec_fft_num, VAR_0->nbits >= 6); VAR_5 = 1 << VAR_3; VAR_9=&VAR_1[0]; VAR_4=(VAR_5 >> 1); do { BF(VAR_9[0].re, VAR_9[0].im, VAR_9[1].re, VAR_9[1].im,
VAR_9[0].re, VAR_9[0].im, VAR_9[1].re, VAR_9[1].im); VAR_9+=2; } while (--VAR_4 != 0); VAR_9=&VAR_1[0]; VAR_4=VAR_5 >> 2; if (VAR_0->inverse) { do { BF(VAR_9[0].re, VAR_9[0].im, VAR_9[2].re, VAR_9[2].im,
VAR_9[0].re, VAR_9[0].im, VAR_9[2].re, VAR_9[2].im); BF(VAR_9[1].re, VAR_9[1].im, VAR_9[3].re, VAR_9[3].im,
VAR_9[1].re, VAR_9[1].im, -VAR_9[3].im, VAR_9[3].re); VAR_9+=4; } while (--VAR_4 != 0); } else { do { BF(VAR_9[0].re, VAR_9[0].im, VAR_9[2].re, VAR_9[2].im,
VAR_9[0].re, VAR_9[0].im, VAR_9[2].re, VAR_9[2].im); BF(VAR_9[1].re, VAR_9[1].im, VAR_9[3].re, VAR_9[3].im,
VAR_9[1].re, VAR_9[1].im, VAR_9[3].im, -VAR_9[3].re); VAR_9+=4; } while (--VAR_4 != 0); } VAR_7 = VAR_5 >> 3; VAR_8 = 1 << 2; VAR_6 = VAR_5 >> 1; do { VAR_9 = VAR_1; q = VAR_1 + VAR_8; for (VAR_4 = 0; VAR_4 < VAR_7; ++VAR_4) { BF(VAR_9->re, VAR_9->im, q->re, q->im,
VAR_9->re, VAR_9->im, q->re, q->im); VAR_9++; q++; for(l = VAR_7; l < VAR_6; l += VAR_7) { CMUL(tmp_re, tmp_im, exptab[l].re, exptab[l].im, q->re, q->im); BF(VAR_9->re, VAR_9->im, q->re, q->im,
VAR_9->re, VAR_9->im, tmp_re, tmp_im); VAR_9++; q++; } VAR_9 += VAR_8; q += VAR_8; } VAR_7 = VAR_7 >> 1; VAR_8 = VAR_8 << 1; } while (VAR_7 != 0); POWERPC_TBL_STOP_COUNT(altivec_fft_num, VAR_0->nbits >= 6); #else
#ifdef CONFIG_DARWIN
register const vector float VAR_2 = (const vector float)(0.); #else
register const vector float VAR_2 = (const vector float){0.,0.,0.,0.}; #endif
int VAR_3 = VAR_0->nbits; int	VAR_4, VAR_5, VAR_6; int	VAR_7, VAR_8; register FFTComplex *VAR_9, *q; FFTComplex *cptr, *cptr1; int VAR_10; POWERPC_TBL_START_COUNT(altivec_fft_num, VAR_0->nbits >= 6); VAR_5 = 1 << VAR_3; { vector float *r, a, b, a1, c1, c2; r = (vector float *)&VAR_1[0]; c1 = vcii(VAR_9,VAR_9,n,n); if (VAR_0->inverse)
{ c2 = vcii(VAR_9,VAR_9,n,VAR_9); } else
{ c2 = vcii(VAR_9,VAR_9,VAR_9,n); } VAR_4 = (VAR_5 >> 2); do { a = vec_ld(0, r); a1 = vec_ld(sizeof(vector float), r); b = vec_perm(a,a,vcprmle(1,0,3,2)); a = vec_madd(a,c1,b); b = vec_perm(a1,a1,vcprmle(1,0,3,2)); b = vec_madd(a1,c1,b); b = vec_perm(b,b,vcprmle(2,3,1,0)); vec_st(vec_madd(b,c2,a), 0, r); vec_st(vec_nmsub(b,c2,a), sizeof(vector float), r); r += 2; } while (--VAR_4 != 0); } VAR_7 = VAR_5 >> 3; VAR_8 = 1 << 2; VAR_6 = VAR_5 >> 1; cptr1 = VAR_0->exptab1; do { VAR_9 = VAR_1; q = VAR_1 + VAR_8; VAR_4 = VAR_7; do { cptr = cptr1; VAR_10 = VAR_8 >> 1; do { vector float a,b,c,t1; a = vec_ld(0, (float*)VAR_9); b = vec_ld(0, (float*)q); c = vec_ld(0, (float*)cptr); t1 = vec_madd(c, vec_perm(b,b,vcprmle(2,2,0,0)),VAR_2); c = vec_ld(sizeof(vector float), (float*)cptr); b = vec_madd(c, vec_perm(b,b,vcprmle(3,3,1,1)),t1); vec_st(vec_add(a,b), 0, (float*)VAR_9); vec_st(vec_sub(a,b), 0, (float*)q); VAR_9 += 2; q += 2; cptr += 4; } while (--VAR_10); VAR_9 += VAR_8; q += VAR_8; } while (--VAR_4); cptr1 += VAR_8 * 2; VAR_7 = VAR_7 >> 1; VAR_8 = VAR_8 << 1; } while (VAR_7 != 0); POWERPC_TBL_STOP_COUNT(altivec_fft_num, VAR_0->nbits >= 6); #endif
}