[cvs] / xvidcore / src / image / interpolate8x8.h Repository:
ViewVC logotype

Diff of /xvidcore/src/image/interpolate8x8.h

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

revision 1.3, Wed Jun 12 20:38:40 2002 UTC revision 1.15, Wed Jan 5 23:02:15 2005 UTC
# Line 1  Line 1 
1    /*****************************************************************************
2     *
3     *  XVID MPEG-4 VIDEO CODEC
4     *  - Interpolation related header  -
5     *
6     *  Copyright(C) 2001-2003 Peter Ross <pross@xvid.org>
7     *
8     *  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
10     *  the Free Software Foundation ; either version 2 of the License, or
11     *  (at your option) any later version.
12     *
13     *  This program is distributed in the hope that it will be useful,
14     *  but WITHOUT ANY WARRANTY ; without even the implied warranty of
15     *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
16     *  GNU General Public License for more details.
17     *
18     *  You should have received a copy of the GNU General Public License
19     *  along with this program ; if not, write to the Free Software
20     *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
21     *
22     * $Id$
23     *
24     ****************************************************************************/
25    
26    #ifndef _INTERPOLATE8X8_H_
27    #define _INTERPOLATE8X8_H_
28    
29  #include "../utils/mem_transfer.h"  #include "../utils/mem_transfer.h"
30    
31  typedef void (INTERPOLATE8X8) (uint8_t * const dst,  typedef void (INTERPOLATE8X8) (uint8_t * const dst,
# Line 6  Line 34 
34                                                             const uint32_t rounding);                                                             const uint32_t rounding);
35  typedef INTERPOLATE8X8 *INTERPOLATE8X8_PTR;  typedef INTERPOLATE8X8 *INTERPOLATE8X8_PTR;
36    
37    typedef void (INTERPOLATE8X8_AVG2) (uint8_t *dst,
38                                                                            const uint8_t *src1,
39                                                                            const uint8_t *src2,
40                                                                            const uint32_t stride,
41                                                                            const uint32_t rounding,
42                                                                            const uint32_t height);
43    typedef INTERPOLATE8X8_AVG2 *INTERPOLATE8X8_AVG2_PTR;
44    
45    typedef void (INTERPOLATE8X8_AVG4) (uint8_t *dst,
46                                                                            const uint8_t *src1,
47                                                                            const uint8_t *src2,
48                                                                            const uint8_t *src3,
49                                                                            const uint8_t *src4,
50                                                                            const uint32_t stride,
51                                                                            const uint32_t rounding);
52    typedef INTERPOLATE8X8_AVG4 *INTERPOLATE8X8_AVG4_PTR;
53    
54    typedef void (INTERPOLATE_LOWPASS) (uint8_t *dst,
55                                                                               uint8_t *src,
56                                                                               int32_t stride,
57                                                                               int32_t rounding);
58    
59    typedef INTERPOLATE_LOWPASS *INTERPOLATE_LOWPASS_PTR;
60    
61    typedef void (INTERPOLATE_LOWPASS_HV) (uint8_t *dst1,
62                                                                                      uint8_t *dst2,
63                                                                                      uint8_t *src,
64                                                                                      int32_t stride,
65                                                                                      int32_t rounding);
66    
67    typedef INTERPOLATE_LOWPASS_HV *INTERPOLATE_LOWPASS_HV_PTR;
68    
69    typedef void (INTERPOLATE8X8_6TAP_LOWPASS) (uint8_t *dst,
70                                                                                    uint8_t *src,
71                                                                                    int32_t stride,
72                                                                                    int32_t rounding);
73    
74    typedef INTERPOLATE8X8_6TAP_LOWPASS *INTERPOLATE8X8_6TAP_LOWPASS_PTR;
75    
76    /* These function do: dst = interpolate(src) */
77  extern INTERPOLATE8X8_PTR interpolate8x8_halfpel_h;  extern INTERPOLATE8X8_PTR interpolate8x8_halfpel_h;
78  extern INTERPOLATE8X8_PTR interpolate8x8_halfpel_v;  extern INTERPOLATE8X8_PTR interpolate8x8_halfpel_v;
79  extern INTERPOLATE8X8_PTR interpolate8x8_halfpel_hv;  extern INTERPOLATE8X8_PTR interpolate8x8_halfpel_hv;
80    
81    /* These functions do: dst = (dst+interpolate(src) + 1)/2
82     * Suitable for direct/interpolated bvop prediction block
83     * building w/o the need for intermediate interpolated result
84     * storing/reading
85     * NB: the rounding applies to the interpolation, but not
86     *     the averaging step which will always use rounding=0 */
87    extern INTERPOLATE8X8_PTR interpolate8x8_halfpel_add;
88    extern INTERPOLATE8X8_PTR interpolate8x8_halfpel_h_add;
89    extern INTERPOLATE8X8_PTR interpolate8x8_halfpel_v_add;
90    extern INTERPOLATE8X8_PTR interpolate8x8_halfpel_hv_add;
91    
92    extern INTERPOLATE8X8_AVG2_PTR interpolate8x8_avg2;
93    extern INTERPOLATE8X8_AVG4_PTR interpolate8x8_avg4;
94    
95    extern INTERPOLATE_LOWPASS_PTR interpolate8x8_lowpass_h;
96    extern INTERPOLATE_LOWPASS_PTR interpolate8x8_lowpass_v;
97    
98    extern INTERPOLATE_LOWPASS_PTR interpolate16x16_lowpass_h;
99    extern INTERPOLATE_LOWPASS_PTR interpolate16x16_lowpass_v;
100    
101    extern INTERPOLATE_LOWPASS_HV_PTR interpolate8x8_lowpass_hv;
102    extern INTERPOLATE_LOWPASS_HV_PTR interpolate16x16_lowpass_hv;
103    
104    extern INTERPOLATE8X8_6TAP_LOWPASS_PTR interpolate8x8_6tap_lowpass_h;
105    extern INTERPOLATE8X8_6TAP_LOWPASS_PTR interpolate8x8_6tap_lowpass_v;
106    
107  INTERPOLATE8X8 interpolate8x8_halfpel_h_c;  INTERPOLATE8X8 interpolate8x8_halfpel_h_c;
108  INTERPOLATE8X8 interpolate8x8_halfpel_v_c;  INTERPOLATE8X8 interpolate8x8_halfpel_v_c;
109  INTERPOLATE8X8 interpolate8x8_halfpel_hv_c;  INTERPOLATE8X8 interpolate8x8_halfpel_hv_c;
110    INTERPOLATE8X8 interpolate8x8_halfpel_add_c;
111    INTERPOLATE8X8 interpolate8x8_halfpel_h_add_c;
112    INTERPOLATE8X8 interpolate8x8_halfpel_v_add_c;
113    INTERPOLATE8X8 interpolate8x8_halfpel_hv_add_c;
114    
115    #ifdef ARCH_IS_IA32
116  INTERPOLATE8X8 interpolate8x8_halfpel_h_mmx;  INTERPOLATE8X8 interpolate8x8_halfpel_h_mmx;
117  INTERPOLATE8X8 interpolate8x8_halfpel_v_mmx;  INTERPOLATE8X8 interpolate8x8_halfpel_v_mmx;
118  INTERPOLATE8X8 interpolate8x8_halfpel_hv_mmx;  INTERPOLATE8X8 interpolate8x8_halfpel_hv_mmx;
119    
120    INTERPOLATE8X8 interpolate8x8_halfpel_add_mmx;
121    INTERPOLATE8X8 interpolate8x8_halfpel_h_add_mmx;
122    INTERPOLATE8X8 interpolate8x8_halfpel_v_add_mmx;
123    INTERPOLATE8X8 interpolate8x8_halfpel_hv_add_mmx;
124    
125  INTERPOLATE8X8 interpolate8x8_halfpel_h_xmm;  INTERPOLATE8X8 interpolate8x8_halfpel_h_xmm;
126  INTERPOLATE8X8 interpolate8x8_halfpel_v_xmm;  INTERPOLATE8X8 interpolate8x8_halfpel_v_xmm;
127  INTERPOLATE8X8 interpolate8x8_halfpel_hv_xmm;  INTERPOLATE8X8 interpolate8x8_halfpel_hv_xmm;
128    
129    INTERPOLATE8X8 interpolate8x8_halfpel_add_xmm;
130    INTERPOLATE8X8 interpolate8x8_halfpel_h_add_xmm;
131    INTERPOLATE8X8 interpolate8x8_halfpel_v_add_xmm;
132    INTERPOLATE8X8 interpolate8x8_halfpel_hv_add_xmm;
133    
134  INTERPOLATE8X8 interpolate8x8_halfpel_h_3dn;  INTERPOLATE8X8 interpolate8x8_halfpel_h_3dn;
135  INTERPOLATE8X8 interpolate8x8_halfpel_v_3dn;  INTERPOLATE8X8 interpolate8x8_halfpel_v_3dn;
136  INTERPOLATE8X8 interpolate8x8_halfpel_hv_3dn;  INTERPOLATE8X8 interpolate8x8_halfpel_hv_3dn;
137    
138    INTERPOLATE8X8 interpolate8x8_halfpel_h_3dne;
139    INTERPOLATE8X8 interpolate8x8_halfpel_v_3dne;
140    INTERPOLATE8X8 interpolate8x8_halfpel_hv_3dne;
141    #endif
142    
143    #ifdef ARCH_IS_IA64
144    INTERPOLATE8X8 interpolate8x8_halfpel_h_ia64;
145    INTERPOLATE8X8 interpolate8x8_halfpel_v_ia64;
146    INTERPOLATE8X8 interpolate8x8_halfpel_hv_ia64;
147    #endif
148    
149    #ifdef ARCH_IS_PPC
150    INTERPOLATE8X8 interpolate8x8_halfpel_h_altivec_c;
151    INTERPOLATE8X8 interpolate8x8_halfpel_v_altivec_c;
152    INTERPOLATE8X8 interpolate8x8_halfpel_hv_altivec_c;
153    
154    INTERPOLATE8X8 interpolate8x8_halfpel_add_altivec_c;
155    INTERPOLATE8X8 interpolate8x8_halfpel_h_add_altivec_c;
156    INTERPOLATE8X8 interpolate8x8_halfpel_v_add_altivec_c;
157    INTERPOLATE8X8 interpolate8x8_halfpel_hv_add_altivec_c;
158    #endif
159    
160    #ifdef ARCH_IS_X86_64
161    INTERPOLATE8X8 interpolate8x8_halfpel_h_x86_64;
162    INTERPOLATE8X8 interpolate8x8_halfpel_v_x86_64;
163    INTERPOLATE8X8 interpolate8x8_halfpel_hv_x86_64;
164    
165    INTERPOLATE8X8 interpolate8x8_halfpel_add_x86_64;
166    INTERPOLATE8X8 interpolate8x8_halfpel_h_add_x86_64;
167    INTERPOLATE8X8 interpolate8x8_halfpel_v_add_x86_64;
168    INTERPOLATE8X8 interpolate8x8_halfpel_hv_add_x86_64;
169    #endif
170    
171    INTERPOLATE8X8_AVG2 interpolate8x8_avg2_c;
172    INTERPOLATE8X8_AVG4 interpolate8x8_avg4_c;
173    
174    #ifdef ARCH_IS_IA32
175    INTERPOLATE8X8_AVG2 interpolate8x8_avg2_mmx;
176    INTERPOLATE8X8_AVG4 interpolate8x8_avg4_mmx;
177    #endif
178    
179    #ifdef ARCH_IS_PPC
180    INTERPOLATE8X8_AVG2 interpolate8x8_avg2_altivec_c;
181    INTERPOLATE8X8_AVG4 interpolate8x8_avg4_altivec_c;
182    #endif
183    
184    #ifdef ARCH_IS_X86_64
185    INTERPOLATE8X8_AVG2 interpolate8x8_avg2_x86_64;
186    INTERPOLATE8X8_AVG4 interpolate8x8_avg4_x86_64;
187    #endif
188    
189    INTERPOLATE_LOWPASS interpolate8x8_lowpass_h_c;
190    INTERPOLATE_LOWPASS interpolate8x8_lowpass_v_c;
191    
192    INTERPOLATE_LOWPASS interpolate16x16_lowpass_h_c;
193    INTERPOLATE_LOWPASS interpolate16x16_lowpass_v_c;
194    
195    INTERPOLATE_LOWPASS_HV interpolate8x8_lowpass_hv_c;
196    INTERPOLATE_LOWPASS_HV interpolate16x16_lowpass_hv_c;
197    
198    INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_h_c;
199    INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_v_c;
200    
201    #ifdef ARCH_IS_IA32
202    INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_h_mmx;
203    INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_v_mmx;
204    #endif
205    
206    #ifdef ARCH_IS_PPC
207    INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_h_altivec_c;
208    #endif
209    
210    #ifdef ARCH_IS_X86_64
211    INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_h_x86_64;
212    INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_v_x86_64;
213    #endif
214    
215  static __inline void  static __inline void
216  interpolate8x8_switch(uint8_t * const cur,  interpolate8x8_switch(uint8_t * const cur,
217                                            const uint8_t * const refn,                                            const uint8_t * const refn,
# Line 36  Line 222 
222                                            const uint32_t stride,                                            const uint32_t stride,
223                                            const uint32_t rounding)                                            const uint32_t rounding)
224  {  {
         int32_t ddx, ddy;  
225    
226          switch (((dx & 1) << 1) + (dy & 1))     // ((dx%2)?2:0)+((dy%2)?1:0)          const uint8_t * const src = refn + (int)((y + (dy>>1)) * stride + x + (dx>>1));
227          {          uint8_t * const dst = cur + (int)(y * stride + x);
228    
229            switch (((dx & 1) << 1) + (dy & 1))     { /* ((dx%2)?2:0)+((dy%2)?1:0) */
230          case 0:          case 0:
231                  ddx = dx / 2;                  transfer8x8_copy(dst, src, stride);
                 ddy = dy / 2;  
                 transfer8x8_copy(cur + y * stride + x,  
                                                  refn + (y + ddy) * stride + x + ddx, stride);  
232                  break;                  break;
   
233          case 1:          case 1:
234                  ddx = dx / 2;                  interpolate8x8_halfpel_v(dst, src, stride, rounding);
                 ddy = (dy - 1) / 2;  
                 interpolate8x8_halfpel_v(cur + y * stride + x,  
                                                                  refn + (y + ddy) * stride + x + ddx, stride,  
                                                                  rounding);  
235                  break;                  break;
   
236          case 2:          case 2:
237                  ddx = (dx - 1) / 2;                  interpolate8x8_halfpel_h(dst, src, stride, rounding);
                 ddy = dy / 2;  
                 interpolate8x8_halfpel_h(cur + y * stride + x,  
                                                                  refn + (y + ddy) * stride + x + ddx, stride,  
                                                                  rounding);  
238                  break;                  break;
239            default:
240                    interpolate8x8_halfpel_hv(dst, src, stride, rounding);
241                    break;
242            }
243    }
244    
245    static __inline void
246    interpolate8x8_add_switch(uint8_t * const cur,
247                                              const uint8_t * const refn,
248                                              const uint32_t x,
249                                              const uint32_t y,
250                                              const int32_t dx,
251                                              const int dy,
252                                              const uint32_t stride,
253                                              const uint32_t rounding)
254    {
255    
256            const uint8_t * const src = refn + (int)((y + (dy>>1)) * stride + x + (dx>>1));
257            uint8_t * const dst = cur + (int)(y * stride + x);
258    
259            switch (((dx & 1) << 1) + (dy & 1))     { /* ((dx%2)?2:0)+((dy%2)?1:0) */
260            case 0:
261                    interpolate8x8_halfpel_add(dst, src, stride, rounding);
262                    break;
263            case 1:
264                    interpolate8x8_halfpel_v_add(dst, src, stride, rounding);
265                    break;
266            case 2:
267                    interpolate8x8_halfpel_h_add(dst, src, stride, rounding);
268                    break;
269          default:          default:
270                  ddx = (dx - 1) / 2;                  interpolate8x8_halfpel_hv_add(dst, src, stride, rounding);
                 ddy = (dy - 1) / 2;  
                 interpolate8x8_halfpel_hv(cur + y * stride + x,  
                                                                   refn + (y + ddy) * stride + x + ddx, stride,  
                                                                   rounding);  
271                  break;                  break;
272          }          }
273  }  }
274    
275    static __inline void
276    interpolate16x16_switch(uint8_t * const cur,
277                                              const uint8_t * const refn,
278                                              const uint32_t x,
279                                              const uint32_t y,
280                                              const int32_t dx,
281                                              const int dy,
282                                              const uint32_t stride,
283                                              const uint32_t rounding)
284    {
285            interpolate8x8_switch(cur, refn, x,   y,   dx, dy, stride, rounding);
286            interpolate8x8_switch(cur, refn, x+8, y,   dx, dy, stride, rounding);
287            interpolate8x8_switch(cur, refn, x,   y+8, dx, dy, stride, rounding);
288            interpolate8x8_switch(cur, refn, x+8, y+8, dx, dy, stride, rounding);
289    }
290    
291  void interpolate8x8_c(uint8_t * const dst,  static __inline void
292                                            const uint8_t * const src,  interpolate16x16_add_switch(uint8_t * const cur,
293                                              const uint8_t * const refn,
294                                              const uint32_t x,
295                                              const uint32_t y,
296                                              const int32_t dx,
297                                              const int dy,
298                                              const uint32_t stride,
299                                              const uint32_t rounding)
300    {
301            interpolate8x8_add_switch(cur, refn, x,   y,   dx, dy, stride, rounding);
302            interpolate8x8_add_switch(cur, refn, x+8, y,   dx, dy, stride, rounding);
303            interpolate8x8_add_switch(cur, refn, x,   y+8, dx, dy, stride, rounding);
304            interpolate8x8_add_switch(cur, refn, x+8, y+8, dx, dy, stride, rounding);
305    }
306    
307    static __inline void
308    interpolate32x32_switch(uint8_t * const cur,
309                                              const uint8_t * const refn,
310                                            const uint32_t x,                                            const uint32_t x,
311                                            const uint32_t y,                                            const uint32_t y,
312                                            const uint32_t stride);                                            const int32_t dx,
313                                              const int dy,
314                                              const uint32_t stride,
315                                              const uint32_t rounding)
316    {
317            interpolate16x16_switch(cur, refn, x,    y,    dx, dy, stride, rounding);
318            interpolate16x16_switch(cur, refn, x+16, y,    dx, dy, stride, rounding);
319            interpolate16x16_switch(cur, refn, x,    y+16, dx, dy, stride, rounding);
320            interpolate16x16_switch(cur, refn, x+16, y+16, dx, dy, stride, rounding);
321    }
322    
323    static __inline void
324    interpolate32x32_add_switch(uint8_t * const cur,
325                                              const uint8_t * const refn,
326                                              const uint32_t x,
327                                              const uint32_t y,
328                                              const int32_t dx,
329                                              const int dy,
330                                              const uint32_t stride,
331                                              const uint32_t rounding)
332    {
333            interpolate16x16_add_switch(cur, refn, x,    y,    dx, dy, stride, rounding);
334            interpolate16x16_add_switch(cur, refn, x+16, y,    dx, dy, stride, rounding);
335            interpolate16x16_add_switch(cur, refn, x,    y+16, dx, dy, stride, rounding);
336            interpolate16x16_add_switch(cur, refn, x+16, y+16, dx, dy, stride, rounding);
337    }
338    
339    static __inline uint8_t *
340    interpolate8x8_switch2(uint8_t * const buffer,
341                                              const uint8_t * const refn,
342                                              const int x,
343                                              const int y,
344                                              const int dx,
345                                              const int dy,
346                                              const uint32_t stride,
347                                              const uint32_t rounding)
348    {
349    
350            const uint8_t * const src = refn + (int)((y + (dy>>1)) * stride + x + (dx>>1));
351    
352            switch (((dx & 1) << 1) + (dy & 1))     { /* ((dx%2)?2:0)+((dy%2)?1:0) */
353            case 0:
354                    return (uint8_t *)src;
355            case 1:
356                    interpolate8x8_halfpel_v(buffer, src, stride, rounding);
357                    break;
358            case 2:
359                    interpolate8x8_halfpel_h(buffer, src, stride, rounding);
360                    break;
361            default:
362                    interpolate8x8_halfpel_hv(buffer, src, stride, rounding);
363                    break;
364            }
365            return buffer;
366    }
367    
368    #endif

Legend:
Removed from v.1.3  
changed lines
  Added in v.1.15

No admin address has been configured
ViewVC Help
Powered by ViewVC 1.0.4