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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 嵌入式/单片机编程 > Windows CE > 查看源码
resample_mmx.c
资源名称：tcpmp.rar　[点击查看]
上传用户：wstnjxml
上传日期：2014-04-03
资源大小：7248k
文件大小：15k
源码类别：

Windows CE

开发平台：

C/C++

resample_mmx.c：源码内容
  2. // MMX optimizations from Michael Niedermayer (michaelni@gmx.at) (under GPL)
  4. /* optimization TODO / NOTES 
  5.     movntq is slightly faster (0.5% with the current test.c benchmark) 
  6. (but thats just test.c so that needs to be testd in reallity)
  7. and it would mean (C / MMX2 / MMX / 3DNOW) versions 
  8. */
  10. static uint64_t __attribute__((aligned(8))) magicF2W= 0x43c0000043c00000LL;
  11. static uint64_t __attribute__((aligned(8))) wm1010= 0xFFFF0000FFFF0000LL;
  12. static uint64_t __attribute__((aligned(8))) wm0101= 0x0000FFFF0000FFFFLL;
  13. static uint64_t __attribute__((aligned(8))) wm1100= 0xFFFFFFFF00000000LL;
  15. static int a52_resample_MONO_to_5_MMX(float * _f, int16_t * s16){
  16.     int32_t * f = (int32_t *) _f;
  17. asm volatile(
  18. "movl $-512, %%esi nt"
  19. "movq "MANGLE(magicF2W)", %%mm7 nt"
  20. "movq "MANGLE(wm1100)", %%mm3 nt"
  21. "movq "MANGLE(wm0101)", %%mm4 nt"
  22. "movq "MANGLE(wm1010)", %%mm5 nt"
  23. "pxor %%mm6, %%mm6 nt"
  24. "1: nt"
  25. "movq (%1, %%esi, 2), %%mm0 nt"
  26. "movq 8(%1, %%esi, 2), %%mm1 nt"
  27. "leal (%%esi, %%esi, 4), %%edi nt"
  28. "psubd %%mm7, %%mm0 nt"
  29. "psubd %%mm7, %%mm1 nt"
  30. "packssdw %%mm1, %%mm0 nt"
  31. "movq %%mm0, %%mm1 nt"
  32. "pand %%mm4, %%mm0 nt"
  33. "pand %%mm5, %%mm1 nt"
  34. "movq %%mm6, (%0, %%edi) nt" // 0 0 0 0
  35. "movd %%mm0, 8(%0, %%edi) nt" // A 0
  36. "pand %%mm3, %%mm0 nt"
  37. "movd %%mm6, 12(%0, %%edi) nt" // 0 0
  38. "movd %%mm1, 16(%0, %%edi) nt" // 0 B
  39. "pand %%mm3, %%mm1 nt"
  40. "movd %%mm6, 20(%0, %%edi) nt" // 0 0
  41. "movq %%mm0, 24(%0, %%edi) nt" // 0 0 C 0
  42. "movq %%mm1, 32(%0, %%edi) nt" // 0 0 0 B
  43. "addl $8, %%esi nt"
  44. " jnz 1b nt"
  45. "emms nt"
  46. :: "r" (s16+1280), "r" (f+256)
  47. :"%esi", "%edi", "memory"
  48. );
  49.     return 5*256;
  50. }
  52. static int a52_resample_STEREO_to_2_MMX(float * _f, int16_t * s16){
  53.     int32_t * f = (int32_t *) _f;
  54. /* benchmark scores are 0.3% better with SSE but we would need to set bias=0 and premultiply it
  55. #ifdef HAVE_SSE
  56. asm volatile(
  57. "movl $-1024, %%esi nt"
  58. "1: nt"
  59. "cvtps2pi (%1, %%esi), %%mm0 nt"
  60. "cvtps2pi 1024(%1, %%esi), %%mm2nt"
  61. "movq %%mm0, %%mm1 nt"
  62. "punpcklwd %%mm2, %%mm0 nt"
  63. "punpckhwd %%mm2, %%mm1 nt"
  64. "movq %%mm0, (%0, %%esi) nt"
  65. "movq %%mm1, 8(%0, %%esi) nt"
  66. "addl $16, %%esi nt"
  67. " jnz 1b nt"
  68. "emms nt"
  69. :: "r" (s16+512), "r" (f+256)
  70. :"%esi", "memory"
  71. );*/
  72. asm volatile(
  73. "movl $-1024, %%esi nt"
  74. "movq "MANGLE(magicF2W)", %%mm7 nt"
  75. "1: nt"
  76. "movq (%1, %%esi), %%mm0 nt"
  77. "movq 8(%1, %%esi), %%mm1 nt"
  78. "movq 1024(%1, %%esi), %%mm2 nt"
  79. "movq 1032(%1, %%esi), %%mm3 nt"
  80. "psubd %%mm7, %%mm0 nt"
  81. "psubd %%mm7, %%mm1 nt"
  82. "psubd %%mm7, %%mm2 nt"
  83. "psubd %%mm7, %%mm3 nt"
  84. "packssdw %%mm1, %%mm0 nt"
  85. "packssdw %%mm3, %%mm2 nt"
  86. "movq %%mm0, %%mm1 nt"
  87. "punpcklwd %%mm2, %%mm0 nt"
  88. "punpckhwd %%mm2, %%mm1 nt"
  89. "movq %%mm0, (%0, %%esi) nt"
  90. "movq %%mm1, 8(%0, %%esi) nt"
  91. "addl $16, %%esi nt"
  92. " jnz 1b nt"
  93. "emms nt"
  94. :: "r" (s16+512), "r" (f+256)
  95. :"%esi", "memory"
  96. );
  97.     return 2*256;
  98. }
  100. static int a52_resample_3F_to_5_MMX(float * _f, int16_t * s16){
  101.     int32_t * f = (int32_t *) _f;
  102. asm volatile(
  103. "movl $-1024, %%esi nt"
  104. "movq "MANGLE(magicF2W)", %%mm7 nt"
  105. "pxor %%mm6, %%mm6 nt"
  106. "movq %%mm7, %%mm5 nt"
  107. "punpckldq %%mm6, %%mm5 nt"
  108. "1: nt"
  109. "movd (%1, %%esi), %%mm0 nt"
  110. "punpckldq 2048(%1, %%esi), %%mm0nt"
  111. "movd 1024(%1, %%esi), %%mm1 nt"
  112. "punpckldq 4(%1, %%esi), %%mm1 nt"
  113. "movd 2052(%1, %%esi), %%mm2 nt"
  114. "movq %%mm7, %%mm3 nt"
  115. "punpckldq 1028(%1, %%esi), %%mm3nt"
  116. "movd 8(%1, %%esi), %%mm4 nt"
  117. "punpckldq 2056(%1, %%esi), %%mm4nt"
  118. "leal (%%esi, %%esi, 4), %%edi nt"
  119. "sarl $1, %%edi nt"
  120. "psubd %%mm7, %%mm0 nt"
  121. "psubd %%mm7, %%mm1 nt"
  122. "psubd %%mm5, %%mm2 nt"
  123. "psubd %%mm7, %%mm3 nt"
  124. "psubd %%mm7, %%mm4 nt"
  125. "packssdw %%mm6, %%mm0 nt"
  126. "packssdw %%mm2, %%mm1 nt"
  127. "packssdw %%mm4, %%mm3 nt"
  128. "movq %%mm0, (%0, %%edi) nt"
  129. "movq %%mm1, 8(%0, %%edi) nt"
  130. "movq %%mm3, 16(%0, %%edi) nt"
  132. "movd 1032(%1, %%esi), %%mm1 nt"
  133. "punpckldq 12(%1, %%esi), %%mm1nt"
  134. "movd 2060(%1, %%esi), %%mm2 nt"
  135. "movq %%mm7, %%mm3 nt"
  136. "punpckldq 1036(%1, %%esi), %%mm3nt"
  137. "pxor %%mm0, %%mm0 nt"
  138. "psubd %%mm7, %%mm1 nt"
  139. "psubd %%mm5, %%mm2 nt"
  140. "psubd %%mm7, %%mm3 nt"
  141. "packssdw %%mm1, %%mm0 nt"
  142. "packssdw %%mm3, %%mm2 nt"
  143. "movq %%mm0, 24(%0, %%edi) nt"
  144. "movq %%mm2, 32(%0, %%edi) nt"
  146. "addl $16, %%esi nt"
  147. " jnz 1b nt"
  148. "emms nt"
  149. :: "r" (s16+1280), "r" (f+256)
  150. :"%esi", "%edi", "memory"
  151. );
  152.     return 5*256;
  153. }
  155. static int a52_resample_2F_2R_to_4_MMX(float * _f, int16_t * s16){
  156.     int32_t * f = (int32_t *) _f;
  157. asm volatile(
  158. "movl $-1024, %%esi nt"
  159. "movq "MANGLE(magicF2W)", %%mm7 nt"
  160. "1: nt"
  161. "movq (%1, %%esi), %%mm0 nt"
  162. "movq 8(%1, %%esi), %%mm1 nt"
  163. "movq 1024(%1, %%esi), %%mm2 nt"
  164. "movq 1032(%1, %%esi), %%mm3 nt"
  165. "psubd %%mm7, %%mm0 nt"
  166. "psubd %%mm7, %%mm1 nt"
  167. "psubd %%mm7, %%mm2 nt"
  168. "psubd %%mm7, %%mm3 nt"
  169. "packssdw %%mm1, %%mm0 nt"
  170. "packssdw %%mm3, %%mm2 nt"
  171. "movq 2048(%1, %%esi), %%mm3 nt"
  172. "movq 2056(%1, %%esi), %%mm4 nt"
  173. "movq 3072(%1, %%esi), %%mm5 nt"
  174. "movq 3080(%1, %%esi), %%mm6 nt"
  175. "psubd %%mm7, %%mm3 nt"
  176. "psubd %%mm7, %%mm4 nt"
  177. "psubd %%mm7, %%mm5 nt"
  178. "psubd %%mm7, %%mm6 nt"
  179. "packssdw %%mm4, %%mm3 nt"
  180. "packssdw %%mm6, %%mm5 nt"
  181. "movq %%mm0, %%mm1 nt"
  182. "movq %%mm3, %%mm4 nt"
  183. "punpcklwd %%mm2, %%mm0 nt"
  184. "punpckhwd %%mm2, %%mm1 nt"
  185. "punpcklwd %%mm5, %%mm3 nt"
  186. "punpckhwd %%mm5, %%mm4 nt"
  187. "movq %%mm0, %%mm2 nt"
  188. "movq %%mm1, %%mm5 nt"
  189. "punpckldq %%mm3, %%mm0 nt"
  190. "punpckhdq %%mm3, %%mm2 nt"
  191. "punpckldq %%mm4, %%mm1 nt"
  192. "punpckhdq %%mm4, %%mm5 nt"
  193. "movq %%mm0, (%0, %%esi,2) nt"
  194. "movq %%mm2, 8(%0, %%esi,2) nt"
  195. "movq %%mm1, 16(%0, %%esi,2) nt"
  196. "movq %%mm5, 24(%0, %%esi,2) nt"
  197. "addl $16, %%esi nt"
  198. " jnz 1b nt"
  199. "emms nt"
  200. :: "r" (s16+1024), "r" (f+256)
  201. :"%esi", "memory"
  202. );
  203.     return 4*256;
  204. }
  206. static int a52_resample_3F_2R_to_5_MMX(float * _f, int16_t * s16){
  207.     int32_t * f = (int32_t *) _f;
  208. asm volatile(
  209. "movl $-1024, %%esi nt"
  210. "movq "MANGLE(magicF2W)", %%mm7 nt"
  211. "1: nt"
  212. "movd (%1, %%esi), %%mm0 nt"
  213. "punpckldq 2048(%1, %%esi), %%mm0nt"
  214. "movd 3072(%1, %%esi), %%mm1 nt"
  215. "punpckldq 4096(%1, %%esi), %%mm1nt"
  216. "movd 1024(%1, %%esi), %%mm2 nt"
  217. "punpckldq 4(%1, %%esi), %%mm2 nt"
  218. "movd 2052(%1, %%esi), %%mm3 nt"
  219. "punpckldq 3076(%1, %%esi), %%mm3nt"
  220. "movd 4100(%1, %%esi), %%mm4 nt"
  221. "punpckldq 1028(%1, %%esi), %%mm4nt"
  222. "movd 8(%1, %%esi), %%mm5 nt"
  223. "punpckldq 2056(%1, %%esi), %%mm5nt"
  224. "leal (%%esi, %%esi, 4), %%edi nt"
  225. "sarl $1, %%edi nt"
  226. "psubd %%mm7, %%mm0 nt"
  227. "psubd %%mm7, %%mm1 nt"
  228. "psubd %%mm7, %%mm2 nt"
  229. "psubd %%mm7, %%mm3 nt"
  230. "psubd %%mm7, %%mm4 nt"
  231. "psubd %%mm7, %%mm5 nt"
  232. "packssdw %%mm1, %%mm0 nt"
  233. "packssdw %%mm3, %%mm2 nt"
  234. "packssdw %%mm5, %%mm4 nt"
  235. "movq %%mm0, (%0, %%edi) nt"
  236. "movq %%mm2, 8(%0, %%edi) nt"
  237. "movq %%mm4, 16(%0, %%edi) nt"
  239. "movd 3080(%1, %%esi), %%mm0 nt"
  240. "punpckldq 4104(%1, %%esi), %%mm0nt"
  241. "movd 1032(%1, %%esi), %%mm1 nt"
  242. "punpckldq 12(%1, %%esi), %%mm1nt"
  243. "movd 2060(%1, %%esi), %%mm2 nt"
  244. "punpckldq 3084(%1, %%esi), %%mm2nt"
  245. "movd 4108(%1, %%esi), %%mm3 nt"
  246. "punpckldq 1036(%1, %%esi), %%mm3nt"
  247. "psubd %%mm7, %%mm0 nt"
  248. "psubd %%mm7, %%mm1 nt"
  249. "psubd %%mm7, %%mm2 nt"
  250. "psubd %%mm7, %%mm3 nt"
  251. "packssdw %%mm1, %%mm0 nt"
  252. "packssdw %%mm3, %%mm2 nt"
  253. "movq %%mm0, 24(%0, %%edi) nt"
  254. "movq %%mm2, 32(%0, %%edi) nt"
  256. "addl $16, %%esi nt"
  257. " jnz 1b nt"
  258. "emms nt"
  259. :: "r" (s16+1280), "r" (f+256)
  260. :"%esi", "%edi", "memory"
  261. );
  262.     return 5*256;
  263. }
  265. static int a52_resample_MONO_LFE_to_6_MMX(float * _f, int16_t * s16){
  266.     int32_t * f = (int32_t *) _f;
  267. asm volatile(
  268. "movl $-1024, %%esi nt"
  269. "movq "MANGLE(magicF2W)", %%mm7 nt"
  270. "pxor %%mm6, %%mm6 nt"
  271. "1: nt"
  272. "movq 1024(%1, %%esi), %%mm0 nt"
  273. "movq 1032(%1, %%esi), %%mm1 nt"
  274. "movq (%1, %%esi), %%mm2 nt"
  275. "movq 8(%1, %%esi), %%mm3 nt"
  276. "psubd %%mm7, %%mm0 nt"
  277. "psubd %%mm7, %%mm1 nt"
  278. "psubd %%mm7, %%mm2 nt"
  279. "psubd %%mm7, %%mm3 nt"
  280. "packssdw %%mm1, %%mm0 nt"
  281. "packssdw %%mm3, %%mm2 nt"
  282. "movq %%mm0, %%mm1 nt"
  283. "punpcklwd %%mm2, %%mm0 nt"
  284. "punpckhwd %%mm2, %%mm1 nt"
  285. "leal (%%esi, %%esi, 2), %%edi nt"
  286. "movq %%mm6, (%0, %%edi) nt"
  287. "movd %%mm0, 8(%0, %%edi) nt"
  288. "punpckhdq %%mm0, %%mm0 nt"
  289. "movq %%mm6, 12(%0, %%edi) nt"
  290. "movd %%mm0, 20(%0, %%edi) nt"
  291. "movq %%mm6, 24(%0, %%edi) nt"
  292. "movd %%mm1, 32(%0, %%edi) nt"
  293. "punpckhdq %%mm1, %%mm1 nt"
  294. "movq %%mm6, 36(%0, %%edi) nt"
  295. "movd %%mm1, 44(%0, %%edi) nt"
  296. "addl $16, %%esi nt"
  297. " jnz 1b nt"
  298. "emms nt"
  299. :: "r" (s16+1536), "r" (f+256)
  300. :"%esi", "%edi", "memory"
  301. );
  302.     return 6*256;
  303. }
  305. static int a52_resample_STEREO_LFE_to_6_MMX(float * _f, int16_t * s16){
  306.     int32_t * f = (int32_t *) _f;
  307. asm volatile(
  308. "movl $-1024, %%esi nt"
  309. "movq "MANGLE(magicF2W)", %%mm7 nt"
  310. "pxor %%mm6, %%mm6 nt"
  311. "1: nt"
  312. "movq 1024(%1, %%esi), %%mm0 nt"
  313. "movq 2048(%1, %%esi), %%mm1 nt"
  314. "movq (%1, %%esi), %%mm5 nt" 
  315. "psubd %%mm7, %%mm0 nt"
  316. "psubd %%mm7, %%mm1 nt"
  317. "psubd %%mm7, %%mm5 nt"
  318. "leal (%%esi, %%esi, 2), %%edi nt"
  320. "pxor %%mm4, %%mm4 nt"
  321. "packssdw %%mm5, %%mm0 nt" // FfAa
  322. "packssdw %%mm4, %%mm1 nt" // 00Bb
  323. "punpckhwd %%mm0, %%mm4 nt" // F0f0
  324. "punpcklwd %%mm1, %%mm0 nt" // BAba
  325. "movq %%mm0, %%mm1 nt" // BAba
  326. "punpckldq %%mm4, %%mm3 nt" // f0XX
  327. "punpckldq %%mm6, %%mm0 nt" // 00ba
  328. "punpckhdq %%mm1, %%mm3 nt" // BAf0
  330. "movq %%mm0, (%0, %%edi) nt" // 00ba
  331. "punpckhdq %%mm4, %%mm0 nt" // F000
  332. "movq %%mm3, 8(%0, %%edi) nt" // BAf0
  333. "movq %%mm0, 16(%0, %%edi) nt" // F000
  334. "addl $8, %%esi nt"
  335. " jnz 1b nt"
  336. "emms nt"
  337. :: "r" (s16+1536), "r" (f+256)
  338. :"%esi", "%edi", "memory"
  339. );
  340.     return 6*256;
  341. }
  343. static int a52_resample_3F_LFE_to_6_MMX(float * _f, int16_t * s16){
  344.     int32_t * f = (int32_t *) _f;
  345. asm volatile(
  346. "movl $-1024, %%esi nt"
  347. "movq "MANGLE(magicF2W)", %%mm7 nt"
  348. "pxor %%mm6, %%mm6 nt"
  349. "1: nt"
  350. "movq 1024(%1, %%esi), %%mm0 nt"
  351. "movq 3072(%1, %%esi), %%mm1 nt"
  352. "movq 2048(%1, %%esi), %%mm4 nt"
  353. "movq (%1, %%esi), %%mm5 nt" 
  354. "psubd %%mm7, %%mm0 nt"
  355. "psubd %%mm7, %%mm1 nt"
  356. "psubd %%mm7, %%mm4 nt"
  357. "psubd %%mm7, %%mm5 nt"
  358. "leal (%%esi, %%esi, 2), %%edi nt"
  360. "packssdw %%mm4, %%mm0 nt" // EeAa
  361. "packssdw %%mm5, %%mm1 nt" // FfBb
  362. "movq %%mm0, %%mm2 nt" // EeAa
  363. "punpcklwd %%mm1, %%mm0 nt" // BAba
  364. "punpckhwd %%mm1, %%mm2 nt" // FEfe
  365. "movq %%mm0, %%mm1 nt" // BAba
  366. "punpckldq %%mm6, %%mm0 nt" // 00ba
  367. "punpckhdq %%mm1, %%mm1 nt" // BABA
  369. "movq %%mm0, (%0, %%edi) nt"
  370. "punpckhdq %%mm2, %%mm0 nt" // FE00
  371. "punpckldq %%mm1, %%mm2 nt" // BAfe
  372. "movq %%mm2, 8(%0, %%edi) nt"
  373. "movq %%mm0, 16(%0, %%edi) nt"
  374. "addl $8, %%esi nt"
  375. " jnz 1b nt"
  376. "emms nt"
  377. :: "r" (s16+1536), "r" (f+256)
  378. :"%esi", "%edi", "memory"
  379. );
  380.     return 6*256;
  381. }
  383. static int a52_resample_2F_2R_LFE_to_6_MMX(float * _f, int16_t * s16){
  384.     int32_t * f = (int32_t *) _f;
  385. asm volatile(
  386. "movl $-1024, %%esi nt"
  387. "movq "MANGLE(magicF2W)", %%mm7 nt"
  388. // "pxor %%mm6, %%mm6 nt"
  389. "1: nt"
  390. "movq 1024(%1, %%esi), %%mm0 nt"
  391. "movq 2048(%1, %%esi), %%mm1 nt"
  392. "movq 3072(%1, %%esi), %%mm2 nt"
  393. "movq 4096(%1, %%esi), %%mm3 nt"
  394. "movq (%1, %%esi), %%mm5 nt" 
  395. "psubd %%mm7, %%mm0 nt"
  396. "psubd %%mm7, %%mm1 nt"
  397. "psubd %%mm7, %%mm2 nt"
  398. "psubd %%mm7, %%mm3 nt"
  399. "psubd %%mm7, %%mm5 nt"
  400. "leal (%%esi, %%esi, 2), %%edi nt"
  402. "packssdw %%mm2, %%mm0 nt" // CcAa
  403. "packssdw %%mm3, %%mm1 nt" // DdBb
  404. "packssdw %%mm5, %%mm5 nt" // FfFf
  405. "movq %%mm0, %%mm2 nt" // CcAa
  406. "punpcklwd %%mm1, %%mm0 nt" // BAba
  407. "punpckhwd %%mm1, %%mm2 nt" // DCdc
  408. "pxor %%mm4, %%mm4 nt" // 0000
  409. "punpcklwd %%mm5, %%mm4 nt" // F0f0
  410. "movq %%mm0, %%mm1 nt" // BAba
  411. "movq %%mm4, %%mm3 nt" // F0f0
  412. "punpckldq %%mm2, %%mm0 nt" // dcba
  413. "punpckhdq %%mm1, %%mm1 nt" // BABA
  414. "punpckldq %%mm1, %%mm4 nt" // BAf0
  415. "punpckhdq %%mm3, %%mm2 nt" // F0DC
  417. "movq %%mm0, (%0, %%edi) nt"
  418. "movq %%mm4, 8(%0, %%edi) nt"
  419. "movq %%mm2, 16(%0, %%edi) nt"
  420. "addl $8, %%esi nt"
  421. " jnz 1b nt"
  422. "emms nt"
  423. :: "r" (s16+1536), "r" (f+256)
  424. :"%esi", "%edi", "memory"
  425. );
  426.     return 6*256;
  427. }
  429. static int a52_resample_3F_2R_LFE_to_6_MMX(float * _f, int16_t * s16){
  430.     int32_t * f = (int32_t *) _f;
  431. asm volatile(
  432. "movl $-1024, %%esi nt"
  433. "movq "MANGLE(magicF2W)", %%mm7 nt"
  434. // "pxor %%mm6, %%mm6 nt"
  435. "1: nt"
  436. "movq 1024(%1, %%esi), %%mm0 nt"
  437. "movq 3072(%1, %%esi), %%mm1 nt"
  438. "movq 4096(%1, %%esi), %%mm2 nt"
  439. "movq 5120(%1, %%esi), %%mm3 nt"
  440. "movq 2048(%1, %%esi), %%mm4 nt"
  441. "movq (%1, %%esi), %%mm5 nt" 
  442. "psubd %%mm7, %%mm0 nt"
  443. "psubd %%mm7, %%mm1 nt"
  444. "psubd %%mm7, %%mm2 nt"
  445. "psubd %%mm7, %%mm3 nt"
  446. "psubd %%mm7, %%mm4 nt"
  447. "psubd %%mm7, %%mm5 nt"
  448. "leal (%%esi, %%esi, 2), %%edi nt"
  450. "packssdw %%mm2, %%mm0 nt" // CcAa
  451. "packssdw %%mm3, %%mm1 nt" // DdBb
  452. "packssdw %%mm4, %%mm4 nt" // EeEe
  453. "packssdw %%mm5, %%mm5 nt" // FfFf
  454. "movq %%mm0, %%mm2 nt" // CcAa
  455. "punpcklwd %%mm1, %%mm0 nt" // BAba
  456. "punpckhwd %%mm1, %%mm2 nt" // DCdc
  457. "punpcklwd %%mm5, %%mm4 nt" // FEfe
  458. "movq %%mm0, %%mm1 nt" // BAba
  459. "movq %%mm4, %%mm3 nt" // FEfe
  460. "punpckldq %%mm2, %%mm0 nt" // dcba
  461. "punpckhdq %%mm1, %%mm1 nt" // BABA
  462. "punpckldq %%mm1, %%mm4 nt" // BAfe
  463. "punpckhdq %%mm3, %%mm2 nt" // FEDC
  465. "movq %%mm0, (%0, %%edi) nt"
  466. "movq %%mm4, 8(%0, %%edi) nt"
  467. "movq %%mm2, 16(%0, %%edi) nt"
  468. "addl $8, %%esi nt"
  469. " jnz 1b nt"
  470. "emms nt"
  471. :: "r" (s16+1536), "r" (f+256)
  472. :"%esi", "%edi", "memory"
  473. );
  474.     return 6*256;
  475. }
  478. static void* a52_resample_MMX(int flags, int ch){
  479.     switch (flags) {
  480.     case A52_MONO:
  481. if(ch==5) return a52_resample_MONO_to_5_MMX;
  482. break;
  483.     case A52_CHANNEL:
  484.     case A52_STEREO:
  485.     case A52_DOLBY:
  486. if(ch==2) return a52_resample_STEREO_to_2_MMX;
  487. break;
  488.     case A52_3F:
  489. if(ch==5) return a52_resample_3F_to_5_MMX;
  490. break;
  491.     case A52_2F2R:
  492. if(ch==4) return a52_resample_2F_2R_to_4_MMX;
  493. break;
  494.     case A52_3F2R:
  495. if(ch==5) return a52_resample_3F_2R_to_5_MMX;
  496. break;
  497.     case A52_MONO | A52_LFE:
  498. if(ch==6) return a52_resample_MONO_LFE_to_6_MMX;
  499. break;
  500.     case A52_CHANNEL | A52_LFE:
  501.     case A52_STEREO | A52_LFE:
  502.     case A52_DOLBY | A52_LFE:
  503. if(ch==6) return a52_resample_STEREO_LFE_to_6_MMX;
  504. break;
  505.     case A52_3F | A52_LFE:
  506. if(ch==6) return a52_resample_3F_LFE_to_6_MMX;
  507. break;
  508.     case A52_2F2R | A52_LFE:
  509. if(ch==6) return a52_resample_2F_2R_LFE_to_6_MMX;
  510. break;
  511.     case A52_3F2R | A52_LFE:
  512. if(ch==6) return a52_resample_3F_2R_LFE_to_6_MMX;
  513. break;
  514.     }
  515.     return NULL;
  516. }