Fix 8x8dct+slices+no sliced threads+cavlc+deblock
[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( pixel, tmp0,[16*16] );
98     ALIGNED_ARRAY_16( pixel, tmp1,[16*16] );
99     pixel *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->mb.i_mb_width;
220     h->mb.i_b8_stride = h->mb.i_mb_width * 2;
221     h->mb.i_b4_stride = h->mb.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         int align = h->param.cpu&X264_CPU_CACHELINE_64 ? 64 : h->param.cpu&X264_CPU_CACHELINE_32 ? 32 : 16;
261         int i_stride, luma_plane_size = 0;
262         int numweightbuf;
263
264         if( h->param.analyse.i_weighted_pred == X264_WEIGHTP_FAKE )
265         {
266             // only need buffer for lookahead
267             if( !h->param.i_sync_lookahead || h == h->thread[h->param.i_threads] )
268             {
269                 // Fake analysis only works on lowres
270                 i_stride = ALIGN( h->mb.i_mb_width*8 + 2*PADH, align );
271                 luma_plane_size = i_stride * (h->mb.i_mb_height*8+2*i_padv);
272                 // Only need 1 buffer for analysis
273                 numweightbuf = 1;
274             }
275             else
276                 numweightbuf = 0;
277         }
278         else
279         {
280             i_stride = ALIGN( h->mb.i_mb_width*16 + 2*PADH, align );
281             luma_plane_size = i_stride * (h->mb.i_mb_height*16+2*i_padv);
282
283             if( h->param.analyse.i_weighted_pred == X264_WEIGHTP_SMART )
284                 //SMART can weight one ref and one offset -1
285                 numweightbuf = 2;
286             else
287                 //blind only has one weighted copy (offset -1)
288                 numweightbuf = 1;
289         }
290
291         for( int i = 0; i < numweightbuf; i++ )
292             CHECKED_MALLOC( h->mb.p_weight_buf[i], luma_plane_size * sizeof(pixel) );
293     }
294
295     return 0;
296 fail:
297     return -1;
298 }
299 void x264_macroblock_cache_free( x264_t *h )
300 {
301     for( int i = 0; i < 2; i++ )
302         for( int j = !i; j < 32; j++ )
303             x264_free( h->mb.mvr[i][j] );
304     for( int i = 0; i < 16; i++ )
305         x264_free( h->mb.p_weight_buf[i] );
306
307     if( h->param.b_cabac )
308     {
309         x264_free( h->mb.chroma_pred_mode );
310         x264_free( h->mb.mvd[0] );
311         x264_free( h->mb.mvd[1] );
312     }
313     x264_free( h->mb.slice_table );
314     x264_free( h->mb.intra4x4_pred_mode );
315     x264_free( h->mb.non_zero_count );
316     x264_free( h->mb.mb_transform_size );
317     x264_free( h->mb.skipbp );
318     x264_free( h->mb.cbp );
319     x264_free( h->mb.qp );
320 }
321
322 int x264_macroblock_thread_allocate( x264_t *h, int b_lookahead )
323 {
324     if( !b_lookahead )
325         for( int i = 0; i <= h->param.b_interlaced; i++ )
326         {
327             for( int j = 0; j < 3; j++ )
328             {
329                 /* shouldn't really be initialized, just silences a valgrind false-positive in predict_8x8_filter_mmx */
330                 CHECKED_MALLOCZERO( h->intra_border_backup[i][j], ((h->mb.i_mb_width*16+32)>>!!j) * sizeof(pixel) );
331                 h->intra_border_backup[i][j] += 8;
332             }
333             CHECKED_MALLOC( h->deblock_strength[i], sizeof(**h->deblock_strength) * h->mb.i_mb_width );
334         }
335
336     /* Allocate scratch buffer */
337     int scratch_size = 0;
338     if( !b_lookahead )
339     {
340         int buf_hpel = (h->thread[0]->fdec->i_width[0]+48) * sizeof(dctcoef);
341         int buf_ssim = h->param.analyse.b_ssim * 8 * (h->param.i_width/4+3) * sizeof(int);
342         int me_range = X264_MIN(h->param.analyse.i_me_range, h->param.analyse.i_mv_range);
343         int buf_tesa = (h->param.analyse.i_me_method >= X264_ME_ESA) *
344             ((me_range*2+24) * sizeof(int16_t) + (me_range+4) * (me_range+1) * 4 * sizeof(mvsad_t));
345         scratch_size = X264_MAX3( buf_hpel, buf_ssim, buf_tesa );
346     }
347     int buf_mbtree = h->param.rc.b_mb_tree * ((h->mb.i_mb_width+3)&~3) * sizeof(int);
348     scratch_size = X264_MAX( scratch_size, buf_mbtree );
349     CHECKED_MALLOC( h->scratch_buffer, scratch_size );
350
351     return 0;
352 fail:
353     return -1;
354 }
355
356 void x264_macroblock_thread_free( x264_t *h, int b_lookahead )
357 {
358     if( !b_lookahead )
359         for( int i = 0; i <= h->param.b_interlaced; i++ )
360         {
361             x264_free( h->deblock_strength[i] );
362             for( int j = 0; j < 3; j++ )
363                 x264_free( h->intra_border_backup[i][j] - 8 );
364         }
365     x264_free( h->scratch_buffer );
366 }
367
368 void x264_macroblock_slice_init( x264_t *h )
369 {
370     h->mb.mv[0] = h->fdec->mv[0];
371     h->mb.mv[1] = h->fdec->mv[1];
372     h->mb.mvr[0][0] = h->fdec->mv16x16;
373     h->mb.ref[0] = h->fdec->ref[0];
374     h->mb.ref[1] = h->fdec->ref[1];
375     h->mb.type = h->fdec->mb_type;
376     h->mb.partition = h->fdec->mb_partition;
377
378     h->fdec->i_ref[0] = h->i_ref0;
379     h->fdec->i_ref[1] = h->i_ref1;
380     for( int i = 0; i < h->i_ref0; i++ )
381         h->fdec->ref_poc[0][i] = h->fref0[i]->i_poc;
382     if( h->sh.i_type == SLICE_TYPE_B )
383     {
384         for( int i = 0; i < h->i_ref1; i++ )
385             h->fdec->ref_poc[1][i] = h->fref1[i]->i_poc;
386
387         map_col_to_list0(-1) = -1;
388         map_col_to_list0(-2) = -2;
389         for( int i = 0; i < h->fref1[0]->i_ref[0]; i++ )
390         {
391             int poc = h->fref1[0]->ref_poc[0][i];
392             map_col_to_list0(i) = -2;
393             for( int j = 0; j < h->i_ref0; j++ )
394                 if( h->fref0[j]->i_poc == poc )
395                 {
396                     map_col_to_list0(i) = j;
397                     break;
398                 }
399         }
400     }
401     else if( h->sh.i_type == SLICE_TYPE_P )
402     {
403         memset( h->mb.cache.skip, 0, sizeof( h->mb.cache.skip ) );
404
405         if( h->sh.i_disable_deblocking_filter_idc != 1 && h->param.analyse.i_weighted_pred )
406         {
407             deblock_ref_table(-2) = -2;
408             deblock_ref_table(-1) = -1;
409             for( int i = 0; i < h->i_ref0 << h->sh.b_mbaff; i++ )
410             {
411                 /* Mask off high bits to avoid frame num collisions with -1/-2.
412                  * In current x264 frame num values don't cover a range of more
413                  * than 32, so 6 bits is enough for uniqueness. */
414                 if( !h->mb.b_interlaced )
415                     deblock_ref_table(i) = h->fref0[i]->i_frame_num&63;
416                 else
417                     deblock_ref_table(i) = ((h->fref0[i>>1]->i_frame_num&63)<<1) + (i&1);
418             }
419         }
420     }
421
422     /* init with not available (for top right idx=7,15) */
423     memset( h->mb.cache.ref, -2, sizeof( h->mb.cache.ref ) );
424
425     if( h->i_ref0 > 0 )
426         for( int field = 0; field <= h->sh.b_mbaff; field++ )
427         {
428             int curpoc = h->fdec->i_poc + field*h->sh.i_delta_poc_bottom;
429             int refpoc = h->fref0[0]->i_poc;
430             if( h->sh.b_mbaff && field )
431                 refpoc += h->sh.i_delta_poc_bottom;
432             int delta = curpoc - refpoc;
433
434             h->fdec->inv_ref_poc[field] = (256 + delta/2) / delta;
435         }
436
437     h->mb.i_neighbour4[6] =
438     h->mb.i_neighbour4[9] =
439     h->mb.i_neighbour4[12] =
440     h->mb.i_neighbour4[14] = MB_LEFT|MB_TOP|MB_TOPLEFT|MB_TOPRIGHT;
441     h->mb.i_neighbour4[3] =
442     h->mb.i_neighbour4[7] =
443     h->mb.i_neighbour4[11] =
444     h->mb.i_neighbour4[13] =
445     h->mb.i_neighbour4[15] =
446     h->mb.i_neighbour8[3] = MB_LEFT|MB_TOP|MB_TOPLEFT;
447 }
448
449 void x264_macroblock_thread_init( x264_t *h )
450 {
451     h->mb.i_me_method = h->param.analyse.i_me_method;
452     h->mb.i_subpel_refine = h->param.analyse.i_subpel_refine;
453     if( h->sh.i_type == SLICE_TYPE_B && (h->mb.i_subpel_refine == 6 || h->mb.i_subpel_refine == 8) )
454         h->mb.i_subpel_refine--;
455     h->mb.b_chroma_me = h->param.analyse.b_chroma_me && h->sh.i_type == SLICE_TYPE_P
456                         && h->mb.i_subpel_refine >= 5;
457     h->mb.b_dct_decimate = h->sh.i_type == SLICE_TYPE_B ||
458                           (h->param.analyse.b_dct_decimate && h->sh.i_type != SLICE_TYPE_I);
459
460
461     /* fdec:      fenc:
462      * yyyyyyy
463      * yYYYY      YYYY
464      * yYYYY      YYYY
465      * yYYYY      YYYY
466      * yYYYY      YYYY
467      * uuu vvv    UUVV
468      * uUU vVV    UUVV
469      * uUU vVV
470      */
471     h->mb.pic.p_fenc[0] = h->mb.pic.fenc_buf;
472     h->mb.pic.p_fenc[1] = h->mb.pic.fenc_buf + 16*FENC_STRIDE;
473     h->mb.pic.p_fenc[2] = h->mb.pic.fenc_buf + 16*FENC_STRIDE + 8;
474     h->mb.pic.p_fdec[0] = h->mb.pic.fdec_buf + 2*FDEC_STRIDE;
475     h->mb.pic.p_fdec[1] = h->mb.pic.fdec_buf + 19*FDEC_STRIDE;
476     h->mb.pic.p_fdec[2] = h->mb.pic.fdec_buf + 19*FDEC_STRIDE + 16;
477 }
478
479 void x264_prefetch_fenc( x264_t *h, x264_frame_t *fenc, int i_mb_x, int i_mb_y )
480 {
481     int stride_y  = fenc->i_stride[0];
482     int stride_uv = fenc->i_stride[1];
483     int off_y = 16 * (i_mb_x + i_mb_y * stride_y);
484     int off_uv = 8 * (i_mb_x + i_mb_y * stride_uv);
485     h->mc.prefetch_fenc( fenc->plane[0]+off_y, stride_y,
486                          fenc->plane[1+(i_mb_x&1)]+off_uv, stride_uv, i_mb_x );
487 }
488
489 static NOINLINE void copy_column8( pixel *dst, pixel *src )
490 {
491     // input pointers are offset by 4 rows because that's faster (smaller instruction size on x86)
492     for( int i = -4; i < 4; i++ )
493         dst[i*FDEC_STRIDE] = src[i*FDEC_STRIDE];
494 }
495
496 static void ALWAYS_INLINE x264_macroblock_load_pic_pointers( x264_t *h, int mb_x, int mb_y, int i, int b_interlaced )
497 {
498     const int w = (i == 0 ? 16 : 8);
499     const int i_stride = h->fdec->i_stride[!!i];
500     const int i_stride2 = i_stride << b_interlaced;
501     const int i_pix_offset = b_interlaced
502                            ? w * (mb_x + (mb_y&~1) * i_stride) + (mb_y&1) * i_stride
503                            : w * (mb_x + mb_y * i_stride);
504     const pixel *plane_fdec = &h->fdec->plane[i][i_pix_offset];
505     const pixel *intra_fdec = &h->intra_border_backup[mb_y & h->sh.b_mbaff][i][mb_x*16>>!!i];
506     int ref_pix_offset[2] = { i_pix_offset, i_pix_offset };
507     x264_frame_t **fref[2] = { h->fref0, h->fref1 };
508     if( b_interlaced )
509         ref_pix_offset[1] += (1-2*(mb_y&1)) * i_stride;
510     h->mb.pic.i_stride[i] = i_stride2;
511     h->mb.pic.p_fenc_plane[i] = &h->fenc->plane[i][i_pix_offset];
512     h->mc.copy[i?PIXEL_8x8:PIXEL_16x16]( h->mb.pic.p_fenc[i], FENC_STRIDE,
513         h->mb.pic.p_fenc_plane[i], i_stride2, w );
514     memcpy( &h->mb.pic.p_fdec[i][-1-FDEC_STRIDE], intra_fdec-1, (w*3/2+1) * sizeof(pixel) );
515     if( b_interlaced )
516         for( int j = 0; j < w; j++ )
517             h->mb.pic.p_fdec[i][-1+j*FDEC_STRIDE] = plane_fdec[-1+j*i_stride2];
518     for( int j = 0; j < h->mb.pic.i_fref[0]; j++ )
519     {
520         h->mb.pic.p_fref[0][j][i==0 ? 0:i+3] = &fref[0][j >> b_interlaced]->plane[i][ref_pix_offset[j&1]];
521         if( i == 0 )
522         {
523             for( int k = 1; k < 4; k++ )
524                 h->mb.pic.p_fref[0][j][k] = &fref[0][j >> b_interlaced]->filtered[k][ref_pix_offset[j&1]];
525             if( h->sh.weight[j][0].weightfn )
526                 h->mb.pic.p_fref_w[j] = &h->fenc->weighted[j >> b_interlaced][ref_pix_offset[j&1]];
527             else
528                 h->mb.pic.p_fref_w[j] = h->mb.pic.p_fref[0][j][0];
529         }
530     }
531     if( h->sh.i_type == SLICE_TYPE_B )
532         for( int j = 0; j < h->mb.pic.i_fref[1]; j++ )
533         {
534             h->mb.pic.p_fref[1][j][i==0 ? 0:i+3] = &fref[1][j >> b_interlaced]->plane[i][ref_pix_offset[j&1]];
535             if( i == 0 )
536                 for( int k = 1; k < 4; k++ )
537                     h->mb.pic.p_fref[1][j][k] = &fref[1][j >> b_interlaced]->filtered[k][ref_pix_offset[j&1]];
538         }
539 }
540
541 static void inline x264_macroblock_cache_load_neighbours( x264_t *h, int mb_x, int mb_y )
542 {
543     int top = (mb_y - (1 << h->mb.b_interlaced)) * h->mb.i_mb_stride + mb_x;
544
545     h->mb.i_mb_x = mb_x;
546     h->mb.i_mb_y = mb_y;
547     h->mb.i_mb_xy = mb_y * h->mb.i_mb_stride + mb_x;
548     h->mb.i_b8_xy = 2*(mb_y * h->mb.i_b8_stride + mb_x);
549     h->mb.i_b4_xy = 4*(mb_y * h->mb.i_b4_stride + mb_x);
550     h->mb.i_neighbour = 0;
551     h->mb.i_neighbour_intra = 0;
552     h->mb.i_neighbour_frame = 0;
553     h->mb.i_mb_top_xy = -1;
554     h->mb.i_mb_left_xy = -1;
555     h->mb.i_mb_topleft_xy = -1;
556     h->mb.i_mb_topright_xy = -1;
557     h->mb.i_mb_type_top = -1;
558     h->mb.i_mb_type_left = -1;
559     h->mb.i_mb_type_topleft = -1;
560     h->mb.i_mb_type_topright = -1;
561
562     if( mb_x > 0 )
563     {
564         h->mb.i_neighbour_frame |= MB_LEFT;
565         h->mb.i_mb_left_xy = h->mb.i_mb_xy - 1;
566         h->mb.i_mb_type_left = h->mb.type[h->mb.i_mb_left_xy];
567         if( h->mb.i_mb_xy > h->sh.i_first_mb )
568         {
569             h->mb.i_neighbour |= MB_LEFT;
570
571             if( !h->param.b_constrained_intra || IS_INTRA( h->mb.i_mb_type_left ) )
572                 h->mb.i_neighbour_intra |= MB_LEFT;
573         }
574     }
575
576     /* We can't predict from the previous threadslice since it hasn't been encoded yet. */
577     if( (h->i_threadslice_start >> h->mb.b_interlaced) != (mb_y >> h->mb.b_interlaced) )
578     {
579         if( top >= 0 )
580         {
581             h->mb.i_neighbour_frame |= MB_TOP;
582             h->mb.i_mb_top_xy = top;
583             h->mb.i_mb_type_top = h->mb.type[h->mb.i_mb_top_xy];
584             if( top >= h->sh.i_first_mb )
585             {
586                 h->mb.i_neighbour |= MB_TOP;
587
588                 if( !h->param.b_constrained_intra || IS_INTRA( h->mb.i_mb_type_top ) )
589                     h->mb.i_neighbour_intra |= MB_TOP;
590
591                 /* We only need to prefetch the top blocks because the left was just written
592                  * to as part of the previous cache_save.  Since most target CPUs use write-allocate
593                  * caches, left blocks are near-guaranteed to be in L1 cache.  Top--not so much. */
594                 x264_prefetch( &h->mb.cbp[top] );
595                 x264_prefetch( h->mb.intra4x4_pred_mode[top] );
596                 x264_prefetch( &h->mb.non_zero_count[top][12] );
597                 /* These aren't always allocated, but prefetching an invalid address can't hurt. */
598                 x264_prefetch( &h->mb.mb_transform_size[top] );
599                 x264_prefetch( &h->mb.skipbp[top] );
600             }
601         }
602
603         if( mb_x > 0 && top - 1 >= 0  )
604         {
605             h->mb.i_neighbour_frame |= MB_TOPLEFT;
606             h->mb.i_mb_topleft_xy = top - 1;
607             h->mb.i_mb_type_topleft = h->mb.type[h->mb.i_mb_topleft_xy];
608             if( top - 1 >= h->sh.i_first_mb )
609             {
610                 h->mb.i_neighbour |= MB_TOPLEFT;
611
612                 if( !h->param.b_constrained_intra || IS_INTRA( h->mb.i_mb_type_topleft ) )
613                     h->mb.i_neighbour_intra |= MB_TOPLEFT;
614             }
615         }
616
617         if( mb_x < h->mb.i_mb_width - 1 && top + 1 >= 0 )
618         {
619             h->mb.i_neighbour_frame |= MB_TOPRIGHT;
620             h->mb.i_mb_topright_xy = top + 1;
621             h->mb.i_mb_type_topright = h->mb.type[h->mb.i_mb_topright_xy];
622             if( top + 1 >= h->sh.i_first_mb )
623             {
624                 h->mb.i_neighbour |= MB_TOPRIGHT;
625
626                 if( !h->param.b_constrained_intra || IS_INTRA( h->mb.i_mb_type_topright ) )
627                     h->mb.i_neighbour_intra |= MB_TOPRIGHT;
628             }
629         }
630     }
631 }
632
633 void x264_macroblock_cache_load( x264_t *h, int mb_x, int mb_y )
634 {
635     x264_macroblock_cache_load_neighbours( h, mb_x, mb_y );
636
637     int left = h->mb.i_mb_left_xy;
638     int top  = h->mb.i_mb_top_xy;
639     int top_y = mb_y - (1 << h->mb.b_interlaced);
640     int s8x8 = h->mb.i_b8_stride;
641     int s4x4 = h->mb.i_b4_stride;
642     int top_8x8 = (2*top_y+1) * s8x8 + 2*mb_x;
643     int top_4x4 = (4*top_y+3) * s4x4 + 4*mb_x;
644     int lists = (1 << h->sh.i_type) & 3;
645
646     /* GCC pessimizes direct loads from heap-allocated arrays due to aliasing. */
647     /* By only dereferencing them once, we avoid this issue. */
648     int8_t (*i4x4)[8] = h->mb.intra4x4_pred_mode;
649     uint8_t (*nnz)[24] = h->mb.non_zero_count;
650     int16_t *cbp = h->mb.cbp;
651
652     /* load cache */
653     if( h->mb.i_neighbour & MB_TOP )
654     {
655         h->mb.cache.i_cbp_top = cbp[top];
656         /* load intra4x4 */
657         CP32( &h->mb.cache.intra4x4_pred_mode[x264_scan8[0] - 8], &i4x4[top][0] );
658
659         /* load non_zero_count */
660         CP32( &h->mb.cache.non_zero_count[x264_scan8[0] - 8], &nnz[top][12] );
661         /* shift because x264_scan8[16] is misaligned */
662         M32( &h->mb.cache.non_zero_count[x264_scan8[16+0] - 9] ) = M16( &nnz[top][18] ) << 8;
663         M32( &h->mb.cache.non_zero_count[x264_scan8[16+4] - 9] ) = M16( &nnz[top][22] ) << 8;
664     }
665     else
666     {
667         h->mb.cache.i_cbp_top = -1;
668
669         /* load intra4x4 */
670         M32( &h->mb.cache.intra4x4_pred_mode[x264_scan8[0] - 8] ) = 0xFFFFFFFFU;
671
672         /* load non_zero_count */
673         M32( &h->mb.cache.non_zero_count[x264_scan8[   0] - 8] ) = 0x80808080U;
674         M32( &h->mb.cache.non_zero_count[x264_scan8[16+0] - 9] ) = 0x80808080U;
675         M32( &h->mb.cache.non_zero_count[x264_scan8[16+4] - 9] ) = 0x80808080U;
676     }
677
678     if( h->mb.i_neighbour & MB_LEFT )
679     {
680         h->mb.cache.i_cbp_left = cbp[left];
681
682         /* load intra4x4 */
683         h->mb.cache.intra4x4_pred_mode[x264_scan8[0 ] - 1] = i4x4[left][4];
684         h->mb.cache.intra4x4_pred_mode[x264_scan8[2 ] - 1] = i4x4[left][5];
685         h->mb.cache.intra4x4_pred_mode[x264_scan8[8 ] - 1] = i4x4[left][6];
686         h->mb.cache.intra4x4_pred_mode[x264_scan8[10] - 1] = i4x4[left][3];
687
688         /* load non_zero_count */
689         h->mb.cache.non_zero_count[x264_scan8[0 ] - 1] = nnz[left][3];
690         h->mb.cache.non_zero_count[x264_scan8[2 ] - 1] = nnz[left][7];
691         h->mb.cache.non_zero_count[x264_scan8[8 ] - 1] = nnz[left][11];
692         h->mb.cache.non_zero_count[x264_scan8[10] - 1] = nnz[left][15];
693
694         h->mb.cache.non_zero_count[x264_scan8[16+0] - 1] = nnz[left][16+1];
695         h->mb.cache.non_zero_count[x264_scan8[16+2] - 1] = nnz[left][16+3];
696
697         h->mb.cache.non_zero_count[x264_scan8[16+4+0] - 1] = nnz[left][16+4+1];
698         h->mb.cache.non_zero_count[x264_scan8[16+4+2] - 1] = nnz[left][16+4+3];
699
700         /* Finish the prefetching */
701         for( int l = 0; l < lists; 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         x264_macroblock_load_pic_pointers( h, mb_x, mb_y, 0, 0 );
754         x264_macroblock_load_pic_pointers( h, mb_x, mb_y, 1, 0 );
755         x264_macroblock_load_pic_pointers( h, mb_x, mb_y, 2, 0 );
756     }
757     else
758     {
759         x264_macroblock_load_pic_pointers( h, mb_x, mb_y, 0, 1 );
760         x264_macroblock_load_pic_pointers( h, mb_x, mb_y, 1, 1 );
761         x264_macroblock_load_pic_pointers( h, mb_x, mb_y, 2, 1 );
762     }
763
764     if( h->fdec->integral )
765     {
766         int offset = 16 * (mb_x + mb_y * h->fdec->i_stride[0]);
767         for( int i = 0; i < h->mb.pic.i_fref[0]; i++ )
768             h->mb.pic.p_integral[0][i] = &h->fref0[i]->integral[offset];
769         for( int i = 0; i < h->mb.pic.i_fref[1]; i++ )
770             h->mb.pic.p_integral[1][i] = &h->fref1[i]->integral[offset];
771     }
772
773     x264_prefetch_fenc( h, h->fenc, mb_x, mb_y );
774
775     /* load ref/mv/mvd */
776     for( int l = 0; l < lists; l++ )
777     {
778         int16_t (*mv)[2] = h->mb.mv[l];
779         int8_t *ref = h->mb.ref[l];
780
781         int i8 = x264_scan8[0] - 1 - 1*8;
782         if( h->mb.i_neighbour & MB_TOPLEFT )
783         {
784             h->mb.cache.ref[l][i8] = ref[top_8x8 - 1];
785             CP32( h->mb.cache.mv[l][i8], mv[top_4x4 - 1] );
786         }
787         else
788         {
789             h->mb.cache.ref[l][i8] = -2;
790             M32( h->mb.cache.mv[l][i8] ) = 0;
791         }
792
793         i8 = x264_scan8[0] - 8;
794         if( h->mb.i_neighbour & MB_TOP )
795         {
796             h->mb.cache.ref[l][i8+0] =
797             h->mb.cache.ref[l][i8+1] = ref[top_8x8 + 0];
798             h->mb.cache.ref[l][i8+2] =
799             h->mb.cache.ref[l][i8+3] = ref[top_8x8 + 1];
800             CP128( h->mb.cache.mv[l][i8], mv[top_4x4] );
801         }
802         else
803         {
804             M128( h->mb.cache.mv[l][i8] ) = M128_ZERO;
805             M32( &h->mb.cache.ref[l][i8] ) = (uint8_t)(-2) * 0x01010101U;
806         }
807
808         i8 = x264_scan8[0] + 4 - 1*8;
809         if( h->mb.i_neighbour & MB_TOPRIGHT )
810         {
811             h->mb.cache.ref[l][i8] = ref[top_8x8 + 2];
812             CP32( h->mb.cache.mv[l][i8], mv[top_4x4 + 4] );
813         }
814         else
815              h->mb.cache.ref[l][i8] = -2;
816
817         i8 = x264_scan8[0] - 1;
818         if( h->mb.i_neighbour & MB_LEFT )
819         {
820             const int ir = h->mb.i_b8_xy - 1;
821             const int iv = h->mb.i_b4_xy - 1;
822             h->mb.cache.ref[l][i8+0*8] =
823             h->mb.cache.ref[l][i8+1*8] = ref[ir + 0*s8x8];
824             h->mb.cache.ref[l][i8+2*8] =
825             h->mb.cache.ref[l][i8+3*8] = ref[ir + 1*s8x8];
826
827             CP32( h->mb.cache.mv[l][i8+0*8], mv[iv + 0*s4x4] );
828             CP32( h->mb.cache.mv[l][i8+1*8], mv[iv + 1*s4x4] );
829             CP32( h->mb.cache.mv[l][i8+2*8], mv[iv + 2*s4x4] );
830             CP32( h->mb.cache.mv[l][i8+3*8], mv[iv + 3*s4x4] );
831         }
832         else
833         {
834             for( int i = 0; i < 4; i++ )
835             {
836                 h->mb.cache.ref[l][i8+i*8] = -2;
837                 M32( h->mb.cache.mv[l][i8+i*8] ) = 0;
838             }
839         }
840
841         if( h->param.b_cabac )
842         {
843             uint8_t (*mvd)[8][2] = h->mb.mvd[l];
844             if( h->mb.i_neighbour & MB_TOP )
845                 CP64( h->mb.cache.mvd[l][x264_scan8[0] - 8], mvd[top][0] );
846             else
847                 M64( h->mb.cache.mvd[l][x264_scan8[0] - 8] ) = 0;
848
849             if( h->mb.i_neighbour & MB_LEFT )
850             {
851                 CP16( h->mb.cache.mvd[l][x264_scan8[0 ] - 1], mvd[left][4] );
852                 CP16( h->mb.cache.mvd[l][x264_scan8[2 ] - 1], mvd[left][5] );
853                 CP16( h->mb.cache.mvd[l][x264_scan8[8 ] - 1], mvd[left][6] );
854                 CP16( h->mb.cache.mvd[l][x264_scan8[10] - 1], mvd[left][3] );
855             }
856             else
857                 for( int i = 0; i < 4; i++ )
858                     M16( h->mb.cache.mvd[l][x264_scan8[0]-1+i*8] ) = 0;
859         }
860     }
861
862     /* load skip */
863     if( h->sh.i_type == SLICE_TYPE_B )
864     {
865         h->mb.bipred_weight = h->mb.bipred_weight_buf[h->mb.b_interlaced&(mb_y&1)];
866         h->mb.dist_scale_factor = h->mb.dist_scale_factor_buf[h->mb.b_interlaced&(mb_y&1)];
867         if( h->param.b_cabac )
868         {
869             uint8_t skipbp;
870             x264_macroblock_cache_skip( h, 0, 0, 4, 4, 0 );
871             skipbp = (h->mb.i_neighbour & MB_LEFT) ? h->mb.skipbp[left] : 0;
872             h->mb.cache.skip[x264_scan8[0] - 1] = skipbp & 0x2;
873             h->mb.cache.skip[x264_scan8[8] - 1] = skipbp & 0x8;
874             skipbp = (h->mb.i_neighbour & MB_TOP) ? h->mb.skipbp[top] : 0;
875             h->mb.cache.skip[x264_scan8[0] - 8] = skipbp & 0x4;
876             h->mb.cache.skip[x264_scan8[4] - 8] = skipbp & 0x8;
877         }
878     }
879
880     if( h->sh.i_type == SLICE_TYPE_P )
881         x264_mb_predict_mv_pskip( h, h->mb.cache.pskip_mv );
882
883     h->mb.i_neighbour4[0] =
884     h->mb.i_neighbour8[0] = (h->mb.i_neighbour_intra & (MB_TOP|MB_LEFT|MB_TOPLEFT))
885                             | ((h->mb.i_neighbour_intra & MB_TOP) ? MB_TOPRIGHT : 0);
886     h->mb.i_neighbour4[4] =
887     h->mb.i_neighbour4[1] = MB_LEFT | ((h->mb.i_neighbour_intra & MB_TOP) ? (MB_TOP|MB_TOPLEFT|MB_TOPRIGHT) : 0);
888     h->mb.i_neighbour4[2] =
889     h->mb.i_neighbour4[8] =
890     h->mb.i_neighbour4[10] =
891     h->mb.i_neighbour8[2] = MB_TOP|MB_TOPRIGHT | ((h->mb.i_neighbour_intra & MB_LEFT) ? (MB_LEFT|MB_TOPLEFT) : 0);
892     h->mb.i_neighbour4[5] =
893     h->mb.i_neighbour8[1] = MB_LEFT | (h->mb.i_neighbour_intra & MB_TOPRIGHT)
894                             | ((h->mb.i_neighbour_intra & MB_TOP) ? MB_TOP|MB_TOPLEFT : 0);
895 }
896
897 void x264_macroblock_cache_load_neighbours_deblock( x264_t *h, int mb_x, int mb_y )
898 {
899     int deblock_on_slice_edges = h->sh.i_disable_deblocking_filter_idc != 2;
900
901     h->mb.i_neighbour = 0;
902     h->mb.i_mb_xy = mb_y * h->mb.i_mb_stride + mb_x;
903
904     if( mb_x > 0 )
905     {
906         h->mb.i_mb_left_xy = h->mb.i_mb_xy - 1;
907         if( deblock_on_slice_edges || h->mb.slice_table[h->mb.i_mb_left_xy] == h->mb.slice_table[h->mb.i_mb_xy] )
908             h->mb.i_neighbour |= MB_LEFT;
909     }
910
911     if( mb_y > h->mb.b_interlaced )
912     {
913         h->mb.i_mb_top_xy = h->mb.i_mb_xy - (h->mb.i_mb_stride << h->mb.b_interlaced);
914         if( deblock_on_slice_edges || h->mb.slice_table[h->mb.i_mb_top_xy] == h->mb.slice_table[h->mb.i_mb_xy] )
915             h->mb.i_neighbour |= MB_TOP;
916     }
917 }
918
919 void x264_macroblock_cache_load_deblock( x264_t *h )
920 {
921     if( IS_INTRA( h->mb.type[h->mb.i_mb_xy] ) )
922         return;
923
924     /* If we have multiple slices and we're deblocking on slice edges, we
925      * have to reload neighbour data. */
926     if( h->sh.i_first_mb && h->sh.i_disable_deblocking_filter_idc != 2 )
927     {
928         int old_neighbour = h->mb.i_neighbour;
929         int mb_x = h->mb.i_mb_x;
930         int mb_y = h->mb.i_mb_y;
931         x264_macroblock_cache_load_neighbours_deblock( h, mb_x, mb_y );
932         int new_neighbour = h->mb.i_neighbour;
933         h->mb.i_neighbour &= ~old_neighbour;
934         if( h->mb.i_neighbour )
935         {
936             int top_y = mb_y - (1 << h->mb.b_interlaced);
937             int top_8x8 = (2*top_y+1) * h->mb.i_b8_stride + 2*mb_x;
938             int top_4x4 = (4*top_y+3) * h->mb.i_b4_stride + 4*mb_x;
939             int s8x8 = h->mb.i_b8_stride;
940             int s4x4 = h->mb.i_b4_stride;
941
942             uint8_t (*nnz)[24] = h->mb.non_zero_count;
943
944             if( h->mb.i_neighbour & MB_TOP )
945                 CP32( &h->mb.cache.non_zero_count[x264_scan8[0] - 8], &nnz[h->mb.i_mb_top_xy][12] );
946
947             if( h->mb.i_neighbour & MB_LEFT )
948             {
949                 int left = h->mb.i_mb_left_xy;
950                 h->mb.cache.non_zero_count[x264_scan8[0 ] - 1] = nnz[left][3];
951                 h->mb.cache.non_zero_count[x264_scan8[2 ] - 1] = nnz[left][7];
952                 h->mb.cache.non_zero_count[x264_scan8[8 ] - 1] = nnz[left][11];
953                 h->mb.cache.non_zero_count[x264_scan8[10] - 1] = nnz[left][15];
954             }
955
956             for( int l = 0; l <= (h->sh.i_type == SLICE_TYPE_B); l++ )
957             {
958                 int16_t (*mv)[2] = h->mb.mv[l];
959                 int8_t *ref = h->mb.ref[l];
960
961                 int i8 = x264_scan8[0] - 8;
962                 if( h->mb.i_neighbour & MB_TOP )
963                 {
964                     h->mb.cache.ref[l][i8+0] =
965                     h->mb.cache.ref[l][i8+1] = ref[top_8x8 + 0];
966                     h->mb.cache.ref[l][i8+2] =
967                     h->mb.cache.ref[l][i8+3] = ref[top_8x8 + 1];
968                     CP128( h->mb.cache.mv[l][i8], mv[top_4x4] );
969                 }
970
971                 i8 = x264_scan8[0] - 1;
972                 if( h->mb.i_neighbour & MB_LEFT )
973                 {
974                     int ir = h->mb.i_b8_xy - 1;
975                     int iv = h->mb.i_b4_xy - 1;
976                     h->mb.cache.ref[l][i8+0*8] =
977                     h->mb.cache.ref[l][i8+1*8] = ref[ir + 0*s8x8];
978                     h->mb.cache.ref[l][i8+2*8] =
979                     h->mb.cache.ref[l][i8+3*8] = ref[ir + 1*s8x8];
980
981                     CP32( h->mb.cache.mv[l][i8+0*8], mv[iv + 0*s4x4] );
982                     CP32( h->mb.cache.mv[l][i8+1*8], mv[iv + 1*s4x4] );
983                     CP32( h->mb.cache.mv[l][i8+2*8], mv[iv + 2*s4x4] );
984                     CP32( h->mb.cache.mv[l][i8+3*8], mv[iv + 3*s4x4] );
985                 }
986             }
987         }
988         h->mb.i_neighbour = new_neighbour;
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( (0x3FF30 >> i_mb_type) & 1 ) /* !INTRA && !SKIP && !DIRECT */
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