1 /*****************************************************************************
2 * frame.c: frame handling
3 *****************************************************************************
4 * Copyright (C) 2003-2010 x264 project
6 * Authors: Laurent Aimar <fenrir@via.ecp.fr>
7 * Loren Merritt <lorenm@u.washington.edu>
8 * Fiona Glaser <fiona@x264.com>
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.
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.
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.
24 * This program is also available under a commercial proprietary license.
25 * For more information, contact us at licensing@x264.com.
26 *****************************************************************************/
30 static int align_stride( int x, int align, int disalign )
32 x = ALIGN( x, align );
33 if( !(x&(disalign-1)) )
38 static int align_plane_size( int x, int disalign )
40 if( !(x&(disalign-1)) )
45 x264_frame_t *x264_frame_new( x264_t *h, int b_fdec )
49 int i_mb_count = h->mb.i_mb_count;
50 int i_stride, i_width, i_lines;
51 int i_padv = PADV << h->param.b_interlaced;
52 int luma_plane_size, chroma_plane_size;
53 int align = h->param.cpu&X264_CPU_CACHELINE_64 ? 64 : h->param.cpu&X264_CPU_CACHELINE_32 ? 32 : 16;
54 int disalign = h->param.cpu&X264_CPU_ALTIVEC ? 1<<9 : 1<<10;
56 CHECKED_MALLOCZERO( frame, sizeof(x264_frame_t) );
58 /* allocate frame data (+64 for extra data for me) */
59 i_width = h->mb.i_mb_width*16;
60 i_lines = h->mb.i_mb_height*16;
61 i_stride = align_stride( i_width + 2*PADH, align, disalign );
64 for( int i = 0; i < 2; i++ )
66 frame->i_width[i] = i_width >> i;
67 frame->i_lines[i] = i_lines >> i;
68 frame->i_stride[i] = i_stride;
71 frame->i_width_lowres = frame->i_width[0]/2;
72 frame->i_lines_lowres = frame->i_lines[0]/2;
73 frame->i_stride_lowres = align_stride( frame->i_width_lowres + 2*PADH, align, disalign<<1 );
75 for( int i = 0; i < h->param.i_bframe + 2; i++ )
76 for( int j = 0; j < h->param.i_bframe + 2; j++ )
77 CHECKED_MALLOC( frame->i_row_satds[i][j], i_lines/16 * sizeof(int) );
80 frame->i_type = X264_TYPE_AUTO;
84 frame->i_frame_num = -1;
85 frame->i_lines_completed = -1;
86 frame->b_fdec = b_fdec;
87 frame->i_pic_struct = PIC_STRUCT_AUTO;
88 frame->i_field_cnt = -1;
90 frame->i_cpb_duration =
91 frame->i_dpb_output_delay =
92 frame->i_cpb_delay = 0;
93 frame->i_coded_fields_lookahead =
94 frame->i_cpb_delay_lookahead = -1;
98 luma_plane_size = align_plane_size( frame->i_stride[0] * (frame->i_lines[0] + 2*i_padv), disalign );
99 chroma_plane_size = (frame->i_stride[1] * (frame->i_lines[1] + i_padv));
101 CHECKED_MALLOC( frame->buffer[1], chroma_plane_size * sizeof(pixel) );
102 frame->plane[1] = frame->buffer[1] + frame->i_stride[1] * i_padv/2 + PADH;
104 /* all 4 luma planes allocated together, since the cacheline split code
105 * requires them to be in-phase wrt cacheline alignment. */
106 if( h->param.analyse.i_subpel_refine && b_fdec )
108 CHECKED_MALLOC( frame->buffer[0], 4*luma_plane_size * sizeof(pixel) );
109 for( int i = 0; i < 4; i++ )
110 frame->filtered[i] = frame->buffer[0] + i*luma_plane_size + frame->i_stride[0] * i_padv + PADH;
111 frame->plane[0] = frame->filtered[0];
115 CHECKED_MALLOC( frame->buffer[0], luma_plane_size * sizeof(pixel) );
116 frame->filtered[0] = frame->plane[0] = frame->buffer[0] + frame->i_stride[0] * i_padv + PADH;
119 frame->b_duplicate = 0;
121 if( b_fdec ) /* fdec frame */
123 CHECKED_MALLOC( frame->mb_type, i_mb_count * sizeof(int8_t));
124 CHECKED_MALLOC( frame->mb_partition, i_mb_count * sizeof(uint8_t));
125 CHECKED_MALLOC( frame->mv[0], 2*16 * i_mb_count * sizeof(int16_t) );
126 CHECKED_MALLOC( frame->mv16x16, 2*(i_mb_count+1) * sizeof(int16_t) );
127 M32( frame->mv16x16[0] ) = 0;
129 CHECKED_MALLOC( frame->ref[0], 4 * i_mb_count * sizeof(int8_t) );
130 if( h->param.i_bframe )
132 CHECKED_MALLOC( frame->mv[1], 2*16 * i_mb_count * sizeof(int16_t) );
133 CHECKED_MALLOC( frame->ref[1], 4 * i_mb_count * sizeof(int8_t) );
138 frame->ref[1] = NULL;
140 CHECKED_MALLOC( frame->i_row_bits, i_lines/16 * sizeof(int) );
141 CHECKED_MALLOC( frame->f_row_qp, i_lines/16 * sizeof(float) );
142 if( h->param.analyse.i_me_method >= X264_ME_ESA )
144 CHECKED_MALLOC( frame->buffer[3],
145 frame->i_stride[0] * (frame->i_lines[0] + 2*i_padv) * sizeof(uint16_t) << h->frames.b_have_sub8x8_esa );
146 frame->integral = (uint16_t*)frame->buffer[3] + frame->i_stride[0] * i_padv + PADH;
149 else /* fenc frame */
151 if( h->frames.b_have_lowres )
153 luma_plane_size = align_plane_size( frame->i_stride_lowres * (frame->i_lines[0]/2 + 2*PADV), disalign );
155 CHECKED_MALLOC( frame->buffer_lowres[0], 4 * luma_plane_size * sizeof(pixel) );
156 for( int i = 0; i < 4; i++ )
157 frame->lowres[i] = frame->buffer_lowres[0] + (frame->i_stride_lowres * PADV + PADH) + i * luma_plane_size;
159 for( int j = 0; j <= !!h->param.i_bframe; j++ )
160 for( int i = 0; i <= h->param.i_bframe; i++ )
162 CHECKED_MALLOCZERO( frame->lowres_mvs[j][i], 2*h->mb.i_mb_count*sizeof(int16_t) );
163 CHECKED_MALLOC( frame->lowres_mv_costs[j][i], h->mb.i_mb_count*sizeof(int) );
165 CHECKED_MALLOC( frame->i_propagate_cost, (i_mb_count+3) * sizeof(uint16_t) );
166 for( int j = 0; j <= h->param.i_bframe+1; j++ )
167 for( int i = 0; i <= h->param.i_bframe+1; i++ )
168 CHECKED_MALLOC( frame->lowres_costs[j][i], (i_mb_count+3) * sizeof(uint16_t) );
169 frame->i_intra_cost = frame->lowres_costs[0][0];
170 memset( frame->i_intra_cost, -1, (i_mb_count+3) * sizeof(uint16_t) );
172 if( h->param.rc.i_aq_mode )
174 CHECKED_MALLOC( frame->f_qp_offset, h->mb.i_mb_count * sizeof(float) );
175 CHECKED_MALLOC( frame->f_qp_offset_aq, h->mb.i_mb_count * sizeof(float) );
176 if( h->frames.b_have_lowres )
177 /* shouldn't really be initialized, just silences a valgrind false-positive in x264_mbtree_propagate_cost_sse2 */
178 CHECKED_MALLOCZERO( frame->i_inv_qscale_factor, (h->mb.i_mb_count+3) * sizeof(uint16_t) );
182 if( x264_pthread_mutex_init( &frame->mutex, NULL ) )
184 if( x264_pthread_cond_init( &frame->cv, NULL ) )
194 void x264_frame_delete( x264_frame_t *frame )
196 /* Duplicate frames are blank copies of real frames (including pointers),
197 * so freeing those pointers would cause a double free later. */
198 if( !frame->b_duplicate )
200 for( int i = 0; i < 4; i++ )
201 x264_free( frame->buffer[i] );
202 for( int i = 0; i < 4; i++ )
203 x264_free( frame->buffer_lowres[i] );
204 for( int i = 0; i < X264_BFRAME_MAX+2; i++ )
205 for( int j = 0; j < X264_BFRAME_MAX+2; j++ )
206 x264_free( frame->i_row_satds[i][j] );
207 for( int j = 0; j < 2; j++ )
208 for( int i = 0; i <= X264_BFRAME_MAX; i++ )
210 x264_free( frame->lowres_mvs[j][i] );
211 x264_free( frame->lowres_mv_costs[j][i] );
213 x264_free( frame->i_propagate_cost );
214 for( int j = 0; j <= X264_BFRAME_MAX+1; j++ )
215 for( int i = 0; i <= X264_BFRAME_MAX+1; i++ )
216 x264_free( frame->lowres_costs[j][i] );
217 x264_free( frame->f_qp_offset );
218 x264_free( frame->f_qp_offset_aq );
219 x264_free( frame->i_inv_qscale_factor );
220 x264_free( frame->i_row_bits );
221 x264_free( frame->f_row_qp );
222 x264_free( frame->mb_type );
223 x264_free( frame->mb_partition );
224 x264_free( frame->mv[0] );
225 x264_free( frame->mv[1] );
227 x264_free( frame->mv16x16-1 );
228 x264_free( frame->ref[0] );
229 x264_free( frame->ref[1] );
230 x264_pthread_mutex_destroy( &frame->mutex );
231 x264_pthread_cond_destroy( &frame->cv );
236 static int get_plane_ptr( x264_t *h, x264_picture_t *src, uint8_t **pix, int *stride, int plane, int xshift, int yshift )
238 int width = h->param.i_width >> xshift;
239 int height = h->param.i_height >> yshift;
240 *pix = src->img.plane[plane];
241 *stride = src->img.i_stride[plane];
242 if( src->img.i_csp & X264_CSP_VFLIP )
244 *pix += (height-1) * *stride;
247 if( width > abs(*stride) )
249 x264_log( h, X264_LOG_ERROR, "Input picture width (%d) is greater than stride (%d)\n", width, *stride );
255 #define get_plane_ptr(...) do{ if( get_plane_ptr(__VA_ARGS__) < 0 ) return -1; }while(0)
257 int x264_frame_copy_picture( x264_t *h, x264_frame_t *dst, x264_picture_t *src )
259 int i_csp = src->img.i_csp & X264_CSP_MASK;
260 if( i_csp <= X264_CSP_NONE || i_csp >= X264_CSP_MAX )
262 x264_log( h, X264_LOG_ERROR, "Invalid input colorspace\n" );
266 dst->i_type = src->i_type;
267 dst->i_qpplus1 = src->i_qpplus1;
268 dst->i_pts = dst->i_reordered_pts = src->i_pts;
269 dst->param = src->param;
270 dst->i_pic_struct = src->i_pic_struct;
271 dst->extra_sei = src->extra_sei;
275 get_plane_ptr( h, src, &pix[0], &stride[0], 0, 0, 0 );
276 h->mc.plane_copy( dst->plane[0], dst->i_stride[0], pix[0], stride[0],
277 h->param.i_width, h->param.i_height );
278 if( i_csp == X264_CSP_NV12 )
280 get_plane_ptr( h, src, &pix[1], &stride[1], 1, 0, 1 );
281 h->mc.plane_copy( dst->plane[1], dst->i_stride[1], pix[1], stride[1],
282 h->param.i_width, h->param.i_height>>1 );
286 get_plane_ptr( h, src, &pix[1], &stride[1], i_csp==X264_CSP_I420 ? 1 : 2, 1, 1 );
287 get_plane_ptr( h, src, &pix[2], &stride[2], i_csp==X264_CSP_I420 ? 2 : 1, 1, 1 );
288 h->mc.plane_copy_interleave( dst->plane[1], dst->i_stride[1],
289 pix[1], stride[1], pix[2], stride[2],
290 h->param.i_width>>1, h->param.i_height>>1 );
295 static void ALWAYS_INLINE pixel_memset( pixel *dst, pixel *src, int len, int size )
297 uint8_t *dstp = (uint8_t*)dst;
299 memset(dst, *src, len);
300 } else if(size == 2) {
302 for(int i=0; i<len; i++)
304 } else if(size == 4) {
306 for(int i=0; i<len; i++)
311 static void plane_expand_border( pixel *pix, int i_stride, int i_width, int i_height, int i_padh, int i_padv, int b_pad_top, int b_pad_bottom, int b_chroma )
313 #define PPIXEL(x, y) ( pix + (x) + (y)*i_stride )
314 for( int y = 0; y < i_height; y++ )
317 pixel_memset( PPIXEL(-i_padh, y), PPIXEL(0, y), i_padh>>b_chroma, sizeof(pixel)<<b_chroma );
319 pixel_memset( PPIXEL(i_width, y), PPIXEL(i_width-1-b_chroma, y), i_padh>>b_chroma, sizeof(pixel)<<b_chroma );
323 for( int y = 0; y < i_padv; y++ )
324 memcpy( PPIXEL(-i_padh, -y-1), PPIXEL(-i_padh, 0), (i_width+2*i_padh) * sizeof(pixel) );
327 for( int y = 0; y < i_padv; y++ )
328 memcpy( PPIXEL(-i_padh, i_height+y), PPIXEL(-i_padh, i_height-1), (i_width+2*i_padh) * sizeof(pixel) );
332 void x264_frame_expand_border( x264_t *h, x264_frame_t *frame, int mb_y, int b_end )
335 if( mb_y & h->sh.b_mbaff )
337 for( int i = 0; i < frame->i_plane; i++ )
339 int stride = frame->i_stride[i];
340 int width = 16*h->sps->i_mb_width;
341 int height = (b_end ? 16*(h->mb.i_mb_height - mb_y) >> h->sh.b_mbaff : 16) >> !!i;
343 int padv = PADV >> !!i;
344 // buffer: 2 chroma, 3 luma (rounded to 4) because deblocking goes beyond the top of the mb
345 pixel *pix = frame->plane[i] + X264_MAX(0, (16*mb_y-4)*stride >> !!i);
346 if( b_end && !b_start )
347 height += 4 >> (!!i + h->sh.b_mbaff);
350 plane_expand_border( pix, stride*2, width, height, padh, padv, b_start, b_end, i );
351 plane_expand_border( pix+stride, stride*2, width, height, padh, padv, b_start, b_end, i );
355 plane_expand_border( pix, stride, width, height, padh, padv, b_start, b_end, i );
360 void x264_frame_expand_border_filtered( x264_t *h, x264_frame_t *frame, int mb_y, int b_end )
362 /* during filtering, 8 extra pixels were filtered on each edge,
363 * but up to 3 of the horizontal ones may be wrong.
364 we want to expand border from the last filtered pixel */
366 int stride = frame->i_stride[0];
367 int width = 16*h->mb.i_mb_width + 8;
368 int height = b_end ? (16*(h->mb.i_mb_height - mb_y) >> h->sh.b_mbaff) + 16 : 16;
371 for( int i = 1; i < 4; i++ )
373 // buffer: 8 luma, to match the hpel filter
374 pixel *pix = frame->filtered[i] + (16*mb_y - (8 << h->sh.b_mbaff)) * stride - 4;
377 plane_expand_border( pix, stride*2, width, height, padh, padv, b_start, b_end, 0 );
378 plane_expand_border( pix+stride, stride*2, width, height, padh, padv, b_start, b_end, 0 );
381 plane_expand_border( pix, stride, width, height, padh, padv, b_start, b_end, 0 );
385 void x264_frame_expand_border_lowres( x264_frame_t *frame )
387 for( int i = 0; i < 4; i++ )
388 plane_expand_border( frame->lowres[i], frame->i_stride_lowres, frame->i_width_lowres, frame->i_lines_lowres, PADH, PADV, 1, 1, 0 );
391 void x264_frame_expand_border_mod16( x264_t *h, x264_frame_t *frame )
393 for( int i = 0; i < frame->i_plane; i++ )
395 int i_width = h->param.i_width;
396 int i_height = h->param.i_height >> !!i;
397 int i_padx = (h->mb.i_mb_width * 16 - h->param.i_width);
398 int i_pady = (h->mb.i_mb_height * 16 - h->param.i_height) >> !!i;
402 for( int y = 0; y < i_height; y++ )
403 pixel_memset( &frame->plane[i][y*frame->i_stride[i] + i_width],
404 &frame->plane[i][y*frame->i_stride[i] + i_width - 1-i],
405 i_padx>>i, sizeof(pixel)<<i );
409 for( int y = i_height; y < i_height + i_pady; y++ )
410 memcpy( &frame->plane[i][y*frame->i_stride[i]],
411 &frame->plane[i][(i_height-(~y&h->param.b_interlaced)-1)*frame->i_stride[i]],
412 (i_width + i_padx) * sizeof(pixel) );
418 void x264_frame_cond_broadcast( x264_frame_t *frame, int i_lines_completed )
420 x264_pthread_mutex_lock( &frame->mutex );
421 frame->i_lines_completed = i_lines_completed;
422 x264_pthread_cond_broadcast( &frame->cv );
423 x264_pthread_mutex_unlock( &frame->mutex );
426 void x264_frame_cond_wait( x264_frame_t *frame, int i_lines_completed )
428 x264_pthread_mutex_lock( &frame->mutex );
429 while( frame->i_lines_completed < i_lines_completed )
430 x264_pthread_cond_wait( &frame->cv, &frame->mutex );
431 x264_pthread_mutex_unlock( &frame->mutex );
436 void x264_frame_push( x264_frame_t **list, x264_frame_t *frame )
439 while( list[i] ) i++;
443 x264_frame_t *x264_frame_pop( x264_frame_t **list )
448 while( list[i+1] ) i++;
454 void x264_frame_unshift( x264_frame_t **list, x264_frame_t *frame )
457 while( list[i] ) i++;
463 x264_frame_t *x264_frame_shift( x264_frame_t **list )
465 x264_frame_t *frame = list[0];
467 for( i = 0; list[i]; i++ )
473 void x264_frame_push_unused( x264_t *h, x264_frame_t *frame )
475 assert( frame->i_reference_count > 0 );
476 frame->i_reference_count--;
477 if( frame->i_reference_count == 0 )
478 x264_frame_push( h->frames.unused[frame->b_fdec], frame );
481 x264_frame_t *x264_frame_pop_unused( x264_t *h, int b_fdec )
484 if( h->frames.unused[b_fdec][0] )
485 frame = x264_frame_pop( h->frames.unused[b_fdec] );
487 frame = x264_frame_new( h, b_fdec );
490 frame->b_last_minigop_bframe = 0;
491 frame->i_reference_count = 1;
492 frame->b_intra_calculated = 0;
493 frame->b_scenecut = 1;
494 frame->b_keyframe = 0;
495 frame->b_corrupt = 0;
497 memset( frame->weight, 0, sizeof(frame->weight) );
498 memset( frame->f_weighted_cost_delta, 0, sizeof(frame->f_weighted_cost_delta) );
503 void x264_frame_push_blank_unused( x264_t *h, x264_frame_t *frame )
505 assert( frame->i_reference_count > 0 );
506 frame->i_reference_count--;
507 if( frame->i_reference_count == 0 )
508 x264_frame_push( h->frames.blank_unused, frame );
511 x264_frame_t *x264_frame_pop_blank_unused( x264_t *h )
514 if( h->frames.blank_unused[0] )
515 frame = x264_frame_pop( h->frames.blank_unused );
517 frame = x264_malloc( sizeof(x264_frame_t) );
520 frame->b_duplicate = 1;
521 frame->i_reference_count = 1;
525 void x264_frame_sort( x264_frame_t **list, int b_dts )
530 for( int i = 0; list[i+1]; i++ )
532 int dtype = list[i]->i_type - list[i+1]->i_type;
533 int dtime = list[i]->i_frame - list[i+1]->i_frame;
534 int swap = b_dts ? dtype > 0 || ( dtype == 0 && dtime > 0 )
538 XCHG( x264_frame_t*, list[i], list[i+1] );
545 void x264_weight_scale_plane( x264_t *h, pixel *dst, int i_dst_stride, pixel *src, int i_src_stride,
546 int i_width, int i_height, x264_weight_t *w )
548 /* Weight horizontal strips of height 16. This was found to be the optimal height
549 * in terms of the cache loads. */
550 while( i_height > 0 )
552 for( int x = 0; x < i_width; x += 16 )
553 w->weightfn[16>>2]( dst+x, i_dst_stride, src+x, i_src_stride, w, X264_MIN( i_height, 16 ) );
555 dst += 16 * i_dst_stride;
556 src += 16 * i_src_stride;
560 void x264_frame_delete_list( x264_frame_t **list )
566 x264_frame_delete( list[i++] );
570 int x264_sync_frame_list_init( x264_sync_frame_list_t *slist, int max_size )
574 slist->i_max_size = max_size;
576 CHECKED_MALLOCZERO( slist->list, (max_size+1) * sizeof(x264_frame_t*) );
577 if( x264_pthread_mutex_init( &slist->mutex, NULL ) ||
578 x264_pthread_cond_init( &slist->cv_fill, NULL ) ||
579 x264_pthread_cond_init( &slist->cv_empty, NULL ) )
586 void x264_sync_frame_list_delete( x264_sync_frame_list_t *slist )
588 x264_pthread_mutex_destroy( &slist->mutex );
589 x264_pthread_cond_destroy( &slist->cv_fill );
590 x264_pthread_cond_destroy( &slist->cv_empty );
591 x264_frame_delete_list( slist->list );
594 void x264_sync_frame_list_push( x264_sync_frame_list_t *slist, x264_frame_t *frame )
596 x264_pthread_mutex_lock( &slist->mutex );
597 while( slist->i_size == slist->i_max_size )
598 x264_pthread_cond_wait( &slist->cv_empty, &slist->mutex );
599 slist->list[ slist->i_size++ ] = frame;
600 x264_pthread_mutex_unlock( &slist->mutex );
601 x264_pthread_cond_broadcast( &slist->cv_fill );
604 x264_frame_t *x264_sync_frame_list_pop( x264_sync_frame_list_t *slist )
607 x264_pthread_mutex_lock( &slist->mutex );
608 while( !slist->i_size )
609 x264_pthread_cond_wait( &slist->cv_fill, &slist->mutex );
610 frame = slist->list[ --slist->i_size ];
611 slist->list[ slist->i_size ] = NULL;
612 x264_pthread_cond_broadcast( &slist->cv_empty );
613 x264_pthread_mutex_unlock( &slist->mutex );