register const vector float vczero = (vector float)( 0., 0., 0., 0.);
register const vector float vczero = (const vector float)(0.);
void fft_calc_altivec(FFTContext *s, FFTComplex *z){    register const vector float vczero = (vector float)( 0., 0., 0., 0.);        int ln = s->nbits;    int	j, np, np2;    int	nblocks, nloops;    register FFTComplex *p, *q;    FFTComplex *cptr, *cptr1;    int k;    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);}