开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > 多媒体编程 > 查看源码
IDCT_MMX.C
资源名称：Visual C++视频音频开发实用工程案例精选.rar　[点击查看]
上传用户：tuheem
上传日期：2007-05-01
资源大小：21889k
文件大小：14k
源码类别：

多媒体编程

开发平台：

Visual C++

IDCT_MMX.C：源码内容
  2. #define DEFAULT_FUNCTION_NAME       idct
  6. #ifndef MMX_IDCT_NAME
  7. #define MMX_IDCT_NAME DEFAULT_FUNCTION_NAME
  8. #endif
  12. #define int16_t short
  13. #define int32_t int
  14. #define uint64_t unsigned __int64
  16. #define MPTR qword ptr 
  19. #define SHORT4_TO_QWORD(A, B, C, D)   ( 
  20.     ( ( ((uint64_t)(A)) & ((uint64_t)(0xffff)) )  <<  0 ) | 
  21.     ( ( ((uint64_t)(B)) & ((uint64_t)(0xffff)) )  << 16 ) | 
  22.     ( ( ((uint64_t)(C)) & ((uint64_t)(0xffff)) )  << 32 ) | 
  23.     ( ( ((uint64_t)(D)) & ((uint64_t)(0xffff)) )  << 48 ) ) 
  25. #define INT2_TO_QWORD(A, B)    ( 
  26.     ( ( ((uint64_t)(A)) & ((uint64_t)(0xffffffff)) )  <<  0 ) | 
  27.     ( ( ((uint64_t)(B)) & ((uint64_t)(0xffffffff)) )  << 32 ) )
  32. #define BITS_INV_ACC      5                             
  33. #define SHIFT_INV_ROW     (16 - BITS_INV_ACC)
  34. #define SHIFT_INV_COL     (1 + BITS_INV_ACC)
  35. #define RND_INV_ROW       (1024 * (6 - BITS_INV_ACC))   
  36. #define RND_INV_COL       (16 * (BITS_INV_ACC - 3))     
  37. #define RND_INV_CORR      (RND_INV_COL - 1)             
  38. #define BITS_FRW_ACC      3                             
  39. #define SHIFT_FRW_COL     BITS_FRW_ACC
  40. #define SHIFT_FRW_ROW     (BITS_FRW_ACC + 17)
  41. #define RND_FRW_ROW       (262144 * (BITS_FRW_ACC - 1)) 
  48. const uint64_t   tg_1_16 = SHORT4_TO_QWORD( 13036,  13036,  13036,  13036 );
  49. const uint64_t   tg_2_16 = SHORT4_TO_QWORD( 27146,  27146,  27146,  27146 );
  50. const uint64_t   tg_3_16 = SHORT4_TO_QWORD(-21746, -21746, -21746, -21746 );
  51. const uint64_t ocos_4_16 = SHORT4_TO_QWORD( 23170,  23170,  23170,  23170 );
  59. #if   SHIFT_INV_ROW == 12
  60. const uint64_t rounder[8] = {
  61. INT2_TO_QWORD( 65536, 65536), /* rounder_0 */       
  62. INT2_TO_QWORD(  7195,  7195), /* rounder_1 */ 
  63. INT2_TO_QWORD(  4520,  4520), /* rounder_2 */ 
  64. INT2_TO_QWORD(  2407,  2407), /* rounder_3 */ 
  65. INT2_TO_QWORD(     0,     0), /* rounder_4 */ 
  66. INT2_TO_QWORD(   240,   240), /* rounder_5 */
  67. INT2_TO_QWORD(  1024,  1024), /* rounder_6 */
  68. INT2_TO_QWORD(  1024,  1024)  /* rounder_7 */
  69. };
  70. #elif SHIFT_INV_ROW == 11
  71. const uint64_t rounder[2*8] = {
  72. INT2_TO_QWORD( 65536, 65536), /* rounder_0 */
  73. INT2_TO_QWORD(  3597,  3597), /* rounder_1 */
  74. INT2_TO_QWORD(  2260,  2260), /* rounder_2 */
  75. INT2_TO_QWORD(  1203,  1203), /* rounder_3 */
  76. INT2_TO_QWORD(   0,   0), /* rounder_4 */
  77. INT2_TO_QWORD(   120, 120), /* rounder_5 */
  78. INT2_TO_QWORD(   512, 512), /* rounder_6 */
  79. INT2_TO_QWORD(   512, 512)  /* rounder_7 */
  80. };
  81. #endif
  86. const uint64_t tab_i_04[8*8] = {
  88.   SHORT4_TO_QWORD( 16384,  16384,  16384, -16384 ),   /* movq-> w06 w04 w02 w00 */
  89. SHORT4_TO_QWORD( 21407,   8867,   8867, -21407 ),   /*        w07 w05 w03 w01 */
  90. SHORT4_TO_QWORD( 16384, -16384,  16384,  16384 ),   /*        w14 w12 w10 w08 */
  91. SHORT4_TO_QWORD( -8867,  21407, -21407,  -8867 ),   /*        w15 w13 w11 w09 */
  92. SHORT4_TO_QWORD( 22725,  12873,  19266, -22725 ),   /*        w22 w20 w18 w16 */
  93. SHORT4_TO_QWORD( 19266,   4520,  -4520, -12873 ),   /*        w23 w21 w19 w17 */
  94. SHORT4_TO_QWORD( 12873,   4520,   4520,  19266 ),   /*        w30 w28 w26 w24 */
  95. SHORT4_TO_QWORD(-22725,  19266, -12873, -22725 ),   /*        w31 w29 w27 w25 */
  97. SHORT4_TO_QWORD( 22725,  22725,  22725, -22725 ),   /* movq-> w06 w04 w02 w00 */
  98. SHORT4_TO_QWORD( 29692,  12299,  12299, -29692 ),   /*        w07 w05 w03 w01 */
  99. SHORT4_TO_QWORD( 22725, -22725,  22725,  22725 ),   /*        w14 w12 w10 w08 */
  100. SHORT4_TO_QWORD(-12299,  29692, -29692, -12299 ),   /*        w15 w13 w11 w09 */
  101. SHORT4_TO_QWORD( 31521,  17855,  26722, -31521 ),   /*        w22 w20 w18 w16 */
  102. SHORT4_TO_QWORD( 26722,   6270,  -6270, -17855 ),   /*        w23 w21 w19 w17 */
  103. SHORT4_TO_QWORD( 17855,   6270,   6270,  26722 ),   /*        w30 w28 w26 w24 */
  104. SHORT4_TO_QWORD(-31521,  26722, -17855, -31521 ),   /*        w31 w29 w27 w25 */
  106. SHORT4_TO_QWORD( 21407,  21407,  21407, -21407 ),   /* movq-> w06 w04 w02 w00 */
  107. SHORT4_TO_QWORD( 27969,  11585,  11585, -27969 ),   /*        w07 w05 w03 w01 */
  108. SHORT4_TO_QWORD( 21407, -21407,  21407,  21407 ),   /*        w14 w12 w10 w08 */
  109. SHORT4_TO_QWORD(-11585,  27969, -27969, -11585 ),   /*        w15 w13 w11 w09 */
  110. SHORT4_TO_QWORD( 29692,  16819,  25172, -29692 ),   /*        w22 w20 w18 w16 */
  111. SHORT4_TO_QWORD( 25172,   5906,  -5906, -16819 ),   /*        w23 w21 w19 w17 */
  112. SHORT4_TO_QWORD( 16819,   5906,   5906,  25172 ),   /*        w30 w28 w26 w24 */
  113. SHORT4_TO_QWORD(-29692,  25172, -16819, -29692 ),   /*        w31 w29 w27 w25 */
  115. SHORT4_TO_QWORD( 19266,  19266,  19266, -19266 ),   /* movq-> w06 w04 w02 w00 */
  116. SHORT4_TO_QWORD( 25172,  10426,  10426, -25172 ),   /*        w07 w05 w03 w01 */
  117. SHORT4_TO_QWORD( 19266, -19266,  19266,  19266 ),   /*        w14 w12 w10 w08 */
  118. SHORT4_TO_QWORD(-10426,  25172, -25172, -10426 ),   /*        w15 w13 w11 w09 */
  119. SHORT4_TO_QWORD( 26722,  15137,  22654, -26722 ),   /*        w22 w20 w18 w16 */
  120. SHORT4_TO_QWORD( 22654,   5315,  -5315, -15137 ),   /*        w23 w21 w19 w17 */
  121. SHORT4_TO_QWORD( 15137,   5315,   5315,  22654 ),   /*        w30 w28 w26 w24 */
  122. SHORT4_TO_QWORD(-26722,  22654, -15137, -26722 ),   /*        w31 w29 w27 w25 */
  124.   SHORT4_TO_QWORD( 16384,  16384,  16384, -16384 ),   /* movq-> w06 w04 w02 w00 */
  125. SHORT4_TO_QWORD( 21407,   8867,   8867, -21407 ),   /*        w07 w05 w03 w01 */
  126. SHORT4_TO_QWORD( 16384, -16384,  16384,  16384 ),   /*        w14 w12 w10 w08 */
  127. SHORT4_TO_QWORD(-8867,  21407, -21407,  -8867  ),   /*        w15 w13 w11 w09 */
  128. SHORT4_TO_QWORD( 22725,  12873,  19266, -22725 ),   /*        w22 w20 w18 w16 */
  129. SHORT4_TO_QWORD( 19266,   4520,  -4520, -12873 ),   /*        w23 w21 w19 w17 */
  130. SHORT4_TO_QWORD( 12873,   4520,   4520,  19266 ),   /*        w30 w28 w26 w24 */
  131. SHORT4_TO_QWORD(-22725,  19266, -12873, -22725 ),   /*        w31 w29 w27 w25 */
  133. SHORT4_TO_QWORD( 19266,  19266,  19266, -19266 ),   /* movq-> w06 w04 w02 w00 */
  134. SHORT4_TO_QWORD( 25172,  10426,  10426, -25172 ),   /*        w07 w05 w03 w01 */
  135. SHORT4_TO_QWORD( 19266, -19266,  19266,  19266 ),   /*        w14 w12 w10 w08 */
  136. SHORT4_TO_QWORD(-10426,  25172, -25172, -10426 ),   /*        w15 w13 w11 w09 */
  137. SHORT4_TO_QWORD( 26722,  15137,  22654, -26722 ),   /*        w22 w20 w18 w16 */
  138. SHORT4_TO_QWORD( 22654,   5315,  -5315, -15137 ),   /*        w23 w21 w19 w17 */
  139. SHORT4_TO_QWORD( 15137,   5315,   5315,  22654 ),   /*        w30 w28 w26 w24 */
  140. SHORT4_TO_QWORD(-26722,  22654, -15137, -26722 ),   /*        w31 w29 w27 w25 */
  142. SHORT4_TO_QWORD( 21407,  21407,  21407, -21407 ),   /* movq-> w06 w04 w02 w00 */
  143. SHORT4_TO_QWORD( 27969,  11585,  11585, -27969 ),   /*        w07 w05 w03 w01 */
  144. SHORT4_TO_QWORD( 21407, -21407,  21407,  21407 ),   /*        w14 w12 w10 w08 */
  145. SHORT4_TO_QWORD(-11585,  27969, -27969, -11585 ),   /*        w15 w13 w11 w09 */
  146. SHORT4_TO_QWORD( 29692,  16819,  25172, -29692 ),   /*        w22 w20 w18 w16 */
  147. SHORT4_TO_QWORD( 25172,   5906,  -5906, -16819 ),   /*        w23 w21 w19 w17 */
  148. SHORT4_TO_QWORD( 16819,   5906,   5906,  25172 ),   /*        w30 w28 w26 w24 */
  149. SHORT4_TO_QWORD(-29692,  25172, -16819, -29692 ),   /*        w31 w29 w27 w25 */
  151. SHORT4_TO_QWORD( 22725,  22725,  22725, -22725 ),   /* movq-> w06 w04 w02 w00 */
  152. SHORT4_TO_QWORD( 29692,  12299,  12299, -29692 ),   /*        w07 w05 w03 w01 */
  153. SHORT4_TO_QWORD( 22725, -22725,  22725,  22725 ),   /*        w14 w12 w10 w08 */
  154. SHORT4_TO_QWORD(-12299,  29692, -29692, -12299 ),   /*        w15 w13 w11 w09 */
  155. SHORT4_TO_QWORD( 31521,  17855,  26722, -31521 ),   /*        w22 w20 w18 w16 */
  156. SHORT4_TO_QWORD( 26722,   6270,  -6270, -17855 ),   /*        w23 w21 w19 w17 */
  157. SHORT4_TO_QWORD( 17855,   6270,   6270,  26722 ),   /*        w30 w28 w26 w24 */
  158. SHORT4_TO_QWORD(-31521,  26722, -17855, -31521 ),   /*        w31 w29 w27 w25 */
  159. };
  163. uint64_t tempMatrix[64/4];
  166. #define INP      eax
  167. #define ROUNDER  ebx
  168. #define TABLE    ecx
  169. #define OUTP     edx
  176. void MMX_IDCT_NAME
  177.               (int16_t *block) {
  178. __asm {
  180. push eax
  181. push ebx 
  182. push ecx
  183. push edx 
  184. push edi
  186. mov INP, block
  187. lea OUTP,  [tempMatrix]
  188. lea TABLE, [tab_i_04]
  189. lea ROUNDER, [rounder]
  191. mov edi, -8;                 ; i=-8
  193. align 16
  195. rowloop:
  197. movq mm0, MPTR [INP]  ; 0 ; x3 x2 x1 x0
  199. movq mm1, MPTR [INP+8] ; 1 ; x7 x6 x5 x4
  200. movq mm2, mm0  ; 2 ; x3 x2 x1 x0
  202. movq mm3, MPTR [TABLE] ; 3 ; w06 w04 w02 w00
  203. punpcklwd mm0, mm1  ; x5 x1 x4 x0
  205. movq mm5, mm0  ; 5 ; x5 x1 x4 x0
  206. punpckldq mm0, mm0  ; x4 x0 x4 x0
  208. movq mm4, MPTR [TABLE+8]  ; 4 ; w07 w05 w03 w01
  209. punpckhwd mm2, mm1 ; 1 ; x7 x3 x6 x2
  211. pmaddwd mm3, mm0  ; x4*w06+x0*w04 x4*w02+x0*w00
  212. movq mm6, mm2  ; 6  ; x7 x3 x6 x2
  214. movq mm1, MPTR [TABLE+32]  ; 1  ; w22 w20 w18 w16
  215. punpckldq mm2, mm2  ; x6 x2 x6 x2
  217. pmaddwd mm4, mm2  ; x6*w07+x2*w05 x6*w03+x2*w01
  218. punpckhdq mm5, mm5  ; x5 x1 x5 x1
  220. pmaddwd mm0, MPTR [TABLE+16]  ; x4*w14+x0*w12 x4*w10+x0*w08
  221. punpckhdq mm6, mm6  ; x7 x3 x7 x3
  223. movq mm7, MPTR [TABLE+40]  ; 7  ; w23 w21 w19 w17
  224. pmaddwd mm1, mm5  ; x5*w22+x1*w20 x5*w18+x1*w16
  226. paddd mm3, MPTR [ROUNDER]  ; +rounder
  227. pmaddwd mm7, mm6  ; x7*w23+x3*w21 x7*w19+x3*w17
  229. pmaddwd mm2, MPTR [TABLE+24]  ; x6*w15+x2*w13 x6*w11+x2*w09
  230. paddd mm3, mm4  ; 4  ; a1=sum(even1) a0=sum(even0)
  232. pmaddwd mm5, MPTR [TABLE+48]  ; x5*w30+x1*w28 x5*w26+x1*w24
  233. movq mm4, mm3  ; 4  ; a1 a0
  235. pmaddwd mm6, MPTR [TABLE+56]  ; x7*w31+x3*w29 x7*w27+x3*w25
  236. paddd mm1, mm7  ; 7  ; b1=sum(odd1) b0=sum(odd0)
  238. paddd mm0, MPTR [ROUNDER] ; +rounder
  239. psubd mm3, mm1  ; a1-b1 a0-b0
  241. psrad mm3, SHIFT_INV_ROW  ; y6=a1-b1 y7=a0-b0
  242. paddd mm1, mm4  ; 4  ; a1+b1 a0+b0
  244. paddd mm0, mm2  ; 2  ; a3=sum(even3) a2=sum(even2)
  245. psrad mm1, SHIFT_INV_ROW  ; y1=a1+b1 y0=a0+b0
  247. paddd mm5, mm6  ; 6  ; b3=sum(odd3) b2=sum(odd2)
  248. movq mm4, mm0  ; 4  ; a3 a2
  250. paddd mm0, mm5  ; a3+b3 a2+b2
  251. psubd mm4, mm5  ; 5  ; a3-b3 a2-b2
  253. psrad mm0, SHIFT_INV_ROW  ; y3=a3+b3 y2=a2+b2
  254. psrad mm4, SHIFT_INV_ROW  ; y4=a3-b3 y5=a2-b2
  256. packssdw mm1, mm0  ; 0  ; y3 y2 y1 y0
  257. packssdw mm4, mm3  ; 3  ; y6 y7 y4 y5
  259. movq mm7, mm4  ; 7  ; y6 y7 y4 y5
  260. psrld mm4, 16  ; 0 y6 0 y4
  262. pslld mm7, 16  ; y7 0 y5 0
  263. movq MPTR [OUTP], mm1  ; 1  ; save y3 y2 y1 y0
  264.                              
  265. por mm7, mm4  ; 4  ; y7 y6 y5 y4
  266. movq MPTR [OUTP+8], mm7  ; 7  ; save y7 y6 y5 y4
  269. add INP,  16                    ; add 1 row to input pointer
  270. add ROUNDER, 8                  ; go to next rounding values
  271. add OUTP, 16                    ; add 1 row to output pointer
  272. add TABLE,64                    ; move to next section of table
  274. add edi, 1
  275. jne rowloop;
  281. lea INP,   [tempMatrix] 
  282. mov OUTP,  block
  283. mov edi, -2;                 ; i=-2
  285. align 16
  286. colloop:
  288. movq mm0, tg_3_16
  290. movq mm3, MPTR [INP+16*3]
  291. movq mm1, mm0 ; tg_3_16
  293. movq mm5, MPTR [INP+16*5]
  294. pmulhw mm0, mm3 ; x3*(tg_3_16-1)
  297. movq mm4, tg_1_16
  298. pmulhw mm1, mm5 ; x5*(tg_3_16-1)
  300. movq mm7, MPTR [INP+16*7]
  301. movq mm2, mm4 ; tg_1_16
  303. movq mm6, MPTR [INP+16*1]
  304. pmulhw mm4, mm7 ; x7*tg_1_16
  306. paddsw mm0, mm3 ; x3*tg_3_16
  307. pmulhw mm2, mm6 ; x1*tg_1_16
  309. paddsw mm1, mm3 ; x3+x5*(tg_3_16-1)
  310. psubsw mm0, mm5 ; x3*tg_3_16-x5 = tm35
  312. movq mm3, ocos_4_16
  313. paddsw mm1, mm5 ; x3+x5*tg_3_16 = tp35
  315. paddsw mm4, mm6 ; x1+tg_1_16*x7 = tp17
  316. psubsw mm2, mm7 ; x1*tg_1_16-x7 = tm17
  318. movq mm5, mm4 ; tp17
  319. movq mm6, mm2 ; tm17
  321. paddsw mm5, mm1 ; tp17+tp35 = b0
  322. psubsw mm6, mm0 ; tm17-tm35 = b3
  324. psubsw mm4, mm1 ; tp17-tp35 = t1
  325. paddsw mm2, mm0 ; tm17+tm35 = t2
  327. movq mm7, tg_2_16
  328. movq mm1, mm4 ; t1
  330. ; movq MPTR [SCRATCH+0], mm5 ; save b0
  331. movq MPTR [OUTP+3*16], mm5 ; save b0
  332. paddsw mm1, mm2 ; t1+t2
  334. ; movq MPTR [SCRATCH+8], mm6 ; save b3
  335. movq MPTR [OUTP+5*16], mm6 ; save b3
  336. psubsw mm4, mm2 ; t1-t2
  338. movq mm5, MPTR [INP+2*16]
  339. movq mm0, mm7 ; tg_2_16
  341. movq mm6, MPTR [INP+6*16]
  342. pmulhw mm0, mm5 ; x2*tg_2_16
  344. pmulhw mm7, mm6 ; x6*tg_2_16
  345. ; slot
  346. pmulhw mm1, mm3 ; ocos_4_16*(t1+t2) = b1/2
  347. ; slot
  348. movq mm2, MPTR [INP+0*16]
  349. pmulhw mm4, mm3 ; ocos_4_16*(t1-t2) = b2/2
  351. psubsw mm0, mm6 ; t2*tg_2_16-x6 = tm26
  352. movq mm3, mm2 ; x0
  354. movq mm6, MPTR [INP+4*16]
  355. paddsw mm7, mm5 ; x2+x6*tg_2_16 = tp26
  357. paddsw mm2, mm6 ; x0+x4 = tp04
  358. psubsw mm3, mm6 ; x0-x4 = tm04
  360. movq mm5, mm2 ; tp04
  361. movq mm6, mm3 ; tm04
  363. psubsw mm2, mm7 ; tp04-tp26 = a3
  364. paddsw mm3, mm0 ; tm04+tm26 = a1
  366. paddsw mm1, mm1 ; b1
  367. paddsw mm4, mm4 ; b2
  369. paddsw mm5, mm7 ; tp04+tp26 = a0
  370. psubsw mm6, mm0 ; tm04-tm26 = a2
  372. movq mm7, mm3 ; a1
  373. movq mm0, mm6 ; a2
  375. paddsw mm3, mm1 ; a1+b1
  376. paddsw mm6, mm4 ; a2+b2
  378. psraw mm3, SHIFT_INV_COL ; dst1
  379. psubsw mm7, mm1 ; a1-b1
  381. psraw mm6, SHIFT_INV_COL ; dst2
  382. psubsw mm0, mm4 ; a2-b2
  384. ; movq mm1, MPTR [SCRATCH+0] ; load b0
  385. movq mm1, MPTR [OUTP+3*16] ; load b0
  386. psraw mm7, SHIFT_INV_COL ; dst6
  388. movq mm4, mm5 ; a0
  389. psraw mm0, SHIFT_INV_COL ; dst5
  391. movq MPTR [OUTP+1*16], mm3
  392. paddsw mm5, mm1 ; a0+b0
  394. movq MPTR [OUTP+2*16], mm6
  395. psubsw mm4, mm1 ; a0-b0
  397. ; movq mm3, MPTR [SCRATCH+8] ; load b3
  398. movq mm3, MPTR [OUTP+5*16] ; load b3
  399. psraw mm5, SHIFT_INV_COL ; dst0
  401. movq mm6, mm2 ; a3
  402. psraw mm4, SHIFT_INV_COL ; dst7
  404. movq MPTR [OUTP+5*16], mm0
  405. paddsw mm2, mm3 ; a3+b3
  407. movq MPTR [OUTP+6*16], mm7
  408. psubsw mm6, mm3 ; a3-b3
  410. movq MPTR [OUTP+0*16], mm5
  411. psraw mm2, SHIFT_INV_COL ; dst3
  413. movq MPTR [OUTP+7*16], mm4
  414. psraw mm6, SHIFT_INV_COL ; dst4
  416. movq MPTR [OUTP+3*16], mm2
  418. movq MPTR [OUTP+4*16], mm6
  420. lea INP,   [tempMatrix] 
  421. mov OUTP,  block
  422. add INP,  8                    ; add 4 cols to input pointer
  423. add OUTP, 8                    ; add 4 cols to output pointer
  425. add edi, 1
  426. jne colloop
  428. pop edi
  429. pop edx 
  430. pop ecx
  431. pop ebx
  432. pop eax
  434. emms 
  436. }; 
  438. };