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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > 图片显示 > 查看源码
idct.cpp
资源名称：jpeg_decode　[点击查看]
上传用户：liulu_102
上传日期：2007-01-14
资源大小：193k
文件大小：38k
源码类别：

图片显示

开发平台：

Visual C++

idct.cpp：源码内容
  1. //********************************************************************************************
  2. //idctInt()和idctFloat()函数摘自“The Independent JPEG Group's JPEG software”文档，由下面
  3. //地址下载：ftp://ftp.simtel.net/pub/simtelnet/msdos/graphics/jpegsr6a.zip 并做了少量修改
  4. //domidct8x8llmW()函数摘自http://developer.intel.com/ 并做了少量修改
  5. //*******************************************************************************************
  7. #include "cjpeg.h"
  9. #define NO_ZERO_ROW_TEST
  10. #define BITS_IN_JSAMPLE 8
  12. #define MAXJSAMPLE 255
  13. #define CENTERJSAMPLE 128
  15. #define CONST_BITS  13
  16. #define PASS1_BITS  2
  18. #define FIX_0_298631336  ((INT32)  2446) /* FIX(0.298631336) */
  19. #define FIX_0_390180644  ((INT32)  3196) /* FIX(0.390180644) */
  20. #define FIX_0_541196100  ((INT32)  4433) /* FIX(0.541196100) */
  21. #define FIX_0_765366865  ((INT32)  6270) /* FIX(0.765366865) */
  22. #define FIX_0_899976223  ((INT32)  7373) /* FIX(0.899976223) */
  23. #define FIX_1_175875602  ((INT32)  9633) /* FIX(1.175875602) */
  24. #define FIX_1_501321110  ((INT32)  12299) /* FIX(1.501321110) */
  25. #define FIX_1_847759065  ((INT32)  15137) /* FIX(1.847759065) */
  26. #define FIX_1_961570560  ((INT32)  16069) /* FIX(1.961570560) */
  27. #define FIX_2_053119869  ((INT32)  16819) /* FIX(2.053119869) */
  28. #define FIX_2_562915447  ((INT32)  20995) /* FIX(2.562915447) */
  29. #define FIX_3_072711026  ((INT32)  25172) /* FIX(3.072711026) */
  31. #define ONE ((INT32) 1)
  32. #define RIGHT_SHIFT(x,shft)  
  33. ((shift_temp = (x)) < 0 ? 
  34.  (shift_temp >> (shft)) | ((~((INT32) 0)) << (32-(shft))) : 
  35.  (shift_temp >> (shft)))
  36. #define DESCALE(x,n)  RIGHT_SHIFT((x) + (ONE << ((n)-1)), n)
  37. #define RANGE_MASK  (MAXJSAMPLE * 4 + 3) /* 2 bits wider than legal samples */
  38. #define RANGE_LIMIT(i) 
  39. if((i) > MAXJSAMPLE)
  40. (i) = MAXJSAMPLE;
  41. else if((i) < 0)
  42. (i) = 0
  44. void CJpeg::idctInt(short* coef_block,BYTE *output_buf)
  45. { INT32 tmp0, tmp1, tmp2, tmp3;
  46. INT32 tmp10, tmp11, tmp12, tmp13;
  47. INT32 z1, z2, z3, z4, z5;
  48. int itp;
  49. short *inptr;
  50. int * wsptr;
  51. char *outptr;
  52. int ctr;
  53. int workspace[64];
  54. INT32 shift_temp;
  56. inptr = coef_block;
  57. wsptr = workspace;
  58. for (ctr = 8; ctr > 0; ctr--)
  59. { if ((inptr[8] | inptr[16] | inptr[24] | inptr[32] | inptr[40] | inptr[48] |
  60. inptr[56]) == 0) 
  61. { int dcval = inptr[0] << PASS1_BITS;
  62.     wsptr[0] = dcval;
  63.     wsptr[8] = dcval;
  64.     wsptr[16] = dcval;
  65.     wsptr[24] = dcval;
  66.     wsptr[32] = dcval;
  67.     wsptr[40] = dcval;
  68.     wsptr[48] = dcval;
  69.     wsptr[56] = dcval;
  70.     inptr++;
  71.     wsptr++;
  72.     continue;
  73.     }
  74.     z2 = inptr[16];
  75.     z3 = inptr[48];
  76.     z1 = (z2 + z3)*FIX_0_541196100;
  77.     tmp2 = z1 + z3*(- FIX_1_847759065);
  78.     tmp3 = z1 + z2*FIX_0_765366865;
  79.     z2 = inptr[0];
  80.     z3 = inptr[32];
  81.     tmp0 = (z2 + z3) << CONST_BITS;
  82.     tmp1 = (z2 - z3) << CONST_BITS;
  83.     tmp10 = tmp0 + tmp3;
  84.     tmp13 = tmp0 - tmp3;
  85.     tmp11 = tmp1 + tmp2;
  86.     tmp12 = tmp1 - tmp2;
  87.     tmp0 = inptr[56];
  88.     tmp1 = inptr[40];
  89.     tmp2 = inptr[24];
  90.     tmp3 = inptr[8];
  91.     z1 = tmp0 + tmp3;
  92.     z2 = tmp1 + tmp2;
  93.     z3 = tmp0 + tmp2;
  94.     z4 = tmp1 + tmp3;
  95.     z5 = (z3 + z4) * FIX_1_175875602; 
  96.     tmp0 = tmp0 * FIX_0_298631336; 
  97.     tmp1 = tmp1 * FIX_2_053119869;
  98.     tmp2 = tmp2 * FIX_3_072711026;
  99.     tmp3 = tmp3 * FIX_1_501321110;
  100.     z1 = z1 * (- FIX_0_899976223);
  101.     z2 = z2 * (- FIX_2_562915447);
  102.     z3 = z3 * (- FIX_1_961570560);
  103.     z4 = z4 * (- FIX_0_390180644);
  104.     z3 += z5;
  105.     z4 += z5;
  106.     tmp0 += z1 + z3;
  107.     tmp1 += z2 + z4;
  108.     tmp2 += z2 + z3;
  109.     tmp3 += z1 + z4;
  110.     wsptr[0] = (int) DESCALE(tmp10 + tmp3, CONST_BITS-PASS1_BITS);
  111.     wsptr[56] = (int) DESCALE(tmp10 - tmp3, CONST_BITS-PASS1_BITS);
  112.     wsptr[8] = (int) DESCALE(tmp11 + tmp2, CONST_BITS-PASS1_BITS);
  113.     wsptr[48] = (int) DESCALE(tmp11 - tmp2, CONST_BITS-PASS1_BITS);
  114.     wsptr[16] = (int) DESCALE(tmp12 + tmp1, CONST_BITS-PASS1_BITS);
  115.     wsptr[40] = (int) DESCALE(tmp12 - tmp1, CONST_BITS-PASS1_BITS);
  116.     wsptr[24] = (int) DESCALE(tmp13 + tmp0, CONST_BITS-PASS1_BITS);
  117.     wsptr[32] = (int) DESCALE(tmp13 - tmp0, CONST_BITS-PASS1_BITS);
  118.     inptr++;
  119.     wsptr++;
  120. }
  121. wsptr = workspace;
  122. for (ctr = 0; ctr < 8; ctr++) {
  123.     outptr = (char*)(output_buf + ctr*8);
  124. #ifndef NO_ZERO_ROW_TEST
  125.     if ((wsptr[1] | wsptr[2] | wsptr[3] | wsptr[4] | wsptr[5] | wsptr[6] | wsptr[7]) == 0)
  126. { char dcval;
  127. itp = (int) DESCALE((INT32) wsptr[0], PASS1_BITS+3)+128;
  128. if(itp > MAXJSAMPLE)
  129. itp = MAXJSAMPLE;
  130. else if(itp < 0)
  131. itp = 0;
  132. dcval = (char)itp;
  133.     outptr[0] = dcval;
  134. outptr[1] = dcval;
  135.     outptr[2] = dcval;
  136. outptr[3] = dcval;
  137.     outptr[4] = dcval;
  138. outptr[5] = dcval;
  139.     outptr[6] = dcval;
  140. outptr[7] = dcval;
  141.     wsptr += 8;
  142. continue;
  143. }
  144. #endif
  145.     z2 = (INT32) wsptr[2];
  146.     z3 = (INT32) wsptr[6];
  147.     z1 = (z2 + z3) * FIX_0_541196100;
  148.     tmp2 = z1 + z3 * (- FIX_1_847759065);
  149.     tmp3 = z1 + z2 * FIX_0_765366865;
  150.     tmp0 = ((INT32) wsptr[0] + (INT32) wsptr[4]) << CONST_BITS;
  151.     tmp1 = ((INT32) wsptr[0] - (INT32) wsptr[4]) << CONST_BITS;
  152.     tmp10 = tmp0 + tmp3;
  153.     tmp13 = tmp0 - tmp3;
  154.     tmp11 = tmp1 + tmp2;
  155.     tmp12 = tmp1 - tmp2;
  156.     tmp0 = (INT32) wsptr[7];
  157.     tmp1 = (INT32) wsptr[5];
  158.     tmp2 = (INT32) wsptr[3];
  159.     tmp3 = (INT32) wsptr[1];
  160.     z1 = tmp0 + tmp3;
  161.     z2 = tmp1 + tmp2;
  162.     z3 = tmp0 + tmp2;
  163.     z4 = tmp1 + tmp3;
  164.     z5 = (z3 + z4) * FIX_1_175875602;
  165.     tmp0 = tmp0 * FIX_0_298631336;
  166.     tmp1 = tmp1 * FIX_2_053119869;
  167.     tmp2 = tmp2 * FIX_3_072711026;
  168.     tmp3 = tmp3 * FIX_1_501321110;
  169.     z1 = z1 * (- FIX_0_899976223);
  170.     z2 = z2 * (- FIX_2_562915447);
  171.     z3 = z3 * (- FIX_1_961570560);
  172.     z4 = z4 * (- FIX_0_390180644);
  173.     z3 += z5;
  174.     z4 += z5;
  175.     tmp0 += z1 + z3;
  176.     tmp1 += z2 + z4;
  177.     tmp2 += z2 + z3;
  178.     tmp3 += z1 + z4;
  179.     itp = (int) DESCALE(tmp10 + tmp3,CONST_BITS+PASS1_BITS+3)+128;
  180. RANGE_LIMIT(itp);
  181. outptr[0] = (char)itp;
  182.     itp = (int) DESCALE(tmp10 - tmp3,CONST_BITS+PASS1_BITS+3)+128;
  183. RANGE_LIMIT(itp);
  184. outptr[7] = (char)itp;
  185.     itp = (int) DESCALE(tmp11 + tmp2,CONST_BITS+PASS1_BITS+3)+128;
  186. RANGE_LIMIT(itp);
  187. outptr[1] = (char)itp;
  188.     itp = (int) DESCALE(tmp11 - tmp2,CONST_BITS+PASS1_BITS+3)+128;
  189. RANGE_LIMIT(itp);
  190. outptr[6] = (char)itp;
  191.     itp = (int) DESCALE(tmp12 + tmp1,CONST_BITS+PASS1_BITS+3)+128;
  192. RANGE_LIMIT(itp);
  193. outptr[2] = (char)itp;
  194.     itp = (int) DESCALE(tmp12 - tmp1,CONST_BITS+PASS1_BITS+3)+128;
  195. RANGE_LIMIT(itp);
  196. outptr[5] = (char)itp;
  197.     itp = (int) DESCALE(tmp13 + tmp0,CONST_BITS+PASS1_BITS+3)+128;
  198. RANGE_LIMIT(itp);
  199. outptr[3] = (char)itp;
  200.     itp = (int) DESCALE(tmp13 - tmp0,CONST_BITS+PASS1_BITS+3)+128;
  201. RANGE_LIMIT(itp);
  202. outptr[4] = (char)itp;
  203.     wsptr += 8;
  204.   }
  205. }
  207. void CJpeg::idctFloat(short* coef_block,BYTE *output_buf)
  208. { float tmp0, tmp1, tmp2, tmp3;
  209. float tmp10, tmp11, tmp12, tmp13;
  210. float z1, z2, z3, z4, z5;
  211. int itp;
  212. short *inptr;
  213. float * wsptr;
  214. char *outptr;
  215. int ctr;
  216. float workspace[64];
  217. INT32 shift_temp;
  219. inptr = coef_block;
  220. wsptr = workspace;
  221. for (ctr = 8; ctr > 0; ctr--)
  222. { if ((inptr[8] | inptr[16] | inptr[24] | inptr[32] | inptr[40] | inptr[48] |
  223. inptr[56]) == 0) 
  224. { float dcval = inptr[0];
  225.     wsptr[0] = dcval;
  226.     wsptr[8] = dcval;
  227.     wsptr[16] = dcval;
  228.     wsptr[24] = dcval;
  229.     wsptr[32] = dcval;
  230.     wsptr[40] = dcval;
  231.     wsptr[48] = dcval;
  232.     wsptr[56] = dcval;
  233.     inptr++;
  234.     wsptr++;
  235.     continue;
  236.     }
  237.     z2 = inptr[16];
  238.     z3 = inptr[48];
  239.     z1 = (z2 + z3)*(float)0.541196100;
  240.     tmp2 = z1 + z3*(float)(-1.847759065);
  241.     tmp3 = z1 + z2*(float)0.765366865;
  242.     z2 = inptr[0];
  243.     z3 = inptr[32];
  244.     tmp0 = (z2 + z3);
  245.     tmp1 = (z2 - z3);
  246.     tmp10 = tmp0 + tmp3;
  247.     tmp13 = tmp0 - tmp3;
  248.     tmp11 = tmp1 + tmp2;
  249.     tmp12 = tmp1 - tmp2;
  250.     tmp0 = inptr[56];
  251.     tmp1 = inptr[40];
  252.     tmp2 = inptr[24];
  253.     tmp3 = inptr[8];
  254.     z1 = tmp0 + tmp3;
  255.     z2 = tmp1 + tmp2;
  256.     z3 = tmp0 + tmp2;
  257.     z4 = tmp1 + tmp3;
  258.     z5 = (z3 + z4) * (float)1.175875602; 
  259.     tmp0 = tmp0 * (float)0.298631336; 
  260.     tmp1 = tmp1 * (float)2.053119869;
  261.     tmp2 = tmp2 * (float)3.072711026;
  262.     tmp3 = tmp3 * (float)1.501321110;
  263.     z1 = z1 * (float)(-0.899976223);
  264.     z2 = z2 * (float)(-2.562915447);
  265.     z3 = z3 * (float)(-1.961570560);
  266.     z4 = z4 * (float)(-0.390180644);
  267.     z3 += z5;
  268.     z4 += z5;
  269.     tmp0 += z1 + z3;
  270.     tmp1 += z2 + z4;
  271.     tmp2 += z2 + z3;
  272.     tmp3 += z1 + z4;
  273.     wsptr[0] = tmp10 + tmp3;
  274.     wsptr[56] = tmp10 - tmp3;
  275.     wsptr[8] = tmp11 + tmp2;
  276.     wsptr[48] = tmp11 - tmp2;
  277.     wsptr[16] = tmp12 + tmp1;
  278.     wsptr[40] = tmp12 - tmp1;
  279.     wsptr[24] = tmp13 + tmp0;
  280.     wsptr[32] = tmp13 - tmp0;
  281.     inptr++;
  282.     wsptr++;
  283. }
  284. wsptr = workspace;
  285. for (ctr = 0; ctr < 8; ctr++) {
  286.     outptr = (char*)(output_buf + ctr*8);
  287.     z2 = wsptr[2];
  288.     z3 = wsptr[6];
  289.     z1 = (z2 + z3) * (float)0.541196100;
  290.     tmp2 = z1 + z3 * (float)(-1.847759065);
  291.     tmp3 = z1 + z2 * (float)0.765366865;
  292.     tmp0 = wsptr[0] + wsptr[4];
  293.     tmp1 = wsptr[0] - wsptr[4];
  294.     tmp10 = tmp0 + tmp3;
  295.     tmp13 = tmp0 - tmp3;
  296.     tmp11 = tmp1 + tmp2;
  297.     tmp12 = tmp1 - tmp2;
  298.     tmp0 = wsptr[7];
  299.     tmp1 = wsptr[5];
  300.     tmp2 = wsptr[3];
  301.     tmp3 = wsptr[1];
  302.     z1 = tmp0 + tmp3;
  303.     z2 = tmp1 + tmp2;
  304.     z3 = tmp0 + tmp2;
  305.     z4 = tmp1 + tmp3;
  306.     z5 = (z3 + z4) * (float)1.175875602;
  307.     tmp0 = tmp0 * (float)0.298631336;
  308.     tmp1 = tmp1 * (float)2.053119869;
  309.     tmp2 = tmp2 * (float)3.072711026;
  310.     tmp3 = tmp3 * (float)1.501321110;
  311.     z1 = z1 * (float)(-0.899976223);
  312.     z2 = z2 * (float)(-2.562915447);
  313.     z3 = z3 * (float)(-1.961570560);
  314.     z4 = z4 * (float)(-0.390180644);
  315.     z3 += z5;
  316.     z4 += z5;
  317.     tmp0 += z1 + z3;
  318.     tmp1 += z2 + z4;
  319.     tmp2 += z2 + z3;
  320.     tmp3 += z1 + z4;
  321.     itp = (int) DESCALE((INT32)(tmp10 + tmp3),3)+128;
  322. RANGE_LIMIT(itp);
  323. outptr[0] = (char)itp;
  324.     itp = (int) DESCALE((INT32)(tmp10 - tmp3),3)+128;
  325. RANGE_LIMIT(itp);
  326. outptr[7] = (char)itp;
  327.     itp = (int) DESCALE((INT32)(tmp11 + tmp2),3)+128;
  328. RANGE_LIMIT(itp);
  329. outptr[1] = (char)itp;
  330.     itp = (int) DESCALE((INT32)(tmp11 - tmp2),3)+128;
  331. RANGE_LIMIT(itp);
  332. outptr[6] = (char)itp;
  333.     itp = (int) DESCALE((INT32)(tmp12 + tmp1),3)+128;
  334. RANGE_LIMIT(itp);
  335. outptr[2] = (char)itp;
  336.     itp = (int) DESCALE((INT32)(tmp12 - tmp1),3)+128;
  337. RANGE_LIMIT(itp);
  338. outptr[5] = (char)itp;
  339.     itp = (int) DESCALE((INT32)(tmp13 + tmp0),3)+128;
  340. RANGE_LIMIT(itp);
  341. outptr[3] = (char)itp;
  342.     itp = (int) DESCALE((INT32)(tmp13 - tmp0),3)+128;
  343. RANGE_LIMIT(itp);
  344. outptr[4] = (char)itp;
  345.     wsptr += 8;
  346.   }
  347. }
  349. void CJpeg::idctMMX(short* coef_block,BYTE *output_buf)
  350. { int workspace[64+4];
  351. domidct8x8llmW(coef_block, workspace, output_buf);
  352. }
  354. static __int64 fix_029_n089n196 = 0x098ea46e098ea46e;
  355. static __int64 fix_n196_n089 = 0xc13be333c13be333;
  356. static __int64 fix_205_n256n039 = 0x41b3a18141b3a181;
  357. static __int64 fix_n039_n256 = 0xf384adfdf384adfd;
  358. static __int64 fix_307n256_n196 = 0x1051c13b1051c13b;
  359. static __int64 fix_n256_n196 = 0xadfdc13badfdc13b;
  360. static __int64 fix_150_n089n039 = 0x300bd6b7300bd6b7;
  361. static __int64 fix_n039_n089 = 0xf384e333f384e333;
  362. static __int64 fix_117_117 = 0x25a125a125a125a1;
  363. static __int64 fix_054_054p076 = 0x115129cf115129cf;
  364. static __int64 fix_054n184_054 = 0xd6301151d6301151;
  365. static __int64 fix_054n184  = 0xd630d630d630d630;
  366. static __int64 fix_054 = 0x1151115111511151;
  367. static __int64 fix_054p076 = 0x29cf29cf29cf29cf;
  368. static __int64 fix_n196p307n256 = 0xd18cd18cd18cd18c;
  369. static __int64 fix_n089n039p150 = 0x06c206c206c206c2;
  370. static __int64 fix_n256 = 0xadfdadfdadfdadfd;
  371. static __int64 fix_n039 = 0xf384f384f384f384;
  372. static __int64 fix_n256n039p205 = 0xe334e334e334e334;
  373. static __int64 fix_n196 = 0xc13bc13bc13bc13b;
  374. static __int64 fix_n089 = 0xe333e333e333e333;
  375. static __int64 fixn089n196p029 = 0xadfcadfcadfcadfc;
  376. static  __int64 const_0x2xx8 = 0x0000010000000100;
  377. static  __int64 const_0x0808 = 0x0808080808080808;
  378. void CJpeg::domidct8x8llmW(short *inptr,int *wsptr,BYTE* outptr)
  379. {
  380. __asm{
  381. mov ebx, inptr
  382. mov esi, wsptr
  383. add esi, 0x07 ;align wsptr to qword
  384. and esi, 0xfffffff8 ;align wsptr to qword
  385. mov eax, esi
  386. /* Pass 1. */
  387. movq mm0, [ebx + 8*4] ;p1(1,0)
  388.     movq mm1, [ebx + 8*12] ;p1(2,0)
  389. movq mm6, [ebx + 8*0] ;p1(5,0)
  390. movq mm2, mm0 ;p1(3,0)
  391. movq mm7, [ebx + 8*8] ;p1(6,0)
  392. punpcklwd mm0, mm1 ;p1(3,1)
  393. movq mm4, mm0 ;p1(3,2)
  394. punpckhwd mm2, mm1 ;p1(3,4)
  395. pmaddwd mm0, fix_054n184_054;p1(3,3)
  396. movq mm5, mm2 ;p1(3,5)
  397. pmaddwd mm2, fix_054n184_054;p1(3,6)
  398. pxor mm1, mm1 ;p1(7,0)
  399. pmaddwd mm4, fix_054_054p076;p1(4,0)
  400. punpcklwd   mm1, mm6 ;p1(7,1)
  401. pmaddwd mm5, fix_054_054p076;p1(4,1)
  402. psrad mm1, 3 ;p1(7,2)
  403. pxor mm3, mm3 ;p1(7,3)
  404. punpcklwd mm3, mm7 ;p1(7,4)
  405. psrad mm3, 3 ;p1(7,5)
  406. paddd mm1, mm3 ;p1(7,6)
  407. movq mm3, mm1 ;p1(7,7)
  408. paddd mm1, mm4 ;p1(7,8)
  409. psubd mm3, mm4 ;p1(7,9)
  410. movq [esi + 8*16], mm1 ;p1(7,10)
  411. pxor mm4, mm4 ;p1(7,12)
  412. movq [esi + 8*22], mm3 ;p1(7,11)
  413. punpckhwd mm4, mm6 ;p1(7,13)
  414. psrad mm4, 3 ;p1(7,14)
  415. pxor mm1, mm1 ;p1(7,15)
  416. punpckhwd mm1, mm7 ;p1(7,16)
  417. psrad mm1, 3 ;p1(7,17)
  418. paddd mm4, mm1 ;p1(7,18)
  419. movq mm3, mm4 ;p1(7,19)
  420. pxor mm1, mm1 ;p1(8,0)
  421. paddd mm3, mm5 ;p1(7,20)
  422. punpcklwd mm1, mm6 ;p1(8,1)
  423. psubd mm4, mm5 ;p1(7,21)
  424. psrad mm1, 3 ;p1(8,2)
  425. movq [esi + 8*17], mm3 ;p1(7,22)
  426. pxor mm5, mm5 ;p1(8,3)
  427. movq [esi + 8*23], mm4 ;p1(7,23)
  428. punpcklwd mm5, mm7 ;p1(8,4)
  429. psrad mm5, 3 ;p1(8,5)
  430. pxor mm4, mm4 ;p1(8,12)
  431. psubd mm1, mm5 ;p1(8,6)
  432. punpckhwd mm4, mm6 ;p1(8,13)
  433. movq mm3, mm1 ;p1(8,7)
  434. psrad mm4, 3 ;p1(8,14)
  435. paddd mm1, mm0 ;p1(8,8)
  436. pxor mm5, mm5 ;p1(8,15)
  437. psubd mm3, mm0 ;p1(8,9)
  438. movq mm0, [ebx + 8*14] ;p1(9,0)
  439. punpckhwd mm5, mm7 ;p1(8,16)
  440. movq [esi + 8*18], mm1 ;p1(8,10)
  441. psrad mm5, 3 ;p1(8,17)
  442. movq [esi + 8*20], mm3 ;p1(8,11)
  443. psubd mm4, mm5 ;p1(8,18)
  444. movq mm3, mm4 ;p1(8,19)
  445. movq mm1, [ebx + 8*6] ;p1(10,0)
  446. paddd mm3, mm2 ;p1(8,20)
  447. psubd mm4, mm2 ;p1(8,21)
  448. movq mm5, mm0 ;p1(11,1)
  449. movq [esi + 8*21], mm4 ;p1(8,23)
  450. movq [esi + 8*19], mm3 ;p1(8,22)
  451. movq mm4, mm0 ;p1(11,0)
  452. punpcklwd mm4, mm1 ;p1(11,2)
  453. movq mm2, [ebx + 8*10] ;p1(12,0)
  454. punpckhwd mm5, mm1 ;p1(11,4)
  455. movq mm3, [ebx + 8*2] ;p1(13,0)
  456. movq mm6, mm2 ;p1(14,0)
  457. pmaddwd mm4, fix_117_117 ;p1(11,3)
  458. movq mm7, mm2 ;p1(14,1)
  459. pmaddwd mm5, fix_117_117 ;p1(11,5)
  460. punpcklwd mm6, mm3 ;p1(14,2)
  461. pmaddwd mm6, fix_117_117 ;p1(14,3)
  462. punpckhwd mm7, mm3 ;p1(14,4)
  463. pmaddwd mm7, fix_117_117 ;p1(14,5)
  464. paddd mm4, mm6 ;p1(15,0)
  465. paddd mm5, mm7 ;p1(15,1)
  466.     movq [esi+8*24], mm4 ;p1(15,2)
  467. movq [esi+8*25], mm5 ;p1(15,3)
  468. movq mm6, mm0 ;p1(16,0)
  469. movq mm7, mm3 ;p1(16,3)
  470. punpcklwd mm6, mm2 ;p1(16,1)
  471. punpcklwd mm7, mm3 ;p1(16,4)
  472. pmaddwd mm6, fix_n039_n089 ;p1(16,2)
  473. pmaddwd mm7, fix_150_n089n039 ;p1(16,5)
  474. movq mm4, mm0 ;p1(16,12)
  475. paddd mm6, [esi+8*24] ;p1(16,6)
  476. punpckhwd mm4, mm2 ;p1(16,13)
  477. paddd mm6, mm7 ;p1(16,7)
  478. pmaddwd mm4, fix_n039_n089 ;p1(16,14)
  479. movq mm7, mm6 ;p1(16,8)
  480. paddd mm4, [esi+8*25] ;p1(16,18)
  481. movq mm5, mm3 ;p1(16,15)
  482. paddd mm6, [esi + 8*16] ;p1(16,9)
  483. punpckhwd mm5, mm3 ;p1(16,16)
  484. paddd mm6, const_0x2xx8 ;p1(16,10)
  485. psrad mm6, 9 ;p1(16,11)
  486. pmaddwd mm5, fix_150_n089n039 ;p1(16,17)
  487. paddd mm4, mm5 ;p1(16,19)
  488. movq mm5, mm4 ;p1(16,20)
  489. paddd mm4, [esi + 8*17] ;p1(16,21)
  490. paddd mm4, const_0x2xx8 ;p1(16,22)
  491. psrad mm4, 9 ;p1(16,23)
  492. packssdw mm6, mm4 ;p1(16,24)
  493. movq [esi + 8*0], mm6 ;p1(16,25)
  494. movq mm4, [esi + 8*16] ;p1(16,26)
  495. psubd mm4, mm7 ;p1(16,27)
  496. movq mm6, [esi + 8*17] ;p1(16,30)
  497. paddd mm4, const_0x2xx8 ;p1(16,28)
  498. movq mm7, mm1 ;p1(17,3)
  499. psrad mm4, 9 ;p1(16,29)
  500. psubd mm6, mm5 ;p1(16,31)
  501. paddd mm6, const_0x2xx8 ;p1(16,32)
  502. punpcklwd mm7, mm1 ;p1(17,4)
  503. pmaddwd mm7, fix_307n256_n196 ;p1(17,5)
  504. psrad mm6, 9 ;p1(16,33)
  505. packssdw mm4, mm6 ;p1(16,34)
  506. movq [esi + 8*14], mm4 ;p1(16,35)
  507. movq mm6, mm0 ;p1(17,0)
  508. movq mm4, mm0 ;p1(17,12)
  509. punpcklwd mm6, mm2 ;p1(17,1)
  510. punpckhwd mm4, mm2 ;p1(17,13)
  511. pmaddwd mm6, fix_n256_n196 ;p1(17,2)
  512. movq mm5, mm1 ;p1(17,15)
  513. paddd mm6, [esi+8*24] ;p1(17,6)
  514. punpckhwd mm5, mm1 ;p1(17,16)
  515. paddd mm6, mm7 ;p1(17,7)
  516. pmaddwd mm4, fix_n256_n196 ;p1(17,14)
  517. movq mm7, mm6 ;p1(17,8)
  518. pmaddwd mm5, fix_307n256_n196 ;p1(17,17)
  519. paddd mm6, [esi + 8*18] ;p1(17,9)
  520. paddd mm6, const_0x2xx8 ;p1(17,10)
  521. psrad mm6, 9 ;p1(17,11)
  522. paddd mm4, [esi+8*25] ;p1(17,18)
  523. paddd mm4, mm5 ;p1(17,19)
  524. movq mm5, mm4 ;p1(17,20)
  525. paddd mm4, [esi + 8*19] ;p1(17,21)
  526. paddd mm4, const_0x2xx8 ;p1(17,22)
  527. psrad mm4, 9 ;p1(17,23)
  528. packssdw mm6, mm4 ;p1(17,24)
  529. movq [esi + 8*2], mm6 ;p1(17,25)
  530. movq mm4, [esi + 8*18] ;p1(17,26)
  531. movq mm6, [esi + 8*19] ;p1(17,30)
  532. psubd mm4, mm7 ;p1(17,27)
  533. paddd mm4, const_0x2xx8 ;p1(17,28)
  534. psubd mm6, mm5 ;p1(17,31)
  535. psrad mm4, 9 ;p1(17,29)
  536. paddd mm6, const_0x2xx8 ;p1(17,32)
  537. psrad mm6, 9 ;p1(17,33)
  538. movq mm7, mm2 ;p1(18,3)
  539. packssdw mm4, mm6 ;p1(17,34)
  540. movq [esi + 8*12], mm4 ;p1(17,35)
  541. movq mm6, mm1 ;p1(18,0)
  542. punpcklwd mm7, mm2 ;p1(18,4)
  543. punpcklwd mm6, mm3 ;p1(18,1)
  544. pmaddwd mm7, fix_205_n256n039 ;p1(18,5)
  545. pmaddwd mm6, fix_n039_n256 ;p1(18,2)
  546. movq mm4, mm1 ;p1(18,12)
  547. paddd mm6, [esi+8*24] ;p1(18,6)
  548. punpckhwd mm4, mm3 ;p1(18,13)
  549. paddd mm6, mm7 ;p1(18,7)
  550. pmaddwd mm4, fix_n039_n256 ;p1(18,14)
  551. movq mm7, mm6 ;p1(18,8)
  552. movq mm5, mm2 ;p1(18,15)
  553. paddd mm6, [esi + 8*20] ;p1(18,9)
  554. punpckhwd mm5, mm2 ;p1(18,16)
  555. paddd mm6, const_0x2xx8 ;p1(18,10)
  556. psrad mm6, 9 ;p1(18,11)
  557. pmaddwd mm5, fix_205_n256n039 ;p1(18,17)
  558. paddd mm4, [esi+8*25] ;p1(18,18)
  559. paddd mm4, mm5 ;p1(18,19)
  560. movq mm5, mm4 ;p1(18,20)
  561. paddd mm4, [esi + 8*21] ;p1(18,21)
  562. paddd mm4, const_0x2xx8 ;p1(18,22)
  563. psrad mm4, 9 ;p1(18,23)
  564. packssdw mm6, mm4 ;p1(18,24)
  565. movq [esi + 8*4], mm6 ;p1(18,25)
  566. movq mm4, [esi + 8*20] ;p1(18,26)
  567. psubd mm4, mm7 ;p1(18,27)
  568. paddd mm4, const_0x2xx8 ;p1(18,28)
  569. movq mm7, mm0 ;p1(19,3)
  570. psrad mm4, 9 ;p1(18,29)
  571. movq mm6, [esi + 8*21] ;p1(18,30)
  572. psubd mm6, mm5 ;p1(18,31)
  573. punpcklwd mm7, mm0 ;p1(19,4)
  574. paddd mm6, const_0x2xx8 ;p1(18,32)
  575. psrad mm6, 9 ;p1(18,33)
  576. pmaddwd mm7, fix_029_n089n196 ;p1(19,5)
  577. packssdw mm4, mm6 ;p1(18,34)
  578. movq [esi + 8*10], mm4 ;p1(18,35)
  579. movq mm6, mm3 ;p1(19,0)
  580. punpcklwd mm6, mm1 ;p1(19,1)
  581. movq mm5, mm0 ;p1(19,15)
  582. pmaddwd mm6, fix_n196_n089 ;p1(19,2)
  583. punpckhwd mm5, mm0 ;p1(19,16)
  584. paddd mm6, [esi+8*24] ;p1(19,6)
  585. movq mm4, mm3 ;p1(19,12)
  586. paddd mm6, mm7 ;p1(19,7)
  587. punpckhwd mm4, mm1 ;p1(19,13)
  588. movq mm7, mm6 ;p1(19,8)
  589. pmaddwd mm4, fix_n196_n089   ;p1(19,14)
  590. paddd mm6, [esi + 8*22] ;p1(19,9)
  591. pmaddwd mm5, fix_029_n089n196 ;p1(19,17)
  592. paddd mm6, const_0x2xx8 ;p1(19,10)
  593. psrad mm6, 9 ;p1(19,11)
  594. paddd mm4, [esi+8*25] ;p1(19,18)
  595. paddd mm4, mm5 ;p1(19,19)
  596. movq mm5, mm4 ;p1(19,20)
  597. paddd mm4, [esi + 8*23] ;p1(19,21)
  598. paddd mm4, const_0x2xx8 ;p1(19,22)
  599. psrad mm4, 9 ;p1(19,23)
  600. packssdw mm6, mm4 ;p1(19,24)
  601. movq [esi + 8*6], mm6 ;p1(19,25)
  602. movq mm4, [esi + 8*22] ;p1(19,26)
  603. psubd mm4, mm7 ;p1(19,27)
  604. movq mm6, [esi + 8*23] ;p1(19,30)
  605. paddd mm4, const_0x2xx8 ;p1(19,28)
  606. psubd mm6, mm5 ;p1(19,31)
  607. psrad mm4, 9 ;p1(19,29)
  608. paddd mm6, const_0x2xx8 ;p1(19,32)
  609. psrad mm6, 9 ;p1(19,33)
  610. packssdw mm4, mm6 ;p1(19,34)
  611. movq [esi + 8*8], mm4 ;p1(19,35)
  612. //************************************************************  
  613. // add edi, 8
  614. add ebx, 8
  615. add esi, 8
  616. //************************************************************  
  617. /* Pass 1. */
  618. movq mm0, [ebx + 8*4] ;p1(1,0)
  619.     movq mm1, [ebx + 8*12] ;p1(2,0)
  620. movq mm6, [ebx + 8*0] ;p1(5,0)
  621. movq mm2, mm0 ;p1(3,0)
  622. movq mm7, [ebx + 8*8] ;p1(6,0)
  623. punpcklwd mm0, mm1 ;p1(3,1)
  624. movq mm4, mm0 ;p1(3,2)
  625. punpckhwd mm2, mm1 ;p1(3,4)
  626. pmaddwd mm0, fix_054n184_054 ;p1(3,3)
  627. movq mm5, mm2 ;p1(3,5)
  628. pmaddwd mm2, fix_054n184_054 ;p1(3,6)
  629. pxor mm1, mm1 ;p1(7,0)
  630. pmaddwd mm4, fix_054_054p076 ;p1(4,0)
  631. punpcklwd   mm1, mm6 ;p1(7,1)
  632. pmaddwd mm5, fix_054_054p076 ;p1(4,1)
  633. psrad mm1, 3 ;p1(7,2)
  634. pxor mm3, mm3 ;p1(7,3)
  635. punpcklwd mm3, mm7 ;p1(7,4)
  636. psrad mm3, 3 ;p1(7,5)
  637. paddd mm1, mm3 ;p1(7,6)
  638. movq mm3, mm1 ;p1(7,7)
  639. paddd mm1, mm4 ;p1(7,8)
  640. psubd mm3, mm4 ;p1(7,9)
  641. movq [esi + 8*16], mm1 ;p1(7,10)
  642. pxor mm4, mm4 ;p1(7,12)
  643. movq [esi + 8*22], mm3 ;p1(7,11)
  644. punpckhwd mm4, mm6 ;p1(7,13)
  645. psrad mm4, 3 ;p1(7,14)
  646. pxor mm1, mm1 ;p1(7,15)
  647. punpckhwd mm1, mm7 ;p1(7,16)
  648. psrad mm1, 3 ;p1(7,17)
  649. paddd mm4, mm1 ;p1(7,18)
  650. movq mm3, mm4 ;p1(7,19)
  651. pxor mm1, mm1 ;p1(8,0)
  652. paddd mm3, mm5 ;p1(7,20)
  653. punpcklwd mm1, mm6 ;p1(8,1)
  654. psubd mm4, mm5 ;p1(7,21)
  655. psrad mm1, 3 ;p1(8,2)
  656. movq [esi + 8*17], mm3 ;p1(7,22)
  657. pxor mm5, mm5 ;p1(8,3)
  658. movq [esi + 8*23], mm4 ;p1(7,23)
  659. punpcklwd mm5, mm7 ;p1(8,4)
  660. psrad mm5, 3 ;p1(8,5)
  661. pxor mm4, mm4 ;p1(8,12)
  662. psubd mm1, mm5 ;p1(8,6)
  663. punpckhwd mm4, mm6 ;p1(8,13)
  664. movq mm3, mm1 ;p1(8,7)
  665. psrad mm4, 3 ;p1(8,14)
  666. paddd mm1, mm0 ;p1(8,8)
  667. pxor mm5, mm5 ;p1(8,15)
  668. psubd mm3, mm0 ;p1(8,9)
  669. movq mm0, [ebx + 8*14] ;p1(9,0)
  670. punpckhwd mm5, mm7 ;p1(8,16)
  671. movq [esi + 8*18], mm1 ;p1(8,10)
  672. psrad mm5, 3 ;p1(8,17)
  673. movq [esi + 8*20], mm3 ;p1(8,11)
  674. psubd mm4, mm5 ;p1(8,18)
  675. movq mm3, mm4 ;p1(8,19)
  676. movq mm1, [ebx + 8*6] ;p1(10,0)
  677. paddd mm3, mm2 ;p1(8,20)
  678. psubd mm4, mm2 ;p1(8,21)
  679. movq mm5, mm0 ;p1(11,1)
  680. movq [esi + 8*21], mm4 ;p1(8,23)
  681. movq [esi + 8*19], mm3 ;p1(8,22)
  682. movq mm4, mm0 ;p1(11,0)
  683. punpcklwd mm4, mm1 ;p1(11,2)
  684. movq mm2, [ebx + 8*10] ;p1(12,0)
  685. punpckhwd mm5, mm1 ;p1(11,4)
  686. movq mm3, [ebx + 8*2] ;p1(13,0)
  687. movq mm6, mm2 ;p1(14,0)
  688. pmaddwd mm4, fix_117_117 ;p1(11,3)
  689. movq mm7, mm2 ;p1(14,1)
  690. pmaddwd mm5, fix_117_117 ;p1(11,5)
  691. punpcklwd mm6, mm3 ;p1(14,2)
  692. pmaddwd mm6, fix_117_117 ;p1(14,3)
  693. punpckhwd mm7, mm3 ;p1(14,4)
  694. pmaddwd mm7, fix_117_117 ;p1(14,5)
  695. paddd mm4, mm6 ;p1(15,0)
  696. paddd mm5, mm7 ;p1(15,1)
  697.     movq [esi+8*24], mm4 ;p1(15,2)
  698. movq [esi+8*25], mm5 ;p1(15,3)
  699. movq mm6, mm0 ;p1(16,0)
  700. movq mm7, mm3 ;p1(16,3)
  701. punpcklwd mm6, mm2 ;p1(16,1)
  702. punpcklwd mm7, mm3 ;p1(16,4)
  703. pmaddwd mm6, fix_n039_n089 ;p1(16,2)
  704. pmaddwd mm7, fix_150_n089n039 ;p1(16,5)
  705. movq mm4, mm0 ;p1(16,12)
  706. paddd mm6, [esi+8*24] ;p1(16,6)
  707. punpckhwd mm4, mm2 ;p1(16,13)
  708. paddd mm6, mm7 ;p1(16,7)
  709. pmaddwd mm4, fix_n039_n089 ;p1(16,14)
  710. movq mm7, mm6 ;p1(16,8)
  711. paddd mm4, [esi+8*25] ;p1(16,18)
  712. movq mm5, mm3 ;p1(16,15)
  713. paddd mm6, [esi + 8*16] ;p1(16,9)
  714. punpckhwd mm5, mm3 ;p1(16,16)
  715. paddd mm6, const_0x2xx8 ;p1(16,10)
  716. psrad mm6, 9 ;p1(16,11)
  717. pmaddwd mm5, fix_150_n089n039 ;p1(16,17)
  718. paddd mm4, mm5 ;p1(16,19)
  719. movq mm5, mm4 ;p1(16,20)
  720. paddd mm4, [esi + 8*17] ;p1(16,21)
  721. paddd mm4, const_0x2xx8 ;p1(16,22)
  722. psrad mm4, 9 ;p1(16,23)
  723. packssdw mm6, mm4 ;p1(16,24)
  724. movq [esi + 8*0], mm6 ;p1(16,25)
  725. movq mm4, [esi + 8*16] ;p1(16,26)
  726. psubd mm4, mm7 ;p1(16,27)
  727. movq mm6, [esi + 8*17] ;p1(16,30)
  728. paddd mm4, const_0x2xx8 ;p1(16,28)
  729. movq mm7, mm1 ;p1(17,3)
  730. psrad mm4, 9 ;p1(16,29)
  731. psubd mm6, mm5 ;p1(16,31)
  732. paddd mm6, const_0x2xx8 ;p1(16,32)
  733. punpcklwd mm7, mm1 ;p1(17,4)
  734. pmaddwd mm7, fix_307n256_n196 ;p1(17,5)
  735. psrad mm6, 9 ;p1(16,33)
  736. packssdw mm4, mm6 ;p1(16,34)
  737. movq [esi + 8*14], mm4 ;p1(16,35)
  738. movq mm6, mm0 ;p1(17,0)
  739. movq mm4, mm0 ;p1(17,12)
  740. punpcklwd mm6, mm2 ;p1(17,1)
  741. punpckhwd mm4, mm2 ;p1(17,13)
  742. pmaddwd mm6, fix_n256_n196 ;p1(17,2)
  743. movq mm5, mm1 ;p1(17,15)
  744. paddd mm6, [esi+8*24] ;p1(17,6)
  745. punpckhwd mm5, mm1 ;p1(17,16)
  746. paddd mm6, mm7 ;p1(17,7)
  747. pmaddwd mm4, fix_n256_n196 ;p1(17,14)
  748. movq mm7, mm6 ;p1(17,8)
  749. pmaddwd mm5, fix_307n256_n196 ;p1(17,17)
  750. paddd mm6, [esi + 8*18] ;p1(17,9)
  751. paddd mm6, const_0x2xx8 ;p1(17,10)
  752. psrad mm6, 9 ;p1(17,11)
  753. paddd mm4, [esi+8*25] ;p1(17,18)
  754. paddd mm4, mm5 ;p1(17,19)
  755. movq mm5, mm4 ;p1(17,20)
  756. paddd mm4, [esi + 8*19] ;p1(17,21)
  757. paddd mm4, const_0x2xx8 ;p1(17,22)
  758. psrad mm4, 9 ;p1(17,23)
  759. packssdw mm6, mm4 ;p1(17,24)
  760. movq [esi + 8*2], mm6 ;p1(17,25)
  761. movq mm4, [esi + 8*18] ;p1(17,26)
  762. movq mm6, [esi + 8*19] ;p1(17,30)
  763. psubd mm4, mm7 ;p1(17,27)
  764. paddd mm4, const_0x2xx8 ;p1(17,28)
  765. psubd mm6, mm5 ;p1(17,31)
  766. psrad mm4, 9 ;p1(17,29)
  767. paddd mm6, const_0x2xx8 ;p1(17,32)
  768. psrad mm6, 9 ;p1(17,33)
  769. movq mm7, mm2 ;p1(18,3)
  770. packssdw mm4, mm6 ;p1(17,34)
  771. movq [esi + 8*12], mm4 ;p1(17,35)
  772. movq mm6, mm1 ;p1(18,0)
  773. punpcklwd mm7, mm2 ;p1(18,4)
  774. punpcklwd mm6, mm3 ;p1(18,1)
  775. pmaddwd mm7, fix_205_n256n039 ;p1(18,5)
  776. pmaddwd mm6, fix_n039_n256 ;p1(18,2)
  777. movq mm4, mm1 ;p1(18,12)
  778. paddd mm6, [esi+8*24] ;p1(18,6)
  779. punpckhwd mm4, mm3 ;p1(18,13)
  780. paddd mm6, mm7 ;p1(18,7)
  781. pmaddwd mm4, fix_n039_n256 ;p1(18,14)
  782. movq mm7, mm6 ;p1(18,8)
  783. movq mm5, mm2 ;p1(18,15)
  784. paddd mm6, [esi + 8*20] ;p1(18,9)
  785. punpckhwd mm5, mm2 ;p1(18,16)
  786. paddd mm6, const_0x2xx8 ;p1(18,10)
  787. psrad mm6, 9 ;p1(18,11)
  788. pmaddwd mm5, fix_205_n256n039 ;p1(18,17)
  789. paddd mm4, [esi+8*25] ;p1(18,18)
  790. paddd mm4, mm5 ;p1(18,19)
  791. movq mm5, mm4 ;p1(18,20)
  792. paddd mm4, [esi + 8*21] ;p1(18,21)
  793. paddd mm4, const_0x2xx8 ;p1(18,22)
  794. psrad mm4, 9 ;p1(18,23)
  795. packssdw mm6, mm4 ;p1(18,24)
  796. movq [esi + 8*4], mm6 ;p1(18,25)
  797. movq mm4, [esi + 8*20] ;p1(18,26)
  798. psubd mm4, mm7 ;p1(18,27)
  799. paddd mm4, const_0x2xx8 ;p1(18,28)
  800. movq mm7, mm0 ;p1(19,3)
  801. psrad mm4, 9 ;p1(18,29)
  802. movq mm6, [esi + 8*21] ;p1(18,30)
  803. psubd mm6, mm5 ;p1(18,31)
  804. punpcklwd mm7, mm0 ;p1(19,4)
  805. paddd mm6, const_0x2xx8 ;p1(18,32)
  806. psrad mm6, 9 ;p1(18,33)
  807. pmaddwd mm7, fix_029_n089n196 ;p1(19,5)
  808. packssdw mm4, mm6 ;p1(18,34)
  809. movq [esi + 8*10], mm4 ;p1(18,35)
  810. movq mm6, mm3 ;p1(19,0)
  811. punpcklwd mm6, mm1 ;p1(19,1)
  812. movq mm5, mm0 ;p1(19,15)
  813. pmaddwd mm6, fix_n196_n089 ;p1(19,2)
  814. punpckhwd mm5, mm0 ;p1(19,16)
  815. paddd mm6, [esi+8*24] ;p1(19,6)
  816. movq mm4, mm3 ;p1(19,12)
  817. paddd mm6, mm7 ;p1(19,7)
  818. punpckhwd mm4, mm1 ;p1(19,13)
  819. movq mm7, mm6 ;p1(19,8)
  820. pmaddwd mm4, fix_n196_n089   ;p1(19,14)
  821. paddd mm6, [esi + 8*22] ;p1(19,9)
  822. pmaddwd mm5, fix_029_n089n196 ;p1(19,17)
  823. paddd mm6, const_0x2xx8 ;p1(19,10)
  824. psrad mm6, 9 ;p1(19,11)
  825. paddd mm4, [esi+8*25] ;p1(19,18)
  826. paddd mm4, mm5 ;p1(19,19)
  827. movq mm5, mm4 ;p1(19,20)
  828. paddd mm4, [esi + 8*23] ;p1(19,21)
  829. paddd mm4, const_0x2xx8 ;p1(19,22)
  830. psrad mm4, 9 ;p1(19,23)
  831. packssdw mm6, mm4 ;p1(19,24)
  832. movq [esi + 8*6], mm6 ;p1(19,25)
  833. movq mm4, [esi + 8*22] ;p1(19,26)
  834. psubd mm4, mm7 ;p1(19,27)
  835. movq mm6, [esi + 8*23] ;p1(19,30)
  836. paddd mm4, const_0x2xx8 ;p1(19,28)
  837. psubd mm6, mm5 ;p1(19,31)
  838. psrad mm4, 9 ;p1(19,29)
  839. paddd mm6, const_0x2xx8 ;p1(19,32)
  840. psrad mm6, 9 ;p1(19,33)
  841. packssdw mm4, mm6 ;p1(19,34)
  842. movq [esi + 8*8], mm4 ;p1(19,35)
  843. //************************************************************  
  844. mov esi, eax
  845. mov edi, outptr
  846. //transpose 4 rows of wsptr
  847. movq mm0, [esi+8*0] ;tran(0)
  848. movq mm1, mm0 ;tran(1)
  849. movq mm2, [esi+8*2] ;tran(2)
  850. punpcklwd mm0, mm2 ;tran(3)
  851. movq mm3, [esi+8*4] ;tran(5)
  852. punpckhwd mm1, mm2 ;tran(4)
  853. movq mm5, [esi+8*6] ;tran(7)
  854. movq mm4, mm3 ;tran(6)
  855. movq mm6, mm0 ;tran(10)
  856. punpcklwd mm3, mm5 ;tran(8)
  857. movq mm7, mm1 ;tran(11)
  858. punpckldq mm0, mm3 ;tran(12)
  859. punpckhwd mm4, mm5 ;tran(9)
  860. movq [esi+8*0], mm0 ;tran(16)
  861. punpckhdq mm6, mm3 ;tran(13)
  862. movq mm0, [esi+8*1] ;tran(20)
  863. punpckldq mm1, mm4 ;tran(14)
  864. movq [esi+8*2], mm6 ;tran(17)
  865. punpckhdq mm7, mm4 ;tran(15)
  866. movq [esi+8*4], mm1 ;tran(18)
  867. movq mm1, mm0 ;tran(21)
  868. movq mm3, [esi+8*5] ;tran(25)
  869. movq mm2, [esi+8*3] ;tran(22)
  870. movq mm4, mm3 ;tran(26)
  871. punpcklwd mm0, mm2 ;tran(23)
  872. movq [esi+8*6], mm7 ;tran(19)
  873. punpckhwd mm1, mm2 ;tran(24)
  874. movq mm5, [esi+8*7] ;tran(27)
  875. punpcklwd mm3, mm5 ;tran(28)
  876. movq mm6, mm0 ;tran(30)
  877. movq mm7, mm1 ;tran(31)
  878. punpckhdq mm6, mm3 ;tran(33)
  879. punpckhwd mm4, mm5 ;tran(29)
  880. movq mm2, mm6 ;p2(1,0)
  881. punpckhdq mm7, mm4 ;tran(35)
  882. movq mm5, [esi + 8*2] ;p2(1,2)
  883. paddw mm2, mm7 ;p2(1,1)
  884. paddw mm5, [esi + 8*6] ;p2(1,3)
  885. punpckldq mm0, mm3 ;tran(32)
  886. paddw mm2, mm5 ;p2(1,4)
  887. punpckldq mm1, mm4 ;tran(34)
  888. movq mm5, [esi + 8*2] ;p2(3,0)
  889. pmulhw mm2, fix_117_117 ;p2(1,5)
  890. movq mm4, mm7 ;p2(2,0)
  891. pmulhw mm4, fixn089n196p029 ;p2(2,1)
  892. movq mm3, mm6 ;p2(6,0)
  893. pmulhw mm3, fix_n256n039p205 ;p2(6,1)
  894. pmulhw mm5, fix_n089 ;p2(3,1)
  895. movq [eax + 8*24], mm2 ;p2(1,6)
  896. movq mm2, [esi + 8*6] ;p2(4,0)
  897. pmulhw mm2, fix_n196 ;p2(4,1)
  898. paddw mm4, [eax + 8*24] ;p2(5,0)
  899. paddw mm3, [eax + 8*24] ;p2(9,0)
  900. paddw mm5, mm2 ;p2(5,1)
  901. movq mm2, [esi + 8*2] ;p2(7,0)
  902. paddw mm5, mm4 ;p2(5,2)
  903. pmulhw mm2, fix_n039 ;p2(7,1)
  904. movq [esi + 8*1], mm5 ;p2(5,3)
  905. movq mm4, [esi + 8*6] ;p2(8,0)
  906. movq mm5, mm6 ;p2(10,0)
  907. pmulhw mm4, fix_n256 ;p2(8,1)
  908. pmulhw mm5, fix_n039 ;p2(10,1)
  909. pmulhw mm6, fix_n256 ;p2(15,0)
  910. paddw mm2, mm4 ;p2(9,1)
  911. movq mm4, mm7 ;p2(11,0)
  912. pmulhw mm4, fix_n089 ;p2(11,1)
  913. paddw mm2, mm3 ;p2(9,2)
  914. movq [esi + 8*3], mm2 ;p2(9,3)
  915. movq mm3, [esi + 8*2] ;p2(13,0)
  916. pmulhw mm7, fix_n196 ;p2(16,0)
  917. pmulhw mm3, fix_n089n039p150 ;p2(13,1)
  918. paddw mm5, mm4 ;p2(12,0)
  919. paddw mm5, [eax + 8*24] ;p2(14,0)
  920. movq mm2, [esi + 8*6] ;p2(18,0)
  921. pmulhw mm2, fix_n196p307n256 ;p2(18,1)
  922. paddw mm5, mm3 ;p2(14,1)
  923. movq [esi + 8*5], mm5 ;p2(14,2)
  924. paddw mm6, mm7 ;p2(17,0)
  925. paddw mm6, [eax + 8*24] ;p2(19,0)
  926. movq mm3, mm1 ;p2(21,0)
  927. movq mm4, [esi + 8*4] ;p2(20,0)
  928. paddw mm6, mm2 ;p2(19,1)
  929. movq [esi + 8*7], mm6 ;p2(19,2)
  930. movq mm5, mm4 ;p2(20,1)
  931. movq mm7, [esi + 8*0] ;p2(26,0)
  932. pmulhw mm4, fix_054p076 ;p2(20,2)
  933. psubw mm7, mm0 ;p2(27,0)
  934. pmulhw mm3, fix_054 ;p2(21,1)
  935. movq mm2, mm0 ;p2(26,1)
  936. pmulhw mm5, fix_054 ;p2(23,0)
  937. psraw mm7, 3 ;p2(27,1)
  938. paddw mm2, [esi + 8*0] ;p2(26,2)
  939. movq mm6, mm7 ;p2(28,0)
  940. pmulhw mm1, fix_054n184 ;p2(24,0)
  941. psraw mm2, 3 ;p2(26,3)
  942. paddw mm4, mm3 ;p2(22,0)
  943. paddw mm5, mm1 ;p2(25,0)
  944. psubw mm6, mm5 ;p2(29,0)
  945. movq mm3, mm2 ;p2(30,0)
  946. paddw mm2, mm4 ;p2(30,1)
  947. paddw mm7, mm5 ;p2(28,1)
  948. movq mm1, mm2 ;p2(32,0)
  949. psubw mm3, mm4 ;p2(31,0)
  950. paddw mm1, [esi + 8*5] ;p2(32,1)
  951. movq mm0, mm7 ;p2(33,0)
  952. psubw mm2, [esi + 8*5] ;p2(32,2)
  953. movq mm4, mm6 ;p2(34,0)
  954. paddw mm1, const_0x0808 ;p2(32,3)
  955. paddw mm2, const_0x0808 ;p2(32,4)
  956. psraw mm1, 4 ;p2(32,5)
  957. psraw mm2, 4 ;p2(32,6)
  958. paddw mm7, [esi + 8*7] ;p2(33,1)
  959. packuswb mm1, mm2 ;p2(32,7)
  960. psubw mm0, [esi + 8*7] ;p2(33,2)
  961. paddw mm7, const_0x0808 ;p2(33,3)
  962. paddw mm0, const_0x0808 ;p2(33,4)
  963. psraw mm7, 4 ;p2(33,5)
  964. psraw mm0, 4 ;p2(33,6)
  965. paddw mm4, [esi + 8*3] ;p2(34,1)
  966. packuswb mm7, mm0 ;p2(33,7)
  967. psubw mm6, [esi + 8*3] ;p2(34,2)
  968. paddw mm4, const_0x0808 ;p2(34,3)
  969. movq mm5, mm3 ;p2(35,0)
  970. paddw mm6, const_0x0808 ;p2(34,4)
  971. psraw mm4, 4 ;p2(34,5)
  972. psraw mm6, 4 ;p2(34,6)
  973. paddw mm3, [esi + 8*1] ;p2(35,1)
  974. packuswb mm4, mm6 ;p2(34,7)
  975. psubw mm5, [esi + 8*1] ;p2(35,2)
  976. movq mm0, mm1 ;p2(36,0)
  977. paddw mm3, const_0x0808 ;p2(35,3)
  978. paddw mm5, const_0x0808 ;p2(35,4)
  979. punpcklbw mm0, mm7 ;p2(36,1)
  980. psraw mm3, 4 ;p2(35,5)
  981. movq mm2, mm4 ;p2(37,0)
  982. psraw mm5, 4 ;p2(35,6)
  983. movq mm6, mm0 ;p2(38,0)
  984. packuswb mm3, mm5 ;p2(35,7)
  985. lea ebx, [edi] ;p2(42,0)
  986. punpckhbw mm7, mm1 ;p2(36,2)
  987. lea ecx, [edi+8] ;p2(42,1)
  988. punpcklbw mm2, mm3 ;p2(37,1)
  989. lea edx, [edi+16] ;p2(42,2)
  990. punpckhbw mm3, mm4 ;p2(37,2)
  991. punpcklwd mm0, mm2 ;p2(38,1)
  992. movq mm5, mm3 ;p2(39,0)
  993. punpckhwd mm6, mm2 ;p2(38,2)
  994. movq mm1, mm0 ;p2(40,0)
  995. punpcklwd mm3, mm7 ;p2(39,1)
  996. punpckldq mm0, mm3 ;p2(40,1)
  997. punpckhdq mm1, mm3 ;p2(40,2)
  998. movq [ebx], mm0 ;p2(43,0)
  999. punpckhwd mm5, mm7 ;p2(39,2)
  1000. movq [ecx], mm1 ;p2(43,1)
  1001. movq mm4, mm6 ;p2(41,0)
  1002. lea ebx, [edi+24] ;p2(43,3)
  1003. punpckldq mm4, mm5 ;p2(41,1)
  1004. punpckhdq mm6, mm5 ;p2(41,2)
  1005. movq [edx], mm4 ;p2(43,2)
  1006. movq [ebx], mm6 ;p2(43,5)
  1007. //************************************************************
  1008. // Process next 4 rows
  1009. add esi, 64
  1010. add edi, 32
  1011. //transpose next 4 rows of wsptr
  1012. movq mm0, [esi+8*0] ;tran(0)
  1013. movq mm1, mm0 ;tran(1)
  1014. movq mm2, [esi+8*2] ;tran(2)
  1015. punpcklwd mm0, mm2 ;tran(3)
  1016. movq mm3, [esi+8*4] ;tran(5)
  1017. punpckhwd mm1, mm2 ;tran(4)
  1018. movq mm5, [esi+8*6] ;tran(7)
  1019. movq mm4, mm3 ;tran(6)
  1020. movq mm6, mm0 ;tran(10)
  1021. punpcklwd mm3, mm5 ;tran(8)
  1022. movq mm7, mm1 ;tran(11)
  1023. punpckldq mm0, mm3 ;tran(12)
  1024. punpckhwd mm4, mm5 ;tran(9)
  1025. movq [esi+8*0], mm0 ;tran(16)
  1026. punpckhdq mm6, mm3 ;tran(13)
  1027. movq mm0, [esi+8*1] ;tran(20)
  1028. punpckldq mm1, mm4 ;tran(14)
  1029. movq [esi+8*2], mm6 ;tran(17)
  1030. punpckhdq mm7, mm4 ;tran(15)
  1031. movq [esi+8*4], mm1 ;tran(18)
  1032. movq mm1, mm0 ;tran(21)
  1033. movq mm3, [esi+8*5] ;tran(25)
  1034. movq mm2, [esi+8*3] ;tran(22)
  1035. movq mm4, mm3 ;tran(26)
  1036. punpcklwd mm0, mm2 ;tran(23)
  1037. movq [esi+8*6], mm7 ;tran(19)
  1038. punpckhwd mm1, mm2 ;tran(24)
  1039. movq mm5, [esi+8*7] ;tran(27)
  1040. punpcklwd mm3, mm5 ;tran(28)
  1041. movq mm6, mm0 ;tran(30)
  1042. movq mm7, mm1 ;tran(31)
  1043. punpckhdq mm6, mm3 ;tran(33)
  1044. punpckhwd mm4, mm5 ;tran(29)
  1045. movq mm2, mm6 ;p2(1,0)
  1046. punpckhdq mm7, mm4 ;tran(35)
  1047. movq mm5, [esi + 8*2] ;p2(1,2)
  1048. paddw mm2, mm7 ;p2(1,1)
  1049. paddw mm5, [esi + 8*6] ;p2(1,3)
  1050. punpckldq mm0, mm3 ;tran(32)
  1051. paddw mm2, mm5 ;p2(1,4)
  1052. punpckldq mm1, mm4 ;tran(34)
  1053. movq mm5, [esi + 8*2] ;p2(3,0)
  1054. pmulhw mm2, fix_117_117 ;p2(1,5)
  1055. movq mm4, mm7 ;p2(2,0)
  1056. pmulhw mm4, fixn089n196p029 ;p2(2,1)
  1057. movq mm3, mm6 ;p2(6,0)
  1058. pmulhw mm3, fix_n256n039p205 ;p2(6,1)
  1059. pmulhw mm5, fix_n089 ;p2(3,1)
  1060. movq [eax + 8*24], mm2 ;p2(1,6)
  1061. movq mm2, [esi + 8*6] ;p2(4,0)
  1062. pmulhw mm2, fix_n196 ;p2(4,1)
  1063. paddw mm4, [eax + 8*24] ;p2(5,0)
  1064. paddw mm3, [eax + 8*24] ;p2(9,0)
  1065. paddw mm5, mm2 ;p2(5,1)
  1066. movq mm2, [esi + 8*2] ;p2(7,0)
  1067. paddw mm5, mm4 ;p2(5,2)
  1068. pmulhw mm2, fix_n039 ;p2(7,1)
  1069. movq [esi + 8*1], mm5 ;p2(5,3)
  1070. movq mm4, [esi + 8*6] ;p2(8,0)
  1071. movq mm5, mm6 ;p2(10,0)
  1072. pmulhw mm4, fix_n256 ;p2(8,1)
  1073. pmulhw mm5, fix_n039 ;p2(10,1)
  1074. pmulhw mm6, fix_n256 ;p2(15,0)
  1075. paddw mm2, mm4 ;p2(9,1)
  1076. movq mm4, mm7 ;p2(11,0)
  1077. pmulhw mm4, fix_n089 ;p2(11,1)
  1078. paddw mm2, mm3 ;p2(9,2)
  1079. movq [esi + 8*3], mm2 ;p2(9,3)
  1080. movq mm3, [esi + 8*2] ;p2(13,0)
  1081. pmulhw mm7, fix_n196 ;p2(16,0)
  1082. pmulhw mm3, fix_n089n039p150 ;p2(13,1)
  1083. paddw mm5, mm4 ;p2(12,0)
  1084. paddw mm5, [eax + 8*24] ;p2(14,0)
  1085. movq mm2, [esi + 8*6] ;p2(18,0)
  1086. pmulhw mm2, fix_n196p307n256 ;p2(18,1)
  1087. paddw mm5, mm3 ;p2(14,1)
  1088. movq [esi + 8*5], mm5 ;p2(14,2)
  1089. paddw mm6, mm7 ;p2(17,0)
  1090. paddw mm6, [eax + 8*24] ;p2(19,0)
  1091. movq mm3, mm1 ;p2(21,0)
  1092. movq mm4, [esi + 8*4] ;p2(20,0)
  1093. paddw mm6, mm2 ;p2(19,1)
  1094. movq [esi + 8*7], mm6 ;p2(19,2)
  1095. movq mm5, mm4 ;p2(20,1)
  1096. movq mm7, [esi + 8*0] ;p2(26,0)
  1097. pmulhw mm4, fix_054p076 ;p2(20,2)
  1098. psubw mm7, mm0 ;p2(27,0)
  1099. pmulhw mm3, fix_054 ;p2(21,1)
  1100. movq mm2, mm0 ;p2(26,1)
  1101. pmulhw mm5, fix_054 ;p2(23,0)
  1102. psraw mm7, 3 ;p2(27,1)
  1103. paddw mm2, [esi + 8*0] ;p2(26,2)
  1104. movq mm6, mm7 ;p2(28,0)
  1105. pmulhw mm1, fix_054n184 ;p2(24,0)
  1106. psraw mm2, 3 ;p2(26,3)
  1107. paddw mm4, mm3 ;p2(22,0)
  1108. paddw mm5, mm1 ;p2(25,0)
  1109. psubw mm6, mm5 ;p2(29,0)
  1110. movq mm3, mm2 ;p2(30,0)
  1111. paddw mm2, mm4 ;p2(30,1)
  1112. paddw mm7, mm5 ;p2(28,1)
  1113. movq mm1, mm2 ;p2(32,0)
  1114. psubw mm3, mm4 ;p2(31,0)
  1115. paddw mm1, [esi + 8*5] ;p2(32,1)
  1116. movq mm0, mm7 ;p2(33,0)
  1117. psubw mm2, [esi + 8*5] ;p2(32,2)
  1118. movq mm4, mm6 ;p2(34,0)
  1119. paddw mm1, const_0x0808 ;p2(32,3)
  1120. paddw mm2, const_0x0808 ;p2(32,4)
  1121. psraw mm1, 4 ;p2(32,5)
  1122. psraw mm2, 4 ;p2(32,6)
  1123. paddw mm7, [esi + 8*7] ;p2(33,1)
  1124. packuswb mm1, mm2 ;p2(32,7)
  1125. psubw mm0, [esi + 8*7] ;p2(33,2)
  1126. paddw mm7, const_0x0808 ;p2(33,3)
  1127. paddw mm0, const_0x0808 ;p2(33,4)
  1128. psraw mm7, 4 ;p2(33,5)
  1129. psraw mm0, 4 ;p2(33,6)
  1130. paddw mm4, [esi + 8*3] ;p2(34,1)
  1131. packuswb mm7, mm0 ;p2(33,7)
  1132. psubw mm6, [esi + 8*3] ;p2(34,2)
  1133. paddw mm4, const_0x0808 ;p2(34,3)
  1134. movq mm5, mm3 ;p2(35,0)
  1135. paddw mm6, const_0x0808 ;p2(34,4)
  1136. psraw mm4, 4 ;p2(34,5)
  1137. psraw mm6, 4 ;p2(34,6)
  1138. paddw mm3, [esi + 8*1] ;p2(35,1)
  1139. packuswb mm4, mm6 ;p2(34,7)
  1140. psubw mm5, [esi + 8*1] ;p2(35,2)
  1141. movq mm0, mm1 ;p2(36,0)
  1142. paddw mm3, const_0x0808 ;p2(35,3)
  1143. paddw mm5, const_0x0808 ;p2(35,4)
  1144. punpcklbw mm0, mm7 ;p2(36,1)
  1145. psraw mm3, 4 ;p2(35,5)
  1146. movq mm2, mm4 ;p2(37,0)
  1147. psraw mm5, 4 ;p2(35,6)
  1148. movq mm6, mm0 ;p2(38,0)
  1149. packuswb mm3, mm5 ;p2(35,7)
  1150. lea ebx, [edi] ;p2(42,0)
  1151. punpckhbw mm7, mm1 ;p2(36,2)
  1152. lea ecx, [edi+8] ;p2(42,1)
  1153. punpcklbw mm2, mm3 ;p2(37,1)
  1154. lea edx, [edi+16] ;p2(42,2)
  1155. punpckhbw mm3, mm4 ;p2(37,2)
  1156. punpcklwd mm0, mm2 ;p2(38,1)
  1157. movq mm5, mm3 ;p2(39,0)
  1158. punpckhwd mm6, mm2 ;p2(38,2)
  1159. movq mm1, mm0 ;p2(40,0)
  1160. punpcklwd mm3, mm7 ;p2(39,1)
  1161. punpckldq mm0, mm3 ;p2(40,1)
  1162. punpckhdq mm1, mm3 ;p2(40,2)
  1163. movq [ebx], mm0 ;p2(43,0)
  1164. punpckhwd mm5, mm7 ;p2(39,2)
  1165. movq [ecx], mm1 ;p2(43,1)
  1166. movq mm4, mm6 ;p2(41,0)
  1167. lea ebx, [edi+24] ;p2(43,3)
  1168. punpckldq mm4, mm5 ;p2(41,1)
  1169. punpckhdq mm6, mm5 ;p2(41,2)
  1170. movq [edx], mm4 ;p2(43,2)
  1171. movq [ebx], mm6 ;p2(43,5)
  1172. emms
  1173. }
  1174. }