40 #if HAVE_GNU_AS && HAVE_ALTIVEC
48 const uint16_t *revtabj = s->
revtab;
49 const uint16_t *revtabk = s->
revtab+n4;
58 vec_f cos,sin,cos0,sin0,cos1,sin1,
re,
im,r0,i0,r1,i1,
a,
b,
c,d;
59 #define CMULA(p,o0,o1,o2,o3)\
62 re = vec_perm(a, b, vcprm(0,2,s0,s2)); \
63 im = vec_perm(a, b, vcprm(s3,s1,3,1)); \
64 cos = vec_perm(cos0, cos1, vcprm(o0,o1,s##o2,s##o3)); \
65 sin = vec_perm(sin0, sin1, vcprm(o0,o1,s##o2,s##o3));\
66 r##p = im*cos - re*sin;\
67 i##p = re*cos + im*sin;
68 #define STORE2(v,dst)\
70 vec_ste(v, 0, output+j*2);\
71 vec_ste(v, 4, output+j*2);
73 a = vec_perm(r##p, i##p, vcprm(0,s0,0,s0));\
74 b = vec_perm(r##p, i##p, vcprm(1,s1,1,s1));\
75 c = vec_perm(r##p, i##p, vcprm(2,s2,2,s2));\
76 d = vec_perm(r##p, i##p, vcprm(3,s3,3,s3));\
77 STORE2(a, revtabk[ p*2-4]);\
78 STORE2(b, revtabk[ p*2-3]);\
79 STORE2(c, revtabj[-p*2+2]);\
80 STORE2(d, revtabj[-p*2+3]);
102 #define CMULB(d0,d1,o)\
107 d0 = im*sin - re*cos;\
108 d1 = re*sin + im*cos;
112 pout[2*j] = vec_perm(a, d, vcprm(0,
s3,1,
s2));
113 pout[2*j+1] = vec_perm(a, d, vcprm(2,
s1,3,
s0));
114 pout[2*k] = vec_perm(c, b, vcprm(0,
s3,1,
s2));
115 pout[2*k+1] = vec_perm(c, b, vcprm(2,
s1,3,
s0));
127 vec_u32 sign = {1
U<<31,1
U<<31,1U<<31,1U<<31};
131 imdct_half_altivec(s, output + n4, input);
133 for (k = 0; k < n16; k++) {
136 p0[-k-1] = vec_perm(a, a, vcprm(3,2,1,0));
137 p1[k] = vec_perm(b, b, vcprm(3,2,1,0));
144 #if HAVE_GNU_AS && HAVE_ALTIVEC