|
0
|
1 //////////////////////////////////////////////////////////////////////////////
|
|
|
2 //
|
|
|
3 // fdctam32.c - AP922 MMX(3D-Now) forward-DCT
|
|
|
4 // ----------
|
|
|
5 // Intel Application Note AP-922 - fast, precise implementation of DCT
|
|
|
6 // http://developer.intel.com/vtune/cbts/appnotes.htm
|
|
|
7 // ----------
|
|
|
8 //
|
|
|
9 // This routine uses a 3D-Now/MMX enhancement to increase the
|
|
|
10 // accuracy of the fdct_col_4 macro. The dct_col function uses 3D-Now's
|
|
|
11 // PMHULHRW instead of MMX's PMHULHW(and POR). The substitution improves
|
|
|
12 // accuracy very slightly with performance penalty. If the target CPU
|
|
|
13 // does not support 3D-Now, then this function cannot be executed.
|
|
|
14 // fdctmm32.c contains the standard MMX implementation of AP-922.
|
|
|
15 //
|
|
|
16 // For a fast, precise MMX implementation of inverse-DCT
|
|
|
17 // visit http://www.elecard.com/peter
|
|
|
18 //
|
|
|
19 // v1.0 07/22/2000 (initial release)
|
|
|
20 // Initial release of AP922 MMX(3D-Now) forward_DCT.
|
|
|
21 // This code was tested with Visual C++ 6.0Pro + service_pack4 +
|
|
|
22 // processor_pack_beta! If you have the processor_pack_beta, you can
|
|
|
23 // remove the #include for amd3dx.h, and substitute the 'normal'
|
|
|
24 // assembly lines for the macro'd versions. Otherwise, this
|
|
|
25 // code should compile 'as is', under Visual C++ 6.0 Pro.
|
|
|
26 //
|
|
|
27 // liaor@iname.com http://members.tripod.com/~liaor
|
|
|
28 //////////////////////////////////////////////////////////////////////////////
|
|
|
29
|
|
|
30 #include <inttypes.h>
|
|
|
31
|
|
|
32 //////////////////////////////////////////////////////////////////////
|
|
|
33 //
|
|
|
34 // constants for the forward DCT
|
|
|
35 // -----------------------------
|
|
|
36 //
|
|
|
37 // Be sure to check that your compiler is aligning all constants to QWORD
|
|
|
38 // (8-byte) memory boundaries! Otherwise the unaligned memory access will
|
|
|
39 // severely stall MMX execution.
|
|
|
40 //
|
|
|
41 //////////////////////////////////////////////////////////////////////
|
|
|
42
|
|
|
43 #define BITS_FRW_ACC 3 //; 2 or 3 for accuracy
|
|
|
44 #define SHIFT_FRW_COL BITS_FRW_ACC
|
|
|
45 #define SHIFT_FRW_ROW (BITS_FRW_ACC + 17)
|
|
|
46 //#define RND_FRW_ROW (262144 * (BITS_FRW_ACC - 1)) //; 1 << (SHIFT_FRW_ROW-1)
|
|
|
47 #define RND_FRW_ROW (1 << (SHIFT_FRW_ROW-1))
|
|
|
48 //#define RND_FRW_COL (2 * (BITS_FRW_ACC - 1)) //; 1 << (SHIFT_FRW_COL-1)
|
|
|
49 #define RND_FRW_COL (1 << (SHIFT_FRW_COL-1))
|
|
|
50
|
|
|
51 //concatenated table, for forward DCT transformation
|
|
|
52 const int16_t fdct_tg_all_16[] = {
|
|
|
53 13036, 13036, 13036, 13036, // tg * (2<<16) + 0.5
|
|
|
54 27146, 27146, 27146, 27146, // tg * (2<<16) + 0.5
|
|
|
55 -21746, -21746, -21746, -21746, // tg * (2<<16) + 0.5
|
|
|
56 -19195, -19195, -19195, -19195, //cos * (2<<16) + 0.5
|
|
|
57 23170, 23170, 23170, 23170 }; //cos * (2<<15) + 0.5
|
|
|
58 const long long fdct_one_corr = 0x0001000100010001LL;
|
|
|
59 const long fdct_r_row[2] = {RND_FRW_ROW, RND_FRW_ROW };
|
|
|
60
|
|
|
61 const int16_t tab_frw_01234567[] = { // forward_dct coeff table
|
|
|
62 //row0
|
|
|
63 16384, 16384, 21407, -8867, // w09 w01 w08 w00
|
|
|
64 16384, 16384, 8867, -21407, // w13 w05 w12 w04
|
|
|
65 16384, -16384, 8867, 21407, // w11 w03 w10 w02
|
|
|
66 -16384, 16384, -21407, -8867, // w15 w07 w14 w06
|
|
|
67 22725, 12873, 19266, -22725, // w22 w20 w18 w16
|
|
|
68 19266, 4520, -4520, -12873, // w23 w21 w19 w17
|
|
|
69 12873, 4520, 4520, 19266, // w30 w28 w26 w24
|
|
|
70 -22725, 19266, -12873, -22725, // w31 w29 w27 w25
|
|
|
71
|
|
|
72 //row1
|
|
|
73 22725, 22725, 29692, -12299, // w09 w01 w08 w00
|
|
|
74 22725, 22725, 12299, -29692, // w13 w05 w12 w04
|
|
|
75 22725, -22725, 12299, 29692, // w11 w03 w10 w02
|
|
|
76 -22725, 22725, -29692, -12299, // w15 w07 w14 w06
|
|
|
77 31521, 17855, 26722, -31521, // w22 w20 w18 w16
|
|
|
78 26722, 6270, -6270, -17855, // w23 w21 w19 w17
|
|
|
79 17855, 6270, 6270, 26722, // w30 w28 w26 w24
|
|
|
80 -31521, 26722, -17855, -31521, // w31 w29 w27 w25
|
|
|
81
|
|
|
82 //row2
|
|
|
83 21407, 21407, 27969, -11585, // w09 w01 w08 w00
|
|
|
84 21407, 21407, 11585, -27969, // w13 w05 w12 w04
|
|
|
85 21407, -21407, 11585, 27969, // w11 w03 w10 w02
|
|
|
86 -21407, 21407, -27969, -11585, // w15 w07 w14 w06
|
|
|
87 29692, 16819, 25172, -29692, // w22 w20 w18 w16
|
|
|
88 25172, 5906, -5906, -16819, // w23 w21 w19 w17
|
|
|
89 16819, 5906, 5906, 25172, // w30 w28 w26 w24
|
|
|
90 -29692, 25172, -16819, -29692, // w31 w29 w27 w25
|
|
|
91
|
|
|
92 //row3
|
|
|
93 19266, 19266, 25172, -10426, // w09 w01 w08 w00
|
|
|
94 19266, 19266, 10426, -25172, // w13 w05 w12 w04
|
|
|
95 19266, -19266, 10426, 25172, // w11 w03 w10 w02
|
|
|
96 -19266, 19266, -25172, -10426, // w15 w07 w14 w06,
|
|
|
97 26722, 15137, 22654, -26722, // w22 w20 w18 w16
|
|
|
98 22654, 5315, -5315, -15137, // w23 w21 w19 w17
|
|
|
99 15137, 5315, 5315, 22654, // w30 w28 w26 w24
|
|
|
100 -26722, 22654, -15137, -26722, // w31 w29 w27 w25,
|
|
|
101
|
|
|
102 //row4
|
|
|
103 16384, 16384, 21407, -8867, // w09 w01 w08 w00
|
|
|
104 16384, 16384, 8867, -21407, // w13 w05 w12 w04
|
|
|
105 16384, -16384, 8867, 21407, // w11 w03 w10 w02
|
|
|
106 -16384, 16384, -21407, -8867, // w15 w07 w14 w06
|
|
|
107 22725, 12873, 19266, -22725, // w22 w20 w18 w16
|
|
|
108 19266, 4520, -4520, -12873, // w23 w21 w19 w17
|
|
|
109 12873, 4520, 4520, 19266, // w30 w28 w26 w24
|
|
|
110 -22725, 19266, -12873, -22725, // w31 w29 w27 w25
|
|
|
111
|
|
|
112 //row5
|
|
|
113 19266, 19266, 25172, -10426, // w09 w01 w08 w00
|
|
|
114 19266, 19266, 10426, -25172, // w13 w05 w12 w04
|
|
|
115 19266, -19266, 10426, 25172, // w11 w03 w10 w02
|
|
|
116 -19266, 19266, -25172, -10426, // w15 w07 w14 w06
|
|
|
117 26722, 15137, 22654, -26722, // w22 w20 w18 w16
|
|
|
118 22654, 5315, -5315, -15137, // w23 w21 w19 w17
|
|
|
119 15137, 5315, 5315, 22654, // w30 w28 w26 w24
|
|
|
120 -26722, 22654, -15137, -26722, // w31 w29 w27 w25
|
|
|
121
|
|
|
122 //row6
|
|
|
123 21407, 21407, 27969, -11585, // w09 w01 w08 w00
|
|
|
124 21407, 21407, 11585, -27969, // w13 w05 w12 w04
|
|
|
125 21407, -21407, 11585, 27969, // w11 w03 w10 w02
|
|
|
126 -21407, 21407, -27969, -11585, // w15 w07 w14 w06,
|
|
|
127 29692, 16819, 25172, -29692, // w22 w20 w18 w16
|
|
|
128 25172, 5906, -5906, -16819, // w23 w21 w19 w17
|
|
|
129 16819, 5906, 5906, 25172, // w30 w28 w26 w24
|
|
|
130 -29692, 25172, -16819, -29692, // w31 w29 w27 w25,
|
|
|
131
|
|
|
132 //row7
|
|
|
133 22725, 22725, 29692, -12299, // w09 w01 w08 w00
|
|
|
134 22725, 22725, 12299, -29692, // w13 w05 w12 w04
|
|
|
135 22725, -22725, 12299, 29692, // w11 w03 w10 w02
|
|
|
136 -22725, 22725, -29692, -12299, // w15 w07 w14 w06,
|
|
|
137 31521, 17855, 26722, -31521, // w22 w20 w18 w16
|
|
|
138 26722, 6270, -6270, -17855, // w23 w21 w19 w17
|
|
|
139 17855, 6270, 6270, 26722, // w30 w28 w26 w24
|
|
|
140 -31521, 26722, -17855, -31521 // w31 w29 w27 w25
|
|
|
141 };
|
|
|
142
|
|
|
143
|
