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

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

流媒体/Mpeg4/MP4

开发平台：

Visual C++

sad_mmx.asm：源码内容
  1. ;/**************************************************************************
  2. ; *
  3. ; * XVID MPEG-4 VIDEO CODEC
  4. ; * mmx/xmm sum of absolute difference
  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. ; * 17.11.2001  bugfix and small improvement for dev16_xmm,
  36. ; *             removed terminate early in sad16_xmm 
  37. ; * 12.11.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
  57. mmx_one times 4 dw 1
  60. section .text
  63. ;===========================================================================
  64. ;
  65. ; uint32_t sad16_mmx(const uint8_t * const cur,
  66. ; const uint8_t * const ref,
  67. ; const uint32_t stride,
  68. ; const uint32_t best_sad);
  69. ;
  70. ; (early termination ignore; slows this down)
  71. ;
  72. ;===========================================================================
  74. align 16
  75. cglobal sad16_mmx
  76. sad16_mmx
  77. push esi
  78. push edi
  80. mov esi, [esp + 8 + 4] ; ref
  81. mov edi, [esp + 8 + 8] ; cur
  82. mov ecx, [esp + 8 + 12] ; stride
  83. mov edx, 16
  85. pxor mm6, mm6 ; mm6 = sum = 0
  86. pxor mm7, mm7 ; mm7 = 0
  87. .loop
  88. movq mm0, [esi] ; ref
  89. movq mm1, [edi] ; cur
  91. movq mm2, [esi+8] ; ref2
  92. movq mm3, [edi+8] ; cur2
  94. movq mm4, mm0 
  95. movq mm5, mm2
  97. psubusb mm0, mm1
  98. psubusb mm2, mm3
  100. psubusb mm1, mm4
  101. psubusb mm3, mm5
  103. por mm0, mm1 ; mm0 = |ref - cur|
  104. por mm2, mm3 ; mm2 = |ref2 - cur2|
  106. movq mm1,mm0
  107. movq mm3,mm2
  109. punpcklbw mm0,mm7
  110. punpcklbw mm2,mm7
  112. punpckhbw mm1,mm7
  113. punpckhbw mm3,mm7
  115. paddusw mm0,mm1
  116. paddusw mm2,mm3
  118. paddusw mm6,mm0 ; sum += mm01
  119. paddusw mm6,mm2 ; sum += mm23
  121. add esi, ecx
  122. add edi, ecx
  123. dec edx
  124. jnz .loop
  126. pmaddwd mm6, [mmx_one] ; merge sum
  127. movq mm7, mm6
  128. psrlq mm7, 32 
  129. paddd mm6, mm7
  130. movd eax, mm6
  132. pop  edi
  133. pop  esi
  135. ret
  138. ;===========================================================================
  139. ;
  140. ; uint32_t sad16_xmm(const uint8_t * const cur,
  141. ; const uint8_t * const ref,
  142. ; const uint32_t stride,
  143. ; const uint32_t best_sad);
  144. ;
  145. ; experimental!
  146. ;
  147. ;===========================================================================
  149. align 16
  150. cglobal sad16_xmm
  151. sad16_xmm
  152. push esi
  153. push edi
  154. push ebx
  156. mov esi, [esp + 12 + 4] ; ref
  157. mov edi, [esp + 12 + 8] ; cur
  158. mov ecx, [esp + 12 + 12] ; stride
  159. mov ebx, [esp + 12 + 16] ; best_sad
  160. ; mov edx, 16
  162. pxor mm6, mm6 ; mm6 = sum = 0
  163. ;.loop
  164. movq mm0, [esi] ; ref
  165. movq mm2, [esi+8] ; ref2
  167. psadbw mm0, [edi] ; mm0 = |ref - cur|
  168. psadbw mm2, [edi+8] ; mm0 = |ref2 - cur2|
  170. paddusw mm6,mm0 ; sum += mm01
  171. paddusw mm6,mm2 ; sum += mm23
  173. add esi, ecx
  174. add edi, ecx
  175. ; dec edx
  176. ; jnz .loop
  178. movq mm0, [esi]
  179. movq mm2, [esi+8]
  181. psadbw mm0, [edi]
  182. psadbw mm2, [edi+8]
  184. paddusw mm6,mm0
  185. paddusw mm6,mm2
  187. add esi, ecx
  188. add edi, ecx
  190. movq mm0, [esi]
  191. movq mm2, [esi+8]
  193. psadbw mm0, [edi]
  194. psadbw mm2, [edi+8]
  196. paddusw mm6,mm0
  197. paddusw mm6,mm2
  199. add esi, ecx
  200. add edi, ecx
  202. movq mm0, [esi]
  203. movq mm2, [esi+8]
  205. psadbw mm0, [edi]
  206. psadbw mm2, [edi+8]
  208. paddusw mm6,mm0
  209. paddusw mm6,mm2
  211. add esi, ecx
  212. add edi, ecx
  214. movq mm0, [esi]
  215. movq mm2, [esi+8]
  217. psadbw mm0, [edi]
  218. psadbw mm2, [edi+8]
  220. paddusw mm6,mm0
  221. paddusw mm6,mm2
  223. add esi, ecx
  224. add edi, ecx
  226. movq mm0, [esi]
  227. movq mm2, [esi+8]
  229. psadbw mm0, [edi]
  230. psadbw mm2, [edi+8]
  232. paddusw mm6,mm0
  233. paddusw mm6,mm2
  235. add esi, ecx
  236. add edi, ecx
  238. movq mm0, [esi]
  239. movq mm2, [esi+8]
  241. psadbw mm0, [edi]
  242. psadbw mm2, [edi+8]
  244. paddusw mm6,mm0
  245. paddusw mm6,mm2
  247. add esi, ecx
  248. add edi, ecx
  250. movq mm0, [esi]
  251. movq mm2, [esi+8]
  253. psadbw mm0, [edi]
  254. psadbw mm2, [edi+8]
  256. paddusw mm6,mm0
  257. paddusw mm6,mm2
  259. add esi, ecx
  260. add edi, ecx
  262. movq mm0, [esi]
  263. movq mm2, [esi+8]
  265. psadbw mm0, [edi]
  266. psadbw mm2, [edi+8]
  268. paddusw mm6,mm0
  269. paddusw mm6,mm2
  271. add esi, ecx
  272. add edi, ecx
  274. movq mm0, [esi]
  275. movq mm2, [esi+8]
  277. psadbw mm0, [edi]
  278. psadbw mm2, [edi+8]
  280. paddusw mm6,mm0
  281. paddusw mm6,mm2
  283. add esi, ecx
  284. add edi, ecx
  286. movq mm0, [esi]
  287. movq mm2, [esi+8]
  289. psadbw mm0, [edi]
  290. psadbw mm2, [edi+8]
  292. paddusw mm6,mm0
  293. paddusw mm6,mm2
  295. add esi, ecx
  296. add edi, ecx
  298. movq mm0, [esi]
  299. movq mm2, [esi+8]
  301. psadbw mm0, [edi]
  302. psadbw mm2, [edi+8]
  304. paddusw mm6,mm0
  305. paddusw mm6,mm2
  307. add esi, ecx
  308. add edi, ecx
  310. movq mm0, [esi]
  311. movq mm2, [esi+8]
  313. psadbw mm0, [edi]
  314. psadbw mm2, [edi+8]
  316. paddusw mm6,mm0
  317. paddusw mm6,mm2
  319. add esi, ecx
  320. add edi, ecx
  322. movq mm0, [esi]
  323. movq mm2, [esi+8]
  325. psadbw mm0, [edi]
  326. psadbw mm2, [edi+8]
  328. paddusw mm6,mm0
  329. paddusw mm6,mm2
  331. add esi, ecx
  332. add edi, ecx
  334. movq mm0, [esi]
  335. movq mm2, [esi+8]
  337. psadbw mm0, [edi]
  338. psadbw mm2, [edi+8]
  340. paddusw mm6,mm0
  341. paddusw mm6,mm2
  343. add esi, ecx
  344. add edi, ecx
  346. movq mm0, [esi]
  347. movq mm2, [esi+8]
  349. psadbw mm0, [edi]
  350. psadbw mm2, [edi+8]
  352. paddusw mm6,mm0
  353. paddusw mm6,mm2
  355. movd eax, mm6
  357. .ret pop ebx
  358. pop  edi
  359. pop  esi
  361. ret
  363. ;===========================================================================
  364. ;
  365. ; uint32_t sad16_sse2(const uint8_t * const cur,
  366. ; const uint8_t * const ref,
  367. ; const uint32_t stride,
  368. ; const uint32_t best_sad);
  369. ;
  370. ;===========================================================================
  372. align 16
  373. cglobal sad16_sse2
  374. sad16_sse2
  375. push esi
  376. push edi
  377. push ebx
  379. mov esi, [esp + 12 + 4] ; ref
  380. mov edi, [esp + 12 + 8] ; cur
  381. mov ecx, [esp + 12 + 12] ; stride
  382. mov ebx, [esp + 12 + 16] ; best_sad
  383. mov edx, 16
  385. pxor xmm2, xmm2 ; xmm2 = sum = 0
  387. ; 0
  388. movdqu xmm0, [esi] ; ref
  389. movdqu xmm1, [edi] ; cur
  390. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  391. paddusw xmm2,xmm0 ; sum += xmm0
  392. add esi, ecx
  393. add edi, ecx
  395. ; 1
  396. movdqu xmm0, [esi] ; ref
  397. movdqu xmm1, [edi] ; cur
  398. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  399. paddusw xmm2,xmm0 ; sum += xmm0
  400. add esi, ecx
  401. add edi, ecx
  403. ; 2
  404. movdqu xmm0, [esi] ; ref
  405. movdqu xmm1, [edi] ; cur
  406. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  407. paddusw xmm2,xmm0 ; sum += xmm0
  408. add esi, ecx
  409. add edi, ecx
  411. ; 3
  412. movdqu xmm0, [esi] ; ref
  413. movdqu xmm1, [edi] ; cur
  414. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  415. paddusw xmm2,xmm0 ; sum += xmm0
  416. add esi, ecx
  417. add edi, ecx
  419. ; 4
  420. movdqu xmm0, [esi] ; ref
  421. movdqu xmm1, [edi] ; cur
  422. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  423. paddusw xmm2,xmm0 ; sum += xmm0
  424. add esi, ecx
  425. add edi, ecx
  427. ; 5
  428. movdqu xmm0, [esi] ; ref
  429. movdqu xmm1, [edi] ; cur
  430. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  431. paddusw xmm2,xmm0 ; sum += xmm0
  432. add esi, ecx
  433. add edi, ecx
  435. ; 6
  436. movdqu xmm0, [esi] ; ref
  437. movdqu xmm1, [edi] ; cur
  438. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  439. paddusw xmm2,xmm0 ; sum += xmm0
  440. add esi, ecx
  441. add edi, ecx
  443. ; 7
  444. movdqu xmm0, [esi] ; ref
  445. movdqu xmm1, [edi] ; cur
  446. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  447. paddusw xmm2,xmm0 ; sum += xmm0
  448. add esi, ecx
  449. add edi, ecx
  451. ; 8
  452. movdqu xmm0, [esi] ; ref
  453. movdqu xmm1, [edi] ; cur
  454. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  455. paddusw xmm2,xmm0 ; sum += xmm0
  456. add esi, ecx
  457. add edi, ecx
  459. ; 9
  460. movdqu xmm0, [esi] ; ref
  461. movdqu xmm1, [edi] ; cur
  462. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  463. paddusw xmm2,xmm0 ; sum += xmm0
  464. add esi, ecx
  465. add edi, ecx
  467. ; 10
  468. movdqu xmm0, [esi] ; ref
  469. movdqu xmm1, [edi] ; cur
  470. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  471. paddusw xmm2,xmm0 ; sum += xmm0
  472. add esi, ecx
  473. add edi, ecx
  475. ; 11
  476. movdqu xmm0, [esi] ; ref
  477. movdqu xmm1, [edi] ; cur
  478. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  479. paddusw xmm2,xmm0 ; sum += xmm0
  480. add esi, ecx
  481. add edi, ecx
  483. ; 12
  484. movdqu xmm0, [esi] ; ref
  485. movdqu xmm1, [edi] ; cur
  486. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  487. paddusw xmm2,xmm0 ; sum += xmm0
  488. add esi, ecx
  489. add edi, ecx
  491. ; 13
  492. movdqu xmm0, [esi] ; ref
  493. movdqu xmm1, [edi] ; cur
  494. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  495. paddusw xmm2,xmm0 ; sum += xmm0
  496. add esi, ecx
  497. add edi, ecx
  499. ; 14
  500. movdqu xmm0, [esi] ; ref
  501. movdqu xmm1, [edi] ; cur
  502. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  503. paddusw xmm2,xmm0 ; sum += xmm0
  504. add esi, ecx
  505. add edi, ecx
  507. ; 15
  508. movdqu xmm0, [esi] ; ref
  509. movdqu xmm1, [edi] ; cur
  510. psadbw xmm0, xmm1 ; xmm0 = |ref - cur|
  511. paddusw xmm2,xmm0 ; sum += xmm0
  513. movd eax, xmm2
  514. psrldq xmm2, 8
  515. movd ecx, xmm2
  516. add eax, ecx
  518. .ret pop ebx
  519. pop  edi
  520. pop  esi
  522. ret
  526. ;===========================================================================
  527. ;
  528. ; uint32_t sad8_mmx(const uint8_t * const cur,
  529. ; const uint8_t * const ref,
  530. ; const uint32_t stride);
  531. ;
  532. ;===========================================================================
  533. align 16
  534. cglobal sad8_mmx
  535. sad8_mmx
  536. push  esi
  537. push  edi
  539. mov  esi, [esp + 8 + 4] ; ref
  540. mov  edi, [esp + 8 + 8] ; cur
  541. mov  ecx, [esp + 8 + 12] ; stride
  542. mov eax, 4
  544. pxor mm6, mm6 ; mm6 = sum = 0
  545. pxor mm7, mm7 ; mm7 = 0
  546. .loop
  547. movq mm0, [esi] ; ref
  548. movq mm1, [edi] ; cur
  550. movq mm2, [esi+ecx] ; ref2
  551. movq mm3, [edi+ecx] ; cur2
  553. movq mm4, mm0 
  554. movq mm5, mm2
  556. psubusb mm0, mm1
  557. psubusb mm2, mm3
  559. psubusb mm1, mm4
  560. psubusb mm3, mm5
  562. por mm0, mm1 ; mm0 = |ref - cur|
  563. por mm2, mm3 ; mm2 = |ref2 - cur2|
  565. movq mm1,mm0
  566. movq mm3,mm2
  568. punpcklbw mm0,mm7
  569. punpcklbw mm2,mm7
  571. punpckhbw mm1,mm7
  572. punpckhbw mm3,mm7
  574. paddusw mm0,mm1
  575. paddusw mm2,mm3
  577. paddusw mm6,mm0 ; sum += mm01
  578. paddusw mm6,mm2 ; sum += mm23
  580. add esi, ecx
  581. add edi, ecx
  582. add esi, ecx
  583. add edi, ecx
  584. dec eax
  585. jnz .loop
  587. pmaddwd mm6, [mmx_one] ; merge sum
  588. movq mm7, mm6
  589. psrlq mm7, 32 
  590. paddd mm6, mm7
  591. movd eax, mm6
  593. pop  edi
  594. pop  esi
  596. ret
  601. ;===========================================================================
  602. ;
  603. ; uint32_t sad8_xmm(const uint8_t * const cur,
  604. ; const uint8_t * const ref,
  605. ; const uint32_t stride);
  606. ;
  607. ; experimental!
  608. ;
  609. ;===========================================================================
  610. align 16
  611. cglobal sad8_xmm
  612. sad8_xmm
  613. push  esi
  614. push  edi
  616. mov  esi, [esp + 8 + 4] ; ref
  617. mov  edi, [esp + 8 + 8] ; cur
  618. mov  ecx, [esp + 8 + 12] ; stride
  619. mov     edx, ecx
  620. shl     edx, 1
  622. ; mov eax, 4
  624. pxor mm6, mm6 ; mm6 = sum = 0
  625. ;.loop
  626. movq mm0, [esi] ; ref
  627. movq mm2, [esi+ecx] ; ref2
  629. psadbw mm0, [edi] ; mm0 = |ref - cur|
  630. psadbw mm2, [edi+ecx] ; mm0 = |ref2 - cur2|
  632. paddusw mm6,mm0 ; sum += mm01
  633. paddusw mm6,mm2 ; sum += mm23
  635. add esi, edx
  636. add edi, edx
  637. ; dec eax
  638. ; jnz .loop
  640. movq mm0, [esi]
  641. movq mm2, [esi+ecx]
  643. psadbw mm0, [edi]
  644. psadbw mm2, [edi+ecx]
  646. paddusw mm6,mm0
  647. paddusw mm6,mm2
  649. add esi, edx
  650. add edi, edx
  652. movq mm0, [esi]
  653. movq mm2, [esi+ecx]
  655. psadbw mm0, [edi]
  656. psadbw mm2, [edi+ecx]
  658. paddusw mm6,mm0
  659. paddusw mm6,mm2
  661. add esi, edx
  662. add edi, edx
  664. movq mm0, [esi]
  665. movq mm2, [esi+ecx]
  667. psadbw mm0, [edi]
  668. psadbw mm2, [edi+ecx]
  670. paddusw mm6,mm0
  671. paddusw mm6,mm2
  673. movd eax, mm6
  675. pop  edi
  676. pop  esi
  678. ret
  682. ;===========================================================================
  683. ;
  684. ; uint32_t dev16_mmx(const uint8_t * const cur,
  685. ; const uint32_t stride);
  686. ;
  687. ;===========================================================================
  689. align 16
  690. cglobal dev16_mmx
  691. dev16_mmx
  693. push  esi
  694. push  edi
  696. pxor mm4, mm4 ; mm23 = sum = 0
  697. pxor mm5, mm5
  699. mov  esi, [esp + 8 + 4] ; cur
  700. mov  ecx, [esp + 8 + 8] ; stride
  701. mov     edi, esi
  703. mov eax, 16
  704. pxor mm7, mm7 ; mm7 = 0
  705. .loop1
  706. movq mm0, [esi]
  707. movq mm2, [esi + 8]
  709. movq mm1, mm0
  710. movq mm3, mm2
  712. punpcklbw mm0, mm7
  713. punpcklbw mm2, mm7
  715. punpckhbw mm1, mm7
  716. punpckhbw mm3, mm7
  718. paddw mm0, mm1
  719. paddw mm2, mm3
  721. paddw mm4, mm0
  722. paddw mm5, mm2
  724. add esi, ecx
  725. dec eax
  726. jnz .loop1
  728. paddusw mm4, mm5
  729. pmaddwd mm4, [mmx_one] ; merge sum
  730. movq mm5, mm4
  731. psrlq mm5, 32 
  732. paddd mm4, mm5
  734. psllq mm4, 32 ; blank upper dword
  735. psrlq mm4, 32 + 8 ; mm4 /= (16*16)
  737. punpckldq mm4, mm4
  738. packssdw mm4, mm4 ; mm4 = mean
  740. pxor mm6, mm6 ; mm6 = dev = 0
  741. mov eax, 16
  742. .loop2
  743. movq mm0, [edi]
  744. movq mm2, [edi + 8]
  746. movq mm1, mm0
  747. movq mm3, mm2
  749. punpcklbw mm0, mm7
  750. punpcklbw mm2, mm7
  752. punpckhbw mm1, mm7 ; mm01 = cur
  753. punpckhbw mm3, mm7 ; mm23 = cur2
  755. movq mm5, mm4 ;
  756. psubusw mm5, mm0 ;
  757. psubusw mm0, mm4 ;
  758. por mm0, mm5 ;
  759. movq mm5, mm4 ;
  760. psubusw mm5, mm1 ;
  761. psubusw mm1, mm4 ;
  762. por mm1, mm5 ; mm01 = |mm01 - mm4|
  765. movq mm5, mm4 ;
  766. psubusw mm5, mm2 ;
  767. psubusw mm2, mm4 ;
  768. por mm2, mm5 ;
  770. movq mm5, mm4 ;
  771. psubusw mm5, mm3 ;
  772. psubusw mm3, mm4 ;
  773. por mm3, mm5 ; mm23 = |mm23 - mm4|
  775. paddw mm0, mm1
  776. paddw mm2, mm3
  778. paddw mm6, mm0
  779. paddw mm6, mm2 ; dev += mm01 + mm23
  781. add edi, ecx
  782. dec eax
  783. jnz .loop2
  785. pmaddwd mm6, [mmx_one] ; merge dev
  786. movq mm7, mm6
  787. psrlq mm7, 32 
  788. paddd mm6, mm7
  789. movd eax, mm6
  791. pop  edi
  792. pop  esi
  794. ret
  798. ;===========================================================================
  799. ;
  800. ; uint32_t dev16_xmm(const uint8_t * const cur,
  801. ; const uint32_t stride);
  802. ;
  803. ; experimental!
  804. ;
  805. ;===========================================================================
  807. align 16
  808. cglobal dev16_xmm
  809. dev16_xmm
  811. push  esi
  812. push  edi
  814. pxor mm4, mm4 ; mm23 = sum = 0
  816. mov  esi, [esp + 8 + 4] ; cur
  817. mov  ecx, [esp + 8 + 8] ; stride
  818. mov     edi, esi
  820. ; mov eax, 16
  821. pxor mm7, mm7 ; mm7 = 0
  822. ;.loop1
  823. movq mm0, [esi]
  824. movq mm2, [esi + 8]
  826. psadbw mm0, mm7 ; abs(cur0 - 0) + abs(cur1 - 0) + ... + abs(cur7 - 0) -> mm0
  827. psadbw mm2, mm7 ; abs(cur8 - 0) + abs(cur9 - 0) + ... + abs(cur15 - 0) -> mm2
  829. paddw mm4,mm0 ; mean += mm0
  830. paddw mm4,mm2 ; mean += mm2
  832. add esi, ecx
  833. ; dec eax
  834. ; jnz .loop1
  836. movq mm0, [esi]
  837. movq mm2, [esi + 8]
  839. psadbw mm0, mm7
  840. psadbw mm2, mm7
  842. paddw mm4,mm0
  843. paddw mm4,mm2
  845. add esi, ecx
  847. movq mm0, [esi]
  848. movq mm2, [esi + 8]
  850. psadbw mm0, mm7
  851. psadbw mm2, mm7
  853. paddw mm4,mm0
  854. paddw mm4,mm2
  856. add esi, ecx
  858. movq mm0, [esi]
  859. movq mm2, [esi + 8]
  861. psadbw mm0, mm7
  862. psadbw mm2, mm7
  864. paddw mm4,mm0
  865. paddw mm4,mm2
  867. add esi, ecx
  869. movq mm0, [esi]
  870. movq mm2, [esi + 8]
  872. psadbw mm0, mm7
  873. psadbw mm2, mm7
  875. paddw mm4,mm0
  876. paddw mm4,mm2
  878. add esi, ecx
  880. movq mm0, [esi]
  881. movq mm2, [esi + 8]
  883. psadbw mm0, mm7
  884. psadbw mm2, mm7
  886. paddw mm4,mm0
  887. paddw mm4,mm2
  889. add esi, ecx
  891. movq mm0, [esi]
  892. movq mm2, [esi + 8]
  894. psadbw mm0, mm7
  895. psadbw mm2, mm7
  897. paddw mm4,mm0
  898. paddw mm4,mm2
  900. add esi, ecx
  902. movq mm0, [esi]
  903. movq mm2, [esi + 8]
  905. psadbw mm0, mm7
  906. psadbw mm2, mm7
  908. paddw mm4,mm0
  909. paddw mm4,mm2
  911. add esi, ecx
  913. movq mm0, [esi]
  914. movq mm2, [esi + 8]
  916. psadbw mm0, mm7
  917. psadbw mm2, mm7
  919. paddw mm4,mm0
  920. paddw mm4,mm2
  922. add esi, ecx
  924. movq mm0, [esi]
  925. movq mm2, [esi + 8]
  927. psadbw mm0, mm7
  928. psadbw mm2, mm7
  930. paddw mm4,mm0
  931. paddw mm4,mm2
  933. add esi, ecx
  935. movq mm0, [esi]
  936. movq mm2, [esi + 8]
  938. psadbw mm0, mm7
  939. psadbw mm2, mm7
  941. paddw mm4,mm0
  942. paddw mm4,mm2
  944. add esi, ecx
  946. movq mm0, [esi]
  947. movq mm2, [esi + 8]
  949. psadbw mm0, mm7
  950. psadbw mm2, mm7
  952. paddw mm4,mm0
  953. paddw mm4,mm2
  955. add esi, ecx
  957. movq mm0, [esi]
  958. movq mm2, [esi + 8]
  960. psadbw mm0, mm7
  961. psadbw mm2, mm7
  963. paddw mm4,mm0
  964. paddw mm4,mm2
  966. add esi, ecx
  968. movq mm0, [esi]
  969. movq mm2, [esi + 8]
  971. psadbw mm0, mm7
  972. psadbw mm2, mm7
  974. paddw mm4,mm0
  975. paddw mm4,mm2
  977. add esi, ecx
  979. movq mm0, [esi]
  980. movq mm2, [esi + 8]
  982. psadbw mm0, mm7
  983. psadbw mm2, mm7
  985. paddw mm4,mm0
  986. paddw mm4,mm2
  988. add esi, ecx
  990. movq mm0, [esi]
  991. movq mm2, [esi + 8]
  993. psadbw mm0, mm7
  994. psadbw mm2, mm7
  996. paddw mm4,mm0
  997. paddw mm4,mm2
  999. movq mm5, mm4
  1000. psllq mm5, 32 
  1001. paddd mm4, mm5
  1003. psrld mm4, 8
  1004. packssdw mm4, mm4
  1005. packuswb mm4, mm4
  1007. pxor mm6, mm6 ; mm6 = dev = 0
  1008. ; mov eax, 16
  1009. ;.loop2
  1010. movq mm0, [edi]
  1011. movq mm2, [edi + 8]
  1013. psadbw mm0, mm4 ; mm0 = |cur - mean|
  1014. psadbw mm2, mm4 ; mm0 = |cur2 - mean|
  1016. paddw mm6,mm0 ; dev += mm01
  1017. paddw mm6,mm2 ; dev += mm23
  1019. add edi, ecx
  1020. ; dec eax
  1021. ; jnz .loop2
  1023. movq mm0, [edi]
  1024. movq mm2, [edi + 8]
  1026. psadbw mm0, mm4
  1027. psadbw mm2, mm4
  1029. paddw mm6,mm0
  1030. paddw mm6,mm2
  1032. add edi, ecx
  1034. movq mm0, [edi]
  1035. movq mm2, [edi + 8]
  1037. psadbw mm0, mm4
  1038. psadbw mm2, mm4
  1040. paddw mm6,mm0
  1041. paddw mm6,mm2
  1043. add edi, ecx
  1045. movq mm0, [edi]
  1046. movq mm2, [edi + 8]
  1048. psadbw mm0, mm4
  1049. psadbw mm2, mm4
  1051. paddw mm6,mm0
  1052. paddw mm6,mm2
  1054. add edi, ecx
  1056. movq mm0, [edi]
  1057. movq mm2, [edi + 8]
  1059. psadbw mm0, mm4
  1060. psadbw mm2, mm4
  1062. paddw mm6,mm0
  1063. paddw mm6,mm2
  1065. add edi, ecx
  1067. movq mm0, [edi]
  1068. movq mm2, [edi + 8]
  1070. psadbw mm0, mm4
  1071. psadbw mm2, mm4
  1073. paddw mm6,mm0
  1074. paddw mm6,mm2
  1076. add edi, ecx
  1078. movq mm0, [edi]
  1079. movq mm2, [edi + 8]
  1081. psadbw mm0, mm4
  1082. psadbw mm2, mm4
  1084. paddw mm6,mm0
  1085. paddw mm6,mm2
  1087. add edi, ecx
  1089. movq mm0, [edi]
  1090. movq mm2, [edi + 8]
  1092. psadbw mm0, mm4
  1093. psadbw mm2, mm4
  1095. paddw mm6,mm0
  1096. paddw mm6,mm2
  1098. add edi, ecx
  1100. movq mm0, [edi]
  1101. movq mm2, [edi + 8]
  1103. psadbw mm0, mm4
  1104. psadbw mm2, mm4
  1106. paddw mm6,mm0
  1107. paddw mm6,mm2
  1109. add edi, ecx
  1111. movq mm0, [edi]
  1112. movq mm2, [edi + 8]
  1114. psadbw mm0, mm4
  1115. psadbw mm2, mm4
  1117. paddw mm6,mm0
  1118. paddw mm6,mm2
  1120. add edi, ecx
  1122. movq mm0, [edi]
  1123. movq mm2, [edi + 8]
  1125. psadbw mm0, mm4
  1126. psadbw mm2, mm4
  1128. paddw mm6,mm0
  1129. paddw mm6,mm2
  1131. add edi, ecx
  1133. movq mm0, [edi]
  1134. movq mm2, [edi + 8]
  1136. psadbw mm0, mm4
  1137. psadbw mm2, mm4
  1139. paddw mm6,mm0
  1140. paddw mm6,mm2
  1142. add edi, ecx
  1144. movq mm0, [edi]
  1145. movq mm2, [edi + 8]
  1147. psadbw mm0, mm4
  1148. psadbw mm2, mm4
  1150. paddw mm6,mm0
  1151. paddw mm6,mm2
  1153. add edi, ecx
  1155. movq mm0, [edi]
  1156. movq mm2, [edi + 8]
  1158. psadbw mm0, mm4
  1159. psadbw mm2, mm4
  1161. paddw mm6,mm0
  1162. paddw mm6,mm2
  1164. add edi, ecx
  1166. movq mm0, [edi]
  1167. movq mm2, [edi + 8]
  1169. psadbw mm0, mm4
  1170. psadbw mm2, mm4
  1172. paddw mm6,mm0
  1173. paddw mm6,mm2
  1175. add edi, ecx
  1177. movq mm0, [edi]
  1178. movq mm2, [edi + 8]
  1180. psadbw mm0, mm4
  1181. psadbw mm2, mm4
  1183. paddw mm6,mm0
  1184. paddw mm6,mm2
  1186. movq mm7, mm6
  1187. psllq mm7, 32 
  1188. paddd mm6, mm7
  1189. movd eax, mm6
  1191. pop  edi
  1192. pop  esi
  1194. ret