[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.5.2.4, Sun Oct 13 13:52:01 2002 UTC revision 1.14, Sun Oct 17 10:20:15 2004 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_  #ifndef _INTERPOLATE8X8_H_
27  #define _INTERPOLATE8X8_H_  #define _INTERPOLATE8X8_H_
28    
# Line 13  Line 38 
38                                                                          const uint8_t *src1,                                                                          const uint8_t *src1,
39                                                                          const uint8_t *src2,                                                                          const uint8_t *src2,
40                                                                          const uint32_t stride,                                                                          const uint32_t stride,
41                                                                          const uint32_t rounding);                                                                          const uint32_t rounding,
42                                                                            const uint32_t height);
43  typedef INTERPOLATE8X8_AVG2 *INTERPOLATE8X8_AVG2_PTR;  typedef INTERPOLATE8X8_AVG2 *INTERPOLATE8X8_AVG2_PTR;
44    
45  typedef void (INTERPOLATE8X8_AVG4) (uint8_t *dst,  typedef void (INTERPOLATE8X8_AVG4) (uint8_t *dst,
# Line 47  Line 73 
73    
74  typedef INTERPOLATE8X8_6TAP_LOWPASS *INTERPOLATE8X8_6TAP_LOWPASS_PTR;  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;  extern INTERPOLATE8X8_AVG2_PTR interpolate8x8_avg2;
93  extern INTERPOLATE8X8_AVG4_PTR interpolate8x8_avg4;  extern INTERPOLATE8X8_AVG4_PTR interpolate8x8_avg4;
94    
# Line 69  Line 107 
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;  INTERPOLATE8X8 interpolate8x8_halfpel_h_ia64;
145  INTERPOLATE8X8 interpolate8x8_halfpel_v_ia64;  INTERPOLATE8X8 interpolate8x8_halfpel_v_ia64;
146  INTERPOLATE8X8 interpolate8x8_halfpel_hv_ia64;  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  INTERPOLATE8X8_AVG2 interpolate8x8_avg2_c;  INTERPOLATE8X8_AVG2 interpolate8x8_avg2_c;
161  INTERPOLATE8X8_AVG4 interpolate8x8_avg4_c;  INTERPOLATE8X8_AVG4 interpolate8x8_avg4_c;
162    
163    #ifdef ARCH_IS_IA32
164  INTERPOLATE8X8_AVG2 interpolate8x8_avg2_mmx;  INTERPOLATE8X8_AVG2 interpolate8x8_avg2_mmx;
165  INTERPOLATE8X8_AVG4 interpolate8x8_avg4_mmx;  INTERPOLATE8X8_AVG4 interpolate8x8_avg4_mmx;
166    #endif
167    
168    #ifdef ARCH_IS_PPC
169    INTERPOLATE8X8_AVG2 interpolate8x8_avg2_altivec_c;
170    INTERPOLATE8X8_AVG4 interpolate8x8_avg4_altivec_c;
171    #endif
172    
173  INTERPOLATE_LOWPASS interpolate8x8_lowpass_h_c;  INTERPOLATE_LOWPASS interpolate8x8_lowpass_h_c;
174  INTERPOLATE_LOWPASS interpolate8x8_lowpass_v_c;  INTERPOLATE_LOWPASS interpolate8x8_lowpass_v_c;
# Line 104  Line 182 
182  INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_h_c;  INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_h_c;
183  INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_v_c;  INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_v_c;
184    
185    #ifdef ARCH_IS_IA32
186  INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_h_mmx;  INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_h_mmx;
187  INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_v_mmx;  INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_v_mmx;
188    #endif
189    
190    #ifdef ARCH_IS_PPC
191    INTERPOLATE8X8_6TAP_LOWPASS interpolate8x8_6tap_lowpass_h_altivec_c;
192    #endif
193    
194  static __inline void  static __inline void
195  interpolate8x8_switch(uint8_t * const cur,  interpolate8x8_switch(uint8_t * const cur,
# Line 117  Line 201 
201                                            const uint32_t stride,                                            const uint32_t stride,
202                                            const uint32_t rounding)                                            const uint32_t rounding)
203  {  {
         int32_t ddx, ddy;  
204    
205          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));
206          {          uint8_t * const dst = cur + (int)(y * stride + x);
207    
208            switch (((dx & 1) << 1) + (dy & 1))     { /* ((dx%2)?2:0)+((dy%2)?1:0) */
209          case 0:          case 0:
210                  ddx = dx / 2;                  transfer8x8_copy(dst, src, stride);
                 ddy = dy / 2;  
                 transfer8x8_copy(cur + y * stride + x,  
                                                  refn + (int)((y + ddy) * stride + x + ddx), stride);  
211                  break;                  break;
   
212          case 1:          case 1:
213                  ddx = dx / 2;                  interpolate8x8_halfpel_v(dst, src, stride, rounding);
                 ddy = (dy - 1) / 2;  
                 interpolate8x8_halfpel_v(cur + y * stride + x,  
                                                                  refn + (int)((y + ddy) * stride + x + ddx), stride,  
                                                                  rounding);  
214                  break;                  break;
   
215          case 2:          case 2:
216                  ddx = (dx - 1) / 2;                  interpolate8x8_halfpel_h(dst, src, stride, rounding);
                 ddy = dy / 2;  
                 interpolate8x8_halfpel_h(cur + y * stride + x,  
                                                                  refn + (int)((y + ddy) * stride + x + ddx), stride,  
                                                                  rounding);  
217                  break;                  break;
   
218          default:          default:
219                  ddx = (dx - 1) / 2;                  interpolate8x8_halfpel_hv(dst, src, stride, rounding);
                 ddy = (dy - 1) / 2;  
                 interpolate8x8_halfpel_hv(cur + y * stride + x,  
                                                                  refn + (int)((y + ddy) * stride + x + ddx), stride,  
                                                                   rounding);  
220                  break;                  break;
221          }          }
222  }  }
223    
224  static __inline uint8_t *  static __inline void
225  interpolate8x8_switch2(uint8_t * const buffer,  interpolate8x8_add_switch(uint8_t * const cur,
226                                            const uint8_t * const refn,                                            const uint8_t * const refn,
227                                            const uint32_t x,                                            const uint32_t x,
228                                            const uint32_t y,                                            const uint32_t y,
# Line 164  Line 231 
231                                            const uint32_t stride,                                            const uint32_t stride,
232                                            const uint32_t rounding)                                            const uint32_t rounding)
233  {  {
         int32_t ddx, ddy;  
234    
235          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));
236          {          uint8_t * const dst = cur + (int)(y * stride + x);
         case 0:  
                 return (uint8_t *)refn + (int)((y + dy/2) * stride + x + dx/2);  
237    
238            switch (((dx & 1) << 1) + (dy & 1))     { /* ((dx%2)?2:0)+((dy%2)?1:0) */
239            case 0:
240                    interpolate8x8_halfpel_add(dst, src, stride, rounding);
241                    break;
242          case 1:          case 1:
243                  ddx = dx / 2;                  interpolate8x8_halfpel_v_add(dst, src, stride, rounding);
                 ddy = (dy - 1) / 2;  
                 interpolate8x8_halfpel_v(buffer,  
                                                                  refn + (int)((y + ddy) * stride + x + ddx), stride,  
                                                                  rounding);  
244                  break;                  break;
   
245          case 2:          case 2:
246                  ddx = (dx - 1) / 2;                  interpolate8x8_halfpel_h_add(dst, src, stride, rounding);
                 ddy = dy / 2;  
                 interpolate8x8_halfpel_h(buffer,  
                                                                  refn + (int)((y + ddy) * stride + x + ddx), stride,  
                                                                  rounding);  
247                  break;                  break;
   
248          default:          default:
249                  ddx = (dx - 1) / 2;                  interpolate8x8_halfpel_hv_add(dst, src, stride, rounding);
                 ddy = (dy - 1) / 2;  
                 interpolate8x8_halfpel_hv(buffer,  
                                                                  refn + (int)((y + ddy) * stride + x + ddx), stride,  
                                                                   rounding);  
250                  break;                  break;
251          }          }
         return buffer;  
