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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > 流媒体/Mpeg4/MP4 > 查看源码
interpolate8x8_mmx.asm
资源名称：NETVIDEO.rar　[点击查看]
上传用户：sun1608
上传日期：2007-02-02
资源大小：6116k
文件大小：18k
源码类别：

流媒体/Mpeg4/MP4

开发平台：

Visual C++

interpolate8x8_mmx.asm：源码内容
  1. ;/**************************************************************************
  2. ; *
  3. ; * XVID MPEG-4 VIDEO CODEC
  4. ; * mmx/xmm/3dnow 8x8 block-based halfpel interpolation
  5. ; *
  6. ; * This program is an implementation of a part of one or more MPEG-4
  7. ; * Video tools as specified in ISO/IEC 14496-2 standard.  Those intending
  8. ; * to use this software module in hardware or software products are
  9. ; * advised that its use may infringe existing patents or copyrights, and
  10. ; * any such use would be at such party's own risk.  The original
  11. ; * developer of this software module and his/her company, and subsequent
  12. ; * editors and their companies, will have no liability for use of this
  13. ; * software or modifications or derivatives thereof.
  14. ; *
  15. ; * This program is free software; you can redistribute it and/or modify
  16. ; * it under the terms of the GNU General Public License as published by
  17. ; * the Free Software Foundation; either version 2 of the License, or
  18. ; * (at your option) any later version.
  19. ; *
  20. ; * This program is distributed in the hope that it will be useful,
  21. ; * but WITHOUT ANY WARRANTY; without even the implied warranty of
  22. ; * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  23. ; * GNU General Public License for more details.
  24. ; *
  25. ; * You should have received a copy of the GNU General Public License
  26. ; * along with this program; if not, write to the Free Software
  27. ; * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  28. ; *
  29. ; *************************************************************************/
  31. ;/**************************************************************************
  32. ; *
  33. ; * History:
  34. ; *
  35. ; * 04.02.2002  added xmm and additional 3dnow optimizations (Isibaar)
  36. ; * 27.12.2001 mofified "compensate_halfpel"
  37. ; * 22.12.2001 inital version; (c)2001 peter ross <pross@cs.rmit.edu.au>
  38. ; *
  39. ; *************************************************************************/
  42. bits 32
  44. %macro cglobal 1 
  45. %ifdef PREFIX
  46. global _%1 
  47. %define %1 _%1
  48. %else
  49. global %1
  50. %endif
  51. %endmacro
  53. section .data
  56. align 16
  59. ;===========================================================================
  60. ; (1 - r) rounding table
  61. ;===========================================================================
  63. rounding1_mmx
  64. times 4 dw 1
  65. times 4 dw 0
  67. ;===========================================================================
  68. ; (2 - r) rounding table  
  69. ;===========================================================================
  71. rounding2_mmx
  72. times 4 dw 2
  73. times 4 dw 1
  75. mmx_one
  76. times 8 db 1
  78. section .text
  80. %macro  CALC_AVG 6
  81. punpcklbw %3, %6
  82. punpckhbw %4, %6
  84. paddusw %1, %3 ; mm01 += mm23
  85. paddusw %2, %4
  86. paddusw %1, %5 ; mm01 += rounding
  87. paddusw %2, %5
  89. psrlw %1, 1 ; mm01 >>= 1
  90. psrlw %2, 1
  92. %endmacro
  96. ;===========================================================================
  97. ;
  98. ; void interpolate8x8_halfpel_h_xmm(uint8_t * const dst,
  99. ; const uint8_t * const src,
  100. ; const uint32_t stride,
  101. ; const uint32_t rounding);
  102. ;
  103. ;===========================================================================
  105. align 16
  106. cglobal interpolate8x8_halfpel_h_xmm
  107. interpolate8x8_halfpel_h_xmm
  108. push esi
  109. push edi
  111. mov edi, [esp + 8 + 4] ; dst
  112. mov esi, [esp + 8 + 8] ; src
  113. mov edx, [esp + 8 + 12] ; stride
  115. mov ecx, edx
  116. shl edx, 1
  118. mov eax, [esp + 8 + 16] ; rounding
  120. or al,al
  121. jnz near halfpel_h_xmm_rounding1 ; needs extra hack, still faster than mmx
  123. halfpel_h_xmm_rounding0
  124. ; mov eax, 4
  126. ;.loop
  127. movq mm0, [esi]
  128. movq mm1, [esi + ecx]
  129. pavgb mm0, [esi + 1] ; mm0 = avg([src], [src+1])
  130. pavgb mm1, [esi + ecx + 1] ; mm1 = avg([src+stride], [src+stride+1])
  132. movq [edi], mm0
  133. movq [edi + ecx], mm1
  135. add esi, edx ; src += 2*stride
  136. add edi, edx ; dst += 2*stride
  138. movq mm0, [esi]
  139. movq mm1, [esi + ecx]
  140. pavgb mm0, [esi + 1]
  141. pavgb mm1, [esi + ecx + 1]
  143. movq [edi], mm0
  144. movq [edi + ecx], mm1
  146. add esi, edx
  147. add edi, edx
  149. movq mm0, [esi]
  150. movq mm1, [esi + ecx]
  151. pavgb mm0, [esi + 1]
  152. pavgb mm1, [esi + ecx + 1]
  154. movq [edi], mm0
  155. movq [edi + ecx], mm1
  157. add esi, edx
  158. add edi, edx
  160. movq mm0, [esi]
  161. movq mm1, [esi + ecx]
  162. pavgb mm0, [esi + 1]
  163. pavgb mm1, [esi + ecx + 1]
  165. movq [edi], mm0
  166. movq [edi + ecx], mm1
  168. ; dec eax
  169. ; jnz .loop
  171. pop edi
  172. pop esi
  174. ret
  176. halfpel_h_xmm_rounding1
  177. movq mm7, [mmx_one]
  178. mov eax, 4
  180. .loop
  181. movq mm0, [esi]
  182. movq mm1, [esi + ecx]
  183. movq mm2, [esi + 1]
  184. movq mm3, [esi + ecx + 1]
  185. movq mm4, mm0 ; backup mm01
  186. movq mm5, mm1
  188. pavgb mm0, mm2 ; mm0 = avg([src], [src+1])
  189. pavgb mm1, mm3 ; mm1 = avg([src+stride], [src+stride+1])
  191. paddb mm2, mm4 ; add values, see if lsb is set
  192. paddb mm3, mm5 ; (to patch pavgb's (a+b+1)/2 rounding)
  193. pand mm2, mm7 ; isolate lsb's
  194. pand mm3, mm7
  195. psubusb mm0, mm2 ; mm23 = 0 if lsb wasn't set
  196. psubusb mm1, mm3
  198. movq [edi], mm0
  199. movq [edi + ecx], mm1
  201. add esi, edx ; src += 2*stride
  202. add edi, edx ; dst += 2*stride
  204. dec eax
  205. jnz .loop
  207. pop edi
  208. pop esi
  210. ret
  213. ;===========================================================================
  214. ;
  215. ; void interpolate8x8_halfpel_h_3dn(uint8_t * const dst,
  216. ; const uint8_t * const src,
  217. ; const uint32_t stride,
  218. ; const uint32_t rounding);
  219. ;
  220. ;===========================================================================
  222. align 16
  223. cglobal interpolate8x8_halfpel_h_3dn
  224. interpolate8x8_halfpel_h_3dn
  225. push esi
  226. push edi
  228. mov edi, [esp + 8 + 4] ; dst
  229. mov esi, [esp + 8 + 8] ; src
  230. mov edx, [esp + 8 + 12] ; stride
  232. mov ecx, edx
  233. shl edx, 1
  235. mov eax, [esp + 8 + 16] ; rounding
  237. or al,al ; perform hack for pavgusb rounding issue
  238. jnz near halfpel_h_3dn_rounding1
  240. halfpel_h_3dn_rounding0
  241. ; mov eax, 4
  243. ;.loop
  244. movq mm0, [esi]
  245. movq mm1, [esi + ecx]
  246. pavgusb mm0, [esi + 1] ; mm0 = avg([src], [src+1])
  247. pavgusb mm1, [esi + ecx + 1] ; mm0 = avg([src+stride], [src+stride+1])
  249. movq [edi], mm0
  250. movq [edi + ecx], mm1
  252. add esi, edx ; src += 2*stride
  253. add edi, edx ; dst += 2*stride
  255. movq mm0, [esi]
  256. movq mm1, [esi + ecx]
  257. pavgusb mm0, [esi + 1]
  258. pavgusb mm1, [esi + ecx + 1]
  260. movq [edi], mm0
  261. movq [edi + ecx], mm1
  263. add esi, edx
  264. add edi, edx
  266. movq mm0, [esi]
  267. movq mm1, [esi + ecx]
  268. pavgusb mm0, [esi + 1]
  269. pavgusb mm1, [esi + ecx + 1]
  271. movq [edi], mm0
  272. movq [edi + ecx], mm1
  274. add esi, edx
  275. add edi, edx
  277. movq mm0, [esi]
  278. movq mm1, [esi + ecx]
  279. pavgusb mm0, [esi + 1]
  280. pavgusb mm1, [esi + ecx + 1]
  282. movq [edi], mm0
  283. movq [edi + ecx], mm1
  285. ; dec eax
  286. ; jnz .loop
  288. pop edi
  289. pop esi
  291. ret
  293. halfpel_h_3dn_rounding1
  294. movq mm7, [mmx_one]
  295. mov eax, 4
  297. .loop
  298. movq mm0, [esi]
  299. movq mm1, [esi + ecx]
  300. movq mm2, [esi + 1]
  301. movq mm3, [esi + ecx + 1]
  302. movq mm4, mm0 ; backup mm01
  303. movq mm5, mm1
  305. pavgusb mm0, mm2 ; mm0 = avg([src], [src+1])
  306. pavgusb mm1, mm3 ; mm1 = avg([src+stride], [src+stride+1])
  308. paddb mm2, mm4 ; add values, see if lsb is set
  309. paddb mm3, mm5 ; (to patch pavgb's (a+b+1)/2 rounding)
  310. pand mm2, mm7 ; isolate lsb's
  311. pand mm3, mm7
  312. psubusb mm0, mm2 ; mm23 = 0 if lsb wasn't set
  313. psubusb mm1, mm3
  315. movq [edi], mm0
  316. movq [edi + ecx], mm1
  318. add esi, edx ; src += 2*stride
  319. add edi, edx ; dst += 2*stride
  321. dec eax
  322. jnz .loop
  324. pop edi
  325. pop esi
  327. ret
  332. ;===========================================================================
  333. ;
  334. ; void interpolate8x8_halfpel_h_mmx(uint8_t * const dst,
  335. ; const uint8_t * const src,
  336. ; const uint32_t stride,
  337. ; const uint32_t rounding);
  338. ;
  339. ;===========================================================================
  341. align 16
  342. cglobal interpolate8x8_halfpel_h_mmx
  343. interpolate8x8_halfpel_h_mmx
  345. push esi
  346. push edi
  348. mov eax, [esp + 8 + 16] ; rounding
  350. interpolate8x8_halfpel_h_mmx.start
  351. movq mm7, [rounding1_mmx + eax * 8]
  353. mov edi, [esp + 8 + 4] ; dst
  354. mov esi, [esp + 8 + 8] ; src
  355. mov edx, [esp + 8 + 12] ; stride
  357. mov eax, 8
  358. pxor mm6, mm6 ; zero
  360. .loop
  361. movq mm0, [esi]
  362. movq mm2, [esi + 1]
  363. movq mm1, mm0
  364. movq mm3, mm2
  366. punpcklbw mm0, mm6 ; mm01 = [src]
  367. punpckhbw mm1, mm6 ; mm23 = [src + 1]
  369. CALC_AVG mm0, mm1, mm2, mm3, mm7, mm6
  371. packuswb mm0, mm1
  372. movq [edi], mm0 ; [dst] = mm01
  374. add esi, edx ; src += stride
  375. add edi, edx ; dst += stride
  377. dec eax
  378. jnz .loop
  380. pop edi
  381. pop esi
  383. ret
  386. ;===========================================================================
  387. ;
  388. ; void interpolate8x8_halfpel_v_xmm(uint8_t * const dst,
  389. ; const uint8_t * const src,
  390. ; const uint32_t stride,
  391. ; const uint32_t rounding);
  392. ;
  393. ;===========================================================================
  395. align 16
  396. cglobal interpolate8x8_halfpel_v_xmm
  397. interpolate8x8_halfpel_v_xmm
  398. push esi
  399. push edi
  401. mov edi, [esp + 8 + 4] ; dst
  402. mov esi, [esp + 8 + 8] ; src
  403. mov edx, [esp + 8 + 12] ; stride
  405. mov eax, [esp + 8 + 16] ; rounding
  407. or al,al
  408. jnz near halfpel_v_xmm_rounding1
  410. halfpel_v_xmm_rounding0
  411. ; mov eax, 8
  413. ;.loop
  414. movq mm0, [esi]
  415. pavgb mm0, [esi + edx] ; mm0 = avg([src], [src+stride])
  416. movq [edi], mm0
  418. add edi, edx ; dst += stride
  419. add esi, edx ; src += stride
  421. movq mm0, [esi]
  422. pavgb mm0, [esi + edx]
  423. movq [edi], mm0
  425. add edi, edx
  426. add esi, edx
  428. movq mm0, [esi]
  429. pavgb mm0, [esi + edx]
  430. movq [edi], mm0
  432. add edi, edx
  433. add esi, edx
  435. movq mm0, [esi]
  436. pavgb mm0, [esi + edx]
  437. movq [edi], mm0
  439. add edi, edx
  440. add esi, edx
  442. movq mm0, [esi]
  443. pavgb mm0, [esi + edx]
  444. movq [edi], mm0
  446. add edi, edx
  447. add esi, edx
  449. movq mm0, [esi]
  450. pavgb mm0, [esi + edx]
  451. movq [edi], mm0
  453. add edi, edx
  454. add esi, edx
  456. movq mm0, [esi]
  457. pavgb mm0, [esi + edx]
  458. movq [edi], mm0
  460. add edi, edx
  461. add esi, edx
  463. movq mm0, [esi]
  464. pavgb mm0, [esi + edx]
  465. movq [edi], mm0
  467. ; dec eax
  468. ; jnz .loop
  470. pop edi
  471. pop esi
  473. ret
  475. halfpel_v_xmm_rounding1
  476. movq mm7, [mmx_one]
  477. mov eax, 8
  479. .loop
  480. movq mm0, [esi]
  481. movq mm1, [esi + edx]
  482. movq mm2, mm0 ; backup mm0
  484. pavgb mm0, mm1 ; mm0 = avg([src], [src+stride])
  486. paddb mm1, mm2 ; isolate and subtract lsb's from original values
  487. pand mm1, mm7 ; (compensate for pavgb rounding)
  488. psubusb mm0, mm1
  490. movq [edi], mm0
  492. add edi, edx ; dst += stride
  493. add esi, edx ; src += stride
  495. dec eax
  496. jnz .loop
  498. pop edi
  499. pop esi
  501. ret
  504. ;===========================================================================
  505. ;
  506. ; void interpolate8x8_halfpel_v_3dn(uint8_t * const dst,
  507. ; const uint8_t * const src,
  508. ; const uint32_t stride,
  509. ; const uint32_t rounding);
  510. ;
  511. ;===========================================================================
  513. align 16
  514. cglobal interpolate8x8_halfpel_v_3dn
  515. interpolate8x8_halfpel_v_3dn
  516. push esi
  517. push edi
  519. mov edi, [esp + 8 + 4] ; dst
  520. mov esi, [esp + 8 + 8] ; src
  521. mov edx, [esp + 8 + 12] ; stride
  523. mov eax, [esp + 8 + 16] ; rounding
  525. or al,al ; hack for pavgusb rounding
  526. jnz halfpel_v_3dn_rounding1
  528. halfpel_v_3dn_rounding0
  529. ; mov eax, 8
  531. ;.loop
  532. movq mm0, [esi]
  533. pavgusb mm0, [esi + edx] ; mm0 = avg([src], [src+stride])
  534. movq [edi], mm0
  536. add edi, edx ; dst += stride
  537. add esi, edx ; src += stride
  539. movq mm0, [esi]
  540. pavgusb mm0, [esi + edx]
  541. movq [edi], mm0
  543. add edi, edx
  544. add esi, edx
  546. movq mm0, [esi]
  547. pavgusb mm0, [esi + edx]
  548. movq [edi], mm0
  550. add edi, edx
  551. add esi, edx
  553. movq mm0, [esi]
  554. pavgusb mm0, [esi + edx]
  555. movq [edi], mm0
  557. add edi, edx
  558. add esi, edx
  560. movq mm0, [esi]
  561. pavgusb mm0, [esi + edx]
  562. movq [edi], mm0
  564. add edi, edx
  565. add esi, edx
  567. movq mm0, [esi]
  568. pavgusb mm0, [esi + edx]
  569. movq [edi], mm0
  571. add edi, edx
  572. add esi, edx
  574. movq mm0, [esi]
  575. pavgusb mm0, [esi + edx]
  576. movq [edi], mm0
  578. add edi, edx
  579. add esi, edx
  581. movq mm0, [esi]
  582. pavgusb mm0, [esi + edx]
  583. movq [edi], mm0
  585. ; dec eax
  586. ; jnz .loop
  588. pop edi
  589. pop esi
  591. ret
  593. halfpel_v_3dn_rounding1
  594. movq mm7, [mmx_one]
  595. mov eax, 8
  597. .loop
  598. movq mm0, [esi]
  599. movq mm1, [esi + edx]
  600. movq mm2, mm0 ; backup mm0
  602. pavgusb mm0, mm1 ; mm0 = avg([src], [src+stride])
  604. paddb mm1, mm2 ; isolate and subtract lsb's from original values
  605. pand mm1, mm7 ; (compensate for pavgb rounding)
  606. psubusb mm0, mm1
  608. movq [edi], mm0
  610. add edi, edx ; dst += stride
  611. add esi, edx ; src += stride
  613. dec eax
  614. jnz .loop
  616. pop edi
  617. pop esi
  619. ret
  624. ;===========================================================================
  625. ;
  626. ; void interpolate8x8_halfpel_v_mmx(uint8_t * const dst,
  627. ; const uint8_t * const src,
  628. ; const uint32_t stride,
  629. ; const uint32_t rounding);
  630. ;
  631. ;===========================================================================
  633. align 16
  634. cglobal interpolate8x8_halfpel_v_mmx
  635. interpolate8x8_halfpel_v_mmx
  637. push esi
  638. push edi
  640. mov eax, [esp + 8 + 16] ; rounding
  642. interpolate8x8_halfpel_v_mmx.start
  643. movq mm7, [rounding1_mmx + eax * 8]
  645. mov edi, [esp + 8 + 4] ; dst
  646. mov esi, [esp + 8 + 8] ; src
  647. mov edx, [esp + 8 + 12] ; stride
  649. mov eax, 8
  651. pxor mm6, mm6 ; zero
  654. .loop
  655. movq mm0, [esi]
  656. movq mm2, [esi + edx]
  657. movq mm1, mm0
  658. movq mm3, mm2
  660. punpcklbw mm0, mm6 ; mm01 = [src]
  661. punpckhbw mm1, mm6 ; mm23 = [src + 1]
  663. CALC_AVG mm0, mm1, mm2, mm3, mm7, mm6
  665. packuswb mm0, mm1
  666. movq [edi], mm0 ; [dst] = mm01
  668. add esi, edx ; src += stride
  669. add edi, edx ; dst += stride
  671. dec eax
  672. jnz .loop
  674. pop edi
  675. pop esi
  677. ret
  680. ;===========================================================================
  681. ;
  682. ; void interpolate8x8_halfpel_hv_xmm(uint8_t * const dst,
  683. ; const uint8_t * const src,
  684. ; const uint32_t stride, 
  685. ; const uint32_t rounding);
  686. ;
  687. ;
  688. ;===========================================================================
  690. align 16
  691. cglobal interpolate8x8_halfpel_hv_xmm
  692. interpolate8x8_halfpel_hv_xmm
  694. push esi
  695. push edi
  697. mov eax, [esp + 8 + 16] ; rounding
  698. or al,al
  699. jnz interpolate8x8_halfpel_hv_mmx.start
  701. mov edi, [esp + 8 + 4] ; dst
  702. mov esi, [esp + 8 + 8] ; src
  703. mov edx, [esp + 8 + 12] ; stride
  705. movq mm7, [mmx_one]
  706. mov eax, 8
  708. .loop
  709. ; current row
  711. movq mm0, [esi]
  712. movq mm1, [esi + 1]
  714. movq mm2, mm0
  716. pavgb mm0, mm1
  718. pxor mm1, mm2
  720. ; next row
  722. movq mm2, [esi + edx]
  723. movq mm3, [esi + edx + 1]
  725. movq mm4, mm2
  727. pavgb mm2, mm3
  729. pxor mm3, mm4
  731. ; add current + next row
  733. por mm1, mm3
  734. movq mm3, mm0
  736. pxor mm3, mm2
  737. pand mm1, mm3
  739. pavgb mm0, mm2
  741. pand mm1, mm7
  742. psubusb mm0, mm1
  744. movq [edi], mm0 ; [dst] = mm01
  746. add esi, edx ; src += stride
  747. add edi, edx ; dst += stride
  749. dec eax
  750. jnz .loop
  752. pop edi
  753. pop esi
  755. ret
  759. ;===========================================================================
  760. ;
  761. ; void interpolate8x8_halfpel_hv_mmx(uint8_t * const dst,
  762. ; const uint8_t * const src,
  763. ; const uint32_t stride, 
  764. ; const uint32_t rounding);
  765. ;
  766. ;
  767. ;===========================================================================
  769. align 16
  770. cglobal interpolate8x8_halfpel_hv_mmx
  771. interpolate8x8_halfpel_hv_mmx
  773. push esi
  774. push edi
  776. mov eax, [esp + 8 + 16] ; rounding
  777. interpolate8x8_halfpel_hv_mmx.start
  779. movq mm7, [rounding2_mmx + eax * 8]
  781. mov edi, [esp + 8 + 4] ; dst
  782. mov esi, [esp + 8 + 8] ; src
  784. mov eax, 8
  786. pxor mm6, mm6 ; zero
  788. mov edx, [esp + 8 + 12] ; stride
  790. .loop
  791. ; current row
  793. movq mm0, [esi]
  794. movq mm2, [esi + 1]
  796. movq mm1, mm0
  797. movq mm3, mm2
  799. punpcklbw mm0, mm6 ; mm01 = [src]
  800. punpcklbw mm2, mm6 ; mm23 = [src + 1]
  801. punpckhbw mm1, mm6
  802. punpckhbw mm3, mm6
  804. paddusw mm0, mm2 ; mm01 += mm23
  805. paddusw mm1, mm3
  807. ; next row
  809. movq mm4, [esi + edx]
  810. movq mm2, [esi + edx + 1]
  812. movq mm5, mm4
  813. movq mm3, mm2
  815. punpcklbw mm4, mm6 ; mm45 = [src + stride]
  816. punpcklbw mm2, mm6 ; mm23 = [src + stride + 1]
  817. punpckhbw mm5, mm6
  818. punpckhbw mm3, mm6
  820. paddusw mm4, mm2 ; mm45 += mm23
  821. paddusw mm5, mm3
  823. ; add current + next row
  825. paddusw mm0, mm4 ; mm01 += mm45
  826. paddusw mm1, mm5
  827. paddusw mm0, mm7 ; mm01 += rounding2
  828. paddusw mm1, mm7
  830. psrlw mm0, 2 ; mm01 >>= 2
  831. psrlw mm1, 2
  833. packuswb mm0, mm1
  834. movq [edi], mm0 ; [dst] = mm01
  836. add esi, edx ; src += stride
  837. add edi, edx ; dst += stride
  839. dec eax
  840. jnz .loop
  842. pop edi
  843. pop esi
  845. ret
  847. ;===========================================================================
  848. ;
  849. ; void interpolate8x8_halfpel_hv_3dn(uint8_t * const dst,
  850. ; const uint8_t * const src,
  851. ; const uint32_t stride, 
  852. ; const uint32_t rounding);
  853. ;
  854. ;
  855. ;===========================================================================
  857. align 16
  858. cglobal interpolate8x8_halfpel_hv_3dn
  859. interpolate8x8_halfpel_hv_3dn
  861. push esi
  862. push edi
  864. mov eax, [esp + 8 + 16] ; rounding
  865.                 or      al,al
  866.                 jnz     near interpolate8x8_halfpel_hv_mmx.start
  868. mov edi, [esp + 8 + 4] ; dst
  869. mov esi, [esp + 8 + 8] ; src
  870. mov edx, [esp + 8 + 12] ; stride
  872. movq mm7, [mmx_one]
  873. mov eax, 8
  875. .loop
  876. ; current row
  878. movq mm0, [esi]
  879. movq mm1, [esi + 1]
  881. movq mm2, mm0
  883. pavgusb mm0, mm1
  885. pxor mm1, mm2
  887. ; next row
  889. movq mm2, [esi + edx]
  890. movq mm3, [esi + edx + 1]
  892. movq mm4, mm2
  894. pavgusb mm2, mm3
  896. pxor mm3, mm4
  898. ; add current + next row
  900. por mm1, mm3
  901. movq mm3, mm0
  903. pxor mm3, mm2
  904. pand mm1, mm3
  906. pavgusb mm0, mm2
  908. pand mm1, mm7
  909. psubusb mm0, mm1
  911. movq [edi], mm0 ; [dst] = mm01
  913. add esi, edx ; src += stride
  914. add edi, edx ; dst += stride
  916. dec eax
  917. jnz .loop
  919. pop edi
  920. pop esi
  922. ret