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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > 流媒体/Mpeg4/MP4 > 查看源码
text_idct_mmx.c
资源名称：encore50src.zip　[点击查看]
上传用户：enenge
上传日期：2007-01-08
资源大小：96k
文件大小：20k
源码类别：

流媒体/Mpeg4/MP4

开发平台：

Visual C++

text_idct_mmx.c：源码内容
  2. /**************************************************************************
  3.  *                                                                        *
  4.  * This code is developed by John Funnell.  This software is an           *
  5.  * an implementation of a part of one or more MPEG-4 Video tools as       *
  6.  * specified in ISO/IEC 14496-2 standard.  Those intending to use this    *
  7.  * software module in hardware or software products are advised that its  *
  8.  * use may infringe existing patents or copyrights, and any such use      *
  9.  * would be at such party's own risk.  The original developer of this     *
  10.  * software module and his/her company, and subsequent editors and their  *
  11.  * companies (including Project Mayo), will have no liability for use of  *
  12.  * this software or modifications or derivatives thereof.                 *
  13.  *                                                                        *
  14.  * Project Mayo gives users of the Codec a license to this software       *
  15.  * module or modifications thereof for use in hardware or software        *
  16.  * products claiming conformance to the MPEG-4 Video Standard as          *
  17.  * described in the Open DivX license.                                    *
  18.  *                                                                        *
  19.  * The complete Open DivX license can be found at                         *
  20.  * http://www.projectmayo.com/opendivx/license.php .                      *
  21.  *                                                                        *
  22.  **************************************************************************/
  24. /**************************************************************************
  25.  *
  26.  *  text_idct_mmx.c
  27.  *
  28.  *  Copyright (C) 2001  Project Mayo
  29.  *
  30.  *  John Funnell
  31.  *
  32.  *  DivX Advance Research Center <darc@projectmayo.com>
  33.  *
  34.  **************************************************************************/
  36. /* This file contains idct functions with MMX optimization.               */
  37. /* Please see acknowledgement below.                                      */
  39. /*
  40. ;
  41. ; MMX iDCT
  42. ;
  43. ; Originally provided by Intel at AP-922
  44. ; http://developer.intel.com/vtune/cbts/strmsimd/922down.htm
  45. ; (See more app notes at http://developer.intel.com/vtune/cbts/strmsimd/appnotes.htm)
  46. ; but in a limited edition.
  47. ; New macro implements a column part for precise iDCT
  48. ; The routine precision now satisfies IEEE standard 1180-1990. 
  49. ;
  50. ; Copyright (c) 2000-2001 Peter Gubanov <peter@elecard.net.ru>
  51. ; Rounding trick Copyright (c) 2000 Michel Lespinasse <walken@zoy.org>
  52. ;
  53. ; http://www.elecard.com/peter/idct.html
  54. ; http://www.linuxvideo.org/mpeg2dec/
  55. ;
  56. ; Conversion to MS-style inline C format (c) 2001 Project Mayo Inc. John Funnell
  57. ;
  58. ; Version v1.0.0 10 January 2001
  59. ;
  60. */
  63. /*
  64.  if you are not setting the name of you MMX idct function with compile flag MMX_IDCT_NAME,
  65.  then you need to choose a name to insert in here:
  66. */
  67. #define DEFAULT_FUNCTION_NAME       Fast_IDCT
  71. #ifndef MMX_IDCT_NAME
  72. #define MMX_IDCT_NAME DEFAULT_FUNCTION_NAME
  73. #endif
  77. /* some types */
  78. #define int16_t short
  79. #define int32_t int
  80. #define uint64_t unsigned __int64
  82. #define MPTR qword ptr 
  85. /* some conversions from readable integer constants to qwords */
  86. #define SHORT4_TO_QWORD(A, B, C, D)   ( 
  87.     ( ( ((uint64_t)(A)) & ((uint64_t)(0xffff)) )  <<  0 ) | 
  88.     ( ( ((uint64_t)(B)) & ((uint64_t)(0xffff)) )  << 16 ) | 
  89.     ( ( ((uint64_t)(C)) & ((uint64_t)(0xffff)) )  << 32 ) | 
  90.     ( ( ((uint64_t)(D)) & ((uint64_t)(0xffff)) )  << 48 ) ) 
  92. #define INT2_TO_QWORD(A, B)    ( 
  93.     ( ( ((uint64_t)(A)) & ((uint64_t)(0xffffffff)) )  <<  0 ) | 
  94.     ( ( ((uint64_t)(B)) & ((uint64_t)(0xffffffff)) )  << 32 ) )
  99. #define BITS_INV_ACC      5                             /* 4 or 5 for IEEE */
  100. #define SHIFT_INV_ROW     (16 - BITS_INV_ACC)
  101. #define SHIFT_INV_COL     (1 + BITS_INV_ACC)
  102. #define RND_INV_ROW       (1024 * (6 - BITS_INV_ACC))   /* 1 << (SHIFT_INV_ROW-1) */
  103. #define RND_INV_COL       (16 * (BITS_INV_ACC - 3))     /* 1 << (SHIFT_INV_COL-1) */
  104. #define RND_INV_CORR      (RND_INV_COL - 1)             /* correction -1.0 and round */
  105. #define BITS_FRW_ACC      3                             /* 2 or 3 for accuracy */
  106. #define SHIFT_FRW_COL     BITS_FRW_ACC
  107. #define SHIFT_FRW_ROW     (BITS_FRW_ACC + 17)
  108. #define RND_FRW_ROW       (262144 * (BITS_FRW_ACC - 1)) /* 1 << (SHIFT_FRW_ROW-1) */
  115. static const uint64_t   tg_1_16 = SHORT4_TO_QWORD( 13036,  13036,  13036,  13036 );
  116. static const uint64_t   tg_2_16 = SHORT4_TO_QWORD( 27146,  27146,  27146,  27146 );
  117. static const uint64_t   tg_3_16 = SHORT4_TO_QWORD(-21746, -21746, -21746, -21746 );
  118. static const uint64_t ocos_4_16 = SHORT4_TO_QWORD( 23170,  23170,  23170,  23170 );
  123. /* rounding value that is added before the row transform's rounding */
  126. #if   SHIFT_INV_ROW == 12
  127. static const uint64_t rounder[8] = {
  128. INT2_TO_QWORD( 65536, 65536), /* rounder_0 */       
  129. INT2_TO_QWORD(  7195,  7195), /* rounder_1 */ 
  130. INT2_TO_QWORD(  4520,  4520), /* rounder_2 */ 
  131. INT2_TO_QWORD(  2407,  2407), /* rounder_3 */ 
  132. INT2_TO_QWORD(     0,     0), /* rounder_4 */ 
  133. INT2_TO_QWORD(   240,   240), /* rounder_5 */
  134. INT2_TO_QWORD(  1024,  1024), /* rounder_6 */
  135. INT2_TO_QWORD(  1024,  1024)  /* rounder_7 */
  136. };
  137. #elif SHIFT_INV_ROW == 11
  138. static const uint64_t rounder[2*8] = {
  139. INT2_TO_QWORD( 65536, 65536), /* rounder_0 */
  140. INT2_TO_QWORD(  3597,  3597), /* rounder_1 */
  141. INT2_TO_QWORD(  2260,  2260), /* rounder_2 */
  142. INT2_TO_QWORD(  1203,  1203), /* rounder_3 */
  143. INT2_TO_QWORD(   0,   0), /* rounder_4 */
  144. INT2_TO_QWORD(   120, 120), /* rounder_5 */
  145. INT2_TO_QWORD(   512, 512), /* rounder_6 */
  146. INT2_TO_QWORD(   512, 512)  /* rounder_7 */
  147. };
  148. #endif
  152. /*
  153. ;=============================================================================
  154. ;
  155. ; The first stage iDCT 8x8 - inverse DCTs of rows
  156. ;
  157. ;-----------------------------------------------------------------------------
  158. ; The 8-point inverse DCT direct algorithm
  159. ;-----------------------------------------------------------------------------
  160. ;
  161. ; static const short w[32] = {
  162. ; FIX(cos_4_16), FIX(cos_2_16), FIX(cos_4_16), FIX(cos_6_16),
  163. ; FIX(cos_4_16), FIX(cos_6_16), -FIX(cos_4_16), -FIX(cos_2_16),
  164. ; FIX(cos_4_16), -FIX(cos_6_16), -FIX(cos_4_16), FIX(cos_2_16),
  165. ; FIX(cos_4_16), -FIX(cos_2_16), FIX(cos_4_16), -FIX(cos_6_16),
  166. ; FIX(cos_1_16), FIX(cos_3_16), FIX(cos_5_16), FIX(cos_7_16),
  167. ; FIX(cos_3_16), -FIX(cos_7_16), -FIX(cos_1_16), -FIX(cos_5_16),
  168. ; FIX(cos_5_16), -FIX(cos_1_16), FIX(cos_7_16), FIX(cos_3_16),
  169. ; FIX(cos_7_16), -FIX(cos_5_16), FIX(cos_3_16), -FIX(cos_1_16) };
  170. ;
  171. ; #define DCT_8_INV_ROW(x, y)
  172. ;{
  173. ; int a0, a1, a2, a3, b0, b1, b2, b3;
  174. ;
  175. ; a0 =x[0]*w[0]+x[2]*w[1]+x[4]*w[2]+x[6]*w[3];
  176. ; a1 =x[0]*w[4]+x[2]*w[5]+x[4]*w[6]+x[6]*w[7];
  177. ; a2 = x[0] * w[ 8] + x[2] * w[ 9] + x[4] * w[10] + x[6] * w[11];
  178. ; a3 = x[0] * w[12] + x[2] * w[13] + x[4] * w[14] + x[6] * w[15];
  179. ; b0 = x[1] * w[16] + x[3] * w[17] + x[5] * w[18] + x[7] * w[19];
  180. ; b1 = x[1] * w[20] + x[3] * w[21] + x[5] * w[22] + x[7] * w[23];
  181. ; b2 = x[1] * w[24] + x[3] * w[25] + x[5] * w[26] + x[7] * w[27];
  182. ; b3 = x[1] * w[28] + x[3] * w[29] + x[5] * w[30] + x[7] * w[31];
  183. ;
  184. ; y[0] = SHIFT_ROUND ( a0 + b0 );
  185. ; y[1] = SHIFT_ROUND ( a1 + b1 );
  186. ; y[2] = SHIFT_ROUND ( a2 + b2 );
  187. ; y[3] = SHIFT_ROUND ( a3 + b3 );
  188. ; y[4] = SHIFT_ROUND ( a3 - b3 );
  189. ; y[5] = SHIFT_ROUND ( a2 - b2 );
  190. ; y[6] = SHIFT_ROUND ( a1 - b1 );
  191. ; y[7] = SHIFT_ROUND ( a0 - b0 );
  192. ;}
  193. ;
  194. ;-----------------------------------------------------------------------------
  195. ;
  196. ; In this implementation the outputs of the iDCT-1D are multiplied
  197. ; for rows 0,4 - by cos_4_16,
  198. ; for rows 1,7 - by cos_1_16,
  199. ; for rows 2,6 - by cos_2_16,
  200. ; for rows 3,5 - by cos_3_16
  201. ; and are shifted to the left for better accuracy
  202. ;
  203. ; For the constants used,
  204. ; FIX(float_const) = (short) (float_const * (1<<15) + 0.5)
  205. ;
  206. ;=============================================================================
  207. */
  210. /*
  211. ;=============================================================================
  212. ; code for MMX
  213. ;=============================================================================
  214. ; Table for rows 0,4 - constants are multiplied by cos_4_16
  215. */
  216. static const uint64_t tab_i_04[8*8] = {
  217. /* Table for rows 0,4 - constants are multiplied by cos_4_16 */
  218.   SHORT4_TO_QWORD( 16384,  16384,  16384, -16384 ),   /* movq-> w06 w04 w02 w00 */
  219. SHORT4_TO_QWORD( 21407,   8867,   8867, -21407 ),   /*        w07 w05 w03 w01 */
  220. SHORT4_TO_QWORD( 16384, -16384,  16384,  16384 ),   /*        w14 w12 w10 w08 */
  221. SHORT4_TO_QWORD( -8867,  21407, -21407,  -8867 ),   /*        w15 w13 w11 w09 */
  222. SHORT4_TO_QWORD( 22725,  12873,  19266, -22725 ),   /*        w22 w20 w18 w16 */
  223. SHORT4_TO_QWORD( 19266,   4520,  -4520, -12873 ),   /*        w23 w21 w19 w17 */
  224. SHORT4_TO_QWORD( 12873,   4520,   4520,  19266 ),   /*        w30 w28 w26 w24 */
  225. SHORT4_TO_QWORD(-22725,  19266, -12873, -22725 ),   /*        w31 w29 w27 w25 */
  226. /* Table for rows 1,7 - constants are multiplied by cos_1_16 */
  227. SHORT4_TO_QWORD( 22725,  22725,  22725, -22725 ),   /* movq-> w06 w04 w02 w00 */
  228. SHORT4_TO_QWORD( 29692,  12299,  12299, -29692 ),   /*        w07 w05 w03 w01 */
  229. SHORT4_TO_QWORD( 22725, -22725,  22725,  22725 ),   /*        w14 w12 w10 w08 */
  230. SHORT4_TO_QWORD(-12299,  29692, -29692, -12299 ),   /*        w15 w13 w11 w09 */
  231. SHORT4_TO_QWORD( 31521,  17855,  26722, -31521 ),   /*        w22 w20 w18 w16 */
  232. SHORT4_TO_QWORD( 26722,   6270,  -6270, -17855 ),   /*        w23 w21 w19 w17 */
  233. SHORT4_TO_QWORD( 17855,   6270,   6270,  26722 ),   /*        w30 w28 w26 w24 */
  234. SHORT4_TO_QWORD(-31521,  26722, -17855, -31521 ),   /*        w31 w29 w27 w25 */
  235. /* Table for rows 2,6 - constants are multiplied by cos_2_16 */
  236. SHORT4_TO_QWORD( 21407,  21407,  21407, -21407 ),   /* movq-> w06 w04 w02 w00 */
  237. SHORT4_TO_QWORD( 27969,  11585,  11585, -27969 ),   /*        w07 w05 w03 w01 */
  238. SHORT4_TO_QWORD( 21407, -21407,  21407,  21407 ),   /*        w14 w12 w10 w08 */
  239. SHORT4_TO_QWORD(-11585,  27969, -27969, -11585 ),   /*        w15 w13 w11 w09 */
  240. SHORT4_TO_QWORD( 29692,  16819,  25172, -29692 ),   /*        w22 w20 w18 w16 */
  241. SHORT4_TO_QWORD( 25172,   5906,  -5906, -16819 ),   /*        w23 w21 w19 w17 */
  242. SHORT4_TO_QWORD( 16819,   5906,   5906,  25172 ),   /*        w30 w28 w26 w24 */
  243. SHORT4_TO_QWORD(-29692,  25172, -16819, -29692 ),   /*        w31 w29 w27 w25 */
  244. /* Table for rows 3,5 - constants are multiplied by cos_3_16 */
  245. SHORT4_TO_QWORD( 19266,  19266,  19266, -19266 ),   /* movq-> w06 w04 w02 w00 */
  246. SHORT4_TO_QWORD( 25172,  10426,  10426, -25172 ),   /*        w07 w05 w03 w01 */
  247. SHORT4_TO_QWORD( 19266, -19266,  19266,  19266 ),   /*        w14 w12 w10 w08 */
  248. SHORT4_TO_QWORD(-10426,  25172, -25172, -10426 ),   /*        w15 w13 w11 w09 */
  249. SHORT4_TO_QWORD( 26722,  15137,  22654, -26722 ),   /*        w22 w20 w18 w16 */
  250. SHORT4_TO_QWORD( 22654,   5315,  -5315, -15137 ),   /*        w23 w21 w19 w17 */
  251. SHORT4_TO_QWORD( 15137,   5315,   5315,  22654 ),   /*        w30 w28 w26 w24 */
  252. SHORT4_TO_QWORD(-26722,  22654, -15137, -26722 ),   /*        w31 w29 w27 w25 */
  253. /* Table for rows 0,4 - constants are multiplied by cos_4_16 */
  254.   SHORT4_TO_QWORD( 16384,  16384,  16384, -16384 ),   /* movq-> w06 w04 w02 w00 */
  255. SHORT4_TO_QWORD( 21407,   8867,   8867, -21407 ),   /*        w07 w05 w03 w01 */
  256. SHORT4_TO_QWORD( 16384, -16384,  16384,  16384 ),   /*        w14 w12 w10 w08 */
  257. SHORT4_TO_QWORD(-8867,  21407, -21407,  -8867  ),   /*        w15 w13 w11 w09 */
  258. SHORT4_TO_QWORD( 22725,  12873,  19266, -22725 ),   /*        w22 w20 w18 w16 */
  259. SHORT4_TO_QWORD( 19266,   4520,  -4520, -12873 ),   /*        w23 w21 w19 w17 */
  260. SHORT4_TO_QWORD( 12873,   4520,   4520,  19266 ),   /*        w30 w28 w26 w24 */
  261. SHORT4_TO_QWORD(-22725,  19266, -12873, -22725 ),   /*        w31 w29 w27 w25 */
  262. /* Table for rows 3,5 - constants are multiplied by cos_3_16 */
  263. SHORT4_TO_QWORD( 19266,  19266,  19266, -19266 ),   /* movq-> w06 w04 w02 w00 */
  264. SHORT4_TO_QWORD( 25172,  10426,  10426, -25172 ),   /*        w07 w05 w03 w01 */
  265. SHORT4_TO_QWORD( 19266, -19266,  19266,  19266 ),   /*        w14 w12 w10 w08 */
  266. SHORT4_TO_QWORD(-10426,  25172, -25172, -10426 ),   /*        w15 w13 w11 w09 */
  267. SHORT4_TO_QWORD( 26722,  15137,  22654, -26722 ),   /*        w22 w20 w18 w16 */
  268. SHORT4_TO_QWORD( 22654,   5315,  -5315, -15137 ),   /*        w23 w21 w19 w17 */
  269. SHORT4_TO_QWORD( 15137,   5315,   5315,  22654 ),   /*        w30 w28 w26 w24 */
  270. SHORT4_TO_QWORD(-26722,  22654, -15137, -26722 ),   /*        w31 w29 w27 w25 */
  271. /* Table for rows 2,6 - constants are multiplied by cos_2_16 */
  272. SHORT4_TO_QWORD( 21407,  21407,  21407, -21407 ),   /* movq-> w06 w04 w02 w00 */
  273. SHORT4_TO_QWORD( 27969,  11585,  11585, -27969 ),   /*        w07 w05 w03 w01 */
  274. SHORT4_TO_QWORD( 21407, -21407,  21407,  21407 ),   /*        w14 w12 w10 w08 */
  275. SHORT4_TO_QWORD(-11585,  27969, -27969, -11585 ),   /*        w15 w13 w11 w09 */
  276. SHORT4_TO_QWORD( 29692,  16819,  25172, -29692 ),   /*        w22 w20 w18 w16 */
  277. SHORT4_TO_QWORD( 25172,   5906,  -5906, -16819 ),   /*        w23 w21 w19 w17 */
  278. SHORT4_TO_QWORD( 16819,   5906,   5906,  25172 ),   /*        w30 w28 w26 w24 */
  279. SHORT4_TO_QWORD(-29692,  25172, -16819, -29692 ),   /*        w31 w29 w27 w25 */
  280. /* Table for rows 1,7 - constants are multiplied by cos_1_16 */
  281. SHORT4_TO_QWORD( 22725,  22725,  22725, -22725 ),   /* movq-> w06 w04 w02 w00 */
  282. SHORT4_TO_QWORD( 29692,  12299,  12299, -29692 ),   /*        w07 w05 w03 w01 */
  283. SHORT4_TO_QWORD( 22725, -22725,  22725,  22725 ),   /*        w14 w12 w10 w08 */
  284. SHORT4_TO_QWORD(-12299,  29692, -29692, -12299 ),   /*        w15 w13 w11 w09 */
  285. SHORT4_TO_QWORD( 31521,  17855,  26722, -31521 ),   /*        w22 w20 w18 w16 */
  286. SHORT4_TO_QWORD( 26722,   6270,  -6270, -17855 ),   /*        w23 w21 w19 w17 */
  287. SHORT4_TO_QWORD( 17855,   6270,   6270,  26722 ),   /*        w30 w28 w26 w24 */
  288. SHORT4_TO_QWORD(-31521,  26722, -17855, -31521 ),   /*        w31 w29 w27 w25 */
  289. };
  293. uint64_t tempMatrix[64/4];
  296. #define INP      eax
  297. #define ROUNDER  ebx
  298. #define TABLE    ecx
  299. #define OUTP     edx
  305. /* the function itself */
  306. void MMX_IDCT_NAME
  307.               (int16_t *block) {
  308. __asm {
  310. push eax
  311. push ebx 
  312. push ecx
  313. push edx 
  314. push edi
  316. mov INP, block
  317. lea OUTP,  [tempMatrix]
  318. lea TABLE, [tab_i_04]
  319. lea ROUNDER, [rounder]
  321. mov edi, -8;                 ; i=-8
  323. align 16
  325. rowloop:
  327. movq mm0, MPTR [INP]  ; 0 ; x3 x2 x1 x0
  329. movq mm1, MPTR [INP+8] ; 1 ; x7 x6 x5 x4
  330. movq mm2, mm0  ; 2 ; x3 x2 x1 x0
  332. movq mm3, MPTR [TABLE] ; 3 ; w06 w04 w02 w00
  333. punpcklwd mm0, mm1  ; x5 x1 x4 x0
  335. movq mm5, mm0  ; 5 ; x5 x1 x4 x0
  336. punpckldq mm0, mm0  ; x4 x0 x4 x0
  338. movq mm4, MPTR [TABLE+8]  ; 4 ; w07 w05 w03 w01
  339. punpckhwd mm2, mm1 ; 1 ; x7 x3 x6 x2
  341. pmaddwd mm3, mm0  ; x4*w06+x0*w04 x4*w02+x0*w00
  342. movq mm6, mm2  ; 6  ; x7 x3 x6 x2
  344. movq mm1, MPTR [TABLE+32]  ; 1  ; w22 w20 w18 w16
  345. punpckldq mm2, mm2  ; x6 x2 x6 x2
  347. pmaddwd mm4, mm2  ; x6*w07+x2*w05 x6*w03+x2*w01
  348. punpckhdq mm5, mm5  ; x5 x1 x5 x1
  350. pmaddwd mm0, MPTR [TABLE+16]  ; x4*w14+x0*w12 x4*w10+x0*w08
  351. punpckhdq mm6, mm6  ; x7 x3 x7 x3
  353. movq mm7, MPTR [TABLE+40]  ; 7  ; w23 w21 w19 w17
  354. pmaddwd mm1, mm5  ; x5*w22+x1*w20 x5*w18+x1*w16
  356. paddd mm3, MPTR [ROUNDER]  ; +rounder
  357. pmaddwd mm7, mm6  ; x7*w23+x3*w21 x7*w19+x3*w17
  359. pmaddwd mm2, MPTR [TABLE+24]  ; x6*w15+x2*w13 x6*w11+x2*w09
  360. paddd mm3, mm4  ; 4  ; a1=sum(even1) a0=sum(even0)
  362. pmaddwd mm5, MPTR [TABLE+48]  ; x5*w30+x1*w28 x5*w26+x1*w24
  363. movq mm4, mm3  ; 4  ; a1 a0
  365. pmaddwd mm6, MPTR [TABLE+56]  ; x7*w31+x3*w29 x7*w27+x3*w25
  366. paddd mm1, mm7  ; 7  ; b1=sum(odd1) b0=sum(odd0)
  368. paddd mm0, MPTR [ROUNDER] ; +rounder
  369. psubd mm3, mm1  ; a1-b1 a0-b0
  371. psrad mm3, SHIFT_INV_ROW  ; y6=a1-b1 y7=a0-b0
  372. paddd mm1, mm4  ; 4  ; a1+b1 a0+b0
  374. paddd mm0, mm2  ; 2  ; a3=sum(even3) a2=sum(even2)
  375. psrad mm1, SHIFT_INV_ROW  ; y1=a1+b1 y0=a0+b0
  377. paddd mm5, mm6  ; 6  ; b3=sum(odd3) b2=sum(odd2)
  378. movq mm4, mm0  ; 4  ; a3 a2
  380. paddd mm0, mm5  ; a3+b3 a2+b2
  381. psubd mm4, mm5  ; 5  ; a3-b3 a2-b2
  383. psrad mm0, SHIFT_INV_ROW  ; y3=a3+b3 y2=a2+b2
  384. psrad mm4, SHIFT_INV_ROW  ; y4=a3-b3 y5=a2-b2
  386. packssdw mm1, mm0  ; 0  ; y3 y2 y1 y0
  387. packssdw mm4, mm3  ; 3  ; y6 y7 y4 y5
  389. movq mm7, mm4  ; 7  ; y6 y7 y4 y5
  390. psrld mm4, 16  ; 0 y6 0 y4
  392. pslld mm7, 16  ; y7 0 y5 0
  393. movq MPTR [OUTP], mm1  ; 1  ; save y3 y2 y1 y0
  394.                              
  395. por mm7, mm4  ; 4  ; y7 y6 y5 y4
  396. movq MPTR [OUTP+8], mm7  ; 7  ; save y7 y6 y5 y4
  399. add INP,  16                    ; add 1 row to input pointer
  400. add ROUNDER, 8                  ; go to next rounding values
  401. add OUTP, 16                    ; add 1 row to output pointer
  402. add TABLE,64                    ; move to next section of table
  404. add edi, 1
  405. jne rowloop;
  408. /* column code starts here... */
  411. lea INP,   [tempMatrix] 
  412. mov OUTP,  block
  413. mov edi, -2;                 ; i=-2
  415. align 16
  416. colloop:
  418. // movq mm0, MPTR [tg_3_16]
  419. movq mm0, tg_3_16
  421. movq mm3, MPTR [INP+16*3]
  422. movq mm1, mm0 ; tg_3_16
  424. movq mm5, MPTR [INP+16*5]
  425. pmulhw mm0, mm3 ; x3*(tg_3_16-1)
  427. // movq mm4, MPTR [tg_1_16]
  428. movq mm4, tg_1_16
  429. pmulhw mm1, mm5 ; x5*(tg_3_16-1)
  431. movq mm7, MPTR [INP+16*7]
  432. movq mm2, mm4 ; tg_1_16
  434. movq mm6, MPTR [INP+16*1]
  435. pmulhw mm4, mm7 ; x7*tg_1_16
  437. paddsw mm0, mm3 ; x3*tg_3_16
  438. pmulhw mm2, mm6 ; x1*tg_1_16
  440. paddsw mm1, mm3 ; x3+x5*(tg_3_16-1)
  441. psubsw mm0, mm5 ; x3*tg_3_16-x5 = tm35
  443. // movq mm3, MPTR [ocos_4_16]
  444. movq mm3, ocos_4_16
  445. paddsw mm1, mm5 ; x3+x5*tg_3_16 = tp35
  447. paddsw mm4, mm6 ; x1+tg_1_16*x7 = tp17
  448. psubsw mm2, mm7 ; x1*tg_1_16-x7 = tm17
  450. movq mm5, mm4 ; tp17
  451. movq mm6, mm2 ; tm17
  453. paddsw mm5, mm1 ; tp17+tp35 = b0
  454. psubsw mm6, mm0 ; tm17-tm35 = b3
  456. psubsw mm4, mm1 ; tp17-tp35 = t1
  457. paddsw mm2, mm0 ; tm17+tm35 = t2
  459. // movq mm7, MPTR [tg_2_16]
  460. movq mm7, tg_2_16
  461. movq mm1, mm4 ; t1
  463. ; movq MPTR [SCRATCH+0], mm5 ; save b0
  464. movq MPTR [OUTP+3*16], mm5 ; save b0
  465. paddsw mm1, mm2 ; t1+t2
  467. ; movq MPTR [SCRATCH+8], mm6 ; save b3
  468. movq MPTR [OUTP+5*16], mm6 ; save b3
  469. psubsw mm4, mm2 ; t1-t2
  471. movq mm5, MPTR [INP+2*16]
  472. movq mm0, mm7 ; tg_2_16
  474. movq mm6, MPTR [INP+6*16]
  475. pmulhw mm0, mm5 ; x2*tg_2_16
  477. pmulhw mm7, mm6 ; x6*tg_2_16
  478. ; slot
  479. pmulhw mm1, mm3 ; ocos_4_16*(t1+t2) = b1/2
  480. ; slot
  481. movq mm2, MPTR [INP+0*16]
  482. pmulhw mm4, mm3 ; ocos_4_16*(t1-t2) = b2/2
  484. psubsw mm0, mm6 ; t2*tg_2_16-x6 = tm26
  485. movq mm3, mm2 ; x0
  487. movq mm6, MPTR [INP+4*16]
  488. paddsw mm7, mm5 ; x2+x6*tg_2_16 = tp26
  490. paddsw mm2, mm6 ; x0+x4 = tp04
  491. psubsw mm3, mm6 ; x0-x4 = tm04
  493. movq mm5, mm2 ; tp04
  494. movq mm6, mm3 ; tm04
  496. psubsw mm2, mm7 ; tp04-tp26 = a3
  497. paddsw mm3, mm0 ; tm04+tm26 = a1
  499. paddsw mm1, mm1 ; b1
  500. paddsw mm4, mm4 ; b2
  502. paddsw mm5, mm7 ; tp04+tp26 = a0
  503. psubsw mm6, mm0 ; tm04-tm26 = a2
  505. movq mm7, mm3 ; a1
  506. movq mm0, mm6 ; a2
  508. paddsw mm3, mm1 ; a1+b1
  509. paddsw mm6, mm4 ; a2+b2
  511. psraw mm3, SHIFT_INV_COL ; dst1
  512. psubsw mm7, mm1 ; a1-b1
  514. psraw mm6, SHIFT_INV_COL ; dst2
  515. psubsw mm0, mm4 ; a2-b2
  517. ; movq mm1, MPTR [SCRATCH+0] ; load b0
  518. movq mm1, MPTR [OUTP+3*16] ; load b0
  519. psraw mm7, SHIFT_INV_COL ; dst6
  521. movq mm4, mm5 ; a0
  522. psraw mm0, SHIFT_INV_COL ; dst5
  524. movq MPTR [OUTP+1*16], mm3
  525. paddsw mm5, mm1 ; a0+b0
  527. movq MPTR [OUTP+2*16], mm6
  528. psubsw mm4, mm1 ; a0-b0
  530. ; movq mm3, MPTR [SCRATCH+8] ; load b3
  531. movq mm3, MPTR [OUTP+5*16] ; load b3
  532. psraw mm5, SHIFT_INV_COL ; dst0
  534. movq mm6, mm2 ; a3
  535. psraw mm4, SHIFT_INV_COL ; dst7
  537. movq MPTR [OUTP+5*16], mm0
  538. paddsw mm2, mm3 ; a3+b3
  540. movq MPTR [OUTP+6*16], mm7
  541. psubsw mm6, mm3 ; a3-b3
  543. movq MPTR [OUTP+0*16], mm5
  544. psraw mm2, SHIFT_INV_COL ; dst3
  546. movq MPTR [OUTP+7*16], mm4
  547. psraw mm6, SHIFT_INV_COL ; dst4
  549. movq MPTR [OUTP+3*16], mm2
  551. movq MPTR [OUTP+4*16], mm6
  553. lea INP,   [tempMatrix] 
  554. mov OUTP,  block
  555. add INP,  8                    ; add 4 cols to input pointer
  556. add OUTP, 8                    ; add 4 cols to output pointer
  558. add edi, 1
  559. jne colloop
  561. pop edi
  562. pop edx 
  563. pop ecx
  564. pop ebx
  565. pop eax
  567. emms 
  569. }; /* __asm */
  571. };