[libav.git] / libavcodec / i386 / fdctdata.c

//////////////////////////////////////////////////////////////////////////////
//
//  fdctam32.c - AP922 MMX(3D-Now) forward-DCT
//  ----------
//  Intel Application Note AP-922 - fast, precise implementation of DCT
//        http://developer.intel.com/vtune/cbts/appnotes.htm
//  ----------
//  
//       This routine uses a 3D-Now/MMX enhancement to increase the
//  accuracy of the fdct_col_4 macro.  The dct_col function uses 3D-Now's
//  PMHULHRW instead of MMX's PMHULHW(and POR).  The substitution improves
//  accuracy very slightly with performance penalty.  If the target CPU
//  does not support 3D-Now, then this function cannot be executed.  
//  fdctmm32.c contains the standard MMX implementation of AP-922.
//  
//  For a fast, precise MMX implementation of inverse-DCT 
//              visit http://www.elecard.com/peter
//
//  v1.0 07/22/2000 (initial release)
//       Initial release of AP922 MMX(3D-Now) forward_DCT.
//       This code was tested with Visual C++ 6.0Pro + service_pack4 + 
//       processor_pack_beta!  If you have the processor_pack_beta, you can
//       remove the #include for amd3dx.h, and substitute the 'normal'
//       assembly lines for the macro'd versions.  Otherwise, this
//       code should compile 'as is', under Visual C++ 6.0 Pro.
//     
//  liaor@iname.com  http://members.tripod.com/~liaor  
//////////////////////////////////////////////////////////////////////////////

#include <inttypes.h>

//////////////////////////////////////////////////////////////////////
//
// constants for the forward DCT
// -----------------------------
//
// Be sure to check that your compiler is aligning all constants to QWORD
// (8-byte) memory boundaries!  Otherwise the unaligned memory access will
// severely stall MMX execution.
//
//////////////////////////////////////////////////////////////////////

#define BITS_FRW_ACC	3 //; 2 or 3 for accuracy
#define SHIFT_FRW_COL	BITS_FRW_ACC
#define SHIFT_FRW_ROW	(BITS_FRW_ACC + 17)
//#define RND_FRW_ROW		(262144 * (BITS_FRW_ACC - 1)) //; 1 << (SHIFT_FRW_ROW-1)
#define RND_FRW_ROW		(1 << (SHIFT_FRW_ROW-1))
//#define RND_FRW_COL		(2 * (BITS_FRW_ACC - 1)) //; 1 << (SHIFT_FRW_COL-1)
#define RND_FRW_COL		(1 << (SHIFT_FRW_COL-1))

//concatenated table, for forward DCT transformation
const int16_t fdct_tg_all_16[] = {
	13036, 13036, 13036, 13036,		// tg * (2<<16) + 0.5
	27146, 27146, 27146, 27146,		// tg * (2<<16) + 0.5
	-21746, -21746, -21746, -21746,	// tg * (2<<16) + 0.5
	-19195, -19195, -19195, -19195,	//cos * (2<<16) + 0.5
	23170, 23170, 23170, 23170 };	//cos * (2<<15) + 0.5
const long long  fdct_one_corr = 0x0001000100010001LL;
const long fdct_r_row[2] = {RND_FRW_ROW, RND_FRW_ROW };

