[svn] / trunk / xvidcore / src / image / x86_asm / interpolate8x8_3dne.asm Repository:
ViewVC logotype

Diff of /trunk/xvidcore/src/image/x86_asm/interpolate8x8_3dne.asm

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 851, Sat Feb 15 15:22:19 2003 UTC revision 1793, Tue Nov 11 20:46:24 2008 UTC
# Line 1  Line 1 
1  ;/**************************************************************************  ;/*****************************************************************************
2  ; *  ; *
3  ; *     XVID MPEG-4 VIDEO CODEC  ; *     XVID MPEG-4 VIDEO CODEC
4  ; *     xmm 8x8 block-based halfpel interpolation  ; *  - 3dne pipeline optimized  8x8 block-based halfpel interpolation -
5  ; *  ; *
6  ; *     This program is an implementation of a part of one or more MPEG-4  ; *  Copyright(C) 2002 Jaan Kalda
 ; *     Video tools as specified in ISO/IEC 14496-2 standard.  Those intending  
 ; *     to use this software module in hardware or software products are  
 ; *     advised that its use may infringe existing patents or copyrights, and  
 ; *     any such use would be at such party's own risk.  The original  
 ; *     developer of this software module and his/her company, and subsequent  
 ; *     editors and their companies, will have no liability for use of this  
 ; *     software or modifications or derivatives thereof.  
7  ; *  ; *
8  ; *     This program is free software; you can redistribute it and/or modify  ; *     This program is free software; you can redistribute it and/or modify
9  ; *     it under the terms of the GNU General Public License as published by  ; *     it under the terms of the GNU General Public License as published by
# Line 24  Line 17 
17  ; *  ; *
18  ; *     You should have received a copy of the GNU General Public License  ; *     You should have received a copy of the GNU General Public License
19  ; *     along with this program; if not, write to the Free Software  ; *     along with this program; if not, write to the Free Software
20  ; *     Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.  ; *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
21  ; *  ; *
22  ; *************************************************************************/  ; ****************************************************************************/
23    
24  ; these 3dne functions are compatible with iSSE, but are optimized specifically for  ; these 3dne functions are compatible with iSSE, but are optimized specifically
25  ; K7 pipelines  ; for K7 pipelines
 ;  
 ;------------------------------------------------------------------------------  
 ; 09.12.2002  Athlon optimizations contributed by Jaan Kalda  
 ;------------------------------------------------------------------------------  
26    
27  bits 32  BITS 32
28    
29  %macro cglobal 1  %macro cglobal 1
30          %ifdef PREFIX          %ifdef PREFIX
31                    %ifdef MARK_FUNCS
32                            global _%1:function %1.endfunc-%1
33                            %define %1 _%1:function %1.endfunc-%1
34                            %define ENDFUNC .endfunc
35                    %else
36                  global _%1                  global _%1
37                  %define %1 _%1                  %define %1 _%1
38                            %define ENDFUNC
39                    %endif
40            %else
41                    %ifdef MARK_FUNCS
42                            global %1:function %1.endfunc-%1
43                            %define ENDFUNC .endfunc
44          %else          %else
45                  global %1                  global %1
46                            %define ENDFUNC
47                    %endif
48          %endif          %endif
49  %endmacro  %endmacro
50  %macro nop4 0  
51  DB 08Dh,074h,026h,0  ;=============================================================================
52  %endmacro  ; Read only data
53    ;=============================================================================
54    
55  %ifdef FORMAT_COFF  %ifdef FORMAT_COFF
56  section .data data  SECTION .rodata
57  %else  %else
58  section .data data align=16  SECTION .rodata align=16
59  %endif  %endif
60    
61    ALIGN 16
62  align 16  mmx_one:
 mmx_one  
