You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

2528 lines
63KB

  1. /*
  2. *
  3. * rgb2rgb.c, Software RGB to RGB convertor
  4. * pluralize by Software PAL8 to RGB convertor
  5. * Software YUV to YUV convertor
  6. * Software YUV to RGB convertor
  7. * Written by Nick Kurshev.
  8. * palette & yuv & runtime cpu stuff by Michael (michaelni@gmx.at) (under GPL)
  9. */
  10. #include <stddef.h>
  11. #include <inttypes.h> /* for __WORDSIZE */
  12. #ifndef __WORDSIZE
  13. // #warning You have misconfigured system and probably will lose performance!
  14. #define __WORDSIZE MP_WORDSIZE
  15. #endif
  16. #undef PREFETCH
  17. #undef MOVNTQ
  18. #undef EMMS
  19. #undef SFENCE
  20. #undef MMREG_SIZE
  21. #undef PREFETCHW
  22. #undef PAVGB
  23. #ifdef HAVE_SSE2
  24. #define MMREG_SIZE 16
  25. #else
  26. #define MMREG_SIZE 8
  27. #endif
  28. #ifdef HAVE_3DNOW
  29. #define PREFETCH "prefetch"
  30. #define PREFETCHW "prefetchw"
  31. #define PAVGB "pavgusb"
  32. #elif defined ( HAVE_MMX2 )
  33. #define PREFETCH "prefetchnta"
  34. #define PREFETCHW "prefetcht0"
  35. #define PAVGB "pavgb"
  36. #else
  37. #define PREFETCH "/nop"
  38. #define PREFETCHW "/nop"
  39. #endif
  40. #ifdef HAVE_3DNOW
  41. /* On K6 femms is faster of emms. On K7 femms is directly mapped on emms. */
  42. #define EMMS "femms"
  43. #else
  44. #define EMMS "emms"
  45. #endif
  46. #ifdef HAVE_MMX2
  47. #define MOVNTQ "movntq"
  48. #define SFENCE "sfence"
  49. #else
  50. #define MOVNTQ "movq"
  51. #define SFENCE "/nop"
  52. #endif
  53. static inline void RENAME(rgb24to32)(const uint8_t *src,uint8_t *dst,unsigned src_size)
  54. {
  55. uint8_t *dest = dst;
  56. const uint8_t *s = src;
  57. const uint8_t *end;
  58. #ifdef HAVE_MMX
  59. const uint8_t *mm_end;
  60. #endif
  61. end = s + src_size;
  62. #ifdef HAVE_MMX
  63. __asm __volatile(PREFETCH" %0"::"m"(*s):"memory");
  64. mm_end = end - 23;
  65. __asm __volatile("movq %0, %%mm7"::"m"(mask32):"memory");
  66. while(s < mm_end)
  67. {
  68. __asm __volatile(
  69. PREFETCH" 32%1\n\t"
  70. "movd %1, %%mm0\n\t"
  71. "punpckldq 3%1, %%mm0\n\t"
  72. "movd 6%1, %%mm1\n\t"
  73. "punpckldq 9%1, %%mm1\n\t"
  74. "movd 12%1, %%mm2\n\t"
  75. "punpckldq 15%1, %%mm2\n\t"
  76. "movd 18%1, %%mm3\n\t"
  77. "punpckldq 21%1, %%mm3\n\t"
  78. "pand %%mm7, %%mm0\n\t"
  79. "pand %%mm7, %%mm1\n\t"
  80. "pand %%mm7, %%mm2\n\t"
  81. "pand %%mm7, %%mm3\n\t"
  82. MOVNTQ" %%mm0, %0\n\t"
  83. MOVNTQ" %%mm1, 8%0\n\t"
  84. MOVNTQ" %%mm2, 16%0\n\t"
  85. MOVNTQ" %%mm3, 24%0"
  86. :"=m"(*dest)
  87. :"m"(*s)
  88. :"memory");
  89. dest += 32;
  90. s += 24;
  91. }
  92. __asm __volatile(SFENCE:::"memory");
  93. __asm __volatile(EMMS:::"memory");
  94. #endif
  95. while(s < end)
  96. {
  97. *dest++ = *s++;
  98. *dest++ = *s++;
  99. *dest++ = *s++;
  100. *dest++ = 0;
  101. }
  102. }
  103. static inline void RENAME(rgb32to24)(const uint8_t *src,uint8_t *dst,unsigned src_size)
  104. {
  105. uint8_t *dest = dst;
  106. const uint8_t *s = src;
  107. const uint8_t *end;
  108. #ifdef HAVE_MMX
  109. const uint8_t *mm_end;
  110. #endif
  111. end = s + src_size;
  112. #ifdef HAVE_MMX
  113. __asm __volatile(PREFETCH" %0"::"m"(*s):"memory");
  114. mm_end = end - 31;
  115. while(s < mm_end)
  116. {
  117. __asm __volatile(
  118. PREFETCH" 32%1\n\t"
  119. "movq %1, %%mm0\n\t"
  120. "movq 8%1, %%mm1\n\t"
  121. "movq 16%1, %%mm4\n\t"
  122. "movq 24%1, %%mm5\n\t"
  123. "movq %%mm0, %%mm2\n\t"
  124. "movq %%mm1, %%mm3\n\t"
  125. "movq %%mm4, %%mm6\n\t"
  126. "movq %%mm5, %%mm7\n\t"
  127. "psrlq $8, %%mm2\n\t"
  128. "psrlq $8, %%mm3\n\t"
  129. "psrlq $8, %%mm6\n\t"
  130. "psrlq $8, %%mm7\n\t"
  131. "pand %2, %%mm0\n\t"
  132. "pand %2, %%mm1\n\t"
  133. "pand %2, %%mm4\n\t"
  134. "pand %2, %%mm5\n\t"
  135. "pand %3, %%mm2\n\t"
  136. "pand %3, %%mm3\n\t"
  137. "pand %3, %%mm6\n\t"
  138. "pand %3, %%mm7\n\t"
  139. "por %%mm2, %%mm0\n\t"
  140. "por %%mm3, %%mm1\n\t"
  141. "por %%mm6, %%mm4\n\t"
  142. "por %%mm7, %%mm5\n\t"
  143. "movq %%mm1, %%mm2\n\t"
  144. "movq %%mm4, %%mm3\n\t"
  145. "psllq $48, %%mm2\n\t"
  146. "psllq $32, %%mm3\n\t"
  147. "pand %4, %%mm2\n\t"
  148. "pand %5, %%mm3\n\t"
  149. "por %%mm2, %%mm0\n\t"
  150. "psrlq $16, %%mm1\n\t"
  151. "psrlq $32, %%mm4\n\t"
  152. "psllq $16, %%mm5\n\t"
  153. "por %%mm3, %%mm1\n\t"
  154. "pand %6, %%mm5\n\t"
  155. "por %%mm5, %%mm4\n\t"
  156. MOVNTQ" %%mm0, %0\n\t"
  157. MOVNTQ" %%mm1, 8%0\n\t"
  158. MOVNTQ" %%mm4, 16%0"
  159. :"=m"(*dest)
  160. :"m"(*s),"m"(mask24l),
  161. "m"(mask24h),"m"(mask24hh),"m"(mask24hhh),"m"(mask24hhhh)
  162. :"memory");
  163. dest += 24;
  164. s += 32;
  165. }
  166. __asm __volatile(SFENCE:::"memory");
  167. __asm __volatile(EMMS:::"memory");
  168. #endif
  169. while(s < end)
  170. {
  171. *dest++ = *s++;
  172. *dest++ = *s++;
  173. *dest++ = *s++;
  174. s++;
  175. }
  176. }
  177. /*
  178. Original by Strepto/Astral
  179. ported to gcc & bugfixed : A'rpi
  180. MMX2, 3DNOW optimization by Nick Kurshev
  181. 32bit c version, and and&add trick by Michael Niedermayer
  182. */
  183. static inline void RENAME(rgb15to16)(const uint8_t *src,uint8_t *dst,unsigned src_size)
  184. {
  185. register const uint8_t* s=src;
  186. register uint8_t* d=dst;
  187. register const uint8_t *end;
  188. const uint8_t *mm_end;
  189. end = s + src_size;
  190. #ifdef HAVE_MMX
  191. __asm __volatile(PREFETCH" %0"::"m"(*s));
  192. __asm __volatile("movq %0, %%mm4"::"m"(mask15s));
  193. mm_end = end - 15;
  194. while(s<mm_end)
  195. {
  196. __asm __volatile(
  197. PREFETCH" 32%1\n\t"
  198. "movq %1, %%mm0\n\t"
  199. "movq 8%1, %%mm2\n\t"
  200. "movq %%mm0, %%mm1\n\t"
  201. "movq %%mm2, %%mm3\n\t"
  202. "pand %%mm4, %%mm0\n\t"
  203. "pand %%mm4, %%mm2\n\t"
  204. "paddw %%mm1, %%mm0\n\t"
  205. "paddw %%mm3, %%mm2\n\t"
  206. MOVNTQ" %%mm0, %0\n\t"
  207. MOVNTQ" %%mm2, 8%0"
  208. :"=m"(*d)
  209. :"m"(*s)
  210. );
  211. d+=16;
  212. s+=16;
  213. }
  214. __asm __volatile(SFENCE:::"memory");
  215. __asm __volatile(EMMS:::"memory");
  216. #endif
  217. mm_end = end - 3;
  218. while(s < mm_end)
  219. {
  220. register unsigned x= *((uint32_t *)s);
  221. *((uint32_t *)d) = (x&0x7FFF7FFF) + (x&0x7FE07FE0);
  222. d+=4;
  223. s+=4;
  224. }
  225. if(s < end)
  226. {
  227. register unsigned short x= *((uint16_t *)s);
  228. *((uint16_t *)d) = (x&0x7FFF) + (x&0x7FE0);
  229. }
  230. }
  231. static inline void RENAME(bgr24torgb24)(const uint8_t *src, uint8_t *dst, unsigned src_size)
  232. {
  233. unsigned j,i,num_pixels=src_size/3;
  234. for(i=0,j=0; j<num_pixels; i+=3,j+=3)
  235. {
  236. dst[j+0] = src[i+2];
  237. dst[j+1] = src[i+1];
  238. dst[j+2] = src[i+0];
  239. }
  240. }
  241. static inline void RENAME(rgb16to15)(const uint8_t *src,uint8_t *dst,unsigned src_size)
  242. {
  243. register const uint8_t* s=src;
  244. register uint8_t* d=dst;
  245. register const uint8_t *end;
  246. const uint8_t *mm_end;
  247. end = s + src_size;
  248. #ifdef HAVE_MMX
  249. __asm __volatile(PREFETCH" %0"::"m"(*s));
  250. __asm __volatile("movq %0, %%mm7"::"m"(mask15rg));
  251. __asm __volatile("movq %0, %%mm6"::"m"(mask15b));
  252. mm_end = end - 15;
  253. while(s<mm_end)
  254. {
  255. __asm __volatile(
  256. PREFETCH" 32%1\n\t"
  257. "movq %1, %%mm0\n\t"
  258. "movq 8%1, %%mm2\n\t"
  259. "movq %%mm0, %%mm1\n\t"
  260. "movq %%mm2, %%mm3\n\t"
  261. "psrlq $1, %%mm0\n\t"
  262. "psrlq $1, %%mm2\n\t"
  263. "pand %%mm7, %%mm0\n\t"
  264. "pand %%mm7, %%mm2\n\t"
  265. "pand %%mm6, %%mm1\n\t"
  266. "pand %%mm6, %%mm3\n\t"
  267. "por %%mm1, %%mm0\n\t"
  268. "por %%mm3, %%mm2\n\t"
  269. MOVNTQ" %%mm0, %0\n\t"
  270. MOVNTQ" %%mm2, 8%0"
  271. :"=m"(*d)
  272. :"m"(*s)
  273. );
  274. d+=16;
  275. s+=16;
  276. }
  277. __asm __volatile(SFENCE:::"memory");
  278. __asm __volatile(EMMS:::"memory");
  279. #endif
  280. mm_end = end - 3;
  281. while(s < mm_end)
  282. {
  283. register uint32_t x= *((uint32_t *)s);
  284. *((uint32_t *)d) = ((x>>1)&0x7FE07FE0) | (x&0x001F001F);
  285. s+=4;
  286. d+=4;
  287. }
  288. if(s < end)
  289. {
  290. register uint16_t x= *((uint16_t *)s);
  291. *((uint16_t *)d) = ((x>>1)&0x7FE0) | (x&0x001F);
  292. s+=2;
  293. d+=2;
  294. }
  295. }
  296. static inline void RENAME(rgb32to16)(const uint8_t *src, uint8_t *dst, unsigned src_size)
  297. {
  298. const uint8_t *s = src;
  299. const uint8_t *end;
  300. #ifdef HAVE_MMX
  301. const uint8_t *mm_end;
  302. #endif
  303. uint16_t *d = (uint16_t *)dst;
  304. end = s + src_size;
  305. #ifdef HAVE_MMX
  306. mm_end = end - 15;
  307. #if 1 //is faster only if multiplies are reasonable fast (FIXME figure out on which cpus this is faster, on Athlon its slightly faster)
  308. asm volatile(
  309. "movq %3, %%mm5 \n\t"
  310. "movq %4, %%mm6 \n\t"
  311. "movq %5, %%mm7 \n\t"
  312. ".balign 16 \n\t"
  313. "1: \n\t"
  314. PREFETCH" 32(%1) \n\t"
  315. "movd (%1), %%mm0 \n\t"
  316. "movd 4(%1), %%mm3 \n\t"
  317. "punpckldq 8(%1), %%mm0 \n\t"
  318. "punpckldq 12(%1), %%mm3 \n\t"
  319. "movq %%mm0, %%mm1 \n\t"
  320. "movq %%mm3, %%mm4 \n\t"
  321. "pand %%mm6, %%mm0 \n\t"
  322. "pand %%mm6, %%mm3 \n\t"
  323. "pmaddwd %%mm7, %%mm0 \n\t"
  324. "pmaddwd %%mm7, %%mm3 \n\t"
  325. "pand %%mm5, %%mm1 \n\t"
  326. "pand %%mm5, %%mm4 \n\t"
  327. "por %%mm1, %%mm0 \n\t"
  328. "por %%mm4, %%mm3 \n\t"
  329. "psrld $5, %%mm0 \n\t"
  330. "pslld $11, %%mm3 \n\t"
  331. "por %%mm3, %%mm0 \n\t"
  332. MOVNTQ" %%mm0, (%0) \n\t"
  333. "addl $16, %1 \n\t"
  334. "addl $8, %0 \n\t"
  335. "cmpl %2, %1 \n\t"
  336. " jb 1b \n\t"
  337. : "+r" (d), "+r"(s)
  338. : "r" (mm_end), "m" (mask3216g), "m" (mask3216br), "m" (mul3216)
  339. );
  340. #else
  341. __asm __volatile(PREFETCH" %0"::"m"(*src):"memory");
  342. __asm __volatile(
  343. "movq %0, %%mm7\n\t"
  344. "movq %1, %%mm6\n\t"
  345. ::"m"(red_16mask),"m"(green_16mask));
  346. while(s < mm_end)
  347. {
  348. __asm __volatile(
  349. PREFETCH" 32%1\n\t"
  350. "movd %1, %%mm0\n\t"
  351. "movd 4%1, %%mm3\n\t"
  352. "punpckldq 8%1, %%mm0\n\t"
  353. "punpckldq 12%1, %%mm3\n\t"
  354. "movq %%mm0, %%mm1\n\t"
  355. "movq %%mm0, %%mm2\n\t"
  356. "movq %%mm3, %%mm4\n\t"
  357. "movq %%mm3, %%mm5\n\t"
  358. "psrlq $3, %%mm0\n\t"
  359. "psrlq $3, %%mm3\n\t"
  360. "pand %2, %%mm0\n\t"
  361. "pand %2, %%mm3\n\t"
  362. "psrlq $5, %%mm1\n\t"
  363. "psrlq $5, %%mm4\n\t"
  364. "pand %%mm6, %%mm1\n\t"
  365. "pand %%mm6, %%mm4\n\t"
  366. "psrlq $8, %%mm2\n\t"
  367. "psrlq $8, %%mm5\n\t"
  368. "pand %%mm7, %%mm2\n\t"
  369. "pand %%mm7, %%mm5\n\t"
  370. "por %%mm1, %%mm0\n\t"
  371. "por %%mm4, %%mm3\n\t"
  372. "por %%mm2, %%mm0\n\t"
  373. "por %%mm5, %%mm3\n\t"
  374. "psllq $16, %%mm3\n\t"
  375. "por %%mm3, %%mm0\n\t"
  376. MOVNTQ" %%mm0, %0\n\t"
  377. :"=m"(*d):"m"(*s),"m"(blue_16mask):"memory");
  378. d += 4;
  379. s += 16;
  380. }
  381. #endif
  382. __asm __volatile(SFENCE:::"memory");
  383. __asm __volatile(EMMS:::"memory");
  384. #endif
  385. while(s < end)
  386. {
  387. const int src= *((uint32_t*)s)++;
  388. *d++ = ((src&0xFF)>>3) + ((src&0xFC00)>>5) + ((src&0xF80000)>>8);
  389. // *d++ = ((src>>3)&0x1F) + ((src>>5)&0x7E0) + ((src>>8)&0xF800);
  390. }
  391. }
  392. static inline void RENAME(rgb32tobgr16)(const uint8_t *src, uint8_t *dst, unsigned int src_size)
  393. {
  394. const uint8_t *s = src;
  395. const uint8_t *end;
  396. #ifdef HAVE_MMX
  397. const uint8_t *mm_end;
  398. #endif
  399. uint16_t *d = (uint16_t *)dst;
  400. end = s + src_size;
  401. #ifdef HAVE_MMX
  402. __asm __volatile(PREFETCH" %0"::"m"(*src):"memory");
  403. __asm __volatile(
  404. "movq %0, %%mm7\n\t"
  405. "movq %1, %%mm6\n\t"
  406. ::"m"(red_16mask),"m"(green_16mask));
  407. mm_end = end - 15;
  408. while(s < mm_end)
  409. {
  410. __asm __volatile(
  411. PREFETCH" 32%1\n\t"
  412. "movd %1, %%mm0\n\t"
  413. "movd 4%1, %%mm3\n\t"
  414. "punpckldq 8%1, %%mm0\n\t"
  415. "punpckldq 12%1, %%mm3\n\t"
  416. "movq %%mm0, %%mm1\n\t"
  417. "movq %%mm0, %%mm2\n\t"
  418. "movq %%mm3, %%mm4\n\t"
  419. "movq %%mm3, %%mm5\n\t"
  420. "psllq $8, %%mm0\n\t"
  421. "psllq $8, %%mm3\n\t"
  422. "pand %%mm7, %%mm0\n\t"
  423. "pand %%mm7, %%mm3\n\t"
  424. "psrlq $5, %%mm1\n\t"
  425. "psrlq $5, %%mm4\n\t"
  426. "pand %%mm6, %%mm1\n\t"
  427. "pand %%mm6, %%mm4\n\t"
  428. "psrlq $19, %%mm2\n\t"
  429. "psrlq $19, %%mm5\n\t"
  430. "pand %2, %%mm2\n\t"
  431. "pand %2, %%mm5\n\t"
  432. "por %%mm1, %%mm0\n\t"
  433. "por %%mm4, %%mm3\n\t"
  434. "por %%mm2, %%mm0\n\t"
  435. "por %%mm5, %%mm3\n\t"
  436. "psllq $16, %%mm3\n\t"
  437. "por %%mm3, %%mm0\n\t"
  438. MOVNTQ" %%mm0, %0\n\t"
  439. :"=m"(*d):"m"(*s),"m"(blue_16mask):"memory");
  440. d += 4;
  441. s += 16;
  442. }
  443. __asm __volatile(SFENCE:::"memory");
  444. __asm __volatile(EMMS:::"memory");
  445. #endif
  446. while(s < end)
  447. {
  448. const int src= *((uint32_t*)s)++;
  449. *d++ = ((src&0xF8)<<8) + ((src&0xFC00)>>5) + ((src&0xF80000)>>19);
  450. }
  451. }
  452. static inline void RENAME(rgb32to15)(const uint8_t *src, uint8_t *dst, unsigned src_size)
  453. {
  454. const uint8_t *s = src;
  455. const uint8_t *end;
  456. #ifdef HAVE_MMX
  457. const uint8_t *mm_end;
  458. #endif
  459. uint16_t *d = (uint16_t *)dst;
  460. end = s + src_size;
  461. #ifdef HAVE_MMX
  462. mm_end = end - 15;
  463. #if 1 //is faster only if multiplies are reasonable fast (FIXME figure out on which cpus this is faster, on Athlon its slightly faster)
  464. asm volatile(
  465. "movq %3, %%mm5 \n\t"
  466. "movq %4, %%mm6 \n\t"
  467. "movq %5, %%mm7 \n\t"
  468. ".balign 16 \n\t"
  469. "1: \n\t"
  470. PREFETCH" 32(%1) \n\t"
  471. "movd (%1), %%mm0 \n\t"
  472. "movd 4(%1), %%mm3 \n\t"
  473. "punpckldq 8(%1), %%mm0 \n\t"
  474. "punpckldq 12(%1), %%mm3 \n\t"
  475. "movq %%mm0, %%mm1 \n\t"
  476. "movq %%mm3, %%mm4 \n\t"
  477. "pand %%mm6, %%mm0 \n\t"
  478. "pand %%mm6, %%mm3 \n\t"
  479. "pmaddwd %%mm7, %%mm0 \n\t"
  480. "pmaddwd %%mm7, %%mm3 \n\t"
  481. "pand %%mm5, %%mm1 \n\t"
  482. "pand %%mm5, %%mm4 \n\t"
  483. "por %%mm1, %%mm0 \n\t"
  484. "por %%mm4, %%mm3 \n\t"
  485. "psrld $6, %%mm0 \n\t"
  486. "pslld $10, %%mm3 \n\t"
  487. "por %%mm3, %%mm0 \n\t"
  488. MOVNTQ" %%mm0, (%0) \n\t"
  489. "addl $16, %1 \n\t"
  490. "addl $8, %0 \n\t"
  491. "cmpl %2, %1 \n\t"
  492. " jb 1b \n\t"
  493. : "+r" (d), "+r"(s)
  494. : "r" (mm_end), "m" (mask3215g), "m" (mask3216br), "m" (mul3215)
  495. );
  496. #else
  497. __asm __volatile(PREFETCH" %0"::"m"(*src):"memory");
  498. __asm __volatile(
  499. "movq %0, %%mm7\n\t"
  500. "movq %1, %%mm6\n\t"
  501. ::"m"(red_15mask),"m"(green_15mask));
  502. while(s < mm_end)
  503. {
  504. __asm __volatile(
  505. PREFETCH" 32%1\n\t"
  506. "movd %1, %%mm0\n\t"
  507. "movd 4%1, %%mm3\n\t"
  508. "punpckldq 8%1, %%mm0\n\t"
  509. "punpckldq 12%1, %%mm3\n\t"
  510. "movq %%mm0, %%mm1\n\t"
  511. "movq %%mm0, %%mm2\n\t"
  512. "movq %%mm3, %%mm4\n\t"
  513. "movq %%mm3, %%mm5\n\t"
  514. "psrlq $3, %%mm0\n\t"
  515. "psrlq $3, %%mm3\n\t"
  516. "pand %2, %%mm0\n\t"
  517. "pand %2, %%mm3\n\t"
  518. "psrlq $6, %%mm1\n\t"
  519. "psrlq $6, %%mm4\n\t"
  520. "pand %%mm6, %%mm1\n\t"
  521. "pand %%mm6, %%mm4\n\t"
  522. "psrlq $9, %%mm2\n\t"
  523. "psrlq $9, %%mm5\n\t"
  524. "pand %%mm7, %%mm2\n\t"
  525. "pand %%mm7, %%mm5\n\t"
  526. "por %%mm1, %%mm0\n\t"
  527. "por %%mm4, %%mm3\n\t"
  528. "por %%mm2, %%mm0\n\t"
  529. "por %%mm5, %%mm3\n\t"
  530. "psllq $16, %%mm3\n\t"
  531. "por %%mm3, %%mm0\n\t"
  532. MOVNTQ" %%mm0, %0\n\t"
  533. :"=m"(*d):"m"(*s),"m"(blue_15mask):"memory");
  534. d += 4;
  535. s += 16;
  536. }
  537. #endif
  538. __asm __volatile(SFENCE:::"memory");
  539. __asm __volatile(EMMS:::"memory");
  540. #endif
  541. while(s < end)
  542. {
  543. const int src= *((uint32_t*)s)++;
  544. *d++ = ((src&0xFF)>>3) + ((src&0xF800)>>6) + ((src&0xF80000)>>9);
  545. }
  546. }
  547. static inline void RENAME(rgb32tobgr15)(const uint8_t *src, uint8_t *dst, unsigned src_size)
  548. {
  549. const uint8_t *s = src;
  550. const uint8_t *end;
  551. #ifdef HAVE_MMX
  552. const uint8_t *mm_end;
  553. #endif
  554. uint16_t *d = (uint16_t *)dst;
  555. end = s + src_size;
  556. #ifdef HAVE_MMX
  557. __asm __volatile(PREFETCH" %0"::"m"(*src):"memory");
  558. __asm __volatile(
  559. "movq %0, %%mm7\n\t"
  560. "movq %1, %%mm6\n\t"
  561. ::"m"(red_15mask),"m"(green_15mask));
  562. mm_end = end - 15;
  563. while(s < mm_end)
  564. {
  565. __asm __volatile(
  566. PREFETCH" 32%1\n\t"
  567. "movd %1, %%mm0\n\t"
  568. "movd 4%1, %%mm3\n\t"
  569. "punpckldq 8%1, %%mm0\n\t"
  570. "punpckldq 12%1, %%mm3\n\t"
  571. "movq %%mm0, %%mm1\n\t"
  572. "movq %%mm0, %%mm2\n\t"
  573. "movq %%mm3, %%mm4\n\t"
  574. "movq %%mm3, %%mm5\n\t"
  575. "psllq $7, %%mm0\n\t"
  576. "psllq $7, %%mm3\n\t"
  577. "pand %%mm7, %%mm0\n\t"
  578. "pand %%mm7, %%mm3\n\t"
  579. "psrlq $6, %%mm1\n\t"
  580. "psrlq $6, %%mm4\n\t"
  581. "pand %%mm6, %%mm1\n\t"
  582. "pand %%mm6, %%mm4\n\t"
  583. "psrlq $19, %%mm2\n\t"
  584. "psrlq $19, %%mm5\n\t"
  585. "pand %2, %%mm2\n\t"
  586. "pand %2, %%mm5\n\t"
  587. "por %%mm1, %%mm0\n\t"
  588. "por %%mm4, %%mm3\n\t"
  589. "por %%mm2, %%mm0\n\t"
  590. "por %%mm5, %%mm3\n\t"
  591. "psllq $16, %%mm3\n\t"
  592. "por %%mm3, %%mm0\n\t"
  593. MOVNTQ" %%mm0, %0\n\t"
  594. :"=m"(*d):"m"(*s),"m"(blue_15mask):"memory");
  595. d += 4;
  596. s += 16;
  597. }
  598. __asm __volatile(SFENCE:::"memory");
  599. __asm __volatile(EMMS:::"memory");
  600. #endif
  601. while(s < end)
  602. {
  603. const int src= *((uint32_t*)s)++;
  604. *d++ = ((src&0xF8)<<7) + ((src&0xF800)>>6) + ((src&0xF80000)>>19);
  605. }
  606. }
  607. static inline void RENAME(rgb24to16)(const uint8_t *src, uint8_t *dst, unsigned src_size)
  608. {
  609. const uint8_t *s = src;
  610. const uint8_t *end;
  611. #ifdef HAVE_MMX
  612. const uint8_t *mm_end;
  613. #endif
  614. uint16_t *d = (uint16_t *)dst;
  615. end = s + src_size;
  616. #ifdef HAVE_MMX
  617. __asm __volatile(PREFETCH" %0"::"m"(*src):"memory");
  618. __asm __volatile(
  619. "movq %0, %%mm7\n\t"
  620. "movq %1, %%mm6\n\t"
  621. ::"m"(red_16mask),"m"(green_16mask));
  622. mm_end = end - 11;
  623. while(s < mm_end)
  624. {
  625. __asm __volatile(
  626. PREFETCH" 32%1\n\t"
  627. "movd %1, %%mm0\n\t"
  628. "movd 3%1, %%mm3\n\t"
  629. "punpckldq 6%1, %%mm0\n\t"
  630. "punpckldq 9%1, %%mm3\n\t"
  631. "movq %%mm0, %%mm1\n\t"
  632. "movq %%mm0, %%mm2\n\t"
  633. "movq %%mm3, %%mm4\n\t"
  634. "movq %%mm3, %%mm5\n\t"
  635. "psrlq $3, %%mm0\n\t"
  636. "psrlq $3, %%mm3\n\t"
  637. "pand %2, %%mm0\n\t"
  638. "pand %2, %%mm3\n\t"
  639. "psrlq $5, %%mm1\n\t"
  640. "psrlq $5, %%mm4\n\t"
  641. "pand %%mm6, %%mm1\n\t"
  642. "pand %%mm6, %%mm4\n\t"
  643. "psrlq $8, %%mm2\n\t"
  644. "psrlq $8, %%mm5\n\t"
  645. "pand %%mm7, %%mm2\n\t"
  646. "pand %%mm7, %%mm5\n\t"
  647. "por %%mm1, %%mm0\n\t"
  648. "por %%mm4, %%mm3\n\t"
  649. "por %%mm2, %%mm0\n\t"
  650. "por %%mm5, %%mm3\n\t"
  651. "psllq $16, %%mm3\n\t"
  652. "por %%mm3, %%mm0\n\t"
  653. MOVNTQ" %%mm0, %0\n\t"
  654. :"=m"(*d):"m"(*s),"m"(blue_16mask):"memory");
  655. d += 4;
  656. s += 12;
  657. }
  658. __asm __volatile(SFENCE:::"memory");
  659. __asm __volatile(EMMS:::"memory");
  660. #endif
  661. while(s < end)
  662. {
  663. const int b= *s++;
  664. const int g= *s++;
  665. const int r= *s++;
  666. *d++ = (b>>3) | ((g&0xFC)<<3) | ((r&0xF8)<<8);
  667. }
  668. }
  669. static inline void RENAME(rgb24tobgr16)(const uint8_t *src, uint8_t *dst, unsigned int src_size)
  670. {
  671. const uint8_t *s = src;
  672. const uint8_t *end;
  673. #ifdef HAVE_MMX
  674. const uint8_t *mm_end;
  675. #endif
  676. uint16_t *d = (uint16_t *)dst;
  677. end = s + src_size;
  678. #ifdef HAVE_MMX
  679. __asm __volatile(PREFETCH" %0"::"m"(*src):"memory");
  680. __asm __volatile(
  681. "movq %0, %%mm7\n\t"
  682. "movq %1, %%mm6\n\t"
  683. ::"m"(red_16mask),"m"(green_16mask));
  684. mm_end = end - 15;
  685. while(s < mm_end)
  686. {
  687. __asm __volatile(
  688. PREFETCH" 32%1\n\t"
  689. "movd %1, %%mm0\n\t"
  690. "movd 3%1, %%mm3\n\t"
  691. "punpckldq 6%1, %%mm0\n\t"
  692. "punpckldq 9%1, %%mm3\n\t"
  693. "movq %%mm0, %%mm1\n\t"
  694. "movq %%mm0, %%mm2\n\t"
  695. "movq %%mm3, %%mm4\n\t"
  696. "movq %%mm3, %%mm5\n\t"
  697. "psllq $8, %%mm0\n\t"
  698. "psllq $8, %%mm3\n\t"
  699. "pand %%mm7, %%mm0\n\t"
  700. "pand %%mm7, %%mm3\n\t"
  701. "psrlq $5, %%mm1\n\t"
  702. "psrlq $5, %%mm4\n\t"
  703. "pand %%mm6, %%mm1\n\t"
  704. "pand %%mm6, %%mm4\n\t"
  705. "psrlq $19, %%mm2\n\t"
  706. "psrlq $19, %%mm5\n\t"
  707. "pand %2, %%mm2\n\t"
  708. "pand %2, %%mm5\n\t"
  709. "por %%mm1, %%mm0\n\t"
  710. "por %%mm4, %%mm3\n\t"
  711. "por %%mm2, %%mm0\n\t"
  712. "por %%mm5, %%mm3\n\t"
  713. "psllq $16, %%mm3\n\t"
  714. "por %%mm3, %%mm0\n\t"
  715. MOVNTQ" %%mm0, %0\n\t"
  716. :"=m"(*d):"m"(*s),"m"(blue_16mask):"memory");
  717. d += 4;
  718. s += 12;
  719. }
  720. __asm __volatile(SFENCE:::"memory");
  721. __asm __volatile(EMMS:::"memory");
  722. #endif
  723. while(s < end)
  724. {
  725. const int r= *s++;
  726. const int g= *s++;
  727. const int b= *s++;
  728. *d++ = (b>>3) | ((g&0xFC)<<3) | ((r&0xF8)<<8);
  729. }
  730. }
  731. static inline void RENAME(rgb24to15)(const uint8_t *src, uint8_t *dst, unsigned src_size)
  732. {
  733. const uint8_t *s = src;
  734. const uint8_t *end;
  735. #ifdef HAVE_MMX
  736. const uint8_t *mm_end;
  737. #endif
  738. uint16_t *d = (uint16_t *)dst;
  739. end = s + src_size;
  740. #ifdef HAVE_MMX
  741. __asm __volatile(PREFETCH" %0"::"m"(*src):"memory");
  742. __asm __volatile(
  743. "movq %0, %%mm7\n\t"
  744. "movq %1, %%mm6\n\t"
  745. ::"m"(red_15mask),"m"(green_15mask));
  746. mm_end = end - 11;
  747. while(s < mm_end)
  748. {
  749. __asm __volatile(
  750. PREFETCH" 32%1\n\t"
  751. "movd %1, %%mm0\n\t"
  752. "movd 3%1, %%mm3\n\t"
  753. "punpckldq 6%1, %%mm0\n\t"
  754. "punpckldq 9%1, %%mm3\n\t"
  755. "movq %%mm0, %%mm1\n\t"
  756. "movq %%mm0, %%mm2\n\t"
  757. "movq %%mm3, %%mm4\n\t"
  758. "movq %%mm3, %%mm5\n\t"
  759. "psrlq $3, %%mm0\n\t"
  760. "psrlq $3, %%mm3\n\t"
  761. "pand %2, %%mm0\n\t"
  762. "pand %2, %%mm3\n\t"
  763. "psrlq $6, %%mm1\n\t"
  764. "psrlq $6, %%mm4\n\t"
  765. "pand %%mm6, %%mm1\n\t"
  766. "pand %%mm6, %%mm4\n\t"
  767. "psrlq $9, %%mm2\n\t"
  768. "psrlq $9, %%mm5\n\t"
  769. "pand %%mm7, %%mm2\n\t"
  770. "pand %%mm7, %%mm5\n\t"
  771. "por %%mm1, %%mm0\n\t"
  772. "por %%mm4, %%mm3\n\t"
  773. "por %%mm2, %%mm0\n\t"
  774. "por %%mm5, %%mm3\n\t"
  775. "psllq $16, %%mm3\n\t"
  776. "por %%mm3, %%mm0\n\t"
  777. MOVNTQ" %%mm0, %0\n\t"
  778. :"=m"(*d):"m"(*s),"m"(blue_15mask):"memory");
  779. d += 4;
  780. s += 12;
  781. }
  782. __asm __volatile(SFENCE:::"memory");
  783. __asm __volatile(EMMS:::"memory");
  784. #endif
  785. while(s < end)
  786. {
  787. const int b= *s++;
  788. const int g= *s++;
  789. const int r= *s++;
  790. *d++ = (b>>3) | ((g&0xF8)<<2) | ((r&0xF8)<<7);
  791. }
  792. }
  793. static inline void RENAME(rgb24tobgr15)(const uint8_t *src, uint8_t *dst, unsigned src_size)
  794. {
  795. const uint8_t *s = src;
  796. const uint8_t *end;
  797. #ifdef HAVE_MMX
  798. const uint8_t *mm_end;
  799. #endif
  800. uint16_t *d = (uint16_t *)dst;
  801. end = s + src_size;
  802. #ifdef HAVE_MMX
  803. __asm __volatile(PREFETCH" %0"::"m"(*src):"memory");
  804. __asm __volatile(
  805. "movq %0, %%mm7\n\t"
  806. "movq %1, %%mm6\n\t"
  807. ::"m"(red_15mask),"m"(green_15mask));
  808. mm_end = end - 15;
  809. while(s < mm_end)
  810. {
  811. __asm __volatile(
  812. PREFETCH" 32%1\n\t"
  813. "movd %1, %%mm0\n\t"
  814. "movd 3%1, %%mm3\n\t"
  815. "punpckldq 6%1, %%mm0\n\t"
  816. "punpckldq 9%1, %%mm3\n\t"
  817. "movq %%mm0, %%mm1\n\t"
  818. "movq %%mm0, %%mm2\n\t"
  819. "movq %%mm3, %%mm4\n\t"
  820. "movq %%mm3, %%mm5\n\t"
  821. "psllq $7, %%mm0\n\t"
  822. "psllq $7, %%mm3\n\t"
  823. "pand %%mm7, %%mm0\n\t"
  824. "pand %%mm7, %%mm3\n\t"
  825. "psrlq $6, %%mm1\n\t"
  826. "psrlq $6, %%mm4\n\t"
  827. "pand %%mm6, %%mm1\n\t"
  828. "pand %%mm6, %%mm4\n\t"
  829. "psrlq $19, %%mm2\n\t"
  830. "psrlq $19, %%mm5\n\t"
  831. "pand %2, %%mm2\n\t"
  832. "pand %2, %%mm5\n\t"
  833. "por %%mm1, %%mm0\n\t"
  834. "por %%mm4, %%mm3\n\t"
  835. "por %%mm2, %%mm0\n\t"
  836. "por %%mm5, %%mm3\n\t"
  837. "psllq $16, %%mm3\n\t"
  838. "por %%mm3, %%mm0\n\t"
  839. MOVNTQ" %%mm0, %0\n\t"
  840. :"=m"(*d):"m"(*s),"m"(blue_15mask):"memory");
  841. d += 4;
  842. s += 12;
  843. }
  844. __asm __volatile(SFENCE:::"memory");
  845. __asm __volatile(EMMS:::"memory");
  846. #endif
  847. while(s < end)
  848. {
  849. const int r= *s++;
  850. const int g= *s++;
  851. const int b= *s++;
  852. *d++ = (b>>3) | ((g&0xF8)<<2) | ((r&0xF8)<<7);
  853. }
  854. }
  855. /*
  856. I use here less accurate approximation by simply
  857. left-shifting the input
  858. value and filling the low order bits with
  859. zeroes. This method improves png's
  860. compression but this scheme cannot reproduce white exactly, since it does not
  861. generate an all-ones maximum value; the net effect is to darken the
  862. image slightly.
  863. The better method should be "left bit replication":
  864. 4 3 2 1 0
  865. ---------
  866. 1 1 0 1 1
  867. 7 6 5 4 3 2 1 0
  868. ----------------
  869. 1 1 0 1 1 1 1 0
  870. |=======| |===|
  871. | Leftmost Bits Repeated to Fill Open Bits
  872. |
  873. Original Bits
  874. */
  875. static inline void RENAME(rgb15to24)(const uint8_t *src, uint8_t *dst, unsigned src_size)
  876. {
  877. const uint16_t *end;
  878. #ifdef HAVE_MMX
  879. const uint16_t *mm_end;
  880. #endif
  881. uint8_t *d = (uint8_t *)dst;
  882. const uint16_t *s = (uint16_t *)src;
  883. end = s + src_size/2;
  884. #ifdef HAVE_MMX
  885. __asm __volatile(PREFETCH" %0"::"m"(*s):"memory");
  886. mm_end = end - 7;
  887. while(s < mm_end)
  888. {
  889. __asm __volatile(
  890. PREFETCH" 32%1\n\t"
  891. "movq %1, %%mm0\n\t"
  892. "movq %1, %%mm1\n\t"
  893. "movq %1, %%mm2\n\t"
  894. "pand %2, %%mm0\n\t"
  895. "pand %3, %%mm1\n\t"
  896. "pand %4, %%mm2\n\t"
  897. "psllq $3, %%mm0\n\t"
  898. "psrlq $2, %%mm1\n\t"
  899. "psrlq $7, %%mm2\n\t"
  900. "movq %%mm0, %%mm3\n\t"
  901. "movq %%mm1, %%mm4\n\t"
  902. "movq %%mm2, %%mm5\n\t"
  903. "punpcklwd %5, %%mm0\n\t"
  904. "punpcklwd %5, %%mm1\n\t"
  905. "punpcklwd %5, %%mm2\n\t"
  906. "punpckhwd %5, %%mm3\n\t"
  907. "punpckhwd %5, %%mm4\n\t"
  908. "punpckhwd %5, %%mm5\n\t"
  909. "psllq $8, %%mm1\n\t"
  910. "psllq $16, %%mm2\n\t"
  911. "por %%mm1, %%mm0\n\t"
  912. "por %%mm2, %%mm0\n\t"
  913. "psllq $8, %%mm4\n\t"
  914. "psllq $16, %%mm5\n\t"
  915. "por %%mm4, %%mm3\n\t"
  916. "por %%mm5, %%mm3\n\t"
  917. "movq %%mm0, %%mm6\n\t"
  918. "movq %%mm3, %%mm7\n\t"
  919. "movq 8%1, %%mm0\n\t"
  920. "movq 8%1, %%mm1\n\t"
  921. "movq 8%1, %%mm2\n\t"
  922. "pand %2, %%mm0\n\t"
  923. "pand %3, %%mm1\n\t"
  924. "pand %4, %%mm2\n\t"
  925. "psllq $3, %%mm0\n\t"
  926. "psrlq $2, %%mm1\n\t"
  927. "psrlq $7, %%mm2\n\t"
  928. "movq %%mm0, %%mm3\n\t"
  929. "movq %%mm1, %%mm4\n\t"
  930. "movq %%mm2, %%mm5\n\t"
  931. "punpcklwd %5, %%mm0\n\t"
  932. "punpcklwd %5, %%mm1\n\t"
  933. "punpcklwd %5, %%mm2\n\t"
  934. "punpckhwd %5, %%mm3\n\t"
  935. "punpckhwd %5, %%mm4\n\t"
  936. "punpckhwd %5, %%mm5\n\t"
  937. "psllq $8, %%mm1\n\t"
  938. "psllq $16, %%mm2\n\t"
  939. "por %%mm1, %%mm0\n\t"
  940. "por %%mm2, %%mm0\n\t"
  941. "psllq $8, %%mm4\n\t"
  942. "psllq $16, %%mm5\n\t"
  943. "por %%mm4, %%mm3\n\t"
  944. "por %%mm5, %%mm3\n\t"
  945. :"=m"(*d)
  946. :"m"(*s),"m"(mask15b),"m"(mask15g),"m"(mask15r), "m"(mmx_null)
  947. :"memory");
  948. /* Borrowed 32 to 24 */
  949. __asm __volatile(
  950. "movq %%mm0, %%mm4\n\t"
  951. "movq %%mm3, %%mm5\n\t"
  952. "movq %%mm6, %%mm0\n\t"
  953. "movq %%mm7, %%mm1\n\t"
  954. "movq %%mm4, %%mm6\n\t"
  955. "movq %%mm5, %%mm7\n\t"
  956. "movq %%mm0, %%mm2\n\t"
  957. "movq %%mm1, %%mm3\n\t"
  958. "psrlq $8, %%mm2\n\t"
  959. "psrlq $8, %%mm3\n\t"
  960. "psrlq $8, %%mm6\n\t"
  961. "psrlq $8, %%mm7\n\t"
  962. "pand %2, %%mm0\n\t"
  963. "pand %2, %%mm1\n\t"
  964. "pand %2, %%mm4\n\t"
  965. "pand %2, %%mm5\n\t"
  966. "pand %3, %%mm2\n\t"
  967. "pand %3, %%mm3\n\t"
  968. "pand %3, %%mm6\n\t"
  969. "pand %3, %%mm7\n\t"
  970. "por %%mm2, %%mm0\n\t"
  971. "por %%mm3, %%mm1\n\t"
  972. "por %%mm6, %%mm4\n\t"
  973. "por %%mm7, %%mm5\n\t"
  974. "movq %%mm1, %%mm2\n\t"
  975. "movq %%mm4, %%mm3\n\t"
  976. "psllq $48, %%mm2\n\t"
  977. "psllq $32, %%mm3\n\t"
  978. "pand %4, %%mm2\n\t"
  979. "pand %5, %%mm3\n\t"
  980. "por %%mm2, %%mm0\n\t"
  981. "psrlq $16, %%mm1\n\t"
  982. "psrlq $32, %%mm4\n\t"
  983. "psllq $16, %%mm5\n\t"
  984. "por %%mm3, %%mm1\n\t"
  985. "pand %6, %%mm5\n\t"
  986. "por %%mm5, %%mm4\n\t"
  987. MOVNTQ" %%mm0, %0\n\t"
  988. MOVNTQ" %%mm1, 8%0\n\t"
  989. MOVNTQ" %%mm4, 16%0"
  990. :"=m"(*d)
  991. :"m"(*s),"m"(mask24l),"m"(mask24h),"m"(mask24hh),"m"(mask24hhh),"m"(mask24hhhh)
  992. :"memory");
  993. d += 24;
  994. s += 8;
  995. }
  996. __asm __volatile(SFENCE:::"memory");
  997. __asm __volatile(EMMS:::"memory");
  998. #endif
  999. while(s < end)
  1000. {
  1001. register uint16_t bgr;
  1002. bgr = *s++;
  1003. *d++ = (bgr&0x1F)<<3;
  1004. *d++ = (bgr&0x3E0)>>2;
  1005. *d++ = (bgr&0x7C00)>>7;
  1006. }
  1007. }
  1008. static inline void RENAME(rgb16to24)(const uint8_t *src, uint8_t *dst, unsigned src_size)
  1009. {
  1010. const uint16_t *end;
  1011. #ifdef HAVE_MMX
  1012. const uint16_t *mm_end;
  1013. #endif
  1014. uint8_t *d = (uint8_t *)dst;
  1015. const uint16_t *s = (const uint16_t *)src;
  1016. end = s + src_size/2;
  1017. #ifdef HAVE_MMX
  1018. __asm __volatile(PREFETCH" %0"::"m"(*s):"memory");
  1019. mm_end = end - 7;
  1020. while(s < mm_end)
  1021. {
  1022. __asm __volatile(
  1023. PREFETCH" 32%1\n\t"
  1024. "movq %1, %%mm0\n\t"
  1025. "movq %1, %%mm1\n\t"
  1026. "movq %1, %%mm2\n\t"
  1027. "pand %2, %%mm0\n\t"
  1028. "pand %3, %%mm1\n\t"
  1029. "pand %4, %%mm2\n\t"
  1030. "psllq $3, %%mm0\n\t"
  1031. "psrlq $3, %%mm1\n\t"
  1032. "psrlq $8, %%mm2\n\t"
  1033. "movq %%mm0, %%mm3\n\t"
  1034. "movq %%mm1, %%mm4\n\t"
  1035. "movq %%mm2, %%mm5\n\t"
  1036. "punpcklwd %5, %%mm0\n\t"
  1037. "punpcklwd %5, %%mm1\n\t"
  1038. "punpcklwd %5, %%mm2\n\t"
  1039. "punpckhwd %5, %%mm3\n\t"
  1040. "punpckhwd %5, %%mm4\n\t"
  1041. "punpckhwd %5, %%mm5\n\t"
  1042. "psllq $8, %%mm1\n\t"
  1043. "psllq $16, %%mm2\n\t"
  1044. "por %%mm1, %%mm0\n\t"
  1045. "por %%mm2, %%mm0\n\t"
  1046. "psllq $8, %%mm4\n\t"
  1047. "psllq $16, %%mm5\n\t"
  1048. "por %%mm4, %%mm3\n\t"
  1049. "por %%mm5, %%mm3\n\t"
  1050. "movq %%mm0, %%mm6\n\t"
  1051. "movq %%mm3, %%mm7\n\t"
  1052. "movq 8%1, %%mm0\n\t"
  1053. "movq 8%1, %%mm1\n\t"
  1054. "movq 8%1, %%mm2\n\t"
  1055. "pand %2, %%mm0\n\t"
  1056. "pand %3, %%mm1\n\t"
  1057. "pand %4, %%mm2\n\t"
  1058. "psllq $3, %%mm0\n\t"
  1059. "psrlq $3, %%mm1\n\t"
  1060. "psrlq $8, %%mm2\n\t"
  1061. "movq %%mm0, %%mm3\n\t"
  1062. "movq %%mm1, %%mm4\n\t"
  1063. "movq %%mm2, %%mm5\n\t"
  1064. "punpcklwd %5, %%mm0\n\t"
  1065. "punpcklwd %5, %%mm1\n\t"
  1066. "punpcklwd %5, %%mm2\n\t"
  1067. "punpckhwd %5, %%mm3\n\t"
  1068. "punpckhwd %5, %%mm4\n\t"
  1069. "punpckhwd %5, %%mm5\n\t"
  1070. "psllq $8, %%mm1\n\t"
  1071. "psllq $16, %%mm2\n\t"
  1072. "por %%mm1, %%mm0\n\t"
  1073. "por %%mm2, %%mm0\n\t"
  1074. "psllq $8, %%mm4\n\t"
  1075. "psllq $16, %%mm5\n\t"
  1076. "por %%mm4, %%mm3\n\t"
  1077. "por %%mm5, %%mm3\n\t"
  1078. :"=m"(*d)
  1079. :"m"(*s),"m"(mask16b),"m"(mask16g),"m"(mask16r),"m"(mmx_null)
  1080. :"memory");
  1081. /* Borrowed 32 to 24 */
  1082. __asm __volatile(
  1083. "movq %%mm0, %%mm4\n\t"
  1084. "movq %%mm3, %%mm5\n\t"
  1085. "movq %%mm6, %%mm0\n\t"
  1086. "movq %%mm7, %%mm1\n\t"
  1087. "movq %%mm4, %%mm6\n\t"
  1088. "movq %%mm5, %%mm7\n\t"
  1089. "movq %%mm0, %%mm2\n\t"
  1090. "movq %%mm1, %%mm3\n\t"
  1091. "psrlq $8, %%mm2\n\t"
  1092. "psrlq $8, %%mm3\n\t"
  1093. "psrlq $8, %%mm6\n\t"
  1094. "psrlq $8, %%mm7\n\t"
  1095. "pand %2, %%mm0\n\t"
  1096. "pand %2, %%mm1\n\t"
  1097. "pand %2, %%mm4\n\t"
  1098. "pand %2, %%mm5\n\t"
  1099. "pand %3, %%mm2\n\t"
  1100. "pand %3, %%mm3\n\t"
  1101. "pand %3, %%mm6\n\t"
  1102. "pand %3, %%mm7\n\t"
  1103. "por %%mm2, %%mm0\n\t"
  1104. "por %%mm3, %%mm1\n\t"
  1105. "por %%mm6, %%mm4\n\t"
  1106. "por %%mm7, %%mm5\n\t"
  1107. "movq %%mm1, %%mm2\n\t"
  1108. "movq %%mm4, %%mm3\n\t"
  1109. "psllq $48, %%mm2\n\t"
  1110. "psllq $32, %%mm3\n\t"
  1111. "pand %4, %%mm2\n\t"
  1112. "pand %5, %%mm3\n\t"
  1113. "por %%mm2, %%mm0\n\t"
  1114. "psrlq $16, %%mm1\n\t"
  1115. "psrlq $32, %%mm4\n\t"
  1116. "psllq $16, %%mm5\n\t"
  1117. "por %%mm3, %%mm1\n\t"
  1118. "pand %6, %%mm5\n\t"
  1119. "por %%mm5, %%mm4\n\t"
  1120. MOVNTQ" %%mm0, %0\n\t"
  1121. MOVNTQ" %%mm1, 8%0\n\t"
  1122. MOVNTQ" %%mm4, 16%0"
  1123. :"=m"(*d)
  1124. :"m"(*s),"m"(mask24l),"m"(mask24h),"m"(mask24hh),"m"(mask24hhh),"m"(mask24hhhh)
  1125. :"memory");
  1126. d += 24;
  1127. s += 8;
  1128. }
  1129. __asm __volatile(SFENCE:::"memory");
  1130. __asm __volatile(EMMS:::"memory");
  1131. #endif
  1132. while(s < end)
  1133. {
  1134. register uint16_t bgr;
  1135. bgr = *s++;
  1136. *d++ = (bgr&0x1F)<<3;
  1137. *d++ = (bgr&0x7E0)>>3;
  1138. *d++ = (bgr&0xF800)>>8;
  1139. }
  1140. }
  1141. static inline void RENAME(rgb15to32)(const uint8_t *src, uint8_t *dst, unsigned src_size)
  1142. {
  1143. const uint16_t *end;
  1144. #ifdef HAVE_MMX
  1145. const uint16_t *mm_end;
  1146. #endif
  1147. uint8_t *d = (uint8_t *)dst;
  1148. const uint16_t *s = (const uint16_t *)src;
  1149. end = s + src_size/2;
  1150. #ifdef HAVE_MMX
  1151. __asm __volatile(PREFETCH" %0"::"m"(*s):"memory");
  1152. __asm __volatile("pxor %%mm7,%%mm7\n\t":::"memory");
  1153. mm_end = end - 3;
  1154. while(s < mm_end)
  1155. {
  1156. __asm __volatile(
  1157. PREFETCH" 32%1\n\t"
  1158. "movq %1, %%mm0\n\t"
  1159. "movq %1, %%mm1\n\t"
  1160. "movq %1, %%mm2\n\t"
  1161. "pand %2, %%mm0\n\t"
  1162. "pand %3, %%mm1\n\t"
  1163. "pand %4, %%mm2\n\t"
  1164. "psllq $3, %%mm0\n\t"
  1165. "psrlq $2, %%mm1\n\t"
  1166. "psrlq $7, %%mm2\n\t"
  1167. "movq %%mm0, %%mm3\n\t"
  1168. "movq %%mm1, %%mm4\n\t"
  1169. "movq %%mm2, %%mm5\n\t"
  1170. "punpcklwd %%mm7, %%mm0\n\t"
  1171. "punpcklwd %%mm7, %%mm1\n\t"
  1172. "punpcklwd %%mm7, %%mm2\n\t"
  1173. "punpckhwd %%mm7, %%mm3\n\t"
  1174. "punpckhwd %%mm7, %%mm4\n\t"
  1175. "punpckhwd %%mm7, %%mm5\n\t"
  1176. "psllq $8, %%mm1\n\t"
  1177. "psllq $16, %%mm2\n\t"
  1178. "por %%mm1, %%mm0\n\t"
  1179. "por %%mm2, %%mm0\n\t"
  1180. "psllq $8, %%mm4\n\t"
  1181. "psllq $16, %%mm5\n\t"
  1182. "por %%mm4, %%mm3\n\t"
  1183. "por %%mm5, %%mm3\n\t"
  1184. MOVNTQ" %%mm0, %0\n\t"
  1185. MOVNTQ" %%mm3, 8%0\n\t"
  1186. :"=m"(*d)
  1187. :"m"(*s),"m"(mask15b),"m"(mask15g),"m"(mask15r)
  1188. :"memory");
  1189. d += 16;
  1190. s += 4;
  1191. }
  1192. __asm __volatile(SFENCE:::"memory");
  1193. __asm __volatile(EMMS:::"memory");
  1194. #endif
  1195. while(s < end)
  1196. {
  1197. #if 0 //slightly slower on athlon
  1198. int bgr= *s++;
  1199. *((uint32_t*)d)++ = ((bgr&0x1F)<<3) + ((bgr&0x3E0)<<6) + ((bgr&0x7C00)<<9);
  1200. #else
  1201. //FIXME this is very likely wrong for bigendian (and the following converters too)
  1202. register uint16_t bgr;
  1203. bgr = *s++;
  1204. *d++ = (bgr&0x1F)<<3;
  1205. *d++ = (bgr&0x3E0)>>2;
  1206. *d++ = (bgr&0x7C00)>>7;
  1207. *d++ = 0;
  1208. #endif
  1209. }
  1210. }
  1211. static inline void RENAME(rgb16to32)(const uint8_t *src, uint8_t *dst, unsigned src_size)
  1212. {
  1213. const uint16_t *end;
  1214. #ifdef HAVE_MMX
  1215. const uint16_t *mm_end;
  1216. #endif
  1217. uint8_t *d = (uint8_t *)dst;
  1218. const uint16_t *s = (uint16_t *)src;
  1219. end = s + src_size/2;
  1220. #ifdef HAVE_MMX
  1221. __asm __volatile(PREFETCH" %0"::"m"(*s):"memory");
  1222. __asm __volatile("pxor %%mm7,%%mm7\n\t":::"memory");
  1223. mm_end = end - 3;
  1224. while(s < mm_end)
  1225. {
  1226. __asm __volatile(
  1227. PREFETCH" 32%1\n\t"
  1228. "movq %1, %%mm0\n\t"
  1229. "movq %1, %%mm1\n\t"
  1230. "movq %1, %%mm2\n\t"
  1231. "pand %2, %%mm0\n\t"
  1232. "pand %3, %%mm1\n\t"
  1233. "pand %4, %%mm2\n\t"
  1234. "psllq $3, %%mm0\n\t"
  1235. "psrlq $3, %%mm1\n\t"
  1236. "psrlq $8, %%mm2\n\t"
  1237. "movq %%mm0, %%mm3\n\t"
  1238. "movq %%mm1, %%mm4\n\t"
  1239. "movq %%mm2, %%mm5\n\t"
  1240. "punpcklwd %%mm7, %%mm0\n\t"
  1241. "punpcklwd %%mm7, %%mm1\n\t"
  1242. "punpcklwd %%mm7, %%mm2\n\t"
  1243. "punpckhwd %%mm7, %%mm3\n\t"
  1244. "punpckhwd %%mm7, %%mm4\n\t"
  1245. "punpckhwd %%mm7, %%mm5\n\t"
  1246. "psllq $8, %%mm1\n\t"
  1247. "psllq $16, %%mm2\n\t"
  1248. "por %%mm1, %%mm0\n\t"
  1249. "por %%mm2, %%mm0\n\t"
  1250. "psllq $8, %%mm4\n\t"
  1251. "psllq $16, %%mm5\n\t"
  1252. "por %%mm4, %%mm3\n\t"
  1253. "por %%mm5, %%mm3\n\t"
  1254. MOVNTQ" %%mm0, %0\n\t"
  1255. MOVNTQ" %%mm3, 8%0\n\t"
  1256. :"=m"(*d)
  1257. :"m"(*s),"m"(mask16b),"m"(mask16g),"m"(mask16r)
  1258. :"memory");
  1259. d += 16;
  1260. s += 4;
  1261. }
  1262. __asm __volatile(SFENCE:::"memory");
  1263. __asm __volatile(EMMS:::"memory");
  1264. #endif
  1265. while(s < end)
  1266. {
  1267. register uint16_t bgr;
  1268. bgr = *s++;
  1269. *d++ = (bgr&0x1F)<<3;
  1270. *d++ = (bgr&0x7E0)>>3;
  1271. *d++ = (bgr&0xF800)>>8;
  1272. *d++ = 0;
  1273. }
  1274. }
  1275. static inline void RENAME(rgb32tobgr32)(const uint8_t *src, uint8_t *dst, unsigned int src_size)
  1276. {
  1277. #ifdef HAVE_MMX
  1278. /* TODO: unroll this loop */
  1279. asm volatile (
  1280. "xorl %%eax, %%eax \n\t"
  1281. ".balign 16 \n\t"
  1282. "1: \n\t"
  1283. PREFETCH" 32(%0, %%eax) \n\t"
  1284. "movq (%0, %%eax), %%mm0 \n\t"
  1285. "movq %%mm0, %%mm1 \n\t"
  1286. "movq %%mm0, %%mm2 \n\t"
  1287. "pslld $16, %%mm0 \n\t"
  1288. "psrld $16, %%mm1 \n\t"
  1289. "pand "MANGLE(mask32r)", %%mm0 \n\t"
  1290. "pand "MANGLE(mask32g)", %%mm2 \n\t"
  1291. "pand "MANGLE(mask32b)", %%mm1 \n\t"
  1292. "por %%mm0, %%mm2 \n\t"
  1293. "por %%mm1, %%mm2 \n\t"
  1294. MOVNTQ" %%mm2, (%1, %%eax) \n\t"
  1295. "addl $8, %%eax \n\t"
  1296. "cmpl %2, %%eax \n\t"
  1297. " jb 1b \n\t"
  1298. :: "r" (src), "r"(dst), "r" (src_size-7)
  1299. : "%eax"
  1300. );
  1301. __asm __volatile(SFENCE:::"memory");
  1302. __asm __volatile(EMMS:::"memory");
  1303. #else
  1304. unsigned i;
  1305. unsigned num_pixels = src_size >> 2;
  1306. for(i=0; i<num_pixels; i++)
  1307. {
  1308. dst[4*i + 0] = src[4*i + 2];
  1309. dst[4*i + 1] = src[4*i + 1];
  1310. dst[4*i + 2] = src[4*i + 0];
  1311. }
  1312. #endif
  1313. }
  1314. static inline void RENAME(rgb24tobgr24)(const uint8_t *src, uint8_t *dst, unsigned int src_size)
  1315. {
  1316. unsigned i;
  1317. #ifdef HAVE_MMX
  1318. int mmx_size= 23 - src_size;
  1319. asm volatile (
  1320. "movq "MANGLE(mask24r)", %%mm5 \n\t"
  1321. "movq "MANGLE(mask24g)", %%mm6 \n\t"
  1322. "movq "MANGLE(mask24b)", %%mm7 \n\t"
  1323. ".balign 16 \n\t"
  1324. "1: \n\t"
  1325. PREFETCH" 32(%1, %%eax) \n\t"
  1326. "movq (%1, %%eax), %%mm0 \n\t" // BGR BGR BG
  1327. "movq (%1, %%eax), %%mm1 \n\t" // BGR BGR BG
  1328. "movq 2(%1, %%eax), %%mm2 \n\t" // R BGR BGR B
  1329. "psllq $16, %%mm0 \n\t" // 00 BGR BGR
  1330. "pand %%mm5, %%mm0 \n\t"
  1331. "pand %%mm6, %%mm1 \n\t"
  1332. "pand %%mm7, %%mm2 \n\t"
  1333. "por %%mm0, %%mm1 \n\t"
  1334. "por %%mm2, %%mm1 \n\t"
  1335. "movq 6(%1, %%eax), %%mm0 \n\t" // BGR BGR BG
  1336. MOVNTQ" %%mm1, (%2, %%eax) \n\t" // RGB RGB RG
  1337. "movq 8(%1, %%eax), %%mm1 \n\t" // R BGR BGR B
  1338. "movq 10(%1, %%eax), %%mm2 \n\t" // GR BGR BGR
  1339. "pand %%mm7, %%mm0 \n\t"
  1340. "pand %%mm5, %%mm1 \n\t"
  1341. "pand %%mm6, %%mm2 \n\t"
  1342. "por %%mm0, %%mm1 \n\t"
  1343. "por %%mm2, %%mm1 \n\t"
  1344. "movq 14(%1, %%eax), %%mm0 \n\t" // R BGR BGR B
  1345. MOVNTQ" %%mm1, 8(%2, %%eax) \n\t" // B RGB RGB R
  1346. "movq 16(%1, %%eax), %%mm1 \n\t" // GR BGR BGR
  1347. "movq 18(%1, %%eax), %%mm2 \n\t" // BGR BGR BG
  1348. "pand %%mm6, %%mm0 \n\t"
  1349. "pand %%mm7, %%mm1 \n\t"
  1350. "pand %%mm5, %%mm2 \n\t"
  1351. "por %%mm0, %%mm1 \n\t"
  1352. "por %%mm2, %%mm1 \n\t"
  1353. MOVNTQ" %%mm1, 16(%2, %%eax) \n\t"
  1354. "addl $24, %%eax \n\t"
  1355. " js 1b \n\t"
  1356. : "+a" (mmx_size)
  1357. : "r" (src-mmx_size), "r"(dst-mmx_size)
  1358. );
  1359. __asm __volatile(SFENCE:::"memory");
  1360. __asm __volatile(EMMS:::"memory");
  1361. if(mmx_size==23) return; //finihsed, was multiple of 8
  1362. src+= src_size;
  1363. dst+= src_size;
  1364. src_size= 23-mmx_size;
  1365. src-= src_size;
  1366. dst-= src_size;
  1367. #endif
  1368. for(i=0; i<src_size; i+=3)
  1369. {
  1370. register uint8_t x;
  1371. x = src[i + 2];
  1372. dst[i + 1] = src[i + 1];
  1373. dst[i + 2] = src[i + 0];
  1374. dst[i + 0] = x;
  1375. }
  1376. }
  1377. static inline void RENAME(yuvPlanartoyuy2)(const uint8_t *ysrc, const uint8_t *usrc, const uint8_t *vsrc, uint8_t *dst,
  1378. unsigned int width, unsigned int height,
  1379. int lumStride, int chromStride, int dstStride, int vertLumPerChroma)
  1380. {
  1381. unsigned y;
  1382. const unsigned chromWidth= width>>1;
  1383. for(y=0; y<height; y++)
  1384. {
  1385. #ifdef HAVE_MMX
  1386. //FIXME handle 2 lines a once (fewer prefetch, reuse some chrom, but very likely limited by mem anyway)
  1387. asm volatile(
  1388. "xorl %%eax, %%eax \n\t"
  1389. ".balign 16 \n\t"
  1390. "1: \n\t"
  1391. PREFETCH" 32(%1, %%eax, 2) \n\t"
  1392. PREFETCH" 32(%2, %%eax) \n\t"
  1393. PREFETCH" 32(%3, %%eax) \n\t"
  1394. "movq (%2, %%eax), %%mm0 \n\t" // U(0)
  1395. "movq %%mm0, %%mm2 \n\t" // U(0)
  1396. "movq (%3, %%eax), %%mm1 \n\t" // V(0)
  1397. "punpcklbw %%mm1, %%mm0 \n\t" // UVUV UVUV(0)
  1398. "punpckhbw %%mm1, %%mm2 \n\t" // UVUV UVUV(8)
  1399. "movq (%1, %%eax,2), %%mm3 \n\t" // Y(0)
  1400. "movq 8(%1, %%eax,2), %%mm5 \n\t" // Y(8)
  1401. "movq %%mm3, %%mm4 \n\t" // Y(0)
  1402. "movq %%mm5, %%mm6 \n\t" // Y(8)
  1403. "punpcklbw %%mm0, %%mm3 \n\t" // YUYV YUYV(0)
  1404. "punpckhbw %%mm0, %%mm4 \n\t" // YUYV YUYV(4)
  1405. "punpcklbw %%mm2, %%mm5 \n\t" // YUYV YUYV(8)
  1406. "punpckhbw %%mm2, %%mm6 \n\t" // YUYV YUYV(12)
  1407. MOVNTQ" %%mm3, (%0, %%eax, 4) \n\t"
  1408. MOVNTQ" %%mm4, 8(%0, %%eax, 4) \n\t"
  1409. MOVNTQ" %%mm5, 16(%0, %%eax, 4) \n\t"
  1410. MOVNTQ" %%mm6, 24(%0, %%eax, 4) \n\t"
  1411. "addl $8, %%eax \n\t"
  1412. "cmpl %4, %%eax \n\t"
  1413. " jb 1b \n\t"
  1414. ::"r"(dst), "r"(ysrc), "r"(usrc), "r"(vsrc), "g" (chromWidth)
  1415. : "%eax"
  1416. );
  1417. #else
  1418. #if defined ARCH_ALPHA && defined HAVE_MVI
  1419. #define pl2yuy2(n) \
  1420. y1 = yc[n]; \
  1421. y2 = yc2[n]; \
  1422. u = uc[n]; \
  1423. v = vc[n]; \
  1424. asm("unpkbw %1, %0" : "=r"(y1) : "r"(y1)); \
  1425. asm("unpkbw %1, %0" : "=r"(y2) : "r"(y2)); \
  1426. asm("unpkbl %1, %0" : "=r"(u) : "r"(u)); \
  1427. asm("unpkbl %1, %0" : "=r"(v) : "r"(v)); \
  1428. yuv1 = (u << 8) + (v << 24); \
  1429. yuv2 = yuv1 + y2; \
  1430. yuv1 += y1; \
  1431. qdst[n] = yuv1; \
  1432. qdst2[n] = yuv2;
  1433. int i;
  1434. uint64_t *qdst = (uint64_t *) dst;
  1435. uint64_t *qdst2 = (uint64_t *) (dst + dstStride);
  1436. const uint32_t *yc = (uint32_t *) ysrc;
  1437. const uint32_t *yc2 = (uint32_t *) (ysrc + lumStride);
  1438. const uint16_t *uc = (uint16_t*) usrc, *vc = (uint16_t*) vsrc;
  1439. for(i = 0; i < chromWidth; i += 8){
  1440. uint64_t y1, y2, yuv1, yuv2;
  1441. uint64_t u, v;
  1442. /* Prefetch */
  1443. asm("ldq $31,64(%0)" :: "r"(yc));
  1444. asm("ldq $31,64(%0)" :: "r"(yc2));
  1445. asm("ldq $31,64(%0)" :: "r"(uc));
  1446. asm("ldq $31,64(%0)" :: "r"(vc));
  1447. pl2yuy2(0);
  1448. pl2yuy2(1);
  1449. pl2yuy2(2);
  1450. pl2yuy2(3);
  1451. yc += 4;
  1452. yc2 += 4;
  1453. uc += 4;
  1454. vc += 4;
  1455. qdst += 4;
  1456. qdst2 += 4;
  1457. }
  1458. y++;
  1459. ysrc += lumStride;
  1460. dst += dstStride;
  1461. #elif __WORDSIZE >= 64
  1462. int i;
  1463. uint64_t *ldst = (uint64_t *) dst;
  1464. const uint8_t *yc = ysrc, *uc = usrc, *vc = vsrc;
  1465. for(i = 0; i < chromWidth; i += 2){
  1466. uint64_t k, l;
  1467. k = yc[0] + (uc[0] << 8) +
  1468. (yc[1] << 16) + (vc[0] << 24);
  1469. l = yc[2] + (uc[1] << 8) +
  1470. (yc[3] << 16) + (vc[1] << 24);
  1471. *ldst++ = k + (l << 32);
  1472. yc += 4;
  1473. uc += 2;
  1474. vc += 2;
  1475. }
  1476. #else
  1477. int i, *idst = (int32_t *) dst;
  1478. const uint8_t *yc = ysrc, *uc = usrc, *vc = vsrc;
  1479. for(i = 0; i < chromWidth; i++){
  1480. *idst++ = yc[0] + (uc[0] << 8) +
  1481. (yc[1] << 16) + (vc[0] << 24);
  1482. yc += 2;
  1483. uc++;
  1484. vc++;
  1485. }
  1486. #endif
  1487. #endif
  1488. if((y&(vertLumPerChroma-1))==(vertLumPerChroma-1) )
  1489. {
  1490. usrc += chromStride;
  1491. vsrc += chromStride;
  1492. }
  1493. ysrc += lumStride;
  1494. dst += dstStride;
  1495. }
  1496. #ifdef HAVE_MMX
  1497. asm( EMMS" \n\t"
  1498. SFENCE" \n\t"
  1499. :::"memory");
  1500. #endif
  1501. }
  1502. /**
  1503. *
  1504. * height should be a multiple of 2 and width should be a multiple of 16 (if this is a
  1505. * problem for anyone then tell me, and ill fix it)
  1506. */
  1507. static inline void RENAME(yv12toyuy2)(const uint8_t *ysrc, const uint8_t *usrc, const uint8_t *vsrc, uint8_t *dst,
  1508. unsigned int width, unsigned int height,
  1509. int lumStride, int chromStride, int dstStride)
  1510. {
  1511. //FIXME interpolate chroma
  1512. RENAME(yuvPlanartoyuy2)(ysrc, usrc, vsrc, dst, width, height, lumStride, chromStride, dstStride, 2);
  1513. }
  1514. /**
  1515. *
  1516. * width should be a multiple of 16
  1517. */
  1518. static inline void RENAME(yuv422ptoyuy2)(const uint8_t *ysrc, const uint8_t *usrc, const uint8_t *vsrc, uint8_t *dst,
  1519. unsigned int width, unsigned int height,
  1520. int lumStride, int chromStride, int dstStride)
  1521. {
  1522. RENAME(yuvPlanartoyuy2)(ysrc, usrc, vsrc, dst, width, height, lumStride, chromStride, dstStride, 1);
  1523. }
  1524. /**
  1525. *
  1526. * height should be a multiple of 2 and width should be a multiple of 16 (if this is a
  1527. * problem for anyone then tell me, and ill fix it)
  1528. */
  1529. static inline void RENAME(yuy2toyv12)(const uint8_t *src, uint8_t *ydst, uint8_t *udst, uint8_t *vdst,
  1530. unsigned int width, unsigned int height,
  1531. int lumStride, int chromStride, int srcStride)
  1532. {
  1533. unsigned y;
  1534. const unsigned chromWidth= width>>1;
  1535. for(y=0; y<height; y+=2)
  1536. {
  1537. #ifdef HAVE_MMX
  1538. asm volatile(
  1539. "xorl %%eax, %%eax \n\t"
  1540. "pcmpeqw %%mm7, %%mm7 \n\t"
  1541. "psrlw $8, %%mm7 \n\t" // FF,00,FF,00...
  1542. ".balign 16 \n\t"
  1543. "1: \n\t"
  1544. PREFETCH" 64(%0, %%eax, 4) \n\t"
  1545. "movq (%0, %%eax, 4), %%mm0 \n\t" // YUYV YUYV(0)
  1546. "movq 8(%0, %%eax, 4), %%mm1 \n\t" // YUYV YUYV(4)
  1547. "movq %%mm0, %%mm2 \n\t" // YUYV YUYV(0)
  1548. "movq %%mm1, %%mm3 \n\t" // YUYV YUYV(4)
  1549. "psrlw $8, %%mm0 \n\t" // U0V0 U0V0(0)
  1550. "psrlw $8, %%mm1 \n\t" // U0V0 U0V0(4)
  1551. "pand %%mm7, %%mm2 \n\t" // Y0Y0 Y0Y0(0)
  1552. "pand %%mm7, %%mm3 \n\t" // Y0Y0 Y0Y0(4)
  1553. "packuswb %%mm1, %%mm0 \n\t" // UVUV UVUV(0)
  1554. "packuswb %%mm3, %%mm2 \n\t" // YYYY YYYY(0)
  1555. MOVNTQ" %%mm2, (%1, %%eax, 2) \n\t"
  1556. "movq 16(%0, %%eax, 4), %%mm1 \n\t" // YUYV YUYV(8)
  1557. "movq 24(%0, %%eax, 4), %%mm2 \n\t" // YUYV YUYV(12)
  1558. "movq %%mm1, %%mm3 \n\t" // YUYV YUYV(8)
  1559. "movq %%mm2, %%mm4 \n\t" // YUYV YUYV(12)
  1560. "psrlw $8, %%mm1 \n\t" // U0V0 U0V0(8)
  1561. "psrlw $8, %%mm2 \n\t" // U0V0 U0V0(12)
  1562. "pand %%mm7, %%mm3 \n\t" // Y0Y0 Y0Y0(8)
  1563. "pand %%mm7, %%mm4 \n\t" // Y0Y0 Y0Y0(12)
  1564. "packuswb %%mm2, %%mm1 \n\t" // UVUV UVUV(8)
  1565. "packuswb %%mm4, %%mm3 \n\t" // YYYY YYYY(8)
  1566. MOVNTQ" %%mm3, 8(%1, %%eax, 2) \n\t"
  1567. "movq %%mm0, %%mm2 \n\t" // UVUV UVUV(0)
  1568. "movq %%mm1, %%mm3 \n\t" // UVUV UVUV(8)
  1569. "psrlw $8, %%mm0 \n\t" // V0V0 V0V0(0)
  1570. "psrlw $8, %%mm1 \n\t" // V0V0 V0V0(8)
  1571. "pand %%mm7, %%mm2 \n\t" // U0U0 U0U0(0)
  1572. "pand %%mm7, %%mm3 \n\t" // U0U0 U0U0(8)
  1573. "packuswb %%mm1, %%mm0 \n\t" // VVVV VVVV(0)
  1574. "packuswb %%mm3, %%mm2 \n\t" // UUUU UUUU(0)
  1575. MOVNTQ" %%mm0, (%3, %%eax) \n\t"
  1576. MOVNTQ" %%mm2, (%2, %%eax) \n\t"
  1577. "addl $8, %%eax \n\t"
  1578. "cmpl %4, %%eax \n\t"
  1579. " jb 1b \n\t"
  1580. ::"r"(src), "r"(ydst), "r"(udst), "r"(vdst), "g" (chromWidth)
  1581. : "memory", "%eax"
  1582. );
  1583. ydst += lumStride;
  1584. src += srcStride;
  1585. asm volatile(
  1586. "xorl %%eax, %%eax \n\t"
  1587. ".balign 16 \n\t"
  1588. "1: \n\t"
  1589. PREFETCH" 64(%0, %%eax, 4) \n\t"
  1590. "movq (%0, %%eax, 4), %%mm0 \n\t" // YUYV YUYV(0)
  1591. "movq 8(%0, %%eax, 4), %%mm1 \n\t" // YUYV YUYV(4)
  1592. "movq 16(%0, %%eax, 4), %%mm2 \n\t" // YUYV YUYV(8)
  1593. "movq 24(%0, %%eax, 4), %%mm3 \n\t" // YUYV YUYV(12)
  1594. "pand %%mm7, %%mm0 \n\t" // Y0Y0 Y0Y0(0)
  1595. "pand %%mm7, %%mm1 \n\t" // Y0Y0 Y0Y0(4)
  1596. "pand %%mm7, %%mm2 \n\t" // Y0Y0 Y0Y0(8)
  1597. "pand %%mm7, %%mm3 \n\t" // Y0Y0 Y0Y0(12)
  1598. "packuswb %%mm1, %%mm0 \n\t" // YYYY YYYY(0)
  1599. "packuswb %%mm3, %%mm2 \n\t" // YYYY YYYY(8)
  1600. MOVNTQ" %%mm0, (%1, %%eax, 2) \n\t"
  1601. MOVNTQ" %%mm2, 8(%1, %%eax, 2) \n\t"
  1602. "addl $8, %%eax \n\t"
  1603. "cmpl %4, %%eax \n\t"
  1604. " jb 1b \n\t"
  1605. ::"r"(src), "r"(ydst), "r"(udst), "r"(vdst), "g" (chromWidth)
  1606. : "memory", "%eax"
  1607. );
  1608. #else
  1609. unsigned i;
  1610. for(i=0; i<chromWidth; i++)
  1611. {
  1612. ydst[2*i+0] = src[4*i+0];
  1613. udst[i] = src[4*i+1];
  1614. ydst[2*i+1] = src[4*i+2];
  1615. vdst[i] = src[4*i+3];
  1616. }
  1617. ydst += lumStride;
  1618. src += srcStride;
  1619. for(i=0; i<chromWidth; i++)
  1620. {
  1621. ydst[2*i+0] = src[4*i+0];
  1622. ydst[2*i+1] = src[4*i+2];
  1623. }
  1624. #endif
  1625. udst += chromStride;
  1626. vdst += chromStride;
  1627. ydst += lumStride;
  1628. src += srcStride;
  1629. }
  1630. #ifdef HAVE_MMX
  1631. asm volatile( EMMS" \n\t"
  1632. SFENCE" \n\t"
  1633. :::"memory");
  1634. #endif
  1635. }
  1636. static inline void RENAME(yvu9toyv12)(const uint8_t *ysrc, const uint8_t *usrc, const uint8_t *vsrc,
  1637. uint8_t *ydst, uint8_t *udst, uint8_t *vdst,
  1638. unsigned int width, unsigned int height, int lumStride, int chromStride)
  1639. {
  1640. /* Y Plane */
  1641. memcpy(ydst, ysrc, width*height);
  1642. /* XXX: implement upscaling for U,V */
  1643. }
  1644. static inline void RENAME(planar2x)(const uint8_t *src, uint8_t *dst, int srcWidth, int srcHeight, int srcStride, int dstStride)
  1645. {
  1646. int x,y;
  1647. dst[0]= src[0];
  1648. // first line
  1649. for(x=0; x<srcWidth-1; x++){
  1650. dst[2*x+1]= (3*src[x] + src[x+1])>>2;
  1651. dst[2*x+2]= ( src[x] + 3*src[x+1])>>2;
  1652. }
  1653. dst[2*srcWidth-1]= src[srcWidth-1];
  1654. dst+= dstStride;
  1655. for(y=1; y<srcHeight; y++){
  1656. #if defined (HAVE_MMX2) || defined (HAVE_3DNOW)
  1657. const int mmxSize= srcWidth&~15;
  1658. asm volatile(
  1659. "movl %4, %%eax \n\t"
  1660. "1: \n\t"
  1661. "movq (%0, %%eax), %%mm0 \n\t"
  1662. "movq (%1, %%eax), %%mm1 \n\t"
  1663. "movq 1(%0, %%eax), %%mm2 \n\t"
  1664. "movq 1(%1, %%eax), %%mm3 \n\t"
  1665. "movq -1(%0, %%eax), %%mm4 \n\t"
  1666. "movq -1(%1, %%eax), %%mm5 \n\t"
  1667. PAVGB" %%mm0, %%mm5 \n\t"
  1668. PAVGB" %%mm0, %%mm3 \n\t"
  1669. PAVGB" %%mm0, %%mm5 \n\t"
  1670. PAVGB" %%mm0, %%mm3 \n\t"
  1671. PAVGB" %%mm1, %%mm4 \n\t"
  1672. PAVGB" %%mm1, %%mm2 \n\t"
  1673. PAVGB" %%mm1, %%mm4 \n\t"
  1674. PAVGB" %%mm1, %%mm2 \n\t"
  1675. "movq %%mm5, %%mm7 \n\t"
  1676. "movq %%mm4, %%mm6 \n\t"
  1677. "punpcklbw %%mm3, %%mm5 \n\t"
  1678. "punpckhbw %%mm3, %%mm7 \n\t"
  1679. "punpcklbw %%mm2, %%mm4 \n\t"
  1680. "punpckhbw %%mm2, %%mm6 \n\t"
  1681. #if 1
  1682. MOVNTQ" %%mm5, (%2, %%eax, 2) \n\t"
  1683. MOVNTQ" %%mm7, 8(%2, %%eax, 2) \n\t"
  1684. MOVNTQ" %%mm4, (%3, %%eax, 2) \n\t"
  1685. MOVNTQ" %%mm6, 8(%3, %%eax, 2) \n\t"
  1686. #else
  1687. "movq %%mm5, (%2, %%eax, 2) \n\t"
  1688. "movq %%mm7, 8(%2, %%eax, 2) \n\t"
  1689. "movq %%mm4, (%3, %%eax, 2) \n\t"
  1690. "movq %%mm6, 8(%3, %%eax, 2) \n\t"
  1691. #endif
  1692. "addl $8, %%eax \n\t"
  1693. " js 1b \n\t"
  1694. :: "r" (src + mmxSize ), "r" (src + srcStride + mmxSize ),
  1695. "r" (dst + mmxSize*2), "r" (dst + dstStride + mmxSize*2),
  1696. "g" (-mmxSize)
  1697. : "%eax"
  1698. );
  1699. #else
  1700. const int mmxSize=1;
  1701. #endif
  1702. dst[0 ]= (3*src[0] + src[srcStride])>>2;
  1703. dst[dstStride]= ( src[0] + 3*src[srcStride])>>2;
  1704. for(x=mmxSize-1; x<srcWidth-1; x++){
  1705. dst[2*x +1]= (3*src[x+0] + src[x+srcStride+1])>>2;
  1706. dst[2*x+dstStride+2]= ( src[x+0] + 3*src[x+srcStride+1])>>2;
  1707. dst[2*x+dstStride+1]= ( src[x+1] + 3*src[x+srcStride ])>>2;
  1708. dst[2*x +2]= (3*src[x+1] + src[x+srcStride ])>>2;
  1709. }
  1710. dst[srcWidth*2 -1 ]= (3*src[srcWidth-1] + src[srcWidth-1 + srcStride])>>2;
  1711. dst[srcWidth*2 -1 + dstStride]= ( src[srcWidth-1] + 3*src[srcWidth-1 + srcStride])>>2;
  1712. dst+=dstStride*2;
  1713. src+=srcStride;
  1714. }
  1715. // last line
  1716. #if 1
  1717. dst[0]= src[0];
  1718. for(x=0; x<srcWidth-1; x++){
  1719. dst[2*x+1]= (3*src[x] + src[x+1])>>2;
  1720. dst[2*x+2]= ( src[x] + 3*src[x+1])>>2;
  1721. }
  1722. dst[2*srcWidth-1]= src[srcWidth-1];
  1723. #else
  1724. for(x=0; x<srcWidth; x++){
  1725. dst[2*x+0]=
  1726. dst[2*x+1]= src[x];
  1727. }
  1728. #endif
  1729. #ifdef HAVE_MMX
  1730. asm volatile( EMMS" \n\t"
  1731. SFENCE" \n\t"
  1732. :::"memory");
  1733. #endif
  1734. }
  1735. /**
  1736. *
  1737. * height should be a multiple of 2 and width should be a multiple of 16 (if this is a
  1738. * problem for anyone then tell me, and ill fix it)
  1739. * chrominance data is only taken from every secound line others are ignored FIXME write HQ version
  1740. */
  1741. static inline void RENAME(uyvytoyv12)(const uint8_t *src, uint8_t *ydst, uint8_t *udst, uint8_t *vdst,
  1742. unsigned int width, unsigned int height,
  1743. int lumStride, int chromStride, int srcStride)
  1744. {
  1745. unsigned y;
  1746. const unsigned chromWidth= width>>1;
  1747. for(y=0; y<height; y+=2)
  1748. {
  1749. #ifdef HAVE_MMX
  1750. asm volatile(
  1751. "xorl %%eax, %%eax \n\t"
  1752. "pcmpeqw %%mm7, %%mm7 \n\t"
  1753. "psrlw $8, %%mm7 \n\t" // FF,00,FF,00...
  1754. ".balign 16 \n\t"
  1755. "1: \n\t"
  1756. PREFETCH" 64(%0, %%eax, 4) \n\t"
  1757. "movq (%0, %%eax, 4), %%mm0 \n\t" // UYVY UYVY(0)
  1758. "movq 8(%0, %%eax, 4), %%mm1 \n\t" // UYVY UYVY(4)
  1759. "movq %%mm0, %%mm2 \n\t" // UYVY UYVY(0)
  1760. "movq %%mm1, %%mm3 \n\t" // UYVY UYVY(4)
  1761. "pand %%mm7, %%mm0 \n\t" // U0V0 U0V0(0)
  1762. "pand %%mm7, %%mm1 \n\t" // U0V0 U0V0(4)
  1763. "psrlw $8, %%mm2 \n\t" // Y0Y0 Y0Y0(0)
  1764. "psrlw $8, %%mm3 \n\t" // Y0Y0 Y0Y0(4)
  1765. "packuswb %%mm1, %%mm0 \n\t" // UVUV UVUV(0)
  1766. "packuswb %%mm3, %%mm2 \n\t" // YYYY YYYY(0)
  1767. MOVNTQ" %%mm2, (%1, %%eax, 2) \n\t"
  1768. "movq 16(%0, %%eax, 4), %%mm1 \n\t" // UYVY UYVY(8)
  1769. "movq 24(%0, %%eax, 4), %%mm2 \n\t" // UYVY UYVY(12)
  1770. "movq %%mm1, %%mm3 \n\t" // UYVY UYVY(8)
  1771. "movq %%mm2, %%mm4 \n\t" // UYVY UYVY(12)
  1772. "pand %%mm7, %%mm1 \n\t" // U0V0 U0V0(8)
  1773. "pand %%mm7, %%mm2 \n\t" // U0V0 U0V0(12)
  1774. "psrlw $8, %%mm3 \n\t" // Y0Y0 Y0Y0(8)
  1775. "psrlw $8, %%mm4 \n\t" // Y0Y0 Y0Y0(12)
  1776. "packuswb %%mm2, %%mm1 \n\t" // UVUV UVUV(8)
  1777. "packuswb %%mm4, %%mm3 \n\t" // YYYY YYYY(8)
  1778. MOVNTQ" %%mm3, 8(%1, %%eax, 2) \n\t"
  1779. "movq %%mm0, %%mm2 \n\t" // UVUV UVUV(0)
  1780. "movq %%mm1, %%mm3 \n\t" // UVUV UVUV(8)
  1781. "psrlw $8, %%mm0 \n\t" // V0V0 V0V0(0)
  1782. "psrlw $8, %%mm1 \n\t" // V0V0 V0V0(8)
  1783. "pand %%mm7, %%mm2 \n\t" // U0U0 U0U0(0)
  1784. "pand %%mm7, %%mm3 \n\t" // U0U0 U0U0(8)
  1785. "packuswb %%mm1, %%mm0 \n\t" // VVVV VVVV(0)
  1786. "packuswb %%mm3, %%mm2 \n\t" // UUUU UUUU(0)
  1787. MOVNTQ" %%mm0, (%3, %%eax) \n\t"
  1788. MOVNTQ" %%mm2, (%2, %%eax) \n\t"
  1789. "addl $8, %%eax \n\t"
  1790. "cmpl %4, %%eax \n\t"
  1791. " jb 1b \n\t"
  1792. ::"r"(src), "r"(ydst), "r"(udst), "r"(vdst), "g" (chromWidth)
  1793. : "memory", "%eax"
  1794. );
  1795. ydst += lumStride;
  1796. src += srcStride;
  1797. asm volatile(
  1798. "xorl %%eax, %%eax \n\t"
  1799. ".balign 16 \n\t"
  1800. "1: \n\t"
  1801. PREFETCH" 64(%0, %%eax, 4) \n\t"
  1802. "movq (%0, %%eax, 4), %%mm0 \n\t" // YUYV YUYV(0)
  1803. "movq 8(%0, %%eax, 4), %%mm1 \n\t" // YUYV YUYV(4)
  1804. "movq 16(%0, %%eax, 4), %%mm2 \n\t" // YUYV YUYV(8)
  1805. "movq 24(%0, %%eax, 4), %%mm3 \n\t" // YUYV YUYV(12)
  1806. "psrlw $8, %%mm0 \n\t" // Y0Y0 Y0Y0(0)
  1807. "psrlw $8, %%mm1 \n\t" // Y0Y0 Y0Y0(4)
  1808. "psrlw $8, %%mm2 \n\t" // Y0Y0 Y0Y0(8)
  1809. "psrlw $8, %%mm3 \n\t" // Y0Y0 Y0Y0(12)
  1810. "packuswb %%mm1, %%mm0 \n\t" // YYYY YYYY(0)
  1811. "packuswb %%mm3, %%mm2 \n\t" // YYYY YYYY(8)
  1812. MOVNTQ" %%mm0, (%1, %%eax, 2) \n\t"
  1813. MOVNTQ" %%mm2, 8(%1, %%eax, 2) \n\t"
  1814. "addl $8, %%eax \n\t"
  1815. "cmpl %4, %%eax \n\t"
  1816. " jb 1b \n\t"
  1817. ::"r"(src), "r"(ydst), "r"(udst), "r"(vdst), "g" (chromWidth)
  1818. : "memory", "%eax"
  1819. );
  1820. #else
  1821. unsigned i;
  1822. for(i=0; i<chromWidth; i++)
  1823. {
  1824. udst[i] = src[4*i+0];
  1825. ydst[2*i+0] = src[4*i+1];
  1826. vdst[i] = src[4*i+2];
  1827. ydst[2*i+1] = src[4*i+3];
  1828. }
  1829. ydst += lumStride;
  1830. src += srcStride;
  1831. for(i=0; i<chromWidth; i++)
  1832. {
  1833. ydst[2*i+0] = src[4*i+1];
  1834. ydst[2*i+1] = src[4*i+3];
  1835. }
  1836. #endif
  1837. udst += chromStride;
  1838. vdst += chromStride;
  1839. ydst += lumStride;
  1840. src += srcStride;
  1841. }
  1842. #ifdef HAVE_MMX
  1843. asm volatile( EMMS" \n\t"
  1844. SFENCE" \n\t"
  1845. :::"memory");
  1846. #endif
  1847. }
  1848. /**
  1849. *
  1850. * height should be a multiple of 2 and width should be a multiple of 2 (if this is a
  1851. * problem for anyone then tell me, and ill fix it)
  1852. * chrominance data is only taken from every secound line others are ignored in the C version FIXME write HQ version
  1853. */
  1854. static inline void RENAME(rgb24toyv12)(const uint8_t *src, uint8_t *ydst, uint8_t *udst, uint8_t *vdst,
  1855. unsigned int width, unsigned int height,
  1856. int lumStride, int chromStride, int srcStride)
  1857. {
  1858. unsigned y;
  1859. const unsigned chromWidth= width>>1;
  1860. #ifdef HAVE_MMX
  1861. for(y=0; y<height-2; y+=2)
  1862. {
  1863. unsigned i;
  1864. for(i=0; i<2; i++)
  1865. {
  1866. asm volatile(
  1867. "movl %2, %%eax \n\t"
  1868. "movq "MANGLE(bgr2YCoeff)", %%mm6 \n\t"
  1869. "movq "MANGLE(w1111)", %%mm5 \n\t"
  1870. "pxor %%mm7, %%mm7 \n\t"
  1871. "leal (%%eax, %%eax, 2), %%ebx \n\t"
  1872. ".balign 16 \n\t"
  1873. "1: \n\t"
  1874. PREFETCH" 64(%0, %%ebx) \n\t"
  1875. "movd (%0, %%ebx), %%mm0 \n\t"
  1876. "movd 3(%0, %%ebx), %%mm1 \n\t"
  1877. "punpcklbw %%mm7, %%mm0 \n\t"
  1878. "punpcklbw %%mm7, %%mm1 \n\t"
  1879. "movd 6(%0, %%ebx), %%mm2 \n\t"
  1880. "movd 9(%0, %%ebx), %%mm3 \n\t"
  1881. "punpcklbw %%mm7, %%mm2 \n\t"
  1882. "punpcklbw %%mm7, %%mm3 \n\t"
  1883. "pmaddwd %%mm6, %%mm0 \n\t"
  1884. "pmaddwd %%mm6, %%mm1 \n\t"
  1885. "pmaddwd %%mm6, %%mm2 \n\t"
  1886. "pmaddwd %%mm6, %%mm3 \n\t"
  1887. #ifndef FAST_BGR2YV12
  1888. "psrad $8, %%mm0 \n\t"
  1889. "psrad $8, %%mm1 \n\t"
  1890. "psrad $8, %%mm2 \n\t"
  1891. "psrad $8, %%mm3 \n\t"
  1892. #endif
  1893. "packssdw %%mm1, %%mm0 \n\t"
  1894. "packssdw %%mm3, %%mm2 \n\t"
  1895. "pmaddwd %%mm5, %%mm0 \n\t"
  1896. "pmaddwd %%mm5, %%mm2 \n\t"
  1897. "packssdw %%mm2, %%mm0 \n\t"
  1898. "psraw $7, %%mm0 \n\t"
  1899. "movd 12(%0, %%ebx), %%mm4 \n\t"
  1900. "movd 15(%0, %%ebx), %%mm1 \n\t"
  1901. "punpcklbw %%mm7, %%mm4 \n\t"
  1902. "punpcklbw %%mm7, %%mm1 \n\t"
  1903. "movd 18(%0, %%ebx), %%mm2 \n\t"
  1904. "movd 21(%0, %%ebx), %%mm3 \n\t"
  1905. "punpcklbw %%mm7, %%mm2 \n\t"
  1906. "punpcklbw %%mm7, %%mm3 \n\t"
  1907. "pmaddwd %%mm6, %%mm4 \n\t"
  1908. "pmaddwd %%mm6, %%mm1 \n\t"
  1909. "pmaddwd %%mm6, %%mm2 \n\t"
  1910. "pmaddwd %%mm6, %%mm3 \n\t"
  1911. #ifndef FAST_BGR2YV12
  1912. "psrad $8, %%mm4 \n\t"
  1913. "psrad $8, %%mm1 \n\t"
  1914. "psrad $8, %%mm2 \n\t"
  1915. "psrad $8, %%mm3 \n\t"
  1916. #endif
  1917. "packssdw %%mm1, %%mm4 \n\t"
  1918. "packssdw %%mm3, %%mm2 \n\t"
  1919. "pmaddwd %%mm5, %%mm4 \n\t"
  1920. "pmaddwd %%mm5, %%mm2 \n\t"
  1921. "addl $24, %%ebx \n\t"
  1922. "packssdw %%mm2, %%mm4 \n\t"
  1923. "psraw $7, %%mm4 \n\t"
  1924. "packuswb %%mm4, %%mm0 \n\t"
  1925. "paddusb "MANGLE(bgr2YOffset)", %%mm0 \n\t"
  1926. MOVNTQ" %%mm0, (%1, %%eax) \n\t"
  1927. "addl $8, %%eax \n\t"
  1928. " js 1b \n\t"
  1929. : : "r" (src+width*3), "r" (ydst+width), "g" (-width)
  1930. : "%eax", "%ebx"
  1931. );
  1932. ydst += lumStride;
  1933. src += srcStride;
  1934. }
  1935. src -= srcStride*2;
  1936. asm volatile(
  1937. "movl %4, %%eax \n\t"
  1938. "movq "MANGLE(w1111)", %%mm5 \n\t"
  1939. "movq "MANGLE(bgr2UCoeff)", %%mm6 \n\t"
  1940. "pxor %%mm7, %%mm7 \n\t"
  1941. "leal (%%eax, %%eax, 2), %%ebx \n\t"
  1942. "addl %%ebx, %%ebx \n\t"
  1943. ".balign 16 \n\t"
  1944. "1: \n\t"
  1945. PREFETCH" 64(%0, %%ebx) \n\t"
  1946. PREFETCH" 64(%1, %%ebx) \n\t"
  1947. #if defined (HAVE_MMX2) || defined (HAVE_3DNOW)
  1948. "movq (%0, %%ebx), %%mm0 \n\t"
  1949. "movq (%1, %%ebx), %%mm1 \n\t"
  1950. "movq 6(%0, %%ebx), %%mm2 \n\t"
  1951. "movq 6(%1, %%ebx), %%mm3 \n\t"
  1952. PAVGB" %%mm1, %%mm0 \n\t"
  1953. PAVGB" %%mm3, %%mm2 \n\t"
  1954. "movq %%mm0, %%mm1 \n\t"
  1955. "movq %%mm2, %%mm3 \n\t"
  1956. "psrlq $24, %%mm0 \n\t"
  1957. "psrlq $24, %%mm2 \n\t"
  1958. PAVGB" %%mm1, %%mm0 \n\t"
  1959. PAVGB" %%mm3, %%mm2 \n\t"
  1960. "punpcklbw %%mm7, %%mm0 \n\t"
  1961. "punpcklbw %%mm7, %%mm2 \n\t"
  1962. #else
  1963. "movd (%0, %%ebx), %%mm0 \n\t"
  1964. "movd (%1, %%ebx), %%mm1 \n\t"
  1965. "movd 3(%0, %%ebx), %%mm2 \n\t"
  1966. "movd 3(%1, %%ebx), %%mm3 \n\t"
  1967. "punpcklbw %%mm7, %%mm0 \n\t"
  1968. "punpcklbw %%mm7, %%mm1 \n\t"
  1969. "punpcklbw %%mm7, %%mm2 \n\t"
  1970. "punpcklbw %%mm7, %%mm3 \n\t"
  1971. "paddw %%mm1, %%mm0 \n\t"
  1972. "paddw %%mm3, %%mm2 \n\t"
  1973. "paddw %%mm2, %%mm0 \n\t"
  1974. "movd 6(%0, %%ebx), %%mm4 \n\t"
  1975. "movd 6(%1, %%ebx), %%mm1 \n\t"
  1976. "movd 9(%0, %%ebx), %%mm2 \n\t"
  1977. "movd 9(%1, %%ebx), %%mm3 \n\t"
  1978. "punpcklbw %%mm7, %%mm4 \n\t"
  1979. "punpcklbw %%mm7, %%mm1 \n\t"
  1980. "punpcklbw %%mm7, %%mm2 \n\t"
  1981. "punpcklbw %%mm7, %%mm3 \n\t"
  1982. "paddw %%mm1, %%mm4 \n\t"
  1983. "paddw %%mm3, %%mm2 \n\t"
  1984. "paddw %%mm4, %%mm2 \n\t"
  1985. "psrlw $2, %%mm0 \n\t"
  1986. "psrlw $2, %%mm2 \n\t"
  1987. #endif
  1988. "movq "MANGLE(bgr2VCoeff)", %%mm1 \n\t"
  1989. "movq "MANGLE(bgr2VCoeff)", %%mm3 \n\t"
  1990. "pmaddwd %%mm0, %%mm1 \n\t"
  1991. "pmaddwd %%mm2, %%mm3 \n\t"
  1992. "pmaddwd %%mm6, %%mm0 \n\t"
  1993. "pmaddwd %%mm6, %%mm2 \n\t"
  1994. #ifndef FAST_BGR2YV12
  1995. "psrad $8, %%mm0 \n\t"
  1996. "psrad $8, %%mm1 \n\t"
  1997. "psrad $8, %%mm2 \n\t"
  1998. "psrad $8, %%mm3 \n\t"
  1999. #endif
  2000. "packssdw %%mm2, %%mm0 \n\t"
  2001. "packssdw %%mm3, %%mm1 \n\t"
  2002. "pmaddwd %%mm5, %%mm0 \n\t"
  2003. "pmaddwd %%mm5, %%mm1 \n\t"
  2004. "packssdw %%mm1, %%mm0 \n\t" // V1 V0 U1 U0
  2005. "psraw $7, %%mm0 \n\t"
  2006. #if defined (HAVE_MMX2) || defined (HAVE_3DNOW)
  2007. "movq 12(%0, %%ebx), %%mm4 \n\t"
  2008. "movq 12(%1, %%ebx), %%mm1 \n\t"
  2009. "movq 18(%0, %%ebx), %%mm2 \n\t"
  2010. "movq 18(%1, %%ebx), %%mm3 \n\t"
  2011. PAVGB" %%mm1, %%mm4 \n\t"
  2012. PAVGB" %%mm3, %%mm2 \n\t"
  2013. "movq %%mm4, %%mm1 \n\t"
  2014. "movq %%mm2, %%mm3 \n\t"
  2015. "psrlq $24, %%mm4 \n\t"
  2016. "psrlq $24, %%mm2 \n\t"
  2017. PAVGB" %%mm1, %%mm4 \n\t"
  2018. PAVGB" %%mm3, %%mm2 \n\t"
  2019. "punpcklbw %%mm7, %%mm4 \n\t"
  2020. "punpcklbw %%mm7, %%mm2 \n\t"
  2021. #else
  2022. "movd 12(%0, %%ebx), %%mm4 \n\t"
  2023. "movd 12(%1, %%ebx), %%mm1 \n\t"
  2024. "movd 15(%0, %%ebx), %%mm2 \n\t"
  2025. "movd 15(%1, %%ebx), %%mm3 \n\t"
  2026. "punpcklbw %%mm7, %%mm4 \n\t"
  2027. "punpcklbw %%mm7, %%mm1 \n\t"
  2028. "punpcklbw %%mm7, %%mm2 \n\t"
  2029. "punpcklbw %%mm7, %%mm3 \n\t"
  2030. "paddw %%mm1, %%mm4 \n\t"
  2031. "paddw %%mm3, %%mm2 \n\t"
  2032. "paddw %%mm2, %%mm4 \n\t"
  2033. "movd 18(%0, %%ebx), %%mm5 \n\t"
  2034. "movd 18(%1, %%ebx), %%mm1 \n\t"
  2035. "movd 21(%0, %%ebx), %%mm2 \n\t"
  2036. "movd 21(%1, %%ebx), %%mm3 \n\t"
  2037. "punpcklbw %%mm7, %%mm5 \n\t"
  2038. "punpcklbw %%mm7, %%mm1 \n\t"
  2039. "punpcklbw %%mm7, %%mm2 \n\t"
  2040. "punpcklbw %%mm7, %%mm3 \n\t"
  2041. "paddw %%mm1, %%mm5 \n\t"
  2042. "paddw %%mm3, %%mm2 \n\t"
  2043. "paddw %%mm5, %%mm2 \n\t"
  2044. "movq "MANGLE(w1111)", %%mm5 \n\t"
  2045. "psrlw $2, %%mm4 \n\t"
  2046. "psrlw $2, %%mm2 \n\t"
  2047. #endif
  2048. "movq "MANGLE(bgr2VCoeff)", %%mm1 \n\t"
  2049. "movq "MANGLE(bgr2VCoeff)", %%mm3 \n\t"
  2050. "pmaddwd %%mm4, %%mm1 \n\t"
  2051. "pmaddwd %%mm2, %%mm3 \n\t"
  2052. "pmaddwd %%mm6, %%mm4 \n\t"
  2053. "pmaddwd %%mm6, %%mm2 \n\t"
  2054. #ifndef FAST_BGR2YV12
  2055. "psrad $8, %%mm4 \n\t"
  2056. "psrad $8, %%mm1 \n\t"
  2057. "psrad $8, %%mm2 \n\t"
  2058. "psrad $8, %%mm3 \n\t"
  2059. #endif
  2060. "packssdw %%mm2, %%mm4 \n\t"
  2061. "packssdw %%mm3, %%mm1 \n\t"
  2062. "pmaddwd %%mm5, %%mm4 \n\t"
  2063. "pmaddwd %%mm5, %%mm1 \n\t"
  2064. "addl $24, %%ebx \n\t"
  2065. "packssdw %%mm1, %%mm4 \n\t" // V3 V2 U3 U2
  2066. "psraw $7, %%mm4 \n\t"
  2067. "movq %%mm0, %%mm1 \n\t"
  2068. "punpckldq %%mm4, %%mm0 \n\t"
  2069. "punpckhdq %%mm4, %%mm1 \n\t"
  2070. "packsswb %%mm1, %%mm0 \n\t"
  2071. "paddb "MANGLE(bgr2UVOffset)", %%mm0 \n\t"
  2072. "movd %%mm0, (%2, %%eax) \n\t"
  2073. "punpckhdq %%mm0, %%mm0 \n\t"
  2074. "movd %%mm0, (%3, %%eax) \n\t"
  2075. "addl $4, %%eax \n\t"
  2076. " js 1b \n\t"
  2077. : : "r" (src+width*6), "r" (src+srcStride+width*6), "r" (udst+width), "r" (vdst+width), "g" (-width)
  2078. : "%eax", "%ebx"
  2079. );
  2080. udst += chromStride;
  2081. vdst += chromStride;
  2082. src += srcStride*2;
  2083. }
  2084. asm volatile( EMMS" \n\t"
  2085. SFENCE" \n\t"
  2086. :::"memory");
  2087. #else
  2088. y=0;
  2089. #endif
  2090. for(; y<height; y+=2)
  2091. {
  2092. unsigned i;
  2093. for(i=0; i<chromWidth; i++)
  2094. {
  2095. unsigned int b= src[6*i+0];
  2096. unsigned int g= src[6*i+1];
  2097. unsigned int r= src[6*i+2];
  2098. unsigned int Y = ((RY*r + GY*g + BY*b)>>RGB2YUV_SHIFT) + 16;
  2099. unsigned int V = ((RV*r + GV*g + BV*b)>>RGB2YUV_SHIFT) + 128;
  2100. unsigned int U = ((RU*r + GU*g + BU*b)>>RGB2YUV_SHIFT) + 128;
  2101. udst[i] = U;
  2102. vdst[i] = V;
  2103. ydst[2*i] = Y;
  2104. b= src[6*i+3];
  2105. g= src[6*i+4];
  2106. r= src[6*i+5];
  2107. Y = ((RY*r + GY*g + BY*b)>>RGB2YUV_SHIFT) + 16;
  2108. ydst[2*i+1] = Y;
  2109. }
  2110. ydst += lumStride;
  2111. src += srcStride;
  2112. for(i=0; i<chromWidth; i++)
  2113. {
  2114. unsigned int b= src[6*i+0];
  2115. unsigned int g= src[6*i+1];
  2116. unsigned int r= src[6*i+2];
  2117. unsigned int Y = ((RY*r + GY*g + BY*b)>>RGB2YUV_SHIFT) + 16;
  2118. ydst[2*i] = Y;
  2119. b= src[6*i+3];
  2120. g= src[6*i+4];
  2121. r= src[6*i+5];
  2122. Y = ((RY*r + GY*g + BY*b)>>RGB2YUV_SHIFT) + 16;
  2123. ydst[2*i+1] = Y;
  2124. }
  2125. udst += chromStride;
  2126. vdst += chromStride;
  2127. ydst += lumStride;
  2128. src += srcStride;
  2129. }
  2130. }
  2131. void RENAME(interleaveBytes)(uint8_t *src1, uint8_t *src2, uint8_t *dest,
  2132. unsigned width, unsigned height, int src1Stride,
  2133. int src2Stride, int dstStride){
  2134. unsigned h;
  2135. for(h=0; h < height; h++)
  2136. {
  2137. unsigned w;
  2138. #ifdef HAVE_MMX
  2139. #ifdef HAVE_SSE2
  2140. asm(
  2141. "xorl %%eax, %%eax \n\t"
  2142. "1: \n\t"
  2143. PREFETCH" 64(%1, %%eax) \n\t"
  2144. PREFETCH" 64(%2, %%eax) \n\t"
  2145. "movdqa (%1, %%eax), %%xmm0 \n\t"
  2146. "movdqa (%1, %%eax), %%xmm1 \n\t"
  2147. "movdqa (%2, %%eax), %%xmm2 \n\t"
  2148. "punpcklbw %%xmm2, %%xmm0 \n\t"
  2149. "punpckhbw %%xmm2, %%xmm1 \n\t"
  2150. "movntdq %%xmm0, (%0, %%eax, 2) \n\t"
  2151. "movntdq %%xmm1, 16(%0, %%eax, 2)\n\t"
  2152. "addl $16, %%eax \n\t"
  2153. "cmpl %3, %%eax \n\t"
  2154. " jb 1b \n\t"
  2155. ::"r"(dest), "r"(src1), "r"(src2), "r" (width-15)
  2156. : "memory", "%eax"
  2157. );
  2158. #else
  2159. asm(
  2160. "xorl %%eax, %%eax \n\t"
  2161. "1: \n\t"
  2162. PREFETCH" 64(%1, %%eax) \n\t"
  2163. PREFETCH" 64(%2, %%eax) \n\t"
  2164. "movq (%1, %%eax), %%mm0 \n\t"
  2165. "movq 8(%1, %%eax), %%mm2 \n\t"
  2166. "movq %%mm0, %%mm1 \n\t"
  2167. "movq %%mm2, %%mm3 \n\t"
  2168. "movq (%2, %%eax), %%mm4 \n\t"
  2169. "movq 8(%2, %%eax), %%mm5 \n\t"
  2170. "punpcklbw %%mm4, %%mm0 \n\t"
  2171. "punpckhbw %%mm4, %%mm1 \n\t"
  2172. "punpcklbw %%mm5, %%mm2 \n\t"
  2173. "punpckhbw %%mm5, %%mm3 \n\t"
  2174. MOVNTQ" %%mm0, (%0, %%eax, 2) \n\t"
  2175. MOVNTQ" %%mm1, 8(%0, %%eax, 2) \n\t"
  2176. MOVNTQ" %%mm2, 16(%0, %%eax, 2) \n\t"
  2177. MOVNTQ" %%mm3, 24(%0, %%eax, 2) \n\t"
  2178. "addl $16, %%eax \n\t"
  2179. "cmpl %3, %%eax \n\t"
  2180. " jb 1b \n\t"
  2181. ::"r"(dest), "r"(src1), "r"(src2), "r" (width-15)
  2182. : "memory", "%eax"
  2183. );
  2184. #endif
  2185. for(w= (width&(~15)); w < width; w++)
  2186. {
  2187. dest[2*w+0] = src1[w];
  2188. dest[2*w+1] = src2[w];
  2189. }
  2190. #else
  2191. for(w=0; w < width; w++)
  2192. {
  2193. dest[2*w+0] = src1[w];
  2194. dest[2*w+1] = src2[w];
  2195. }
  2196. #endif
  2197. dest += dstStride;
  2198. src1 += src1Stride;
  2199. src2 += src2Stride;
  2200. }
  2201. #ifdef HAVE_MMX
  2202. asm(
  2203. EMMS" \n\t"
  2204. SFENCE" \n\t"
  2205. ::: "memory"
  2206. );
  2207. #endif
  2208. }
  2209. static inline void RENAME(vu9_to_vu12)(const uint8_t *src1, const uint8_t *src2,
  2210. uint8_t *dst1, uint8_t *dst2,
  2211. unsigned width, unsigned height,
  2212. int srcStride1, int srcStride2,
  2213. int dstStride1, int dstStride2)
  2214. {
  2215. unsigned int y,x,h;
  2216. int w;
  2217. w=width/2; h=height/2;
  2218. #ifdef HAVE_MMX
  2219. asm volatile(
  2220. PREFETCH" %0\n\t"
  2221. PREFETCH" %1\n\t"
  2222. ::"m"(*(src1+srcStride1)),"m"(*(src2+srcStride2)):"memory");
  2223. #endif
  2224. for(y=0;y<h;y++){
  2225. const uint8_t* s1=src1+srcStride1*(y>>1);
  2226. uint8_t* d=dst1+dstStride1*y;
  2227. x=0;
  2228. #ifdef HAVE_MMX
  2229. for(;x<w-31;x+=32)
  2230. {
  2231. asm volatile(
  2232. PREFETCH" 32%1\n\t"
  2233. "movq %1, %%mm0\n\t"
  2234. "movq 8%1, %%mm2\n\t"
  2235. "movq 16%1, %%mm4\n\t"
  2236. "movq 24%1, %%mm6\n\t"
  2237. "movq %%mm0, %%mm1\n\t"
  2238. "movq %%mm2, %%mm3\n\t"
  2239. "movq %%mm4, %%mm5\n\t"
  2240. "movq %%mm6, %%mm7\n\t"
  2241. "punpcklbw %%mm0, %%mm0\n\t"
  2242. "punpckhbw %%mm1, %%mm1\n\t"
  2243. "punpcklbw %%mm2, %%mm2\n\t"
  2244. "punpckhbw %%mm3, %%mm3\n\t"
  2245. "punpcklbw %%mm4, %%mm4\n\t"
  2246. "punpckhbw %%mm5, %%mm5\n\t"
  2247. "punpcklbw %%mm6, %%mm6\n\t"
  2248. "punpckhbw %%mm7, %%mm7\n\t"
  2249. MOVNTQ" %%mm0, %0\n\t"
  2250. MOVNTQ" %%mm1, 8%0\n\t"
  2251. MOVNTQ" %%mm2, 16%0\n\t"
  2252. MOVNTQ" %%mm3, 24%0\n\t"
  2253. MOVNTQ" %%mm4, 32%0\n\t"
  2254. MOVNTQ" %%mm5, 40%0\n\t"
  2255. MOVNTQ" %%mm6, 48%0\n\t"
  2256. MOVNTQ" %%mm7, 56%0"
  2257. :"=m"(d[2*x])
  2258. :"m"(s1[x])
  2259. :"memory");
  2260. }
  2261. #endif
  2262. for(;x<w;x++) d[2*x]=d[2*x+1]=s1[x];
  2263. }
  2264. for(y=0;y<h;y++){
  2265. const uint8_t* s2=src2+srcStride2*(y>>1);
  2266. uint8_t* d=dst2+dstStride2*y;
  2267. x=0;
  2268. #ifdef HAVE_MMX
  2269. for(;x<w-31;x+=32)
  2270. {
  2271. asm volatile(
  2272. PREFETCH" 32%1\n\t"
  2273. "movq %1, %%mm0\n\t"
  2274. "movq 8%1, %%mm2\n\t"
  2275. "movq 16%1, %%mm4\n\t"
  2276. "movq 24%1, %%mm6\n\t"
  2277. "movq %%mm0, %%mm1\n\t"
  2278. "movq %%mm2, %%mm3\n\t"
  2279. "movq %%mm4, %%mm5\n\t"
  2280. "movq %%mm6, %%mm7\n\t"
  2281. "punpcklbw %%mm0, %%mm0\n\t"
  2282. "punpckhbw %%mm1, %%mm1\n\t"
  2283. "punpcklbw %%mm2, %%mm2\n\t"
  2284. "punpckhbw %%mm3, %%mm3\n\t"
  2285. "punpcklbw %%mm4, %%mm4\n\t"
  2286. "punpckhbw %%mm5, %%mm5\n\t"
  2287. "punpcklbw %%mm6, %%mm6\n\t"
  2288. "punpckhbw %%mm7, %%mm7\n\t"
  2289. MOVNTQ" %%mm0, %0\n\t"
  2290. MOVNTQ" %%mm1, 8%0\n\t"
  2291. MOVNTQ" %%mm2, 16%0\n\t"
  2292. MOVNTQ" %%mm3, 24%0\n\t"
  2293. MOVNTQ" %%mm4, 32%0\n\t"
  2294. MOVNTQ" %%mm5, 40%0\n\t"
  2295. MOVNTQ" %%mm6, 48%0\n\t"
  2296. MOVNTQ" %%mm7, 56%0"
  2297. :"=m"(d[2*x])
  2298. :"m"(s2[x])
  2299. :"memory");
  2300. }
  2301. #endif
  2302. for(;x<w;x++) d[2*x]=d[2*x+1]=s2[x];
  2303. }
  2304. #ifdef HAVE_MMX
  2305. asm(
  2306. EMMS" \n\t"
  2307. SFENCE" \n\t"
  2308. ::: "memory"
  2309. );
  2310. #endif
  2311. }
  2312. static inline void RENAME(yvu9_to_yuy2)(const uint8_t *src1, const uint8_t *src2, const uint8_t *src3,
  2313. uint8_t *dst,
  2314. unsigned width, unsigned height,
  2315. int srcStride1, int srcStride2,
  2316. int srcStride3, int dstStride)
  2317. {
  2318. unsigned y,x,w,h;
  2319. w=width/2; h=height;
  2320. for(y=0;y<h;y++){
  2321. const uint8_t* yp=src1+srcStride1*y;
  2322. const uint8_t* up=src2+srcStride2*(y>>2);
  2323. const uint8_t* vp=src3+srcStride3*(y>>2);
  2324. uint8_t* d=dst+dstStride*y;
  2325. x=0;
  2326. #ifdef HAVE_MMX
  2327. for(;x<w-7;x+=8)
  2328. {
  2329. asm volatile(
  2330. PREFETCH" 32(%1, %0)\n\t"
  2331. PREFETCH" 32(%2, %0)\n\t"
  2332. PREFETCH" 32(%3, %0)\n\t"
  2333. "movq (%1, %0, 4), %%mm0\n\t" /* Y0Y1Y2Y3Y4Y5Y6Y7 */
  2334. "movq (%2, %0), %%mm1\n\t" /* U0U1U2U3U4U5U6U7 */
  2335. "movq (%3, %0), %%mm2\n\t" /* V0V1V2V3V4V5V6V7 */
  2336. "movq %%mm0, %%mm3\n\t" /* Y0Y1Y2Y3Y4Y5Y6Y7 */
  2337. "movq %%mm1, %%mm4\n\t" /* U0U1U2U3U4U5U6U7 */
  2338. "movq %%mm2, %%mm5\n\t" /* V0V1V2V3V4V5V6V7 */
  2339. "punpcklbw %%mm1, %%mm1\n\t" /* U0U0 U1U1 U2U2 U3U3 */
  2340. "punpcklbw %%mm2, %%mm2\n\t" /* V0V0 V1V1 V2V2 V3V3 */
  2341. "punpckhbw %%mm4, %%mm4\n\t" /* U4U4 U5U5 U6U6 U7U7 */
  2342. "punpckhbw %%mm5, %%mm5\n\t" /* V4V4 V5V5 V6V6 V7V7 */
  2343. "movq %%mm1, %%mm6\n\t"
  2344. "punpcklbw %%mm2, %%mm1\n\t" /* U0V0 U0V0 U1V1 U1V1*/
  2345. "punpcklbw %%mm1, %%mm0\n\t" /* Y0U0 Y1V0 Y2U0 Y3V0*/
  2346. "punpckhbw %%mm1, %%mm3\n\t" /* Y4U1 Y5V1 Y6U1 Y7V1*/
  2347. MOVNTQ" %%mm0, (%4, %0, 8)\n\t"
  2348. MOVNTQ" %%mm3, 8(%4, %0, 8)\n\t"
  2349. "punpckhbw %%mm2, %%mm6\n\t" /* U2V2 U2V2 U3V3 U3V3*/
  2350. "movq 8(%1, %0, 4), %%mm0\n\t"
  2351. "movq %%mm0, %%mm3\n\t"
  2352. "punpcklbw %%mm6, %%mm0\n\t" /* Y U2 Y V2 Y U2 Y V2*/
  2353. "punpckhbw %%mm6, %%mm3\n\t" /* Y U3 Y V3 Y U3 Y V3*/
  2354. MOVNTQ" %%mm0, 16(%4, %0, 8)\n\t"
  2355. MOVNTQ" %%mm3, 24(%4, %0, 8)\n\t"
  2356. "movq %%mm4, %%mm6\n\t"
  2357. "movq 16(%1, %0, 4), %%mm0\n\t"
  2358. "movq %%mm0, %%mm3\n\t"
  2359. "punpcklbw %%mm5, %%mm4\n\t"
  2360. "punpcklbw %%mm4, %%mm0\n\t" /* Y U4 Y V4 Y U4 Y V4*/
  2361. "punpckhbw %%mm4, %%mm3\n\t" /* Y U5 Y V5 Y U5 Y V5*/
  2362. MOVNTQ" %%mm0, 32(%4, %0, 8)\n\t"
  2363. MOVNTQ" %%mm3, 40(%4, %0, 8)\n\t"
  2364. "punpckhbw %%mm5, %%mm6\n\t"
  2365. "movq 24(%1, %0, 4), %%mm0\n\t"
  2366. "movq %%mm0, %%mm3\n\t"
  2367. "punpcklbw %%mm6, %%mm0\n\t" /* Y U6 Y V6 Y U6 Y V6*/
  2368. "punpckhbw %%mm6, %%mm3\n\t" /* Y U7 Y V7 Y U7 Y V7*/
  2369. MOVNTQ" %%mm0, 48(%4, %0, 8)\n\t"
  2370. MOVNTQ" %%mm3, 56(%4, %0, 8)\n\t"
  2371. : "+r" (x)
  2372. : "r"(yp), "r" (up), "r"(vp), "r"(d)
  2373. :"memory");
  2374. }
  2375. #endif
  2376. for(; x<w; x++)
  2377. {
  2378. const int x2= x<<2;
  2379. d[8*x+0]=yp[x2];
  2380. d[8*x+1]=up[x];
  2381. d[8*x+2]=yp[x2+1];
  2382. d[8*x+3]=vp[x];
  2383. d[8*x+4]=yp[x2+2];
  2384. d[8*x+5]=up[x];
  2385. d[8*x+6]=yp[x2+3];
  2386. d[8*x+7]=vp[x];
  2387. }
  2388. }
  2389. #ifdef HAVE_MMX
  2390. asm(
  2391. EMMS" \n\t"
  2392. SFENCE" \n\t"
  2393. ::: "memory"
  2394. );
  2395. #endif
  2396. }