252  }  }
253    
254    static __inline void
255  static __inline void interpolate8x8_quarterpel(uint8_t * const cur,  interpolate16x16_switch(uint8_t * const cur,
256                                       uint8_t * const refn,                                            const uint8_t * const refn,
257                                           uint8_t * const refh,                                            const uint32_t x,
258                                           uint8_t * const refv,                                            const uint32_t y,
259                                           uint8_t * const refhv,                                            const int32_t dx,
260                                       const uint32_t x, const uint32_t y,                                            const int dy,
                                          const int32_t dx,  const int dy,  
261                                           const uint32_t stride,                                           const uint32_t stride,
262                                           const uint32_t rounding)                                           const uint32_t rounding)
263  {  {
264          const int32_t xRef = x*4 + dx;          interpolate8x8_switch(cur, refn, x,   y,   dx, dy, stride, rounding);
265          const int32_t yRef = y*4 + dy;          interpolate8x8_switch(cur, refn, x+8, y,   dx, dy, stride, rounding);
266            interpolate8x8_switch(cur, refn, x,   y+8, dx, dy, stride, rounding);
267          uint8_t *src, *dst;          interpolate8x8_switch(cur, refn, x+8, y+8, dx, dy, stride, rounding);
         uint8_t *halfpel_h, *halfpel_v, *halfpel_hv;  
         int32_t x_int, y_int, x_frac, y_frac;  
   
         x_int = xRef/4;  
         if (xRef < 0 && xRef % 4)  
                 x_int--;  
   
         x_frac = xRef - (4*x_int);  
   
         y_int  = yRef/4;  
         if (yRef < 0 && yRef % 4)  
                 y_int--;  
   
         y_frac = yRef - (4*y_int);  
   
         src = refn + y_int * stride + x_int;  
         halfpel_h = refh;  
         halfpel_v = refv;  
         halfpel_hv = refhv;  
   
         dst = cur + y * stride + x;  
   
         switch((y_frac << 2) | (x_frac)) {  
   
         case 0:  
                 transfer8x8_copy(dst, src, stride);  
                 break;  
   
         case 1:  
                 interpolate8x8_lowpass_h(halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg2(dst, src, halfpel_h, stride, rounding);  
                 break;  
   
         case 2:  
             interpolate8x8_lowpass_h(dst, src, stride, rounding);  
                 break;  
   
         case 3:  
                 interpolate8x8_lowpass_h(halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg2(dst, src + 1, halfpel_h, stride, rounding);  
                 break;  
   
         case 4:  
                 interpolate8x8_lowpass_v(halfpel_v, src, stride, rounding);  
                 interpolate8x8_avg2(dst, src, halfpel_v, stride, rounding);  
                 break;  
   
         case 5:  
                 interpolate8x8_lowpass_v(halfpel_v, src, stride, rounding);  
                 interpolate8x8_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg4(dst, src, halfpel_h, halfpel_v, halfpel_hv, stride, rounding);  
                 break;  
   
         case 6:  
                 interpolate8x8_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg2(dst, halfpel_h, halfpel_hv, stride, rounding);  
                 break;  
   
         case 7:  
                 interpolate8x8_lowpass_v(halfpel_v, src+1, stride, rounding);  
                 interpolate8x8_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg4(dst, src+1, halfpel_h, halfpel_v, halfpel_hv, stride, rounding);  
                 break;  
   
         case 8:  
             interpolate8x8_lowpass_v(dst, src, stride, rounding);  
                 break;  
   
         case 9:  
                 interpolate8x8_lowpass_v(halfpel_v, src, stride, rounding);  
                 interpolate8x8_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg2(dst, halfpel_v, halfpel_hv, stride, rounding);  
                 break;  
   
         case 10:  
                 interpolate8x8_lowpass_hv(dst, halfpel_h, src, stride, rounding);  
                 break;  
   
         case 11:  
                 interpolate8x8_lowpass_v(halfpel_v, src+1, stride, rounding);  
                 interpolate8x8_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg2(dst, halfpel_v, halfpel_hv, stride, rounding);  
                 break;  
   
         case 12:  
                 interpolate8x8_lowpass_v(halfpel_v, src, stride, rounding);  
                 interpolate8x8_avg2(dst, src+stride, halfpel_v, stride, rounding);  
                 break;  
   
         case 13:  
                 interpolate8x8_lowpass_v(halfpel_v, src, stride, rounding);  
                 interpolate8x8_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg4(dst, src+stride, halfpel_h+stride, halfpel_v, halfpel_hv, stride, rounding);  
                 break;  
   
         case 14:  
                 interpolate8x8_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg2(dst, halfpel_h+stride, halfpel_hv, stride, rounding);  
                 break;  
   
         case 15:  
                 interpolate8x8_lowpass_v(halfpel_v, src+1, stride, rounding);  
                 interpolate8x8_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg4(dst, src+stride+1, halfpel_h+stride, halfpel_v, halfpel_hv, stride, rounding);  
                 break;  
         }  
268  }  }
269    
270    static __inline void
271  static __inline void interpolate16x16_quarterpel(uint8_t * const cur,  interpolate16x16_add_switch(uint8_t * const cur,
272                                       uint8_t * const refn,                                            const uint8_t * const refn,
273                                           uint8_t * const refh,                                            const uint32_t x,
274                                           uint8_t * const refv,                                            const uint32_t y,
275                                           uint8_t * const refhv,                                            const int32_t dx,
276                                       const uint32_t x, const uint32_t y,                                            const int dy,
                                          const int32_t dx,  const int dy,  
277                                           const uint32_t stride,                                           const uint32_t stride,
278                                           const uint32_t rounding)                                           const uint32_t rounding)
279  {  {
280          const int32_t xRef = x*4 + dx;          interpolate8x8_add_switch(cur, refn, x,   y,   dx, dy, stride, rounding);
281          const int32_t yRef = y*4 + dy;          interpolate8x8_add_switch(cur, refn, x+8, y,   dx, dy, stride, rounding);
282            interpolate8x8_add_switch(cur, refn, x,   y+8, dx, dy, stride, rounding);
283          uint8_t *src, *dst;          interpolate8x8_add_switch(cur, refn, x+8, y+8, dx, dy, stride, rounding);
284          uint8_t *halfpel_h, *halfpel_v, *halfpel_hv;  }
         int32_t x_int, y_int, x_frac, y_frac;  
   
         x_int = xRef/4;  
         if (xRef < 0 && xRef % 4)  
                 x_int--;  
   
         x_frac = xRef - (4*x_int);  
   
         y_int  = yRef/4;  
         if (yRef < 0 && yRef % 4)  
                 y_int--;  
285    
286          y_frac = yRef - (4*y_int);  static __inline void
287    interpolate32x32_switch(uint8_t * const cur,
288                                              const uint8_t * const refn,
289                                              const uint32_t x,
290                                              const uint32_t y,
291                                              const int32_t dx,
292                                              const int dy,
293                                              const uint32_t stride,
294                                              const uint32_t rounding)
295    {
296            interpolate16x16_switch(cur, refn, x,    y,    dx, dy, stride, rounding);
297            interpolate16x16_switch(cur, refn, x+16, y,    dx, dy, stride, rounding);
298            interpolate16x16_switch(cur, refn, x,    y+16, dx, dy, stride, rounding);
299            interpolate16x16_switch(cur, refn, x+16, y+16, dx, dy, stride, rounding);
300    }
301    
302          src = refn + y_int * stride + x_int;  static __inline void
303          halfpel_h = refh;  interpolate32x32_add_switch(uint8_t * const cur,
304          halfpel_v = refv;                                            const uint8_t * const refn,
305          halfpel_hv = refhv;                                            const uint32_t x,
306                                              const uint32_t y,
307                                              const int32_t dx,
308                                              const int dy,
309                                              const uint32_t stride,
310                                              const uint32_t rounding)
311    {
312            interpolate16x16_add_switch(cur, refn, x,    y,    dx, dy, stride, rounding);
313            interpolate16x16_add_switch(cur, refn, x+16, y,    dx, dy, stride, rounding);
314            interpolate16x16_add_switch(cur, refn, x,    y+16, dx, dy, stride, rounding);
315            interpolate16x16_add_switch(cur, refn, x+16, y+16, dx, dy, stride, rounding);
316    }
317    
318          dst = cur + y * stride + x;  static __inline uint8_t *
319    interpolate8x8_switch2(uint8_t * const buffer,
320                                              const uint8_t * const refn,
321                                              const int x,
322                                              const int y,
323                                              const int dx,
324                                              const int dy,
325                                              const uint32_t stride,
326                                              const uint32_t rounding)
327    {
328    
329          switch((y_frac << 2) | (x_frac)) {          const uint8_t * const src = refn + (int)((y + (dy>>1)) * stride + x + (dx>>1));
330    
331            switch (((dx & 1) << 1) + (dy & 1))     { /* ((dx%2)?2:0)+((dy%2)?1:0) */
332          case 0:          case 0:
333                  transfer8x8_copy(dst, src, stride);                  return (uint8_t *)src;
                 transfer8x8_copy(dst+8, src+8, stride);  
                 transfer8x8_copy(dst+8*stride, src+8*stride, stride);  
                 transfer8x8_copy(dst+8*stride+8, src+8*stride+8, stride);  
                 break;  
   
334          case 1:          case 1:
335                  interpolate16x16_lowpass_h(halfpel_h, src, stride, rounding);                  interpolate8x8_halfpel_v(buffer, src, stride, rounding);
                 interpolate8x8_avg2(dst, src, halfpel_h, stride, rounding);  
                 interpolate8x8_avg2(dst+8, src+8, halfpel_h+8, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride, src+8*stride, halfpel_h+8*stride, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride+8, src+8*stride+8, halfpel_h+8*stride+8, stride, rounding);  
336                  break;                  break;
   
337          case 2:          case 2:
338              interpolate16x16_lowpass_h(dst, src, stride, rounding);                  interpolate8x8_halfpel_h(buffer, src, stride, rounding);
                 break;  
   
         case 3:  
                 interpolate16x16_lowpass_h(halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg2(dst, src + 1, halfpel_h, stride, rounding);  
                 interpolate8x8_avg2(dst+8, src + 8 + 1, halfpel_h+8, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride, src + 8*stride + 1, halfpel_h+8*stride, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride+8, src+8*stride+8 + 1, halfpel_h+8*stride+8, stride, rounding);  
                 break;  
   
         case 4:  
                 interpolate16x16_lowpass_v(halfpel_v, src, stride, rounding);  
                 interpolate8x8_avg2(dst, src, halfpel_v, stride, rounding);  
                 interpolate8x8_avg2(dst+8, src+8, halfpel_v+8, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride, src+8*stride, halfpel_v+8*stride, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride+8, src+8*stride+8, halfpel_v+8*stride+8, stride, rounding);  
                 break;  
   
         case 5:  
                 interpolate16x16_lowpass_v(halfpel_v, src, stride, rounding);  
                 interpolate16x16_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg4(dst, src, halfpel_h, halfpel_v, halfpel_hv, stride, rounding);  
                 interpolate8x8_avg4(dst+8, src+8, halfpel_h+8, halfpel_v+8, halfpel_hv+8, stride, rounding);  
                 interpolate8x8_avg4(dst+8*stride, src+8*stride, halfpel_h+8*stride, halfpel_v+8*stride, halfpel_hv+8*stride, stride, rounding);  
                 interpolate8x8_avg4(dst+8*stride+8, src+8*stride+8, halfpel_h+8*stride+8, halfpel_v+8*stride+8, halfpel_hv+8*stride+8, stride, rounding);  
                 break;  
   
         case 6:  
                 interpolate16x16_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg2(dst, halfpel_h, halfpel_hv, stride, rounding);  
                 interpolate8x8_avg2(dst+8, halfpel_h+8, halfpel_hv+8, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride, halfpel_h+8*stride, halfpel_hv+8*stride, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride+8, halfpel_h+8*stride+8, halfpel_hv+8*stride+8, stride, rounding);  
                 break;  
   
         case 7:  
                 interpolate16x16_lowpass_v(halfpel_v, src+1, stride, rounding);  
                 interpolate16x16_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg4(dst, src+1, halfpel_h, halfpel_v, halfpel_hv, stride, rounding);  
                 interpolate8x8_avg4(dst+8, src+8+1, halfpel_h+8, halfpel_v+8, halfpel_hv+8, stride, rounding);  
                 interpolate8x8_avg4(dst+8*stride, src+1+8*stride, halfpel_h+8*stride, halfpel_v+8*stride, halfpel_hv+8*stride, stride, rounding);  
                 interpolate8x8_avg4(dst+8*stride+8, src+1+8*stride+8, halfpel_h+8*stride+8, halfpel_v+8*stride+8, halfpel_hv+8*stride+8, stride, rounding);  
                 break;  
   
         case 8:  
             interpolate16x16_lowpass_v(dst, src, stride, rounding);  
                 break;  
   
         case 9:  
                 interpolate16x16_lowpass_v(halfpel_v, src, stride, rounding);  
                 interpolate16x16_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg2(dst, halfpel_v, halfpel_hv, stride, rounding);  
                 interpolate8x8_avg2(dst+8, halfpel_v+8, halfpel_hv+8, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride, halfpel_v+8*stride, halfpel_hv+8*stride, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride+8, halfpel_v+8*stride+8, halfpel_hv+8*stride+8, stride, rounding);  
                 break;  
   
         case 10:  
                 interpolate16x16_lowpass_hv(dst, halfpel_h, src, stride, rounding);  
                 break;  
   
         case 11:  
                 interpolate16x16_lowpass_v(halfpel_v, src+1, stride, rounding);  
                 interpolate16x16_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg2(dst, halfpel_v, halfpel_hv, stride, rounding);  
                 interpolate8x8_avg2(dst+8, halfpel_v+8, halfpel_hv+8, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride, halfpel_v+8*stride, halfpel_hv+8*stride, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride+8, halfpel_v+8*stride+8, halfpel_hv+8*stride+8, stride, rounding);  
                 break;  
   
         case 12:  
                 interpolate16x16_lowpass_v(halfpel_v, src, stride, rounding);  
                 interpolate8x8_avg2(dst, src+stride, halfpel_v, stride, rounding);  
                 interpolate8x8_avg2(dst+8, src+stride+8, halfpel_v+8, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride, src+stride+8*stride, halfpel_v+8*stride, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride+8, src+stride+8*stride+8, halfpel_v+8*stride+8, stride, rounding);  
                 break;  
   
         case 13:  
                 interpolate16x16_lowpass_v(halfpel_v, src, stride, rounding);  
                 interpolate16x16_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg4(dst, src+stride, halfpel_h+stride, halfpel_v, halfpel_hv, stride, rounding);  
                 interpolate8x8_avg4(dst+8, src+stride+8, halfpel_h+stride+8, halfpel_v+8, halfpel_hv+8, stride, rounding);  
                 interpolate8x8_avg4(dst+8*stride, src+stride+8*stride, halfpel_h+stride+8*stride, halfpel_v+8*stride, halfpel_hv+8*stride, stride, rounding);  
                 interpolate8x8_avg4(dst+8*stride+8, src+stride+8*stride+8, halfpel_h+stride+8*stride+8, halfpel_v+8*stride+8, halfpel_hv+8*stride+8, stride, rounding);  
                 break;  
   
         case 14:  
                 interpolate16x16_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg2(dst, halfpel_h+stride, halfpel_hv, stride, rounding);  
                 interpolate8x8_avg2(dst+8, halfpel_h+stride+8, halfpel_hv+8, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride, halfpel_h+stride+8*stride, halfpel_hv+8*stride, stride, rounding);  
                 interpolate8x8_avg2(dst+8*stride+8, halfpel_h+stride+8*stride+8, halfpel_hv+8*stride+8, stride, rounding);  
                 break;  
   
         case 15:  
                 interpolate16x16_lowpass_v(halfpel_v, src+1, stride, rounding);  
                 interpolate16x16_lowpass_hv(halfpel_hv, halfpel_h, src, stride, rounding);  
                 interpolate8x8_avg4(dst, src+stride+1, halfpel_h+stride, halfpel_v, halfpel_hv, stride, rounding);  
                 interpolate8x8_avg4(dst+8, src+stride+1+8, halfpel_h+stride+8, halfpel_v+8, halfpel_hv+8, stride, rounding);  
                 interpolate8x8_avg4(dst+8*stride, src+stride+1+8*stride, halfpel_h+stride+8*stride, halfpel_v+8*stride, halfpel_hv+8*stride, stride, rounding);  
                 interpolate8x8_avg4(dst+8*stride+8, src+stride+1+8*stride+8, halfpel_h+stride+8*stride+8, halfpel_v+8*stride+8, halfpel_hv+8*stride+8, stride, rounding);  
339                  break;                  break;
340          default:          default:
341                  interpolate8x8_quarterpel(cur, refn, refh, refv, refhv, x, y, dx, dy, stride, rounding);                  interpolate8x8_halfpel_hv(buffer, src, stride, rounding);
                 interpolate8x8_quarterpel(cur, refn, refh, refv, refhv, x+8, y, dx, dy, stride, rounding);  
                 interpolate8x8_quarterpel(cur, refn, refh, refv, refhv, x, y+8, dx, dy, stride, rounding);  
                 interpolate8x8_quarterpel(cur, refn, refh, refv, refhv, x+8, y+8, dx, dy, stride, rounding);  
                 return;  
342                  break;                  break;
343          }          }
344            return buffer;
345  }  }
346    
347  #endif  #endif

Legend:
Removed from v.1.5.2.4  
changed lines
  Added in v.1.14

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