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

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

多媒体编程

开发平台：

Visual C++

deblock_vert_lpf9.c：源码内容
  3. #include "postprocess_mmx.h"
  8. const static uint64_t mm_fours  = 0x0004000400040004;
  11. INLINE void deblock_vert_lpf9(uint64_t *v_local, uint64_t *p1p2, uint8_t *v, int stride) {
  13. #ifdef PP_SELF_CHECK
  14. int j, k;
  15. uint8_t selfcheck[64], *vv;
  16. int p1, p2, psum;
  17. int l1 = 1 * stride;
  18. int l2 = 2 * stride;
  19. int l3 = 3 * stride;
  20. int l4 = 4 * stride;
  21. int l5 = 5 * stride;
  22. int l6 = 6 * stride;
  23. int l7 = 7 * stride;
  24. int l8 = 8 * stride;
  25. #endif
  28. #ifdef PP_SELF_CHECK
  29. for (j=0; j<8; j++) { 
  30. vv = &(v[j]);
  31. p1 = ((uint16_t *)(&(p1p2[0+j/4])))[j%4]; 
  32. p2 = ((uint16_t *)(&(p1p2[2+j/4])))[j%4]; 
  33. psum = p1 + p1 + p1 + vv[l1] + vv[l2] + vv[l3] + vv[l4] + 4; 
  34. selfcheck[j+8*0] = (((psum + vv[l1]) << 1) - (vv[l4] - vv[l5])) >> 4; 
  35. psum += vv[l5] - p1; 
  36. selfcheck[j+8*1] = (((psum + vv[l2]) << 1) - (vv[l5] - vv[l6])) >> 4; 
  37. psum += vv[l6] - p1; 
  38. selfcheck[j+8*2] = (((psum + vv[l3]) << 1) - (vv[l6] - vv[l7])) >> 4; 
  39. psum += vv[l7] - p1; 
  40. selfcheck[j+8*3] = (((psum + vv[l4]) << 1) + p1 - vv[l1] - (vv[l7] - vv[l8])) >> 4; 
  41. psum += vv[l8] - vv[l1];  
  42. selfcheck[j+8*4] = (((psum + vv[l5]) << 1) + (vv[l1] - vv[l2]) - vv[l8] + p2) >> 4; 
  43. psum += p2 - vv[l2];  
  44. selfcheck[j+8*5] = (((psum + vv[l6]) << 1) + (vv[l2] - vv[l3])) >> 4; 
  45. psum += p2 - vv[l3]; 
  46. selfcheck[j+8*6] = (((psum + vv[l7]) << 1) + (vv[l3] - vv[l4])) >> 4; 
  47. psum += p2 - vv[l4]; 
  48. selfcheck[j+8*7] = (((psum + vv[l8]) << 1) + (vv[l4] - vv[l5])) >> 4; 
  49. }
  50. #endif
  53. __asm {
  54. push eax
  55. push ebx
  56. push ecx
  58. mov eax, p1p2
  59. mov ebx, v_local
  60. mov ecx, v
  62. movq   mm2, [eax]                  /* mm2 = p1p2[0]               0 1 2w3 4 5 6 7    */
  63. add   ecx, stride                    /* ecx points at v[1*stride]   0 1 2 3 4 5 6 7    */     
  65. movq   mm3, 8[eax]                 /* mm3 = p1p2[1]               0 1 2 3w4 5 6 7    */
  67. movq   mm4, mm_fours                /* mm4 = 0x0004000400040004    0 1 2 3 4w5 6 7    */
  68.         */
  69. movq   mm5, mm4                     /* mm5 = 0x0004000400040004    0 1 2 3 4 5w6 7    */
  71. paddsw mm4, 16[ebx]                 /* mm4 += vv[1] left           0 1 2 3 4m5 6 7    */
  72. paddw  mm5, mm3                     /* mm5 += p2 left              0 1 2 3r4 5m6 7    */
  74. paddsw mm4, 32[ebx]                 /* mm4 += vv[2] left           0 1 2 3 4m5 6 7    */
  75. paddw  mm5, mm3                     /* mm5 += p2 left              0 1 2 3r4 5m6 7    */
  77. paddsw mm4, 48[ebx]                 /* mm4 += vv[3] left           0 1 2 3 4m5 6 7    */
  78. paddw  mm5, mm3                     /* mm5 += p2 left              0 1 2 3r4 5m6 7    */
  80. paddsw mm5, 24[ebx]                 /* mm5 += vv[1] right          0 1 2 3 4 5m6 7    */
  81. paddw  mm4, mm2                     /* mm4 += p1 left              0 1 2r3 4m5 6 7    */
  83. paddsw mm5, 40[ebx]                 /* mm5 += vv[2] right          0 1 2 3 4 5m6 7    */
  84. paddw  mm4, mm2                     /* mm4 += p1 left              0 1 2r3 4m5 6 7    */
  86. paddsw mm5, 56[ebx]                 /* mm5 += vv[3] right          0 1 2 3 4 5m6 7    */
  87. paddw  mm4, mm2                     /* mm4 += p1 left              0 1 2r3 4m5 6 7    */
  89. paddsw mm4, 64[ebx]                 /* mm4 += vv[4] left           0 1 2 3 4m5 6 7    */
  91. paddsw mm5, 72[ebx]                 /* mm5 += vv[4] right          0 1 2 3 4 5m6 7    */
  93. movq   mm0, mm4                     /* mm0 = psum left             0w1 2 3 4 5 6 7    */ 
  95. paddsw mm0, 16[ebx]                 /* mm0 += vv[1] left           0m1 2 3 4 5 6 7    */
  96. movq   mm1, mm5                     /* mm1 = psum right            0 1w2 3 4 5r6 7    */ 
  98. paddsw mm1, 24[ebx]                 /* mm1 += vv[1] right          0 1 2 3 4 5 6 7    */
  99. psllw  mm0, 1                       /* mm0 <<= 1                   0m1 2 3 4 5 6 7    */
  101. psubsw mm0, 64[ebx]                 /* mm0 -= vv[4] left           0m1 2 3 4 5 6 7    */
  102. psllw  mm1, 1                       /* mm1 <<= 1                   0 1 2 3 4 5 6 7    */
  104. psubsw mm1, 72[ebx]                 /* mm1 -= vv[4] right          0 1m2 3 4 5 6 7    */
  106. paddsw mm0, 80[ebx]                 /* mm0 += vv[5] left           0m1 2 3 4 5 6 7    */
  108. paddsw mm1, 88[ebx]                 /* mm1 += vv[5] right          0 1m2 3 4 5 6 7    */
  109. psrlw  mm0, 4                       /* mm0 >>= 4                   0m1 2 3 4 5 6 7    */
  111. paddsw mm4, 80[ebx]                 /* mm4 += vv[5] left           0 1 2 3 4m5 6 7    */
  112. psrlw  mm1, 4                       /* mm1 >>= 4                   0 1m2 3 4 5 6 7    */
  114. paddsw mm5, 88[ebx]                 /* mm5 += vv[5] right          0 1 2 3 4 5 6 7    */
  115. psubsw mm4, [eax]                  /* mm4 -= p1 left              0 1 2 3 4 5 6 7    */
  117. packuswb mm0, mm1                   /* pack mm1, mm0 to mm0        0m1 2 3 4 5 6 7    */
  118. psubsw mm5, 8[eax]                 /* mm5 -= p1 right             0 1 2 3 4 5 6 7    */
  121. movq   mm6, mm4                     /* mm6 = psum left             0 1 2 3 4 5 6 7    */ 
  123. paddsw mm6, 32[ebx]                 /* mm6 += vv[2] left           0 1 2 3 4 5 6 7    */
  124. movq   mm7, mm5                     /* mm7 = psum right            0 1 2 3 4 5 6 7    */ 
  126. paddsw mm7, 40[ebx]                 /* mm7 += vv[2] right          0 1 2 3 4 5 6 7    */
  127. psllw  mm6, 1                       /* mm6 <<= 1                   0 1 2 3 4 5 6 7    */
  129. psubsw mm6, 80[ebx]                 /* mm6 -= vv[5] left           0 1 2 3 4 5 6 7    */
  130. psllw  mm7, 1                       /* mm7 <<= 1                   0 1 2 3 4 5 6 7    */
  132. psubsw mm7, 88[ebx]                 /* mm7 -= vv[5] right          0 1 2 3 4 5 6 7    */
  134. movq   [ecx], mm0                     /* v[1*stride] = mm0           0 1 2 3 4 5 6 7    */
  136. paddsw mm6, 96[ebx]                 /* mm6 += vv[6] left           0 1 2 3 4 5 6 7    */
  137. add   ecx, stride                    /* ecx points at v[2*stride]   0 1 2 3 4 5 6 7    */     
  139. paddsw mm7, 104[ebx]                /* mm7 += vv[6] right          0 1 2 3 4 5 6 7    */
  142. paddsw mm4, 96[ebx]                 /* mm4 += vv[6] left           0 1 2 3 4 5 6 7    */
  143. psrlw  mm6, 4                       /* mm6 >>= 4                   0 1 2 3 4 5 6 7    */
  145. paddsw mm5, 104[ebx]                /* mm5 += vv[6] right          0 1 2 3 4 5 6 7    */
  146. psrlw  mm7, 4                       /* mm7 >>= 4                   0 1 2 3 4 5 6 7    */
  148. psubsw mm4, [eax]                  /* mm4 -= p1 left              0 1 2 3 4 5 6 7    */
  149. packuswb mm6, mm7                   /* pack mm7, mm6 to mm6        0 1 2 3 4 5 6 7    */
  151. psubsw mm5, 8[eax]                 /* mm5 -= p1 right             0 1 2 3 4 5 6 7    */
  154. movq   mm0, mm4                     /* mm0 = psum left             0 1 2 3 4 5 6 7    */ 
  156. paddsw mm0, 48[ebx]                 /* mm0 += vv[3] left           0 1 2 3 4 5 6 7    */
  157. movq   mm1, mm5                     /* mm1 = psum right            0 1 2 3 4 5 6 7    */ 
  159. paddsw mm1, 56[ebx]                 /* mm1 += vv[3] right          0 1 2 3 4 5 6 7    */
  160. psllw  mm0, 1                       /* mm0 <<= 1                   0 1 2 3 4 5 6 7    */
  162. psubsw mm0, 96[ebx]                 /* mm0 -= vv[6] left           0 1 2 3 4 5 6 7    */
  163. psllw  mm1, 1                       /* mm1 <<= 1                   0 1 2 3 4 5 6 7    */
  165. psubsw mm1, 104[ebx]                /* mm1 -= vv[6] right          0 1 2 3 4 5 6 7    */
  167. movq   [ecx], mm6                     /* v[2*stride] = mm6           0 1 2 3 4 5 6 7    */
  168. paddsw mm0, 112[ebx]                /* mm0 += vv[7] left           0 1 2 3 4 5 6 7    */
  170. paddsw mm1, 120[ebx]                /* mm1 += vv[7] right          0 1 2 3 4 5 6 7    */
  171. add   ecx, stride                    /* ecx points at v[3*stride]   0 1 2 3 4 5 6 7    */     
  174. paddsw mm4, 112[ebx]                /* mm4 += vv[5] left           0 1 2 3 4 5 6 7    */
  175. psrlw  mm0, 4                       /* mm0 >>= 4                   0 1 2 3 4 5 6 7    */
  177. paddsw mm5, 120[ebx]                /* mm5 += vv[5] right          0 1 2 3 4 5 6 7    */
  178. psrlw  mm1, 4                       /* mm1 >>= 4                   0 1 2 3 4 5 6 7    */
  180. psubsw mm4, [eax]                  /* mm4 -= p1 left              0 1 2 3 4 5 6 7    */
  181. packuswb mm0, mm1                   /* pack mm1, mm0 to mm0        0 1 2 3 4 5 6 7    */
  183. psubsw mm5, 8[eax]                 /* mm5 -= p1 right             0 1 2 3 4 5 6 7    */
  185. movq   [ecx], mm0                     /* v[3*stride] = mm0           0 1 2 3 4 5 6 7    */
  186. movq   mm6, mm4                     /* mm6 = psum left             0 1 2 3 4 5 6 7    */ 
  188. paddsw mm6, 64[ebx]                 /* mm6 += vv[4] left           0 1 2 3 4 5 6 7    */
  189. movq   mm7, mm5                     /* mm7 = psum right            0 1 2 3 4 5 6 7    */ 
  191. paddsw mm7, 72[ebx]                 /* mm7 += vv[4] right          0 1 2 3 4 5 6 7    */
  192. psllw  mm6, 1                       /* mm6 <<= 1                   0 1 2 3 4 5 6 7    */
  194. paddsw mm6, [eax]                  /* mm6 += p1 left              0 1 2 3 4 5 6 7    */
  195. psllw  mm7, 1                       /* mm7 <<= 1                   0 1 2 3 4 5 6 7    */
  197. paddsw mm7, 8[eax]                 /* mm7 += p1 right             0 1 2 3 4 5 6 7    */
  199. psubsw mm6, 16[ebx]                 /* mm6 -= vv[1] left           0 1 2 3 4 5 6 7    */
  201. psubsw mm7, 24[ebx]                 /* mm7 -= vv[1] right          0 1 2 3 4 5 6 7    */
  203. psubsw mm6, 112[ebx]                /* mm6 -= vv[7] left           0 1 2 3 4 5 6 7    */
  205. psubsw mm7, 120[ebx]                /* mm7 -= vv[7] right          0 1 2 3 4 5 6 7    */
  207. paddsw mm6, 128[ebx]                /* mm6 += vv[8] left           0 1 2 3 4 5 6 7    */
  208. add   ecx, stride                    /* ecx points at v[4*stride]   0 1 2 3 4 5 6 7    */     
  210. paddsw mm7, 136[ebx]                /* mm7 += vv[8] right          0 1 2 3 4 5 6 7    */
  213. paddsw mm4, 128[ebx]                /* mm4 += vv[5] left           0 1 2 3 4 5 6 7    */
  214. psrlw  mm6, 4                       /* mm6 >>= 4                   0 1 2 3 4 5 6 7    */
  216. paddsw mm5, 136[ebx]                /* mm5 += vv[5] right          0 1 2 3 4 5 6 7    */
  217. psrlw  mm7, 4                       /* mm7 >>= 4                   0 1 2 3 4 5 6 7    */
  219. psubsw mm4, 16[ebx]                 /* mm4 -= vv[1] left           0 1 2 3 4 5 6 7    */
  220. packuswb mm6, mm7                   /* pack mm7, mm6 to mm6        0 1 2 3 4 5 6 7    */
  222. psubsw mm5, 24[ebx]                 /* mm5 -= vv[1] right          0 1 2 3 4 5 6 7    */
  224. movq   mm0, mm4                     /* mm0 = psum left             0 1 2 3 4 5 6 7    */ 
  226. paddsw mm0, 80[ebx]                 /* mm0 += vv[5] left           0 1 2 3 4 5 6 7    */
  227. movq   mm1, mm5                     /* mm1 = psum right            0 1 2 3 4 5 6 7    */ 
  229. paddsw mm1, 88[ebx]                 /* mm1 += vv[5] right          0 1 2 3 4 5 6 7    */
  230. psllw  mm0, 1                       /* mm0 <<= 1                   0 1 2 3 4 5 6 7    */
  232. paddsw mm0, 16[eax]                /* mm0 += p2 left              0 1 2 3 4 5 6 7    */
  233. psllw  mm1, 1                       /* mm1 <<= 1                   0 1 2 3 4 5 6 7    */
  235. paddsw mm1, 24[eax]                /* mm1 += p2 right             0 1 2 3 4 5 6 7    */
  237. paddsw mm0, 16[ebx]                 /* mm0 += vv[1] left           0 1 2 3 4 5 6 7    */
  238. movq   [ecx], mm6                     /* v[4*stride] = mm6           0 1 2 3 4 5 6 7    */
  240. paddsw mm1, 24[ebx]                 /* mm1 += vv[1] right          0 1 2 3 4 5 6 7    */
  242. psubsw mm0, 32[ebx]                 /* mm0 -= vv[2] left           0 1 2 3 4 5 6 7    */
  244. psubsw mm1, 40[ebx]                 /* mm1 -= vv[2] right          0 1 2 3 4 5 6 7    */
  246. psubsw mm0, 128[ebx]                /* mm0 -= vv[8] left           0 1 2 3 4 5 6 7    */
  248. psubsw mm1, 136[ebx]                /* mm1 -= vv[8] right          0 1 2 3 4 5 6 7    */
  250. paddsw mm4, 16[eax]                /* mm4 += p2 left              0 1 2 3 4 5 6 7    */
  251. add   ecx, stride                   /* ecx points at v[5*stride]   0 1 2 3 4 5 6 7    */     
  253. paddsw mm5, 24[eax]                /* mm5 += p2 right             0 1 2 3 4 5 6 7    */
  255. psubsw mm4, 32[ebx]                 /* mm4 -= vv[2] left           0 1 2 3 4 5 6 7    */
  257. psubsw mm5, 40[ebx]                 /* mm5 -= vv[2] right          0 1 2 3 4 5 6 7    */
  259. movq   mm6, mm4                     /* mm6 = psum left             0 1 2 3 4 5 6 7    */ 
  261. paddsw mm6, 96[ebx]                 /* mm6 += vv[6] left           0 1 2 3 4 5 6 7    */
  262. movq   mm7, mm5                     /* mm7 = psum right            0 1 2 3 4 5 6 7    */ 
  264. paddsw mm7, 104[ebx]                /* mm7 += vv[6] right          0 1 2 3 4 5 6 7    */
  265. psllw  mm6, 1                       /* mm6 <<= 1                   0 1 2 3 4 5 6 7    */
  267. paddsw mm6, 32[ebx]                 /* mm6 += vv[2] left           0 1 2 3 4 5 6 7    */
  268. psllw  mm7, 1                       /* mm7 <<= 1                   0 1 2 3 4 5 6 7    */
  270. paddsw mm7, 40[ebx]                 /* mm7 += vv[2] right          0 1 2 3 4 5 6 7    */
  271. psrlw  mm0, 4                       /* mm0 >>= 4                   0 1 2 3 4 5 6 7    */
  273. psubsw mm6, 48[ebx]                 /* mm6 -= vv[3] left           0 1 2 3 4 5 6 7    */
  274. psrlw  mm1, 4                       /* mm1 >>= 4                   0 1 2 3 4 5 6 7    */
  276. psubsw mm7, 56[ebx]                 /* mm7 -= vv[3] right          0 1 2 3 4 5 6 7    */
  277. packuswb mm0, mm1                   /* pack mm1, mm0 to mm0        0 1 2 3 4 5 6 7    */
  279. movq   [ecx], mm0                     /* v[5*stride] = mm0           0 1 2 3 4 5 6 7    */
  281. paddsw mm4, 16[eax]                /* mm4 += p2 left              0 1 2 3 4 5 6 7    */
  282. psrlw  mm6, 4                       /* mm6 >>= 4                   0 1 2 3 4 5 6 7    */
  284. paddsw mm5, 24[eax]                /* mm5 += p2 right             0 1 2 3 4 5 6 7    */
  285. psrlw  mm7, 4                       /* mm7 >>= 4                   0 1 2 3 4 5 6 7    */
  287. psubsw mm4, 48[ebx]                 /* mm4 -= vv[3] left           0 1 2 3 4 5 6 7    */
  288. add   ecx, stride                    /* ecx points at v[6*stride]   0 1 2 3 4 5 6 7    */     
  290. psubsw mm5, 56[ebx]                 /* mm5 -= vv[3] right           0 1 2 3 4 5 6 7    */
  292. movq   mm0, mm4                     /* mm0 = psum left             0 1 2 3 4 5 6 7    */ 
  294. paddsw mm0, 112[ebx]                /* mm0 += vv[7] left           0 1 2 3 4 5 6 7    */
  295. movq   mm1, mm5                     /* mm1 = psum right            0 1 2 3 4 5 6 7    */ 
  297. paddsw mm1, 120[ebx]                /* mm1 += vv[7] right          0 1 2 3 4 5 6 7    */
  298. psllw  mm0, 1                       /* mm0 <<= 1                   0 1 2 3 4 5 6 7    */
  300. paddsw mm0, 48[ebx]                 /* mm0 += vv[3] left           0 1 2 3 4 5 6 7    */
  301. psllw  mm1, 1                       /* mm1 <<= 1                   0 1 2 3 4 5 6 7    */
  303. paddsw mm1, 56[ebx]                 /* mm1 += vv[3] right          0 1 2 3 4 5 6 7    */
  304. packuswb mm6, mm7                   /* pack mm7, mm6 to mm6        0 1 2 3 4 5 6 7    */
  306. psubsw mm0, 64[ebx]                 /* mm0 -= vv[4] left           0 1 2 3 4 5 6 7    */
  308. psubsw mm1, 72[ebx]                 /* mm1 -= vv[4] right          0 1 2 3 4 5 6 7    */
  309. psrlw  mm0, 4                       /* mm0 >>= 4                   0 1 2 3 4 5 6 7    */
  311. movq   [ecx], mm6                     /* v[6*stride] = mm6           0 1 2 3 4 5 6 7    */
  313. paddsw mm4, 16[eax]                /* mm4 += p2 left               0 1 2 3 4 5 6 7    */
  315. paddsw mm5, 24[eax]                /* mm5 += p2 right              0 1 2 3 4 5 6 7    */
  316. add    ecx, stride                    /* ecx points at v[7*stride]   0 1 2 3 4 5 6 7    */     
  318. psubsw mm4, 64[ebx]                 /* mm4 -= vv[4] left            0 1 2 3 4 5 6 7    */
  319. psrlw  mm1, 4                       /* mm1 >>= 4                   0 1 2 3 4 5 6 7    */
  321. psubsw mm5, 72[ebx]                 /* mm5 -= vv[4] right           0 1 2 3 4 5 6 7    */
  323. movq   mm6, mm4                     /* mm6 = psum left             0 1 2 3 4 5 6 7    */ 
  325. paddsw mm6, 128[ebx]                /* mm6 += vv[8] left           0 1 2 3 4 5 6 7    */
  326. movq   mm7, mm5                     /* mm7 = psum right            0 1 2 3 4 5 6 7    */ 
  328. paddsw mm7, 136[ebx]                /* mm7 += vv[8] right          0 1 2 3 4 5 6 7    */
  329. psllw  mm6, 1                       /* mm6 <<= 1                   0 1 2 3 4 5 6 7    */
  331. paddsw mm6, 64[ebx]                 /* mm6 += vv[4] left           0 1 2 3 4 5 6 7    */
  332. psllw  mm7, 1                       /* mm7 <<= 1                   0 1 2 3 4 5 6 7    */
  334. paddsw mm7, 72[ebx]                 /* mm7 += vv[4] right          0 1 2 3 4 5 6 7    */
  335. packuswb mm0, mm1                   /* pack mm1, mm0 to mm0        0 1 2 3 4 5 6 7    */
  337. psubsw mm6, 80[ebx]                 /* mm6 -= vv[5] left           0 1 2 3 4 5 6 7    */
  339. psubsw mm7, 88[ebx]                 /* mm7 -= vv[5] right          0 1 2 3 4 5 6 7    */
  340. psrlw  mm6, 4                       /* mm6 >>= 4                   0 1 2 3 4 5 6 7    */
  342. movq   [ecx], mm0                   /* v[7*stride] = mm0           0 1 2 3 4 5 6 7    */
  343. psrlw  mm7, 4                       /* mm7 >>= 4                   0 1 2 3 4 5 6 7    */
  345. packuswb mm6, mm7                   /* pack mm7, mm6 to mm6        0 1 2 3 4 5 6 7    */
  347. add   ecx, stride                   /* ecx points at v[8*stride]   0 1 2 3 4 5 6 7    */     
  349. nop                                 /*                             0 1 2 3 4 5 6 7    */     
  351. movq   [ecx], mm6                   /* v[8*stride] = mm6           0 1 2 3 4 5 6 7    */
  354. pop ecx
  355. pop ebx
  356. pop eax
  359. };
  361. #ifdef PP_SELF_CHECK
  363. for (k=0; k<8; k++) { /* loop top->bottom */
  364. for (j=0; j<8; j++) { /* loop left->right */
  365. vv = &(v[(k+1)*stride + j]);
  366. if (*vv != selfcheck[j+8*k]) {
  367. printf("ERROR: problem with vertical LPF9 filter in row %dn", k+1);
  368. }
  369. }
  370. }
  371. #endif
  373. }