Some deblocking-related optimizations
[x262.git] / common / macroblock.c
1 /*****************************************************************************
2  * macroblock.c: h264 encoder library
3  *****************************************************************************
4  * Copyright (C) 2003-2008 x264 project
5  *
6  * Authors: Jason Garrett-Glaser <darkshikari@gmail.com>
7  *          Laurent Aimar <fenrir@via.ecp.fr>
8  *          Loren Merritt <lorenm@u.washington.edu>
9  *
10  * This program is free software; you can redistribute it and/or modify
11  * it under the terms of the GNU General Public License as published by
12  * the Free Software Foundation; either version 2 of the License, or
13  * (at your option) any later version.
14  *
15  * This program is distributed in the hope that it will be useful,
16  * but WITHOUT ANY WARRANTY; without even the implied warranty of
17  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
18  * GNU General Public License for more details.
19  *
20  * You should have received a copy of the GNU General Public License
21  * along with this program; if not, write to the Free Software
22  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02111, USA.
23  *****************************************************************************/
24
25 #include "common.h"
26 #include "encoder/me.h"
27
28 static NOINLINE void x264_mb_mc_0xywh( x264_t *h, int x, int y, int width, int height )
29 {
30     int i8    = x264_scan8[0]+x+8*y;
31     int i_ref = h->mb.cache.ref[0][i8];
32     int mvx   = x264_clip3( h->mb.cache.mv[0][i8][0], h->mb.mv_min[0], h->mb.mv_max[0] ) + 4*4*x;
33     int mvy   = x264_clip3( h->mb.cache.mv[0][i8][1], h->mb.mv_min[1], h->mb.mv_max[1] ) + 4*4*y;
34
35     h->mc.mc_luma( &h->mb.pic.p_fdec[0][4*y*FDEC_STRIDE+4*x], FDEC_STRIDE,
36                    h->mb.pic.p_fref[0][i_ref], h->mb.pic.i_stride[0],
37                    mvx, mvy, 4*width, 4*height, &h->sh.weight[i_ref][0] );
38
39     // chroma is offset if MCing from a field of opposite parity
40     if( h->mb.b_interlaced & i_ref )
41         mvy += (h->mb.i_mb_y & 1)*4 - 2;
42
43     h->mc.mc_chroma( &h->mb.pic.p_fdec[1][2*y*FDEC_STRIDE+2*x], FDEC_STRIDE,
44                      h->mb.pic.p_fref[0][i_ref][4], h->mb.pic.i_stride[1],
45                      mvx, mvy, 2*width, 2*height );
46
47     if( h->sh.weight[i_ref][1].weightfn )
48         h->sh.weight[i_ref][1].weightfn[width>>1]( &h->mb.pic.p_fdec[1][2*y*FDEC_STRIDE+2*x], FDEC_STRIDE,
49                                                    &h->mb.pic.p_fdec[1][2*y*FDEC_STRIDE+2*x], FDEC_STRIDE,
50                                                    &h->sh.weight[i_ref][1], height*2 );
51
52     h->mc.mc_chroma( &h->mb.pic.p_fdec[2][2*y*FDEC_STRIDE+2*x], FDEC_STRIDE,
53                      h->mb.pic.p_fref[0][i_ref][5], h->mb.pic.i_stride[2],
54                      mvx, mvy, 2*width, 2*height );
55
56     if( h->sh.weight[i_ref][2].weightfn )
57         h->sh.weight[i_ref][2].weightfn[width>>1]( &h->mb.pic.p_fdec[2][2*y*FDEC_STRIDE+2*x], FDEC_STRIDE,
58                                                    &h->mb.pic.p_fdec[2][2*y*FDEC_STRIDE+2*x], FDEC_STRIDE,
59                                                    &h->sh.weight[i_ref][2],height*2 );
60
61 }
62 static NOINLINE void x264_mb_mc_1xywh( x264_t *h, int x, int y, int width, int height )
63 {
64     int i8    = x264_scan8[0]+x+8*y;
65     int i_ref = h->mb.cache.ref[1][i8];
66     int mvx   = x264_clip3( h->mb.cache.mv[1][i8][0], h->mb.mv_min[0], h->mb.mv_max[0] ) + 4*4*x;
67     int mvy   = x264_clip3( h->mb.cache.mv[1][i8][1], h->mb.mv_min[1], h->mb.mv_max[1] ) + 4*4*y;
68
69     h->mc.mc_luma( &h->mb.pic.p_fdec[0][4*y*FDEC_STRIDE+4*x], FDEC_STRIDE,
70                    h->mb.pic.p_fref[1][i_ref], h->mb.pic.i_stride[0],
71                    mvx, mvy, 4*width, 4*height, weight_none );
72
73     if( h->mb.b_interlaced & i_ref )
74         mvy += (h->mb.i_mb_y & 1)*4 - 2;
75
76     h->mc.mc_chroma( &h->mb.pic.p_fdec[1][2*y*FDEC_STRIDE+2*x], FDEC_STRIDE,
77                      h->mb.pic.p_fref[1][i_ref][4], h->mb.pic.i_stride[1],
78                      mvx, mvy, 2*width, 2*height );
79
80     h->mc.mc_chroma( &h->mb.pic.p_fdec[2][2*y*FDEC_STRIDE+2*x], FDEC_STRIDE,
81                      h->mb.pic.p_fref[1][i_ref][5], h->mb.pic.i_stride[2],
82                      mvx, mvy, 2*width, 2*height );
83 }
84
85 static NOINLINE void x264_mb_mc_01xywh( x264_t *h, int x, int y, int width, int height )
86 {
87     int i8 = x264_scan8[0]+x+8*y;
88     int i_ref0 = h->mb.cache.ref[0][i8];
89     int i_ref1 = h->mb.cache.ref[1][i8];
90     int weight = h->mb.bipred_weight[i_ref0][i_ref1];
91     int mvx0   = x264_clip3( h->mb.cache.mv[0][i8][0], h->mb.mv_min[0], h->mb.mv_max[0] ) + 4*4*x;
92     int mvx1   = x264_clip3( h->mb.cache.mv[1][i8][0], h->mb.mv_min[0], h->mb.mv_max[0] ) + 4*4*x;
93     int mvy0   = x264_clip3( h->mb.cache.mv[0][i8][1], h->mb.mv_min[1], h->mb.mv_max[1] ) + 4*4*y;
94     int mvy1   = x264_clip3( h->mb.cache.mv[1][i8][1], h->mb.mv_min[1], h->mb.mv_max[1] ) + 4*4*y;
95     int i_mode = x264_size2pixel[height][width];
96     int i_stride0 = 16, i_stride1 = 16;
97     ALIGNED_ARRAY_16( uint8_t, tmp0,[16*16] );
98     ALIGNED_ARRAY_16( uint8_t, tmp1,[16*16] );
99     uint8_t *src0, *src1;
100
101     src0 = h->mc.get_ref( tmp0, &i_stride0, h->mb.pic.p_fref[0][i_ref0], h->mb.pic.i_stride[0],
102                           mvx0, mvy0, 4*width, 4*height, weight_none );
103     src1 = h->mc.get_ref( tmp1, &i_stride1, h->mb.pic.p_fref[1][i_ref1], h->mb.pic.i_stride[0],
104                           mvx1, mvy1, 4*width, 4*height, weight_none );
105     h->mc.avg[i_mode]( &h->mb.pic.p_fdec[0][4*y*FDEC_STRIDE+4*x], FDEC_STRIDE,
106                        src0, i_stride0, src1, i_stride1, weight );
107
108     if( h->mb.b_interlaced & i_ref0 )
109         mvy0 += (h->mb.i_mb_y & 1)*4 - 2;
110     if( h->mb.b_interlaced & i_ref1 )
111         mvy1 += (h->mb.i_mb_y & 1)*4 - 2;
112
113     h->mc.mc_chroma( tmp0, 16, h->mb.pic.p_fref[0][i_ref0][4], h->mb.pic.i_stride[1],
114                      mvx0, mvy0, 2*width, 2*height );
115     h->mc.mc_chroma( tmp1, 16, h->mb.pic.p_fref[1][i_ref1][4], h->mb.pic.i_stride[1],
116                      mvx1, mvy1, 2*width, 2*height );
117     h->mc.avg[i_mode+3]( &h->mb.pic.p_fdec[1][2*y*FDEC_STRIDE+2*x], FDEC_STRIDE, tmp0, 16, tmp1, 16, weight );
118     h->mc.mc_chroma( tmp0, 16, h->mb.pic.p_fref[0][i_ref0][5], h->mb.pic.i_stride[2],
119                      mvx0, mvy0, 2*width, 2*height );
120     h->mc.mc_chroma( tmp1, 16, h->mb.pic.p_fref[1][i_ref1][5], h->mb.pic.i_stride[2],
121                      mvx1, mvy1, 2*width, 2*height );
122     h->mc.avg[i_mode+3]( &h->mb.pic.p_fdec[2][2*y*FDEC_STRIDE+2*x], FDEC_STRIDE, tmp0, 16, tmp1, 16, weight );
123 }
124
125 void x264_mb_mc_8x8( x264_t *h, int i8 )
126 {
127     int x = 2*(i8&1);
128     int y = 2*(i8>>1);
129
130     if( h->sh.i_type == SLICE_TYPE_P )
131     {
132         switch( h->mb.i_sub_partition[i8] )
133         {
134             case D_L0_8x8:
135                 x264_mb_mc_0xywh( h, x, y, 2, 2 );
136                 break;
137             case D_L0_8x4:
138                 x264_mb_mc_0xywh( h, x, y+0, 2, 1 );
139                 x264_mb_mc_0xywh( h, x, y+1, 2, 1 );
140                 break;
141             case D_L0_4x8:
142                 x264_mb_mc_0xywh( h, x+0, y, 1, 2 );
143                 x264_mb_mc_0xywh( h, x+1, y, 1, 2 );
144                 break;
145             case D_L0_4x4:
146                 x264_mb_mc_0xywh( h, x+0, y+0, 1, 1 );
147                 x264_mb_mc_0xywh( h, x+1, y+0, 1, 1 );
148                 x264_mb_mc_0xywh( h, x+0, y+1, 1, 1 );
149                 x264_mb_mc_0xywh( h, x+1, y+1, 1, 1 );
150                 break;
151         }
152     }
153     else
154     {
155         int scan8 = x264_scan8[0] + x + 8*y;
156
157         if( h->mb.cache.ref[0][scan8] >= 0 )
158             if( h->mb.cache.ref[1][scan8] >= 0 )
159                 x264_mb_mc_01xywh( h, x, y, 2, 2 );
160             else
161                 x264_mb_mc_0xywh( h, x, y, 2, 2 );
162         else
163             x264_mb_mc_1xywh( h, x, y, 2, 2 );
164     }
165 }
166
167 void x264_mb_mc( x264_t *h )
168 {
169     if( h->mb.i_partition == D_8x8 )
170     {
171         for( int i = 0; i < 4; i++ )
172             x264_mb_mc_8x8( h, i );
173     }
174     else
175     {
176         int ref0a = h->mb.cache.ref[0][x264_scan8[ 0]];
177         int ref0b = h->mb.cache.ref[0][x264_scan8[12]];
178         int ref1a = h->mb.cache.ref[1][x264_scan8[ 0]];
179         int ref1b = h->mb.cache.ref[1][x264_scan8[12]];
180
181         if( h->mb.i_partition == D_16x16 )
182         {
183             if( ref0a >= 0 )
184                 if( ref1a >= 0 ) x264_mb_mc_01xywh( h, 0, 0, 4, 4 );
185                 else             x264_mb_mc_0xywh ( h, 0, 0, 4, 4 );
186             else                 x264_mb_mc_1xywh ( h, 0, 0, 4, 4 );
187         }
188         else if( h->mb.i_partition == D_16x8 )
189         {
190             if( ref0a >= 0 )
191                 if( ref1a >= 0 ) x264_mb_mc_01xywh( h, 0, 0, 4, 2 );
192                 else             x264_mb_mc_0xywh ( h, 0, 0, 4, 2 );
193             else                 x264_mb_mc_1xywh ( h, 0, 0, 4, 2 );
194
195             if( ref0b >= 0 )
196                 if( ref1b >= 0 ) x264_mb_mc_01xywh( h, 0, 2, 4, 2 );
197                 else             x264_mb_mc_0xywh ( h, 0, 2, 4, 2 );
198             else                 x264_mb_mc_1xywh ( h, 0, 2, 4, 2 );
199         }
200         else if( h->mb.i_partition == D_8x16 )
201         {
202             if( ref0a >= 0 )
203                 if( ref1a >= 0 ) x264_mb_mc_01xywh( h, 0, 0, 2, 4 );
204                 else             x264_mb_mc_0xywh ( h, 0, 0, 2, 4 );
205             else                 x264_mb_mc_1xywh ( h, 0, 0, 2, 4 );
206
207             if( ref0b >= 0 )
208                 if( ref1b >= 0 ) x264_mb_mc_01xywh( h, 2, 0, 2, 4 );
209                 else             x264_mb_mc_0xywh ( h, 2, 0, 2, 4 );
210             else                 x264_mb_mc_1xywh ( h, 2, 0, 2, 4 );
211         }
212     }
213 }
214
215 int x264_macroblock_cache_allocate( x264_t *h )
216 {
217     int i_mb_count = h->mb.i_mb_count;
218
219     h->mb.i_mb_stride = h->sps->i_mb_width;
220     h->mb.i_b8_stride = h->sps->i_mb_width * 2;
221     h->mb.i_b4_stride = h->sps->i_mb_width * 4;
222
223     h->mb.b_interlaced = h->param.b_interlaced;
224
225     CHECKED_MALLOC( h->mb.qp, i_mb_count * sizeof(int8_t) );
226     CHECKED_MALLOC( h->mb.cbp, i_mb_count * sizeof(int16_t) );
227     CHECKED_MALLOC( h->mb.skipbp, i_mb_count * sizeof(int8_t) );
228     CHECKED_MALLOC( h->mb.mb_transform_size, i_mb_count * sizeof(int8_t) );
229     CHECKED_MALLOC( h->mb.slice_table, i_mb_count * sizeof(uint16_t) );
230     memset( h->mb.slice_table, -1, i_mb_count * sizeof(uint16_t) );
231
232     /* 0 -> 3 top(4), 4 -> 6 : left(3) */
233     CHECKED_MALLOC( h->mb.intra4x4_pred_mode, i_mb_count * 8 * sizeof(int8_t) );
234
235     /* all coeffs */
236     CHECKED_MALLOC( h->mb.non_zero_count, i_mb_count * 24 * sizeof(uint8_t) );
237
238     if( h->param.b_cabac )
239     {
240         CHECKED_MALLOC( h->mb.chroma_pred_mode, i_mb_count * sizeof(int8_t) );
241         CHECKED_MALLOC( h->mb.mvd[0], i_mb_count * sizeof( **h->mb.mvd ) );
242         CHECKED_MALLOC( h->mb.mvd[1], i_mb_count * sizeof( **h->mb.mvd ) );
243     }
244
245     for( int i = 0; i < 2; i++ )
246     {
247         int i_refs = X264_MIN(16, (i ? 1 + !!h->param.i_bframe_pyramid : h->param.i_frame_reference) ) << h->param.b_interlaced;
248         if( h->param.analyse.i_weighted_pred == X264_WEIGHTP_SMART )
249             i_refs = X264_MIN(16, i_refs + 2); //smart weights add two duplicate frames
250         else if( h->param.analyse.i_weighted_pred == X264_WEIGHTP_BLIND )
251             i_refs = X264_MIN(16, i_refs + 1); //blind weights add one duplicate frame
252
253         for( int j = !i; j < i_refs; j++ )
254             CHECKED_MALLOC( h->mb.mvr[i][j], 2 * i_mb_count * sizeof(int16_t) );
255     }
256
257     if( h->param.analyse.i_weighted_pred )
258     {
259         int i_padv = PADV << h->param.b_interlaced;
260 #define ALIGN(x,a) (((x)+((a)-1))&~((a)-1))
261         int align = h->param.cpu&X264_CPU_CACHELINE_64 ? 64 : h->param.cpu&X264_CPU_CACHELINE_32 ? 32 : 16;
262         int i_stride, luma_plane_size = 0;
263         int numweightbuf;
264
265         if( h->param.analyse.i_weighted_pred == X264_WEIGHTP_FAKE )
266         {
267             // only need buffer for lookahead
268             if( !h->param.i_sync_lookahead || h == h->thread[h->param.i_threads] )
269             {
270                 // Fake analysis only works on lowres
271                 i_stride = ALIGN( h->sps->i_mb_width*8 + 2*PADH, align );
272                 luma_plane_size = i_stride * (h->sps->i_mb_height*8+2*i_padv);
273                 // Only need 1 buffer for analysis
274                 numweightbuf = 1;
275             }
276             else
277                 numweightbuf = 0;
278         }
279         else
280         {
281             i_stride = ALIGN( h->sps->i_mb_width*16 + 2*PADH, align );
282             luma_plane_size = i_stride * (h->sps->i_mb_height*16+2*i_padv);
283
284             if( h->param.analyse.i_weighted_pred == X264_WEIGHTP_SMART )
285                 //SMART can weight one ref and one offset -1
286                 numweightbuf = 2;
287             else
288                 //blind only has one weighted copy (offset -1)
289                 numweightbuf = 1;
290         }
291
292         for( int i = 0; i < numweightbuf; i++ )
293             CHECKED_MALLOC( h->mb.p_weight_buf[i], luma_plane_size );
294 #undef ALIGN
295     }
296
297     return 0;
298 fail:
299     return -1;
300 }
301 void x264_macroblock_cache_free( x264_t *h )
302 {
303     for( int i = 0; i < 2; i++ )
304         for( int j = !i; j < 32; j++ )
305             x264_free( h->mb.mvr[i][j] );
306     for( int i = 0; i < 16; i++ )
307         x264_free( h->mb.p_weight_buf[i] );
308
309     if( h->param.b_cabac )
310     {
311         x264_free( h->mb.chroma_pred_mode );
312         x264_free( h->mb.mvd[0] );
313         x264_free( h->mb.mvd[1] );
314     }
315     x264_free( h->mb.slice_table );
316     x264_free( h->mb.intra4x4_pred_mode );
317     x264_free( h->mb.non_zero_count );
318     x264_free( h->mb.mb_transform_size );
319     x264_free( h->mb.skipbp );
320     x264_free( h->mb.cbp );
321     x264_free( h->mb.qp );
322 }
323
324 int x264_macroblock_thread_allocate( x264_t *h, int b_lookahead )
325 {
326     if( !b_lookahead )
327         for( int i = 0; i <= h->param.b_interlaced; i++ )
328         {
329             for( int j = 0; j < 3; j++ )
330             {
331                 /* shouldn't really be initialized, just silences a valgrind false-positive in predict_8x8_filter_mmx */
332                 CHECKED_MALLOCZERO( h->intra_border_backup[i][j], (h->sps->i_mb_width*16+32)>>!!j );
333                 h->intra_border_backup[i][j] += 8;
334             }
335             CHECKED_MALLOC( h->deblock_strength[i], sizeof(**h->deblock_strength) * h->sps->i_mb_width );
336         }
337
338     /* Allocate scratch buffer */
339     int scratch_size = 0;
340     if( !b_lookahead )
341     {
342         int buf_hpel = (h->thread[0]->fdec->i_width[0]+48) * sizeof(int16_t);
343         int buf_ssim = h->param.analyse.b_ssim * 8 * (h->param.i_width/4+3) * sizeof(int);
344         int me_range = X264_MIN(h->param.analyse.i_me_range, h->param.analyse.i_mv_range);
345         int buf_tesa = (h->param.analyse.i_me_method >= X264_ME_ESA) *
346             ((me_range*2+18) * sizeof(int16_t) + (me_range+4) * (me_range+1) * 4 * sizeof(mvsad_t));
347         scratch_size = X264_MAX3( buf_hpel, buf_ssim, buf_tesa );
348     }
349     int buf_mbtree = h->param.rc.b_mb_tree * ((h->sps->i_mb_width+3)&~3) * sizeof(int);
350     scratch_size = X264_MAX( scratch_size, buf_mbtree );
351     CHECKED_MALLOC( h->scratch_buffer, scratch_size );
352
353     return 0;
354 fail:
355     return -1;
356 }
357
358 void x264_macroblock_thread_free( x264_t *h, int b_lookahead )
359 {
360     if( !b_lookahead )
361         for( int i = 0; i <= h->param.b_interlaced; i++ )
362         {
363             x264_free( h->deblock_strength[i] );
364             for( int j = 0; j < 3; j++ )
365                 x264_free( h->intra_border_backup[i][j] - 8 );
366         }
367     x264_free( h->scratch_buffer );
368 }
369
370 void x264_macroblock_slice_init( x264_t *h )
371 {
372     h->mb.mv[0] = h->fdec->mv[0];
373     h->mb.mv[1] = h->fdec->mv[1];
374     h->mb.mvr[0][0] = h->fdec->mv16x16;
375     h->mb.ref[0] = h->fdec->ref[0];
376     h->mb.ref[1] = h->fdec->ref[1];
377     h->mb.type = h->fdec->mb_type;
378     h->mb.partition = h->fdec->mb_partition;
379
380     h->fdec->i_ref[0] = h->i_ref0;
381     h->fdec->i_ref[1] = h->i_ref1;
382     for( int i = 0; i < h->i_ref0; i++ )
383         h->fdec->ref_poc[0][i] = h->fref0[i]->i_poc;
384     if( h->sh.i_type == SLICE_TYPE_B )
385     {
386         for( int i = 0; i < h->i_ref1; i++ )
387             h->fdec->ref_poc[1][i] = h->fref1[i]->i_poc;
388
389         map_col_to_list0(-1) = -1;
390         map_col_to_list0(-2) = -2;
391         for( int i = 0; i < h->fref1[0]->i_ref[0]; i++ )
392         {
393             int poc = h->fref1[0]->ref_poc[0][i];
394             map_col_to_list0(i) = -2;
395             for( int j = 0; j < h->i_ref0; j++ )
396                 if( h->fref0[j]->i_poc == poc )
397                 {
398                     map_col_to_list0(i) = j;
399                     break;
400                 }
401         }
402     }
403     else if( h->sh.i_type == SLICE_TYPE_P )
404     {
405         memset( h->mb.cache.skip, 0, sizeof( h->mb.cache.skip ) );
406
407         if( h->sh.i_disable_deblocking_filter_idc != 1 && h->param.analyse.i_weighted_pred )
408         {
409             deblock_ref_table(-2) = -2;
410             deblock_ref_table(-1) = -1;
411             for( int i = 0; i < h->i_ref0 << h->sh.b_mbaff; i++ )
412             {
413                 /* Mask off high bits to avoid frame num collisions with -1/-2.
414                  * In current x264 frame num values don't cover a range of more
415                  * than 32, so 6 bits is enough for uniqueness. */
416                 if( !h->mb.b_interlaced )
417                     deblock_ref_table(i) = h->fref0[i]->i_frame_num&63;
418                 else
419                     deblock_ref_table(i) = ((h->fref0[i>>1]->i_frame_num&63)<<1) + (i&1);
420             }
421         }
422     }
423
424     /* init with not available (for top right idx=7,15) */
425     memset( h->mb.cache.ref, -2, sizeof( h->mb.cache.ref ) );
426
427     if( h->i_ref0 > 0 )
428         for( int field = 0; field <= h->sh.b_mbaff; field++ )
429         {
430             int curpoc = h->fdec->i_poc + field*h->sh.i_delta_poc_bottom;
431             int refpoc = h->fref0[0]->i_poc;
432             if( h->sh.b_mbaff && field )
433                 refpoc += h->sh.i_delta_poc_bottom;
434             int delta = curpoc - refpoc;
435
436             h->fdec->inv_ref_poc[field] = (256 + delta/2) / delta;
437         }
438
439     h->mb.i_neighbour4[6] =
440     h->mb.i_neighbour4[9] =
441     h->mb.i_neighbour4[12] =
442     h->mb.i_neighbour4[14] = MB_LEFT|MB_TOP|MB_TOPLEFT|MB_TOPRIGHT;
443     h->mb.i_neighbour4[3] =
444     h->mb.i_neighbour4[7] =
445     h->mb.i_neighbour4[11] =
446     h->mb.i_neighbour4[13] =
447     h->mb.i_neighbour4[15] =
448     h->mb.i_neighbour8[3] = MB_LEFT|MB_TOP|MB_TOPLEFT;
449 }
450
451 void x264_macroblock_thread_init( x264_t *h )
452 {
453     h->mb.i_me_method = h->param.analyse.i_me_method;
454     h->mb.i_subpel_refine = h->param.analyse.i_subpel_refine;
455     if( h->sh.i_type == SLICE_TYPE_B && (h->mb.i_subpel_refine == 6 || h->mb.i_subpel_refine == 8) )
456         h->mb.i_subpel_refine--;
457     h->mb.b_chroma_me = h->param.analyse.b_chroma_me && h->sh.i_type == SLICE_TYPE_P
458                         && h->mb.i_subpel_refine >= 5;
459     h->mb.b_dct_decimate = h->sh.i_type == SLICE_TYPE_B ||
460                           (h->param.analyse.b_dct_decimate && h->sh.i_type != SLICE_TYPE_I);
461
462
463     /* fdec:      fenc:
464      * yyyyyyy
465      * yYYYY      YYYY
466      * yYYYY      YYYY
467      * yYYYY      YYYY
468      * yYYYY      YYYY
469      * uuu vvv    UUVV
470      * uUU vVV    UUVV
471      * uUU vVV
472      */
473     h->mb.pic.p_fenc[0] = h->mb.pic.fenc_buf;
474     h->mb.pic.p_fenc[1] = h->mb.pic.fenc_buf + 16*FENC_STRIDE;
475     h->mb.pic.p_fenc[2] = h->mb.pic.fenc_buf + 16*FENC_STRIDE + 8;
476     h->mb.pic.p_fdec[0] = h->mb.pic.fdec_buf + 2*FDEC_STRIDE;
477     h->mb.pic.p_fdec[1] = h->mb.pic.fdec_buf + 19*FDEC_STRIDE;
478     h->mb.pic.p_fdec[2] = h->mb.pic.fdec_buf + 19*FDEC_STRIDE + 16;
479 }
480
481 void x264_prefetch_fenc( x264_t *h, x264_frame_t *fenc, int i_mb_x, int i_mb_y )
482 {
483     int stride_y  = fenc->i_stride[0];
484     int stride_uv = fenc->i_stride[1];
485     int off_y = 16 * (i_mb_x + i_mb_y * stride_y);
486     int off_uv = 8 * (i_mb_x + i_mb_y * stride_uv);
487     h->mc.prefetch_fenc( fenc->plane[0]+off_y, stride_y,
488                          fenc->plane[1+(i_mb_x&1)]+off_uv, stride_uv, i_mb_x );
489 }
490
491 static NOINLINE void copy_column8( uint8_t *dst, uint8_t *src )
492 {
493     // input pointers are offset by 4 rows because that's faster (smaller instruction size on x86)
494     for( int i = -4; i < 4; i++ )
495         dst[i*FDEC_STRIDE] = src[i*FDEC_STRIDE];
496 }
497
498 static void ALWAYS_INLINE x264_macroblock_load_pic_pointers( x264_t *h, int mb_x, int mb_y, int i )
499 {
500     const int w = (i == 0 ? 16 : 8);
501     const int i_stride = h->fdec->i_stride[!!i];
502     const int i_stride2 = i_stride << h->mb.b_interlaced;
503     const int i_pix_offset = h->mb.b_interlaced
504                            ? w * (mb_x + (mb_y&~1) * i_stride) + (mb_y&1) * i_stride
505                            : w * (mb_x + mb_y * i_stride);
506     const uint8_t *plane_fdec = &h->fdec->plane[i][i_pix_offset];
507     const uint8_t *intra_fdec = &h->intra_border_backup[mb_y & h->sh.b_mbaff][i][mb_x*16>>!!i];
508     int ref_pix_offset[2] = { i_pix_offset, i_pix_offset };
509     x264_frame_t **fref[2] = { h->fref0, h->fref1 };
510     if( h->mb.b_interlaced )
511         ref_pix_offset[1] += (1-2*(mb_y&1)) * i_stride;
512     h->mb.pic.i_stride[i] = i_stride2;
513     h->mb.pic.p_fenc_plane[i] = &h->fenc->plane[i][i_pix_offset];
514     h->mc.copy[i?PIXEL_8x8:PIXEL_16x16]( h->mb.pic.p_fenc[i], FENC_STRIDE,
515         h->mb.pic.p_fenc_plane[i], i_stride2, w );
516     memcpy( &h->mb.pic.p_fdec[i][-1-FDEC_STRIDE], intra_fdec-1, w*3/2+1 );
517     if( h->mb.b_interlaced )
518         for( int j = 0; j < w; j++ )
519             h->mb.pic.p_fdec[i][-1+j*FDEC_STRIDE] = plane_fdec[-1+j*i_stride2];
520     for( int j = 0; j < h->mb.pic.i_fref[0]; j++ )
521     {
522         h->mb.pic.p_fref[0][j][i==0 ? 0:i+3] = &fref[0][j >> h->mb.b_interlaced]->plane[i][ref_pix_offset[j&1]];
523         if( i == 0 )
524         {
525             for( int k = 1; k < 4; k++ )
526                 h->mb.pic.p_fref[0][j][k] = &fref[0][j >> h->mb.b_interlaced]->filtered[k][ref_pix_offset[j&1]];
527             if( h->sh.weight[j][0].weightfn )
528                 h->mb.pic.p_fref_w[j] = &h->fenc->weighted[j >> h->mb.b_interlaced][ref_pix_offset[j&1]];
529             else
530                 h->mb.pic.p_fref_w[j] = h->mb.pic.p_fref[0][j][0];
531         }
532     }
533     if( h->sh.i_type == SLICE_TYPE_B )
534         for( int j = 0; j < h->mb.pic.i_fref[1]; j++ )
535         {
536             h->mb.pic.p_fref[1][j][i==0 ? 0:i+3] = &fref[1][j >> h->mb.b_interlaced]->plane[i][ref_pix_offset[j&1]];
537             if( i == 0 )
538                 for( int k = 1; k < 4; k++ )
539                     h->mb.pic.p_fref[1][j][k] = &fref[1][j >> h->mb.b_interlaced]->filtered[k][ref_pix_offset[j&1]];
540         }
541 }
542
543 static void inline x264_macroblock_cache_load_neighbours( x264_t *h, int mb_x, int mb_y )
544 {
545     int top = (mb_y - (1 << h->mb.b_interlaced)) * h->mb.i_mb_stride + mb_x;
546
547     h->mb.i_mb_x = mb_x;
548     h->mb.i_mb_y = mb_y;
549     h->mb.i_mb_xy = mb_y * h->mb.i_mb_stride + mb_x;
550     h->mb.i_b8_xy = 2*(mb_y * h->mb.i_b8_stride + mb_x);
551     h->mb.i_b4_xy = 4*(mb_y * h->mb.i_b4_stride + mb_x);
552     h->mb.i_neighbour = 0;
553     h->mb.i_neighbour_intra = 0;
554     h->mb.i_neighbour_frame = 0;
555     h->mb.i_mb_top_xy = -1;
556     h->mb.i_mb_left_xy = -1;
557     h->mb.i_mb_topleft_xy = -1;
558     h->mb.i_mb_topright_xy = -1;
559     h->mb.i_mb_type_top = -1;
560     h->mb.i_mb_type_left = -1;
561     h->mb.i_mb_type_topleft = -1;
562     h->mb.i_mb_type_topright = -1;
563
564     if( mb_x > 0 )
565     {
566         h->mb.i_neighbour_frame |= MB_LEFT;
567         h->mb.i_mb_left_xy = h->mb.i_mb_xy - 1;
568         h->mb.i_mb_type_left = h->mb.type[h->mb.i_mb_left_xy];
569         if( h->mb.i_mb_xy > h->sh.i_first_mb )
570         {
571             h->mb.i_neighbour |= MB_LEFT;
572
573             if( !h->param.b_constrained_intra || IS_INTRA( h->mb.i_mb_type_left ) )
574                 h->mb.i_neighbour_intra |= MB_LEFT;
575         }
576     }
577
578     /* We can't predict from the previous threadslice since it hasn't been encoded yet. */
579     if( (h->i_threadslice_start >> h->mb.b_interlaced) != (mb_y >> h->mb.b_interlaced) )
580     {
581         if( top >= 0 )
582         {
583             h->mb.i_neighbour_frame |= MB_TOP;
584             h->mb.i_mb_top_xy = top;
585             h->mb.i_mb_type_top = h->mb.type[h->mb.i_mb_top_xy];
586             if( top >= h->sh.i_first_mb )
587             {
588                 h->mb.i_neighbour |= MB_TOP;
589
590                 if( !h->param.b_constrained_intra || IS_INTRA( h->mb.i_mb_type_top ) )
591                     h->mb.i_neighbour_intra |= MB_TOP;
592
593                 /* We only need to prefetch the top blocks because the left was just written
594                  * to as part of the previous cache_save.  Since most target CPUs use write-allocate
595                  * caches, left blocks are near-guaranteed to be in L1 cache.  Top--not so much. */
596                 x264_prefetch( &h->mb.cbp[top] );
597                 x264_prefetch( h->mb.intra4x4_pred_mode[top] );
598                 x264_prefetch( &h->mb.non_zero_count[top][12] );
599                 /* These aren't always allocated, but prefetching an invalid address can't hurt. */
600                 x264_prefetch( &h->mb.mb_transform_size[top] );
601                 x264_prefetch( &h->mb.skipbp[top] );
602             }
603         }
604
605         if( mb_x > 0 && top - 1 >= 0  )
606         {
607             h->mb.i_neighbour_frame |= MB_TOPLEFT;
608             h->mb.i_mb_topleft_xy = top - 1;
609             h->mb.i_mb_type_topleft = h->mb.type[h->mb.i_mb_topleft_xy];
610             if( top - 1 >= h->sh.i_first_mb )
611             {
612                 h->mb.i_neighbour |= MB_TOPLEFT;
613
614                 if( !h->param.b_constrained_intra || IS_INTRA( h->mb.i_mb_type_topleft ) )
615                     h->mb.i_neighbour_intra |= MB_TOPLEFT;
616             }
617         }
618
619         if( mb_x < h->sps->i_mb_width - 1 && top + 1 >= 0 )
620         {
621             h->mb.i_neighbour_frame |= MB_TOPRIGHT;
622             h->mb.i_mb_topright_xy = top + 1;
623             h->mb.i_mb_type_topright = h->mb.type[h->mb.i_mb_topright_xy];
624             if( top + 1 >= h->sh.i_first_mb )
625             {
626                 h->mb.i_neighbour |= MB_TOPRIGHT;
627
628                 if( !h->param.b_constrained_intra || IS_INTRA( h->mb.i_mb_type_topright ) )
629                     h->mb.i_neighbour_intra |= MB_TOPRIGHT;
630             }
631         }
632     }
633 }
634
635 void x264_macroblock_cache_load( x264_t *h, int mb_x, int mb_y )
636 {
637     x264_macroblock_cache_load_neighbours( h, mb_x, mb_y );
638
639     int left = h->mb.i_mb_left_xy;
640     int top  = h->mb.i_mb_top_xy;
641     int top_y = mb_y - (1 << h->mb.b_interlaced);
642     int top_8x8 = (2*top_y+1) * h->mb.i_b8_stride + 2*mb_x;
643     int top_4x4 = (4*top_y+3) * h->mb.i_b4_stride + 4*mb_x;
644
645     /* GCC pessimizes direct loads from heap-allocated arrays due to aliasing. */
646     /* By only dereferencing them once, we avoid this issue. */
647     int8_t (*i4x4)[8] = h->mb.intra4x4_pred_mode;
648     uint8_t (*nnz)[24] = h->mb.non_zero_count;
649     int16_t *cbp = h->mb.cbp;
650
651     /* load cache */
652     if( h->mb.i_neighbour & MB_TOP )
653     {
654         h->mb.cache.i_cbp_top = cbp[top];
655         /* load intra4x4 */
656         CP32( &h->mb.cache.intra4x4_pred_mode[x264_scan8[0] - 8], &i4x4[top][0] );
657
658         /* load non_zero_count */
659         CP32( &h->mb.cache.non_zero_count[x264_scan8[0] - 8], &nnz[top][12] );
660         /* shift because x264_scan8[16] is misaligned */
661         M32( &h->mb.cache.non_zero_count[x264_scan8[16+0] - 9] ) = M16( &nnz[top][18] ) << 8;
662         M32( &h->mb.cache.non_zero_count[x264_scan8[16+4] - 9] ) = M16( &nnz[top][22] ) << 8;
663     }
664     else
665     {
666         h->mb.cache.i_cbp_top = -1;
667
668         /* load intra4x4 */
669         M32( &h->mb.cache.intra4x4_pred_mode[x264_scan8[0] - 8] ) = 0xFFFFFFFFU;
670
671         /* load non_zero_count */
672         M32( &h->mb.cache.non_zero_count[x264_scan8[   0] - 8] ) = 0x80808080U;
673         M32( &h->mb.cache.non_zero_count[x264_scan8[16+0] - 9] ) = 0x80808080U;
674         M32( &h->mb.cache.non_zero_count[x264_scan8[16+4] - 9] ) = 0x80808080U;
675     }
676
677     if( h->mb.i_neighbour & MB_LEFT )
678     {
679         h->mb.cache.i_cbp_left = cbp[left];
680
681         /* load intra4x4 */
682         h->mb.cache.intra4x4_pred_mode[x264_scan8[0 ] - 1] = i4x4[left][4];
683         h->mb.cache.intra4x4_pred_mode[x264_scan8[2 ] - 1] = i4x4[left][5];
684         h->mb.cache.intra4x4_pred_mode[x264_scan8[8 ] - 1] = i4x4[left][6];
685         h->mb.cache.intra4x4_pred_mode[x264_scan8[10] - 1] = i4x4[left][3];
686
687         /* load non_zero_count */
688         h->mb.cache.non_zero_count[x264_scan8[0 ] - 1] = nnz[left][3];
689         h->mb.cache.non_zero_count[x264_scan8[2 ] - 1] = nnz[left][7];
690         h->mb.cache.non_zero_count[x264_scan8[8 ] - 1] = nnz[left][11];
691         h->mb.cache.non_zero_count[x264_scan8[10] - 1] = nnz[left][15];
692
693         h->mb.cache.non_zero_count[x264_scan8[16+0] - 1] = nnz[left][16+1];
694         h->mb.cache.non_zero_count[x264_scan8[16+2] - 1] = nnz[left][16+3];
695
696         h->mb.cache.non_zero_count[x264_scan8[16+4+0] - 1] = nnz[left][16+4+1];
697         h->mb.cache.non_zero_count[x264_scan8[16+4+2] - 1] = nnz[left][16+4+3];
698
699         /* Finish the prefetching */
700         if( h->sh.i_type != SLICE_TYPE_I )
701             for( int l = 0; l < (h->sh.i_type == SLICE_TYPE_B) + 1; l++ )
702             {
703                 x264_prefetch( &h->mb.mv[l][top_4x4-1] );
704                 /* Top right being not in the same cacheline as top left will happen
705                  * once every 4 MBs, so one extra prefetch is worthwhile */
706                 x264_prefetch( &h->mb.mv[l][top_4x4+4] );
707                 x264_prefetch( &h->mb.ref[l][top_8x8-1] );
708                 x264_prefetch( &h->mb.mvd[l][top] );
709             }
710     }
711     else
712     {
713         h->mb.cache.i_cbp_left = -1;
714
715         h->mb.cache.intra4x4_pred_mode[x264_scan8[0 ] - 1] =
716         h->mb.cache.intra4x4_pred_mode[x264_scan8[2 ] - 1] =
717         h->mb.cache.intra4x4_pred_mode[x264_scan8[8 ] - 1] =
718         h->mb.cache.intra4x4_pred_mode[x264_scan8[10] - 1] = -1;
719
720         /* load non_zero_count */
721         h->mb.cache.non_zero_count[x264_scan8[0 ] - 1] =
722         h->mb.cache.non_zero_count[x264_scan8[2 ] - 1] =
723         h->mb.cache.non_zero_count[x264_scan8[8 ] - 1] =
724         h->mb.cache.non_zero_count[x264_scan8[10] - 1] =
725         h->mb.cache.non_zero_count[x264_scan8[16+0] - 1] =
726         h->mb.cache.non_zero_count[x264_scan8[16+2] - 1] =
727         h->mb.cache.non_zero_count[x264_scan8[16+4+0] - 1] =
728         h->mb.cache.non_zero_count[x264_scan8[16+4+2] - 1] = 0x80;
729     }
730
731     if( h->pps->b_transform_8x8_mode )
732     {
733         h->mb.cache.i_neighbour_transform_size =
734             ( (h->mb.i_neighbour & MB_LEFT) && h->mb.mb_transform_size[left] )
735           + ( (h->mb.i_neighbour & MB_TOP) && h->mb.mb_transform_size[top]  );
736     }
737
738     if( h->sh.b_mbaff )
739     {
740         h->mb.pic.i_fref[0] = h->i_ref0 << h->mb.b_interlaced;
741         h->mb.pic.i_fref[1] = h->i_ref1 << h->mb.b_interlaced;
742         h->mb.cache.i_neighbour_interlaced =
743             !!(h->mb.i_neighbour & MB_LEFT)
744           + !!(h->mb.i_neighbour & MB_TOP);
745     }
746
747     if( !h->mb.b_interlaced )
748     {
749         copy_column8( h->mb.pic.p_fdec[0]-1+ 4*FDEC_STRIDE, h->mb.pic.p_fdec[0]+15+ 4*FDEC_STRIDE );
750         copy_column8( h->mb.pic.p_fdec[0]-1+12*FDEC_STRIDE, h->mb.pic.p_fdec[0]+15+12*FDEC_STRIDE );
751         copy_column8( h->mb.pic.p_fdec[1]-1+ 4*FDEC_STRIDE, h->mb.pic.p_fdec[1]+ 7+ 4*FDEC_STRIDE );
752         copy_column8( h->mb.pic.p_fdec[2]-1+ 4*FDEC_STRIDE, h->mb.pic.p_fdec[2]+ 7+ 4*FDEC_STRIDE );
753     }
754
755     /* load picture pointers */
756     x264_macroblock_load_pic_pointers( h, mb_x, mb_y, 0 );
757     x264_macroblock_load_pic_pointers( h, mb_x, mb_y, 1 );
758     x264_macroblock_load_pic_pointers( h, mb_x, mb_y, 2 );
759
760     if( h->fdec->integral )
761     {
762         int offset = 16 * (mb_x + mb_y * h->fdec->i_stride[0]);
763         for( int i = 0; i < h->mb.pic.i_fref[0]; i++ )
764             h->mb.pic.p_integral[0][i] = &h->fref0[i]->integral[offset];
765         for( int i = 0; i < h->mb.pic.i_fref[1]; i++ )
766             h->mb.pic.p_integral[1][i] = &h->fref1[i]->integral[offset];
767     }
768
769     x264_prefetch_fenc( h, h->fenc, mb_x, mb_y );
770
771     /* load ref/mv/mvd */
772     if( h->sh.i_type != SLICE_TYPE_I )
773     {
774         int s8x8 = h->mb.i_b8_stride;
775         int s4x4 = h->mb.i_b4_stride;
776
777         for( int l = 0; l < (h->sh.i_type == SLICE_TYPE_B) + 1; l++ )
778         {
779             int16_t (*mv)[2] = h->mb.mv[l];
780             int8_t *ref = h->mb.ref[l];
781
782             int i8 = x264_scan8[0] - 1 - 1*8;
783             if( h->mb.i_neighbour & MB_TOPLEFT )
784             {
785                 h->mb.cache.ref[l][i8] = ref[top_8x8 - 1];
786                 CP32( h->mb.cache.mv[l][i8], mv[top_4x4 - 1] );
787             }
788             else
789             {
790                 h->mb.cache.ref[l][i8] = -2;
791                 M32( h->mb.cache.mv[l][i8] ) = 0;
792             }
793
794             i8 = x264_scan8[0] - 8;
795             if( h->mb.i_neighbour & MB_TOP )
796             {
797                 h->mb.cache.ref[l][i8+0] =
798                 h->mb.cache.ref[l][i8+1] = ref[top_8x8 + 0];
799                 h->mb.cache.ref[l][i8+2] =
800                 h->mb.cache.ref[l][i8+3] = ref[top_8x8 + 1];
801                 CP128( h->mb.cache.mv[l][i8], mv[top_4x4] );
802             }
803             else
804             {
805                 M128( h->mb.cache.mv[l][i8] ) = M128_ZERO;
806                 M32( &h->mb.cache.ref[l][i8] ) = (uint8_t)(-2) * 0x01010101U;
807             }
808
809             i8 = x264_scan8[0] + 4 - 1*8;
810             if( h->mb.i_neighbour & MB_TOPRIGHT )
811             {
812                 h->mb.cache.ref[l][i8] = ref[top_8x8 + 2];
813                 CP32( h->mb.cache.mv[l][i8], mv[top_4x4 + 4] );
814             }
815             else
816                  h->mb.cache.ref[l][i8] = -2;
817
818             i8 = x264_scan8[0] - 1;
819             if( h->mb.i_neighbour & MB_LEFT )
820             {
821                 const int ir = h->mb.i_b8_xy - 1;
822                 const int iv = h->mb.i_b4_xy - 1;
823                 h->mb.cache.ref[l][i8+0*8] =
824                 h->mb.cache.ref[l][i8+1*8] = ref[ir + 0*s8x8];
825                 h->mb.cache.ref[l][i8+2*8] =
826                 h->mb.cache.ref[l][i8+3*8] = ref[ir + 1*s8x8];
827
828                 CP32( h->mb.cache.mv[l][i8+0*8], mv[iv + 0*s4x4] );
829                 CP32( h->mb.cache.mv[l][i8+1*8], mv[iv + 1*s4x4] );
830                 CP32( h->mb.cache.mv[l][i8+2*8], mv[iv + 2*s4x4] );
831                 CP32( h->mb.cache.mv[l][i8+3*8], mv[iv + 3*s4x4] );
832             }
833             else
834             {
835                 for( int i = 0; i < 4; i++ )
836                 {
837                     h->mb.cache.ref[l][i8+i*8] = -2;
838                     M32( h->mb.cache.mv[l][i8+i*8] ) = 0;
839                 }
840             }
841
842             if( h->param.b_cabac )
843             {
844                 uint8_t (*mvd)[8][2] = h->mb.mvd[l];
845                 if( h->mb.i_neighbour & MB_TOP )
846                     CP64( h->mb.cache.mvd[l][x264_scan8[0] - 8], mvd[top][0] );
847                 else
848                     M64( h->mb.cache.mvd[l][x264_scan8[0] - 8] ) = 0;
849
850                 if( h->mb.i_neighbour & MB_LEFT )
851                 {
852                     CP16( h->mb.cache.mvd[l][x264_scan8[0 ] - 1], mvd[left][4] );
853                     CP16( h->mb.cache.mvd[l][x264_scan8[2 ] - 1], mvd[left][5] );
854                     CP16( h->mb.cache.mvd[l][x264_scan8[8 ] - 1], mvd[left][6] );
855                     CP16( h->mb.cache.mvd[l][x264_scan8[10] - 1], mvd[left][3] );
856                 }
857                 else
858                     for( int i = 0; i < 4; i++ )
859                         M16( h->mb.cache.mvd[l][x264_scan8[0]-1+i*8] ) = 0;
860             }
861         }
862
863         /* load skip */
864         if( h->sh.i_type == SLICE_TYPE_B )
865         {
866             h->mb.bipred_weight = h->mb.bipred_weight_buf[h->mb.b_interlaced&(mb_y&1)];
867             h->mb.dist_scale_factor = h->mb.dist_scale_factor_buf[h->mb.b_interlaced&(mb_y&1)];
868             if( h->param.b_cabac )
869             {
870                 uint8_t skipbp;
871                 x264_macroblock_cache_skip( h, 0, 0, 4, 4, 0 );
872                 skipbp = (h->mb.i_neighbour & MB_LEFT) ? h->mb.skipbp[left] : 0;
873                 h->mb.cache.skip[x264_scan8[0] - 1] = skipbp & 0x2;
874                 h->mb.cache.skip[x264_scan8[8] - 1] = skipbp & 0x8;
875                 skipbp = (h->mb.i_neighbour & MB_TOP) ? h->mb.skipbp[top] : 0;
876                 h->mb.cache.skip[x264_scan8[0] - 8] = skipbp & 0x4;
877                 h->mb.cache.skip[x264_scan8[4] - 8] = skipbp & 0x8;
878             }
879         }
880
881         if( h->sh.i_type == SLICE_TYPE_P )
882             x264_mb_predict_mv_pskip( h, h->mb.cache.pskip_mv );
883     }
884
885     h->mb.i_neighbour4[0] =
886     h->mb.i_neighbour8[0] = (h->mb.i_neighbour_intra & (MB_TOP|MB_LEFT|MB_TOPLEFT))
887                             | ((h->mb.i_neighbour_intra & MB_TOP) ? MB_TOPRIGHT : 0);
888     h->mb.i_neighbour4[4] =
889     h->mb.i_neighbour4[1] = MB_LEFT | ((h->mb.i_neighbour_intra & MB_TOP) ? (MB_TOP|MB_TOPLEFT|MB_TOPRIGHT) : 0);
890     h->mb.i_neighbour4[2] =
891     h->mb.i_neighbour4[8] =
892     h->mb.i_neighbour4[10] =
893     h->mb.i_neighbour8[2] = MB_TOP|MB_TOPRIGHT | ((h->mb.i_neighbour_intra & MB_LEFT) ? (MB_LEFT|MB_TOPLEFT) : 0);
894     h->mb.i_neighbour4[5] =
895     h->mb.i_neighbour8[1] = MB_LEFT | (h->mb.i_neighbour_intra & MB_TOPRIGHT)
896                             | ((h->mb.i_neighbour_intra & MB_TOP) ? MB_TOP|MB_TOPLEFT : 0);
897 }
898
899 void x264_macroblock_cache_load_neighbours_deblock( x264_t *h, int mb_x, int mb_y )
900 {
901     int deblock_on_slice_edges = h->sh.i_disable_deblocking_filter_idc != 2;
902
903     h->mb.i_neighbour = 0;
904     h->mb.i_mb_xy = mb_y * h->mb.i_mb_stride + mb_x;
905
906     if( mb_x > 0 )
907     {
908         h->mb.i_mb_left_xy = h->mb.i_mb_xy - 1;
909         if( deblock_on_slice_edges || h->mb.slice_table[h->mb.i_mb_left_xy] == h->mb.slice_table[h->mb.i_mb_xy] )
910             h->mb.i_neighbour |= MB_LEFT;
911     }
912
913     if( mb_y > h->mb.b_interlaced )
914     {
915         h->mb.i_mb_top_xy = h->mb.i_mb_xy - (h->mb.i_mb_stride << h->mb.b_interlaced);
916         if( deblock_on_slice_edges || h->mb.slice_table[h->mb.i_mb_top_xy] == h->mb.slice_table[h->mb.i_mb_xy] )
917             h->mb.i_neighbour |= MB_TOP;
918     }
919 }
920
921 void x264_macroblock_cache_load_deblock( x264_t *h )
922 {
923     if( IS_INTRA( h->mb.type[h->mb.i_mb_xy] ) )
924         return;
925
926     /* If we have multiple slices and we're deblocking on slice edges, we
927      * have to reload neighbour data. */
928     if( h->sh.i_first_mb && h->sh.i_disable_deblocking_filter_idc != 2 )
929     {
930         int old_neighbour = h->mb.i_neighbour;
931         int mb_x = h->mb.i_mb_x;
932         int mb_y = h->mb.i_mb_y;
933         x264_macroblock_cache_load_neighbours_deblock( h, mb_x, mb_y );
934         h->mb.i_neighbour &= ~old_neighbour;
935         if( h->mb.i_neighbour )
936         {
937             int top_y = mb_y - (1 << h->mb.b_interlaced);
938             int top_8x8 = (2*top_y+1) * h->mb.i_b8_stride + 2*mb_x;
939             int top_4x4 = (4*top_y+3) * h->mb.i_b4_stride + 4*mb_x;
940             int s8x8 = h->mb.i_b8_stride;
941             int s4x4 = h->mb.i_b4_stride;
942
943             uint8_t (*nnz)[24] = h->mb.non_zero_count;
944
945             if( h->mb.i_neighbour & MB_TOP )
946                 CP32( &h->mb.cache.non_zero_count[x264_scan8[0] - 8], &nnz[h->mb.i_mb_top_xy][12] );
947
948             if( h->mb.i_neighbour & MB_LEFT )
949             {
950                 int left = h->mb.i_mb_left_xy;
951                 h->mb.cache.non_zero_count[x264_scan8[0 ] - 1] = nnz[left][3];
952                 h->mb.cache.non_zero_count[x264_scan8[2 ] - 1] = nnz[left][7];
953                 h->mb.cache.non_zero_count[x264_scan8[8 ] - 1] = nnz[left][11];
954                 h->mb.cache.non_zero_count[x264_scan8[10] - 1] = nnz[left][15];
955             }
956
957             for( int l = 0; l <= (h->sh.i_type == SLICE_TYPE_B); l++ )
958             {
959                 int16_t (*mv)[2] = h->mb.mv[l];
960                 int8_t *ref = h->mb.ref[l];
961
962                 int i8 = x264_scan8[0] - 8;
963                 if( h->mb.i_neighbour & MB_TOP )
964                 {
965                     h->mb.cache.ref[l][i8+0] =
966                     h->mb.cache.ref[l][i8+1] = ref[top_8x8 + 0];
967                     h->mb.cache.ref[l][i8+2] =
968                     h->mb.cache.ref[l][i8+3] = ref[top_8x8 + 1];
969                     CP128( h->mb.cache.mv[l][i8], mv[top_4x4] );
970                 }
971
972                 i8 = x264_scan8[0] - 1;
973                 if( h->mb.i_neighbour & MB_LEFT )
974                 {
975                     int ir = h->mb.i_b8_xy - 1;
976                     int iv = h->mb.i_b4_xy - 1;
977                     h->mb.cache.ref[l][i8+0*8] =
978                     h->mb.cache.ref[l][i8+1*8] = ref[ir + 0*s8x8];
979                     h->mb.cache.ref[l][i8+2*8] =
980                     h->mb.cache.ref[l][i8+3*8] = ref[ir + 1*s8x8];
981
982                     CP32( h->mb.cache.mv[l][i8+0*8], mv[iv + 0*s4x4] );
983                     CP32( h->mb.cache.mv[l][i8+1*8], mv[iv + 1*s4x4] );
984                     CP32( h->mb.cache.mv[l][i8+2*8], mv[iv + 2*s4x4] );
985                     CP32( h->mb.cache.mv[l][i8+3*8], mv[iv + 3*s4x4] );
986                 }
987             }
988         }
989     }
990
991     if( h->param.analyse.i_weighted_pred && h->sh.i_type == SLICE_TYPE_P )
992     {
993         /* Handle reference frame duplicates */
994         int i8 = x264_scan8[0] - 8;
995         h->mb.cache.ref[0][i8+0] =
996         h->mb.cache.ref[0][i8+1] = deblock_ref_table(h->mb.cache.ref[0][i8+0]);
997         h->mb.cache.ref[0][i8+2] =
998         h->mb.cache.ref[0][i8+3] = deblock_ref_table(h->mb.cache.ref[0][i8+2]);
999
1000         i8 = x264_scan8[0] - 1;
1001         h->mb.cache.ref[0][i8+0*8] =
1002         h->mb.cache.ref[0][i8+1*8] = deblock_ref_table(h->mb.cache.ref[0][i8+0*8]);
1003         h->mb.cache.ref[0][i8+2*8] =
1004         h->mb.cache.ref[0][i8+3*8] = deblock_ref_table(h->mb.cache.ref[0][i8+2*8]);
1005
1006         int ref0 = deblock_ref_table(h->mb.cache.ref[0][x264_scan8[ 0]]);
1007         int ref1 = deblock_ref_table(h->mb.cache.ref[0][x264_scan8[ 4]]);
1008         int ref2 = deblock_ref_table(h->mb.cache.ref[0][x264_scan8[ 8]]);
1009         int ref3 = deblock_ref_table(h->mb.cache.ref[0][x264_scan8[12]]);
1010         uint32_t reftop = pack16to32( (uint8_t)ref0, (uint8_t)ref1 ) * 0x0101;
1011         uint32_t refbot = pack16to32( (uint8_t)ref2, (uint8_t)ref3 ) * 0x0101;
1012
1013         M32( &h->mb.cache.ref[0][x264_scan8[0]+8*0] ) = reftop;
1014         M32( &h->mb.cache.ref[0][x264_scan8[0]+8*1] ) = reftop;
1015         M32( &h->mb.cache.ref[0][x264_scan8[0]+8*2] ) = refbot;
1016         M32( &h->mb.cache.ref[0][x264_scan8[0]+8*3] ) = refbot;
1017     }
1018
1019     /* Munge NNZ for cavlc + 8x8dct */
1020     if( !h->param.b_cabac && h->pps->b_transform_8x8_mode )
1021     {
1022         uint8_t (*nnz)[24] = h->mb.non_zero_count;
1023         int top = h->mb.i_mb_top_xy;
1024         int left = h->mb.i_mb_left_xy;
1025
1026         if( (h->mb.i_neighbour & MB_TOP) && h->mb.mb_transform_size[top] )
1027         {
1028             int i8 = x264_scan8[0] - 8;
1029             int nnz_top0 = M16( &nnz[top][8] ) | M16( &nnz[top][12] );
1030             int nnz_top1 = M16( &nnz[top][10] ) | M16( &nnz[top][14] );
1031             M16( &h->mb.cache.non_zero_count[i8+0] ) = nnz_top0 ? 0x0101 : 0;
1032             M16( &h->mb.cache.non_zero_count[i8+2] ) = nnz_top1 ? 0x0101 : 0;
1033         }
1034
1035         if( (h->mb.i_neighbour & MB_LEFT) && h->mb.mb_transform_size[left] )
1036         {
1037             int i8 = x264_scan8[0] - 1;
1038             int nnz_left0 = M16( &nnz[left][2] ) | M16( &nnz[left][6] );
1039             int nnz_left1 = M16( &nnz[left][10] ) | M16( &nnz[left][14] );
1040             h->mb.cache.non_zero_count[i8+8*0] = !!nnz_left0;
1041             h->mb.cache.non_zero_count[i8+8*1] = !!nnz_left0;
1042             h->mb.cache.non_zero_count[i8+8*2] = !!nnz_left1;
1043             h->mb.cache.non_zero_count[i8+8*3] = !!nnz_left1;
1044         }
1045
1046         if( h->mb.mb_transform_size[h->mb.i_mb_xy] )
1047         {
1048             int nnz0 = M16( &h->mb.cache.non_zero_count[x264_scan8[ 0]] ) | M16( &h->mb.cache.non_zero_count[x264_scan8[ 2]] );
1049             int nnz1 = M16( &h->mb.cache.non_zero_count[x264_scan8[ 4]] ) | M16( &h->mb.cache.non_zero_count[x264_scan8[ 6]] );
1050             int nnz2 = M16( &h->mb.cache.non_zero_count[x264_scan8[ 8]] ) | M16( &h->mb.cache.non_zero_count[x264_scan8[10]] );
1051             int nnz3 = M16( &h->mb.cache.non_zero_count[x264_scan8[12]] ) | M16( &h->mb.cache.non_zero_count[x264_scan8[14]] );
1052             uint32_t nnztop = pack16to32( !!nnz0, !!nnz1 ) * 0x0101;
1053             uint32_t nnzbot = pack16to32( !!nnz2, !!nnz3 ) * 0x0101;
1054
1055             M32( &h->mb.cache.non_zero_count[x264_scan8[0]+8*0] ) = nnztop;
1056             M32( &h->mb.cache.non_zero_count[x264_scan8[0]+8*1] ) = nnztop;
1057             M32( &h->mb.cache.non_zero_count[x264_scan8[0]+8*2] ) = nnzbot;
1058             M32( &h->mb.cache.non_zero_count[x264_scan8[0]+8*3] ) = nnzbot;
1059         }
1060     }
1061 }
1062
1063 static void ALWAYS_INLINE x264_macroblock_store_pic( x264_t *h, int i )
1064 {
1065     int w = i ? 8 : 16;
1066     int i_stride = h->fdec->i_stride[!!i];
1067     int i_stride2 = i_stride << h->mb.b_interlaced;
1068     int i_pix_offset = h->mb.b_interlaced
1069                      ? w * (h->mb.i_mb_x + (h->mb.i_mb_y&~1) * i_stride) + (h->mb.i_mb_y&1) * i_stride
1070                      : w * (h->mb.i_mb_x + h->mb.i_mb_y * i_stride);
1071     h->mc.copy[i?PIXEL_8x8:PIXEL_16x16]( &h->fdec->plane[i][i_pix_offset], i_stride2,
1072                                          h->mb.pic.p_fdec[i], FDEC_STRIDE, w );
1073 }
1074
1075 void x264_macroblock_cache_save( x264_t *h )
1076 {
1077     const int i_mb_xy = h->mb.i_mb_xy;
1078     const int i_mb_type = x264_mb_type_fix[h->mb.i_type];
1079     const int s8x8 = h->mb.i_b8_stride;
1080     const int s4x4 = h->mb.i_b4_stride;
1081     const int i_mb_4x4 = h->mb.i_b4_xy;
1082     const int i_mb_8x8 = h->mb.i_b8_xy;
1083
1084     /* GCC pessimizes direct stores to heap-allocated arrays due to aliasing. */
1085     /* By only dereferencing them once, we avoid this issue. */
1086     int8_t *i4x4 = h->mb.intra4x4_pred_mode[i_mb_xy];
1087     uint8_t *nnz = h->mb.non_zero_count[i_mb_xy];
1088
1089     x264_macroblock_store_pic( h, 0 );
1090     x264_macroblock_store_pic( h, 1 );
1091     x264_macroblock_store_pic( h, 2 );
1092
1093     x264_prefetch_fenc( h, h->fdec, h->mb.i_mb_x, h->mb.i_mb_y );
1094
1095     h->mb.type[i_mb_xy] = i_mb_type;
1096     h->mb.slice_table[i_mb_xy] = h->sh.i_first_mb;
1097     h->mb.partition[i_mb_xy] = IS_INTRA( i_mb_type ) ? D_16x16 : h->mb.i_partition;
1098     h->mb.i_mb_prev_xy = i_mb_xy;
1099
1100     /* save intra4x4 */
1101     if( i_mb_type == I_4x4 )
1102     {
1103         CP32( &i4x4[0], &h->mb.cache.intra4x4_pred_mode[x264_scan8[10]] );
1104         M32( &i4x4[4] ) = pack8to32( h->mb.cache.intra4x4_pred_mode[x264_scan8[5] ],
1105                                      h->mb.cache.intra4x4_pred_mode[x264_scan8[7] ],
1106                                      h->mb.cache.intra4x4_pred_mode[x264_scan8[13] ], 0);
1107     }
1108     else if( !h->param.b_constrained_intra || IS_INTRA(i_mb_type) )
1109         M64( i4x4 ) = I_PRED_4x4_DC * 0x0101010101010101ULL;
1110     else
1111         M64( i4x4 ) = (uint8_t)(-1) * 0x0101010101010101ULL;
1112
1113
1114     if( i_mb_type == I_PCM )
1115     {
1116         h->mb.qp[i_mb_xy] = 0;
1117         h->mb.i_last_dqp = 0;
1118         h->mb.i_cbp_chroma = 2;
1119         h->mb.i_cbp_luma = 0xf;
1120         h->mb.cbp[i_mb_xy] = 0x72f;   /* all set */
1121         h->mb.b_transform_8x8 = 0;
1122         memset( nnz, 16, sizeof( *h->mb.non_zero_count ) );
1123         for( int i = 0; i < 24; i++ )
1124             h->mb.cache.non_zero_count[x264_scan8[i]] = 16;
1125     }
1126     else
1127     {
1128         /* save non zero count */
1129         CP32( &nnz[0*4], &h->mb.cache.non_zero_count[x264_scan8[0]+0*8] );
1130         CP32( &nnz[1*4], &h->mb.cache.non_zero_count[x264_scan8[0]+1*8] );
1131         CP32( &nnz[2*4], &h->mb.cache.non_zero_count[x264_scan8[0]+2*8] );
1132         CP32( &nnz[3*4], &h->mb.cache.non_zero_count[x264_scan8[0]+3*8] );
1133         M16( &nnz[16+0*2] ) = M32( &h->mb.cache.non_zero_count[x264_scan8[16+0*2]-1] ) >> 8;
1134         M16( &nnz[16+1*2] ) = M32( &h->mb.cache.non_zero_count[x264_scan8[16+1*2]-1] ) >> 8;
1135         M16( &nnz[16+2*2] ) = M32( &h->mb.cache.non_zero_count[x264_scan8[16+2*2]-1] ) >> 8;
1136         M16( &nnz[16+3*2] ) = M32( &h->mb.cache.non_zero_count[x264_scan8[16+3*2]-1] ) >> 8;
1137
1138         if( h->mb.i_type != I_16x16 && h->mb.i_cbp_luma == 0 && h->mb.i_cbp_chroma == 0 )
1139             h->mb.i_qp = h->mb.i_last_qp;
1140         h->mb.qp[i_mb_xy] = h->mb.i_qp;
1141         h->mb.i_last_dqp = h->mb.i_qp - h->mb.i_last_qp;
1142         h->mb.i_last_qp = h->mb.i_qp;
1143     }
1144
1145     if( h->mb.i_cbp_luma == 0 && h->mb.i_type != I_8x8 )
1146         h->mb.b_transform_8x8 = 0;
1147     h->mb.mb_transform_size[i_mb_xy] = h->mb.b_transform_8x8;
1148
1149     if( h->sh.i_type != SLICE_TYPE_I )
1150     {
1151         int16_t (*mv0)[2] = &h->mb.mv[0][i_mb_4x4];
1152         int16_t (*mv1)[2] = &h->mb.mv[1][i_mb_4x4];
1153         int8_t *ref0 = &h->mb.ref[0][i_mb_8x8];
1154         int8_t *ref1 = &h->mb.ref[1][i_mb_8x8];
1155         if( !IS_INTRA( i_mb_type ) )
1156         {
1157             ref0[0+0*s8x8] = h->mb.cache.ref[0][x264_scan8[0]];
1158             ref0[1+0*s8x8] = h->mb.cache.ref[0][x264_scan8[4]];
1159             ref0[0+1*s8x8] = h->mb.cache.ref[0][x264_scan8[8]];
1160             ref0[1+1*s8x8] = h->mb.cache.ref[0][x264_scan8[12]];
1161             CP128( &mv0[0*s4x4], h->mb.cache.mv[0][x264_scan8[0]+8*0] );
1162             CP128( &mv0[1*s4x4], h->mb.cache.mv[0][x264_scan8[0]+8*1] );
1163             CP128( &mv0[2*s4x4], h->mb.cache.mv[0][x264_scan8[0]+8*2] );
1164             CP128( &mv0[3*s4x4], h->mb.cache.mv[0][x264_scan8[0]+8*3] );
1165             if( h->sh.i_type == SLICE_TYPE_B )
1166             {
1167                 ref1[0+0*s8x8] = h->mb.cache.ref[1][x264_scan8[0]];
1168                 ref1[1+0*s8x8] = h->mb.cache.ref[1][x264_scan8[4]];
1169                 ref1[0+1*s8x8] = h->mb.cache.ref[1][x264_scan8[8]];
1170                 ref1[1+1*s8x8] = h->mb.cache.ref[1][x264_scan8[12]];
1171                 CP128( &mv1[0*s4x4], h->mb.cache.mv[1][x264_scan8[0]+8*0] );
1172                 CP128( &mv1[1*s4x4], h->mb.cache.mv[1][x264_scan8[0]+8*1] );
1173                 CP128( &mv1[2*s4x4], h->mb.cache.mv[1][x264_scan8[0]+8*2] );
1174                 CP128( &mv1[3*s4x4], h->mb.cache.mv[1][x264_scan8[0]+8*3] );
1175             }
1176         }
1177         else
1178         {
1179             M16( &ref0[0*s8x8] ) = (uint8_t)(-1) * 0x0101;
1180             M16( &ref0[1*s8x8] ) = (uint8_t)(-1) * 0x0101;
1181             M128( &mv0[0*s4x4] ) = M128_ZERO;
1182             M128( &mv0[1*s4x4] ) = M128_ZERO;
1183             M128( &mv0[2*s4x4] ) = M128_ZERO;
1184             M128( &mv0[3*s4x4] ) = M128_ZERO;
1185             if( h->sh.i_type == SLICE_TYPE_B )
1186             {
1187                 M16( &ref1[0*s8x8] ) = (uint8_t)(-1) * 0x0101;
1188                 M16( &ref1[1*s8x8] ) = (uint8_t)(-1) * 0x0101;
1189                 M128( &mv1[0*s4x4] ) = M128_ZERO;
1190                 M128( &mv1[1*s4x4] ) = M128_ZERO;
1191                 M128( &mv1[2*s4x4] ) = M128_ZERO;
1192                 M128( &mv1[3*s4x4] ) = M128_ZERO;
1193             }
1194         }
1195     }
1196
1197     if( h->param.b_cabac )
1198     {
1199         uint8_t (*mvd0)[2] = h->mb.mvd[0][i_mb_xy];
1200         uint8_t (*mvd1)[2] = h->mb.mvd[1][i_mb_xy];
1201         if( IS_INTRA(i_mb_type) && i_mb_type != I_PCM )
1202             h->mb.chroma_pred_mode[i_mb_xy] = x264_mb_pred_mode8x8c_fix[ h->mb.i_chroma_pred_mode ];
1203         else
1204             h->mb.chroma_pred_mode[i_mb_xy] = I_PRED_CHROMA_DC;
1205
1206         if( !IS_INTRA( i_mb_type ) && !IS_SKIP( i_mb_type ) && !IS_DIRECT( i_mb_type ) )
1207         {
1208             CP64( mvd0[0], h->mb.cache.mvd[0][x264_scan8[10]] );
1209             CP16( mvd0[4], h->mb.cache.mvd[0][x264_scan8[5 ]] );
1210             CP16( mvd0[5], h->mb.cache.mvd[0][x264_scan8[7 ]] );
1211             CP16( mvd0[6], h->mb.cache.mvd[0][x264_scan8[13]] );
1212             if( h->sh.i_type == SLICE_TYPE_B )
1213             {
1214                 CP64( mvd1[0], h->mb.cache.mvd[1][x264_scan8[10]] );
1215                 CP16( mvd1[4], h->mb.cache.mvd[1][x264_scan8[5 ]] );
1216                 CP16( mvd1[5], h->mb.cache.mvd[1][x264_scan8[7 ]] );
1217                 CP16( mvd1[6], h->mb.cache.mvd[1][x264_scan8[13]] );
1218             }
1219         }
1220         else
1221         {
1222             M128( mvd0[0] ) = M128_ZERO;
1223             if( h->sh.i_type == SLICE_TYPE_B )
1224                 M128( mvd1[0] ) = M128_ZERO;
1225         }
1226
1227         if( h->sh.i_type == SLICE_TYPE_B )
1228         {
1229             if( i_mb_type == B_SKIP || i_mb_type == B_DIRECT )
1230                 h->mb.skipbp[i_mb_xy] = 0xf;
1231             else if( i_mb_type == B_8x8 )
1232             {
1233                 int skipbp = ( h->mb.i_sub_partition[0] == D_DIRECT_8x8 ) << 0;
1234                 skipbp    |= ( h->mb.i_sub_partition[1] == D_DIRECT_8x8 ) << 1;
1235                 skipbp    |= ( h->mb.i_sub_partition[2] == D_DIRECT_8x8 ) << 2;
1236                 skipbp    |= ( h->mb.i_sub_partition[3] == D_DIRECT_8x8 ) << 3;
1237                 h->mb.skipbp[i_mb_xy] = skipbp;
1238             }
1239             else
1240                 h->mb.skipbp[i_mb_xy] = 0;
1241         }
1242     }
1243 }
1244
1245
1246 void x264_macroblock_bipred_init( x264_t *h )
1247 {
1248     for( int field = 0; field <= h->sh.b_mbaff; field++ )
1249         for( int i_ref0 = 0; i_ref0 < (h->i_ref0<<h->sh.b_mbaff); i_ref0++ )
1250         {
1251             int poc0 = h->fref0[i_ref0>>h->sh.b_mbaff]->i_poc;
1252             if( h->sh.b_mbaff && field^(i_ref0&1) )
1253                 poc0 += h->sh.i_delta_poc_bottom;
1254             for( int i_ref1 = 0; i_ref1 < (h->i_ref1<<h->sh.b_mbaff); i_ref1++ )
1255             {
1256                 int dist_scale_factor;
1257                 int poc1 = h->fref1[i_ref1>>h->sh.b_mbaff]->i_poc;
1258                 if( h->sh.b_mbaff && field^(i_ref1&1) )
1259                     poc1 += h->sh.i_delta_poc_bottom;
1260                 int cur_poc = h->fdec->i_poc + field*h->sh.i_delta_poc_bottom;
1261                 int td = x264_clip3( poc1 - poc0, -128, 127 );
1262                 if( td == 0 /* || pic0 is a long-term ref */ )
1263                     dist_scale_factor = 256;
1264                 else
1265                 {
1266                     int tb = x264_clip3( cur_poc - poc0, -128, 127 );
1267                     int tx = (16384 + (abs(td) >> 1)) / td;
1268                     dist_scale_factor = x264_clip3( (tb * tx + 32) >> 6, -1024, 1023 );
1269                 }
1270
1271                 h->mb.dist_scale_factor_buf[field][i_ref0][i_ref1] = dist_scale_factor;
1272
1273                 dist_scale_factor >>= 2;
1274                 if( h->param.analyse.b_weighted_bipred
1275                       && dist_scale_factor >= -64
1276                       && dist_scale_factor <= 128 )
1277                 {
1278                     h->mb.bipred_weight_buf[field][i_ref0][i_ref1] = 64 - dist_scale_factor;
1279                     // ssse3 implementation of biweight doesn't support the extrema.
1280                     // if we ever generate them, we'll have to drop that optimization.
1281                     assert( dist_scale_factor >= -63 && dist_scale_factor <= 127 );
1282                 }
1283                 else
1284                     h->mb.bipred_weight_buf[field][i_ref0][i_ref1] = 32;
1285             }
1286         }
1287 }
1288