const int16_t tab_frw_01234567[] = {  // forward_dct coeff table
    //row0
    16384, 16384, 21407, -8867,     //    w09 w01 w08 w00
    16384, 16384, 8867, -21407,     //    w13 w05 w12 w04
    16384, -16384, 8867, 21407,     //    w11 w03 w10 w02
    -16384, 16384, -21407, -8867,   //    w15 w07 w14 w06
    22725, 12873, 19266, -22725,    //    w22 w20 w18 w16
    19266, 4520, -4520, -12873,     //    w23 w21 w19 w17
    12873, 4520, 4520, 19266,       //    w30 w28 w26 w24
    -22725, 19266, -12873, -22725,  //    w31 w29 w27 w25

    //row1
    22725, 22725, 29692, -12299,    //    w09 w01 w08 w00
    22725, 22725, 12299, -29692,    //    w13 w05 w12 w04
    22725, -22725, 12299, 29692,    //    w11 w03 w10 w02
    -22725, 22725, -29692, -12299,  //    w15 w07 w14 w06
    31521, 17855, 26722, -31521,    //    w22 w20 w18 w16
    26722, 6270, -6270, -17855,     //    w23 w21 w19 w17
    17855, 6270, 6270, 26722,       //    w30 w28 w26 w24
    -31521, 26722, -17855, -31521,  //    w31 w29 w27 w25

    //row2
    21407, 21407, 27969, -11585,    //    w09 w01 w08 w00
    21407, 21407, 11585, -27969,    //    w13 w05 w12 w04
    21407, -21407, 11585, 27969,    //    w11 w03 w10 w02
    -21407, 21407, -27969, -11585,  //    w15 w07 w14 w06
    29692, 16819, 25172, -29692,    //    w22 w20 w18 w16
    25172, 5906, -5906, -16819,     //    w23 w21 w19 w17
    16819, 5906, 5906, 25172,       //    w30 w28 w26 w24
    -29692, 25172, -16819, -29692,  //    w31 w29 w27 w25

    //row3
    19266, 19266, 25172, -10426,    //    w09 w01 w08 w00
    19266, 19266, 10426, -25172,    //    w13 w05 w12 w04
    19266, -19266, 10426, 25172,    //    w11 w03 w10 w02
    -19266, 19266, -25172, -10426,  //    w15 w07 w14 w06, 
    26722, 15137, 22654, -26722,    //    w22 w20 w18 w16
    22654, 5315, -5315, -15137,     //    w23 w21 w19 w17
    15137, 5315, 5315, 22654,       //    w30 w28 w26 w24
    -26722, 22654, -15137, -26722,  //    w31 w29 w27 w25, 

    //row4
    16384, 16384, 21407, -8867,     //    w09 w01 w08 w00
    16384, 16384, 8867, -21407,     //    w13 w05 w12 w04
    16384, -16384, 8867, 21407,     //    w11 w03 w10 w02
    -16384, 16384, -21407, -8867,   //    w15 w07 w14 w06
    22725, 12873, 19266, -22725,    //    w22 w20 w18 w16
    19266, 4520, -4520, -12873,     //    w23 w21 w19 w17
    12873, 4520, 4520, 19266,       //    w30 w28 w26 w24
    -22725, 19266, -12873, -22725,  //    w31 w29 w27 w25 

    //row5
    19266, 19266, 25172, -10426,    //    w09 w01 w08 w00
    19266, 19266, 10426, -25172,    //    w13 w05 w12 w04
    19266, -19266, 10426, 25172,    //    w11 w03 w10 w02
    -19266, 19266, -25172, -10426,  //    w15 w07 w14 w06
    26722, 15137, 22654, -26722,    //    w22 w20 w18 w16
    22654, 5315, -5315, -15137,     //    w23 w21 w19 w17
    15137, 5315, 5315, 22654,       //    w30 w28 w26 w24
    -26722, 22654, -15137, -26722,  //    w31 w29 w27 w25

    //row6
    21407, 21407, 27969, -11585,    //    w09 w01 w08 w00
    21407, 21407, 11585, -27969,    //    w13 w05 w12 w04
    21407, -21407, 11585, 27969,    //    w11 w03 w10 w02
    -21407, 21407, -27969, -11585,  //    w15 w07 w14 w06, 
    29692, 16819, 25172, -29692,    //    w22 w20 w18 w16
    25172, 5906, -5906, -16819,     //    w23 w21 w19 w17
    16819, 5906, 5906, 25172,       //    w30 w28 w26 w24
    -29692, 25172, -16819, -29692,  //    w31 w29 w27 w25, 

    //row7
    22725, 22725, 29692, -12299,    //    w09 w01 w08 w00
    22725, 22725, 12299, -29692,    //    w13 w05 w12 w04
    22725, -22725, 12299, 29692,    //    w11 w03 w10 w02
    -22725, 22725, -29692, -12299,  //    w15 w07 w14 w06, 
    31521, 17855, 26722, -31521,    //    w22 w20 w18 w16
    26722, 6270, -6270, -17855,     //    w23 w21 w19 w17
    17855, 6270, 6270, 26722,       //    w30 w28 w26 w24
    -31521, 26722, -17855, -31521   //    w31 w29 w27 w25
};
Commit	Line	Data
de6d9b64 FB	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