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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 游戏 > 模拟服务器 > 查看源码
KJpegIdct.c
资源名称：SwordOnline.rar　[点击查看]
上传用户：dzyhzl
上传日期：2019-04-29
资源大小：56270k
文件大小：24k
源码类别：

模拟服务器

开发平台：

C/C++

KJpegIdct.c：源码内容
  1. //---------------------------------------------------------------------------
  2. // Sword3 Engine (c) 1999-2000 by Kingsoft
  3. //
  4. // File: KJpegIdct.cpp
  5. // Date: 2000.08.08
  6. // Code: Daniel Wang
  7. // Desc: Jpeg 解码 MMX 版本 8x8 AAN IDCT
  8. // From: Cloud Wu's JPEG Decoder
  9. //---------------------------------------------------------------------------
  10. #include <windows.h>
  11. #include "KJpegLib.h"
  12. /****************************************************************************
  13.   本函数摘取自 developer.intel.com
  15.   ; esi - input and output data pointer
  16.   ; the input data is tranposed and each 16 bit element in the 8x8 matrix 
  17.   ; is left aligned:
  18.   ; for example in 11...1110000 format
  19.   ; If the iDCT is of I macroblock then 0.5 needs to be added to the 
  20.   ; DC Component 
  21.   ; (element[0][0] of the matrix)
  22. ****************************************************************************/
  23. static __int64 x0005000200010001 = 0x0005000200010001;
  24. static __int64 x0040000000000000 = 0x40000000000000;
  25. static __int64 x5a825a825a825a82 = 0x5a825a825a825a82; // 23170
  26. static __int64 x539f539f539f539f = 0x539f539f539f539f; // 21407
  27. static __int64 x4546454645464546 = 0x4546454645464546; // 17734
  28. static __int64 x61f861f861f861f8 = 0x61f861f861f861f8; // 25080
  29. //---------------------------------------------------------------------------
  30. static __int64 scratch1 = 0;
  31. static __int64 scratch3 = 0;
  32. static __int64 scratch5 = 0;
  33. static __int64 scratch7 = 0;
  34. //---------------------------------------------------------------------------
  35. static WORD preSC[64] = {
  36. 16384, 22725, 21407, 19266,  16384, 12873, 8867,  4520,
  37. 22725, 31521, 29692, 26722,  22725, 17855, 12299, 6270,
  38. 21407, 29692, 27969, 25172,  21407, 16819, 11585, 5906,
  39. 19266, 26722, 25172, 22654,  19266, 15137, 10426, 5315,
  40. 16384, 22725, 21407, 19266,  16384, 12873, 8867,  4520,
  41. 12873, 17855, 16819, 15137,  25746, 20228, 13933, 7103,
  42. 17734, 24598, 23170, 20853,  17734, 13933, 9597,  4892,
  43. 18081, 25080, 23624, 21261,  18081, 14206, 9785,  4988
  44. };
  45. //---------------------------------------------------------------------------
  46. // 函数: jpeg_idct_mmx
  47. // 功能: MMX 版本 8x8 AAN IDCT
  48. // 参数: buf 解码后的缓存
  49. // 返回: void
  50. //---------------------------------------------------------------------------
  51. void jpeg_idct_mmx(short* buf)
  52. {
  53. __asm
  54. {
  55. lea ecx, [preSC];
  56. mov esi, buf ; source
  58. ;slot
  60. ; column 0: even part
  61. ; use V4, V12, V0, V8 to produce V22..V25
  62. movq mm0, dword ptr [ecx+8*12] ; maybe the first mul can be done together
  63. ; with the dequantization in iHuff module ?
  64. ;slot
  66. pmulhw mm0, dword ptr [esi+8*12] ; V12
  67. ;slot
  69. movq mm1, dword ptr [ecx+8*4]
  70. ;slot
  72. pmulhw mm1, dword ptr [esi+8*4] ; V4
  73. ;slot
  75. movq mm3, dword ptr [ecx+8*0]
  76. psraw mm0, 1 ; t64=t66
  78. pmulhw mm3, dword ptr [esi+8*0] ; V0
  79. ;slot
  81. movq mm5, dword ptr [ecx+8*8] ; duplicate V4
  82. movq mm2, mm1                           ; added 11/1/96
  84. pmulhw mm5, dword ptr [esi+8*8] ; V8
  85. psubsw mm1, mm0 ; V16
  88. pmulhw mm1, dword ptr x5a825a825a825a82 ; 23170 ->V18
  89. paddsw mm2, mm0                          ; V17
  91. movq mm0, mm2 ; duplicate V17
  92. psraw mm2, 1 ; t75=t82
  94. psraw mm0, 2   ; t72
  95. movq mm4, mm3 ; duplicate V0
  97. paddsw mm3, mm5 ; V19
  98. psubsw mm4, mm5 ; V20 ;mm5 free
  100. ;moved from the block below
  101. movq mm7, dword ptr [ecx+8*10]
  102. psraw mm3, 1 ; t74=t81
  104. movq mm6, mm3 ; duplicate t74=t81
  105. psraw mm4, 2 ; t77=t79
  107. psubsw mm1, mm0 ; V21 ; mm0 free
  108. paddsw mm3, mm2 ; V22
  110. movq mm5, mm1 ; duplicate V21
  111. paddsw mm1, mm4 ; V23
  113. movq dword ptr [esi+8*4], mm3 ; V22
  114. psubsw mm4, mm5 ; V24; mm5 free
  116. movq dword ptr [esi+8*12], mm1 ; V23
  117. psubsw mm6, mm2 ; V25; mm2 free
  119. movq dword ptr [esi+8*0], mm4 ; V24
  120. ;slot
  122. ; keep mm6 alive all along the next block
  123. ;movq dword ptr [esi+8*8], mm6 ; V25
  125. ; column 0: odd part
  126. ; use V2, V6, V10, V14 to produce V31, V39, V40, V41
  128. ;moved above
  129. ;movq mm7, dword ptr [ecx+8*10]
  131. pmulhw mm7, dword ptr [esi+8*10] ; V10
  132. ;slot
  134. movq mm0, dword ptr [ecx+8*6]
  135. ;slot
  137. pmulhw mm0, dword ptr [esi+8*6] ; V6
  138. ;slot
  140. movq mm5, dword ptr [ecx+8*2]
  141. movq mm3, mm7 ; duplicate V10
  143. pmulhw mm5, dword ptr [esi+8*2] ; V2
  144. ;slot
  146. movq mm4, dword ptr [ecx+8*14]
  147. psubsw mm7, mm0 ; V26
  149. pmulhw mm4, dword ptr [esi+8*14] ; V14
  150. paddsw mm3, mm0 ; V29 ; free mm0
  152. movq mm1, mm7 ; duplicate V26
  153. psraw mm3, 1 ; t91=t94
  155. pmulhw mm7, dword ptr x539f539f539f539f  ; V33
  156. psraw mm1, 1 ; t96
  158. movq mm0, mm5 ; duplicate V2
  159. psraw mm4, 2 ; t85=t87
  161. paddsw mm5, mm4 ; V27
  162. psubsw mm0, mm4 ; V28 ; free mm4
  164. movq mm2, mm0 ; duplicate V28
  165. psraw mm5, 1 ; t90=t93
  167. pmulhw mm0, dword ptr x4546454645464546  ; V35
  168. psraw mm2, 1 ; t97
  170. movq mm4, mm5 ; duplicate t90=t93
  171. psubsw mm1, mm2 ; V32 ; free mm2
  173. pmulhw mm1, dword ptr x61f861f861f861f8  ; V36
  174. psllw mm7, 1 ; t107
  176. paddsw mm5, mm3 ; V31
  177. psubsw mm4, mm3 ; V30 ; free mm3
  179. pmulhw mm4, dword ptr x5a825a825a825a82  ; V34
  180. nop ;slot
  182. psubsw mm0, mm1 ; V38
  183. psubsw mm1, mm7 ; V37 ; free mm7
  185. psllw mm1, 1 ; t114
  186. ;move from the next block
  187. movq mm3, mm6 ; duplicate V25
  189. ;move from the next block
  190. movq mm7, dword ptr [esi+8*4] ; V22
  191. psllw mm0, 1 ; t110
  193. psubsw mm0, mm5 ; V39 (mm5 still needed for next block)
  194. psllw mm4, 2 ; t112
  196. ;move from the next block
  197. movq mm2, dword ptr [esi+8*12] ; V23
  198. psubsw mm4, mm0 ; V40
  200. paddsw mm1, mm4 ; V41; free mm0
  201. ;move from the next block
  202. psllw mm2, 1 ; t117=t125
  204. ; column 0: output butterfly
  205. ;move above
  206. ;movq mm3, mm6 ; duplicate V25
  207. ;movq mm7, dword ptr [esi+8*4] ; V22
  208. ;movq mm2, dword ptr [esi+8*12]     ; V23
  209. ;psllw mm2, 1 ; t117=t125
  211. psubsw mm6, mm1 ; tm6
  212. paddsw mm3, mm1 ; tm8; free mm1
  214. movq mm1, mm7 ; duplicate V22
  215. paddsw mm7, mm5 ; tm0
  217. movq dword ptr [esi+8*8], mm3      ; tm8; free mm3
  218. psubsw mm1, mm5 ; tm14; free mm5
  220. movq dword ptr [esi+8*6], mm6 ; tm6; free mm6
  221. movq mm3, mm2 ; duplicate t117=t125
  223. movq mm6, dword ptr [esi+8*0] ; V24
  224. paddsw mm2, mm0 ; tm2
  226. movq dword ptr [esi+8*0], mm7       ; tm0; free mm7
  227. psubsw mm3, mm0 ; tm12; free mm0
  229. movq dword ptr [esi+8*14], mm1 ; tm14; free mm1
  230. psllw mm6, 1 ; t119=t123
  232. movq dword ptr [esi+8*2], mm2 ; tm2; free mm2
  233. movq mm0, mm6 ; duplicate t119=t123
  235. movq dword ptr [esi+8*12], mm3       ; tm12; free mm3
  236. paddsw mm6, mm4 ; tm4
  238. ;moved from next block
  239. movq mm1, dword ptr [ecx+8*5]
  240. psubsw mm0, mm4 ; tm10; free mm4
  242. ;moved from next block
  243. pmulhw mm1, dword ptr [esi+8*5] ; V5
  244. ;slot
  246. movq dword ptr [esi+8*4], mm6 ; tm4; free mm6
  247. ;slot
  249. movq dword ptr [esi+8*10], mm0      ; tm10; free mm0
  250. ;slot
  252. ; column 1: even part
  253. ; use V5, V13, V1, V9 to produce V56..V59
  254. ;moved to prev block
  255. ;movq mm1, dword ptr [ecx+8*5]
  256. ;pmulhw mm1, dword ptr [esi+8*5] ; V5
  258. movq mm7, dword ptr [ecx+8*13]
  259. psllw mm1, 1 ; t128=t130
  261. pmulhw mm7, dword ptr [esi+8*13] ; V13
  262. movq mm2, mm1 ; duplicate t128=t130
  264. movq mm3, dword ptr [ecx+8*1]
  265. ;slot
  267. pmulhw mm3, dword ptr [esi+8*1] ; V1
  268. ;slot
  270. movq mm5, dword ptr [ecx+8*9]
  271. psubsw mm1, mm7 ; V50
  273. pmulhw mm5, dword ptr [esi+8*9] ; V9
  274. paddsw mm2, mm7 ; V51
  276. pmulhw mm1, dword ptr x5a825a825a825a82 ; 23170 ->V52
  277. movq mm6, mm2 ; duplicate V51
  279. psraw mm2, 1 ; t138=t144
  280. movq mm4, mm3 ; duplicate V1
  282. psraw mm6, 2   ; t136
  283. paddsw mm3, mm5 ; V53
  285. psubsw mm4, mm5 ; V54 ;mm5 free
  286. movq mm7, mm3 ; duplicate V53
  288. ;moved from next block
  289. movq mm0, dword ptr [ecx+8*11]
  290. psraw mm4, 1 ; t140=t142
  292. psubsw mm1, mm6 ; V55 ; mm6 free
  293. paddsw mm3, mm2 ; V56
  295. movq mm5, mm4 ; duplicate t140=t142
  296. paddsw mm4, mm1 ; V57
  298. movq dword ptr [esi+8*5], mm3 ; V56
  299. psubsw mm5, mm1 ; V58; mm1 free
  301. movq dword ptr [esi+8*13], mm4 ; V57
  302. psubsw mm7, mm2 ; V59; mm2 free
  304. movq dword ptr [esi+8*9], mm5 ; V58
  305. ;slot
  307. ; keep mm7 alive all along the next block
  308. ;movq dword ptr [esi+8*1], mm7 ; V59
  310. ;moved above
  311. ;movq mm0, dword ptr [ecx+8*11]
  313. pmulhw mm0, dword ptr [esi+8*11] ; V11
  314. ;slot
  316. movq mm6, dword ptr [ecx+8*7]
  317. ;slot
  319. pmulhw mm6, dword ptr [esi+8*7] ; V7
  320. ;slot
  322. movq mm4, dword ptr [ecx+8*15]
  323. movq mm3, mm0 ; duplicate V11
  325. pmulhw mm4, dword ptr [esi+8*15] ; V15
  326. ;slot
  328. movq mm5, dword ptr [ecx+8*3]
  329. psllw mm6,1 ; t146=t152
  331. pmulhw mm5, dword ptr [esi+8*3] ; V3
  332. paddsw mm0, mm6 ; V63
  334. ; note that V15 computation has a correction step:
  335. ; this is a 'magic' constant that rebiases the results to be closer to the expected result
  336. ; this magic constant can be refined to reduce the error even more
  337. ; by doing the correction step in a later stage when the number is actually multiplied by 16
  339. paddw mm4, dword ptr x0005000200010001;
  340. psubsw mm3, mm6 ; V60 ; free mm6
  342. psraw mm0, 1 ; t154=t156
  343. movq mm1, mm3 ; duplicate V60
  345. pmulhw mm1, dword ptr x539f539f539f539f  ; V67
  346. movq mm6, mm5 ; duplicate V3
  348. psraw mm4, 2 ; t148=t150
  349. ;slot
  351. paddsw mm5, mm4 ; V61
  352. psubsw mm6, mm4 ; V62 ; free mm4
  354. movq mm4, mm5 ; duplicate V61
  355. psllw mm1, 1 ; t169
  357. paddsw mm5, mm0 ; V65 -> result
  358. psubsw mm4, mm0 ; V64 ; free mm0
  360. pmulhw mm4, dword ptr x5a825a825a825a82  ; V68
  361. psraw mm3, 1 ; t158
  363. psubsw mm3, mm6 ; V66
  364. movq mm2, mm5 ; duplicate V65
  366. pmulhw mm3, dword ptr x61f861f861f861f8  ; V70
  367. psllw mm6, 1 ; t165
  369. pmulhw mm6, dword ptr x4546454645464546  ; V69
  370. psraw mm2, 1 ; t172
  372. ;moved from next block
  373. movq mm0, dword ptr [esi+8*5] ; V56
  374. psllw mm4, 1 ; t174
  376. ;moved from next block
  377. psraw mm0, 1 ; t177=t188
  378. nop ; slot
  380. psubsw mm6, mm3 ; V72
  381. psubsw mm3, mm1 ; V71 ; free mm1
  383. psubsw mm6, mm2 ; V73 ; free mm2
  384. ;moved from next block
  385. psraw mm5, 1 ; t178=t189
  387. psubsw mm4, mm6 ; V74
  388. ;moved from next block
  389. movq mm1, mm0 ; duplicate t177=t188
  391. paddsw mm3, mm4 ; V75
  392. ;moved from next block
  393. paddsw mm0, mm5 ; tm1
  395. ;location 
  396. ;  5 - V56
  397. ; 13 - V57
  398. ;  9 - V58
  399. ;  X - V59, mm7
  400. ;  X - V65, mm5
  401. ;  X - V73, mm6
  402. ;  X - V74, mm4
  403. ;  X - V75, mm3
  404. ; free mm0, mm1 & mm2
  405. ;move above
  406. ;movq mm0, dword ptr [esi+8*5] ; V56
  407. ;psllw mm0, 1 ; t177=t188 ! new !!
  408. ;psllw mm5, 1 ; t178=t189 ! new !!
  409. ;movq mm1, mm0 ; duplicate t177=t188
  410. ;paddsw mm0, mm5 ; tm1
  412. movq mm2, dword ptr [esi+8*13]      ; V57
  413. psubsw mm1, mm5 ; tm15; free mm5
  415. movq dword ptr [esi+8*1], mm0       ; tm1; free mm0
  416. psraw mm7, 1 ; t182=t184 ! new !!
  418. ;save the store as used directly in the transpose
  419. ;movq dword ptr [esi+8*15], mm1 ; tm15; free mm1
  420. movq mm5, mm7                    ; duplicate t182=t184
  421. psubsw mm7, mm3 ; tm7
  423. paddsw mm5, mm3 ; tm9; free mm3
  424. ;slot
  426. movq mm0, dword ptr [esi+8*9] ; V58
  427. movq mm3, mm2 ; duplicate V57
  429. movq dword ptr [esi+8*7], mm7 ; tm7; free mm7
  430. psubsw mm3, mm6 ; tm13
  432. paddsw mm2, mm6 ; tm3 ; free mm6
  433. ; moved up from the transpose 
  434. movq mm7, mm3                  
  436. ; moved up from the transpose 
  437. punpcklwd mm3, mm1
  438. movq mm6, mm0 ; duplicate V58
  440. movq dword ptr [esi+8*3], mm2 ; tm3; free mm2
  441. paddsw mm0, mm4 ; tm5
  443. psubsw mm6, mm4 ; tm11; free mm4
  444. ; moved up from the transpose 
  445. punpckhwd mm7, mm1
  447. movq dword ptr [esi+8*5], mm0 ; tm5; free mm0
  448. ; moved up from the transpose 
  449. movq mm2, mm5 
  451. ; transpose - M4 part
  452. ;  ---------       ---------
  453. ; | M1 | M2 |     | M1'| M3'|
  454. ;  ---------  -->  ---------
  455. ; | M3 | M4 |     | M2'| M4'|
  456. ;  ---------       ---------
  457. ; Two alternatives: use full dword approach so the following code can be 
  458. ; scheduled before the transpose is done without stores, or use the faster
  459. ; half dword stores (when possible)
  461. movd dword ptr [esi+8*9+4], mm3  ; MS part of tmt9
  462. punpcklwd mm5, mm6
  464. movd dword ptr [esi+8*13+4], mm7 ; MS part of tmt13 
  465. punpckhwd mm2, mm6
  467. movd dword ptr [esi+8*9], mm5 ; LS part of tmt9 
  468. punpckhdq mm5, mm3 ; free mm3
  470. movd dword ptr [esi+8*13], mm2 ; LS part of tmt13
  471. punpckhdq mm2, mm7 ; free mm7
  473. ; moved up from the M3 transpose 
  474. movq mm0, dword ptr [esi+8*8] 
  475. ;slot
  477. ; moved up from the M3 transpose 
  478. movq mm1, dword ptr [esi+8*10] 
  479. ; moved up from the M3 transpose 
  480. movq mm3, mm0 
  482. ; shuffle the rest of the data, and write it with 2 dword writes
  483. movq dword ptr [esi+8*11], mm5  ; tmt11
  484. ; moved up from the M3 transpose 
  485. punpcklwd mm0, mm1
  487. movq dword ptr [esi+8*15], mm2 ; tmt15
  488. ; moved up from the M3 transpose 
  489. punpckhwd mm3, mm1
  491. ; transpose - M3 part
  493. ; moved up to previous code section
  494. ;movq mm0, dword ptr [esi+8*8] 
  495. ;movq mm1, dword ptr [esi+8*10] 
  496. ;movq mm3, mm0 
  497. ;punpcklwd mm0, mm1
  498. ;punpckhwd mm3, mm1
  500. movq mm6, dword ptr [esi+8*12] 
  501. ;slot
  503. movq mm4, dword ptr [esi+8*14] 
  504. movq mm2, mm6 
  506. ; shuffle the data and write out the lower parts of the transposed in 4 dwords
  507. punpcklwd mm6, mm4
  508. movq mm1, mm0
  510. punpckhdq mm1, mm6
  511. movq mm7, mm3
  513. punpckhwd mm2, mm4 ; free mm4
  514. ;slot
  516. punpckldq mm0, mm6 ; free mm6
  517. ;slot
  519. ;moved from next block
  520. movq mm4, dword ptr [esi+8*13] ; tmt13
  521. punpckldq mm3, mm2
  523. punpckhdq mm7, mm2 ; free mm2
  524. ;moved from next block
  525. movq mm5, mm3 ; duplicate tmt5
  527. ; column 1: even part (after transpose)
  529. ;moved above
  530. ;movq mm5, mm3 ; duplicate tmt5
  531. ;movq mm4, dword ptr [esi+8*13] ; tmt13
  533. psubsw mm3, mm4 ; V134
  534. ;slot
  536. pmulhw mm3, dword ptr x5a825a825a825a82 ; 23170 ->V136
  537. ;slot
  539. movq mm6, dword ptr [esi+8*9]   ; tmt9
  540. paddsw mm5, mm4 ; V135 ; mm4 free
  542. movq mm4, mm0 ; duplicate tmt1
  543. paddsw mm0, mm6 ; V137
  545. psubsw mm4, mm6 ; V138 ; mm6 free
  546. psllw mm3, 2 ; t290
  547. psubsw mm3, mm5 ; V139
  548. movq mm6, mm0 ; duplicate V137
  550. paddsw mm0, mm5 ; V140
  551. movq mm2, mm4 ; duplicate V138
  553. paddsw mm2, mm3 ; V141
  554. psubsw mm4, mm3 ; V142 ; mm3 free
  556. movq dword ptr [esi+8*9], mm0  ; V140
  557. psubsw mm6, mm5 ; V143 ; mm5 free
  559. ;moved from next block
  560. movq mm0, dword ptr[esi+8*11] ; tmt11
  561. ;slot
  563. movq dword ptr [esi+8*13], mm2  ; V141
  564. ;moved from next block
  565. movq mm2, mm0 ; duplicate tmt11
  567. ; column 1: odd part (after transpose)
  569. ;moved up to the prev block
  570. ;movq mm0, dword ptr[esi+8*11] ; tmt11
  571. ;movq mm2, mm0 ; duplicate tmt11
  573. movq mm5, dword ptr[esi+8*15] ; tmt15
  574. psubsw mm0, mm7 ; V144
  576. movq mm3, mm0 ; duplicate V144
  577. paddsw mm2, mm7 ; V147 ; free mm7
  579. pmulhw mm0, dword ptr x539f539f539f539f ; 21407-> V151
  580. movq mm7, mm1 ; duplicate tmt3
  582. paddsw mm7, mm5 ; V145
  583. psubsw mm1, mm5 ; V146 ; free mm5
  585. psubsw mm3, mm1 ; V150
  586. movq mm5, mm7 ; duplicate V145
  588. pmulhw mm1, dword ptr x4546454645464546 ; 17734-> V153
  589. psubsw mm5, mm2 ; V148
  591. pmulhw mm3, dword ptr x61f861f861f861f8 ; 25080-> V154
  592. psllw mm0, 2 ; t311
  594. pmulhw mm5, dword ptr x5a825a825a825a82 ; 23170-> V152
  595. paddsw mm7, mm2 ; V149 ; free mm2
  597. psllw mm1, 1 ; t313
  598. nop ; slot
  600. ;without the nop above - freeze here for one clock
  601. ;the nop cleans the mess a little bit
  602. movq mm2, mm3 ; duplicate V154
  603. psubsw mm3, mm0 ; V155 ; free mm0
  605. psubsw mm1, mm2 ; V156 ; free mm2
  606. ;moved from the next block
  607. movq mm2, mm6 ; duplicate V143
  609. ;moved from the next block
  610. movq mm0, dword ptr[esi+8*13] ; V141
  611. psllw mm1, 1 ; t315
  613. psubsw mm1, mm7 ; V157 (keep V149)
  614. psllw mm5, 2 ; t317
  616. psubsw mm5, mm1 ; V158
  617. psllw mm3, 1 ; t319
  619. paddsw mm3, mm5 ; V159
  620. ;slot
  622. ; column 1: output butterfly (after transform)
  623. ;moved to the prev block
  624. ;movq mm2, mm6 ; duplicate V143
  625. ;movq mm0, dword ptr[esi+8*13] ; V141
  627. psubsw mm2, mm3 ; V163
  628. paddsw mm6, mm3 ; V164 ; free mm3
  630. movq mm3, mm4 ; duplicate V142
  631. psubsw mm4, mm5 ; V165 ; free mm5
  633. movq dword ptr scratch7, mm2 ; out7
  634. psraw mm6, 4
  636. psraw mm4, 4
  637. paddsw mm3, mm5 ; V162
  639. movq mm2, dword ptr[esi+8*9] ; V140
  640. movq mm5, mm0 ; duplicate V141
  642. ;in order not to perculate this line up, we read [esi+8*9] very near to this location
  643. // (9)
  644. movq dword ptr [esi+8*9], mm6 ; out9
  645. paddsw mm0, mm1 ; V161
  647. movq dword ptr scratch5, mm3 ; out5
  648. psubsw mm5, mm1 ; V166 ; free mm1
  650. // (11)
  651. movq dword ptr[esi+8*11], mm4 ; out11
  652. psraw mm5, 4
  654. movq dword ptr scratch3, mm0 ; out3
  655. movq mm4, mm2 ; duplicate V140
  657. // (13)
  658. movq dword ptr[esi+8*13], mm5 ; out13
  659. paddsw mm2, mm7 ; V160
  661. ;moved from the next block
  662. movq mm0, dword ptr [esi+8*1] 
  663. psubsw mm4, mm7 ; V167 ; free mm7
  665. ;moved from the next block
  666. movq mm7, dword ptr [esi+8*3] 
  667. psraw mm4, 4
  669. movq dword ptr scratch1, mm2 ; out1
  670. ;moved from the next block
  671. movq mm1, mm0 
  673. // (15)
  674. movq dword ptr[esi+8*15], mm4 ; out15
  675. ;moved from the next block
  676. punpcklwd mm0, mm7
  678. ; transpose - M2 parts
  679. ;moved up to the prev block
  680. ;movq mm0, dword ptr [esi+8*1] 
  681. ;movq mm7, dword ptr [esi+8*3] 
  682. ;movq mm1, mm0 
  683. ;punpcklwd mm0, mm7
  685. movq mm5, dword ptr [esi+8*5] 
  686. punpckhwd mm1, mm7
  688. movq mm4, dword ptr [esi+8*7] 
  689. movq mm3, mm5 
  691. ; shuffle the data and write out the lower parts of the trasposed in 4 dwords
  692. movd dword ptr [esi+8*8], mm0  ; LS part of tmt8
  693. punpcklwd mm5, mm4
  695. movd dword ptr [esi+8*12], mm1  ; LS part of tmt12
  696. punpckhwd mm3, mm4
  698. movd dword ptr [esi+8*8+4], mm5  ; MS part of tmt8
  699. punpckhdq mm0, mm5 ; tmt10
  701. movd dword ptr [esi+8*12+4], mm3  ; MS part of tmt12
  702. punpckhdq mm1, mm3 ; tmt14
  705. ; transpose - M1 parts
  706. movq mm7, dword ptr [esi] 
  707. ;slot
  709. movq mm2, dword ptr [esi+8*2] 
  710. movq mm6, mm7 
  712. movq mm5, dword ptr [esi+8*4] 
  713. punpcklwd mm7, mm2
  715. movq mm4, dword ptr [esi+8*6] 
  716. punpckhwd mm6, mm2 ; free mm2
  718. movq mm3, mm5 
  719. punpcklwd mm5, mm4
  721. punpckhwd mm3, mm4 ; free mm4
  722. movq mm2, mm7
  724. movq mm4, mm6
  725. punpckldq mm7, mm5 ; tmt0
  727. punpckhdq mm2, mm5 ; tmt2 ; free mm5
  728. ;slot
  730. ; shuffle the rest of the data, and write it with 2 dword writes
  731. punpckldq mm6, mm3 ; tmt4
  732. ;move from next block
  733. movq mm5, mm2 ; duplicate tmt2
  735. punpckhdq mm4, mm3 ; tmt6 ; free mm3
  736. ;move from next block
  737. movq mm3, mm0 ; duplicate tmt10
  739. ; column 0: odd part (after transpose)
  740. ;moved up to prev block
  741. ;movq mm3, mm0 ; duplicate tmt10
  742. ;movq mm5, mm2 ; duplicate tmt2
  744. psubsw mm0, mm4 ; V110
  745. paddsw mm3, mm4 ; V113 ; free mm4
  747. movq mm4, mm0 ; duplicate V110
  748. paddsw mm2, mm1 ; V111
  750. pmulhw mm0, dword ptr x539f539f539f539f ; 21407-> V117
  751. psubsw mm5, mm1 ; V112 ; free mm1
  753. psubsw mm4, mm5 ; V116
  754. movq mm1, mm2 ; duplicate V111
  756. pmulhw mm5, dword ptr x4546454645464546 ; 17734-> V119
  757. psubsw mm2, mm3 ; V114
  759. pmulhw mm4, dword ptr x61f861f861f861f8 ; 25080-> V120
  760. paddsw mm1, mm3 ; V115 ; free mm3
  762. pmulhw mm2, dword ptr x5a825a825a825a82 ; 23170-> V118
  763. psllw mm0, 2 ; t266
  765. movq dword ptr[esi+8*0], mm1 ; save V115
  766. psllw mm5, 1 ; t268
  768. psubsw mm5, mm4 ; V122
  769. psubsw mm4, mm0 ; V121 ; free mm0
  771. psllw mm5, 1 ; t270
  772. ;slot
  774. psubsw mm5, mm1 ; V123 ; free mm1
  775. psllw mm2, 2 ; t272
  777. psubsw mm2, mm5 ; V124 (keep V123)
  778. psllw mm4, 1 ; t274
  781. movq dword ptr[esi+8*2], mm5 ; save V123 ; free mm5
  782. paddsw mm4, mm2 ; V125 (keep V124)
  784. ; column 0: even part (after transpose)
  785. movq mm0, dword ptr[esi+8*12] ; tmt12
  786. movq mm3, mm6 ; duplicate tmt4
  788. psubsw mm6, mm0 ; V100
  789. paddsw mm3, mm0 ; V101 ; free mm0
  791. pmulhw mm6, dword ptr x5a825a825a825a82  ; 23170 ->V102
  792. movq mm5, mm7 ; duplicate tmt0
  794. movq mm1, dword ptr[esi+8*8] ; tmt8
  795. ;slot
  797. paddsw mm7, mm1 ; V103
  798. psubsw mm5, mm1 ; V104 ; free mm1
  800. movq mm0, mm7 ; duplicate V103
  801. psllw mm6, 2 ; t245
  803. paddsw mm7, mm3 ; V106
  804. movq mm1, mm5 ; duplicate V104
  806. psubsw mm6, mm3 ; V105
  807. psubsw mm0, mm3 ; V109; free mm3   
  809. paddsw mm5, mm6 ; V107
  810. psubsw mm1, mm6 ; V108 ; free mm6
  812. ; column 0: output butterfly (after transform)
  813. movq mm3, mm1 ; duplicate V108
  814. paddsw mm1, mm2 ; out4
  816. psraw mm1, 4
  817. psubsw mm3, mm2 ; out10 ; free mm2
  819. psraw mm3, 4
  820. movq mm6, mm0 ; duplicate V109
  822. // (4)
  823. movq dword ptr[esi+8*4], mm1 ; out4 ; free mm1
  824. psubsw mm0, mm4 ; out6
  826. // (10)
  827. movq dword ptr[esi+8*10], mm3 ; out10 ; free mm3
  828. psraw mm0, 4
  830. paddsw mm6, mm4 ; out8 ; free mm4
  831. movq mm1, mm7 ; duplicate V106
  833. // (6)
  834. movq dword ptr[esi+8*6], mm0 ; out6 ; free mm0
  835. psraw mm6, 4
  837. movq mm4, dword ptr[esi+8*0]  ; V115
  838. ;slot
  840. // (8)
  841. movq dword ptr[esi+8*8], mm6 ; out8 ; free mm6
  842. movq mm2, mm5 ; duplicate V107
  844. movq mm3, dword ptr[esi+8*2]  ; V123
  845. paddsw mm7, mm4 ; out0
  847. ;moved up from next block
  848. movq mm0, dword ptr scratch3
  849. psraw mm7, 4
  851. ;moved up from next block
  852. movq mm6, dword ptr scratch5
  853. psubsw mm1, mm4 ; out14 ; free mm4
  855. paddsw mm5, mm3 ; out2
  856. psraw mm1, 4
  858. // (0)
  859. // movq dword ptr[esi], mm7 ; out0 ; free mm7
  860. psraw mm5, 4
  862. // (14)
  863. movq dword ptr[esi+8*14], mm1 ; out14 ; free mm1
  864. psubsw mm2, mm3 ; out12 ; free mm3
  866. // (2)
  867. // movq dword ptr[esi+8*2], mm5 ; out2 ; free mm5
  868. psraw mm2, 4
  870. ;moved up to the prev block
  871. movq mm4, dword ptr scratch7
  872. ;moved up to the prev block
  873. psraw mm0, 4
  875. // (12)
  876. //movq dword ptr[esi+8*12], mm2 ; out12 ; free mm2
  878. ;moved up to the prev block
  879. movq mm1, dword ptr scratch1
  880. psraw mm6, 4
  882. packsswb mm7,mm7 
  884. psraw mm1, 4
  886. movd [esi],mm7; //out0
  888. psraw mm4, 4
  890. packsswb mm1,mm1;
  893. // (3)
  894. //movq dword ptr [esi+8*3], mm0 ; out3
  896. // (5)
  897. //movq dword ptr [esi+8*5], mm6 ; out5
  899. // (7)
  900. //movq dword ptr [esi+8*7], mm4 ; out7
  902. // (1)
  903. //movq dword ptr [esi+8*1], mm1 ; out1
  905. // 将 16bit 变成 8bit
  906. // mm1=1,mm4=7,mm6=5,mm0=3,mm2=12,mm5=2,mm7=0
  908. packsswb mm7,[esi+8*4]; //mm7=4
  909. packsswb mm6,mm6;
  910. packsswb mm5,mm5;
  911. punpckhdq mm7,mm6;
  912. movd [esi+4*1],mm1; //out1   free mm1
  913. movd [esi+4*2],mm5; //out2   free mm5
  914. packsswb mm0,mm0;
  915. packsswb mm1,[esi+8*6]; //mm1=6
  916. packsswb mm4,mm4;
  917. movd [esi+4*3],mm0; //out3   free mm0
  918. punpckhdq mm1,mm4;
  919. packsswb mm3,[esi+8*8]; //mm3=8
  920. movq [esi+4*4],mm7; //out4,5 free mm6,mm7
  921. packsswb mm0,[esi+8*9]; //mm0=9
  922. packsswb mm6,[esi+8*10];//mm6=10
  923. movq [esi+4*6],mm1;     //out6,7 free mm1,mm4
  924. packsswb mm1,[esi+8*11];//mm1=11
  925. punpckhdq mm3,mm0;
  926. punpckhdq mm6,mm1;
  927. packsswb mm2,mm2;
  928. packsswb mm4,[esi+8*13];//mm4=13
  929. movq [esi+4*8],mm3;     //out8,9 free mm0,mm3
  930. punpckhdq mm2,mm4;
  931. packsswb mm0,[esi+8*14];//mm3=14
  932. packsswb mm5,[esi+8*15];//mm5=15
  933. movq [esi+4*10],mm6; //out10,11  free mm1,mm6
  934. punpckhdq mm0,mm5;
  935. movq [esi+4*12],mm2; //out12,13  free mm3,mm4
  936. movq [esi+4*14],mm0;
  938. emms;
  939. }
  940. }
  941. //---------------------------------------------------------------------------