28 #define SUINT32 int32_t
30 #define SUINT unsigned
31 #define SUINT32 uint32_t
35 # define dct32 ff_dct32_float
36 # define FIXHR(x) ((float)(x))
37 # define MULH3(x, y, s) ((s)*(y)*(x))
38 # define INTFLOAT float
39 # define SUINTFLOAT float
41 # define dct32 ff_dct32_fixed
42 # define FIXHR(a) ((int)((a) * (1LL<<32) + 0.5))
43 # define MULH3(x, y, s) MULH((s)*(x), y)
45 # define SUINTFLOAT SUINT
53 #define COS0_0 FIXHR(0.50060299823519630134/2)
54 #define COS0_1 FIXHR(0.50547095989754365998/2)
55 #define COS0_2 FIXHR(0.51544730992262454697/2)
56 #define COS0_3 FIXHR(0.53104259108978417447/2)
57 #define COS0_4 FIXHR(0.55310389603444452782/2)
58 #define COS0_5 FIXHR(0.58293496820613387367/2)
59 #define COS0_6 FIXHR(0.62250412303566481615/2)
60 #define COS0_7 FIXHR(0.67480834145500574602/2)
61 #define COS0_8 FIXHR(0.74453627100229844977/2)
62 #define COS0_9 FIXHR(0.83934964541552703873/2)
63 #define COS0_10 FIXHR(0.97256823786196069369/2)
64 #define COS0_11 FIXHR(1.16943993343288495515/4)
65 #define COS0_12 FIXHR(1.48416461631416627724/4)
66 #define COS0_13 FIXHR(2.05778100995341155085/8)
67 #define COS0_14 FIXHR(3.40760841846871878570/8)
68 #define COS0_15 FIXHR(10.19000812354805681150/32)
70 #define COS1_0 FIXHR(0.50241928618815570551/2)
71 #define COS1_1 FIXHR(0.52249861493968888062/2)
72 #define COS1_2 FIXHR(0.56694403481635770368/2)
73 #define COS1_3 FIXHR(0.64682178335999012954/2)
74 #define COS1_4 FIXHR(0.78815462345125022473/2)
75 #define COS1_5 FIXHR(1.06067768599034747134/4)
76 #define COS1_6 FIXHR(1.72244709823833392782/4)
77 #define COS1_7 FIXHR(5.10114861868916385802/16)
79 #define COS2_0 FIXHR(0.50979557910415916894/2)
80 #define COS2_1 FIXHR(0.60134488693504528054/2)
81 #define COS2_2 FIXHR(0.89997622313641570463/2)
82 #define COS2_3 FIXHR(2.56291544774150617881/8)
84 #define COS3_0 FIXHR(0.54119610014619698439/2)
85 #define COS3_1 FIXHR(1.30656296487637652785/4)
87 #define COS4_0 FIXHR(M_SQRT1_2/2)
90 #define BF(a, b, c, s)\
92 tmp0 = val##a + val##b;\
93 tmp1 = val##a - val##b;\
95 val##b = MULH3(tmp1, c, 1<<(s));\
98 #define BF0(a, b, c, s)\
100 tmp0 = tab[a] + tab[b];\
101 tmp1 = tab[a] - tab[b];\
103 val##b = MULH3(tmp1, c, 1<<(s));\
106 #define BF1(a, b, c, d)\
108 BF(a, b, COS4_0, 1);\
109 BF(c, d,-COS4_0, 1);\
113 #define BF2(a, b, c, d)\
115 BF(a, b, COS4_0, 1);\
116 BF(c, d,-COS4_0, 1);\
123 #define ADD(a, b) val##a += val##b
130 SUINTFLOAT val0 , val1 , val2 , val3 , val4 , val5 , val6 , val7 ,
131 val8 , val9 , val10, val11, val12, val13, val14, val15,
132 val16, val17, val18, val19, val20, val21, val22, val23,
133 val24, val25, val26, val27, val28, val29, val30, val31;
271 out[ 1] = val16 + val24;
272 out[17] = val17 + val25;
273 out[ 9] = val18 + val26;
274 out[25] = val19 + val27;
275 out[ 5] = val20 + val28;
276 out[21] = val21 + val29;
277 out[13] = val22 + val30;
278 out[29] = val23 + val31;
279 out[ 3] = val24 + val20;
280 out[19] = val25 + val21;
281 out[11] = val26 + val22;
282 out[27] = val27 + val23;
283 out[ 7] = val28 + val18;
284 out[23] = val29 + val19;
285 out[15] = val30 + val17;
common internal API header
static const struct twinvq_data tab