63  times 8 db 1  times 8 db 1
64    
65  align 8  ALIGN 8
66  mm_minusone:  mm_minusone:
67  dd -1,-1  dd -1,-1
68    
69  section .text  ;=============================================================================
70    ; Macros
71    ;=============================================================================
72    
73    %macro nop4 0
74    DB 08Dh,074h,026h,0
75    %endmacro
76    
77    ;=============================================================================
78    ; Macros
79    ;=============================================================================
80    
81    SECTION .text
82    
83  cglobal interpolate8x8_halfpel_h_3dne  cglobal interpolate8x8_halfpel_h_3dne
84  cglobal interpolate8x8_halfpel_v_3dne  cglobal interpolate8x8_halfpel_v_3dne
85  cglobal interpolate8x8_halfpel_hv_3dne  cglobal interpolate8x8_halfpel_hv_3dne
86    
87  ;===========================================================================  cglobal interpolate8x4_halfpel_h_3dne
88    cglobal interpolate8x4_halfpel_v_3dne
89    cglobal interpolate8x4_halfpel_hv_3dne
90    
91    ;-----------------------------------------------------------------------------
92  ;  ;
93  ; void interpolate8x8_halfpel_h_3dne(uint8_t * const dst,  ; void interpolate8x8_halfpel_h_3dne(uint8_t * const dst,
94  ;                                               const uint8_t * const src,  ;                                               const uint8_t * const src,
95  ;                                               const uint32_t stride,  ;                                               const uint32_t stride,
96  ;                                               const uint32_t rounding);  ;                                               const uint32_t rounding);
97  ;  ;
98  ;===========================================================================  ;-----------------------------------------------------------------------------
99    
100  %macro COPY_H_SSE_RND0 1  %macro COPY_H_SSE_RND0 1
101  %if (%1)  %if (%1)
102    movq mm0,  [eax]    movq mm0,  [eax]
103  %else  %else
104    movq mm0,  [dword eax]    movq mm0, [eax+0]
105      ; ---
106      ; nasm >0.99.x rejects the original statement:
107      ;   movq mm0, [dword eax]
108      ; as it is ambiguous. for this statement nasm <0.99.x would
109      ; generate "movq mm0,[eax+0]"
110      ; ---
111  %endif  %endif
112    pavgb mm0, [eax+1]    pavgb mm0, [eax+1]
113    movq mm1,  [eax+edx]    movq mm1,  [eax+edx]
# Line 112  Line 137 
137          movq [ecx+edx], mm1          movq [ecx+edx], mm1
138  %endmacro  %endmacro
139    
140  align 16  ALIGN 16
141  interpolate8x8_halfpel_h_3dne:  interpolate8x8_halfpel_h_3dne:
142    
143    mov eax, [esp+ 8] ; Src    mov eax, [esp+ 8] ; Src
# Line 131  Line 156 
156    COPY_H_SSE_RND0 1    COPY_H_SSE_RND0 1
157    ret    ret
158    
159  .rounding1  .rounding1:
160     ; we use: (i+j)/2 = ( i+j+1 )/2 - (i^j)&1     ; we use: (i+j)/2 = ( i+j+1 )/2 - (i^j)&1
161    mov ecx, [esp+ 4] ; Dst    mov ecx, [esp+ 4] ; Dst
162    movq mm7, [mmx_one]    movq mm7, [mmx_one]
# Line 143  Line 168 
168    lea ecx,[ecx+2*edx]    lea ecx,[ecx+2*edx]
169    COPY_H_SSE_RND1    COPY_H_SSE_RND1
170    ret    ret
171    ENDFUNC
172    
173  ;===========================================================================  ;-----------------------------------------------------------------------------
174  ;  ;
175  ; void interpolate8x8_halfpel_v_3dne(uint8_t * const dst,  ; void interpolate8x8_halfpel_v_3dne(uint8_t * const dst,
176  ;                                               const uint8_t * const src,  ;                                               const uint8_t * const src,
177  ;                                               const uint32_t stride,  ;                                               const uint32_t stride,
178  ;                                               const uint32_t rounding);  ;                                               const uint32_t rounding);
179  ;  ;
180  ;===========================================================================  ;-----------------------------------------------------------------------------
181    
182  align 16  ALIGN 16
183  interpolate8x8_halfpel_v_3dne:  interpolate8x8_halfpel_v_3dne:
184    
185    mov eax, [esp+ 8] ; Src    mov eax, [esp+ 8] ; Src
# Line 201  Line 227 
227    movq [ecx+edx],mm7    movq [ecx+edx],mm7
228    ret    ret
229    
230  align 8  ALIGN 8
231  .rounding1  .rounding1:
232    pcmpeqb mm0,mm0    pcmpeqb mm0,mm0
233    psubusb mm0,[eax]    psubusb mm0,[eax]
234    add eax,edx    add eax,edx
# Line 266  Line 292 
292    movq [ecx], mm4    movq [ecx], mm4
293    movq [ecx+edx], mm5    movq [ecx+edx], mm5
294    ret    ret
295  ;===========================================================================  ENDFUNC
296    
297    ;-----------------------------------------------------------------------------
298  ;  ;
299  ; void interpolate8x8_halfpel_hv_3dne(uint8_t * const dst,  ; void interpolate8x8_halfpel_hv_3dne(uint8_t * const dst,
300  ;                                               const uint8_t * const src,  ;                                               const uint8_t * const src,
# Line 274  Line 302 
302  ;                                               const uint32_t rounding);  ;                                               const uint32_t rounding);
303  ;  ;
304  ;  ;
305  ;===========================================================================  ;-----------------------------------------------------------------------------
306    
307  ; The trick is to correct the result of 'pavgb' with some combination of the  ; The trick is to correct the result of 'pavgb' with some combination of the
308  ; lsb's of the 4 input values i,j,k,l, and their intermediate 'pavgb' (s and t).  ; lsb's of the 4 input values i,j,k,l, and their intermediate 'pavgb' (s and t).
# Line 360  Line 388 
388      movq [ecx+edx], mm0      movq [ecx+edx], mm0
389  %endmacro  %endmacro
390    
391  align 16  ALIGN 16
392  interpolate8x8_halfpel_hv_3dne:  interpolate8x8_halfpel_hv_3dne:
393    mov eax, [esp+ 8] ; Src    mov eax, [esp+ 8] ; Src
394    mov edx, [esp+12] ; stride    mov edx, [esp+12] ; stride
# Line 386  Line 414 
414    COPY_HV_SSE_RND0    COPY_HV_SSE_RND0
415    ret    ret
416    
417  align 16  ALIGN 16
418  .rounding1  .rounding1:
419    COPY_HV_SSE_RND1    COPY_HV_SSE_RND1
420      lea ecx,[ecx+2*edx]      lea ecx,[ecx+2*edx]
421    COPY_HV_SSE_RND1    COPY_HV_SSE_RND1
# Line 396  Line 424 
424      lea ecx,[ecx+2*edx]      lea ecx,[ecx+2*edx]
425    COPY_HV_SSE_RND1    COPY_HV_SSE_RND1
426    ret    ret
427    ENDFUNC
428    
429    ;-----------------------------------------------------------------------------
430    ;
431    ; void interpolate8x4_halfpel_h_3dne(uint8_t * const dst,
432    ;                       const uint8_t * const src,
433    ;                       const uint32_t stride,
434    ;                       const uint32_t rounding);
435    ;
436    ;-----------------------------------------------------------------------------
437    
438    ALIGN 16
439    interpolate8x4_halfpel_h_3dne:
440    
441      mov eax, [esp+ 8] ; Src
442      mov edx, [esp+12] ; stride
443      dec dword [esp+16]; rounding
444    
445      jz .rounding1
446      mov ecx, [esp+ 4] ; Dst
447    
448      COPY_H_SSE_RND0 0
449      lea ecx,[ecx+2*edx]
450      COPY_H_SSE_RND0 1
451      ret
452    
453    .rounding1:
454     ; we use: (i+j)/2 = ( i+j+1 )/2 - (i^j)&1
455      mov ecx, [esp+ 4] ; Dst
456      movq mm7, [mmx_one]
457      COPY_H_SSE_RND1
458      lea ecx, [ecx+2*edx]
459      COPY_H_SSE_RND1
460      ret
461    ENDFUNC
462    
463    ;-----------------------------------------------------------------------------
464    ;
465    ; void interpolate8x4_halfpel_v_3dne(uint8_t * const dst,
466    ;                       const uint8_t * const src,
467    ;                       const uint32_t stride,
468    ;                       const uint32_t rounding);
469    ;
470    ;-----------------------------------------------------------------------------
471    
472    ALIGN 16
473    interpolate8x4_halfpel_v_3dne:
474    
475      mov eax, [esp+ 8] ; Src
476      mov edx, [esp+12] ; stride
477      dec dword [esp+16]; rounding
478    
479        ; we process 2 line at a time
480    
481      jz .rounding1
482      pxor mm2,mm2
483      movq mm0, [eax]
484      movq mm1, [eax+edx]
485      por mm2, [eax+2*edx]      ; Something like preload (pipelining)
486      mov ecx, [esp+ 4] ; Dst
487      lea eax, [eax+2*edx]
488      pxor mm4, mm4
489      pavgb mm0, mm1
490      pavgb mm1, mm2
491      movq [byte ecx], mm0
492      movq [ecx+edx], mm1
493    
494      pxor mm6, mm6
495      add eax, edx
496      lea ecx, [ecx+2*edx]
497      movq mm3, [byte eax]
498      por mm4, [eax+edx]
499      lea eax, [eax+2*edx]
500      pavgb mm2, mm3
501      pavgb mm3, mm4
502      movq [ecx], mm2
503      movq [ecx+edx], mm3
504    
505      ret
506    
507    ALIGN 8
508    .rounding1:
509      pcmpeqb mm0, mm0
510      psubusb mm0, [eax]            ; eax==line0
511      add eax, edx                  ; eax==line1
512      mov ecx, [esp+ 4] ; Dst
513    
514      push esi
515    
516      pcmpeqb mm1, mm1
517      pcmpeqb mm2, mm2
518      mov esi, mm_minusone
519      psubusb mm1, [byte eax]       ; line1
520      psubusb mm2, [eax+edx]        ; line2
521      lea eax, [eax+2*edx]          ; eax==line3
522      movq mm6, [esi]
523      movq mm7, [esi]
524      pavgb mm0, mm1
525      pavgb mm1, mm2
526      psubusb mm6, mm0
527      psubusb mm7, mm1
528      movq [ecx], mm6               ; store line0
529      movq [ecx+edx], mm7           ; store line1
530    
531      lea ecx, [ecx+2*edx]
532      pcmpeqb mm3, mm3
533      pcmpeqb mm4, mm4
534      psubusb mm3, [eax]            ; line3
535      psubusb mm4, [eax+edx]        ; line4
536      lea eax, [eax+2*edx]          ; eax==line 5
537      pavgb mm2, mm3
538      pavgb mm3, mm4
539      movq mm0, [esi]
540      movq mm1, [esi]
541      psubusb mm0, mm2
542      psubusb mm1, mm3
543      movq [ecx], mm0
544      movq [ecx+edx], mm1
545    
546      pop esi
547    
548      ret
549    
550    ENDFUNC
551    
552    ;-----------------------------------------------------------------------------
553    ;
554    ; void interpolate8x4_halfpel_hv_3dne(uint8_t * const dst,
555    ;                       const uint8_t * const src,
556    ;                       const uint32_t stride,
557    ;                       const uint32_t rounding);
558    ;
559    ;
560    ;-----------------------------------------------------------------------------
561    
562    ALIGN 16
563    interpolate8x4_halfpel_hv_3dne:
564      mov eax, [esp+ 8]     ; Src
565      mov edx, [esp+12]     ; stride
566      dec dword [esp+16]    ; rounding
567    
568        ; loop invariants: mm2=(i+j+1)/2  and  mm3= i^j
569      movq mm2, [eax]
570      movq mm3, [eax+1]
571      movq mm6, mm2
572      pavgb mm2, mm3
573      pxor mm3, mm6         ; mm2/mm3 ready
574      mov ecx, [esp+ 4]     ; Dst
575      movq mm7, [mmx_one]
576    
577      jz near .rounding1
578      lea ebp,[byte ebp]
579      COPY_HV_SSE_RND0
580      lea ecx,[ecx+2*edx]
581      COPY_HV_SSE_RND0
582      ret
583    
584    ALIGN 16
585    .rounding1:
586      COPY_HV_SSE_RND1
587      lea ecx,[ecx+2*edx]
588      COPY_HV_SSE_RND1
589      ret
590    ENDFUNC
591    
592    
593    %ifidn __OUTPUT_FORMAT__,elf
594    section ".note.GNU-stack" noalloc noexec nowrite progbits
595    %endif
596    
Legend:
Removed from v.851  
changed lines
  Added in v.1793
No admin address has been configured
 ViewVC Help
Powered by ViewVC 1.0.4